+/*
+ VARNA is a tool for the automated drawing, visualization and annotation of the secondary structure of RNA, designed as a companion software for web servers and databases.
+ Copyright (C) 2008 Kevin Darty, Alain Denise and Yann Ponty.
+ electronic mail : Yann.Ponty@lri.fr
+ paper mail : LRI, bat 490 Université Paris-Sud 91405 Orsay Cedex France
+
+ This file is part of VARNA version 3.1.
+ VARNA version 3.1 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.
+
+ VARNA version 3.1 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 VARNA version 3.1.
+ If not, see http://www.gnu.org/licenses.
+ */
+package fr.orsay.lri.varna.models.rna;
+
+import java.awt.Color;
+import java.awt.Point;
+import java.awt.geom.GeneralPath;
+import java.awt.geom.Point2D;
+import java.awt.geom.Rectangle2D;
+import java.io.ByteArrayInputStream;
+import java.io.ByteArrayOutputStream;
+import java.io.FileOutputStream;
+import java.io.FileWriter;
+import java.io.IOException;
+import java.io.ObjectInputStream;
+import java.io.ObjectOutputStream;
+import java.io.OutputStreamWriter;
+import java.io.Reader;
+import java.io.Serializable;
+import java.io.StreamTokenizer;
+import java.util.ArrayList;
+import java.util.Arrays;
+import java.util.Collection;
+import java.util.Comparator;
+import java.util.HashSet;
+import java.util.Hashtable;
+import java.util.List;
+import java.util.Set;
+import java.util.Stack;
+import java.util.Vector;
+
+import javax.xml.transform.sax.TransformerHandler;
+
+import org.xml.sax.SAXException;
+import org.xml.sax.helpers.AttributesImpl;
+
+import fr.orsay.lri.varna.VARNAPanel;
+import fr.orsay.lri.varna.applications.templateEditor.Couple;
+import fr.orsay.lri.varna.exceptions.ExceptionExportFailed;
+import fr.orsay.lri.varna.exceptions.ExceptionFileFormatOrSyntax;
+import fr.orsay.lri.varna.exceptions.ExceptionNAViewAlgorithm;
+import fr.orsay.lri.varna.exceptions.ExceptionPermissionDenied;
+import fr.orsay.lri.varna.exceptions.ExceptionUnmatchedClosingParentheses;
+import fr.orsay.lri.varna.exceptions.ExceptionWritingForbidden;
+import fr.orsay.lri.varna.factories.RNAFactory;
+import fr.orsay.lri.varna.interfaces.InterfaceVARNAListener;
+import fr.orsay.lri.varna.interfaces.InterfaceVARNAObservable;
+import fr.orsay.lri.varna.models.VARNAConfig;
+import fr.orsay.lri.varna.models.VARNAConfig.BP_STYLE;
+import fr.orsay.lri.varna.models.annotations.ChemProbAnnotation;
+import fr.orsay.lri.varna.models.annotations.HighlightRegionAnnotation;
+import fr.orsay.lri.varna.models.annotations.TextAnnotation;
+import fr.orsay.lri.varna.models.export.PSExport;
+import fr.orsay.lri.varna.models.export.SVGExport;
+import fr.orsay.lri.varna.models.export.SecStrDrawingProducer;
+import fr.orsay.lri.varna.models.export.TikzExport;
+import fr.orsay.lri.varna.models.export.XFIGExport;
+import fr.orsay.lri.varna.models.naView.NAView;
+import fr.orsay.lri.varna.models.rna.ModeleBackboneElement.BackboneType;
+import fr.orsay.lri.varna.models.templates.DrawRNATemplateCurveMethod;
+import fr.orsay.lri.varna.models.templates.DrawRNATemplateMethod;
+import fr.orsay.lri.varna.models.templates.RNATemplate;
+import fr.orsay.lri.varna.models.templates.RNATemplateDrawingAlgorithmException;
+import fr.orsay.lri.varna.models.templates.RNATemplateMapping;
+import fr.orsay.lri.varna.utils.RNAMLParser;
+import fr.orsay.lri.varna.utils.XMLUtils;
+import fr.orsay.lri.varna.views.VueUI;
+
+/**
+ * The RNA model which contain the base list and the draw algorithm mode
+ *
+ * @author darty
+ *
+ */
+public class RNA extends InterfaceVARNAObservable implements Serializable {
+ /**
+ *
+ */
+ private static final long serialVersionUID = 7541274455751497303L;
+
+ /**
+ * Selects the "Feynman diagram" drawing algorithm that places the bases on
+ * a circle and draws the base-pairings as chords of the circle graph.
+ */
+
+ public static final int DRAW_MODE_CIRCULAR = 1;
+ /**
+ * Selects the "tree drawing" algorithm. Draws each loop on a circle whose
+ * radius depends on the number of bases involved in the loop. As some
+ * helices can be overlapping in the result, basic interaction is provided
+ * so that the user can "disentangle" the drawing by spinning the helices
+ * around the axis defined by their multiloop (bulge or internal loop)
+ * origin. This is roughly the initial placement strategy of RNAViz.
+ *
+ * @see <a href="http://rnaviz.sourceforge.net/">RNAViz</a>
+ */
+ public static final int DRAW_MODE_RADIATE = 2;
+
+ /**
+ * Selects the NAView algorithm.
+ */
+ public static final int DRAW_MODE_NAVIEW = 3;
+ /**
+ * Selects the linear algorithm.
+ */
+ public static final int DRAW_MODE_LINEAR = 4;
+
+ public static final int DRAW_MODE_VARNA_VIEW = 5;
+
+ /**
+ * Selects the RNAView algorithm.
+ */
+ public static final int DRAW_MODE_MOTIFVIEW = 6;
+
+ public static final int DRAW_MODE_TEMPLATE = 7;
+
+ public static final int DEFAULT_DRAW_MODE = DRAW_MODE_RADIATE;
+
+ public int BASE_RADIUS = 10;
+ public static final double LOOP_DISTANCE = 40.0; // distance between base
+ // pairs in an helix
+ public static final double BASE_PAIR_DISTANCE = 65.0; // distance between
+ // the two bases of
+ // a pair
+ public static final double MULTILOOP_DISTANCE = 35.0;
+ public static final double VIRTUAL_LOOP_RADIUS = 40.0;
+
+ public double CHEM_PROB_DIST = 14;
+ public double CHEM_PROB_BASE_LENGTH = 30;
+ public double CHEM_PROB_ARROW_HEIGHT = 10;
+ public double CHEM_PROB_ARROW_WIDTH = 5;
+ public double CHEM_PROB_TRIANGLE_WIDTH = 2.5;
+ public double CHEM_PROB_PIN_SEMIDIAG = 6;
+ public double CHEM_PROB_DOT_RADIUS = 6.;
+ public static double CHEM_PROB_ARROW_THICKNESS = 2.0;
+
+ public static ArrayList<String> NormalBases = new ArrayList<String>();
+ {
+ NormalBases.add("a");
+ NormalBases.add("c");
+ NormalBases.add("g");
+ NormalBases.add("u");
+ NormalBases.add("t");
+ }
+
+ public GeneralPath _debugShape = null;
+
+ /**
+ * The draw algorithm mode
+ */
+ private int _drawMode = DRAW_MODE_RADIATE;
+ private boolean _drawn = false;
+ private String _name = "";
+ private String _id = "";
+ public double _bpHeightIncrement = VARNAConfig.DEFAULT_BP_INCREMENT;
+ /**
+ * the base list
+ */
+ private ArrayList<ModeleBase> _listeBases;
+ /**
+ * the strand list
+ */
+ StructureTemp _listStrands = new StructureTemp();
+ /**
+ * Additional bonds and info can be specified here.
+ */
+ private ArrayList<ModeleBP> _structureAux = new ArrayList<ModeleBP>();
+ private ArrayList<TextAnnotation> _listeAnnotations = new ArrayList<TextAnnotation>();
+ private ArrayList<HighlightRegionAnnotation> _listeRegionHighlights = new ArrayList<HighlightRegionAnnotation>();
+ private ArrayList<ChemProbAnnotation> _chemProbAnnotations = new ArrayList<ChemProbAnnotation>();
+ private ModeleBackbone _backbone = new ModeleBackbone();
+
+ public static String XML_ELEMENT_NAME = "RNA";
+ public static String XML_VAR_BASE_SPACING_NAME = "spacing";
+ public static String XML_VAR_DRAWN_NAME = "drawn";
+ public static String XML_VAR_NAME_NAME = "name";
+ public static String XML_VAR_DRAWN_MODE_NAME = "mode";
+ public static String XML_VAR_ID_NAME = "id";
+ public static String XML_VAR_BP_HEIGHT_NAME = "delta";
+ public static String XML_VAR_BASES_NAME = "bases";
+ public static String XML_VAR_BASEPAIRS_NAME = "BPs";
+ public static String XML_VAR_ANNOTATIONS_NAME = "annotations";
+ public static String XML_VAR_BACKBONE_NAME = "backbone";
+
+ public void toXML(TransformerHandler hd) throws SAXException {
+ AttributesImpl atts = new AttributesImpl();
+ atts.addAttribute("", "", XML_VAR_DRAWN_NAME, "CDATA", "" + _drawn);
+ atts.addAttribute("", "", XML_VAR_DRAWN_MODE_NAME, "CDATA", ""
+ + _drawMode);
+ atts.addAttribute("", "", XML_VAR_ID_NAME, "CDATA", "" + _id);
+ atts.addAttribute("", "", XML_VAR_BP_HEIGHT_NAME, "CDATA", ""
+ + _bpHeightIncrement);
+ hd.startElement("", "", XML_ELEMENT_NAME, atts);
+
+ atts.clear();
+ hd.startElement("", "", XML_VAR_NAME_NAME, atts);
+ XMLUtils.exportCDATAString(hd, "" + _name);
+ hd.endElement("", "", XML_VAR_NAME_NAME);
+
+ atts.clear();
+ hd.startElement("", "", XML_VAR_BASES_NAME, atts);
+ for (ModeleBase mb : _listeBases) {
+ mb.toXML(hd);
+ }
+ hd.endElement("", "", XML_VAR_BASES_NAME);
+ atts.clear();
+
+ hd.startElement("", "", XML_VAR_BASEPAIRS_NAME, atts);
+ for (ModeleBP mbp : getSecStrBPs()) {
+ mbp.toXML(hd, true);
+ }
+ for (ModeleBP mbp : _structureAux) {
+ mbp.toXML(hd, false);
+ }
+ hd.endElement("", "", XML_VAR_BASEPAIRS_NAME);
+ atts.clear();
+
+ getBackbone().toXML(hd);
+ atts.clear();
+
+ hd.startElement("", "", XML_VAR_ANNOTATIONS_NAME, atts);
+ for (TextAnnotation ta : _listeAnnotations) {
+ ta.toXML(hd);
+ }
+ for (HighlightRegionAnnotation hra : _listeRegionHighlights) {
+ hra.toXML(hd);
+ }
+ for (ChemProbAnnotation cpa : _chemProbAnnotations) {
+ cpa.toXML(hd);
+ }
+ hd.endElement("", "", XML_VAR_ANNOTATIONS_NAME);
+ hd.endElement("", "", XML_ELEMENT_NAME);
+ }
+
+ public ModeleBackbone getBackbone() {
+ return _backbone;
+ }
+
+ public void setBackbone(ModeleBackbone b) {
+ _backbone = b;
+ }
+
+ transient private ArrayList<InterfaceVARNAListener> _listeVARNAListener = new ArrayList<InterfaceVARNAListener>();
+
+ public RNA() {
+ this("");
+ }
+
+ public RNA(String name) {
+ _name = name;
+ _listeBases = new ArrayList<ModeleBase>();
+ _drawn = false;
+ init();
+ }
+
+ public String toString() {
+ if (_name.equals("")) {
+ return getStructDBN();
+ } else {
+ return _name;
+ }
+ }
+
+ public RNA(RNA r) {
+ _drawMode = r._drawMode;
+ _listeBases.addAll(r._listeBases);
+ _listeVARNAListener = (ArrayList<InterfaceVARNAListener>) r._listeVARNAListener;
+ _drawn = r._drawn;
+ init();
+ }
+
+ public void init() {
+ }
+
+ public void saveRNADBN(String path, String title)
+ throws ExceptionWritingForbidden {
+ try {
+ FileWriter out = new FileWriter(path);
+ if (!title.equals("")) {
+ out.write("> " + title + "\n");
+ }
+ out.write(getListeBasesToString());
+ out.write('\n');
+ String str = "";
+ for (int i = 0; i < _listeBases.size(); i++) {
+ if (_listeBases.get(i).getElementStructure() == -1) {
+ str += '.';
+ } else {
+ if (_listeBases.get(i).getElementStructure() > i) {
+ str += '(';
+ } else {
+ str += ')';
+ }
+ }
+ }
+ out.write(str);
+ out.write('\n');
+ out.close();
+ } catch (IOException e) {
+ throw new ExceptionWritingForbidden(e.getMessage());
+ }
+ }
+
+ public Color getBaseInnerColor(int i, VARNAConfig conf) {
+ Color result = _listeBases.get(i).getStyleBase().getBaseInnerColor();
+ String res = _listeBases.get(i).getContent();
+ if (conf._drawColorMap) {
+ result = conf._cm.getColorForValue(_listeBases.get(i).getValue());
+ } else if ((conf._colorDashBases && (res.contains("-")))) {
+ result = conf._dashBasesColor;
+ } else if ((conf._colorSpecialBases && !NormalBases.contains(res
+ .toLowerCase()))) {
+ result = conf._specialBasesColor;
+ }
+ return result;
+ }
+
+ public Color getBaseOuterColor(int i, VARNAConfig conf) {
+ Color result = _listeBases.get(i).getStyleBase()
+ .getBaseOutlineColor();
+ return result;
+ }
+
+ private static double correctComponent(double c)
+ {
+ c = c / 255.0;
+ if (c <= 0.03928)
+ c = c/12.92;
+ else
+ c = Math.pow(((c+0.055)/1.055) , 2.4);
+ return c;
+ }
+ public static double getLuminance(Color c)
+ {
+ return 0.2126 * correctComponent(c.getRed()) + 0.7152 * correctComponent(c.getGreen()) + 0.0722 * correctComponent(c.getBlue());
+ }
+
+ public static boolean whiteLabelPreferrable(Color c)
+ {
+ if (getLuminance(c) > 0.32)
+ return false;
+ return true;
+ }
+
+
+
+ public Color getBaseNameColor(int i, VARNAConfig conf) {
+ Color result = _listeBases.get(i).getStyleBase().getBaseNameColor();
+ if ( RNA.whiteLabelPreferrable(getBaseInnerColor(i, conf)))
+ {
+ result=Color.white;
+ }
+
+ return result;
+ }
+
+ public Color getBasePairColor(ModeleBP bp, VARNAConfig conf) {
+ Color bondColor = conf._bondColor;
+ if (conf._useBaseColorsForBPs) {
+ bondColor = _listeBases.get(bp.getPartner5().getIndex())
+ .getStyleBase().getBaseInnerColor();
+ }
+ if (bp != null) {
+ bondColor = bp.getStyle().getColor(bondColor);
+ }
+ return bondColor;
+ }
+
+ public double getBasePairThickness(ModeleBP bp, VARNAConfig conf) {
+ double thickness = bp.getStyle().getThickness(conf._bpThickness);
+ return thickness;
+ }
+
+ private void drawSymbol(SecStrDrawingProducer out, double posx,
+ double posy, double normx, double normy, double radius,
+ boolean isCIS, ModeleBP.Edge e, double thickness) {
+ Color bck = out.getCurrentColor();
+ switch (e) {
+ case WC:
+ if (isCIS) {
+ out.fillCircle(posx, posy, (radius / 2.0), thickness, bck);
+ } else {
+ out.fillCircle(posx, posy, (radius / 2.0), thickness,
+ Color.white);
+ out.setColor(bck);
+ out.drawCircle(posx, posy, (radius / 2.0), thickness);
+ }
+ break;
+ case HOOGSTEEN: {
+ double xtab[] = new double[4];
+ double ytab[] = new double[4];
+ xtab[0] = posx - radius * normx / 2.0 - radius * normy / 2.0;
+ ytab[0] = posy - radius * normy / 2.0 + radius * normx / 2.0;
+ xtab[1] = posx + radius * normx / 2.0 - radius * normy / 2.0;
+ ytab[1] = posy + radius * normy / 2.0 + radius * normx / 2.0;
+ xtab[2] = posx + radius * normx / 2.0 + radius * normy / 2.0;
+ ytab[2] = posy + radius * normy / 2.0 - radius * normx / 2.0;
+ xtab[3] = posx - radius * normx / 2.0 + radius * normy / 2.0;
+ ytab[3] = posy - radius * normy / 2.0 - radius * normx / 2.0;
+ if (isCIS) {
+ out.fillPolygon(xtab, ytab, bck);
+ } else {
+ out.fillPolygon(xtab, ytab, Color.white);
+ out.setColor(bck);
+ out.drawPolygon(xtab, ytab, thickness);
+ }
+ }
+ break;
+ case SUGAR: {
+ double ix = radius * normx / 2.0;
+ double iy = radius * normy / 2.0;
+ double jx = radius * normy / 2.0;
+ double jy = -radius * normx / 2.0;
+ double xtab[] = new double[3];
+ double ytab[] = new double[3];
+ xtab[0] = posx - ix + jx;
+ ytab[0] = posy - iy + jy;
+ xtab[1] = posx + ix + jx;
+ ytab[1] = posy + iy + jy;
+ xtab[2] = posx - jx;
+ ytab[2] = posy - jy;
+
+ if (isCIS) {
+ out.fillPolygon(xtab, ytab, bck);
+ } else {
+ out.fillPolygon(xtab, ytab, Color.white);
+ out.setColor(bck);
+ out.drawPolygon(xtab, ytab, thickness);
+ }
+ }
+ break;
+ }
+ out.setColor(bck);
+ }
+
+ private void drawBasePairArc(SecStrDrawingProducer out, int i, int j,
+ Point2D.Double orig, Point2D.Double dest, ModeleBP style,
+ VARNAConfig conf) {
+ double coef;
+ double distance;
+ Point2D.Double center = new Point2D.Double((orig.x + dest.x)/2., (orig.y + dest.y)/2. + BASE_RADIUS);
+ if (j - i == 1)
+ coef = _bpHeightIncrement * 2;
+ else
+ coef = _bpHeightIncrement * 1;
+ distance = (int) Math.round(dest.x - orig.x);
+ if (conf._mainBPStyle != BP_STYLE.LW) {
+ out.drawArc(center, distance, distance * coef, 180, 0);
+ } else {
+ double thickness = getBasePairThickness(style, conf);
+ double radiusCircle = ((BASE_PAIR_DISTANCE - BASE_RADIUS) / 5.0);
+
+ if (style.isCanonical()) {
+ if (style.isCanonicalGC()) {
+ if ((orig.x != dest.x) || (orig.y != dest.y)) {
+ out.drawArc(center, distance - BASE_RADIUS / 2.,
+ distance * coef - BASE_RADIUS / 2, 180, 0);
+ out.drawArc(center, distance + BASE_RADIUS / 2.,
+ distance * coef + BASE_RADIUS / 2, 180, 0);
+ }
+ } else if (!style.isWobbleUG()) {
+ out.drawArc(center, distance, distance * coef, 180, 0);
+ drawSymbol(out, center.x, center.y + distance * coef / 2., 180., 0,
+ radiusCircle, style.isCIS(),
+ style.getEdgePartner5(), thickness);
+ } else {
+ out.drawArc(orig, distance, distance * coef, 180, 0);
+ }
+ } else {
+ ModeleBP.Edge p1 = style.getEdgePartner5();
+ ModeleBP.Edge p2 = style.getEdgePartner3();
+ out.drawArc(center, distance, distance * coef, 180, 0);
+ if (p1 == p2) {
+ drawSymbol(out, center.x, center.y + distance * coef / 2., 1., 0,
+ radiusCircle, style.isCIS(),
+ style.getEdgePartner5(), thickness);
+ } else {
+ drawSymbol(out, center.x - BASE_RADIUS,
+ center.y + distance * coef / 2., 1., 0, radiusCircle,
+ style.isCIS(), p1, thickness);
+ drawSymbol(out, center.x + BASE_RADIUS,
+ center.y + distance * coef / 2., 1., 0, radiusCircle,
+ style.isCIS(), p2, thickness);
+ }
+ }
+ }
+ }
+
+ private void drawBasePair(SecStrDrawingProducer out, Point2D.Double orig,
+ Point2D.Double dest, ModeleBP style, VARNAConfig conf) {
+ double dx = dest.x - orig.x;
+ double dy = dest.y - orig.y;
+ double dist = Math.sqrt((dest.x - orig.x) * (dest.x - orig.x)
+ + (dest.y - orig.y) * (dest.y - orig.y));
+ dx /= dist;
+ dy /= dist;
+ double nx = -dy;
+ double ny = dx;
+ orig = new Point2D.Double(orig.x + BASE_RADIUS * dx, orig.y
+ + BASE_RADIUS * dy);
+ dest = new Point2D.Double(dest.x - BASE_RADIUS * dx, dest.y
+ - BASE_RADIUS * dy);
+ if (conf._mainBPStyle == VARNAConfig.BP_STYLE.LW) {
+ double thickness = getBasePairThickness(style, conf);
+ double radiusCircle = ((BASE_PAIR_DISTANCE - BASE_RADIUS) / 5.0);
+
+ if (style.isCanonical()) {
+ if (style.isCanonicalGC()) {
+ if ((orig.x != dest.x) || (orig.y != dest.y)) {
+ nx *= BASE_RADIUS / 4.0;
+ ny *= BASE_RADIUS / 4.0;
+ out.drawLine((orig.x + nx), (orig.y + ny),
+ (dest.x + nx), (dest.y + ny), conf._bpThickness);
+ out.drawLine((orig.x - nx), (orig.y - ny),
+ (dest.x - nx), (dest.y - ny), conf._bpThickness);
+ }
+ } else if (style.isCanonicalAU()) {
+ out.drawLine(orig.x, orig.y, dest.x, dest.y,
+ conf._bpThickness);
+ } else if (style.isWobbleUG()) {
+ double cx = (dest.x + orig.x) / 2.0;
+ double cy = (dest.y + orig.y) / 2.0;
+ out.drawLine(orig.x, orig.y, dest.x, dest.y,
+ conf._bpThickness);
+ drawSymbol(out, cx, cy, nx, ny, radiusCircle, false,
+ ModeleBP.Edge.WC, thickness);
+ }
+
+ else {
+ double cx = (dest.x + orig.x) / 2.0;
+ double cy = (dest.y + orig.y) / 2.0;
+ out.drawLine(orig.x, orig.y, dest.x, dest.y,
+ conf._bpThickness);
+ drawSymbol(out, cx, cy, nx, ny, radiusCircle,
+ style.isCIS(), style.getEdgePartner5(), thickness);
+ }
+ } else {
+ ModeleBP.Edge p1 = style.getEdgePartner5();
+ ModeleBP.Edge p2 = style.getEdgePartner3();
+ double cx = (dest.x + orig.x) / 2.0;
+ double cy = (dest.y + orig.y) / 2.0;
+ out.drawLine(orig.x, orig.y, dest.x, dest.y, conf._bpThickness);
+ if (p1 == p2) {
+ drawSymbol(out, cx, cy, nx, ny, radiusCircle,
+ style.isCIS(), p1, thickness);
+ } else {
+ double vdx = (dest.x - orig.x);
+ double vdy = (dest.y - orig.y);
+ vdx /= 6.0;
+ vdy /= 6.0;
+ drawSymbol(out, cx + vdx, cy + vdy, nx, ny, radiusCircle,
+ style.isCIS(), p2, thickness);
+ drawSymbol(out, cx - vdx, cy - vdy, nx, ny, radiusCircle,
+ style.isCIS(), p1, thickness);
+ }
+ }
+ } else if (conf._mainBPStyle == VARNAConfig.BP_STYLE.RNAVIZ) {
+ double xcenter = (orig.x + dest.x) / 2.0;
+ double ycenter = (orig.y + dest.y) / 2.0;
+ out.fillCircle(xcenter, ycenter, 3.0 * conf._bpThickness,
+ conf._bpThickness, out.getCurrentColor());
+ } else if (conf._mainBPStyle == VARNAConfig.BP_STYLE.SIMPLE) {
+ out.drawLine(orig.x, orig.y, dest.x, dest.y, conf._bpThickness);
+ }
+ }
+
+ private void drawColorMap(VARNAConfig _conf, SecStrDrawingProducer out) {
+ double v1 = _conf._cm.getMinValue();
+ double v2 = _conf._cm.getMaxValue();
+ int x, y;
+ double xSpaceAvail = 0;
+ double ySpaceAvail = 0;
+ double thickness = 1.0;
+ /*
+ * ySpaceAvail =
+ * Math.min((getHeight()-rnabbox.height*scaleFactor-getTitleHeight
+ * ())/2.0,scaleFactor*(_conf._colorMapHeight+VARNAConfig.
+ * DEFAULT_COLOR_MAP_FONT_SIZE)); if ((int)ySpaceAvail==0) { xSpaceAvail
+ * =
+ * Math.min((getWidth()-rnabbox.width*scaleFactor)/2,scaleFactor*(_conf
+ * ._colorMapWidth)+VARNAConfig.DEFAULT_COLOR_MAP_STRIPE_WIDTH); }
+ */
+ Rectangle2D.Double currentBBox = out.getBoundingBox();
+
+ double xBase = (currentBBox.getMaxX() - _conf._colorMapWidth - _conf._colorMapXOffset);
+ // double yBase = (minY - _conf._colorMapHeight +
+ // _conf._colorMapYOffset);
+ double yBase = (currentBBox.getMinY() - _conf._colorMapHeight - VARNAConfig.DEFAULT_COLOR_MAP_FONT_SIZE);
+
+ for (int i = 0; i < _conf._colorMapWidth; i++) {
+ double ratio = (((double) i) / ((double) _conf._colorMapWidth - 1));
+ double val = v1 + (v2 - v1) * ratio;
+ Color c = _conf._cm.getColorForValue(val);
+ x = (int) (xBase + i);
+ y = (int) yBase;
+ out.fillRectangle(x, y, VARNAConfig.DEFAULT_COLOR_MAP_STRIPE_WIDTH,
+ _conf._colorMapHeight, c);
+ }
+ out.setColor(VARNAConfig.DEFAULT_COLOR_MAP_OUTLINE);
+ out.drawRectangle(xBase, yBase, (double) _conf._colorMapWidth
+ + VARNAConfig.DEFAULT_COLOR_MAP_STRIPE_WIDTH - 1,
+ _conf._colorMapHeight, thickness);
+
+ out.setColor(VARNAConfig.DEFAULT_COLOR_MAP_FONT_COLOR);
+ out.setFont(out.getCurrentFont(),
+ VARNAConfig.DEFAULT_COLOR_MAP_FONT_SIZE / 1.5);
+ out.drawText(xBase, yBase + _conf._colorMapHeight
+ + VARNAConfig.DEFAULT_COLOR_MAP_FONT_SIZE / 1.7,
+ "" + _conf._cm.getMinValue());
+ out.drawText(xBase + VARNAConfig.DEFAULT_COLOR_MAP_STRIPE_WIDTH
+ + _conf._colorMapWidth, yBase + _conf._colorMapHeight
+ + VARNAConfig.DEFAULT_COLOR_MAP_FONT_SIZE / 1.7,
+ "" + _conf._cm.getMaxValue());
+ out.drawText(
+ xBase
+ + (VARNAConfig.DEFAULT_COLOR_MAP_STRIPE_WIDTH + _conf._colorMapWidth)
+ / 2.0, yBase
+ - (VARNAConfig.DEFAULT_COLOR_MAP_FONT_SIZE / 1.7),
+ _conf._colorMapCaption);
+
+ }
+
+ private void renderRegionHighlights(SecStrDrawingProducer out,
+ Point2D.Double[] realCoords, Point2D.Double[] realCenters) {
+ for (HighlightRegionAnnotation r : _listeRegionHighlights) {
+ GeneralPath s = r.getShape(realCoords, realCenters, 1.0);
+ out.setColor(r.getFillColor());
+ out.fillPolygon(s, r.getFillColor());
+ out.setColor(r.getOutlineColor());
+ out.drawPolygon(s, 1l);
+ }
+
+ }
+
+ private void saveRNA(String path, VARNAConfig conf, double scale,
+ SecStrDrawingProducer out) throws ExceptionWritingForbidden {
+ out.setScale(scale);
+ // Computing bounding boxes
+ double EPSMargin = 40;
+ double minX = Double.MAX_VALUE;
+ double maxX = Double.MIN_VALUE;
+ double minY = Double.MAX_VALUE;
+ double maxY = Double.MIN_VALUE;
+
+ double x0, y0, x1, y1, xc, yc, xp, yp, dx, dy, norm;
+
+ for (int i = 0; i < _listeBases.size(); i++) {
+ minX = Math.min(minX, (_listeBases.get(i).getCoords().getX()
+ - BASE_RADIUS - EPSMargin));
+ minY = Math.min(minY, -(_listeBases.get(i).getCoords().getY()
+ - BASE_RADIUS - EPSMargin));
+ maxX = Math.max(maxX, (_listeBases.get(i).getCoords().getX()
+ + BASE_RADIUS + EPSMargin));
+ maxY = Math.max(maxY, -(_listeBases.get(i).getCoords().getY()
+ + BASE_RADIUS + EPSMargin));
+ }
+
+ // Rescaling everything
+ Point2D.Double[] coords = new Point2D.Double[_listeBases.size()];
+ Point2D.Double[] centers = new Point2D.Double[_listeBases.size()];
+ for (int i = 0; i < _listeBases.size(); i++) {
+ xp = (_listeBases.get(i).getCoords().getX() - minX);
+ yp = -(_listeBases.get(i).getCoords().getY() - minY);
+ coords[i] = new Point2D.Double(xp, yp);
+
+ Point2D.Double centerBck = getCenter(i);
+ if (get_drawMode() == RNA.DRAW_MODE_NAVIEW
+ || get_drawMode() == RNA.DRAW_MODE_RADIATE) {
+ if ((_listeBases.get(i).getElementStructure() != -1)
+ && i < _listeBases.size() - 1 && i > 1) {
+ ModeleBase b1 = get_listeBases().get(i - 1);
+ ModeleBase b2 = get_listeBases().get(i + 1);
+ int j1 = b1.getElementStructure();
+ int j2 = b2.getElementStructure();
+ if ((j1 == -1) ^ (j2 == -1)) {
+ // alors la position du nombre associé doit etre
+ // décalé
+ Point2D.Double a1 = b1.getCoords();
+ Point2D.Double a2 = b2.getCoords();
+ Point2D.Double c1 = b1.getCenter();
+ Point2D.Double c2 = b2.getCenter();
+
+ centerBck.x = _listeBases.get(i).getCoords().x
+ + (c1.x - a1.x) / c1.distance(a1)
+ + (c2.x - a2.x) / c2.distance(a2);
+ centerBck.y = _listeBases.get(i).getCoords().y
+ + (c1.y - a1.y) / c1.distance(a1)
+ + (c2.y - a2.y) / c2.distance(a2);
+ }
+ }
+ }
+ xc = (centerBck.getX() - minX);
+ yc = -(centerBck.getY() - minY);
+ centers[i] = new Point2D.Double(xc, yc);
+ }
+
+ // Drawing background
+ if (conf._drawBackground)
+ out.setBackgroundColor(conf._backgroundColor);
+
+ // Drawing region highlights
+ renderRegionHighlights(out, coords, centers);
+
+ // Drawing backbone
+ if (conf._drawBackbone)
+ {
+ for (int i = 1; i < _listeBases.size(); i++) {
+ Point2D.Double p1 = coords[i - 1];
+ Point2D.Double p2 = coords[i];
+ x0 = p1.x;
+ y0 = p1.y;
+ x1 = p2.x;
+ y1 = p2.y;
+ Point2D.Double vn = new Point2D.Double();
+ double dist = p1.distance(p2);
+ int a = _listeBases.get(i - 1).getElementStructure();
+ int b = _listeBases.get(i).getElementStructure();
+ BackboneType bt = _backbone.getTypeBefore(i);
+ boolean consecutivePair = (a == i) && (b == i - 1);
+
+ if (dist > 0) {
+ if (bt != BackboneType.DISCONTINUOUS_TYPE) {
+ Color c = _backbone.getColorBefore(i, conf._backboneColor);
+ if (bt == BackboneType.MISSING_PART_TYPE) {
+ c.brighter();
+ }
+ out.setColor(c);
+
+ vn.x = (x1 - x0) / dist;
+ vn.y = (y1 - y0) / dist;
+ if (consecutivePair
+ &&(getDrawMode() != RNA.DRAW_MODE_LINEAR)
+ && (getDrawMode() != RNA.DRAW_MODE_CIRCULAR)) {
+ int dir = 0;
+ if (i + 1 < coords.length) {
+ dir = (testDirectionality(i - 1, i, i + 1) ? 1 : -1);
+ } else if (i - 2 >= 0) {
+ dir = (testDirectionality(i - 2, i - 1, i) ? 1 : -1);
+ }
+ Point2D.Double centerSeg = new Point2D.Double(
+ (p1.x + p2.x) / 2.0, (p1.y + p2.y) / 2.0);
+ double centerDist = RNA.VIRTUAL_LOOP_RADIUS * scale;
+ Point2D.Double centerLoop = new Point2D.Double(
+ centerSeg.x + centerDist * dir * vn.y,
+ centerSeg.y - centerDist * dir * vn.x);
+ // Debug crosshair
+ //out.drawLine(centerLoop.x - 5, centerLoop.y,
+ // centerLoop.x + 5, centerLoop.y, 2.0);
+ //out.drawLine(centerLoop.x, centerLoop.y - 5,
+ // centerLoop.x, centerLoop.y + 5, 2.0);
+
+ double radius = centerLoop.distance(p1);
+ double a1 = 360.
+ * (Math.atan2(p1.y - centerLoop.y, p1.x - centerLoop.x))
+ / (2. * Math.PI);
+ double a2 = 360.
+ * (Math.atan2(p2.y - centerLoop.y, p2.x - centerLoop.x))
+ / (2. * Math.PI);
+ if (dir>0)
+ {
+ double tmp = a1;
+ a1 = a2;
+ a2 = tmp;
+ }
+ if (a1 < 0) {
+ a1 += 360.;
+ }
+ if (a2 < 0) {
+ a2 += 360.;
+ }
+ out.drawArc(centerLoop, 2. * radius, 2. * radius, a1,
+ a2);
+ } else {
+ out.drawLine((x0 + BASE_RADIUS * vn.x),
+ (y0 + BASE_RADIUS * vn.y), (x1 - BASE_RADIUS
+ * vn.x), (y1 - BASE_RADIUS * vn.y), 1.0);
+ }
+ }
+ }
+ }
+ }
+
+ // Drawing bonds
+ for (int i = 0; i < _listeBases.size(); i++) {
+ if (_listeBases.get(i).getElementStructure() > i) {
+ ModeleBP style = _listeBases.get(i).getStyleBP();
+ if (style.isCanonical() || conf._drawnNonCanonicalBP) {
+ Color bpcol = getBasePairColor(style, conf);
+ out.setColor(bpcol);
+
+ int j = _listeBases.get(i).getElementStructure();
+ x0 = coords[i].x;
+ y0 = coords[i].y;
+ x1 = coords[j].x;
+ y1 = coords[j].y;
+ dx = x1 - x0;
+ dy = y1 - y0;
+ norm = Math.sqrt(dx * dx + dy * dy);
+ dx /= norm;
+ dy /= norm;
+
+ if (_drawMode == DRAW_MODE_CIRCULAR
+ || _drawMode == DRAW_MODE_RADIATE
+ || _drawMode == DRAW_MODE_NAVIEW) {
+ drawBasePair(out, new Point2D.Double(x0, y0),
+ new Point2D.Double(x1, y1), style, conf);
+ } else if (_drawMode == DRAW_MODE_LINEAR) {
+ drawBasePairArc(out, i, j, new Point2D.Double(x0, y0),
+ new Point2D.Double(x1, y1), style, conf);
+ }
+ }
+ }
+ }
+
+ // Drawing additional bonds
+ if (conf._drawnNonPlanarBP) {
+ for (int i = 0; i < _structureAux.size(); i++) {
+ ModeleBP bp = _structureAux.get(i);
+ out.setColor(getBasePairColor(bp, conf));
+
+ int a = bp.getPartner5().getIndex();
+ int b = bp.getPartner3().getIndex();
+
+ if (bp.isCanonical() || conf._drawnNonCanonicalBP) {
+ x0 = coords[a].x;
+ y0 = coords[a].y;
+ x1 = coords[b].x;
+ y1 = coords[b].y;
+ dx = x1 - x0;
+ dy = y1 - y0;
+ norm = Math.sqrt(dx * dx + dy * dy);
+ dx /= norm;
+ dy /= norm;
+ if ((_drawMode == DRAW_MODE_CIRCULAR)
+ || (_drawMode == DRAW_MODE_RADIATE)
+ || _drawMode == DRAW_MODE_NAVIEW) {
+ drawBasePair(out, new Point2D.Double(x0, y0),
+ new Point2D.Double(x1, y1), bp, conf);
+ } else if (_drawMode == DRAW_MODE_LINEAR) {
+ drawBasePairArc(out, a, b, new Point2D.Double(x0, y0),
+ new Point2D.Double(x1, y1), bp, conf);
+ }
+ }
+ }
+ }
+
+ // Drawing Bases
+ double baseFontSize = (1.5 * BASE_RADIUS);
+ out.setFont(PSExport.FONT_HELVETICA_BOLD, baseFontSize);
+
+ for (int i = 0; i < _listeBases.size(); i++) {
+ x0 = coords[i].x;
+ y0 = coords[i].y;
+
+ Color baseInnerColor = getBaseInnerColor(i, conf);
+ Color baseOuterColor = getBaseOuterColor(i, conf);
+ Color baseNameColor = getBaseNameColor(i, conf);
+ if ( RNA.whiteLabelPreferrable(baseInnerColor))
+ {
+ baseNameColor=Color.white;
+ }
+
+
+ if (_listeBases.get(i) instanceof ModeleBasesComparison) {
+ ModeleBasesComparison mb = (ModeleBasesComparison) _listeBases
+ .get(i);
+ if (conf._fillBases) {
+ out.fillRectangle(x0 - 1.5 * BASE_RADIUS, y0 - BASE_RADIUS,
+ 3 * BASE_RADIUS, 2 * BASE_RADIUS,
+ baseInnerColor);
+ }
+ if (conf._drawOutlineBases) {
+ out.setColor(baseOuterColor);
+ out.drawRectangle(x0 - 1.5 * BASE_RADIUS, y0 - BASE_RADIUS,
+ 3 * BASE_RADIUS, 2 * BASE_RADIUS, 1l);
+ out.drawLine(x0, y0 - BASE_RADIUS, x0, y0 + BASE_RADIUS, 1l);
+ }
+
+ out.setColor(baseNameColor);
+ out.drawText(x0 - .75 * BASE_RADIUS, y0, "" + mb.getBase1());
+ out.drawText(x0 + .75 * BASE_RADIUS, y0, "" + mb.getBase2());
+ } else if (_listeBases.get(i) instanceof ModeleBaseNucleotide) {
+ if (conf._fillBases) {
+ out.fillCircle(x0, y0, BASE_RADIUS, 1l,
+ baseInnerColor);
+ }
+ if (conf._drawOutlineBases) {
+ out.setColor(baseOuterColor);
+ out.drawCircle(x0, y0, BASE_RADIUS, 1l);
+ }
+ out.setColor(baseNameColor);
+ out.drawText(x0, y0, _listeBases.get(i).getContent());
+
+ }
+ }
+
+ // Drawing base numbers
+ double numFontSize = (double) (1.5 * BASE_RADIUS);
+ out.setFont(PSExport.FONT_HELVETICA_BOLD, numFontSize);
+
+ for (int i = 0; i < _listeBases.size(); i++) {
+ int basenum = _listeBases.get(i).getBaseNumber();
+ if (basenum == -1) {
+ basenum = i + 1;
+ }
+ ModeleBase mb = _listeBases.get(i);
+ if (this.isNumberDrawn(mb, conf._numPeriod)) {
+ out.setColor(mb.getStyleBase()
+ .getBaseNumberColor());
+ x0 = coords[i].x;
+ y0 = coords[i].y;
+ x1 = centers[i].x;
+ y1 = centers[i].y;
+ dx = x1 - x0;
+ dy = y1 - y0;
+ norm = Math.sqrt(dx * dx + dy * dy);
+ dx /= norm;
+ dy /= norm;
+ Point2D.Double vn = VARNAPanel.computeExcentricUnitVector(i,coords,centers);
+
+ out.drawLine((x0 + 1.5 * BASE_RADIUS * vn.x), (y0 + 1.5
+ * BASE_RADIUS * vn.y), (x0 + 2.5 * BASE_RADIUS * vn.x),
+ (y0 + 2.5 * BASE_RADIUS * vn.y), 1);
+ out.drawText(
+ (x0 + (conf._distNumbers + 1.0) * BASE_RADIUS * vn.x),
+ (y0 + (conf._distNumbers + 1.0) * BASE_RADIUS * vn.y), mb.getLabel());
+ }
+ }
+ renderAnnotations(out, minX, minY, conf);
+
+ // Draw color map
+ if (conf._drawColorMap) {
+ drawColorMap(conf, out);
+ }
+
+ // Drawing Title
+ Rectangle2D.Double currentBBox = out.getBoundingBox();
+ double titleFontSize = (2.0 * conf._titleFont.getSize());
+ out.setColor(conf._titleColor);
+ out.setFont(PSExport.FONT_HELVETICA, titleFontSize);
+ double yTitle = currentBBox.y - titleFontSize / 2.0;
+ if (!getName().equals("")) {
+ out.drawText((maxX - minX) / 2.0, yTitle, getName());
+ }
+
+ OutputStreamWriter fout;
+
+ try {
+ fout = new OutputStreamWriter(new FileOutputStream(path), "UTF-8");
+
+ fout.write(out.export());
+ fout.close();
+ } catch (IOException e) {
+ throw new ExceptionWritingForbidden(e.getMessage());
+ }
+ }
+
+ Point2D.Double buildCaptionPosition(ModeleBase mb, double heightEstimate,
+ VARNAConfig conf) {
+ double radius = 2.0;
+ if (isNumberDrawn(mb, conf._numPeriod)) {
+ radius += (conf._distNumbers + 1.0);
+ }
+ Point2D.Double center = mb.getCenter();
+ Point2D.Double p = mb.getCoords();
+ double realDistance = BASE_RADIUS * radius + heightEstimate;
+ return new Point2D.Double(center.getX() + (p.getX() - center.getX())
+ * ((p.distance(center) + realDistance) / p.distance(center)),
+ center.getY()
+ + (p.getY() - center.getY())
+ * ((p.distance(center) + realDistance) / p
+ .distance(center)));
+ }
+
+ public double getBPHeightIncrement() {
+ return this._bpHeightIncrement;
+ }
+
+ public void setBPHeightIncrement(double d) {
+ _bpHeightIncrement = d;
+ }
+
+ private void drawChemProbAnnotation(SecStrDrawingProducer out,
+ ChemProbAnnotation cpa, Point2D.Double anchor, double minX,
+ double minY) {
+ out.setColor(cpa.getColor());
+ Point2D.Double v = cpa.getDirVector();
+ Point2D.Double vn = cpa.getNormalVector();
+ Point2D.Double base = new Point2D.Double((anchor.x + CHEM_PROB_DIST
+ * v.x), (anchor.y + CHEM_PROB_DIST * v.y));
+ Point2D.Double edge = new Point2D.Double(
+ (base.x + CHEM_PROB_BASE_LENGTH * cpa.getIntensity() * v.x),
+ (base.y + CHEM_PROB_BASE_LENGTH * cpa.getIntensity() * v.y));
+ double thickness = CHEM_PROB_ARROW_THICKNESS * cpa.getIntensity();
+ switch (cpa.getType()) {
+ case ARROW: {
+ Point2D.Double arrowTip1 = new Point2D.Double(
+ (base.x + cpa.getIntensity()
+ * (CHEM_PROB_ARROW_WIDTH * vn.x + CHEM_PROB_ARROW_HEIGHT
+ * v.x)),
+ (base.y + cpa.getIntensity()
+ * (CHEM_PROB_ARROW_WIDTH * vn.y + CHEM_PROB_ARROW_HEIGHT
+ * v.y)));
+ Point2D.Double arrowTip2 = new Point2D.Double(
+ (base.x + cpa.getIntensity()
+ * (-CHEM_PROB_ARROW_WIDTH * vn.x + CHEM_PROB_ARROW_HEIGHT
+ * v.x)),
+ (base.y + cpa.getIntensity()
+ * (-CHEM_PROB_ARROW_WIDTH * vn.y + CHEM_PROB_ARROW_HEIGHT
+ * v.y)));
+ out.drawLine(base.x - minX, minY - base.y, edge.x - minX, minY
+ - edge.y, thickness);
+ out.drawLine(base.x - minX, minY - base.y, arrowTip1.x - minX, minY
+ - arrowTip1.y, thickness);
+ out.drawLine(base.x - minX, minY - base.y, arrowTip2.x - minX, minY
+ - arrowTip2.y, thickness);
+ }
+ break;
+ case PIN: {
+ Point2D.Double side1 = new Point2D.Double(
+ (edge.x - cpa.getIntensity()
+ * (CHEM_PROB_PIN_SEMIDIAG * v.x)),
+ (edge.y - cpa.getIntensity()
+ * (CHEM_PROB_PIN_SEMIDIAG * v.y)));
+ Point2D.Double side2 = new Point2D.Double(
+ (edge.x - cpa.getIntensity()
+ * (CHEM_PROB_PIN_SEMIDIAG * vn.x)),
+ (edge.y - cpa.getIntensity()
+ * (CHEM_PROB_PIN_SEMIDIAG * vn.y)));
+ Point2D.Double side3 = new Point2D.Double(
+ (edge.x + cpa.getIntensity()
+ * (CHEM_PROB_PIN_SEMIDIAG * v.x)),
+ (edge.y + cpa.getIntensity()
+ * (CHEM_PROB_PIN_SEMIDIAG * v.y)));
+ Point2D.Double side4 = new Point2D.Double(
+ (edge.x + cpa.getIntensity()
+ * (CHEM_PROB_PIN_SEMIDIAG * vn.x)),
+ (edge.y + cpa.getIntensity()
+ * (CHEM_PROB_PIN_SEMIDIAG * vn.y)));
+ GeneralPath p2 = new GeneralPath();
+ p2.moveTo((float) (side1.x - minX), (float) (minY - side1.y));
+ p2.lineTo((float) (side2.x - minX), (float) (minY - side2.y));
+ p2.lineTo((float) (side3.x - minX), (float) (minY - side3.y));
+ p2.lineTo((float) (side4.x - minX), (float) (minY - side4.y));
+ p2.closePath();
+ out.fillPolygon(p2, cpa.getColor());
+ out.drawLine(base.x - minX, minY - base.y, edge.x - minX, minY
+ - edge.y, thickness);
+ }
+ break;
+ case TRIANGLE: {
+ Point2D.Double arrowTip1 = new Point2D.Double(
+ (edge.x + cpa.getIntensity()
+ * (CHEM_PROB_TRIANGLE_WIDTH * vn.x)),
+ (edge.y + cpa.getIntensity()
+ * (CHEM_PROB_TRIANGLE_WIDTH * vn.y)));
+ Point2D.Double arrowTip2 = new Point2D.Double(
+ (edge.x + cpa.getIntensity()
+ * (-CHEM_PROB_TRIANGLE_WIDTH * vn.x)),
+ (edge.y + cpa.getIntensity()
+ * (-CHEM_PROB_TRIANGLE_WIDTH * vn.y)));
+ GeneralPath p2 = new GeneralPath();
+ p2.moveTo((float) (base.x - minX), (float) (minY - base.y));
+ p2.lineTo((float) (arrowTip1.x - minX),
+ (float) (minY - arrowTip1.y));
+ p2.lineTo((float) (arrowTip2.x - minX),
+ (float) (minY - arrowTip2.y));
+ p2.closePath();
+ out.fillPolygon(p2, cpa.getColor());
+ }
+ break;
+ case DOT: {
+ Double radius = CHEM_PROB_DOT_RADIUS * cpa.getIntensity();
+ Point2D.Double center = new Point2D.Double((base.x + radius * v.x)
+ - minX, minY - (base.y + radius * v.y));
+ out.fillCircle(center.x, center.y, radius, thickness,
+ cpa.getColor());
+ }
+ break;
+ }
+ }
+
+ private void renderAnnotations(SecStrDrawingProducer out, double minX,
+ double minY, VARNAConfig conf) {
+ for (TextAnnotation textAnnotation : getAnnotations()) {
+ out.setColor(textAnnotation.getColor());
+ out.setFont(PSExport.FONT_HELVETICA_BOLD, 2.0 * textAnnotation
+ .getFont().getSize());
+ Point2D.Double position = textAnnotation.getCenterPosition();
+ if (textAnnotation.getType() == TextAnnotation.AnchorType.BASE) {
+ ModeleBase mb = (ModeleBase) textAnnotation.getAncrage();
+ double fontHeight = Math.ceil(textAnnotation.getFont()
+ .getSize());
+ position = buildCaptionPosition(mb, fontHeight, conf);
+
+ }
+ out.drawText(position.x - minX, -(position.y - minY),
+ textAnnotation.getTexte());
+ }
+ for (ChemProbAnnotation cpa : getChemProbAnnotations()) {
+ Point2D.Double anchor = cpa.getAnchorPosition();
+ drawChemProbAnnotation(out, cpa, anchor, minX, minY);
+ }
+ }
+
+ public boolean isNumberDrawn(ModeleBase mb, int numPeriod) {
+ if (numPeriod <= 0)
+ return false;
+ return ((mb.getIndex() == 0) || ((mb.getBaseNumber()) % numPeriod == 0) || (mb
+ .getIndex() == get_listeBases().size() - 1));
+ }
+
+ public void saveRNAEPS(String path, VARNAConfig conf)
+ throws ExceptionWritingForbidden {
+ PSExport out = new PSExport();
+ saveRNA(path, conf, 0.4, out);
+ }
+
+ public void saveRNAXFIG(String path, VARNAConfig conf)
+ throws ExceptionWritingForbidden {
+ XFIGExport out = new XFIGExport();
+ saveRNA(path, conf, 20, out);
+ }
+
+ public void saveRNATIKZ(String path, VARNAConfig conf)
+ throws ExceptionWritingForbidden {
+ TikzExport out = new TikzExport();
+ saveRNA(path, conf, 0.15, out);
+ }
+
+ public void saveRNASVG(String path, VARNAConfig conf)
+ throws ExceptionWritingForbidden {
+ SVGExport out = new SVGExport();
+ saveRNA(path, conf, 0.5, out);
+ }
+
+ public Rectangle2D.Double getBBox() {
+ Rectangle2D.Double result = new Rectangle2D.Double(10, 10, 10, 10);
+ double minx, maxx, miny, maxy;
+ minx = Double.MAX_VALUE;
+ miny = Double.MAX_VALUE;
+ maxx = -Double.MAX_VALUE;
+ maxy = -Double.MAX_VALUE;
+ for (int i = 0; i < _listeBases.size(); i++) {
+ minx = Math.min(
+ _listeBases.get(i).getCoords().getX() - BASE_RADIUS, minx);
+ miny = Math.min(
+ _listeBases.get(i).getCoords().getY() - BASE_RADIUS, miny);
+ maxx = Math.max(
+ _listeBases.get(i).getCoords().getX() + BASE_RADIUS, maxx);
+ maxy = Math.max(
+ _listeBases.get(i).getCoords().getY() + BASE_RADIUS, maxy);
+ }
+ result.x = minx;
+ result.y = miny;
+ result.width = Math.max(maxx - minx, 1);
+ result.height = Math.max(maxy - miny, 1);
+ if (_drawMode == RNA.DRAW_MODE_LINEAR) {
+ double realHeight = _bpHeightIncrement * result.width / 2.0;
+ result.height += realHeight;
+ result.y -= realHeight;
+ }
+ return result;
+ }
+
+ public void setCoord(int index, Point2D.Double p) {
+ setCoord(index, p.x, p.y);
+ }
+
+ public void setCoord(int index, double x, double y) {
+ if (index < _listeBases.size()) {
+ _listeBases.get(index).setCoords(new Point2D.Double(x, y));
+ }
+ }
+
+ public Point2D.Double getCoords(int i) {
+ if (i < _listeBases.size() && i >= 0) {
+ return _listeBases.get(i).getCoords();
+ }
+ return new Point2D.Double();
+ }
+
+ public String getBaseContent(int i) {
+ if ((i >= 0) && (i < _listeBases.size())) {
+ return _listeBases.get(i).getContent();
+ }
+ return "";
+ }
+
+ public int getBaseNumber(int i) {
+ if ((i >= 0) && (i < _listeBases.size())) {
+ return _listeBases.get(i).getBaseNumber();
+ }
+ return -1;
+ }
+
+ public Point2D.Double getCenter(int i) {
+ if (i < _listeBases.size()) {
+ return _listeBases.get(i).getCenter();
+ }
+
+ return new Point2D.Double();
+ }
+
+ public void setCenter(int i, double x, double y) {
+ setCenter(i, new Point2D.Double(x, y));
+ }
+
+ public void setCenter(int i, Point2D.Double p) {
+ if (i < _listeBases.size()) {
+ _listeBases.get(i).setCenter(p);
+ }
+ }
+
+ public void drawRNACircle(VARNAConfig conf) {
+ _drawn = true;
+ _drawMode = DRAW_MODE_CIRCULAR;
+ int radius = (int) ((3 * (_listeBases.size() + 1) * BASE_RADIUS) / (2 * Math.PI));
+ double angle;
+ for (int i = 0; i < _listeBases.size(); i++) {
+ angle = -((((double) -(i + 1)) * 2.0 * Math.PI)
+ / ((double) (_listeBases.size() + 1)) - Math.PI / 2.0);
+ _listeBases
+ .get(i)
+ .setCoords(
+ new Point2D.Double(
+ (radius * Math.cos(angle) * conf._spaceBetweenBases),
+ (radius * Math.sin(angle) * conf._spaceBetweenBases)));
+ _listeBases.get(i).setCenter(new Point2D.Double(0, 0));
+ }
+ }
+
+ public void drawRNAVARNAView(VARNAConfig conf) {
+ _drawn = true;
+ _drawMode = DRAW_MODE_VARNA_VIEW;
+ VARNASecDraw vs = new VARNASecDraw();
+ vs.drawRNA(1, this);
+ }
+
+ public void drawRNALine(VARNAConfig conf) {
+ _drawn = true;
+ _drawMode = DRAW_MODE_LINEAR;
+ for (int i = 0; i < get_listeBases().size(); i++) {
+ get_listeBases().get(i).setCoords(
+ new Point2D.Double(i * conf._spaceBetweenBases * 20, 0));
+ get_listeBases().get(i).setCenter(
+ new Point2D.Double(i * conf._spaceBetweenBases * 20, -10));
+ }
+ }
+
+ public RNATemplateMapping drawRNATemplate(RNATemplate template, boolean straightBulges,
+ VARNAConfig conf) throws RNATemplateDrawingAlgorithmException {
+ return drawRNATemplate(template, conf,
+ DrawRNATemplateMethod.getDefault(),
+ DrawRNATemplateCurveMethod.getDefault(), straightBulges);
+ }
+
+ public RNATemplateMapping drawRNATemplate(RNATemplate template,
+ VARNAConfig conf, DrawRNATemplateMethod helixLengthAdjustmentMethod,
+ boolean straightBulges)
+ throws RNATemplateDrawingAlgorithmException {
+ return drawRNATemplate(template, conf, helixLengthAdjustmentMethod,
+ DrawRNATemplateCurveMethod.getDefault(),straightBulges);
+ }
+
+ public RNATemplateMapping drawRNATemplate(RNATemplate template,
+ VARNAConfig conf,
+ DrawRNATemplateMethod helixLengthAdjustmentMethod,
+ DrawRNATemplateCurveMethod curveMethod,
+ boolean straightBulges)
+ throws RNATemplateDrawingAlgorithmException {
+ _drawn = true;
+ _drawMode = DRAW_MODE_TEMPLATE;
+
+ DrawRNATemplate drawRNATemplate = new DrawRNATemplate(this);
+ drawRNATemplate.drawRNATemplate(template, conf,
+ helixLengthAdjustmentMethod, curveMethod,
+ straightBulges);
+ return drawRNATemplate.getMapping();
+ }
+
+ private static double objFun(int n1, int n2, double r, double bpdist,
+ double multidist) {
+ return (((double) n1) * 2.0 * Math.asin(((double) bpdist) / (2.0 * r))
+ + ((double) n2) * 2.0
+ * Math.asin(((double) multidist) / (2.0 * r)) - (2.0 * Math.PI));
+ }
+
+ public double determineRadius(int nbHel, int nbUnpaired, double startRadius) {
+ return determineRadius(nbHel, nbUnpaired, startRadius,
+ BASE_PAIR_DISTANCE, MULTILOOP_DISTANCE);
+ }
+
+ public static double determineRadius(int nbHel, int nbUnpaired,
+ double startRadius, double bpdist, double multidist) {
+ double xmin = bpdist / 2.0;
+ double xmax = 3.0 * multidist + 1;
+ double x = (xmin + xmax) / 2.0;
+ double y = 10000.0;
+ double ymin = -1000.0;
+ double ymax = 1000.0;
+ int numIt = 0;
+ double precision = 0.00001;
+ while ((Math.abs(y) > precision) && (numIt < 10000)) {
+ x = (xmin + xmax) / 2.0;
+ y = objFun(nbHel, nbUnpaired, x, bpdist, multidist);
+ ymin = objFun(nbHel, nbUnpaired, xmax, bpdist, multidist);
+ ymax = objFun(nbHel, nbUnpaired, xmin, bpdist, multidist);
+ if (ymin > 0.0) {
+ xmax = xmax + (xmax - xmin);
+ } else if ((y <= 0.0) && (ymax > 0.0)) {
+ xmax = x;
+ } else if ((y >= 0.0) && (ymin < 0.0)) {
+ xmin = x;
+ } else if (ymax < 0.0) {
+ xmin = Math.max(xmin - (x - xmin),
+ Math.max(bpdist / 2.0, multidist / 2.0));
+ xmax = x;
+ }
+ numIt++;
+ }
+ return x;
+ }
+
+ public void drawRNA(VARNAConfig conf) throws ExceptionNAViewAlgorithm {
+ drawRNA(RNA.DEFAULT_DRAW_MODE, conf);
+ }
+
+ public void drawRNA(int mode, VARNAConfig conf)
+ throws ExceptionNAViewAlgorithm {
+ _drawMode = mode;
+ switch (get_drawMode()) {
+ case RNA.DRAW_MODE_RADIATE:
+ drawRNARadiate(conf);
+ break;
+ case RNA.DRAW_MODE_LINEAR:
+ drawRNALine(conf);
+ break;
+ case RNA.DRAW_MODE_CIRCULAR:
+ drawRNACircle(conf);
+ break;
+ case RNA.DRAW_MODE_NAVIEW:
+ drawRNANAView(conf);
+ break;
+ case RNA.DRAW_MODE_VARNA_VIEW:
+ drawRNAVARNAView(conf);
+ break;
+ default:
+ break;
+ }
+
+ }
+
+ public int getDrawMode() {
+ return _drawMode;
+ }
+
+ public static double HYSTERESIS_EPSILON = .15;
+ public static final double[] HYSTERESIS_ATTRACTORS = {0.,Math.PI/4.,Math.PI/2.,3.*Math.PI/4.,Math.PI,5.*(Math.PI)/4.,3.*(Math.PI)/2,7.*(Math.PI)/4.};
+
+ public static double normalizeAngle(double angle)
+ {
+ return normalizeAngle(angle,0.);
+ }
+
+ public static double normalizeAngle(double angle, double fromVal)
+ {
+ double toVal = fromVal +2.*Math.PI;
+ double result = angle;
+ while(result<fromVal)
+ {
+ result += 2.*Math.PI;
+ }
+ while(result >= toVal)
+ {
+ result -= 2.*Math.PI;
+ }
+ return result;
+ }
+
+ public static double correctHysteresis(double angle)
+ {
+ double result = normalizeAngle(angle);
+
+ for (int i=0;i<HYSTERESIS_ATTRACTORS.length;i++)
+ {
+ double att = HYSTERESIS_ATTRACTORS[i];
+ if (Math.abs(normalizeAngle(att-result,-Math.PI))<HYSTERESIS_EPSILON)
+ {
+ result = att;
+ }
+ }
+ return result;
+ }
+
+
+
+ private void distributeUnpaired(
+ double radius,
+ double angle,
+ double pHel,
+ double base,
+ Point2D.Double center,
+ Vector<Integer> bases)
+ {
+ double mydist = Math.abs(radius*(angle / (bases.size() + 1)));
+ double addedRadius= 0.;
+ Point2D.Double PA = new Point2D.Double(center.x + radius * Math.cos(base + pHel),
+ center.y + radius * Math.sin(base + pHel));
+ Point2D.Double PB = new Point2D.Double(center.x + radius * Math.cos(base + pHel+angle),
+ center.y + radius * Math.sin(base + pHel+angle));
+ double dist = PA.distance(PB);
+ Point2D.Double VN = new Point2D.Double((PB.y-PA.y)/dist,(-PB.x+PA.x)/dist);
+ if (mydist<2*BASE_RADIUS)
+ {
+ addedRadius=Math.min(1.0,(2*BASE_RADIUS-mydist)/4)*computeRadius(mydist, 2.29*(bases.size() + 1)*BASE_RADIUS-mydist);
+ }
+
+
+ ArrayList<Point2D.Double> pos = computeNewAngles(bases.size(),center,VN, angle,base + pHel,radius,addedRadius);
+ for (int i = 0; i < bases.size(); i++)
+ {
+ int k = bases.get(i);
+ setCoord(k, pos.get(i));
+ }
+
+ }
+
+ private double computeRadius(double b, double pobj)
+ {
+ double a=b, aL=a, aU=Double.POSITIVE_INFINITY;
+ double h = (a-b)*(a-b)/((a+b)*(a+b));
+ double p = Math.PI*(a+b)*(1+h/4.+h*h/64.+h*h*h/256.+25.*h*h*h*h/16384.)/2.0;
+ double aold = a+1.;
+ while ((Math.abs(p-pobj)>10e-4)&&(aold!=a)){
+ aold = a;
+ if (p<pobj)
+ {
+ aL = a;
+ if (aU==Double.POSITIVE_INFINITY)
+ {a *= 2.;}
+ else
+ { a = (a+aU)/2.; }
+ }
+ else
+ {
+ aU = a;
+ a = (a+aL)/2.0;
+ }
+ h = (a-b)*(a-b)/((a+b)*(a+b));
+ p = (Math.PI*(a+b)*(1+h/4.+h*h/64.+h*h*h/256.+25.*h*h*h*h/16384.))/2.0;
+ }
+ return a;
+ }
+
+ public static double computeAngle(Point2D.Double center, Point2D.Double p) {
+ double dist = center.distance(p);
+ double angle = Math.asin((p.y - center.y) / dist);
+ if (p.x - center.x < 0) {
+ angle = Math.PI - angle;
+ }
+ return angle;
+ }
+
+ private Point2D.Double rotatePoint(Point2D.Double center, Point2D.Double p,
+ double angle) {
+ double dist = p.distance(center);
+ double oldAngle = Math.asin((p.y - center.y) / dist);
+
+ if (p.x - center.x < 0) {
+ oldAngle = Math.PI - oldAngle;
+ }
+
+ double newX = (center.x + dist * Math.cos(oldAngle + angle));
+ double newY = (center.y + dist * Math.sin(oldAngle + angle));
+
+ return new Point2D.Double(newX, newY);
+ }
+
+
+
+ private void rotateHelix(Point2D.Double center, int i, int j, double angle) {
+ for (int k = i; k <= j; k++) {
+ Point2D.Double oldp = getCoords(k);
+ Point2D.Double newp = rotatePoint(center, oldp, angle);
+ setCoord(k, newp);
+ if ((k != i) && (k != j)) {
+
+ Point2D.Double oldc = get_listeBases().get(k)
+ .getCenter();
+ Point2D.Double newc = rotatePoint(center, oldc, angle);
+ setCenter(k,newc);
+ }
+ }
+ }
+
+
+
+ private void fixUnpairedPositions(boolean isDirect,
+ double angleRightPartner,
+ double angleLimitLeft, double angleLimitRight, double angleLeftPartner,
+ double radius, double base, Point2D.Double center,Vector<Integer> prevBases,Vector<Integer> nextBases)
+ {
+ if (isDirect) {
+ double anglePrev = normalizeAngle(angleLimitLeft - angleRightPartner);
+ double angleNext = normalizeAngle(angleLeftPartner - angleLimitRight);
+ distributeUnpaired(radius,anglePrev, angleRightPartner, base,
+ center,prevBases);
+ distributeUnpaired(radius,-angleNext, angleLeftPartner, base,
+ center,nextBases);
+ } else {
+ double anglePrev = normalizeAngle(angleLeftPartner - angleLimitRight);
+ double angleNext = normalizeAngle(angleLimitLeft - angleRightPartner);
+ distributeUnpaired(radius,-anglePrev, angleLeftPartner, base,
+ center,prevBases);
+ distributeUnpaired(radius,angleNext, angleRightPartner, base,
+ center,nextBases);
+ }
+
+ }
+
+ private static Point2D.Double getPoint(double angleLine, double angleBulge, Point2D.Double center,
+ Point2D.Double VN,double radius, double addedRadius, double dirBulge)
+ {
+ return new Point2D.Double(
+ center.x + radius * Math.cos(angleLine)+
+ dirBulge*addedRadius*Math.sin(angleBulge)*VN.x,
+ center.y + radius * Math.sin(angleLine)+
+ dirBulge*addedRadius*Math.sin(angleBulge)*VN.y);
+ }
+
+
+ private ArrayList<Point2D.Double> computeNewAngles(int numPoints, Point2D.Double center,
+ Point2D.Double VN, double angle, double angleBase, double radius, double addedRadius)
+ {
+ ArrayList<Point2D.Double> result = new ArrayList<Point2D.Double>();
+ if (numPoints>0)
+ {
+ ArrayList<Double> factors = new ArrayList<Double>();
+
+
+ Point2D.Double prevP = new Point2D.Double(
+ center.x + radius * Math.cos(angleBase),
+ center.y + radius * Math.sin(angleBase));
+
+
+ double fact = 0.;
+
+ double angleBulge = 0.;
+ double dirBulge = (angle<0)?-1.:1.;
+ double dtarget =2.*BASE_RADIUS;
+
+ for (int i = 0; i < numPoints; i++)
+ {
+ double lbound = fact;
+ double ubound = 1.0;
+ double angleLine = angleBase + angle*fact;
+ angleBulge = Math.PI*fact;
+ Point2D.Double currP = getPoint(angleLine, angleBulge, center,VN,radius, addedRadius, dirBulge);
+
+ int numIter = 0;
+ while ((Math.abs(currP.distance(prevP)-dtarget)>0.01)&& (numIter<100))
+ {
+ if (currP.distance(prevP)> dtarget)
+ {
+ ubound = fact;
+ fact = (fact+lbound)/2.0;
+ }
+ else
+ {
+ lbound = fact;
+ fact = (fact+ubound)/2.0;
+ }
+ angleLine = angleBase + angle*fact;
+ angleBulge = Math.PI*fact;
+ currP = getPoint(angleLine, angleBulge, center,VN,radius, addedRadius, dirBulge);
+ numIter++;
+ }
+ factors.add(fact);
+ prevP = currP;
+ }
+
+
+ double rescale = 1.0/(factors.get(factors.size()-1)+factors.get(0));
+
+ for(int j=0;j<factors.size();j++)
+ {
+ factors.set(j,factors.get(j)*rescale);
+ }
+
+
+ if (addedRadius>0)
+ {
+ prevP = getPoint(angleBase, 0, center,VN,radius, addedRadius, dirBulge);
+ double totDist = 0.0;
+ for(int j=0;j<factors.size();j++)
+ {
+ double newfact = factors.get(j);
+ double angleLine = angleBase + angle*newfact;
+ angleBulge = Math.PI*newfact;
+ Point2D.Double currP = getPoint(angleLine, angleBulge, center,VN,radius, addedRadius, dirBulge);
+ totDist += currP.distance(prevP);
+ prevP = currP;
+ }
+ totDist += getPoint(angleBase+angle, Math.PI, center,VN,radius, addedRadius, dirBulge).distance(prevP);
+ dtarget = totDist/(numPoints+1);
+ fact = 0.0;
+ factors=new ArrayList<Double>();
+ prevP = new Point2D.Double(
+ center.x + radius * Math.cos(angleBase),
+ center.y + radius * Math.sin(angleBase));
+ for (int i = 0; i < numPoints; i++)
+ {
+ double lbound = fact;
+ double ubound = 1.5;
+ double angleLine = angleBase + angle*fact;
+ angleBulge = Math.PI*fact;
+ Point2D.Double currP = getPoint(angleLine, angleBulge, center,VN,radius, addedRadius, dirBulge);
+
+ int numIter = 0;
+ while ((Math.abs(currP.distance(prevP)-dtarget)>0.01)&& (numIter<100))
+ {
+ if (currP.distance(prevP)> dtarget)
+ {
+ ubound = fact;
+ fact = (fact+lbound)/2.0;
+ }
+ else
+ {
+ lbound = fact;
+ fact = (fact+ubound)/2.0;
+ }
+ angleLine = angleBase + angle*fact;
+ angleBulge = Math.PI*fact;
+ currP = getPoint(angleLine, angleBulge, center,VN,radius, addedRadius, dirBulge);
+ numIter++;
+ }
+ factors.add(fact);
+ prevP = currP;
+ }
+ rescale = 1.0/(factors.get(factors.size()-1)+factors.get(0));
+ for(int j=0;j<factors.size();j++)
+ {
+ factors.set(j,factors.get(j)*rescale);
+ }
+ }
+
+ for(int j=0;j<factors.size();j++)
+ {
+ double newfact = factors.get(j);
+ double angleLine = angleBase + angle*newfact;
+ angleBulge = Math.PI*newfact;
+ result.add(getPoint(angleLine, angleBulge, center,VN,radius, addedRadius, dirBulge));
+ }
+ }
+ return result;
+ }
+
+
+
+ void drawLoop(int i, int j, double x, double y, double dirAngle,
+ Point2D.Double[] coords, Point2D.Double[] centers, double[] angles,
+ boolean straightBulges) {
+ if (i > j) {
+ return;
+ }
+
+ // BasePaired
+ if (_listeBases.get(i).getElementStructure() == j) {
+ double normalAngle = Math.PI / 2.0;
+ centers[i] = new Point2D.Double(x, y);
+ centers[j] = new Point2D.Double(x, y);
+ coords[i].x = (x + BASE_PAIR_DISTANCE
+ * Math.cos(dirAngle - normalAngle) / 2.0);
+ coords[i].y = (y + BASE_PAIR_DISTANCE
+ * Math.sin(dirAngle - normalAngle) / 2.0);
+ coords[j].x = (x + BASE_PAIR_DISTANCE
+ * Math.cos(dirAngle + normalAngle) / 2.0);
+ coords[j].y = (y + BASE_PAIR_DISTANCE
+ * Math.sin(dirAngle + normalAngle) / 2.0);
+ drawLoop(i + 1, j - 1, x + LOOP_DISTANCE * Math.cos(dirAngle), y
+ + LOOP_DISTANCE * Math.sin(dirAngle), dirAngle, coords,
+ centers, angles, straightBulges);
+ } else {
+ int k = i;
+ Vector<Integer> basesMultiLoop = new Vector<Integer>();
+ Vector<Integer> helices = new Vector<Integer>();
+ int l;
+ while (k <= j) {
+ l = _listeBases.get(k).getElementStructure();
+ if (l > k) {
+ basesMultiLoop.add(new Integer(k));
+ basesMultiLoop.add(new Integer(l));
+ helices.add(new Integer(k));
+ k = l + 1;
+ } else {
+ basesMultiLoop.add(new Integer(k));
+ k++;
+ }
+ }
+ int mlSize = basesMultiLoop.size() + 2;
+ int numHelices = helices.size() + 1;
+ double totalLength = MULTILOOP_DISTANCE * (mlSize - numHelices)
+ + BASE_PAIR_DISTANCE * numHelices;
+ double multiLoopRadius;
+ double angleIncrementML;
+ double angleIncrementBP;
+ if (mlSize > 3) {
+ multiLoopRadius = determineRadius(numHelices, mlSize
+ - numHelices, (totalLength) / (2.0 * Math.PI),
+ BASE_PAIR_DISTANCE, MULTILOOP_DISTANCE);
+ angleIncrementML = -2.0
+ * Math.asin(((float) MULTILOOP_DISTANCE)
+ / (2.0 * multiLoopRadius));
+ angleIncrementBP = -2.0
+ * Math.asin(((float) BASE_PAIR_DISTANCE)
+ / (2.0 * multiLoopRadius));
+ }
+ else {
+ multiLoopRadius = 35.0;
+ angleIncrementBP = -2.0
+ * Math.asin(((float) BASE_PAIR_DISTANCE)
+ / (2.0 * multiLoopRadius));
+ angleIncrementML = (-2.0 * Math.PI - angleIncrementBP) / 2.0;
+ }
+ // System.out.println("MLr:"+multiLoopRadius+" iBP:"+angleIncrementBP+" iML:"+angleIncrementML);
+
+ double centerDist = Math.sqrt(Math.max(Math.pow(multiLoopRadius, 2)
+ - Math.pow(BASE_PAIR_DISTANCE / 2.0, 2), 0.0))
+ - LOOP_DISTANCE;
+ Point2D.Double mlCenter = new Point2D.Double(
+ (x + (centerDist * Math.cos(dirAngle))),
+ (y + (centerDist * Math.sin(dirAngle))));
+
+ // Base directing angle for (multi|hairpin) loop, from the center's
+ // perspective
+ double baseAngle = dirAngle
+ // U-turn
+ + Math.PI
+ // Account for already drawn supporting base-pair
+ + 0.5 * angleIncrementBP
+ // Base cannot be paired twice, so next base is at
+ // "unpaired base distance"
+ + 1.0 * angleIncrementML;
+
+ ArrayList<Integer> currUnpaired = new ArrayList<Integer>();
+ Couple<Double,Double> currInterval = new Couple<Double,Double>(0.,baseAngle-1.0 * angleIncrementML);
+ ArrayList<Couple<ArrayList<Integer>,Couple<Double,Double>>> intervals = new ArrayList<Couple<ArrayList<Integer>,Couple<Double,Double>>>();
+
+ for (k = basesMultiLoop.size() - 1; k >= 0; k--) {
+ l = basesMultiLoop.get(k).intValue();
+ //System.out.println(l+" ");
+ centers[l] = mlCenter;
+ boolean isPaired = (_listeBases.get(l).getElementStructure() != -1);
+ boolean isPaired3 = isPaired && (_listeBases.get(l).getElementStructure() < l);
+ boolean isPaired5 = isPaired && !isPaired3;
+ if (isPaired3) {
+ if ((numHelices == 2) && straightBulges)
+ {
+ baseAngle = dirAngle-angleIncrementBP/2.;
+ }
+ else
+ {
+ baseAngle = correctHysteresis(baseAngle+angleIncrementBP/2.)-angleIncrementBP/2.;
+ }
+ currInterval.first = baseAngle;
+ intervals.add(new Couple<ArrayList<Integer>,Couple<Double,Double>>(currUnpaired,currInterval));
+ currInterval = new Couple<Double,Double>(-1.,-1.);
+ currUnpaired = new ArrayList<Integer>();
+ }
+ else if (isPaired5)
+ {
+ currInterval.second = baseAngle;
+ }
+ else
+ {
+ currUnpaired.add(l);
+ }
+
+ angles[l] = baseAngle;
+ if (isPaired3)
+ {
+ baseAngle += angleIncrementBP;
+ }
+ else {
+ baseAngle += angleIncrementML;
+ }
+ }
+ currInterval.first = dirAngle
+ - Math.PI
+ - 0.5 * angleIncrementBP;
+ intervals.add(new Couple<ArrayList<Integer>,Couple<Double,Double>>(currUnpaired,currInterval));
+ //System.out.println("Inc. ML:"+angleIncrementML+" BP:"+angleIncrementBP);
+
+ for(Couple<ArrayList<Integer>,Couple<Double,Double>> inter: intervals)
+ {
+ //double mid = inter.second.second;
+ double mina = inter.second.first;
+ double maxa = normalizeAngle(inter.second.second,mina);
+ //System.out.println(""+mina+" " +maxa);
+
+ for (int n=0;n<inter.first.size();n++)
+ {
+ double ratio = (1.+n)/(1.+inter.first.size());
+ int b = inter.first.get(n);
+ angles[b] = mina + (1.-ratio)*(maxa-mina);
+ }
+ }
+
+
+ for (k = basesMultiLoop.size() - 1; k >= 0; k--) {
+ l = basesMultiLoop.get(k).intValue();
+ coords[l].x = mlCenter.x + multiLoopRadius
+ * Math.cos(angles[l]);
+ coords[l].y = mlCenter.y + multiLoopRadius
+ * Math.sin(angles[l]);
+ }
+
+ // System.out.println("n1:"+n1+" n2:"+n2);
+ double newAngle;
+ int m, n;
+ for (k = 0; k < helices.size(); k++) {
+ m = helices.get(k).intValue();
+ n = _listeBases.get(m).getElementStructure();
+ newAngle = (angles[m] + angles[n]) / 2.0;
+ drawLoop(m + 1, n - 1, (LOOP_DISTANCE * Math.cos(newAngle))
+ + (coords[m].x + coords[n].x) / 2.0,
+ (LOOP_DISTANCE * Math.sin(newAngle))
+ + (coords[m].y + coords[n].y) / 2.0, newAngle,
+ coords, centers, angles, straightBulges);
+ }
+ }
+ }
+
+ private Vector<Integer> getPreviousUnpaired(Point h)
+ {
+ Vector<Integer> prevBases = new Vector<Integer>();
+ boolean over = false;
+ int i = h.y + 1;
+ while (!over) {
+ if (i >=get_listeBases().size()) {
+ over = true;
+ } else {
+ if (get_listeBases().get(i)
+ .getElementStructure() == -1) {
+ prevBases.add(new Integer(i));
+ } else {
+ over = true;
+ }
+ }
+ i++;
+ }
+ return prevBases;
+ }
+
+ private Vector<Integer> getNextUnpaired(Point h)
+ {
+ boolean over = false;
+ int i = h.x - 1;
+ Vector<Integer> nextBases = new Vector<Integer>();
+ while (!over) {
+ if (i < 0) {
+ over = true;
+ } else {
+ if (get_listeBases().get(i)
+ .getElementStructure() == -1) {
+ nextBases.add(new Integer(i));
+ } else {
+ over = true;
+ }
+ }
+ i--;
+ }
+ return nextBases;
+ }
+
+
+ public void rotateEverything(double delta, double base, double pLimL, double pLimR, Point h, Point ml, Hashtable<Integer,Point2D.Double> backupPos)
+ {
+ boolean isDirect = testDirectionality(ml.x, ml.y, h.x);
+ Point2D.Double center = get_listeBases().get(h.x).getCenter();
+ for(int k=h.x;k<=h.y;k++)
+ { backupPos.put(k, getBaseAt(k).getCoords()); }
+ rotateHelix(center, h.x, h.y, delta);
+
+ // Re-assigns unpaired atoms
+ Point2D.Double helixStart = getCoords(h.x);
+ Point2D.Double helixStop = getCoords(h.y);
+ double pHelR,pHelL;
+ if (isDirect) {
+ pHelR = computeAngle(center, helixStop) - base;
+ pHelL = computeAngle(center, helixStart) - base;
+ } else {
+ pHelL = computeAngle(center, helixStop) - base;
+ pHelR = computeAngle(center, helixStart) - base;
+ }
+
+ Vector<Integer> prevBases = getPreviousUnpaired(h);
+ Vector<Integer> nextBases = getNextUnpaired(h);
+
+ double radius = center.distance(helixStart);
+
+ for (int j = 0; j < prevBases.size(); j++)
+ {
+ int k = prevBases.get(j);
+ backupPos.put(k, getCoords(k));
+ }
+ for (int j = 0; j < nextBases.size(); j++)
+ {
+ int k = nextBases.get(j);
+ backupPos.put(k, getCoords(k));
+ }
+ fixUnpairedPositions(isDirect, pHelR, pLimL, pLimR, pHelL, radius, base,center,prevBases,nextBases);
+ }
+
+
+
+ public void drawRNARadiate() {
+ drawRNARadiate(-1.0, VARNAConfig.DEFAULT_SPACE_BETWEEN_BASES, true, true);
+ }
+
+ public void drawRNARadiate(VARNAConfig conf) {
+ drawRNARadiate(-1.0, conf._spaceBetweenBases, conf._flatExteriorLoop, false);
+ }
+
+ public static final double FLAT_RECURSIVE_INCREMENT = 20.;
+
+ public void drawRNARadiate(double dirAngle, double _spaceBetweenBases,
+ boolean flatExteriorLoop, boolean straightBulges) {
+ _drawn = true;
+ straightBulges = true;
+ _drawMode = DRAW_MODE_RADIATE;
+ Point2D.Double[] coords = new Point2D.Double[_listeBases.size()];
+ Point2D.Double[] centers = new Point2D.Double[_listeBases.size()];
+ double[] angles = new double[_listeBases.size()];
+ for (int i = 0; i < _listeBases.size(); i++) {
+ coords[i] = new Point2D.Double(0, 0);
+ centers[i] = new Point2D.Double(0, 0);
+ }
+ if (flatExteriorLoop) {
+ dirAngle += 1.0 - Math.PI / 2.0;
+ int i = 0;
+ double x = 0.0;
+ double y = 0.0;
+ double vx = -Math.sin(dirAngle);
+ double vy = Math.cos(dirAngle);
+ while (i < _listeBases.size()) {
+ coords[i].x = x;
+ coords[i].y = y;
+ centers[i].x = x + BASE_PAIR_DISTANCE * vy;
+ centers[i].y = y - BASE_PAIR_DISTANCE * vx;
+ int j = _listeBases.get(i).getElementStructure();
+ if (j > i) {
+ double increment = 0.;
+ if (i+1<_listeBases.size())
+ {
+ if (_listeBases.get(i+1).getElementStructure()==-1)
+ {
+ //increment = -FLAT_RECURSIVE_INCREMENT;
+ }
+ }
+ drawLoop(i, j, x + (BASE_PAIR_DISTANCE * vx / 2.0), y
+ + (BASE_PAIR_DISTANCE * vy / 2.0)+increment, dirAngle,
+ coords, centers, angles, straightBulges);
+ centers[i].x = coords[i].x + BASE_PAIR_DISTANCE * vy;
+ centers[i].y = y - BASE_PAIR_DISTANCE * vx;
+ i = j;
+ x += BASE_PAIR_DISTANCE * vx;
+ y += BASE_PAIR_DISTANCE * vy;
+ centers[i].x = coords[i].x + BASE_PAIR_DISTANCE * vy;
+ centers[i].y = y - BASE_PAIR_DISTANCE * vx;
+ }
+ x += MULTILOOP_DISTANCE * vx;
+ y += MULTILOOP_DISTANCE * vy;
+ i += 1;
+ }
+ } else {
+ drawLoop(0, _listeBases.size() - 1, 0, 0, dirAngle, coords, centers, angles, straightBulges);
+ }
+ for (int i = 0; i < _listeBases.size(); i++) {
+ _listeBases.get(i).setCoords(
+ new Point2D.Double(coords[i].x * _spaceBetweenBases,
+ coords[i].y * _spaceBetweenBases));
+ _listeBases.get(i).setCenter(
+ new Point2D.Double(centers[i].x * _spaceBetweenBases,
+ centers[i].y * _spaceBetweenBases));
+ }
+
+ // TODO
+ // change les centres des bases de la premiere helice vers la boucle la
+ // plus proche
+ }
+
+ public void drawRNANAView(VARNAConfig conf) throws ExceptionNAViewAlgorithm {
+ _drawMode = DRAW_MODE_NAVIEW;
+ _drawn = true;
+
+ ArrayList<Double> X = new ArrayList<Double>(_listeBases.size());
+ ArrayList<Double> Y = new ArrayList<Double>(_listeBases.size());
+ ArrayList<Short> pair_table = new ArrayList<Short>(_listeBases.size());
+
+ for (int i = 0; i < _listeBases.size(); i++) {
+ pair_table.add(Short.valueOf(String.valueOf(_listeBases.get(i)
+ .getElementStructure())));
+ }
+ NAView naView = new NAView();
+ naView.naview_xy_coordinates(pair_table, X, Y);
+
+ // Updating individual base positions
+ for (int i = 0; i < _listeBases.size(); i++) {
+ _listeBases.get(i).setCoords(
+ new Point2D.Double(
+ X.get(i) * 2.5 * conf._spaceBetweenBases, Y.get(i)
+ * 2.5 * conf._spaceBetweenBases));
+ }
+
+ // Updating centers
+ for (int i = 0; i < _listeBases.size(); i++) {
+ int indicePartner = _listeBases.get(i).getElementStructure();
+ if (indicePartner != -1) {
+ Point2D.Double base = _listeBases.get(i).getCoords();
+ Point2D.Double partner = _listeBases.get(indicePartner)
+ .getCoords();
+ _listeBases.get(i).setCenter(
+ new Point2D.Double((base.x + partner.x) / 2.0,
+ (base.y + partner.y) / 2.0));
+ } else {
+ Vector<Integer> loop = getLoopBases(i);
+ double tmpx = 0.0;
+ double tmpy = 0.0;
+ for (int j = 0; j < loop.size(); j++) {
+ int partner = loop.elementAt(j);
+ Point2D.Double loopmember = _listeBases.get(partner)
+ .getCoords();
+ tmpx += loopmember.x;
+ tmpy += loopmember.y;
+ }
+ _listeBases.get(i).setCenter(
+ new Point2D.Double(tmpx / loop.size(), tmpy
+ / loop.size()));
+ }
+ }
+ }
+
+ /*
+ * public void drawMOTIFView() { _drawn = true; _drawMode =
+ * DRAW_MODE_MOTIFVIEW; int spaceBetweenStrand =0; Motif motif = new
+ * Motif(this,get_listeBases()); motif.listStrand(); for (int i = 0; i <
+ * motif.getListStrand().sizeStruct(); i++ ){ for (int j = 0; j <
+ * motif.getListStrand().getStrand(i).sizeStrand(); j++ ){ int indice =
+ * motif.getListStrand().getStrand(i).getMB(j).getIndex();
+ * get_listeBases().get(indice).setCoords( new Point2D.Double(0,0));
+ * get_listeBases().get(indice).setCenter( new Point2D.Double(0, 0));
+ *
+ * } } //Recherche du brin central int centralStrand =
+ * motif.getCentralStrand();
+ *
+ * //Cas o? l'on a un motif en ?toile if(centralStrand!=-1){ //On positionne
+ * le brin central motif.positionneSpecificStrand(centralStrand,
+ * spaceBetweenStrand);
+ *
+ * //On place les autres brins par rapport a ce brin central
+ * motif.orderStrands(centralStrand); }
+ *
+ * else { centralStrand = 0; motif.positionneStrand(); motif.ajusteStrand();
+ * } motif.reajustement(); motif.deviationBasePair();
+ * motif.setCenterMotif(); }
+ */
+
+ public ArrayList<ModeleBase> getAllPartners(int indice) {
+ ArrayList<ModeleBase> result = new ArrayList<ModeleBase>();
+ ModeleBase me = this.getBaseAt(indice);
+ int i = me.getElementStructure();
+ if (i != -1) {
+ result.add(getBaseAt(i));
+ }
+ ArrayList<ModeleBP> msbps = getAuxBPs(indice);
+ for (ModeleBP m : msbps) {
+ result.add(m.getPartner(me));
+ }
+ return result;
+ }
+
+ public int get_drawMode() {
+ return _drawMode;
+ }
+
+ public void setDrawMode(int drawMode) {
+ _drawMode = drawMode;
+ }
+
+ public Set<Integer> getSeparatorPositions(String s) {
+ HashSet<Integer> result = new HashSet<Integer>();
+ int index = s.indexOf(DBNStrandSep);
+ while (index >= 0) {
+ result.add(index);
+ index = s.indexOf(DBNStrandSep, index + 1);
+ }
+ return result;
+ }
+
+ public static String DBNStrandSep = "&";
+
+ public void setRNA(String seq, String str)
+ throws ExceptionFileFormatOrSyntax,
+ ExceptionUnmatchedClosingParentheses {
+ ArrayList<String> al = RNA.explodeSequence(seq);
+ Set<Integer> sepPos = getSeparatorPositions(str);
+ ArrayList<String> alRes = new ArrayList<String>();
+ Set<Integer> resSepPos = new HashSet<Integer>();
+ String strRes = "";
+ for (int i = 0; i < al.size(); i++) {
+ if (sepPos.contains(i) && al.get(i).equals(DBNStrandSep)) {
+ resSepPos.add(alRes.size() - 1);
+ } else {
+ alRes.add(al.get(i));
+ if (i<str.length())
+ {
+ strRes += str.charAt(i);
+ }
+ else
+ {
+ strRes += '.';
+ }
+ }
+ }
+ for (int i = al.size(); i < str.length(); i++) {
+ alRes.add(" ");
+ strRes += str.charAt(i);
+ }
+ setRNA(alRes, strRes);
+ for (int i : resSepPos) {
+ _backbone.addElement(new ModeleBackboneElement(i,
+ BackboneType.DISCONTINUOUS_TYPE));
+ }
+ }
+
+ public void setRNA(String seq) {
+ ArrayList<String> s = RNA.explodeSequence(seq);
+ int[] str = new int[s.size()];
+ for (int i = 0; i < str.length; i++) {
+ str[i] = -1;
+ }
+ try {
+ setRNA(s, str);
+ } catch (ExceptionFileFormatOrSyntax e) {
+ // TODO Auto-generated catch block
+ e.printStackTrace();
+ }
+ }
+
+ public void setRNA(String seq, int[] str)
+ throws ExceptionFileFormatOrSyntax,
+ ExceptionUnmatchedClosingParentheses {
+ setRNA(RNA.explodeSequence(seq), str);
+ }
+
+ public void setRNA(String[] seq, int[] str)
+ throws ExceptionFileFormatOrSyntax {
+ setRNA(seq, str, 1);
+ }
+
+ public void setRNA(List<String> seq, int[] str)
+ throws ExceptionFileFormatOrSyntax {
+ setRNA(seq.toArray(new String[seq.size()]), str, 1);
+ }
+
+ public void setRNA(List<String> seq, int[] str, int baseIndex)
+ throws ExceptionFileFormatOrSyntax {
+ setRNA(seq.toArray(new String[seq.size()]), str, baseIndex);
+ }
+
+ public void setRNA(String[] seq, int[] str, int baseIndex)
+ throws ExceptionFileFormatOrSyntax {
+ clearAnnotations();
+ _listeBases = new ArrayList<ModeleBase>();
+ if (seq.length != str.length) {
+ warningEmition("Sequence length " + seq.length
+ + " differs from that of secondary structure " + str.length
+ + ". \nAdapting sequence length ...");
+ if (seq.length < str.length) {
+ String[] nseq = new String[str.length];
+ for (int i = 0; i < seq.length; i++) {
+ nseq[i] = seq[i];
+ }
+ for (int i = seq.length; i < nseq.length; i++) {
+ nseq[i] = "";
+ }
+ seq = nseq;
+ } else {
+ String[] seqTmp = new String[str.length];
+ for (int i = 0; i < str.length; i++) {
+ seqTmp[i] = seq[i];
+ }
+ seq = seqTmp;
+ }
+ }
+ for (int i = 0; i < str.length; i++) {
+ _listeBases.add(new ModeleBaseNucleotide(seq[i], i, baseIndex + i));
+ }
+ applyStruct(str);
+ }
+
+ /**
+ * Sets the RNA to be drawn. Uses when comparison mode is on. Will draw the
+ * super-structure passed in parameters and apply specials styles to the
+ * bases owning by each RNA alignment and both.
+ *
+ * @param seq
+ * - The sequence of the super-structure This sequence shall be
+ * designed like this:
+ * <code>firstRNA1stBaseSecondRNA1stBaseFirstRNA2ndBaseSecondRNA2ndBase [...]</code>
+ * <br>
+ * <b>Example:</b> <code>AAC-GUAGA--UGG</code>
+ * @param struct
+ * - The super-structure
+ * @param basesOwn
+ * - The RNA owning bases array (each index will be:0 when common
+ * base, 1 when first RNA alignment base, 2 when second RNA
+ * alignment base)
+ * @throws ExceptionUnmatchedClosingParentheses
+ * @throws ExceptionFileFormatOrSyntax
+ */
+ public void setRNA(String seq, String struct, ArrayList<Integer> basesOwn)
+ throws ExceptionUnmatchedClosingParentheses,
+ ExceptionFileFormatOrSyntax {
+ clearAnnotations();
+ _listeBases = new ArrayList<ModeleBase>();
+ // On "parse" la structure (repérage des points, tiret et couples
+ // parentheses ouvrante/fermante)
+ int[] array_struct = parseStruct(struct);
+ int size = struct.length();
+ int j = 0;
+ for (int i = 0; i < size; i++) {
+ ModeleBase mb;
+ if (seq.charAt(j) != seq.charAt(j + 1)) {
+ ModeleBasesComparison mbc = new ModeleBasesComparison(
+ seq.charAt(j), seq.charAt(j + 1), i);
+ mbc.set_appartenance(basesOwn.get(i));
+ mbc.setBaseNumber(i + 1);
+ mb = mbc;
+ } else {
+ mb = new ModeleBaseNucleotide("" + seq.charAt(j), i, i + 1);
+
+ }
+ _listeBases.add(mb);
+ j += 2;
+ }
+ for (int i = 0; i < size; i++) {
+ if (array_struct[i] != -1) {
+ this.addBPNow(i, array_struct[i]);
+ }
+
+ j += 2;
+ }
+ }
+
+ public void setRNA(List<String> seq, String dbnStr)
+ throws ExceptionUnmatchedClosingParentheses,
+ ExceptionFileFormatOrSyntax {
+ clearAnnotations();
+ int[] finStr = RNAFactory.parseSecStr(dbnStr);
+ setRNA(seq, finStr);
+ }
+
+ public static ArrayList<String> explodeSequence(String seq) {
+ ArrayList<String> analyzedSeq = new ArrayList<String>();
+ int i = 0;
+ while (i < seq.length()) {
+ if (seq.charAt(i) == '{') {
+ boolean found = false;
+ String buf = "";
+ i++;
+ while (!found & (i < seq.length())) {
+ if (seq.charAt(i) != '}') {
+ buf += seq.charAt(i);
+ i++;
+ } else {
+ found = true;
+ }
+ }
+ analyzedSeq.add(buf);
+ } else {
+ analyzedSeq.add("" + seq.charAt(i));
+ }
+ i++;
+ }
+ return analyzedSeq;
+ }
+
+ public int[] parseStruct(String str)
+ throws ExceptionUnmatchedClosingParentheses,
+ ExceptionFileFormatOrSyntax {
+ int[] result = new int[str.length()];
+ int unexpectedChar = -1;
+ Stack<Integer> p = new Stack<Integer>();
+ for (int i = 0; i < str.length(); i++) {
+ char c = str.charAt(i);
+ if (c == '(') {
+ p.push(new Integer(i));
+ } else if (c == '.' || c == '-' || c == ':') {
+ result[i] = -1;
+ } else if (c == ')') {
+ if (p.size() == 0) {
+ throw new ExceptionUnmatchedClosingParentheses(i + 1);
+ }
+ int j = p.pop().intValue();
+ result[i] = j;
+ result[j] = i;
+ } else {
+ if (unexpectedChar == -1)
+ unexpectedChar = i;
+ break;
+ }
+ }
+
+ if (unexpectedChar != -1) {
+ // warningEmition("Unexpected Character at index:" +
+ // unexpectedChar);
+ }
+
+ if (p.size() != 0) {
+ throw new ExceptionUnmatchedClosingParentheses(
+ p.pop().intValue() + 1);
+ }
+
+ return result;
+ }
+
+ public Point getHelixInterval(int index) {
+ if ((index < 0) || (index >= _listeBases.size())) {
+ return new Point(index, index);
+ }
+ int j = _listeBases.get(index).getElementStructure();
+ if (j != -1) {
+ int minH = index;
+ int maxH = index;
+ if (j > index) {
+ maxH = j;
+ } else {
+ minH = j;
+ }
+ boolean over = false;
+ while (!over) {
+ if ((minH < 0) || (maxH >= _listeBases.size())) {
+ over = true;
+ } else {
+ if (_listeBases.get(minH).getElementStructure() == maxH) {
+ minH--;
+ maxH++;
+ } else {
+ over = true;
+ }
+ }
+ }
+ minH++;
+ maxH--;
+ return new Point(minH, maxH);
+ }
+ return new Point(0, 0);
+ }
+
+ public Point getExteriorHelix(int index) {
+ Point h = getHelixInterval(index);
+ int a = h.x;
+ int b = h.y;
+ while (!((h.x==0)))
+ {
+ a = h.x;
+ b = h.y;
+ h = getHelixInterval(a-1);
+ }
+ return new Point(a, b);
+ }
+
+
+ public ArrayList<Integer> getHelix(int index) {
+ ArrayList<Integer> result = new ArrayList<Integer>();
+ if ((index < 0) || (index >= _listeBases.size())) {
+ return result;
+ }
+ Point p = getHelixInterval(index);
+ for (int i = p.x; i <= p.y; i++) {
+ result.add(i);
+ result.add(this._listeBases.get(i).getElementStructure());
+ }
+ return result;
+ }
+
+ public Point getMultiLoop(int index) {
+ if ((index < 0) || (index >= _listeBases.size())) {
+ return new Point(index, index);
+ }
+ Point h = getHelixInterval(index);
+ int minH = h.x - 1;
+ int maxH = h.y + 1;
+ boolean over = false;
+ while (!over) {
+ if (minH < 0) {
+ over = true;
+ minH = 0;
+ } else {
+ if (_listeBases.get(minH).getElementStructure() == -1) {
+ minH--;
+ } else if (_listeBases.get(minH).getElementStructure() < minH) {
+ minH = _listeBases.get(minH).getElementStructure() - 1;
+ } else {
+ over = true;
+ }
+ }
+ }
+ over = false;
+ while (!over) {
+ if (maxH > _listeBases.size() - 1) {
+ over = true;
+ maxH = _listeBases.size() - 1;
+ } else {
+ if (_listeBases.get(maxH).getElementStructure() == -1) {
+ maxH++;
+ } else if (_listeBases.get(maxH).getElementStructure() > maxH) {
+ maxH = _listeBases.get(maxH).getElementStructure() + 1;
+ } else {
+ over = true;
+ }
+ }
+ }
+ return new Point(minH, maxH);
+ }
+
+ public Vector<Integer> getLoopBases(int startIndex) {
+ Vector<Integer> result = new Vector<Integer>();
+
+ if ((startIndex < 0) || (startIndex >= _listeBases.size())) {
+ return result;
+ }
+ int index = startIndex;
+ result.add(startIndex);
+ if (_listeBases.get(index).getElementStructure() <= index) {
+ index = (index + 1) % _listeBases.size();
+ } else {
+ index = _listeBases.get(index).getElementStructure();
+ result.add(index);
+ index = (index + 1) % _listeBases.size();
+ }
+
+ while (index != startIndex) {
+ result.add(index);
+ if (_listeBases.get(index).getElementStructure() == -1) {
+ index = (index + 1) % _listeBases.size();
+ } else {
+ index = _listeBases.get(index).getElementStructure();
+ result.add(index);
+ index = (index + 1) % _listeBases.size();
+ }
+ }
+ return result;
+ }
+
+ /**
+ * Returns the RNA secondary structure displayed by this panel as a
+ * well-parenthesized word, accordingly to the DBN format
+ *
+ * @return This panel's secondary structure
+ */
+ public String getStructDBN() {
+ String result = "";
+ for (int i = 0; i < _listeBases.size(); i++) {
+ int j = _listeBases.get(i).getElementStructure();
+ if (j == -1) {
+ result += ".";
+ } else if (i > j) {
+ result += ")";
+ } else {
+ result += "(";
+ }
+ }
+ return addStrandSeparators(result);
+ }
+
+ private ArrayList<ModeleBP> getNonCrossingSubset(
+ ArrayList<ArrayList<ModeleBP>> rankedBPs) {
+ ArrayList<ModeleBP> currentBPs = new ArrayList<ModeleBP>();
+ Stack<Integer> pile = new Stack<Integer>();
+ for (int i = 0; i < rankedBPs.size(); i++) {
+ ArrayList<ModeleBP> lbp = rankedBPs.get(i);
+ if (!lbp.isEmpty()) {
+ ModeleBP bp = lbp.get(0);
+ boolean ok = true;
+ if (!pile.empty()) {
+ int x = pile.peek();
+ if ((bp.getIndex3() >= x)) {
+ ok = false;
+ }
+ }
+ if (ok) {
+ lbp.remove(0);
+ currentBPs.add(bp);
+ pile.add(bp.getIndex3());
+ }
+ }
+ if (!pile.empty() && (i == pile.peek())) {
+ pile.pop();
+ }
+ }
+ return currentBPs;
+ }
+
+ public ArrayList<int[]> paginateStructure() {
+ ArrayList<int[]> result = new ArrayList<int[]>();
+ // Mumbo jumbo to sort the basepair list
+ ArrayList<ModeleBP> bps = this.getAllBPs();
+ ModeleBP[] mt = new ModeleBP[bps.size()];
+ bps.toArray(mt);
+ Arrays.sort(mt, new Comparator<ModeleBP>() {
+ public int compare(ModeleBP arg0, ModeleBP arg1) {
+ if (arg0.getIndex5() != arg1.getIndex5())
+ return arg0.getIndex5() - arg1.getIndex5();
+ else
+ return arg0.getIndex3() - arg1.getIndex3();
+
+ }
+ });
+ ArrayList<ArrayList<ModeleBP>> rankedBps = new ArrayList<ArrayList<ModeleBP>>();
+ for (int i = 0; i < getSize(); i++) {
+ rankedBps.add(new ArrayList<ModeleBP>());
+ }
+ for (int i = 0; i < mt.length; i++) {
+ rankedBps.get(mt[i].getIndex5()).add(mt[i]);
+ }
+
+ while (!bps.isEmpty()) {
+ //System.out.println("Page: " + result.size());
+ ArrayList<ModeleBP> currentBPs = getNonCrossingSubset(rankedBps);
+ int[] ss = new int[this.getSize()];
+ for (int i = 0; i < ss.length; i++) {
+ ss[i] = -1;
+ }
+
+ for (int i = 0; i < currentBPs.size(); i++) {
+ ModeleBP mbp = currentBPs.get(i);
+ ss[mbp.getIndex3()] = mbp.getIndex5();
+ ss[mbp.getIndex5()] = mbp.getIndex3();
+ }
+ bps.removeAll(currentBPs);
+ result.add(ss);
+ }
+ return result;
+ }
+
+ private void showBasic(int[] res) {
+ for (int i = 0; i < res.length; i++) {
+ System.out.print(res[i] + ",");
+ }
+ System.out.println();
+
+ }
+
+ public int[] getStrandShifts()
+ {
+ int[] result = new int[getSize()];
+ int acc = 0;
+ for (int i=0;i<getSize();i++)
+ {
+ if (_backbone.getTypeBefore(i)==BackboneType.DISCONTINUOUS_TYPE)
+ {
+ acc++;
+ }
+ result[i] = acc;
+ }
+ return result;
+
+ }
+
+ public String addStrandSeparators(String s)
+ {
+ String res = "";
+ for (int i=0;i<s.length();i++)
+ {
+ res += s.charAt(i);
+ if (_backbone.getTypeAfter(i)==BackboneType.DISCONTINUOUS_TYPE)
+ {
+ res += DBNStrandSep;
+ }
+ }
+ return res;
+ }
+
+
+ public String getStructDBN(boolean includeMostPKs) {
+ String result = getStructDBN();
+ if (includeMostPKs) {
+ ArrayList<int[]> pages = paginateStructure();
+ char[] res = new char[getSize()];
+ for (int i = 0; i < res.length; i++) {
+ res[i] = '.';
+ }
+ char[] open = { '(', '[', '{', '<', 'A', 'B', 'C', 'D', 'E', 'F',
+ 'G', 'H', 'I', 'J', 'K', 'L', 'M', 'N',
+ 'O', 'P', 'Q', 'R', 'S', 'T', 'U', 'V', 'W', 'X', 'Y', 'Z'};
+
+ char[] close = new char[open.length];
+ close[0] = ')';
+ close[1] = ']';
+ close[2] = '}';
+ close[3] = '>';
+ for (int i=4; i<open.length;i++)
+ {
+ close[i] = Character.toLowerCase(open[i]);
+ }
+
+ for (int p = 0; p < Math.min(pages.size(), open.length); p++) {
+ int[] page = pages.get(p);
+ //showBasic(page);
+ for (int i = 0; i < res.length; i++) {
+ if (page[i] != -1 && page[i] > i && res[i] == '.'
+ && res[page[i]] == '.') {
+ res[i] = open[p];
+ res[page[i]] = close[p];
+ }
+ }
+ }
+ result = "";
+ for (int i = 0; i < res.length; i++) {
+ result += res[i];
+ }
+
+ }
+ return addStrandSeparators(result);
+
+ }
+
+ public String getStructDBN(int[] str) {
+ String result = "";
+ for (int i = 0; i < str.length; i++) {
+ if (str[i] == -1) {
+ result += ".";
+ } else if (str[i] > i) {
+ result += "(";
+ } else {
+ result += ")";
+ }
+ }
+ return addStrandSeparators(result);
+ }
+
+ /**
+ * Returns the raw nucleotides sequence for the displayed RNA
+ *
+ * @return The RNA sequence
+ */
+ public String getSeq() {
+ String result = "";
+ for (int i = 0; i < _listeBases.size(); i++) {
+ result += ((ModeleBase) _listeBases.get(i)).getContent();
+ }
+ return addStrandSeparators(result);
+ }
+
+ public String getStructBPSEQ() {
+ String result = "";
+ int[] str = getNonOverlappingStruct();
+ for (int i = 0; i < _listeBases.size(); i++) {
+ result += (i + 1) + " "
+ + ((ModeleBaseNucleotide) _listeBases.get(i)).getContent()
+ + " " + (str[i] + 1) + "\n";
+ }
+ return result;
+ }
+
+ public int[] getNonCrossingStruct() {
+ int[] result = new int[_listeBases.size()];
+ // Adding "planar" base-pairs
+ for (int i = 0; i < _listeBases.size(); i++) {
+ result[i] = _listeBases.get(i).getElementStructure();
+ }
+ return result;
+ }
+
+ public int[] getNonOverlappingStruct() {
+ int[] result = getNonCrossingStruct();
+ // Adding additional base pairs when possible (No more than one
+ // base-pair per base)
+ for (int i = 0; i < _structureAux.size(); i++) {
+ ModeleBP msbp = _structureAux.get(i);
+ ModeleBase mb5 = msbp.getPartner5();
+ ModeleBase mb3 = msbp.getPartner3();
+ int j5 = mb5.getIndex();
+ int j3 = mb3.getIndex();
+ if ((result[j3] == -1) && (result[j5] == -1)) {
+ result[j3] = j5;
+ result[j5] = j3;
+ }
+ }
+ return result;
+ }
+
+ public String getStructCT() {
+ String result = "";
+ for (int i = 0; i < _listeBases.size(); i++) {
+ result += (i + 1) + " "
+ + ((ModeleBase) _listeBases.get(i)).getContent() + " " + i
+ + " " + (i + 2) + " "
+ + (_listeBases.get(i).getElementStructure() + 1) + " "
+ + (i + 1) + "\n";
+ }
+ return result;
+ }
+
+ public void saveAsBPSEQ(String path, String title)
+ throws ExceptionExportFailed, ExceptionPermissionDenied {
+ try {
+ FileWriter f = new FileWriter(path);
+ f.write("# " + title + "\n");
+ f.write(this.getStructBPSEQ() + "\n");
+ f.close();
+ } catch (IOException e) {
+ throw new ExceptionExportFailed(e.getMessage(), path);
+ }
+ }
+
+ public void saveAsCT(String path, String title)
+ throws ExceptionExportFailed, ExceptionPermissionDenied {
+ try {
+ FileWriter f = new FileWriter(path);
+ f.write("" + _listeBases.size() + " " + title + "\n");
+ f.write(this.getStructCT() + "\n");
+ f.close();
+ } catch (IOException e) {
+ throw new ExceptionExportFailed(e.getMessage(), path);
+ }
+ }
+
+ public void saveAsDBN(String path, String title)
+ throws ExceptionExportFailed, ExceptionPermissionDenied {
+ try {
+ FileWriter f = new FileWriter(path);
+ f.write("> " + title + "\n");
+ f.write(getListeBasesToString() + "\n");
+ f.write(getStructDBN() + "\n");
+ f.close();
+ } catch (IOException e) {
+ throw new ExceptionExportFailed(e.getMessage(), path);
+ }
+ }
+
+ public String getListeBasesToString() {
+ String s = new String();
+ for (int i = 0; i < _listeBases.size(); i++) {
+ s += ((ModeleBaseNucleotide) _listeBases.get(i)).getContent();
+ }
+ return addStrandSeparators(s);
+ }
+
+ public void applyBPs(ArrayList<ModeleBP> allbps) {
+ ArrayList<ModeleBP> planar = new ArrayList<ModeleBP>();
+ ArrayList<ModeleBP> others = new ArrayList<ModeleBP>();
+ // System.err.println("Sequence: "+this.getSeq());
+ RNAMLParser.planarize(allbps, planar, others, getSize());
+ // System.err.println("All:"+allbps);
+ // System.err.println("=> Planar: "+planar);
+ // System.err.println("=> Others: "+others);
+
+ for (ModeleBP mb : planar) {
+ addBPnow(mb.getPartner5().getIndex(), mb.getPartner3().getIndex(),
+ mb);
+ }
+
+ for (ModeleBP mb : others) {
+ addBPAux(mb.getPartner5().getIndex(), mb.getPartner3().getIndex(),
+ mb);
+ }
+ }
+
+ public void set_listeBases(ArrayList<ModeleBase> _liste) {
+ this._listeBases = _liste;
+ }
+
+ public void addVARNAListener(InterfaceVARNAListener rl) {
+ _listeVARNAListener.add(rl);
+ }
+
+ public void warningEmition(String warningMessage) {
+ for (int i = 0; i < _listeVARNAListener.size(); i++) {
+ _listeVARNAListener.get(i).onWarningEmitted(warningMessage);
+ }
+ }
+
+ public void applyStyleOnBases(ArrayList<Integer> basesList,
+ ModelBaseStyle style) {
+ for (int i = 1; i < basesList.size(); i++) {
+ _listeBases.get(basesList.get(i)).setStyleBase(style);
+ }
+ }
+
+ private int[] correctReciprocity(int[] str) {
+ int[] result = new int[str.length];
+ for (int i = 0; i < str.length; i++) {
+ if (str[i] != -1) {
+ if (i == str[str[i]]) {
+ result[i] = str[i];
+ } else {
+ str[str[i]] = i;
+ }
+ } else {
+ result[i] = -1;
+ }
+ }
+ return result;
+ }
+
+ private void applyStruct(int[] str) throws ExceptionFileFormatOrSyntax {
+ str = correctReciprocity(str);
+
+ int[] planarSubset = RNAMLParser.planarize(str);
+ _structureAux.clear();
+
+ for (int i = 0; i < planarSubset.length; i++) {
+ if (str[i] > i) {
+ if (planarSubset[i] > i) {
+ addBPNow(i, planarSubset[i]);
+ } else if ((planarSubset[i] != str[i])) {
+ addBPAux(i, str[i]);
+ }
+ }
+ }
+
+ }
+
+ public ArrayList<ModeleBase> get_listeBases() {
+ return _listeBases;
+ }
+
+ public int getSize() {
+ return _listeBases.size();
+ }
+
+ public ArrayList<Integer> findAll() {
+ ArrayList<Integer> listAll = new ArrayList<Integer>();
+ for (int i = 0; i < get_listeBases().size(); i++) {
+ listAll.add(i);
+ }
+ return listAll;
+ }
+
+ public ArrayList<Integer> findBulge(int index) {
+ ArrayList<Integer> listUp = new ArrayList<Integer>();
+ if (get_listeBases().get(index).getElementStructure() == -1) {
+ int i = index;
+ boolean over = false;
+ while ((i < get_listeBases().size()) && !over) {
+ int j = get_listeBases().get(i).getElementStructure();
+ if (j == -1) {
+ listUp.add(i);
+ i++;
+ } else {
+ over = true;
+ }
+ }
+ i = index - 1;
+ over = false;
+ while ((i >= 0) && !over) {
+ int j = get_listeBases().get(i).getElementStructure();
+ if (j == -1) {
+ listUp.add(i);
+ i--;
+ } else {
+ over = true;
+ }
+ }
+ }
+ return listUp;
+ }
+
+ public ArrayList<Integer> findStem(int index) {
+ ArrayList<Integer> listUp = new ArrayList<Integer>();
+ int i = index;
+ do {
+ listUp.add(i);
+ int j = get_listeBases().get(i).getElementStructure();
+ if (j == -1) {
+ i = (i + 1) % getSize();
+ } else {
+ if ((j < i) && (index <= i) && (j <= index)) {
+ i = j;
+ } else {
+ i = (i + 1) % getSize();
+ }
+ }
+ } while (i != index);
+ return listUp;
+ }
+
+ public int getHelixCountOnLoop(int indice) {
+ int cptHelice = 0;
+ if (indice < 0 || indice >= get_listeBases().size())
+ return cptHelice;
+ int i = indice;
+ int j = get_listeBases().get(i).getElementStructure();
+ // Only way to distinguish "supporting base-pair" from others
+ boolean justJumped = false;
+ if ((j != -1) && (j < i)) {
+ i = j + 1;
+ indice = i;
+ }
+ do {
+ j = get_listeBases().get(i).getElementStructure();
+ if ((j != -1) && (!justJumped)) {
+ i = j;
+ justJumped = true;
+ cptHelice++;
+ } else {
+ i = (i + 1) % get_listeBases().size();
+ justJumped = false;
+ }
+ } while (i != indice);
+ return cptHelice;
+ }
+
+ public ArrayList<Integer> findLoop(int indice) {
+ return findLoopForward(indice);
+ }
+
+ public ArrayList<Integer> findLoopForward(int indice) {
+ ArrayList<Integer> base = new ArrayList<Integer>();
+ if (indice < 0 || indice >= get_listeBases().size())
+ return base;
+ int i = indice;
+ int j = get_listeBases().get(i).getElementStructure();
+ // Only way to distinguish "supporting base-pair" from others
+ boolean justJumped = false;
+ if (j != -1) {
+ i = Math.min(i, j) + 1;
+ indice = i;
+ }
+ do {
+ base.add(i);
+ j = get_listeBases().get(i).getElementStructure();
+ if ((j != -1) && (!justJumped)) {
+ i = j;
+ justJumped = true;
+ } else {
+ i = (i + 1) % get_listeBases().size();
+ justJumped = false;
+ }
+ } while (i != indice);
+ return base;
+ }
+
+ public ArrayList<Integer> findPair(int indice) {
+ ArrayList<Integer> base = new ArrayList<Integer>();
+ int j = get_listeBases().get(indice).getElementStructure();
+ if (j != -1) {
+ base.add(Math.min(indice, j));
+ base.add(Math.max(indice, j));
+ }
+
+ return base;
+
+ }
+
+ public ArrayList<Integer> findLoopBackward(int indice) {
+ ArrayList<Integer> base = new ArrayList<Integer>();
+ if (indice < 0 || indice >= get_listeBases().size())
+ return base;
+ int i = indice;
+ int j = get_listeBases().get(i).getElementStructure();
+ // Only way to distinguish "supporting base-pair" from others
+ boolean justJumped = false;
+ if (j != -1) {
+ i = Math.min(i, j) - 1;
+ indice = i;
+ }
+ if (i < 0) {
+ return base;
+ }
+ do {
+ base.add(i);
+ j = get_listeBases().get(i).getElementStructure();
+ if ((j != -1) && (!justJumped)) {
+ i = j;
+ justJumped = true;
+ } else {
+ i = (i + get_listeBases().size() - 1) % get_listeBases().size();
+ justJumped = false;
+ }
+ } while (i != indice);
+ return base;
+ }
+
+ public ArrayList<Integer> findHelix(int indice) {
+ ArrayList<Integer> list = new ArrayList<Integer>();
+ if (get_listeBases().get(indice).getElementStructure() != -1) {
+ list.add(indice);
+ list.add(get_listeBases().get(indice).getElementStructure());
+ int i = 1, prec = get_listeBases().get(indice)
+ .getElementStructure();
+ while (indice + i < get_listeBases().size()
+ && get_listeBases().get(indice + i).getElementStructure() != -1
+ && get_listeBases().get(indice + i).getElementStructure() == prec - 1) {
+ list.add(indice + i);
+ list.add(get_listeBases().get(indice + i).getElementStructure());
+ prec = get_listeBases().get(indice + i).getElementStructure();
+ i++;
+ }
+ i = -1;
+ prec = get_listeBases().get(indice).getElementStructure();
+ while (indice + i >= 0
+ && get_listeBases().get(indice + i).getElementStructure() != -1
+ && get_listeBases().get(indice + i).getElementStructure() == prec + 1) {
+ list.add(indice + i);
+ list.add(get_listeBases().get(indice + i).getElementStructure());
+ prec = get_listeBases().get(indice + i).getElementStructure();
+ i--;
+ }
+ }
+ return list;
+ }
+
+ public ArrayList<Integer> find3Prime(int indice) {
+ ArrayList<Integer> list = new ArrayList<Integer>();
+ boolean over = false;
+ while ((indice >= 0) && !over) {
+ over = (get_listeBases().get(indice).getElementStructure() != -1);
+ indice--;
+ }
+ indice++;
+ if (over) {
+ indice++;
+ }
+ for (int i = indice; i < get_listeBases().size(); i++) {
+ list.add(i);
+ if (get_listeBases().get(i).getElementStructure() != -1) {
+ return new ArrayList<Integer>();
+ }
+ }
+ return list;
+ }
+
+ public ArrayList<Integer> find5Prime(int indice) {
+ ArrayList<Integer> list = new ArrayList<Integer>();
+ for (int i = 0; i <= indice; i++) {
+ list.add(i);
+ if (get_listeBases().get(i).getElementStructure() != -1) {
+ return new ArrayList<Integer>();
+ }
+ }
+ return list;
+ }
+
+ public static Double angle(Point2D.Double p1, Point2D.Double p2,
+ Point2D.Double p3) {
+ Double alpha = Math.atan2(p1.y - p2.y, p1.x - p2.x);
+ Double beta = Math.atan2(p3.y - p2.y, p3.x - p2.x);
+ Double angle = (beta - alpha);
+
+ // Correction de l'angle pour le resituer entre 0 et 2PI
+ while (angle < 0.0 || angle > 2 * Math.PI) {
+ if (angle < 0.0)
+ angle += 2 * Math.PI;
+ else if (angle > 2 * Math.PI)
+ angle -= 2 * Math.PI;
+ }
+ return angle;
+ }
+
+ public ArrayList<Integer> findNonPairedBaseGroup(Integer get_nearestBase) {
+ // detection 3', 5', bulge
+ ArrayList<Integer> list = new ArrayList<Integer>();
+ int indice = get_nearestBase;
+ boolean nonpairedUp = true, nonpairedDown = true;
+ while (indice < get_listeBases().size() && nonpairedUp) {
+ if (get_listeBases().get(indice).getElementStructure() == -1) {
+ list.add(indice);
+ indice++;
+ } else {
+ nonpairedUp = false;
+ }
+ }
+ indice = get_nearestBase - 1;
+ while (indice >= 0 && nonpairedDown) {
+ if (get_listeBases().get(indice).getElementStructure() == -1) {
+ list.add(indice);
+ indice--;
+ } else {
+ nonpairedDown = false;
+ }
+ }
+ return list;
+ }
+
+ /*
+ * public boolean getDrawn() { return _drawn; }
+ */
+
+ public ArrayList<ModeleBP> getStructureAux() {
+ return _structureAux;
+ }
+
+ /**
+ * Translates a base number into its corresponding index. Although both
+ * should be unique, base numbers are not necessarily contiguous, and
+ * indices should be preferred for any reasonably complex algorithmic
+ * treatment.
+ *
+ * @param num
+ * The base number
+ * @return The first index whose associated Base model has base number
+ * <code>num</code>, <code>-1</code> of no such base model exists.
+ */
+
+ public int getIndexFromBaseNumber(int num) {
+ for (int i = 0; i < this._listeBases.size(); i++) {
+ if (_listeBases.get(i).getBaseNumber() == num) {
+ return i;
+ }
+ }
+ return -1;
+ }
+
+ /**
+ * Adds a base pair to this RNA's structure. Tries to add it to the
+ * secondary structure first, eventually adding it to the 'tertiary'
+ * interactions if it clashes with the current secondary structure.
+ *
+ * @param baseNumber5
+ * - Base number of the origin of this base pair
+ * @param baseNumber3
+ * - Base number of the destination of this base pair
+ */
+
+ public void addBPToStructureUsingNumbers(int baseNumber5, int baseNumber3) {
+ int i = getIndexFromBaseNumber(baseNumber5);
+ int j = getIndexFromBaseNumber(baseNumber3);
+ addBP(i, j);
+ }
+
+ /**
+ * Adds a base pair to this RNA's structure. Tries to add it to the
+ * secondary structure first, possibly adding it to the 'tertiary'
+ * interactions if it clashes with the current secondary structure.
+ *
+ * @param number5
+ * - Base number of the origin of this base pair
+ * @param number3
+ * - Base number of the destination of this base pair
+ */
+
+ public void addBPToStructureUsingNumbers(int number5, int number3,
+ ModeleBP msbp) {
+ addBP(getIndexFromBaseNumber(number5), getIndexFromBaseNumber(number3),
+ msbp);
+ }
+
+ public void addBP(int index5, int index3) {
+ int i = index5;
+ int j = index3;
+ ModeleBase part5 = _listeBases.get(i);
+ ModeleBase part3 = _listeBases.get(j);
+ ModeleBP msbp = new ModeleBP(part5, part3);
+ addBP(i, j, msbp);
+ }
+
+ public void addBP(int index5, int index3, ModeleBP msbp) {
+ int i = index5;
+ int j = index3;
+
+ if (j < i) {
+ int k = j;
+ j = i;
+ i = k;
+ }
+ if (i != -1) {
+ for (int k = i; k <= j; k++) {
+ ModeleBase tmp = _listeBases.get(k);
+ int l = tmp.getElementStructure();
+ if (l != -1) {
+ if ((l <= i) || (l >= j)) {
+ addBPAux(i, j, msbp);
+ return;
+ }
+ }
+ }
+ addBPnow(i, j, msbp);
+ }
+ }
+
+ public void removeBP(ModeleBP ms) {
+ if (_structureAux.contains(ms)) {
+ _structureAux.remove(ms);
+ } else {
+ ModeleBase m5 = ms.getPartner5();
+ ModeleBase m3 = ms.getPartner3();
+ int i = m5.getIndex();
+ int j = m3.getIndex();
+ if ((m5.getElementStructure() == m3.getIndex())
+ && (m3.getElementStructure() == m5.getIndex())) {
+ m5.removeElementStructure();
+ m3.removeElementStructure();
+ }
+ }
+ }
+
+ /**
+ * Register base-pair, no question asked. More precisely, this function will
+ * not try to determine if the base-pairs crosses any other.
+ *
+ * @param i
+ * @param j
+ * @param msbp
+ */
+ private void addBPNow(int i, int j) {
+ if (j < i) {
+ int k = j;
+ j = i;
+ i = k;
+ }
+
+ ModeleBase part5 = _listeBases.get(i);
+ ModeleBase part3 = _listeBases.get(j);
+ ModeleBP msbp = new ModeleBP(part5, part3);
+ addBPnow(i, j, msbp);
+ }
+
+ /**
+ * Register base-pair, no question asked. More precisely, this function will
+ * not try to determine if the base-pairs crosses any other.
+ *
+ * @param i
+ * @param j
+ * @param msbp
+ */
+ private void addBPnow(int i, int j, ModeleBP msbp) {
+ if (j < i) {
+ int k = j;
+ j = i;
+ i = k;
+ }
+ ModeleBase part5 = _listeBases.get(i);
+ ModeleBase part3 = _listeBases.get(j);
+ msbp.setPartner5(part5);
+ msbp.setPartner3(part3);
+ part5.setElementStructure(j, msbp);
+ part3.setElementStructure(i, msbp);
+ }
+
+ public void addBPAux(int i, int j) {
+ ModeleBase part5 = _listeBases.get(i);
+ ModeleBase part3 = _listeBases.get(j);
+ ModeleBP msbp = new ModeleBP(part5, part3);
+ addBPAux(i, j, msbp);
+ }
+
+ public void addBPAux(int i, int j, ModeleBP msbp) {
+ if (j < i) {
+ int k = j;
+ j = i;
+ i = k;
+ }
+ ModeleBase part5 = _listeBases.get(i);
+ ModeleBase part3 = _listeBases.get(j);
+ msbp.setPartner5(part5);
+ msbp.setPartner3(part3);
+ _structureAux.add(msbp);
+ }
+
+ public ArrayList<ModeleBP> getBPsAt(int i) {
+ ArrayList<ModeleBP> result = new ArrayList<ModeleBP>();
+ if (_listeBases.get(i).getElementStructure() != -1) {
+ result.add(_listeBases.get(i).getStyleBP());
+ }
+ for (int k = 0; k < _structureAux.size(); k++) {
+ ModeleBP bp = _structureAux.get(k);
+ if ((bp.getPartner5().getIndex() == i)
+ || (bp.getPartner3().getIndex() == i)) {
+ result.add(bp);
+ }
+ }
+ return result;
+
+ }
+
+ public ModeleBP getBPStyle(int i, int j) {
+ ModeleBP result = null;
+ if (i > j) {
+ int k = j;
+ j = i;
+ i = k;
+ }
+ if (_listeBases.get(i).getElementStructure() == j) {
+ result = _listeBases.get(i).getStyleBP();
+ }
+ for (int k = 0; k < _structureAux.size(); k++) {
+ ModeleBP bp = _structureAux.get(k);
+ if ((bp.getPartner5().getIndex() == i)
+ && (bp.getPartner3().getIndex() == j)) {
+ result = bp;
+ }
+ }
+ return result;
+ }
+
+ public ArrayList<ModeleBP> getSecStrBPs() {
+ ArrayList<ModeleBP> result = new ArrayList<ModeleBP>();
+ for (int i = 0; i < this.getSize(); i++) {
+ ModeleBase mb = _listeBases.get(i);
+ int k = mb.getElementStructure();
+ if ((k != -1) && (k > i)) {
+ result.add(mb.getStyleBP());
+ }
+ }
+ return result;
+ }
+
+ public ArrayList<ModeleBP> getAuxBPs() {
+ ArrayList<ModeleBP> result = new ArrayList<ModeleBP>();
+ for (ModeleBP bp : _structureAux) {
+ result.add(bp);
+ }
+ return result;
+ }
+
+ public ArrayList<ModeleBP> getAllBPs() {
+ ArrayList<ModeleBP> result = new ArrayList<ModeleBP>();
+ result.addAll(getSecStrBPs());
+ result.addAll(getAuxBPs());
+ return result;
+ }
+
+ public ArrayList<ModeleBP> getAuxBPs(int i) {
+ ArrayList<ModeleBP> result = new ArrayList<ModeleBP>();
+ for (ModeleBP bp : _structureAux) {
+ if ((bp.getPartner5().getIndex() == i)
+ || (bp.getPartner3().getIndex() == i)) {
+ result.add(bp);
+ }
+ }
+ return result;
+ }
+
+ public void setBaseInnerColor(Color c) {
+ for (int i = 0; i < _listeBases.size(); i++) {
+ ModeleBase mb = _listeBases.get(i);
+ mb.getStyleBase().setBaseInnerColor(c);
+ }
+ }
+
+ public void setBaseNumbersColor(Color c) {
+ for (int i = 0; i < _listeBases.size(); i++) {
+ ModeleBase mb = _listeBases.get(i);
+ mb.getStyleBase().setBaseNumberColor(c);
+ }
+ }
+
+ public void setBaseNameColor(Color c) {
+ for (int i = 0; i < _listeBases.size(); i++) {
+ ModeleBase mb = _listeBases.get(i);
+ mb.getStyleBase().setBaseNameColor(c);
+ }
+ }
+
+ public void setBaseOutlineColor(Color c) {
+ for (int i = 0; i < _listeBases.size(); i++) {
+ ModeleBase mb = _listeBases.get(i);
+ mb.getStyleBase().setBaseOutlineColor(c);
+ }
+ }
+
+ public String getName() {
+ return _name;
+ }
+
+ public void setName(String n) {
+ _name = n;
+ }
+
+ public ArrayList<TextAnnotation> getAnnotations() {
+ return _listeAnnotations;
+ }
+
+ public boolean removeAnnotation(TextAnnotation t) {
+ return _listeAnnotations.remove(t);
+ }
+
+ public void addAnnotation(TextAnnotation t) {
+ _listeAnnotations.add(t);
+ }
+
+ public void removeAnnotation(String filter) {
+ ArrayList<TextAnnotation> condamne = new ArrayList<TextAnnotation>();
+ for (TextAnnotation t : _listeAnnotations) {
+ if (t.getTexte().contains(filter)) {
+ condamne.add(t);
+ }
+ }
+ for (TextAnnotation t : condamne) {
+ _listeAnnotations.remove(t);
+ }
+ }
+
+ public void clearAnnotations() {
+ _listeAnnotations.clear();
+ }
+
+ private boolean _strandEndsAnnotated = false;
+
+ public void autoAnnotateStrandEnds() {
+ if (!_strandEndsAnnotated) {
+ int tailleListBases = _listeBases.size();
+ boolean endAnnotate = false;
+ addAnnotation(new TextAnnotation("5'", _listeBases.get(0)));
+ for (int i = 0; i < _listeBases.size() - 1; i++) {
+ int realposA = _listeBases.get(i).getBaseNumber();
+ int realposB = _listeBases.get(i + 1).getBaseNumber();
+ if (realposB - realposA != 1) {
+ addAnnotation(new TextAnnotation("3'", _listeBases.get(i)));
+ addAnnotation(new TextAnnotation("5'",
+ _listeBases.get(i + 1)));
+ if (i + 1 == _listeBases.size() - 1) {
+ endAnnotate = true;
+ }
+ }
+ }
+ if (!endAnnotate) {
+ addAnnotation(new TextAnnotation("3'",
+ _listeBases.get(tailleListBases - 1)));
+ }
+ _strandEndsAnnotated = true;
+ } else {
+ removeAnnotation("3'");
+ removeAnnotation("5'");
+ _strandEndsAnnotated = false;
+ }
+ }
+
+ public void autoAnnotateHelices() {
+ Stack<Integer> p = new Stack<Integer>();
+ p.push(0);
+ int nbH = 1;
+ while (!p.empty()) {
+ int i = p.pop();
+ if (i < _listeBases.size()) {
+ ModeleBase mb = _listeBases.get(i);
+ int j = mb.getElementStructure();
+ if (j == -1) {
+ p.push(i + 1);
+ } else {
+ if (j > i) {
+ ModeleBase mbp = _listeBases.get(j);
+ p.push(j + 1);
+ ArrayList<ModeleBase> h = new ArrayList<ModeleBase>();
+ int k = 1;
+ while (mb.getElementStructure() == mbp.getIndex()) {
+ h.add(mb);
+ h.add(mbp);
+ mb = _listeBases.get(i + k);
+ mbp = _listeBases.get(j - k);
+
+ k++;
+ }
+ try {
+ addAnnotation(new TextAnnotation("H" + nbH++, h,
+ TextAnnotation.AnchorType.HELIX));
+ } catch (Exception e) {
+ e.printStackTrace();
+ }
+ p.push(i + k);
+ }
+ }
+ }
+ }
+ }
+
+ public void autoAnnotateTerminalLoops() {
+ Stack<Integer> p = new Stack<Integer>();
+ p.push(0);
+ int nbT = 1;
+ while (!p.empty()) {
+ int i = p.pop();
+ if (i < _listeBases.size()) {
+ ModeleBase mb = _listeBases.get(i);
+ int j = mb.getElementStructure();
+ if (j == -1) {
+ int k = 1;
+ ArrayList<ModeleBase> t = new ArrayList<ModeleBase>();
+ while ((i + k < getSize())
+ && (mb.getElementStructure() == -1)) {
+ t.add(mb);
+ mb = _listeBases.get(i + k);
+ k++;
+ }
+ if (mb.getElementStructure() != -1) {
+ if (mb.getElementStructure() == i - 1) {
+ try {
+ t.add(_listeBases.get(i - 1));
+ t.add(_listeBases.get(i + k - 1));
+ addAnnotation(new TextAnnotation("T" + nbT++,
+ t, TextAnnotation.AnchorType.LOOP));
+ } catch (Exception e) {
+ e.printStackTrace();
+ }
+ }
+ p.push(i + k - 1);
+ }
+
+ } else {
+ if (j > i) {
+ p.push(j + 1);
+ p.push(i + 1);
+ }
+ }
+ }
+ }
+ }
+
+ public void autoAnnotateInteriorLoops() {
+ Stack<Integer> p = new Stack<Integer>();
+ p.push(0);
+ int nbT = 1;
+ while (!p.empty()) {
+ int i = p.pop();
+ if (i < _listeBases.size()) {
+ ModeleBase mb = _listeBases.get(i);
+ int j = mb.getElementStructure();
+ if (j == -1) {
+ int k = i + 1;
+ ArrayList<ModeleBase> t = new ArrayList<ModeleBase>();
+ boolean terminal = true;
+ while ((k < getSize())
+ && ((mb.getElementStructure() >= i) || (mb
+ .getElementStructure() == -1))) {
+ t.add(mb);
+ mb = _listeBases.get(k);
+ if ((mb.getElementStructure() == -1)
+ || (mb.getElementStructure() < k))
+ k++;
+ else {
+ p.push(k);
+ terminal = false;
+ k = mb.getElementStructure();
+ }
+ }
+ if (mb.getElementStructure() != -1) {
+ if ((mb.getElementStructure() == i - 1) && !terminal) {
+ try {
+ t.add(_listeBases.get(i - 1));
+ t.add(_listeBases.get(k - 1));
+ addAnnotation(new TextAnnotation("I" + nbT++,
+ t, TextAnnotation.AnchorType.LOOP));
+ } catch (Exception e) {
+ e.printStackTrace();
+ }
+ p.push(k - 1);
+ }
+ }
+ } else {
+ if (j > i) {
+ p.push(i + 1);
+ }
+ }
+ }
+ }
+ }
+
+ @SuppressWarnings("unchecked")
+ public TextAnnotation getAnnotation(TextAnnotation.AnchorType type,
+ ModeleBase base) {
+ TextAnnotation result = null;
+ for (TextAnnotation t : _listeAnnotations) {
+ if (t.getType() == type) {
+ switch (type) {
+ case BASE:
+ if (base == (ModeleBase) t.getAncrage())
+ return t;
+ break;
+ case HELIX:
+ case LOOP: {
+ ArrayList<ModeleBase> mbl = (ArrayList<ModeleBase>) t
+ .getAncrage();
+ if (mbl.contains(base))
+ return t;
+ }
+ break;
+ }
+ }
+ }
+ return result;
+ }
+
+ public void addChemProbAnnotation(ChemProbAnnotation cpa) {
+ //System.err.println(cpa.isOut());
+ _chemProbAnnotations.add(cpa);
+ }
+
+ public ArrayList<ChemProbAnnotation> getChemProbAnnotations() {
+ return _chemProbAnnotations;
+ }
+
+ public void setColorMapValues(Double[] values, ModeleColorMap cm) {
+ setColorMapValues(values, cm, false);
+ }
+
+ public void adaptColorMapToValues(ModeleColorMap cm) {
+ double min = Double.MAX_VALUE;
+ double max = Double.MIN_VALUE;
+ for (int i = 0; i < Math.min(_listeBases.size(), _listeBases.size()); i++) {
+ ModeleBase mb = _listeBases.get(i);
+ max = Math.max(max, mb.getValue());
+ min = Math.min(min, mb.getValue());
+ }
+ cm.rescale(min, max);
+ }
+
+
+ private ArrayList<Double> loadDotPlot(StreamTokenizer st)
+ {
+ ArrayList<Double> result = new ArrayList<Double>();
+ try {
+ boolean inSeq = false;
+ String sequence = "";
+ ArrayList<Double> accumulator = new ArrayList<Double>();
+ int type = st.nextToken();
+ Hashtable<Couple<Integer,Integer>,Double> BP = new Hashtable<Couple<Integer,Integer>,Double>();
+ while (type != StreamTokenizer.TT_EOF) {
+ switch (type) {
+ case (StreamTokenizer.TT_NUMBER):
+ accumulator.add(st.nval);
+ break;
+ case (StreamTokenizer.TT_EOL):
+ break;
+ case (StreamTokenizer.TT_WORD):
+ if (st.sval.equals("/sequence"))
+ {
+ inSeq = true;
+ }
+ else if (st.sval.equals("ubox"))
+ {
+ int i = accumulator.get(accumulator.size()-3).intValue()-1;
+ int j = accumulator.get(accumulator.size()-2).intValue()-1;
+ double val = accumulator.get(accumulator.size()-1);
+ //System.err.println((char) type);
+ BP.put(new Couple<Integer, Integer>(Math.min(i, j), Math.max(i, j)),val*val);
+ accumulator.clear();
+ }
+ else if (inSeq)
+ {
+ sequence += st.sval;
+ }
+ break;
+ case ')':
+ inSeq = false;
+ break;
+ }
+ type = st.nextToken();
+ }
+ for (int i = 0; i < getSize(); i++) {
+ int j = getBaseAt(i).getElementStructure();
+ if (j != -1) {
+ Couple<Integer, Integer> coor = new Couple<Integer, Integer>(
+ Math.min(i, j), Math.max(i, j));
+ if (BP.containsKey(coor)) {
+ result.add(BP.get(coor));
+ } else {
+ result.add(0.);
+ }
+ } else {
+ double acc = 1.0;
+ for (int k = 0; k < getSize(); k++) {
+ Couple<Integer, Integer> coor = new Couple<Integer, Integer>(
+ Math.min(i, k), Math.max(i, k));
+ if (BP.containsKey(coor)) {
+ acc -= BP.get(coor);
+ }
+ }
+ result.add(acc);
+ }
+ }
+ } catch (IOException e) {
+ // TODO Auto-generated catch block
+ e.printStackTrace();
+ }
+ return result;
+ }
+
+ public void readValues(Reader r, ModeleColorMap cm) {
+ try {
+ StreamTokenizer st = new StreamTokenizer(r);
+ st.eolIsSignificant(true);
+ st.wordChars('/', '/');
+ st.parseNumbers();
+ ArrayList<Double> vals = new ArrayList<Double>();
+ ArrayList<Double> curVals = new ArrayList<Double>();
+ int type = st.nextToken();
+ boolean isDotPlot = false;
+ if (type=='%')
+ {
+ vals = loadDotPlot(st);
+ isDotPlot = true;
+ }
+ else
+ {
+ while (type != StreamTokenizer.TT_EOF) {
+ switch (type) {
+ case (StreamTokenizer.TT_NUMBER):
+ curVals.add(st.nval);
+ break;
+ case (StreamTokenizer.TT_EOL):
+ if (curVals.size() > 0) {
+ vals.add(curVals.get(curVals.size()-1));
+ curVals = new ArrayList<Double>();
+ }
+ break;
+ }
+ type = st.nextToken();
+ }
+ if (curVals.size() > 0)
+ vals.add(curVals.get(curVals.size()-1));
+ }
+
+ Double[] v = new Double[vals.size()];
+ for (int i = 0; i < Math.min(vals.size(), getSize()); i++) {
+ v[i] = vals.get(i);
+ }
+ setColorMapValues(v, cm, true);
+ if (isDotPlot)
+ {
+ cm.setMinValue(0.0);
+ cm.setMaxValue(1.0);
+ }
+ } catch (IOException e) {
+ e.printStackTrace();
+ }
+ }
+
+ public void setColorMapValues(Double[] values, ModeleColorMap cm,
+ boolean rescaleColorMap) {
+ if (values.length > 0) {
+ for (int i = 0; i < Math.min(values.length, _listeBases.size()); i++) {
+ ModeleBase mb = _listeBases.get(i);
+ mb.setValue(values[i]);
+ }
+ if (rescaleColorMap) {
+ adaptColorMapToValues(cm);
+ }
+ }
+ }
+
+ public Double[] getColorMapValues() {
+ Double[] values = new Double[_listeBases.size()];
+ for (int i = 0; i < _listeBases.size(); i++) {
+ values[i] = _listeBases.get(i).getValue();
+ }
+ return values;
+ }
+
+ public void rescaleColorMap(ModeleColorMap cm) {
+ Double max = Double.MIN_VALUE;
+ Double min = Double.MAX_VALUE;
+ for (int i = 0; i < _listeBases.size(); i++) {
+ Double value = _listeBases.get(i).getValue();
+ max = Math.max(max, value);
+ min = Math.min(min, value);
+ }
+ cm.rescale(min, max);
+ }
+
+ public void addBase(ModeleBase mb) {
+ _listeBases.add(mb);
+ }
+
+ public void setSequence(String s) {
+ setSequence(RNA.explodeSequence(s));
+ }
+
+ public void setSequence(List<String> s) {
+ int i = 0;
+ int j = 0;
+ while ((i < s.size()) && (j < _listeBases.size())) {
+ ModeleBase mb = _listeBases.get(j);
+ if (mb instanceof ModeleBaseNucleotide) {
+ ((ModeleBaseNucleotide) mb).setBase(s.get(i));
+ i++;
+ j++;
+ } else if (mb instanceof ModeleBasesComparison) {
+ ((ModeleBasesComparison) mb)
+ .setBase1(((s.get(i).length() > 0) ? s.get(i).charAt(0)
+ : ' '));
+ ((ModeleBasesComparison) mb)
+ .setBase2(((s.get(i + 1).length() > 0) ? s.get(i + 1)
+ .charAt(0) : ' '));
+ i += 2;
+ j++;
+ } else
+ j++;
+ }
+ for (i = _listeBases.size(); i < s.size(); i++) {
+ _listeBases.add(new ModeleBaseNucleotide(s.get(i), i));
+ }
+ }
+
+ public void eraseSequence() {
+ int j = 0;
+ while ((j < _listeBases.size())) {
+ ModeleBase mb = _listeBases.get(j);
+ if (mb instanceof ModeleBaseNucleotide) {
+ ((ModeleBaseNucleotide) mb).setBase("");
+ j++;
+ } else if (mb instanceof ModeleBasesComparison) {
+ ((ModeleBasesComparison) mb).setBase1(' ');
+ ((ModeleBasesComparison) mb).setBase2(' ');
+ j++;
+ } else
+ j++;
+ }
+ }
+
+ public RNA clone() {
+ try {
+ ByteArrayOutputStream out = new ByteArrayOutputStream();
+ ObjectOutputStream oout = new ObjectOutputStream(out);
+ oout.writeObject(this);
+
+ ObjectInputStream in = new ObjectInputStream(
+ new ByteArrayInputStream(out.toByteArray()));
+ return (RNA) in.readObject();
+ } catch (Exception e) {
+ throw new RuntimeException("cannot clone class ["
+ + this.getClass().getName() + "] via serialization: "
+ + e.toString());
+ }
+ }
+
+ /**
+ * Returns the base at index <code>index</code>. Indices are contiguous in
+ * the sequence over an interval <code>[0,this.getSize()-1]</code>, where
+ * <code>n</code> is the length of the sequence.
+ *
+ * @param index
+ * The index, <code>0 ≤ index < this.getSize()</code>, of the
+ * base model
+ * @return The base model of index <code>index</code>
+ */
+ public ModeleBase getBaseAt(int index) {
+ return this._listeBases.get(index);
+ }
+
+ /**
+ * Returns the set of bases of indices in <code>indices</code>. Indices are
+ * contiguous in the sequence, and belong to an interval
+ * <code>[0,n-1]</code>, where <code>n</code> is the length of the sequence.
+ *
+ * @param indices
+ * A Collection of indices <code>i</code>,
+ * <code>0 ≤ index < this.getSize()</code>, where some base
+ * models are found.
+ * @return A list of base model of indices in <code>indices</code>
+ */
+ public ArrayList<ModeleBase> getBasesAt(
+ Collection<? extends Integer> indices) {
+ ArrayList<ModeleBase> mbs = new ArrayList<ModeleBase>();
+ for (int i : indices) {
+ mbs.add(getBaseAt(i));
+ }
+ return mbs;
+ }
+
+ public ArrayList<ModeleBase> getBasesBetween(int from, int to) {
+ ArrayList<ModeleBase> mbs = new ArrayList<ModeleBase>();
+ int bck = Math.min(from, to);
+ to = Math.max(from, to);
+ from = bck;
+ for (int i = from; i <= to; i++) {
+ mbs.add(getBaseAt(i));
+ }
+ return mbs;
+ }
+
+ public void addHighlightRegion(HighlightRegionAnnotation n) {
+ _listeRegionHighlights.add(n);
+ }
+
+ public void removeHighlightRegion(HighlightRegionAnnotation n) {
+ _listeRegionHighlights.remove(n);
+ }
+
+ public void removeChemProbAnnotation(ChemProbAnnotation a) {
+ _chemProbAnnotations.remove(a);
+ }
+
+ public void clearChemProbAnnotations() {
+ _chemProbAnnotations.clear();
+ }
+
+ public void addHighlightRegion(int from, int to, Color fill, Color outline,
+ double radius) {
+ _listeRegionHighlights.add(new HighlightRegionAnnotation(
+ getBasesBetween(from, to), fill, outline, radius));
+ }
+
+ public void addHighlightRegion(int from, int to) {
+ _listeRegionHighlights.add(new HighlightRegionAnnotation(
+ getBasesBetween(from, to)));
+ }
+
+ public ArrayList<HighlightRegionAnnotation> getHighlightRegion() {
+ return _listeRegionHighlights;
+ }
+
+ /**
+ * Rotates the RNA coordinates by a certain angle
+ *
+ * @param angleDegres
+ * Rotation angle, in degrees
+ */
+ public void globalRotation(Double angleDegres) {
+ if (_listeBases.size() > 0) {
+
+ // angle en radian
+ Double angle = angleDegres * Math.PI / 180;
+
+ // initialisation du minimum et dumaximum
+ Double maxX = _listeBases.get(0).getCoords().x;
+ Double maxY = _listeBases.get(0).getCoords().y;
+ Double minX = _listeBases.get(0).getCoords().x;
+ Double minY = _listeBases.get(0).getCoords().y;
+ // mise a jour du minimum et du maximum
+ for (int i = 0; i < _listeBases.size(); i++) {
+ if (_listeBases.get(i).getCoords().getX() < minX)
+ minX = _listeBases.get(i).getCoords().getX();
+ if (_listeBases.get(i).getCoords().getY() < minY)
+ minY = _listeBases.get(i).getCoords().getY();
+ if (_listeBases.get(i).getCoords().getX() > maxX)
+ maxX = _listeBases.get(i).getCoords().getX();
+ if (_listeBases.get(i).getCoords().getX() > maxY)
+ maxY = _listeBases.get(i).getCoords().getY();
+ }
+ // creation du point central
+ Point2D.Double centre = new Point2D.Double((maxX - minX) / 2,
+ (maxY - minY) / 2);
+ Double x, y;
+ for (int i = 0; i < _listeBases.size(); i++) {
+ // application de la rotation au centre de chaque base
+ // x' = cos(theta)*(x-xc) - sin(theta)*(y-yc) + xc
+ x = Math.cos(angle)
+ * (_listeBases.get(i).getCenter().getX() - centre.x)
+ - Math.sin(angle)
+ * (_listeBases.get(i).getCenter().getY() - centre.y)
+ + centre.x;
+ // y' = sin(theta)*(x-xc) + cos(theta)*(y-yc) + yc
+ y = Math.sin(angle)
+ * (_listeBases.get(i).getCenter().getX() - centre.x)
+ + Math.cos(angle)
+ * (_listeBases.get(i).getCenter().getY() - centre.y)
+ + centre.y;
+ _listeBases.get(i).setCenter(new Point2D.Double(x, y));
+
+ // application de la rotation au coordonnees de chaque
+ // base
+ // x' = cos(theta)*(x-xc) - sin(theta)*(y-yc) + xc
+ x = Math.cos(angle)
+ * (_listeBases.get(i).getCoords().getX() - centre.x)
+ - Math.sin(angle)
+ * (_listeBases.get(i).getCoords().getY() - centre.y)
+ + centre.x;
+ // y' = sin(theta)*(x-xc) + cos(theta)*(y-yc) + yc
+ y = Math.sin(angle)
+ * (_listeBases.get(i).getCoords().getX() - centre.x)
+ + Math.cos(angle)
+ * (_listeBases.get(i).getCoords().getY() - centre.y)
+ + centre.y;
+ _listeBases.get(i).setCoords(new Point2D.Double(x, y));
+ }
+ }
+ }
+
+ private static double MIN_DISTANCE = 10.;
+
+
+ /**
+ * Flip an helix around its supporting base
+ */
+ public void flipHelix(Point h) {
+ if (h.x!=-1 && h.y!=-1 && h.x!=h.y)
+ {
+ int hBeg=h.x;
+ int hEnd=h.y;
+ Point2D.Double A = getCoords(hBeg);
+ Point2D.Double B = getCoords(hEnd);
+ Point2D.Double AB = new Point2D.Double(B.x - A.x, B.y - A.y);
+ double normAB = Math.sqrt(AB.x * AB.x + AB.y * AB.y);
+ // Creating a coordinate system centered on A and having
+ // unit x-vector Ox.
+ Point2D.Double O = A;
+ Point2D.Double Ox = new Point2D.Double(AB.x / normAB, AB.y / normAB);
+ Hashtable<Integer,Point2D.Double> old = new Hashtable<Integer,Point2D.Double>();
+ for (int i = hBeg + 1; i < hEnd; i++) {
+ Point2D.Double P = getCoords(i);
+ Point2D.Double nP = project(O, Ox, P);
+ old.put(i, nP);
+ setCoord(i, nP);
+ Point2D.Double Center = getCenter(i);
+ setCenter(i, project(O, Ox, Center));
+ }
+ }
+ }
+
+ public static Point2D.Double project(Point2D.Double O, Point2D.Double Ox,
+ Point2D.Double C) {
+ Point2D.Double OC = new Point2D.Double(C.x - O.x, C.y - O.y);
+ // Projection of OC on OI => OX
+ double normOX = (Ox.x * OC.x + Ox.y * OC.y);
+ Point2D.Double OX = new Point2D.Double((normOX * Ox.x), (normOX * Ox.y));
+ // Portion of OC orthogonal to Ox => XC
+ Point2D.Double XC = new Point2D.Double(OC.x - OX.x, OC.y - OX.y);
+ // Reflexive image of C with respect to Ox => CP
+ Point2D.Double OCP = new Point2D.Double(OX.x - XC.x, OX.y - XC.y);
+ Point2D.Double CP = new Point2D.Double(O.x + OCP.x, O.y + OCP.y);
+ return CP;
+ }
+
+
+ public boolean testDirectionality(int i, int j, int k) {
+
+ // Which direction are we heading toward?
+ Point2D.Double pi = getCoords(i);
+ Point2D.Double pj = getCoords(j);
+ Point2D.Double pk = getCoords(k);
+ return testDirectionality(pi, pj, pk);
+ }
+
+ public static boolean testDirectionality(Point2D.Double pi,
+ Point2D.Double pj, Point2D.Double pk) {
+
+ // Which direction are we heading toward?
+ double test = (pj.x - pi.x) * (pk.y - pj.y) - (pj.y - pi.y)
+ * (pk.x - pj.x);
+ return test < 0.0;
+ }
+
+ public double getOrientation() {
+ double maxDist = Double.MIN_VALUE;
+ double angle = 0;
+ for (int i = 0; i < _listeBases.size(); i++) {
+ ModeleBase b1 = _listeBases.get(i);
+ for (int j = i + 1; j < _listeBases.size(); j++) {
+ ModeleBase b2 = _listeBases.get(j);
+ Point2D.Double p1 = b1._coords.toPoint2D();
+ Point2D.Double p2 = b2._coords.toPoint2D();
+ double dist = p1.distance(p2);
+ if (dist > maxDist) {
+ maxDist = dist;
+ angle = computeAngle(p1, p2);
+ }
+ }
+ }
+ return angle;
+ }
+
+ public boolean hasVirtualLoops() {
+ boolean consecutiveBPs = false;
+ for (int i = 0; i < _listeBases.size(); i++) {
+ int j = _listeBases.get(i).getElementStructure();
+ if (j == i + 1) {
+ consecutiveBPs = true;
+ }
+
+ }
+ return ((_drawMode != DRAW_MODE_LINEAR)
+ && (_drawMode != DRAW_MODE_CIRCULAR) && (consecutiveBPs));
+ }
+
+ public String getHTMLDescription() {
+ String result = "<table>";
+ result += "<tr><td><b>Name:</b></td><td>" + this._name + "</td></tr>";
+ result += "<tr><td><b>Length:</b></td><td>" + this.getSize()
+ + " nts</td></tr>";
+ result += "<tr><td><b>Base-pairs:</b></td><td>"
+ + this.getAllBPs().size() + " </td></tr>";
+ return result + "</table>";
+ }
+
+ public String getID() {
+ return _id;
+ }
+
+ public void setID(String id) {
+ _id = id;
+ }
+
+ public static ArrayList<Integer> getGapPositions(String gapString) {
+ ArrayList<Integer> result = new ArrayList<Integer>();
+ for (int i = 0; i < gapString.length(); i++) {
+ char c = gapString.charAt(i);
+ if (c == '.' || c == ':') {
+ result.add(i);
+ }
+ }
+ return result;
+ }
+
+ public RNA restrictTo(String gapString) {
+ return restrictTo(getGapPositions(gapString));
+ }
+
+ public RNA restrictTo(ArrayList<Integer> positions) {
+ RNA result = new RNA();
+ String oldSeq = this.getSeq();
+ String newSeq = "";
+ HashSet<Integer> removedPos = new HashSet<Integer>(positions);
+ int[] matching = new int[oldSeq.length()];
+ int j = 0;
+ for (int i = 0; i < oldSeq.length(); i++) {
+ matching[i] = j;
+ if (!removedPos.contains(i)) {
+ newSeq += oldSeq.charAt(i);
+ j++;
+ }
+ }
+ result.setRNA(newSeq);
+ for (ModeleBP m : getAllBPs()) {
+ if (removedPos.contains(m.getIndex5())
+ || removedPos.contains(m.getIndex3())) {
+ int i5 = matching[m.getIndex5()];
+ int i3 = matching[m.getIndex3()];
+ ModeleBP msbp = new ModeleBP(result.getBaseAt(i5),
+ result.getBaseAt(i3), m.getEdgePartner5(),
+ m.getEdgePartner3(), m.getStericity());
+ result.addBP(i5, i3, msbp);
+ }
+ }
+ return result;
+ }
+
+ public void rescale(double d) {
+ for (ModeleBase mb : _listeBases) {
+ mb._coords.x *= d;
+ mb._coords.y *= d;
+ mb._center.x *= d;
+ mb._center.y *= d;
+ }
+ }
+
+ /**
+ * Necessary for DrawRNATemplate (which is why the method is
+ * package-visible).
+ */
+ ArrayList<ModeleBase> getListeBases() {
+ return _listeBases;
+ }
+}