import jalview.datamodel.SequenceFeature;
import jalview.datamodel.SequenceI;
import jalview.ext.rbvi.chimera.JalviewChimeraBinding;
+import jalview.gui.AlignmentPanel;
import jalview.gui.Desktop;
import jalview.gui.StructureViewer.ViewerType;
import jalview.io.DataSourceType;
import jalview.schemes.ResidueProperties;
import jalview.structure.AtomSpec;
import jalview.structure.AtomSpecModel;
+import jalview.structure.StructureCommand;
import jalview.structure.StructureCommandI;
import jalview.structure.StructureCommandsI;
import jalview.structure.StructureListener;
public static class SuperposeData
{
public String filename;
-
+
public String pdbId;
-
+
public String chain = "";
-
+
public boolean isRna;
-
+
/*
* The pdb residue number (if any) mapped to columns of the alignment
*/
public int[] pdbResNo; // or use SparseIntArray?
-
+
public String modelId;
-
+
/**
* Constructor
*
chains = newchains;
return chainmaps > 0;
}
+
public StructureSelectionManager getSsm()
{
return ssm;
* @return
*/
protected int findSuperposableResidues(AlignmentI alignment,
- BitSet matched, AAStructureBindingModel.SuperposeData[] structures)
+ BitSet matched,
+ AAStructureBindingModel.SuperposeData[] structures)
{
int refStructure = -1;
String[] files = getStructureFiles();
* Calculate the superposable alignment columns ('matched'), and the
* corresponding structure residue positions (structures.pdbResNo)
*/
- int refStructure = findSuperposableResidues(alignment,
- matched, structures);
+ int refStructure = findSuperposableResidues(alignment, matched,
+ structures);
/*
* require at least 4 positions to be able to execute superposition
int nmatched = matched.cardinality();
if (nmatched < MIN_POS_TO_SUPERPOSE)
{
- String msg = MessageManager.formatMessage("label.insufficient_residues",
- nmatched);
+ String msg = MessageManager
+ .formatMessage("label.insufficient_residues", nmatched);
error += view.getViewName() + ": " + msg + "; ";
continue;
}
return error;
}
- private AtomSpecModel getAtomSpec(AAStructureBindingModel.SuperposeData superposeData,
+ private AtomSpecModel getAtomSpec(
+ AAStructureBindingModel.SuperposeData superposeData,
BitSet matched)
{
AtomSpecModel model = new AtomSpecModel();
{
return;
}
-
+
/*
* build a map of {Residue3LetterCode, Color}
*/
*/
final JalviewStructureDisplayI theViewer = getViewer();
final long handle = msg == null ? 0 : theViewer.startProgressBar(msg);
-
+
SwingUtilities.invokeLater(new Runnable()
{
@Override
*/
final JalviewStructureDisplayI theViewer = getViewer();
final long handle = msg == null ? 0 : theViewer.startProgressBar(msg);
-
+
List<String> response = getReply ? new ArrayList<>() : null;
try
{
AlignmentI al = viewport.getAlignment();
Map<Object, AtomSpecModel> colourMap = new LinkedHashMap<>();
Color lastColour = null;
-
+
for (int pdbfnum = 0; pdbfnum < files.length; pdbfnum++)
{
final String modelId = getModelIdForFile(files[pdbfnum]);
StructureMapping[] mapping = ssm.getMapping(files[pdbfnum]);
-
+
if (mapping == null || mapping.length < 1)
{
continue;
}
-
+
int startPos = -1, lastPos = -1;
String lastChain = "";
for (int s = 0; s < sequence[pdbfnum].length; s++)
continue;
}
int pos = mapping[m].getPDBResNum(asp.findPosition(r));
-
+
if (pos < 1 || pos == lastPos)
{
continue;
}
-
+
Color colour = sr.getResidueColour(seq, r, finder);
-
+
/*
* darker colour for hidden regions
*/
{
colour = Color.GRAY;
}
-
+
final String chain = mapping[m].getChain();
-
+
/*
* Just keep incrementing the end position for this colour range
* _unless_ colour, PDB model or chain has changed, or there is a
{
if (startPos != -1)
{
- addAtomSpecRange(colourMap, lastColour, modelId,
- startPos, lastPos, lastChain);
+ addAtomSpecRange(colourMap, lastColour, modelId, startPos,
+ lastPos, lastChain);
}
startPos = pos;
}
}
/**
- * Helper method to add one contiguous range to the AtomSpec model for the given
- * value (creating the model if necessary). As used by Jalview, {@code value} is
+ * Helper method to add one contiguous range to the AtomSpec model for the
+ * given value (creating the model if necessary). As used by Jalview,
+ * {@code value} is
* <ul>
* <li>a colour, when building a 'colour structure by sequence' command</li>
* <li>a feature value, when building a 'set Chimera attributes from features'
* @param chain
*/
public static final void addAtomSpecRange(Map<Object, AtomSpecModel> map,
- Object value,
- String model, int startPos, int endPos, String chain)
+ Object value, String model, int startPos, int endPos,
+ String chain)
{
/*
* Get/initialize map of data for the colour
atomSpec = new AtomSpecModel();
map.put(value, atomSpec);
}
-
+
atomSpec.addRange(model, startPos, endPos, chain);
}
saveSession(f);
} catch (IOException e)
{
- Cache.log.error(String.format("Error saving %s session: %s",
- prefix, e.toString()));
+ Cache.log.error(String.format("Error saving %s session: %s", prefix,
+ e.toString()));
}
return f;
*/
protected void saveSession(File f)
{
- StructureCommandI cmd = commandGenerator
- .saveSession(f.getPath());
+ StructureCommandI cmd = commandGenerator.saveSession(f.getPath());
if (cmd != null)
{
executeCommand(cmd, false);
}
stopListening();
-
+
if (forceClose)
{
StructureCommandI cmd = getCommandGenerator().closeViewer();
{
return theMap;
}
-
+
AlignViewportI viewport = viewPanel.getAlignViewport();
List<String> visibleFeatures = fr.getDisplayedFeatureTypes();
-
+
/*
* if alignment is showing features from complement, we also transfer
* these features to the corresponding mapped structure residues
{
return theMap;
}
-
+
AlignmentI alignment = viewPanel.getAlignment();
SequenceI[][] seqs = getSequence();
{
String modelId = getModelIdForFile(files[pdbfnum]);
StructureMapping[] mapping = ssm.getMapping(files[pdbfnum]);
-
+
if (mapping == null || mapping.length < 1)
{
continue;
}
-
+
for (int seqNo = 0; seqNo < seqs[pdbfnum].length; seqNo++)
{
for (int m = 0; m < mapping.length; m++)
}
/**
- * Scans visible features in mapped positions of the CDS/peptide complement, and
- * adds any found to the map of attribute values/structure positions
+ * Scans visible features in mapped positions of the CDS/peptide complement,
+ * and adds any found to the map of attribute values/structure positions
*
* @param complementRenderer
* @param structureMapping
for (SequenceFeature sf : mf.features)
{
String type = sf.getType();
-
+
/*
* Don't copy features which originated from Chimera
*/
{
continue;
}
-
+
/*
* record feature 'value' (score/description/type) as at the
* corresponding structure position
*/
List<int[]> mappedRanges = structureMapping
.getPDBResNumRanges(seqPos, seqPos);
-
+
if (!mappedRanges.isEmpty())
{
String value = sf.getDescription();
for (SequenceFeature sf : sfs)
{
String type = sf.getType();
-
+
/*
* Don't copy features which originated from Chimera
*/
{
continue;
}
-
+
List<int[]> mappedRanges = mapping.getPDBResNumRanges(sf.getBegin(),
sf.getEnd());
-
+
if (!mappedRanges.isEmpty())
{
String value = sf.getDescription();
{
externalViewerMonitor = new Thread(new Runnable()
{
-
+
@Override
public void run()
{
executeCommands(commands, false, null);
}
}
+
+ /**
+ * If supported by the structure viewer, queries it for all residue attributes
+ * with the given attribute name, and creates features on corresponding
+ * residues of the alignment. Returns the number of features added.
+ *
+ * @param attName
+ * @param alignmentPanel
+ * @return
+ */
+ public int copyStructureAttributesToFeatures(String attName,
+ AlignmentPanel alignmentPanel)
+ {
+ StructureCommandI cmd = getCommandGenerator()
+ .getResidueAttributes(attName);
+ if (cmd == null)
+ {
+ return 0;
+ }
+ List<String> residueAttributes = executeCommand(cmd, true);
+
+ int featuresAdded = createFeaturesForAttributes(attName,
+ residueAttributes);
+ if (featuresAdded > 0)
+ {
+ alignmentPanel.getFeatureRenderer().featuresAdded();
+ }
+ return featuresAdded;
+ }
+
+ /**
+ * Parses {@code residueAttributes} and creates sequence features on any
+ * mapped alignment residues. Returns the number of features created.
+ * <p>
+ * {@code residueAttributes} is the reply from the structure viewer to a
+ * command to list any residue attributes for the given attribute name. Syntax
+ * and parsing of this is viewer-specific.
+ *
+ * @param attName
+ * @param residueAttributes
+ * @return
+ */
+ protected int createFeaturesForAttributes(String attName,
+ List<String> residueAttributes)
+ {
+ return 0;
+ }
}