/*
* Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$)
* Copyright (C) $$Year-Rel$$ The Jalview Authors
*
* This file is part of Jalview.
*
* Jalview is free software: you can redistribute it and/or
* modify it under the terms of the GNU General Public License
* as published by the Free Software Foundation, either version 3
* of the License, or (at your option) any later version.
*
* Jalview is distributed in the hope that it will be useful, but
* WITHOUT ANY WARRANTY; without even the implied warranty
* of MERCHANTABILITY or FITNESS FOR A PARTICULAR
* PURPOSE. See the GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with Jalview. If not, see .
* The Jalview Authors are detailed in the 'AUTHORS' file.
*/
package jalview.gui;
import jalview.api.structures.JalviewStructureDisplayI;
import jalview.bin.Cache;
import jalview.datamodel.AlignedCodonFrame;
import jalview.datamodel.Alignment;
import jalview.datamodel.AlignmentAnnotation;
import jalview.datamodel.AlignmentI;
import jalview.datamodel.PDBEntry;
import jalview.datamodel.SequenceI;
import jalview.datamodel.StructureViewerModel;
import jalview.datamodel.StructureViewerModel.StructureData;
import jalview.schemabinding.version2.AlcodMap;
import jalview.schemabinding.version2.AlcodonFrame;
import jalview.schemabinding.version2.Annotation;
import jalview.schemabinding.version2.AnnotationColours;
import jalview.schemabinding.version2.AnnotationElement;
import jalview.schemabinding.version2.CalcIdParam;
import jalview.schemabinding.version2.DBRef;
import jalview.schemabinding.version2.Features;
import jalview.schemabinding.version2.Group;
import jalview.schemabinding.version2.HiddenColumns;
import jalview.schemabinding.version2.JGroup;
import jalview.schemabinding.version2.JSeq;
import jalview.schemabinding.version2.JalviewModel;
import jalview.schemabinding.version2.JalviewModelSequence;
import jalview.schemabinding.version2.MapListFrom;
import jalview.schemabinding.version2.MapListTo;
import jalview.schemabinding.version2.Mapping;
import jalview.schemabinding.version2.MappingChoice;
import jalview.schemabinding.version2.OtherData;
import jalview.schemabinding.version2.PdbentryItem;
import jalview.schemabinding.version2.Pdbids;
import jalview.schemabinding.version2.Property;
import jalview.schemabinding.version2.Sequence;
import jalview.schemabinding.version2.SequenceSet;
import jalview.schemabinding.version2.SequenceSetProperties;
import jalview.schemabinding.version2.Setting;
import jalview.schemabinding.version2.StructureState;
import jalview.schemabinding.version2.ThresholdLine;
import jalview.schemabinding.version2.Tree;
import jalview.schemabinding.version2.UserColours;
import jalview.schemabinding.version2.Viewport;
import jalview.schemes.AnnotationColourGradient;
import jalview.schemes.ColourSchemeI;
import jalview.schemes.ColourSchemeProperty;
import jalview.schemes.GraduatedColor;
import jalview.schemes.ResidueColourScheme;
import jalview.schemes.ResidueProperties;
import jalview.schemes.UserColourScheme;
import jalview.structure.StructureSelectionManager;
import jalview.structures.models.AAStructureBindingModel;
import jalview.util.MessageManager;
import jalview.util.Platform;
import jalview.util.jarInputStreamProvider;
import jalview.viewmodel.AlignmentViewport;
import jalview.viewmodel.seqfeatures.FeatureRendererSettings;
import jalview.viewmodel.seqfeatures.FeaturesDisplayed;
import jalview.ws.jws2.Jws2Discoverer;
import jalview.ws.jws2.dm.AAConSettings;
import jalview.ws.jws2.jabaws2.Jws2Instance;
import jalview.ws.params.ArgumentI;
import jalview.ws.params.AutoCalcSetting;
import jalview.ws.params.WsParamSetI;
import java.awt.Rectangle;
import java.io.BufferedReader;
import java.io.DataInputStream;
import java.io.DataOutputStream;
import java.io.File;
import java.io.FileInputStream;
import java.io.FileOutputStream;
import java.io.IOException;
import java.io.InputStreamReader;
import java.io.OutputStreamWriter;
import java.io.PrintWriter;
import java.lang.reflect.InvocationTargetException;
import java.net.MalformedURLException;
import java.net.URL;
import java.util.ArrayList;
import java.util.Enumeration;
import java.util.HashMap;
import java.util.HashSet;
import java.util.Hashtable;
import java.util.IdentityHashMap;
import java.util.Iterator;
import java.util.LinkedHashMap;
import java.util.List;
import java.util.Map;
import java.util.Map.Entry;
import java.util.Set;
import java.util.StringTokenizer;
import java.util.Vector;
import java.util.jar.JarEntry;
import java.util.jar.JarInputStream;
import java.util.jar.JarOutputStream;
import javax.swing.JInternalFrame;
import javax.swing.JOptionPane;
import javax.swing.SwingUtilities;
import org.exolab.castor.xml.Unmarshaller;
/**
* Write out the current jalview desktop state as a Jalview XML stream.
*
* Note: the vamsas objects referred to here are primitive versions of the
* VAMSAS project schema elements - they are not the same and most likely never
* will be :)
*
* @author $author$
* @version $Revision: 1.134 $
*/
public class Jalview2XML
{
/*
* SequenceI reference -> XML ID string in jalview XML. Populated as XML reps
* of sequence objects are created.
*/
IdentityHashMap seqsToIds = null;
/**
* jalview XML Sequence ID to jalview sequence object reference (both dataset
* and alignment sequences. Populated as XML reps of sequence objects are
* created.)
*/
Map seqRefIds = null;
Vector frefedSequence = null;
boolean raiseGUI = true; // whether errors are raised in dialog boxes or not
/*
* Map of reconstructed AlignFrame objects that appear to have come from
* SplitFrame objects (have a dna/protein complement view).
*/
private Map splitFrameCandidates = new HashMap();
/**
* create/return unique hash string for sq
*
* @param sq
* @return new or existing unique string for sq
*/
String seqHash(SequenceI sq)
{
if (seqsToIds == null)
{
initSeqRefs();
}
if (seqsToIds.containsKey(sq))
{
return seqsToIds.get(sq);
}
else
{
// create sequential key
String key = "sq" + (seqsToIds.size() + 1);
key = makeHashCode(sq, key); // check we don't have an external reference
// for it already.
seqsToIds.put(sq, key);
return key;
}
}
void clearSeqRefs()
{
if (_cleartables)
{
if (seqRefIds != null)
{
seqRefIds.clear();
}
if (seqsToIds != null)
{
seqsToIds.clear();
}
// seqRefIds = null;
// seqsToIds = null;
}
else
{
// do nothing
warn("clearSeqRefs called when _cleartables was not set. Doing nothing.");
// seqRefIds = new Hashtable();
// seqsToIds = new IdentityHashMap();
}
}
void initSeqRefs()
{
if (seqsToIds == null)
{
seqsToIds = new IdentityHashMap();
}
if (seqRefIds == null)
{
seqRefIds = new HashMap();
}
}
public Jalview2XML()
{
}
public Jalview2XML(boolean raiseGUI)
{
this.raiseGUI = raiseGUI;
}
public void resolveFrefedSequences()
{
if (frefedSequence.size() > 0)
{
int r = 0, rSize = frefedSequence.size();
while (r < rSize)
{
Object[] ref = (Object[]) frefedSequence.elementAt(r);
if (ref != null)
{
String sref = (String) ref[0];
if (seqRefIds.containsKey(sref))
{
if (ref[1] instanceof jalview.datamodel.Mapping)
{
SequenceI seq = seqRefIds.get(sref);
while (seq.getDatasetSequence() != null)
{
seq = seq.getDatasetSequence();
}
((jalview.datamodel.Mapping) ref[1]).setTo(seq);
}
else
{
if (ref[1] instanceof jalview.datamodel.AlignedCodonFrame)
{
SequenceI seq = seqRefIds.get(sref);
while (seq.getDatasetSequence() != null)
{
seq = seq.getDatasetSequence();
}
if (ref[2] != null
&& ref[2] instanceof jalview.datamodel.Mapping)
{
jalview.datamodel.Mapping mp = (jalview.datamodel.Mapping) ref[2];
((jalview.datamodel.AlignedCodonFrame) ref[1]).addMap(
seq, mp.getTo(), mp.getMap());
}
else
{
System.err
.println("IMPLEMENTATION ERROR: Unimplemented forward sequence references for AlcodonFrames involving "
+ ref[2].getClass() + " type objects.");
}
}
else
{
System.err
.println("IMPLEMENTATION ERROR: Unimplemented forward sequence references for "
+ ref[1].getClass() + " type objects.");
}
}
frefedSequence.remove(r);
rSize--;
}
else
{
System.err
.println("IMPLEMENTATION WARNING: Unresolved forward reference for hash string "
+ ref[0]
+ " with objecttype "
+ ref[1].getClass());
r++;
}
}
else
{
// empty reference
frefedSequence.remove(r);
rSize--;
}
}
}
}
/**
* This maintains a map of viewports, the key being the seqSetId. Important to
* set historyItem and redoList for multiple views
*/
Map viewportsAdded = new HashMap();
Map annotationIds = new HashMap();
String uniqueSetSuffix = "";
/**
* List of pdbfiles added to Jar
*/
List pdbfiles = null;
// SAVES SEVERAL ALIGNMENT WINDOWS TO SAME JARFILE
public void saveState(File statefile)
{
FileOutputStream fos = null;
try
{
fos = new FileOutputStream(statefile);
JarOutputStream jout = new JarOutputStream(fos);
saveState(jout);
} catch (Exception e)
{
// TODO: inform user of the problem - they need to know if their data was
// not saved !
if (errorMessage == null)
{
errorMessage = "Couldn't write Jalview Archive to output file '"
+ statefile + "' - See console error log for details";
}
else
{
errorMessage += "(output file was '" + statefile + "')";
}
e.printStackTrace();
} finally
{
if (fos != null)
{
try
{
fos.close();
} catch (IOException e)
{
// ignore
}
}
}
reportErrors();
}
/**
* Writes a jalview project archive to the given Jar output stream.
*
* @param jout
*/
public void saveState(JarOutputStream jout)
{
AlignFrame[] frames = Desktop.getAlignFrames(); // Desktop.desktop.getAllFrames();
if (frames == null)
{
return;
}
Hashtable dsses = new Hashtable();
try
{
// NOTE UTF-8 MUST BE USED FOR WRITING UNICODE CHARS
// //////////////////////////////////////////////////
Vector shortNames = new Vector();
// REVERSE ORDER
for (int i = frames.length - 1; i > -1; i--)
{
AlignFrame af = frames[i];
// skip ?
if (skipList != null
&& skipList
.containsKey(af.getViewport().getSequenceSetId()))
{
continue;
}
String shortName = af.getTitle();
if (shortName.indexOf(File.separatorChar) > -1)
{
shortName = shortName.substring(shortName
.lastIndexOf(File.separatorChar) + 1);
}
int count = 1;
while (shortNames.contains(shortName))
{
if (shortName.endsWith("_" + (count - 1)))
{
shortName = shortName.substring(0, shortName.lastIndexOf("_"));
}
shortName = shortName.concat("_" + count);
count++;
}
shortNames.addElement(shortName);
if (!shortName.endsWith(".xml"))
{
shortName = shortName + ".xml";
}
int ap, apSize = af.alignPanels.size();
for (ap = 0; ap < apSize; ap++)
{
AlignmentPanel apanel = af.alignPanels.get(ap);
String fileName = apSize == 1 ? shortName : ap + shortName;
if (!fileName.endsWith(".xml"))
{
fileName = fileName + ".xml";
}
saveState(apanel, fileName, jout);
String dssid = getDatasetIdRef(af.getViewport().getAlignment()
.getDataset());
if (!dsses.containsKey(dssid))
{
dsses.put(dssid, af);
}
}
}
writeDatasetFor(dsses, "" + jout.hashCode() + " " + uniqueSetSuffix,
jout);
try
{
jout.flush();
} catch (Exception foo)
{
}
;
jout.close();
} catch (Exception ex)
{
// TODO: inform user of the problem - they need to know if their data was
// not saved !
if (errorMessage == null)
{
errorMessage = "Couldn't write Jalview Archive - see error output for details";
}
ex.printStackTrace();
}
}
// USE THIS METHOD TO SAVE A SINGLE ALIGNMENT WINDOW
public boolean saveAlignment(AlignFrame af, String jarFile,
String fileName)
{
try
{
int ap = 0;
int apSize = af.alignPanels.size();
FileOutputStream fos = new FileOutputStream(jarFile);
JarOutputStream jout = new JarOutputStream(fos);
Hashtable dsses = new Hashtable();
for (AlignmentPanel apanel : af.alignPanels)
{
String jfileName = apSize == 1 ? fileName : fileName + ap;
ap++;
if (!jfileName.endsWith(".xml"))
{
jfileName = jfileName + ".xml";
}
saveState(apanel, jfileName, jout);
String dssid = getDatasetIdRef(af.getViewport().getAlignment()
.getDataset());
if (!dsses.containsKey(dssid))
{
dsses.put(dssid, af);
}
}
writeDatasetFor(dsses, fileName, jout);
try
{
jout.flush();
} catch (Exception foo)
{
}
;
jout.close();
return true;
} catch (Exception ex)
{
errorMessage = "Couldn't Write alignment view to Jalview Archive - see error output for details";
ex.printStackTrace();
return false;
}
}
private void writeDatasetFor(Hashtable dsses,
String fileName, JarOutputStream jout)
{
for (String dssids : dsses.keySet())
{
AlignFrame _af = dsses.get(dssids);
String jfileName = fileName + " Dataset for " + _af.getTitle();
if (!jfileName.endsWith(".xml"))
{
jfileName = jfileName + ".xml";
}
saveState(_af.alignPanel, jfileName, true, jout);
}
}
/**
* create a JalviewModel from an alignment view and marshall it to a
* JarOutputStream
*
* @param ap
* panel to create jalview model for
* @param fileName
* name of alignment panel written to output stream
* @param jout
* jar output stream
* @param out
* jar entry name
*/
public JalviewModel saveState(AlignmentPanel ap, String fileName,
JarOutputStream jout)
{
return saveState(ap, fileName, false, jout);
}
/**
* create a JalviewModel from an alignment view and marshall it to a
* JarOutputStream
*
* @param ap
* panel to create jalview model for
* @param fileName
* name of alignment panel written to output stream
* @param storeDS
* when true, only write the dataset for the alignment, not the data
* associated with the view.
* @param jout
* jar output stream
* @param out
* jar entry name
*/
public JalviewModel saveState(AlignmentPanel ap, String fileName,
boolean storeDS, JarOutputStream jout)
{
initSeqRefs();
List viewIds = new ArrayList();
List userColours = new ArrayList();
AlignViewport av = ap.av;
JalviewModel object = new JalviewModel();
object.setVamsasModel(new jalview.schemabinding.version2.VamsasModel());
object.setCreationDate(new java.util.Date(System.currentTimeMillis()));
object.setVersion(jalview.bin.Cache.getDefault("VERSION",
"Development Build"));
jalview.datamodel.AlignmentI jal = av.getAlignment();
if (av.hasHiddenRows())
{
jal = jal.getHiddenSequences().getFullAlignment();
}
SequenceSet vamsasSet = new SequenceSet();
Sequence vamsasSeq;
JalviewModelSequence jms = new JalviewModelSequence();
vamsasSet.setGapChar(jal.getGapCharacter() + "");
if (jal.getDataset() != null)
{
// dataset id is the dataset's hashcode
vamsasSet.setDatasetId(getDatasetIdRef(jal.getDataset()));
if (storeDS)
{
// switch jal and the dataset
jal = jal.getDataset();
}
}
if (jal.getProperties() != null)
{
Enumeration en = jal.getProperties().keys();
while (en.hasMoreElements())
{
String key = en.nextElement().toString();
SequenceSetProperties ssp = new SequenceSetProperties();
ssp.setKey(key);
ssp.setValue(jal.getProperties().get(key).toString());
vamsasSet.addSequenceSetProperties(ssp);
}
}
JSeq jseq;
Set calcIdSet = new HashSet();
// SAVE SEQUENCES
String id = "";
jalview.datamodel.SequenceI jds, jdatasq;
for (int i = 0; i < jal.getHeight(); i++)
{
jds = jal.getSequenceAt(i);
jdatasq = jds.getDatasetSequence() == null ? jds : jds
.getDatasetSequence();
id = seqHash(jds);
if (seqRefIds.get(id) != null)
{
// This happens for two reasons: 1. multiple views are being serialised.
// 2. the hashCode has collided with another sequence's code. This DOES
// HAPPEN! (PF00072.15.stk does this)
// JBPNote: Uncomment to debug writing out of files that do not read
// back in due to ArrayOutOfBoundExceptions.
// System.err.println("vamsasSeq backref: "+id+"");
// System.err.println(jds.getName()+"
// "+jds.getStart()+"-"+jds.getEnd()+" "+jds.getSequenceAsString());
// System.err.println("Hashcode: "+seqHash(jds));
// SequenceI rsq = (SequenceI) seqRefIds.get(id + "");
// System.err.println(rsq.getName()+"
// "+rsq.getStart()+"-"+rsq.getEnd()+" "+rsq.getSequenceAsString());
// System.err.println("Hashcode: "+seqHash(rsq));
}
else
{
vamsasSeq = createVamsasSequence(id, jds);
vamsasSet.addSequence(vamsasSeq);
seqRefIds.put(id, jds);
}
jseq = new JSeq();
jseq.setStart(jds.getStart());
jseq.setEnd(jds.getEnd());
jseq.setColour(av.getSequenceColour(jds).getRGB());
jseq.setId(id); // jseq id should be a string not a number
if (!storeDS)
{
// Store any sequences this sequence represents
if (av.hasHiddenRows())
{
jseq.setHidden(av.getAlignment().getHiddenSequences()
.isHidden(jds));
if (av.isHiddenRepSequence(jal.getSequenceAt(i)))
{
jalview.datamodel.SequenceI[] reps = av
.getRepresentedSequences(jal.getSequenceAt(i))
.getSequencesInOrder(jal);
for (int h = 0; h < reps.length; h++)
{
if (reps[h] != jal.getSequenceAt(i))
{
jseq.addHiddenSequences(jal.findIndex(reps[h]));
}
}
}
}
}
if (jds.getSequenceFeatures() != null)
{
jalview.datamodel.SequenceFeature[] sf = jds
.getSequenceFeatures();
int index = 0;
while (index < sf.length)
{
Features features = new Features();
features.setBegin(sf[index].getBegin());
features.setEnd(sf[index].getEnd());
features.setDescription(sf[index].getDescription());
features.setType(sf[index].getType());
features.setFeatureGroup(sf[index].getFeatureGroup());
features.setScore(sf[index].getScore());
if (sf[index].links != null)
{
for (int l = 0; l < sf[index].links.size(); l++)
{
OtherData keyValue = new OtherData();
keyValue.setKey("LINK_" + l);
keyValue.setValue(sf[index].links.elementAt(l).toString());
features.addOtherData(keyValue);
}
}
if (sf[index].otherDetails != null)
{
String key;
Enumeration keys = sf[index].otherDetails.keys();
while (keys.hasMoreElements())
{
key = keys.nextElement().toString();
OtherData keyValue = new OtherData();
keyValue.setKey(key);
keyValue.setValue(sf[index].otherDetails.get(key).toString());
features.addOtherData(keyValue);
}
}
jseq.addFeatures(features);
index++;
}
}
if (jdatasq.getPDBId() != null)
{
Enumeration en = jdatasq.getPDBId().elements();
while (en.hasMoreElements())
{
Pdbids pdb = new Pdbids();
jalview.datamodel.PDBEntry entry = (jalview.datamodel.PDBEntry) en
.nextElement();
pdb.setId(entry.getId());
pdb.setType(entry.getType());
/*
* Store any structure views associated with this sequence. This
* section copes with duplicate entries in the project, so a dataset
* only view *should* be coped with sensibly.
*/
// This must have been loaded, is it still visible?
JInternalFrame[] frames = Desktop.desktop.getAllFrames();
String matchedFile = null;
for (int f = frames.length - 1; f > -1; f--)
{
if (frames[f] instanceof StructureViewerBase)
{
StructureViewerBase viewFrame = (StructureViewerBase) frames[f];
matchedFile = saveStructureState(ap, jds, pdb, entry,
viewIds, matchedFile, viewFrame);
}
}
if (matchedFile != null || entry.getFile() != null)
{
if (entry.getFile() != null)
{
// use entry's file
matchedFile = entry.getFile();
}
pdb.setFile(matchedFile); // entry.getFile());
if (pdbfiles == null)
{
pdbfiles = new ArrayList();
}
if (!pdbfiles.contains(entry.getId()))
{
pdbfiles.add(entry.getId());
DataInputStream dis = null;
try
{
File file = new File(matchedFile);
if (file.exists() && jout != null)
{
byte[] data = new byte[(int) file.length()];
jout.putNextEntry(new JarEntry(entry.getId()));
dis = new DataInputStream(new FileInputStream(file));
dis.readFully(data);
DataOutputStream dout = new DataOutputStream(jout);
dout.write(data, 0, data.length);
dout.flush();
jout.closeEntry();
}
} catch (Exception ex)
{
ex.printStackTrace();
} finally
{
if (dis != null)
{
try
{
dis.close();
} catch (IOException e)
{
// ignore
}
}
}
}
}
if (entry.getProperty() != null)
{
PdbentryItem item = new PdbentryItem();
Hashtable properties = entry.getProperty();
Enumeration en2 = properties.keys();
while (en2.hasMoreElements())
{
Property prop = new Property();
String key = en2.nextElement().toString();
prop.setName(key);
prop.setValue(properties.get(key).toString());
item.addProperty(prop);
}
pdb.addPdbentryItem(item);
}
jseq.addPdbids(pdb);
}
}
jms.addJSeq(jseq);
}
if (!storeDS && av.hasHiddenRows())
{
jal = av.getAlignment();
}
// SAVE MAPPINGS
if (jal.getCodonFrames() != null)
{
Set jac = jal.getCodonFrames();
for (AlignedCodonFrame acf : jac)
{
AlcodonFrame alc = new AlcodonFrame();
vamsasSet.addAlcodonFrame(alc);
if (acf.getProtMappings() != null
&& acf.getProtMappings().length > 0)
{
SequenceI[] dnas = acf.getdnaSeqs();
jalview.datamodel.Mapping[] pmaps = acf.getProtMappings();
for (int m = 0; m < pmaps.length; m++)
{
AlcodMap alcmap = new AlcodMap();
alcmap.setDnasq(seqHash(dnas[m]));
alcmap.setMapping(createVamsasMapping(pmaps[m], dnas[m], null,
false));
alc.addAlcodMap(alcmap);
}
}
// {
// AlcodonFrame alc = new AlcodonFrame();
// vamsasSet.addAlcodonFrame(alc);
// for (int p = 0; p < acf.aaWidth; p++)
// {
// Alcodon cmap = new Alcodon();
// if (acf.codons[p] != null)
// {
// // Null codons indicate a gapped column in the translated peptide
// // alignment.
// cmap.setPos1(acf.codons[p][0]);
// cmap.setPos2(acf.codons[p][1]);
// cmap.setPos3(acf.codons[p][2]);
// }
// alc.addAlcodon(cmap);
// }
// if (acf.getProtMappings() != null
// && acf.getProtMappings().length > 0)
// {
// SequenceI[] dnas = acf.getdnaSeqs();
// jalview.datamodel.Mapping[] pmaps = acf.getProtMappings();
// for (int m = 0; m < pmaps.length; m++)
// {
// AlcodMap alcmap = new AlcodMap();
// alcmap.setDnasq(seqHash(dnas[m]));
// alcmap.setMapping(createVamsasMapping(pmaps[m], dnas[m], null,
// false));
// alc.addAlcodMap(alcmap);
// }
// }
}
}
// SAVE TREES
// /////////////////////////////////
if (!storeDS && av.currentTree != null)
{
// FIND ANY ASSOCIATED TREES
// NOT IMPLEMENTED FOR HEADLESS STATE AT PRESENT
if (Desktop.desktop != null)
{
JInternalFrame[] frames = Desktop.desktop.getAllFrames();
for (int t = 0; t < frames.length; t++)
{
if (frames[t] instanceof TreePanel)
{
TreePanel tp = (TreePanel) frames[t];
if (tp.treeCanvas.av.getAlignment() == jal)
{
Tree tree = new Tree();
tree.setTitle(tp.getTitle());
tree.setCurrentTree((av.currentTree == tp.getTree()));
tree.setNewick(tp.getTree().toString());
tree.setThreshold(tp.treeCanvas.threshold);
tree.setFitToWindow(tp.fitToWindow.getState());
tree.setFontName(tp.getTreeFont().getName());
tree.setFontSize(tp.getTreeFont().getSize());
tree.setFontStyle(tp.getTreeFont().getStyle());
tree.setMarkUnlinked(tp.placeholdersMenu.getState());
tree.setShowBootstrap(tp.bootstrapMenu.getState());
tree.setShowDistances(tp.distanceMenu.getState());
tree.setHeight(tp.getHeight());
tree.setWidth(tp.getWidth());
tree.setXpos(tp.getX());
tree.setYpos(tp.getY());
tree.setId(makeHashCode(tp, null));
jms.addTree(tree);
}
}
}
}
}
// SAVE ANNOTATIONS
/**
* store forward refs from an annotationRow to any groups
*/
IdentityHashMap groupRefs = new IdentityHashMap();
if (storeDS)
{
for (SequenceI sq : jal.getSequences())
{
// Store annotation on dataset sequences only
jalview.datamodel.AlignmentAnnotation[] aa = sq.getAnnotation();
if (aa != null && aa.length > 0)
{
storeAlignmentAnnotation(aa, groupRefs, av, calcIdSet, storeDS,
vamsasSet);
}
}
}
else
{
if (jal.getAlignmentAnnotation() != null)
{
// Store the annotation shown on the alignment.
jalview.datamodel.AlignmentAnnotation[] aa = jal
.getAlignmentAnnotation();
storeAlignmentAnnotation(aa, groupRefs, av, calcIdSet, storeDS,
vamsasSet);
}
}
// SAVE GROUPS
if (jal.getGroups() != null)
{
JGroup[] groups = new JGroup[jal.getGroups().size()];
int i = -1;
for (jalview.datamodel.SequenceGroup sg : jal.getGroups())
{
groups[++i] = new JGroup();
groups[i].setStart(sg.getStartRes());
groups[i].setEnd(sg.getEndRes());
groups[i].setName(sg.getName());
if (groupRefs.containsKey(sg))
{
// group has references so set it's ID field
groups[i].setId(groupRefs.get(sg).toString());
}
if (sg.cs != null)
{
if (sg.cs.conservationApplied())
{
groups[i].setConsThreshold(sg.cs.getConservationInc());
if (sg.cs instanceof jalview.schemes.UserColourScheme)
{
groups[i].setColour(setUserColourScheme(sg.cs, userColours,
jms));
}
else
{
groups[i]
.setColour(ColourSchemeProperty.getColourName(sg.cs));
}
}
else if (sg.cs instanceof jalview.schemes.AnnotationColourGradient)
{
groups[i].setColour("AnnotationColourGradient");
groups[i].setAnnotationColours(constructAnnotationColours(
(jalview.schemes.AnnotationColourGradient) sg.cs,
userColours, jms));
}
else if (sg.cs instanceof jalview.schemes.UserColourScheme)
{
groups[i]
.setColour(setUserColourScheme(sg.cs, userColours, jms));
}
else
{
groups[i].setColour(ColourSchemeProperty.getColourName(sg.cs));
}
groups[i].setPidThreshold(sg.cs.getThreshold());
}
groups[i].setOutlineColour(sg.getOutlineColour().getRGB());
groups[i].setDisplayBoxes(sg.getDisplayBoxes());
groups[i].setDisplayText(sg.getDisplayText());
groups[i].setColourText(sg.getColourText());
groups[i].setTextCol1(sg.textColour.getRGB());
groups[i].setTextCol2(sg.textColour2.getRGB());
groups[i].setTextColThreshold(sg.thresholdTextColour);
groups[i].setShowUnconserved(sg.getShowNonconserved());
groups[i].setIgnoreGapsinConsensus(sg.getIgnoreGapsConsensus());
groups[i].setShowConsensusHistogram(sg.isShowConsensusHistogram());
groups[i].setShowSequenceLogo(sg.isShowSequenceLogo());
groups[i].setNormaliseSequenceLogo(sg.isNormaliseSequenceLogo());
for (int s = 0; s < sg.getSize(); s++)
{
jalview.datamodel.Sequence seq = (jalview.datamodel.Sequence) sg
.getSequenceAt(s);
groups[i].addSeq(seqHash(seq));
}
}
jms.setJGroup(groups);
}
if (!storeDS)
{
// /////////SAVE VIEWPORT
Viewport view = new Viewport();
view.setTitle(ap.alignFrame.getTitle());
view.setSequenceSetId(makeHashCode(av.getSequenceSetId(),
av.getSequenceSetId()));
view.setId(av.getViewId());
if (av.getCodingComplement() != null)
{
view.setComplementId(av.getCodingComplement().getViewId());
}
view.setViewName(av.viewName);
view.setGatheredViews(av.isGatherViewsHere());
Rectangle position = ap.av.getExplodedGeometry();
if (position == null)
{
position = ap.alignFrame.getBounds();
}
view.setXpos(position.x);
view.setYpos(position.y);
view.setWidth(position.width);
view.setHeight(position.height);
view.setStartRes(av.startRes);
view.setStartSeq(av.startSeq);
if (av.getGlobalColourScheme() instanceof jalview.schemes.UserColourScheme)
{
view.setBgColour(setUserColourScheme(av.getGlobalColourScheme(),
userColours, jms));
}
else if (av.getGlobalColourScheme() instanceof jalview.schemes.AnnotationColourGradient)
{
AnnotationColours ac = constructAnnotationColours(
(jalview.schemes.AnnotationColourGradient) av
.getGlobalColourScheme(),
userColours, jms);
view.setAnnotationColours(ac);
view.setBgColour("AnnotationColourGradient");
}
else
{
view.setBgColour(ColourSchemeProperty.getColourName(av
.getGlobalColourScheme()));
}
ColourSchemeI cs = av.getGlobalColourScheme();
if (cs != null)
{
if (cs.conservationApplied())
{
view.setConsThreshold(cs.getConservationInc());
if (cs instanceof jalview.schemes.UserColourScheme)
{
view.setBgColour(setUserColourScheme(cs, userColours, jms));
}
}
if (cs instanceof ResidueColourScheme)
{
view.setPidThreshold(cs.getThreshold());
}
}
view.setConservationSelected(av.getConservationSelected());
view.setPidSelected(av.getAbovePIDThreshold());
view.setFontName(av.font.getName());
view.setFontSize(av.font.getSize());
view.setFontStyle(av.font.getStyle());
view.setRenderGaps(av.isRenderGaps());
view.setShowAnnotation(av.isShowAnnotation());
view.setShowBoxes(av.getShowBoxes());
view.setShowColourText(av.getColourText());
view.setShowFullId(av.getShowJVSuffix());
view.setRightAlignIds(av.isRightAlignIds());
view.setShowSequenceFeatures(av.isShowSequenceFeatures());
view.setShowText(av.getShowText());
view.setShowUnconserved(av.getShowUnconserved());
view.setWrapAlignment(av.getWrapAlignment());
view.setTextCol1(av.getTextColour().getRGB());
view.setTextCol2(av.getTextColour2().getRGB());
view.setTextColThreshold(av.getThresholdTextColour());
view.setShowConsensusHistogram(av.isShowConsensusHistogram());
view.setShowSequenceLogo(av.isShowSequenceLogo());
view.setNormaliseSequenceLogo(av.isNormaliseSequenceLogo());
view.setShowGroupConsensus(av.isShowGroupConsensus());
view.setShowGroupConservation(av.isShowGroupConservation());
view.setShowNPfeatureTooltip(av.isShowNPFeats());
view.setShowDbRefTooltip(av.isShowDBRefs());
view.setFollowHighlight(av.isFollowHighlight());
view.setFollowSelection(av.followSelection);
view.setIgnoreGapsinConsensus(av.isIgnoreGapsConsensus());
if (av.getFeaturesDisplayed() != null)
{
jalview.schemabinding.version2.FeatureSettings fs = new jalview.schemabinding.version2.FeatureSettings();
String[] renderOrder = ap.getSeqPanel().seqCanvas
.getFeatureRenderer().getRenderOrder()
.toArray(new String[0]);
Vector settingsAdded = new Vector();
Object gstyle = null;
GraduatedColor gcol = null;
if (renderOrder != null)
{
for (int ro = 0; ro < renderOrder.length; ro++)
{
gstyle = ap.getSeqPanel().seqCanvas.getFeatureRenderer()
.getFeatureStyle(renderOrder[ro]);
Setting setting = new Setting();
setting.setType(renderOrder[ro]);
if (gstyle instanceof GraduatedColor)
{
gcol = (GraduatedColor) gstyle;
setting.setColour(gcol.getMaxColor().getRGB());
setting.setMincolour(gcol.getMinColor().getRGB());
setting.setMin(gcol.getMin());
setting.setMax(gcol.getMax());
setting.setColourByLabel(gcol.isColourByLabel());
setting.setAutoScale(gcol.isAutoScale());
setting.setThreshold(gcol.getThresh());
setting.setThreshstate(gcol.getThreshType());
}
else
{
setting.setColour(ap.getSeqPanel().seqCanvas
.getFeatureRenderer()
.getColour(renderOrder[ro]).getRGB());
}
setting.setDisplay(av.getFeaturesDisplayed().isVisible(
renderOrder[ro]));
float rorder = ap.getSeqPanel().seqCanvas.getFeatureRenderer()
.getOrder(renderOrder[ro]);
if (rorder > -1)
{
setting.setOrder(rorder);
}
fs.addSetting(setting);
settingsAdded.addElement(renderOrder[ro]);
}
}
// Make sure we save none displayed feature settings
Iterator en = ap.getSeqPanel().seqCanvas.getFeatureRenderer()
.getFeatureColours().keySet().iterator();
while (en.hasNext())
{
String key = en.next().toString();
if (settingsAdded.contains(key))
{
continue;
}
Setting setting = new Setting();
setting.setType(key);
setting.setColour(ap.getSeqPanel().seqCanvas.getFeatureRenderer()
.getColour(key).getRGB());
setting.setDisplay(false);
float rorder = ap.getSeqPanel().seqCanvas.getFeatureRenderer()
.getOrder(key);
if (rorder > -1)
{
setting.setOrder(rorder);
}
fs.addSetting(setting);
settingsAdded.addElement(key);
}
// is groups actually supposed to be a map here ?
en = ap.getSeqPanel().seqCanvas.getFeatureRenderer()
.getFeatureGroups().iterator();
Vector groupsAdded = new Vector();
while (en.hasNext())
{
String grp = en.next().toString();
if (groupsAdded.contains(grp))
{
continue;
}
Group g = new Group();
g.setName(grp);
g.setDisplay(((Boolean) ap.getSeqPanel().seqCanvas
.getFeatureRenderer().checkGroupVisibility(grp, false))
.booleanValue());
fs.addGroup(g);
groupsAdded.addElement(grp);
}
jms.setFeatureSettings(fs);
}
if (av.hasHiddenColumns())
{
if (av.getColumnSelection() == null
|| av.getColumnSelection().getHiddenColumns() == null)
{
warn("REPORT BUG: avoided null columnselection bug (DMAM reported). Please contact Jim about this.");
}
else
{
for (int c = 0; c < av.getColumnSelection().getHiddenColumns()
.size(); c++)
{
int[] region = av.getColumnSelection().getHiddenColumns()
.get(c);
HiddenColumns hc = new HiddenColumns();
hc.setStart(region[0]);
hc.setEnd(region[1]);
view.addHiddenColumns(hc);
}
}
}
if (calcIdSet.size() > 0)
{
for (String calcId : calcIdSet)
{
if (calcId.trim().length() > 0)
{
CalcIdParam cidp = createCalcIdParam(calcId, av);
// Some calcIds have no parameters.
if (cidp != null)
{
view.addCalcIdParam(cidp);
}
}
}
}
jms.addViewport(view);
}
object.setJalviewModelSequence(jms);
object.getVamsasModel().addSequenceSet(vamsasSet);
if (jout != null && fileName != null)
{
// We may not want to write the object to disk,
// eg we can copy the alignViewport to a new view object
// using save and then load
try
{
JarEntry entry = new JarEntry(fileName);
jout.putNextEntry(entry);
PrintWriter pout = new PrintWriter(new OutputStreamWriter(jout,
"UTF-8"));
org.exolab.castor.xml.Marshaller marshaller = new org.exolab.castor.xml.Marshaller(
pout);
marshaller.marshal(object);
pout.flush();
jout.closeEntry();
} catch (Exception ex)
{
// TODO: raise error in GUI if marshalling failed.
ex.printStackTrace();
}
}
return object;
}
/**
* Save the state of a structure viewer
*
* @param ap
* @param jds
* @param pdb
* the archive XML element under which to save the state
* @param entry
* @param viewIds
* @param matchedFile
* @param viewFrame
* @return
*/
protected String saveStructureState(AlignmentPanel ap, SequenceI jds,
Pdbids pdb, PDBEntry entry, List viewIds,
String matchedFile, StructureViewerBase viewFrame)
{
final AAStructureBindingModel bindingModel = viewFrame.getBinding();
for (int peid = 0; peid < bindingModel.getPdbCount(); peid++)
{
final PDBEntry pdbentry = bindingModel.getPdbEntry(peid);
final String pdbId = pdbentry.getId();
if (!pdbId.equals(entry.getId())
&& !(entry.getId().length() > 4 && entry.getId()
.toLowerCase().startsWith(pdbId.toLowerCase())))
{
continue;
}
if (matchedFile == null)
{
matchedFile = pdbentry.getFile();
}
else if (!matchedFile.equals(pdbentry.getFile()))
{
Cache.log
.warn("Probably lost some PDB-Sequence mappings for this structure file (which apparently has same PDB Entry code): "
+ pdbentry.getFile());
}
// record the
// file so we
// can get at it if the ID
// match is ambiguous (e.g.
// 1QIP==1qipA)
String statestring = viewFrame.getStateInfo();
for (int smap = 0; smap < viewFrame.getBinding().getSequence()[peid].length; smap++)
{
// if (jal.findIndex(jmol.jmb.sequence[peid][smap]) > -1)
if (jds == viewFrame.getBinding().getSequence()[peid][smap])
{
StructureState state = new StructureState();
state.setVisible(true);
state.setXpos(viewFrame.getX());
state.setYpos(viewFrame.getY());
state.setWidth(viewFrame.getWidth());
state.setHeight(viewFrame.getHeight());
final String viewId = viewFrame.getViewId();
state.setViewId(viewId);
state.setAlignwithAlignPanel(viewFrame.isUsedforaligment(ap));
state.setColourwithAlignPanel(viewFrame.isUsedforcolourby(ap));
state.setColourByJmol(viewFrame.isColouredByViewer());
/*
* Only store each structure viewer's state once in each XML document.
*/
if (!viewIds.contains(viewId))
{
viewIds.add(viewId);
state.setContent(statestring.replaceAll("\n", ""));
}
else
{
state.setContent("# duplicate state");
}
pdb.addStructureState(state);
}
}
}
return matchedFile;
}
private AnnotationColours constructAnnotationColours(
AnnotationColourGradient acg, List userColours,
JalviewModelSequence jms)
{
AnnotationColours ac = new AnnotationColours();
ac.setAboveThreshold(acg.getAboveThreshold());
ac.setThreshold(acg.getAnnotationThreshold());
ac.setAnnotation(acg.getAnnotation());
if (acg.getBaseColour() instanceof jalview.schemes.UserColourScheme)
{
ac.setColourScheme(setUserColourScheme(acg.getBaseColour(),
userColours, jms));
}
else
{
ac.setColourScheme(ColourSchemeProperty.getColourName(acg
.getBaseColour()));
}
ac.setMaxColour(acg.getMaxColour().getRGB());
ac.setMinColour(acg.getMinColour().getRGB());
ac.setPerSequence(acg.isSeqAssociated());
ac.setPredefinedColours(acg.isPredefinedColours());
return ac;
}
private void storeAlignmentAnnotation(AlignmentAnnotation[] aa,
IdentityHashMap groupRefs, AlignmentViewport av,
Set calcIdSet, boolean storeDS, SequenceSet vamsasSet)
{
for (int i = 0; i < aa.length; i++)
{
Annotation an = new Annotation();
if (aa[i].annotationId != null)
{
annotationIds.put(aa[i].annotationId, aa[i]);
}
an.setId(aa[i].annotationId);
an.setVisible(aa[i].visible);
an.setDescription(aa[i].description);
if (aa[i].sequenceRef != null)
{
// TODO later annotation sequenceRef should be the XML ID of the
// sequence rather than its display name
an.setSequenceRef(aa[i].sequenceRef.getName());
}
if (aa[i].groupRef != null)
{
Object groupIdr = groupRefs.get(aa[i].groupRef);
if (groupIdr == null)
{
// make a locally unique String
groupRefs.put(aa[i].groupRef,
groupIdr = ("" + System.currentTimeMillis()
+ aa[i].groupRef.getName() + groupRefs.size()));
}
an.setGroupRef(groupIdr.toString());
}
// store all visualization attributes for annotation
an.setGraphHeight(aa[i].graphHeight);
an.setCentreColLabels(aa[i].centreColLabels);
an.setScaleColLabels(aa[i].scaleColLabel);
an.setShowAllColLabels(aa[i].showAllColLabels);
an.setBelowAlignment(aa[i].belowAlignment);
if (aa[i].graph > 0)
{
an.setGraph(true);
an.setGraphType(aa[i].graph);
an.setGraphGroup(aa[i].graphGroup);
if (aa[i].getThreshold() != null)
{
ThresholdLine line = new ThresholdLine();
line.setLabel(aa[i].getThreshold().label);
line.setValue(aa[i].getThreshold().value);
line.setColour(aa[i].getThreshold().colour.getRGB());
an.setThresholdLine(line);
}
}
else
{
an.setGraph(false);
}
an.setLabel(aa[i].label);
if (aa[i] == av.getAlignmentQualityAnnot()
|| aa[i] == av.getAlignmentConservationAnnotation()
|| aa[i] == av.getAlignmentConsensusAnnotation()
|| aa[i].autoCalculated)
{
// new way of indicating autocalculated annotation -
an.setAutoCalculated(aa[i].autoCalculated);
}
if (aa[i].hasScore())
{
an.setScore(aa[i].getScore());
}
if (aa[i].getCalcId() != null)
{
calcIdSet.add(aa[i].getCalcId());
an.setCalcId(aa[i].getCalcId());
}
if (aa[i].hasProperties())
{
for (String pr : aa[i].getProperties())
{
Property prop = new Property();
prop.setName(pr);
prop.setValue(aa[i].getProperty(pr));
an.addProperty(prop);
}
}
AnnotationElement ae;
if (aa[i].annotations != null)
{
an.setScoreOnly(false);
for (int a = 0; a < aa[i].annotations.length; a++)
{
if ((aa[i] == null) || (aa[i].annotations[a] == null))
{
continue;
}
ae = new AnnotationElement();
if (aa[i].annotations[a].description != null)
{
ae.setDescription(aa[i].annotations[a].description);
}
if (aa[i].annotations[a].displayCharacter != null)
{
ae.setDisplayCharacter(aa[i].annotations[a].displayCharacter);
}
if (!Float.isNaN(aa[i].annotations[a].value))
{
ae.setValue(aa[i].annotations[a].value);
}
ae.setPosition(a);
if (aa[i].annotations[a].secondaryStructure > ' ')
{
ae.setSecondaryStructure(aa[i].annotations[a].secondaryStructure
+ "");
}
if (aa[i].annotations[a].colour != null
&& aa[i].annotations[a].colour != java.awt.Color.black)
{
ae.setColour(aa[i].annotations[a].colour.getRGB());
}
an.addAnnotationElement(ae);
if (aa[i].autoCalculated)
{
// only write one non-null entry into the annotation row -
// sufficient to get the visualization attributes necessary to
// display data
continue;
}
}
}
else
{
an.setScoreOnly(true);
}
if (!storeDS || (storeDS && !aa[i].autoCalculated))
{
// skip autocalculated annotation - these are only provided for
// alignments
vamsasSet.addAnnotation(an);
}
}
}
private CalcIdParam createCalcIdParam(String calcId, AlignViewport av)
{
AutoCalcSetting settings = av.getCalcIdSettingsFor(calcId);
if (settings != null)
{
CalcIdParam vCalcIdParam = new CalcIdParam();
vCalcIdParam.setCalcId(calcId);
vCalcIdParam.addServiceURL(settings.getServiceURI());
// generic URI allowing a third party to resolve another instance of the
// service used for this calculation
for (String urls : settings.getServiceURLs())
{
vCalcIdParam.addServiceURL(urls);
}
vCalcIdParam.setVersion("1.0");
if (settings.getPreset() != null)
{
WsParamSetI setting = settings.getPreset();
vCalcIdParam.setName(setting.getName());
vCalcIdParam.setDescription(setting.getDescription());
}
else
{
vCalcIdParam.setName("");
vCalcIdParam.setDescription("Last used parameters");
}
// need to be able to recover 1) settings 2) user-defined presets or
// recreate settings from preset 3) predefined settings provided by
// service - or settings that can be transferred (or discarded)
vCalcIdParam.setParameters(settings.getWsParamFile().replace("\n",
"|\\n|"));
vCalcIdParam.setAutoUpdate(settings.isAutoUpdate());
// todo - decide if updateImmediately is needed for any projects.
return vCalcIdParam;
}
return null;
}
private boolean recoverCalcIdParam(CalcIdParam calcIdParam,
AlignViewport av)
{
if (calcIdParam.getVersion().equals("1.0"))
{
Jws2Instance service = Jws2Discoverer.getDiscoverer()
.getPreferredServiceFor(calcIdParam.getServiceURL());
if (service != null)
{
WsParamSetI parmSet = null;
try
{
parmSet = service.getParamStore().parseServiceParameterFile(
calcIdParam.getName(), calcIdParam.getDescription(),
calcIdParam.getServiceURL(),
calcIdParam.getParameters().replace("|\\n|", "\n"));
} catch (IOException x)
{
warn("Couldn't parse parameter data for "
+ calcIdParam.getCalcId(), x);
return false;
}
List argList = null;
if (calcIdParam.getName().length() > 0)
{
parmSet = service.getParamStore()
.getPreset(calcIdParam.getName());
if (parmSet != null)
{
// TODO : check we have a good match with settings in AACon -
// otherwise we'll need to create a new preset
}
}
else
{
argList = parmSet.getArguments();
parmSet = null;
}
AAConSettings settings = new AAConSettings(
calcIdParam.isAutoUpdate(), service, parmSet, argList);
av.setCalcIdSettingsFor(calcIdParam.getCalcId(), settings,
calcIdParam.isNeedsUpdate());
return true;
}
else
{
warn("Cannot resolve a service for the parameters used in this project. Try configuring a JABAWS server.");
return false;
}
}
throw new Error(MessageManager.formatMessage(
"error.unsupported_version_calcIdparam", new Object[]
{ calcIdParam.toString() }));
}
/**
* External mapping between jalview objects and objects yielding a valid and
* unique object ID string. This is null for normal Jalview project IO, but
* non-null when a jalview project is being read or written as part of a
* vamsas session.
*/
IdentityHashMap jv2vobj = null;
/**
* Construct a unique ID for jvobj using either existing bindings or if none
* exist, the result of the hashcode call for the object.
*
* @param jvobj
* jalview data object
* @return unique ID for referring to jvobj
*/
private String makeHashCode(Object jvobj, String altCode)
{
if (jv2vobj != null)
{
Object id = jv2vobj.get(jvobj);
if (id != null)
{
return id.toString();
}
// check string ID mappings
if (jvids2vobj != null && jvobj instanceof String)
{
id = jvids2vobj.get(jvobj);
}
if (id != null)
{
return id.toString();
}
// give up and warn that something has gone wrong
warn("Cannot find ID for object in external mapping : " + jvobj);
}
return altCode;
}
/**
* return local jalview object mapped to ID, if it exists
*
* @param idcode
* (may be null)
* @return null or object bound to idcode
*/
private Object retrieveExistingObj(String idcode)
{
if (idcode != null && vobj2jv != null)
{
return vobj2jv.get(idcode);
}
return null;
}
/**
* binding from ID strings from external mapping table to jalview data model
* objects.
*/
private Hashtable vobj2jv;
private Sequence createVamsasSequence(String id, SequenceI jds)
{
return createVamsasSequence(true, id, jds, null);
}
private Sequence createVamsasSequence(boolean recurse, String id,
SequenceI jds, SequenceI parentseq)
{
Sequence vamsasSeq = new Sequence();
vamsasSeq.setId(id);
vamsasSeq.setName(jds.getName());
vamsasSeq.setSequence(jds.getSequenceAsString());
vamsasSeq.setDescription(jds.getDescription());
jalview.datamodel.DBRefEntry[] dbrefs = null;
if (jds.getDatasetSequence() != null)
{
vamsasSeq.setDsseqid(seqHash(jds.getDatasetSequence()));
if (jds.getDatasetSequence().getDBRef() != null)
{
dbrefs = jds.getDatasetSequence().getDBRef();
}
}
else
{
vamsasSeq.setDsseqid(id); // so we can tell which sequences really are
// dataset sequences only
dbrefs = jds.getDBRef();
}
if (dbrefs != null)
{
for (int d = 0; d < dbrefs.length; d++)
{
DBRef dbref = new DBRef();
dbref.setSource(dbrefs[d].getSource());
dbref.setVersion(dbrefs[d].getVersion());
dbref.setAccessionId(dbrefs[d].getAccessionId());
if (dbrefs[d].hasMap())
{
Mapping mp = createVamsasMapping(dbrefs[d].getMap(), parentseq,
jds, recurse);
dbref.setMapping(mp);
}
vamsasSeq.addDBRef(dbref);
}
}
return vamsasSeq;
}
private Mapping createVamsasMapping(jalview.datamodel.Mapping jmp,
SequenceI parentseq, SequenceI jds, boolean recurse)
{
Mapping mp = null;
if (jmp.getMap() != null)
{
mp = new Mapping();
jalview.util.MapList mlst = jmp.getMap();
List r = mlst.getFromRanges();
for (int[] range : r)
{
MapListFrom mfrom = new MapListFrom();
mfrom.setStart(range[0]);
mfrom.setEnd(range[1]);
mp.addMapListFrom(mfrom);
}
r = mlst.getToRanges();
for (int[] range : r)
{
MapListTo mto = new MapListTo();
mto.setStart(range[0]);
mto.setEnd(range[1]);
mp.addMapListTo(mto);
}
mp.setMapFromUnit(mlst.getFromRatio());
mp.setMapToUnit(mlst.getToRatio());
if (jmp.getTo() != null)
{
MappingChoice mpc = new MappingChoice();
if (recurse
&& (parentseq != jmp.getTo() || parentseq
.getDatasetSequence() != jmp.getTo()))
{
mpc.setSequence(createVamsasSequence(false, seqHash(jmp.getTo()),
jmp.getTo(), jds));
}
else
{
String jmpid = "";
SequenceI ps = null;
if (parentseq != jmp.getTo()
&& parentseq.getDatasetSequence() != jmp.getTo())
{
// chaining dbref rather than a handshaking one
jmpid = seqHash(ps = jmp.getTo());
}
else
{
jmpid = seqHash(ps = parentseq);
}
mpc.setDseqFor(jmpid);
if (!seqRefIds.containsKey(mpc.getDseqFor()))
{
jalview.bin.Cache.log.debug("creatign new DseqFor ID");
seqRefIds.put(mpc.getDseqFor(), ps);
}
else
{
jalview.bin.Cache.log.debug("reusing DseqFor ID");
}
}
mp.setMappingChoice(mpc);
}
}
return mp;
}
String setUserColourScheme(jalview.schemes.ColourSchemeI cs,
List userColours, JalviewModelSequence jms)
{
String id = null;
jalview.schemes.UserColourScheme ucs = (jalview.schemes.UserColourScheme) cs;
boolean newucs = false;
if (!userColours.contains(ucs))
{
userColours.add(ucs);
newucs = true;
}
id = "ucs" + userColours.indexOf(ucs);
if (newucs)
{
// actually create the scheme's entry in the XML model
java.awt.Color[] colours = ucs.getColours();
jalview.schemabinding.version2.UserColours uc = new jalview.schemabinding.version2.UserColours();
jalview.schemabinding.version2.UserColourScheme jbucs = new jalview.schemabinding.version2.UserColourScheme();
for (int i = 0; i < colours.length; i++)
{
jalview.schemabinding.version2.Colour col = new jalview.schemabinding.version2.Colour();
col.setName(ResidueProperties.aa[i]);
col.setRGB(jalview.util.Format.getHexString(colours[i]));
jbucs.addColour(col);
}
if (ucs.getLowerCaseColours() != null)
{
colours = ucs.getLowerCaseColours();
for (int i = 0; i < colours.length; i++)
{
jalview.schemabinding.version2.Colour col = new jalview.schemabinding.version2.Colour();
col.setName(ResidueProperties.aa[i].toLowerCase());
col.setRGB(jalview.util.Format.getHexString(colours[i]));
jbucs.addColour(col);
}
}
uc.setId(id);
uc.setUserColourScheme(jbucs);
jms.addUserColours(uc);
}
return id;
}
jalview.schemes.UserColourScheme getUserColourScheme(
JalviewModelSequence jms, String id)
{
UserColours[] uc = jms.getUserColours();
UserColours colours = null;
for (int i = 0; i < uc.length; i++)
{
if (uc[i].getId().equals(id))
{
colours = uc[i];
break;
}
}
java.awt.Color[] newColours = new java.awt.Color[24];
for (int i = 0; i < 24; i++)
{
newColours[i] = new java.awt.Color(Integer.parseInt(colours
.getUserColourScheme().getColour(i).getRGB(), 16));
}
jalview.schemes.UserColourScheme ucs = new jalview.schemes.UserColourScheme(
newColours);
if (colours.getUserColourScheme().getColourCount() > 24)
{
newColours = new java.awt.Color[23];
for (int i = 0; i < 23; i++)
{
newColours[i] = new java.awt.Color(Integer.parseInt(colours
.getUserColourScheme().getColour(i + 24).getRGB(), 16));
}
ucs.setLowerCaseColours(newColours);
}
return ucs;
}
/**
* contains last error message (if any) encountered by XML loader.
*/
String errorMessage = null;
/**
* flag to control whether the Jalview2XML_V1 parser should be deferred to if
* exceptions are raised during project XML parsing
*/
public boolean attemptversion1parse = true;
/**
* Load a jalview project archive from a jar file
*
* @param file
* - HTTP URL or filename
*/
public AlignFrame loadJalviewAlign(final String file)
{
jalview.gui.AlignFrame af = null;
try
{
// create list to store references for any new Jmol viewers created
newStructureViewers = new Vector();
// UNMARSHALLER SEEMS TO CLOSE JARINPUTSTREAM, MOST ANNOYING
// Workaround is to make sure caller implements the JarInputStreamProvider
// interface
// so we can re-open the jar input stream for each entry.
jarInputStreamProvider jprovider = createjarInputStreamProvider(file);
af = loadJalviewAlign(jprovider);
} catch (MalformedURLException e)
{
errorMessage = "Invalid URL format for '" + file + "'";
reportErrors();
} finally
{
try
{
SwingUtilities.invokeAndWait(new Runnable()
{
public void run()
{
setLoadingFinishedForNewStructureViewers();
};
});
} catch (Exception x)
{
}
}
return af;
}
private jarInputStreamProvider createjarInputStreamProvider(
final String file) throws MalformedURLException
{
URL url = null;
errorMessage = null;
uniqueSetSuffix = null;
seqRefIds = null;
viewportsAdded.clear();
frefedSequence = null;
if (file.startsWith("http://"))
{
url = new URL(file);
}
final URL _url = url;
return new jarInputStreamProvider()
{
@Override
public JarInputStream getJarInputStream() throws IOException
{
if (_url != null)
{
return new JarInputStream(_url.openStream());
}
else
{
return new JarInputStream(new FileInputStream(file));
}
}
@Override
public String getFilename()
{
return file;
}
};
}
/**
* Recover jalview session from a jalview project archive. Caller may
* initialise uniqueSetSuffix, seqRefIds, viewportsAdded and frefedSequence
* themselves. Any null fields will be initialised with default values,
* non-null fields are left alone.
*
* @param jprovider
* @return
*/
public AlignFrame loadJalviewAlign(final jarInputStreamProvider jprovider)
{
errorMessage = null;
if (uniqueSetSuffix == null)
{
uniqueSetSuffix = System.currentTimeMillis() % 100000 + "";
}
if (seqRefIds == null)
{
seqRefIds = new HashMap();
}
if (frefedSequence == null)
{
frefedSequence = new Vector();
}
AlignFrame af = null, _af = null;
Map gatherToThisFrame = new HashMap();
final String file = jprovider.getFilename();
try
{
JarInputStream jin = null;
JarEntry jarentry = null;
int entryCount = 1;
do
{
jin = jprovider.getJarInputStream();
for (int i = 0; i < entryCount; i++)
{
jarentry = jin.getNextJarEntry();
}
if (jarentry != null && jarentry.getName().endsWith(".xml"))
{
InputStreamReader in = new InputStreamReader(jin, "UTF-8");
JalviewModel object = new JalviewModel();
Unmarshaller unmar = new Unmarshaller(object);
unmar.setValidation(false);
object = (JalviewModel) unmar.unmarshal(in);
if (true) // !skipViewport(object))
{
_af = loadFromObject(object, file, true, jprovider);
if (object.getJalviewModelSequence().getViewportCount() > 0)
{
af = _af;
if (af.viewport.isGatherViewsHere())
{
gatherToThisFrame.put(af.viewport.getSequenceSetId(), af);
}
}
}
entryCount++;
}
else if (jarentry != null)
{
// Some other file here.
entryCount++;
}
} while (jarentry != null);
resolveFrefedSequences();
} catch (java.io.FileNotFoundException ex)
{
ex.printStackTrace();
errorMessage = "Couldn't locate Jalview XML file : " + file;
System.err.println("Exception whilst loading jalview XML file : "
+ ex + "\n");
} catch (java.net.UnknownHostException ex)
{
ex.printStackTrace();
errorMessage = "Couldn't locate Jalview XML file : " + file;
System.err.println("Exception whilst loading jalview XML file : "
+ ex + "\n");
} catch (Exception ex)
{
System.err.println("Parsing as Jalview Version 2 file failed.");
ex.printStackTrace(System.err);
if (attemptversion1parse)
{
// Is Version 1 Jar file?
try
{
af = new Jalview2XML_V1(raiseGUI).LoadJalviewAlign(jprovider);
} catch (Exception ex2)
{
System.err.println("Exception whilst loading as jalviewXMLV1:");
ex2.printStackTrace();
af = null;
}
}
if (Desktop.instance != null)
{
Desktop.instance.stopLoading();
}
if (af != null)
{
System.out.println("Successfully loaded archive file");
return af;
}
ex.printStackTrace();
System.err.println("Exception whilst loading jalview XML file : "
+ ex + "\n");
} catch (OutOfMemoryError e)
{
// Don't use the OOM Window here
errorMessage = "Out of memory loading jalview XML file";
System.err.println("Out of memory whilst loading jalview XML file");
e.printStackTrace();
}
if (Desktop.instance != null)
{
Desktop.instance.stopLoading();
}
/*
* Regather multiple views (with the same sequence set id) to the frame (if
* any) that is flagged as the one to gather to, i.e. convert them to tabbed
* views instead of separate frames. Note this doesn't restore a state where
* some expanded views in turn have tabbed views - the last "first tab" read
* in will play the role of gatherer for all.
*/
for (AlignFrame fr : gatherToThisFrame.values())
{
Desktop.instance.gatherViews(fr);
}
restoreSplitFrames();
if (errorMessage != null)
{
reportErrors();
}
return af;
}
/**
* Try to reconstruct and display SplitFrame windows, where each contains
* complementary dna and protein alignments. Done by pairing up AlignFrame
* objects (created earlier) which have complementary viewport ids associated.
*/
protected void restoreSplitFrames()
{
List gatherTo = new ArrayList();
List addedToSplitFrames = new ArrayList();
Map dna = new HashMap();
/*
* Identify the DNA alignments
*/
for (Entry candidate : splitFrameCandidates
.entrySet())
{
AlignFrame af = candidate.getValue();
if (af.getViewport().getAlignment().isNucleotide())
{
dna.put(candidate.getKey().getId(), af);
}
}
/*
* Try to match up the protein complements
*/
for (Entry candidate : splitFrameCandidates
.entrySet())
{
AlignFrame af = candidate.getValue();
if (!af.getViewport().getAlignment().isNucleotide())
{
String complementId = candidate.getKey().getComplementId();
// only non-null complements should be in the Map
if (complementId != null && dna.containsKey(complementId))
{
final AlignFrame dnaFrame = dna.get(complementId);
SplitFrame sf = createSplitFrame(dnaFrame, af);
addedToSplitFrames.add(dnaFrame);
addedToSplitFrames.add(af);
if (af.viewport.isGatherViewsHere())
{
gatherTo.add(sf);
}
}
}
}
/*
* Open any that we failed to pair up (which shouldn't happen!) as
* standalone AlignFrame's.
*/
for (Entry candidate : splitFrameCandidates
.entrySet())
{
AlignFrame af = candidate.getValue();
if (!addedToSplitFrames.contains(af)) {
Viewport view = candidate.getKey();
Desktop.addInternalFrame(af, view.getTitle(), view.getWidth(),
view.getHeight());
System.err.println("Failed to restore view " + view.getTitle()
+ " to split frame");
}
}
/*
* Gather back into tabbed views as flagged.
*/
for (SplitFrame sf : gatherTo)
{
Desktop.instance.gatherViews(sf);
}
splitFrameCandidates.clear();
}
/**
* Construct and display one SplitFrame holding DNA and protein alignments.
*
* @param dnaFrame
* @param proteinFrame
* @return
*/
protected SplitFrame createSplitFrame(AlignFrame dnaFrame,
AlignFrame proteinFrame)
{
SplitFrame splitFrame = new SplitFrame(dnaFrame, proteinFrame);
String title = MessageManager.getString("label.linked_view_title");
int width = (int) dnaFrame.getBounds().getWidth();
int height = (int) (dnaFrame.getBounds().getHeight()
+ proteinFrame.getBounds().getHeight() + 50);
Desktop.addInternalFrame(splitFrame, title, width, height);
/*
* And compute cDNA consensus (couldn't do earlier with consensus as
* mappings were not yet present)
*/
proteinFrame.viewport.alignmentChanged(proteinFrame.alignPanel);
return splitFrame;
}
/**
* check errorMessage for a valid error message and raise an error box in the
* GUI or write the current errorMessage to stderr and then clear the error
* state.
*/
protected void reportErrors()
{
reportErrors(false);
}
protected void reportErrors(final boolean saving)
{
if (errorMessage != null)
{
final String finalErrorMessage = errorMessage;
if (raiseGUI)
{
javax.swing.SwingUtilities.invokeLater(new Runnable()
{
@Override
public void run()
{
JOptionPane.showInternalMessageDialog(Desktop.desktop,
finalErrorMessage, "Error "
+ (saving ? "saving" : "loading")
+ " Jalview file", JOptionPane.WARNING_MESSAGE);
}
});
}
else
{
System.err.println("Problem loading Jalview file: " + errorMessage);
}
}
errorMessage = null;
}
Map alreadyLoadedPDB = new HashMap();
/**
* when set, local views will be updated from view stored in JalviewXML
* Currently (28th Sep 2008) things will go horribly wrong in vamsas document
* sync if this is set to true.
*/
private final boolean updateLocalViews = false;
String loadPDBFile(jarInputStreamProvider jprovider, String pdbId)
{
if (alreadyLoadedPDB.containsKey(pdbId))
{
return alreadyLoadedPDB.get(pdbId).toString();
}
try
{
JarInputStream jin = jprovider.getJarInputStream();
/*
* if (jprovider.startsWith("http://")) { jin = new JarInputStream(new
* URL(jprovider).openStream()); } else { jin = new JarInputStream(new
* FileInputStream(jprovider)); }
*/
JarEntry entry = null;
do
{
entry = jin.getNextJarEntry();
} while (entry != null && !entry.getName().equals(pdbId));
if (entry != null)
{
BufferedReader in = new BufferedReader(new InputStreamReader(jin));
File outFile = File.createTempFile("jalview_pdb", ".txt");
outFile.deleteOnExit();
PrintWriter out = new PrintWriter(new FileOutputStream(outFile));
String data;
while ((data = in.readLine()) != null)
{
out.println(data);
}
try
{
out.flush();
} catch (Exception foo)
{
}
;
out.close();
String t = outFile.getAbsolutePath();
alreadyLoadedPDB.put(pdbId, t);
return t;
}
else
{
warn("Couldn't find PDB file entry in Jalview Jar for " + pdbId);
}
} catch (Exception ex)
{
ex.printStackTrace();
}
return null;
}
private class JvAnnotRow
{
public JvAnnotRow(int i, AlignmentAnnotation jaa)
{
order = i;
template = jaa;
}
/**
* persisted version of annotation row from which to take vis properties
*/
public jalview.datamodel.AlignmentAnnotation template;
/**
* original position of the annotation row in the alignment
*/
public int order;
}
/**
* Load alignment frame from jalview XML DOM object
*
* @param object
* DOM
* @param file
* filename source string
* @param loadTreesAndStructures
* when false only create Viewport
* @param jprovider
* data source provider
* @return alignment frame created from view stored in DOM
*/
AlignFrame loadFromObject(JalviewModel object, String file,
boolean loadTreesAndStructures, jarInputStreamProvider jprovider)
{
SequenceSet vamsasSet = object.getVamsasModel().getSequenceSet(0);
Sequence[] vamsasSeq = vamsasSet.getSequence();
JalviewModelSequence jms = object.getJalviewModelSequence();
Viewport view = (jms.getViewportCount() > 0) ? jms.getViewport(0)
: null;
// ////////////////////////////////
// LOAD SEQUENCES
List hiddenSeqs = null;
jalview.datamodel.Sequence jseq;
List tmpseqs = new ArrayList();
boolean multipleView = false;
JSeq[] jseqs = object.getJalviewModelSequence().getJSeq();
int vi = 0; // counter in vamsasSeq array
for (int i = 0; i < jseqs.length; i++)
{
String seqId = jseqs[i].getId();
if (seqRefIds.get(seqId) != null)
{
tmpseqs.add(seqRefIds.get(seqId));
multipleView = true;
}
else
{
jseq = new jalview.datamodel.Sequence(vamsasSeq[vi].getName(),
vamsasSeq[vi].getSequence());
jseq.setDescription(vamsasSeq[vi].getDescription());
jseq.setStart(jseqs[i].getStart());
jseq.setEnd(jseqs[i].getEnd());
jseq.setVamsasId(uniqueSetSuffix + seqId);
seqRefIds.put(vamsasSeq[vi].getId(), jseq);
tmpseqs.add(jseq);
vi++;
}
if (jseqs[i].getHidden())
{
if (hiddenSeqs == null)
{
hiddenSeqs = new ArrayList();
}
hiddenSeqs.add(seqRefIds.get(seqId));
}
}
// /
// Create the alignment object from the sequence set
// ///////////////////////////////
SequenceI[] orderedSeqs = tmpseqs
.toArray(new SequenceI[tmpseqs.size()]);
Alignment al = new Alignment(orderedSeqs);
// / Add the alignment properties
for (int i = 0; i < vamsasSet.getSequenceSetPropertiesCount(); i++)
{
SequenceSetProperties ssp = vamsasSet.getSequenceSetProperties(i);
al.setProperty(ssp.getKey(), ssp.getValue());
}
// /
// SequenceFeatures are added to the DatasetSequence,
// so we must create or recover the dataset before loading features
// ///////////////////////////////
if (vamsasSet.getDatasetId() == null || vamsasSet.getDatasetId() == "")
{
// older jalview projects do not have a dataset id.
al.setDataset(null);
}
else
{
// recover dataset - passing on flag indicating if this a 'viewless'
// sequence set (a.k.a. a stored dataset for the project)
recoverDatasetFor(vamsasSet, al, object.getJalviewModelSequence()
.getViewportCount() == 0);
}
// ///////////////////////////////
Hashtable pdbloaded = new Hashtable(); // TODO nothing writes to this??
if (!multipleView)
{
// load sequence features, database references and any associated PDB
// structures for the alignment
for (int i = 0; i < vamsasSeq.length; i++)
{
if (jseqs[i].getFeaturesCount() > 0)
{
Features[] features = jseqs[i].getFeatures();
for (int f = 0; f < features.length; f++)
{
jalview.datamodel.SequenceFeature sf = new jalview.datamodel.SequenceFeature(
features[f].getType(), features[f].getDescription(),
features[f].getStatus(), features[f].getBegin(),
features[f].getEnd(), features[f].getFeatureGroup());
sf.setScore(features[f].getScore());
for (int od = 0; od < features[f].getOtherDataCount(); od++)
{
OtherData keyValue = features[f].getOtherData(od);
if (keyValue.getKey().startsWith("LINK"))
{
sf.addLink(keyValue.getValue());
}
else
{
sf.setValue(keyValue.getKey(), keyValue.getValue());
}
}
al.getSequenceAt(i).getDatasetSequence().addSequenceFeature(sf);
}
}
if (vamsasSeq[i].getDBRefCount() > 0)
{
addDBRefs(al.getSequenceAt(i).getDatasetSequence(), vamsasSeq[i]);
}
if (jseqs[i].getPdbidsCount() > 0)
{
Pdbids[] ids = jseqs[i].getPdbids();
for (int p = 0; p < ids.length; p++)
{
jalview.datamodel.PDBEntry entry = new jalview.datamodel.PDBEntry();
entry.setId(ids[p].getId());
if (ids[p].getType() != null)
{
if (ids[p].getType().equalsIgnoreCase("PDB"))
{
entry.setType(PDBEntry.Type.PDB);
}
else
{
entry.setType(PDBEntry.Type.FILE);
}
}
if (ids[p].getFile() != null)
{
if (!pdbloaded.containsKey(ids[p].getFile()))
{
entry.setFile(loadPDBFile(jprovider, ids[p].getId()));
}
else
{
entry.setFile(pdbloaded.get(ids[p].getId()).toString());
}
}
StructureSelectionManager.getStructureSelectionManager(
Desktop.instance).registerPDBEntry(entry);
al.getSequenceAt(i).getDatasetSequence().addPDBId(entry);
}
}
}
} // end !multipleview
// ///////////////////////////////
// LOAD SEQUENCE MAPPINGS
if (vamsasSet.getAlcodonFrameCount() > 0)
{
// TODO Potentially this should only be done once for all views of an
// alignment
AlcodonFrame[] alc = vamsasSet.getAlcodonFrame();
for (int i = 0; i < alc.length; i++)
{
AlignedCodonFrame cf = new AlignedCodonFrame();
if (alc[i].getAlcodMapCount() > 0)
{
AlcodMap[] maps = alc[i].getAlcodMap();
for (int m = 0; m < maps.length; m++)
{
SequenceI dnaseq = seqRefIds.get(maps[m].getDnasq());
// Load Mapping
jalview.datamodel.Mapping mapping = null;
// attach to dna sequence reference.
if (maps[m].getMapping() != null)
{
mapping = addMapping(maps[m].getMapping());
}
if (dnaseq != null)
{
cf.addMap(dnaseq, mapping.getTo(), mapping.getMap());
}
else
{
// defer to later
frefedSequence.add(new Object[]
{ maps[m].getDnasq(), cf, mapping });
}
}
}
al.addCodonFrame(cf);
}
}
// ////////////////////////////////
// LOAD ANNOTATIONS
List autoAlan = new ArrayList();
/**
* store any annotations which forward reference a group's ID
*/
Hashtable> groupAnnotRefs = new Hashtable>();
if (vamsasSet.getAnnotationCount() > 0)
{
Annotation[] an = vamsasSet.getAnnotation();
for (int i = 0; i < an.length; i++)
{
/**
* test if annotation is automatically calculated for this view only
*/
boolean autoForView = false;
if (an[i].getLabel().equals("Quality")
|| an[i].getLabel().equals("Conservation")
|| an[i].getLabel().equals("Consensus"))
{
// Kludge for pre 2.5 projects which lacked the autocalculated flag
autoForView = true;
if (!an[i].hasAutoCalculated())
{
an[i].setAutoCalculated(true);
}
}
if (autoForView
|| (an[i].hasAutoCalculated() && an[i].isAutoCalculated()))
{
// remove ID - we don't recover annotation from other views for
// view-specific annotation
an[i].setId(null);
}
// set visiblity for other annotation in this view
if (an[i].getId() != null
&& annotationIds.containsKey(an[i].getId()))
{
AlignmentAnnotation jda = annotationIds.get(an[i].getId());
// in principle Visible should always be true for annotation displayed
// in multiple views
if (an[i].hasVisible())
{
jda.visible = an[i].getVisible();
}
al.addAnnotation(jda);
continue;
}
// Construct new annotation from model.
AnnotationElement[] ae = an[i].getAnnotationElement();
jalview.datamodel.Annotation[] anot = null;
java.awt.Color firstColour = null;
int anpos;
if (!an[i].getScoreOnly())
{
anot = new jalview.datamodel.Annotation[al.getWidth()];
for (int aa = 0; aa < ae.length && aa < anot.length; aa++)
{
anpos = ae[aa].getPosition();
if (anpos >= anot.length)
{
continue;
}
anot[anpos] = new jalview.datamodel.Annotation(
ae[aa].getDisplayCharacter(), ae[aa].getDescription(),
(ae[aa].getSecondaryStructure() == null || ae[aa]
.getSecondaryStructure().length() == 0) ? ' '
: ae[aa].getSecondaryStructure().charAt(0),
ae[aa].getValue()
);
// JBPNote: Consider verifying dataflow for IO of secondary
// structure annotation read from Stockholm files
// this was added to try to ensure that
// if (anot[ae[aa].getPosition()].secondaryStructure>' ')
// {
// anot[ae[aa].getPosition()].displayCharacter = "";
// }
anot[anpos].colour = new java.awt.Color(ae[aa].getColour());
if (firstColour == null)
{
firstColour = anot[anpos].colour;
}
}
}
jalview.datamodel.AlignmentAnnotation jaa = null;
if (an[i].getGraph())
{
float llim = 0, hlim = 0;
// if (autoForView || an[i].isAutoCalculated()) {
// hlim=11f;
// }
jaa = new jalview.datamodel.AlignmentAnnotation(an[i].getLabel(),
an[i].getDescription(), anot, llim, hlim,
an[i].getGraphType());
jaa.graphGroup = an[i].getGraphGroup();
jaa._linecolour = firstColour;
if (an[i].getThresholdLine() != null)
{
jaa.setThreshold(new jalview.datamodel.GraphLine(an[i]
.getThresholdLine().getValue(), an[i]
.getThresholdLine().getLabel(), new java.awt.Color(
an[i].getThresholdLine().getColour())));
}
if (autoForView || an[i].isAutoCalculated())
{
// Hardwire the symbol display line to ensure that labels for
// histograms are displayed
jaa.hasText = true;
}
}
else
{
jaa = new jalview.datamodel.AlignmentAnnotation(an[i].getLabel(),
an[i].getDescription(), anot);
jaa._linecolour = firstColour;
}
// register new annotation
if (an[i].getId() != null)
{
annotationIds.put(an[i].getId(), jaa);
jaa.annotationId = an[i].getId();
}
// recover sequence association
if (an[i].getSequenceRef() != null)
{
if (al.findName(an[i].getSequenceRef()) != null)
{
jaa.createSequenceMapping(al.findName(an[i].getSequenceRef()),
1, true);
al.findName(an[i].getSequenceRef()).addAlignmentAnnotation(jaa);
}
}
// and make a note of any group association
if (an[i].getGroupRef() != null && an[i].getGroupRef().length() > 0)
{
ArrayList aal = groupAnnotRefs
.get(an[i].getGroupRef());
if (aal == null)
{
aal = new ArrayList();
groupAnnotRefs.put(an[i].getGroupRef(), aal);
}
aal.add(jaa);
}
if (an[i].hasScore())
{
jaa.setScore(an[i].getScore());
}
if (an[i].hasVisible())
{
jaa.visible = an[i].getVisible();
}
if (an[i].hasCentreColLabels())
{
jaa.centreColLabels = an[i].getCentreColLabels();
}
if (an[i].hasScaleColLabels())
{
jaa.scaleColLabel = an[i].getScaleColLabels();
}
if (an[i].hasAutoCalculated() && an[i].isAutoCalculated())
{
// newer files have an 'autoCalculated' flag and store calculation
// state in viewport properties
jaa.autoCalculated = true; // means annotation will be marked for
// update at end of load.
}
if (an[i].hasGraphHeight())
{
jaa.graphHeight = an[i].getGraphHeight();
}
if (an[i].hasBelowAlignment())
{
jaa.belowAlignment = an[i].isBelowAlignment();
}
jaa.setCalcId(an[i].getCalcId());
if (an[i].getPropertyCount() > 0)
{
for (jalview.schemabinding.version2.Property prop : an[i]
.getProperty())
{
jaa.setProperty(prop.getName(), prop.getValue());
}
}
if (jaa.autoCalculated)
{
autoAlan.add(new JvAnnotRow(i, jaa));
}
else
// if (!autoForView)
{
// add autocalculated group annotation and any user created annotation
// for the view
al.addAnnotation(jaa);
}
}
}
// ///////////////////////
// LOAD GROUPS
// Create alignment markup and styles for this view
if (jms.getJGroupCount() > 0)
{
JGroup[] groups = jms.getJGroup();
boolean addAnnotSchemeGroup = false;
for (int i = 0; i < groups.length; i++)
{
ColourSchemeI cs = null;
if (groups[i].getColour() != null)
{
if (groups[i].getColour().startsWith("ucs"))
{
cs = getUserColourScheme(jms, groups[i].getColour());
}
else if (groups[i].getColour().equals("AnnotationColourGradient")
&& groups[i].getAnnotationColours() != null)
{
addAnnotSchemeGroup = true;
cs = null;
}
else
{
cs = ColourSchemeProperty.getColour(al, groups[i].getColour());
}
if (cs != null)
{
cs.setThreshold(groups[i].getPidThreshold(), true);
}
}
Vector seqs = new Vector();
for (int s = 0; s < groups[i].getSeqCount(); s++)
{
String seqId = groups[i].getSeq(s) + "";
jalview.datamodel.SequenceI ts = seqRefIds.get(seqId);
if (ts != null)
{
seqs.addElement(ts);
}
}
if (seqs.size() < 1)
{
continue;
}
jalview.datamodel.SequenceGroup sg = new jalview.datamodel.SequenceGroup(
seqs, groups[i].getName(), cs, groups[i].getDisplayBoxes(),
groups[i].getDisplayText(), groups[i].getColourText(),
groups[i].getStart(), groups[i].getEnd());
sg.setOutlineColour(new java.awt.Color(groups[i].getOutlineColour()));
sg.textColour = new java.awt.Color(groups[i].getTextCol1());
sg.textColour2 = new java.awt.Color(groups[i].getTextCol2());
sg.setShowNonconserved(groups[i].hasShowUnconserved() ? groups[i]
.isShowUnconserved() : false);
sg.thresholdTextColour = groups[i].getTextColThreshold();
if (groups[i].hasShowConsensusHistogram())
{
sg.setShowConsensusHistogram(groups[i].isShowConsensusHistogram());
}
;
if (groups[i].hasShowSequenceLogo())
{
sg.setshowSequenceLogo(groups[i].isShowSequenceLogo());
}
if (groups[i].hasNormaliseSequenceLogo())
{
sg.setNormaliseSequenceLogo(groups[i].isNormaliseSequenceLogo());
}
if (groups[i].hasIgnoreGapsinConsensus())
{
sg.setIgnoreGapsConsensus(groups[i].getIgnoreGapsinConsensus());
}
if (groups[i].getConsThreshold() != 0)
{
jalview.analysis.Conservation c = new jalview.analysis.Conservation(
"All", ResidueProperties.propHash, 3,
sg.getSequences(null), 0, sg.getWidth() - 1);
c.calculate();
c.verdict(false, 25);
sg.cs.setConservation(c);
}
if (groups[i].getId() != null && groupAnnotRefs.size() > 0)
{
// re-instate unique group/annotation row reference
ArrayList jaal = groupAnnotRefs
.get(groups[i].getId());
if (jaal != null)
{
for (jalview.datamodel.AlignmentAnnotation jaa : jaal)
{
jaa.groupRef = sg;
if (jaa.autoCalculated)
{
// match up and try to set group autocalc alignment row for this
// annotation
if (jaa.label.startsWith("Consensus for "))
{
sg.setConsensus(jaa);
}
// match up and try to set group autocalc alignment row for this
// annotation
if (jaa.label.startsWith("Conservation for "))
{
sg.setConservationRow(jaa);
}
}
}
}
}
al.addGroup(sg);
if (addAnnotSchemeGroup)
{
// reconstruct the annotation colourscheme
sg.cs = constructAnnotationColour(
groups[i].getAnnotationColours(), null, al, jms, false);
}
}
}
if (view == null)
{
// only dataset in this model, so just return.
return null;
}
// ///////////////////////////////
// LOAD VIEWPORT
// If we just load in the same jar file again, the sequenceSetId
// will be the same, and we end up with multiple references
// to the same sequenceSet. We must modify this id on load
// so that each load of the file gives a unique id
String uniqueSeqSetId = view.getSequenceSetId() + uniqueSetSuffix;
String viewId = (view.getId() == null ? null : view.getId()
+ uniqueSetSuffix);
AlignFrame af = null;
AlignViewport av = null;
// now check to see if we really need to create a new viewport.
if (multipleView && viewportsAdded.size() == 0)
{
// We recovered an alignment for which a viewport already exists.
// TODO: fix up any settings necessary for overlaying stored state onto
// state recovered from another document. (may not be necessary).
// we may need a binding from a viewport in memory to one recovered from
// XML.
// and then recover its containing af to allow the settings to be applied.
// TODO: fix for vamsas demo
System.err
.println("About to recover a viewport for existing alignment: Sequence set ID is "
+ uniqueSeqSetId);
Object seqsetobj = retrieveExistingObj(uniqueSeqSetId);
if (seqsetobj != null)
{
if (seqsetobj instanceof String)
{
uniqueSeqSetId = (String) seqsetobj;
System.err
.println("Recovered extant sequence set ID mapping for ID : New Sequence set ID is "
+ uniqueSeqSetId);
}
else
{
System.err
.println("Warning : Collision between sequence set ID string and existing jalview object mapping.");
}
}
}
/**
* indicate that annotation colours are applied across all groups (pre
* Jalview 2.8.1 behaviour)
*/
boolean doGroupAnnColour = isVersionStringLaterThan("2.8.1",
object.getVersion());
AlignmentPanel ap = null;
boolean isnewview = true;
if (viewId != null)
{
// Check to see if this alignment already has a view id == viewId
jalview.gui.AlignmentPanel views[] = Desktop
.getAlignmentPanels(uniqueSeqSetId);
if (views != null && views.length > 0)
{
for (int v = 0; v < views.length; v++)
{
if (views[v].av.getViewId().equalsIgnoreCase(viewId))
{
// recover the existing alignpanel, alignframe, viewport
af = views[v].alignFrame;
av = views[v].av;
ap = views[v];
// TODO: could even skip resetting view settings if we don't want to
// change the local settings from other jalview processes
isnewview = false;
}
}
}
}
if (isnewview)
{
af = loadViewport(file, jseqs, hiddenSeqs, al, jms, view,
uniqueSeqSetId, viewId, autoAlan);
av = af.viewport;
ap = af.alignPanel;
}
// LOAD TREES
// /////////////////////////////////////
if (loadTreesAndStructures && jms.getTreeCount() > 0)
{
try
{
for (int t = 0; t < jms.getTreeCount(); t++)
{
Tree tree = jms.getTree(t);
TreePanel tp = (TreePanel) retrieveExistingObj(tree.getId());
if (tp == null)
{
tp = af.ShowNewickTree(
new jalview.io.NewickFile(tree.getNewick()),
tree.getTitle(), tree.getWidth(), tree.getHeight(),
tree.getXpos(), tree.getYpos());
if (tree.getId() != null)
{
// perhaps bind the tree id to something ?
}
}
else
{
// update local tree attributes ?
// TODO: should check if tp has been manipulated by user - if so its
// settings shouldn't be modified
tp.setTitle(tree.getTitle());
tp.setBounds(new Rectangle(tree.getXpos(), tree.getYpos(), tree
.getWidth(), tree.getHeight()));
tp.av = av; // af.viewport; // TODO: verify 'associate with all
// views'
// works still
tp.treeCanvas.av = av; // af.viewport;
tp.treeCanvas.ap = ap; // af.alignPanel;
}
if (tp == null)
{
warn("There was a problem recovering stored Newick tree: \n"
+ tree.getNewick());
continue;
}
tp.fitToWindow.setState(tree.getFitToWindow());
tp.fitToWindow_actionPerformed(null);
if (tree.getFontName() != null)
{
tp.setTreeFont(new java.awt.Font(tree.getFontName(), tree
.getFontStyle(), tree.getFontSize()));
}
else
{
tp.setTreeFont(new java.awt.Font(view.getFontName(), view
.getFontStyle(), tree.getFontSize()));
}
tp.showPlaceholders(tree.getMarkUnlinked());
tp.showBootstrap(tree.getShowBootstrap());
tp.showDistances(tree.getShowDistances());
tp.treeCanvas.threshold = tree.getThreshold();
if (tree.getCurrentTree())
{
af.viewport.setCurrentTree(tp.getTree());
}
}
} catch (Exception ex)
{
ex.printStackTrace();
}
}
// //LOAD STRUCTURES
if (loadTreesAndStructures)
{
loadStructures(jprovider, jseqs, af, ap);
}
// and finally return.
return af;
}
/**
* Load and link any saved structure viewers.
*
* @param jprovider
* @param jseqs
* @param af
* @param ap
*/
protected void loadStructures(jarInputStreamProvider jprovider,
JSeq[] jseqs, AlignFrame af, AlignmentPanel ap)
{
/*
* Run through all PDB ids on the alignment, and collect mappings between
* distinct view ids and all sequences referring to that view.
*/
Map structureViewers = new LinkedHashMap();
for (int i = 0; i < jseqs.length; i++)
{
if (jseqs[i].getPdbidsCount() > 0)
{
Pdbids[] ids = jseqs[i].getPdbids();
for (int p = 0; p < ids.length; p++)
{
final int structureStateCount = ids[p].getStructureStateCount();
for (int s = 0; s < structureStateCount; s++)
{
// check to see if we haven't already created this structure view
final StructureState structureState = ids[p]
.getStructureState(s);
String sviewid = (structureState.getViewId() == null) ? null
: structureState.getViewId() + uniqueSetSuffix;
jalview.datamodel.PDBEntry jpdb = new jalview.datamodel.PDBEntry();
// Originally : ids[p].getFile()
// : TODO: verify external PDB file recovery still works in normal
// jalview project load
jpdb.setFile(loadPDBFile(jprovider, ids[p].getId()));
jpdb.setId(ids[p].getId());
int x = structureState.getXpos();
int y = structureState.getYpos();
int width = structureState.getWidth();
int height = structureState.getHeight();
// Probably don't need to do this anymore...
// Desktop.desktop.getComponentAt(x, y);
// TODO: NOW: check that this recovers the PDB file correctly.
String pdbFile = loadPDBFile(jprovider, ids[p].getId());
jalview.datamodel.SequenceI seq = seqRefIds.get(jseqs[i]
.getId() + "");
if (sviewid == null)
{
sviewid = "_jalview_pre2_4_" + x + "," + y + "," + width
+ "," + height;
}
if (!structureViewers.containsKey(sviewid))
{
structureViewers.put(sviewid, new StructureViewerModel(x, y,
width, height, false, false, true));
// Legacy pre-2.7 conversion JAL-823 :
// do not assume any view has to be linked for colour by
// sequence
}
// assemble String[] { pdb files }, String[] { id for each
// file }, orig_fileloc, SequenceI[][] {{ seqs_file 1 }, {
// seqs_file 2}, boolean[] {
// linkAlignPanel,superposeWithAlignpanel}} from hash
StructureViewerModel jmoldat = structureViewers.get(sviewid);
jmoldat.setAlignWithPanel(jmoldat.isAlignWithPanel()
| (structureState.hasAlignwithAlignPanel() ? structureState
.getAlignwithAlignPanel() : false));
/*
* Default colour by linked panel to false if not specified (e.g.
* for pre-2.7 projects)
*/
boolean colourWithAlignPanel = jmoldat.isColourWithAlignPanel();
colourWithAlignPanel |= (structureState
.hasColourwithAlignPanel() ? structureState
.getColourwithAlignPanel() : false);
jmoldat.setColourWithAlignPanel(colourWithAlignPanel);
/*
* Default colour by viewer to true if not specified (e.g. for
* pre-2.7 projects)
*/
boolean colourByViewer = jmoldat.isColourByViewer();
colourByViewer &= structureState.hasColourByJmol() ? structureState
.getColourByJmol() : true;
jmoldat.setColourByViewer(colourByViewer);
if (jmoldat.getStateData().length() < structureState
.getContent().length())
{
{
jmoldat.setStateData(structureState.getContent());
}
}
if (ids[p].getFile() != null)
{
File mapkey = new File(ids[p].getFile());
StructureData seqstrmaps = jmoldat.getFileData().get(mapkey);
if (seqstrmaps == null)
{
jmoldat.getFileData().put(
mapkey,
seqstrmaps = jmoldat.new StructureData(pdbFile,
ids[p].getId()));
}
if (!seqstrmaps.getSeqList().contains(seq))
{
seqstrmaps.getSeqList().add(seq);
// TODO and chains?
}
}
else
{
errorMessage = ("The Jmol views in this project were imported\nfrom an older version of Jalview.\nPlease review the sequence colour associations\nin the Colour by section of the Jmol View menu.\n\nIn the case of problems, see note at\nhttp://issues.jalview.org/browse/JAL-747");
warn(errorMessage);
}
}
}
}
}
// Instantiate the associated structure views
for (Entry entry : structureViewers
.entrySet())
{
createOrLinkStructureViewer(entry, af, ap);
}
}
/**
*
* @param viewerData
* @param af
* @param ap
*/
protected void createOrLinkStructureViewer(
Entry viewerData, AlignFrame af,
AlignmentPanel ap)
{
final StructureViewerModel svattrib = viewerData.getValue();
/*
* Search for any viewer windows already open from other alignment views
* that exactly match the stored structure state
*/
StructureViewerBase comp = findMatchingViewer(viewerData);
if (comp != null)
{
linkStructureViewer(ap, comp, svattrib);
return;
}
/*
* Pending an XML element for ViewerType, just check if stateData contains
* "chimera" (part of the chimera session filename).
*/
if (svattrib.getStateData().indexOf("chimera") > -1)
{
createChimeraViewer(viewerData, af);
}
else
{
createJmolViewer(viewerData, af);
}
}
/**
* Create a new Chimera viewer.
*
* @param viewerData
* @param af
*/
protected void createChimeraViewer(
Entry viewerData, AlignFrame af)
{
final StructureViewerModel data = viewerData.getValue();
String chimeraSession = data.getStateData();
if (new File(chimeraSession).exists())
{
Set> fileData = data.getFileData()
.entrySet();
List pdbs = new ArrayList();
List allseqs = new ArrayList();
for (Entry pdb : fileData)
{
String filePath = pdb.getValue().getFilePath();
String pdbId = pdb.getValue().getPdbId();
// pdbs.add(new PDBEntry(filePath, pdbId));
pdbs.add(new PDBEntry(pdbId, null, PDBEntry.Type.PDB, filePath));
final List seqList = pdb.getValue().getSeqList();
SequenceI[] seqs = seqList.toArray(new SequenceI[seqList.size()]);
allseqs.add(seqs);
}
boolean colourByChimera = data.isColourByViewer();
boolean colourBySequence = data.isColourWithAlignPanel();
// TODO can/should this be done via StructureViewer (like Jmol)?
final PDBEntry[] pdbArray = pdbs.toArray(new PDBEntry[pdbs.size()]);
final SequenceI[][] seqsArray = allseqs.toArray(new SequenceI[allseqs
.size()][]);
new ChimeraViewFrame(chimeraSession, af.alignPanel, pdbArray,
seqsArray, colourByChimera, colourBySequence);
}
else
{
Cache.log.error("Chimera session file " + chimeraSession
+ " not found");
}
}
/**
* Create a new Jmol window. First parse the Jmol state to translate filenames
* loaded into the view, and record the order in which files are shown in the
* Jmol view, so we can add the sequence mappings in same order.
*
* @param viewerData
* @param af
*/
protected void createJmolViewer(
final Entry viewerData,
AlignFrame af)
{
final StructureViewerModel svattrib = viewerData.getValue();
String state = svattrib.getStateData();
List pdbfilenames = new ArrayList();
List seqmaps = new ArrayList();
List pdbids = new ArrayList();
StringBuilder newFileLoc = new StringBuilder(64);
int cp = 0, ncp, ecp;
Map oldFiles = svattrib.getFileData();
while ((ncp = state.indexOf("load ", cp)) > -1)
{
do
{
// look for next filename in load statement
newFileLoc.append(state.substring(cp,
ncp = (state.indexOf("\"", ncp + 1) + 1)));
String oldfilenam = state.substring(ncp,
ecp = state.indexOf("\"", ncp));
// recover the new mapping data for this old filename
// have to normalize filename - since Jmol and jalview do
// filename
// translation differently.
StructureData filedat = oldFiles.get(new File(oldfilenam));
newFileLoc.append(Platform.escapeString(filedat.getFilePath()));
pdbfilenames.add(filedat.getFilePath());
pdbids.add(filedat.getPdbId());
seqmaps.add(filedat.getSeqList().toArray(new SequenceI[0]));
newFileLoc.append("\"");
cp = ecp + 1; // advance beyond last \" and set cursor so we can
// look for next file statement.
} while ((ncp = state.indexOf("/*file*/", cp)) > -1);
}
if (cp > 0)
{
// just append rest of state
newFileLoc.append(state.substring(cp));
}
else
{
System.err.print("Ignoring incomplete Jmol state for PDB ids: ");
newFileLoc = new StringBuilder(state);
newFileLoc.append("; load append ");
for (File id : oldFiles.keySet())
{
// add this and any other pdb files that should be present in
// the viewer
StructureData filedat = oldFiles.get(id);
newFileLoc.append(filedat.getFilePath());
pdbfilenames.add(filedat.getFilePath());
pdbids.add(filedat.getPdbId());
seqmaps.add(filedat.getSeqList().toArray(new SequenceI[0]));
newFileLoc.append(" \"");
newFileLoc.append(filedat.getFilePath());
newFileLoc.append("\"");
}
newFileLoc.append(";");
}
if (newFileLoc.length() > 0)
{
int histbug = newFileLoc.indexOf("history = ");
histbug += 10;
int diff = histbug == -1 ? -1 : newFileLoc.indexOf(";", histbug);
String val = (diff == -1) ? null : newFileLoc
.substring(histbug, diff);
if (val != null && val.length() >= 4)
{
if (val.contains("e"))
{
if (val.trim().equals("true"))
{
val = "1";
}
else
{
val = "0";
}
newFileLoc.replace(histbug, diff, val);
}
}
final String[] pdbf = pdbfilenames.toArray(new String[pdbfilenames
.size()]);
final String[] id = pdbids.toArray(new String[pdbids.size()]);
final SequenceI[][] sq = seqmaps
.toArray(new SequenceI[seqmaps.size()][]);
final String fileloc = newFileLoc.toString();
final String sviewid = viewerData.getKey();
final AlignFrame alf = af;
final Rectangle rect = new Rectangle(svattrib.getX(),
svattrib.getY(), svattrib.getWidth(), svattrib.getHeight());
try
{
javax.swing.SwingUtilities.invokeAndWait(new Runnable()
{
@Override
public void run()
{
JalviewStructureDisplayI sview = null;
try
{
// JAL-1333 note - we probably can't migrate Jmol views to UCSF
// Chimera!
sview = new StructureViewer(alf.alignPanel
.getStructureSelectionManager()).createView(
StructureViewer.ViewerType.JMOL, pdbf, id, sq,
alf.alignPanel, svattrib, fileloc, rect, sviewid);
addNewStructureViewer(sview);
} catch (OutOfMemoryError ex)
{
new OOMWarning("restoring structure view for PDB id " + id,
(OutOfMemoryError) ex.getCause());
if (sview != null && sview.isVisible())
{
sview.closeViewer(false);
sview.setVisible(false);
sview.dispose();
}
}
}
});
} catch (InvocationTargetException ex)
{
warn("Unexpected error when opening Jmol view.", ex);
} catch (InterruptedException e)
{
// e.printStackTrace();
}
}
}
/**
* Returns any open frame that matches given structure viewer data. The match
* is based on the unique viewId, or (for older project versions) the frame's
* geometry.
*
* @param viewerData
* @return
*/
protected StructureViewerBase findMatchingViewer(
Entry viewerData)
{
final String sviewid = viewerData.getKey();
final StructureViewerModel svattrib = viewerData.getValue();
StructureViewerBase comp = null;
JInternalFrame[] frames = getAllFrames();
for (JInternalFrame frame : frames)
{
if (frame instanceof StructureViewerBase)
{
/*
* Post jalview 2.4 schema includes structure view id
*/
if (sviewid != null
&& ((StructureViewerBase) frame).getViewId()
.equals(sviewid))
{
comp = (AppJmol) frame;
// todo: break?
}
/*
* Otherwise test for matching position and size of viewer frame
*/
else if (frame.getX() == svattrib.getX()
&& frame.getY() == svattrib.getY()
&& frame.getHeight() == svattrib.getHeight()
&& frame.getWidth() == svattrib.getWidth())
{
comp = (AppJmol) frame;
// todo: break?
}
}
}
return comp;
}
/**
* Link an AlignmentPanel to an existing structure viewer.
*
* @param ap
* @param viewer
* @param oldFiles
* @param useinViewerSuperpos
* @param usetoColourbyseq
* @param viewerColouring
*/
protected void linkStructureViewer(AlignmentPanel ap,
StructureViewerBase viewer, StructureViewerModel svattrib)
{
// NOTE: if the jalview project is part of a shared session then
// view synchronization should/could be done here.
final boolean useinViewerSuperpos = svattrib.isAlignWithPanel();
final boolean usetoColourbyseq = svattrib.isColourWithAlignPanel();
final boolean viewerColouring = svattrib.isColourByViewer();
Map oldFiles = svattrib.getFileData();
/*
* Add mapping for sequences in this view to an already open viewer
*/
final AAStructureBindingModel binding = viewer.getBinding();
for (File id : oldFiles.keySet())
{
// add this and any other pdb files that should be present in the
// viewer
StructureData filedat = oldFiles.get(id);
String pdbFile = filedat.getFilePath();
SequenceI[] seq = filedat.getSeqList().toArray(new SequenceI[0]);
binding.getSsm().setMapping(seq, null, pdbFile,
jalview.io.AppletFormatAdapter.FILE);
binding.addSequenceForStructFile(pdbFile, seq);
}
// and add the AlignmentPanel's reference to the view panel
viewer.addAlignmentPanel(ap);
if (useinViewerSuperpos)
{
viewer.useAlignmentPanelForSuperposition(ap);
}
else
{
viewer.excludeAlignmentPanelForSuperposition(ap);
}
if (usetoColourbyseq)
{
viewer.useAlignmentPanelForColourbyseq(ap, !viewerColouring);
}
else
{
viewer.excludeAlignmentPanelForColourbyseq(ap);
}
}
/**
* Get all frames within the Desktop.
*
* @return
*/
protected JInternalFrame[] getAllFrames()
{
JInternalFrame[] frames = null;
// TODO is this necessary - is it safe - risk of hanging?
do
{
try
{
frames = Desktop.desktop.getAllFrames();
} catch (ArrayIndexOutOfBoundsException e)
{
// occasional No such child exceptions are thrown here...
try
{
Thread.sleep(10);
} catch (InterruptedException f)
{
}
}
} while (frames == null);
return frames;
}
/**
*
* @param supported
* - minimum version we are comparing against
* @param version
* - version of data being processsed.
* @return true if version is development/null or evaluates to the same or
* later X.Y.Z (where X,Y,Z are like [0-9]+b?[0-9]*)
*/
private boolean isVersionStringLaterThan(String supported, String version)
{
if (version == null || version.equalsIgnoreCase("DEVELOPMENT BUILD")
|| version.equalsIgnoreCase("Test")
|| version.equalsIgnoreCase("AUTOMATED BUILD"))
{
System.err.println("Assuming project file with "
+ (version == null ? "null" : version)
+ " is compatible with Jalview version " + supported);
return true;
}
else
{
StringTokenizer currentV = new StringTokenizer(supported, "."), fileV = new StringTokenizer(
version, ".");
while (currentV.hasMoreTokens() && fileV.hasMoreTokens())
{
// convert b to decimal to catch bugfix releases within a series
String curT = currentV.nextToken().toLowerCase().replace('b', '.');
String fileT = fileV.nextToken().toLowerCase().replace('b', '.');
try
{
if (Float.valueOf(curT) > Float.valueOf(fileT))
{
// current version is newer than the version that wrote the file
return false;
}
} catch (NumberFormatException nfe)
{
System.err
.println("** WARNING: Version comparison failed for tokens ("
+ curT
+ ") and ("
+ fileT
+ ")\n** Current: '"
+ supported + "' and Version: '" + version + "'");
}
}
if (currentV.hasMoreElements())
{
// fileV has no minor version but identical series to current
return false;
}
}
return true;
}
Vector newStructureViewers = null;
protected void addNewStructureViewer(JalviewStructureDisplayI sview)
{
if (newStructureViewers != null)
{
sview.getBinding().setFinishedLoadingFromArchive(false);
newStructureViewers.add(sview);
}
}
protected void setLoadingFinishedForNewStructureViewers()
{
if (newStructureViewers != null)
{
for (JalviewStructureDisplayI sview : newStructureViewers)
{
sview.getBinding().setFinishedLoadingFromArchive(true);
}
newStructureViewers.clear();
newStructureViewers = null;
}
}
AlignFrame loadViewport(String file, JSeq[] JSEQ,
List hiddenSeqs, Alignment al,
JalviewModelSequence jms, Viewport view, String uniqueSeqSetId,
String viewId, List autoAlan)
{
AlignFrame af = null;
af = new AlignFrame(al, view.getWidth(), view.getHeight(),
uniqueSeqSetId, viewId);
af.setFileName(file, "Jalview");
for (int i = 0; i < JSEQ.length; i++)
{
af.viewport.setSequenceColour(af.viewport.getAlignment()
.getSequenceAt(i), new java.awt.Color(JSEQ[i].getColour()));
}
af.viewport.setGatherViewsHere(view.getGatheredViews());
if (view.getSequenceSetId() != null)
{
AlignmentViewport av = viewportsAdded.get(uniqueSeqSetId);
af.viewport.setSequenceSetId(uniqueSeqSetId);
if (av != null)
{
// propagate shared settings to this new view
af.viewport.setHistoryList(av.getHistoryList());
af.viewport.setRedoList(av.getRedoList());
}
else
{
viewportsAdded.put(uniqueSeqSetId, af.viewport);
}
// TODO: check if this method can be called repeatedly without
// side-effects if alignpanel already registered.
PaintRefresher.Register(af.alignPanel, uniqueSeqSetId);
}
// apply Hidden regions to view.
if (hiddenSeqs != null)
{
for (int s = 0; s < JSEQ.length; s++)
{
jalview.datamodel.SequenceGroup hidden = new jalview.datamodel.SequenceGroup();
for (int r = 0; r < JSEQ[s].getHiddenSequencesCount(); r++)
{
hidden.addSequence(
al.getSequenceAt(JSEQ[s].getHiddenSequences(r)), false);
}
af.viewport.hideRepSequences(al.getSequenceAt(s), hidden);
}
// jalview.datamodel.SequenceI[] hseqs = new
// jalview.datamodel.SequenceI[hiddenSeqs
// .size()];
//
// for (int s = 0; s < hiddenSeqs.size(); s++)
// {
// hseqs[s] = (jalview.datamodel.SequenceI) hiddenSeqs.elementAt(s);
// }
SequenceI[] hseqs = hiddenSeqs.toArray(new SequenceI[hiddenSeqs
.size()]);
af.viewport.hideSequence(hseqs);
}
// recover view properties and display parameters
if (view.getViewName() != null)
{
af.viewport.viewName = view.getViewName();
af.setInitialTabVisible();
}
af.setBounds(view.getXpos(), view.getYpos(), view.getWidth(),
view.getHeight());
af.viewport.setShowAnnotation(view.getShowAnnotation());
af.viewport.setAbovePIDThreshold(view.getPidSelected());
af.viewport.setColourText(view.getShowColourText());
af.viewport.setConservationSelected(view.getConservationSelected());
af.viewport.setShowJVSuffix(view.getShowFullId());
af.viewport.setRightAlignIds(view.getRightAlignIds());
af.viewport.setFont(
new java.awt.Font(view.getFontName(), view.getFontStyle(), view
.getFontSize()), true);
// TODO: allow custom charWidth/Heights to be restored by updating them
// after setting font - which means set above to false
af.viewport.setRenderGaps(view.getRenderGaps());
af.viewport.setWrapAlignment(view.getWrapAlignment());
af.viewport.setShowAnnotation(view.getShowAnnotation());
af.viewport.setShowBoxes(view.getShowBoxes());
af.viewport.setShowText(view.getShowText());
af.viewport.setTextColour(new java.awt.Color(view.getTextCol1()));
af.viewport.setTextColour2(new java.awt.Color(view.getTextCol2()));
af.viewport.setThresholdTextColour(view.getTextColThreshold());
af.viewport.setShowUnconserved(view.hasShowUnconserved() ? view
.isShowUnconserved() : false);
af.viewport.setStartRes(view.getStartRes());
af.viewport.setStartSeq(view.getStartSeq());
af.alignPanel.updateLayout();
ColourSchemeI cs = null;
// apply colourschemes
if (view.getBgColour() != null)
{
if (view.getBgColour().startsWith("ucs"))
{
cs = getUserColourScheme(jms, view.getBgColour());
}
else if (view.getBgColour().startsWith("Annotation"))
{
AnnotationColours viewAnnColour = view.getAnnotationColours();
cs = constructAnnotationColour(viewAnnColour, af, al, jms, true);
// annpos
}
else
{
cs = ColourSchemeProperty.getColour(al, view.getBgColour());
}
if (cs != null)
{
cs.setThreshold(view.getPidThreshold(), true);
cs.setConsensus(af.viewport.getSequenceConsensusHash());
}
}
af.viewport.setGlobalColourScheme(cs);
af.viewport.setColourAppliesToAllGroups(false);
if (view.getConservationSelected() && cs != null)
{
cs.setConservationInc(view.getConsThreshold());
}
af.changeColour(cs);
af.viewport.setColourAppliesToAllGroups(true);
af.viewport.setShowSequenceFeatures(view.getShowSequenceFeatures());
if (view.hasCentreColumnLabels())
{
af.viewport.setCentreColumnLabels(view.getCentreColumnLabels());
}
if (view.hasIgnoreGapsinConsensus())
{
af.viewport.setIgnoreGapsConsensus(view.getIgnoreGapsinConsensus(),
null);
}
if (view.hasFollowHighlight())
{
af.viewport.setFollowHighlight(view.getFollowHighlight());
}
if (view.hasFollowSelection())
{
af.viewport.followSelection = view.getFollowSelection();
}
if (view.hasShowConsensusHistogram())
{
af.viewport.setShowConsensusHistogram(view
.getShowConsensusHistogram());
}
else
{
af.viewport.setShowConsensusHistogram(true);
}
if (view.hasShowSequenceLogo())
{
af.viewport.setShowSequenceLogo(view.getShowSequenceLogo());
}
else
{
af.viewport.setShowSequenceLogo(false);
}
if (view.hasNormaliseSequenceLogo())
{
af.viewport.setNormaliseSequenceLogo(view.getNormaliseSequenceLogo());
}
if (view.hasShowDbRefTooltip())
{
af.viewport.setShowDBRefs(view.getShowDbRefTooltip());
}
if (view.hasShowNPfeatureTooltip())
{
af.viewport.setShowNPFeats(view.hasShowNPfeatureTooltip());
}
if (view.hasShowGroupConsensus())
{
af.viewport.setShowGroupConsensus(view.getShowGroupConsensus());
}
else
{
af.viewport.setShowGroupConsensus(false);
}
if (view.hasShowGroupConservation())
{
af.viewport.setShowGroupConservation(view.getShowGroupConservation());
}
else
{
af.viewport.setShowGroupConservation(false);
}
// recover featre settings
if (jms.getFeatureSettings() != null)
{
FeaturesDisplayed fdi;
af.viewport.setFeaturesDisplayed(fdi = new FeaturesDisplayed());
String[] renderOrder = new String[jms.getFeatureSettings()
.getSettingCount()];
Hashtable featureGroups = new Hashtable();
Hashtable featureColours = new Hashtable();
Hashtable featureOrder = new Hashtable();
for (int fs = 0; fs < jms.getFeatureSettings().getSettingCount(); fs++)
{
Setting setting = jms.getFeatureSettings().getSetting(fs);
if (setting.hasMincolour())
{
GraduatedColor gc = setting.hasMin() ? new GraduatedColor(
new java.awt.Color(setting.getMincolour()),
new java.awt.Color(setting.getColour()),
setting.getMin(), setting.getMax()) : new GraduatedColor(
new java.awt.Color(setting.getMincolour()),
new java.awt.Color(setting.getColour()), 0, 1);
if (setting.hasThreshold())
{
gc.setThresh(setting.getThreshold());
gc.setThreshType(setting.getThreshstate());
}
gc.setAutoScaled(true); // default
if (setting.hasAutoScale())
{
gc.setAutoScaled(setting.getAutoScale());
}
if (setting.hasColourByLabel())
{
gc.setColourByLabel(setting.getColourByLabel());
}
// and put in the feature colour table.
featureColours.put(setting.getType(), gc);
}
else
{
featureColours.put(setting.getType(),
new java.awt.Color(setting.getColour()));
}
renderOrder[fs] = setting.getType();
if (setting.hasOrder())
{
featureOrder.put(setting.getType(), setting.getOrder());
}
else
{
featureOrder.put(setting.getType(), new Float(fs
/ jms.getFeatureSettings().getSettingCount()));
}
if (setting.getDisplay())
{
fdi.setVisible(setting.getType());
}
}
Hashtable fgtable = new Hashtable();
for (int gs = 0; gs < jms.getFeatureSettings().getGroupCount(); gs++)
{
Group grp = jms.getFeatureSettings().getGroup(gs);
fgtable.put(grp.getName(), new Boolean(grp.getDisplay()));
}
// FeatureRendererSettings frs = new FeatureRendererSettings(renderOrder,
// fgtable, featureColours, jms.getFeatureSettings().hasTransparency() ?
// jms.getFeatureSettings().getTransparency() : 0.0, featureOrder);
FeatureRendererSettings frs = new FeatureRendererSettings(
renderOrder, fgtable, featureColours, 1.0f, featureOrder);
af.alignPanel.getSeqPanel().seqCanvas.getFeatureRenderer()
.transferSettings(frs);
}
if (view.getHiddenColumnsCount() > 0)
{
for (int c = 0; c < view.getHiddenColumnsCount(); c++)
{
af.viewport.hideColumns(view.getHiddenColumns(c).getStart(), view
.getHiddenColumns(c).getEnd() // +1
);
}
}
if (view.getCalcIdParam() != null)
{
for (CalcIdParam calcIdParam : view.getCalcIdParam())
{
if (calcIdParam != null)
{
if (recoverCalcIdParam(calcIdParam, af.viewport))
{
}
else
{
warn("Couldn't recover parameters for "
+ calcIdParam.getCalcId());
}
}
}
}
af.setMenusFromViewport(af.viewport);
// TODO: we don't need to do this if the viewport is aready visible.
/*
* Add the AlignFrame to the desktop (it may be 'gathered' later), unless it
* has a 'cdna/protein complement' view, in which case save it in order to
* populate a SplitFrame once all views have been read in.
*/
String complementaryViewId = view.getComplementId();
if (complementaryViewId == null)
{
Desktop.addInternalFrame(af, view.getTitle(), view.getWidth(),
view.getHeight());
// recompute any autoannotation
af.alignPanel.updateAnnotation(false, true);
reorderAutoannotation(af, al, autoAlan);
af.alignPanel.alignmentChanged();
}
else
{
splitFrameCandidates.put(view, af);
}
return af;
}
private ColourSchemeI constructAnnotationColour(
AnnotationColours viewAnnColour, AlignFrame af, Alignment al,
JalviewModelSequence jms, boolean checkGroupAnnColour)
{
boolean propagateAnnColour = false;
ColourSchemeI cs = null;
AlignmentI annAlignment = af != null ? af.viewport.getAlignment() : al;
if (checkGroupAnnColour && al.getGroups() != null
&& al.getGroups().size() > 0)
{
// pre 2.8.1 behaviour
// check to see if we should transfer annotation colours
propagateAnnColour = true;
for (jalview.datamodel.SequenceGroup sg : al.getGroups())
{
if (sg.cs instanceof AnnotationColourGradient)
{
propagateAnnColour = false;
}
}
}
// int find annotation
if (annAlignment.getAlignmentAnnotation() != null)
{
for (int i = 0; i < annAlignment.getAlignmentAnnotation().length; i++)
{
if (annAlignment.getAlignmentAnnotation()[i].label
.equals(viewAnnColour.getAnnotation()))
{
if (annAlignment.getAlignmentAnnotation()[i].getThreshold() == null)
{
annAlignment.getAlignmentAnnotation()[i]
.setThreshold(new jalview.datamodel.GraphLine(
viewAnnColour.getThreshold(), "Threshold",
java.awt.Color.black)
);
}
if (viewAnnColour.getColourScheme().equals("None"))
{
cs = new AnnotationColourGradient(
annAlignment.getAlignmentAnnotation()[i],
new java.awt.Color(viewAnnColour.getMinColour()),
new java.awt.Color(viewAnnColour.getMaxColour()),
viewAnnColour.getAboveThreshold());
}
else if (viewAnnColour.getColourScheme().startsWith("ucs"))
{
cs = new AnnotationColourGradient(
annAlignment.getAlignmentAnnotation()[i],
getUserColourScheme(jms,
viewAnnColour.getColourScheme()),
viewAnnColour.getAboveThreshold());
}
else
{
cs = new AnnotationColourGradient(
annAlignment.getAlignmentAnnotation()[i],
ColourSchemeProperty.getColour(al,
viewAnnColour.getColourScheme()),
viewAnnColour.getAboveThreshold());
}
if (viewAnnColour.hasPerSequence())
{
((AnnotationColourGradient) cs).setSeqAssociated(viewAnnColour
.isPerSequence());
}
if (viewAnnColour.hasPredefinedColours())
{
((AnnotationColourGradient) cs)
.setPredefinedColours(viewAnnColour
.isPredefinedColours());
}
if (propagateAnnColour && al.getGroups() != null)
{
// Also use these settings for all the groups
for (int g = 0; g < al.getGroups().size(); g++)
{
jalview.datamodel.SequenceGroup sg = al.getGroups().get(g);
if (sg.cs == null)
{
continue;
}
/*
* if (viewAnnColour.getColourScheme().equals("None" )) { sg.cs =
* new AnnotationColourGradient(
* annAlignment.getAlignmentAnnotation()[i], new
* java.awt.Color(viewAnnColour. getMinColour()), new
* java.awt.Color(viewAnnColour. getMaxColour()),
* viewAnnColour.getAboveThreshold()); } else
*/
{
sg.cs = new AnnotationColourGradient(
annAlignment.getAlignmentAnnotation()[i], sg.cs,
viewAnnColour.getAboveThreshold());
if (cs instanceof AnnotationColourGradient)
{
if (viewAnnColour.hasPerSequence())
{
((AnnotationColourGradient) cs)
.setSeqAssociated(viewAnnColour.isPerSequence());
}
if (viewAnnColour.hasPredefinedColours())
{
((AnnotationColourGradient) cs)
.setPredefinedColours(viewAnnColour
.isPredefinedColours());
}
}
}
}
}
break;
}
}
}
return cs;
}
private void reorderAutoannotation(AlignFrame af, Alignment al,
List autoAlan)
{
// copy over visualization settings for autocalculated annotation in the
// view
if (al.getAlignmentAnnotation() != null)
{
/**
* Kludge for magic autoannotation names (see JAL-811)
*/
String[] magicNames = new String[]
{ "Consensus", "Quality", "Conservation" };
JvAnnotRow nullAnnot = new JvAnnotRow(-1, null);
Hashtable visan = new Hashtable();
for (String nm : magicNames)
{
visan.put(nm, nullAnnot);
}
for (JvAnnotRow auan : autoAlan)
{
visan.put(auan.template.label
+ (auan.template.getCalcId() == null ? "" : "\t"
+ auan.template.getCalcId()), auan);
}
int hSize = al.getAlignmentAnnotation().length;
List reorder = new ArrayList();
// work through any autoCalculated annotation already on the view
// removing it if it should be placed in a different location on the
// annotation panel.
List remains = new ArrayList(visan.keySet());
for (int h = 0; h < hSize; h++)
{
jalview.datamodel.AlignmentAnnotation jalan = al
.getAlignmentAnnotation()[h];
if (jalan.autoCalculated)
{
String k;
JvAnnotRow valan = visan.get(k = jalan.label);
if (jalan.getCalcId() != null)
{
valan = visan.get(k = jalan.label + "\t" + jalan.getCalcId());
}
if (valan != null)
{
// delete the auto calculated row from the alignment
al.deleteAnnotation(jalan, false);
remains.remove(k);
hSize--;
h--;
if (valan != nullAnnot)
{
if (jalan != valan.template)
{
// newly created autoannotation row instance
// so keep a reference to the visible annotation row
// and copy over all relevant attributes
if (valan.template.graphHeight >= 0)
{
jalan.graphHeight = valan.template.graphHeight;
}
jalan.visible = valan.template.visible;
}
reorder.add(new JvAnnotRow(valan.order, jalan));
}
}
}
}
// Add any (possibly stale) autocalculated rows that were not appended to
// the view during construction
for (String other : remains)
{
JvAnnotRow othera = visan.get(other);
if (othera != nullAnnot && othera.template.getCalcId() != null
&& othera.template.getCalcId().length() > 0)
{
reorder.add(othera);
}
}
// now put the automatic annotation in its correct place
int s = 0, srt[] = new int[reorder.size()];
JvAnnotRow[] rws = new JvAnnotRow[reorder.size()];
for (JvAnnotRow jvar : reorder)
{
rws[s] = jvar;
srt[s++] = jvar.order;
}
reorder.clear();
jalview.util.QuickSort.sort(srt, rws);
// and re-insert the annotation at its correct position
for (JvAnnotRow jvar : rws)
{
al.addAnnotation(jvar.template, jvar.order);
}
af.alignPanel.adjustAnnotationHeight();
}
}
Hashtable skipList = null;
/**
* TODO remove this method
*
* @param view
* @return AlignFrame bound to sequenceSetId from view, if one exists. private
* AlignFrame getSkippedFrame(Viewport view) { if (skipList==null) {
* throw new Error("Implementation Error. No skipList defined for this
* Jalview2XML instance."); } return (AlignFrame)
* skipList.get(view.getSequenceSetId()); }
*/
/**
* Check if the Jalview view contained in object should be skipped or not.
*
* @param object
* @return true if view's sequenceSetId is a key in skipList
*/
private boolean skipViewport(JalviewModel object)
{
if (skipList == null)
{
return false;
}
String id;
if (skipList.containsKey(id = object.getJalviewModelSequence()
.getViewport()[0].getSequenceSetId()))
{
if (Cache.log != null && Cache.log.isDebugEnabled())
{
Cache.log.debug("Skipping seuqence set id " + id);
}
return true;
}
return false;
}
public void addToSkipList(AlignFrame af)
{
if (skipList == null)
{
skipList = new Hashtable();
}
skipList.put(af.getViewport().getSequenceSetId(), af);
}
public void clearSkipList()
{
if (skipList != null)
{
skipList.clear();
skipList = null;
}
}
private void recoverDatasetFor(SequenceSet vamsasSet, Alignment al,
boolean ignoreUnrefed)
{
jalview.datamodel.Alignment ds = getDatasetFor(vamsasSet.getDatasetId());
Vector dseqs = null;
if (ds == null)
{
// create a list of new dataset sequences
dseqs = new Vector();
}
for (int i = 0, iSize = vamsasSet.getSequenceCount(); i < iSize; i++)
{
Sequence vamsasSeq = vamsasSet.getSequence(i);
ensureJalviewDatasetSequence(vamsasSeq, ds, dseqs, ignoreUnrefed);
}
// create a new dataset
if (ds == null)
{
SequenceI[] dsseqs = new SequenceI[dseqs.size()];
dseqs.copyInto(dsseqs);
ds = new jalview.datamodel.Alignment(dsseqs);
debug("Created new dataset " + vamsasSet.getDatasetId()
+ " for alignment " + System.identityHashCode(al));
addDatasetRef(vamsasSet.getDatasetId(), ds);
}
// set the dataset for the newly imported alignment.
if (al.getDataset() == null && !ignoreUnrefed)
{
al.setDataset(ds);
}
}
/**
*
* @param vamsasSeq
* sequence definition to create/merge dataset sequence for
* @param ds
* dataset alignment
* @param dseqs
* vector to add new dataset sequence to
*/
private void ensureJalviewDatasetSequence(Sequence vamsasSeq,
AlignmentI ds, Vector dseqs, boolean ignoreUnrefed)
{
// JBP TODO: Check this is called for AlCodonFrames to support recovery of
// xRef Codon Maps
SequenceI sq = seqRefIds.get(vamsasSeq.getId());
SequenceI dsq = null;
if (sq != null && sq.getDatasetSequence() != null)
{
dsq = sq.getDatasetSequence();
}
if (sq == null && ignoreUnrefed)
{
return;
}
String sqid = vamsasSeq.getDsseqid();
if (dsq == null)
{
// need to create or add a new dataset sequence reference to this sequence
if (sqid != null)
{
dsq = seqRefIds.get(sqid);
}
// check again
if (dsq == null)
{
// make a new dataset sequence
dsq = sq.createDatasetSequence();
if (sqid == null)
{
// make up a new dataset reference for this sequence
sqid = seqHash(dsq);
}
dsq.setVamsasId(uniqueSetSuffix + sqid);
seqRefIds.put(sqid, dsq);
if (ds == null)
{
if (dseqs != null)
{
dseqs.addElement(dsq);
}
}
else
{
ds.addSequence(dsq);
}
}
else
{
if (sq != dsq)
{ // make this dataset sequence sq's dataset sequence
sq.setDatasetSequence(dsq);
// and update the current dataset alignment
if (ds == null)
{
if (dseqs != null)
{
if (!dseqs.contains(dsq))
{
dseqs.add(dsq);
}
}
else
{
if (ds.findIndex(dsq) < 0)
{
ds.addSequence(dsq);
}
}
}
}
}
}
// TODO: refactor this as a merge dataset sequence function
// now check that sq (the dataset sequence) sequence really is the union of
// all references to it
// boolean pre = sq.getStart() < dsq.getStart();
// boolean post = sq.getEnd() > dsq.getEnd();
// if (pre || post)
if (sq != dsq)
{
// StringBuffer sb = new StringBuffer();
String newres = jalview.analysis.AlignSeq.extractGaps(
jalview.util.Comparison.GapChars, sq.getSequenceAsString());
if (!newres.equalsIgnoreCase(dsq.getSequenceAsString())
&& newres.length() > dsq.getLength())
{
// Update with the longer sequence.
synchronized (dsq)
{
/*
* if (pre) { sb.insert(0, newres .substring(0, dsq.getStart() -
* sq.getStart())); dsq.setStart(sq.getStart()); } if (post) {
* sb.append(newres.substring(newres.length() - sq.getEnd() -
* dsq.getEnd())); dsq.setEnd(sq.getEnd()); }
*/
dsq.setSequence(newres);
}
// TODO: merges will never happen if we 'know' we have the real dataset
// sequence - this should be detected when id==dssid
System.err
.println("DEBUG Notice: Merged dataset sequence (if you see this often, post at http://issues.jalview.org/browse/JAL-1474)"); // ("
// + (pre ? "prepended" : "") + " "
// + (post ? "appended" : ""));
}
}
}
java.util.Hashtable datasetIds = null;
java.util.IdentityHashMap dataset2Ids = null;
private Alignment getDatasetFor(String datasetId)
{
if (datasetIds == null)
{
datasetIds = new Hashtable();
return null;
}
if (datasetIds.containsKey(datasetId))
{
return (Alignment) datasetIds.get(datasetId);
}
return null;
}
private void addDatasetRef(String datasetId, Alignment dataset)
{
if (datasetIds == null)
{
datasetIds = new Hashtable();
}
datasetIds.put(datasetId, dataset);
}
/**
* make a new dataset ID for this jalview dataset alignment
*
* @param dataset
* @return
*/
private String getDatasetIdRef(jalview.datamodel.Alignment dataset)
{
if (dataset.getDataset() != null)
{
warn("Serious issue! Dataset Object passed to getDatasetIdRef is not a Jalview DATASET alignment...");
}
String datasetId = makeHashCode(dataset, null);
if (datasetId == null)
{
// make a new datasetId and record it
if (dataset2Ids == null)
{
dataset2Ids = new IdentityHashMap();
}
else
{
datasetId = (String) dataset2Ids.get(dataset);
}
if (datasetId == null)
{
datasetId = "ds" + dataset2Ids.size() + 1;
dataset2Ids.put(dataset, datasetId);
}
}
return datasetId;
}
private void addDBRefs(SequenceI datasetSequence, Sequence sequence)
{
for (int d = 0; d < sequence.getDBRefCount(); d++)
{
DBRef dr = sequence.getDBRef(d);
jalview.datamodel.DBRefEntry entry = new jalview.datamodel.DBRefEntry(
sequence.getDBRef(d).getSource(), sequence.getDBRef(d)
.getVersion(), sequence.getDBRef(d).getAccessionId());
if (dr.getMapping() != null)
{
entry.setMap(addMapping(dr.getMapping()));
}
datasetSequence.addDBRef(entry);
}
}
private jalview.datamodel.Mapping addMapping(Mapping m)
{
SequenceI dsto = null;
// Mapping m = dr.getMapping();
int fr[] = new int[m.getMapListFromCount() * 2];
Enumeration f = m.enumerateMapListFrom();
for (int _i = 0; f.hasMoreElements(); _i += 2)
{
MapListFrom mf = (MapListFrom) f.nextElement();
fr[_i] = mf.getStart();
fr[_i + 1] = mf.getEnd();
}
int fto[] = new int[m.getMapListToCount() * 2];
f = m.enumerateMapListTo();
for (int _i = 0; f.hasMoreElements(); _i += 2)
{
MapListTo mf = (MapListTo) f.nextElement();
fto[_i] = mf.getStart();
fto[_i + 1] = mf.getEnd();
}
jalview.datamodel.Mapping jmap = new jalview.datamodel.Mapping(dsto,
fr, fto, (int) m.getMapFromUnit(), (int) m.getMapToUnit());
if (m.getMappingChoice() != null)
{
MappingChoice mc = m.getMappingChoice();
if (mc.getDseqFor() != null)
{
String dsfor = "" + mc.getDseqFor();
if (seqRefIds.containsKey(dsfor))
{
/**
* recover from hash
*/
jmap.setTo(seqRefIds.get(dsfor));
}
else
{
frefedSequence.add(new Object[]
{ dsfor, jmap });
}
}
else
{
/**
* local sequence definition
*/
Sequence ms = mc.getSequence();
SequenceI djs = null;
String sqid = ms.getDsseqid();
if (sqid != null && sqid.length() > 0)
{
/*
* recover dataset sequence
*/
djs = seqRefIds.get(sqid);
}
else
{
System.err
.println("Warning - making up dataset sequence id for DbRef sequence map reference");
sqid = ((Object) ms).toString(); // make up a new hascode for
// undefined dataset sequence hash
// (unlikely to happen)
}
if (djs == null)
{
/**
* make a new dataset sequence and add it to refIds hash
*/
djs = new jalview.datamodel.Sequence(ms.getName(),
ms.getSequence());
djs.setStart(jmap.getMap().getToLowest());
djs.setEnd(jmap.getMap().getToHighest());
djs.setVamsasId(uniqueSetSuffix + sqid);
jmap.setTo(djs);
seqRefIds.put(sqid, djs);
}
jalview.bin.Cache.log.debug("about to recurse on addDBRefs.");
addDBRefs(djs, ms);
}
}
return (jmap);
}
public jalview.gui.AlignmentPanel copyAlignPanel(AlignmentPanel ap,
boolean keepSeqRefs)
{
initSeqRefs();
jalview.schemabinding.version2.JalviewModel jm = saveState(ap, null,
null);
if (!keepSeqRefs)
{
clearSeqRefs();
jm.getJalviewModelSequence().getViewport(0).setSequenceSetId(null);
}
else
{
uniqueSetSuffix = "";
jm.getJalviewModelSequence().getViewport(0).setId(null); // we don't
// overwrite the
// view we just
// copied
}
if (this.frefedSequence == null)
{
frefedSequence = new Vector();
}
viewportsAdded.clear();
AlignFrame af = loadFromObject(jm, null, false, null);
af.alignPanels.clear();
af.closeMenuItem_actionPerformed(true);
/*
* if(ap.av.getAlignment().getAlignmentAnnotation()!=null) { for(int i=0;
* i();
}
if (seqsToIds == null)
{
seqsToIds = new IdentityHashMap();
}
seqRefIds.put(jv2vobj.get(jvobj).toString(), (SequenceI) jvobj);
seqsToIds.put((SequenceI) jvobj, id);
}
else if (jvobj instanceof jalview.datamodel.AlignmentAnnotation)
{
String anid;
AlignmentAnnotation jvann = (AlignmentAnnotation) jvobj;
annotationIds.put(anid = jv2vobj.get(jvobj).toString(), jvann);
if (jvann.annotationId == null)
{
jvann.annotationId = anid;
}
if (!jvann.annotationId.equals(anid))
{
// TODO verify that this is the correct behaviour
this.warn("Overriding Annotation ID for " + anid
+ " from different id : " + jvann.annotationId);
jvann.annotationId = anid;
}
}
else if (jvobj instanceof String)
{
if (jvids2vobj == null)
{
jvids2vobj = new Hashtable();
jvids2vobj.put(jvobj, jv2vobj.get(jvobj).toString());
}
}
else
{
Cache.log.debug("Ignoring " + jvobj.getClass() + " (ID = " + id);
}
}
}
/**
* set the uniqueSetSuffix used to prefix/suffix object IDs for jalview
* objects created from the project archive. If string is null (default for
* construction) then suffix will be set automatically.
*
* @param string
*/
public void setUniqueSetSuffix(String string)
{
uniqueSetSuffix = string;
}
/**
* uses skipList2 as the skipList for skipping views on sequence sets
* associated with keys in the skipList
*
* @param skipList2
*/
public void setSkipList(Hashtable skipList2)
{
skipList = skipList2;
}
}