X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=src%2Fjalview%2Fgui%2FPCAPanel.java;h=772e67e8b8fccb12aa41810cc8e1af3c81e18e63;hb=442b6bc246ee877d5e977baf80fd96299e0e05c4;hp=d7e0cd78d3c7bea344592aa99589c2570127f426;hpb=c6032953b9b83b438ca8a1661d6b52a9a1950e23;p=jalview.git diff --git a/src/jalview/gui/PCAPanel.java b/src/jalview/gui/PCAPanel.java index d7e0cd7..772e67e 100755 --- a/src/jalview/gui/PCAPanel.java +++ b/src/jalview/gui/PCAPanel.java @@ -1,6 +1,6 @@ /* * Jalview - A Sequence Alignment Editor and Viewer - * Copyright (C) 2006 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle + * Copyright (C) 2007 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle * * This program is free software; you can redistribute it and/or * modify it under the terms of the GNU General Public License @@ -18,19 +18,16 @@ */ package jalview.gui; -import jalview.analysis.*; - -import jalview.datamodel.*; - -import jalview.jbgui.*; +import java.util.*; import java.awt.*; import java.awt.event.*; - -import java.util.*; -import javax.swing.*; import java.awt.print.*; +import javax.swing.*; +import jalview.analysis.*; +import jalview.datamodel.*; +import jalview.jbgui.*; /** * DOCUMENT ME! @@ -38,450 +35,450 @@ import java.awt.print.*; * @author $author$ * @version $Revision$ */ -public class PCAPanel extends GPCAPanel implements Runnable +public class PCAPanel + extends GPCAPanel implements Runnable { - PCA pca; - int top; - RotatableCanvas rc; - AlignmentPanel ap; - AlignViewport av; - AlignmentView seqstrings; - SequenceI [] seqs; - - /** - * Creates a new PCAPanel object. - * - * @param av DOCUMENT ME! - * @param s DOCUMENT ME! - */ - public PCAPanel(AlignmentPanel ap) + PCA pca; + int top; + RotatableCanvas rc; + AlignmentPanel ap; + AlignViewport av; + AlignmentView seqstrings; + SequenceI[] seqs; + + /** + * Creates a new PCAPanel object. + * + * @param av DOCUMENT ME! + * @param s DOCUMENT ME! + */ + public PCAPanel(AlignmentPanel ap) + { + this.av = ap.av; + this.ap = ap; + + boolean sameLength = true; + + seqstrings = av.getAlignmentView(av.getSelectionGroup() != null); + if (av.getSelectionGroup() == null) { - this.av = ap.av; - this.ap = ap; - - boolean sameLength = true; - - seqstrings = av.getAlignmentView(av.getSelectionGroup()!=null); - if(av.getSelectionGroup()==null) - { - seqs = av.alignment.getSequencesArray(); - } - else - { - seqs = av.getSelectionGroup().getSequencesInOrder(av.alignment); - } - SeqCigar sq[]=seqstrings.getSequences(); - int length = sq[0].getWidth(); - - for (int i = 0; i < seqs.length; i++) - { - if (sq[i].getWidth() != length) - { - sameLength = false; - break; - } - } - - if (!sameLength) - { - JOptionPane.showMessageDialog(Desktop.desktop, - "The sequences must be aligned before calculating PCA.\n" + - "Try using the Pad function in the edit menu,\n" + - "or one of the multiple sequence alignment web services.", - "Sequences not aligned", - JOptionPane.WARNING_MESSAGE); - - return; - } + seqs = av.alignment.getSequencesArray(); + } + else + { + seqs = av.getSelectionGroup().getSequencesInOrder(av.alignment); + } + SeqCigar sq[] = seqstrings.getSequences(); + int length = sq[0].getWidth(); + for (int i = 0; i < seqs.length; i++) + { + if (sq[i].getWidth() != length) + { + sameLength = false; + break; + } + } - Desktop.addInternalFrame(this, "Principal component analysis", - 400, 400); + if (!sameLength) + { + JOptionPane.showMessageDialog(Desktop.desktop, + "The sequences must be aligned before calculating PCA.\n" + + "Try using the Pad function in the edit menu,\n" + + "or one of the multiple sequence alignment web services.", + "Sequences not aligned", + JOptionPane.WARNING_MESSAGE); + + return; + } + Desktop.addInternalFrame(this, "Principal component analysis", + 400, 400); + PaintRefresher.Register(this, av.getSequenceSetId()); - PaintRefresher.Register(this, av.getSequenceSetId()); + rc = new RotatableCanvas(ap); + this.getContentPane().add(rc, BorderLayout.CENTER); + Thread worker = new Thread(this); + worker.start(); + } - rc = new RotatableCanvas(ap); - this.getContentPane().add(rc, BorderLayout.CENTER); - Thread worker = new Thread(this); - worker.start(); - } + public void bgcolour_actionPerformed(ActionEvent e) + { + Color col = JColorChooser.showDialog(this, "Select Background Colour", + rc.bgColour); - public void bgcolour_actionPerformed(ActionEvent e) + if (col != null) { - Color col = JColorChooser.showDialog(this, "Select Background Colour", - rc.bgColour); - - if(col!=null) - rc.bgColour = col; - rc.repaint(); + rc.bgColour = col; } + rc.repaint(); + } - - - - /** - * DOCUMENT ME! - */ - public void run() + /** + * DOCUMENT ME! + */ + public void run() + { + try { - try{ - pca = new PCA(seqstrings.getSequenceStrings(' ')); - pca.run(); + pca = new PCA(seqstrings.getSequenceStrings(' ')); + pca.run(); - // Now find the component coordinates - int ii = 0; + // Now find the component coordinates + int ii = 0; - while ( (ii < seqs.length) && (seqs[ii] != null)) - { - ii++; - } + while ( (ii < seqs.length) && (seqs[ii] != null)) + { + ii++; + } - double[][] comps = new double[ii][ii]; + double[][] comps = new double[ii][ii]; - for (int i = 0; i < ii; i++) + for (int i = 0; i < ii; i++) + { + if (pca.getEigenvalue(i) > 1e-4) { - if (pca.getEigenvalue(i) > 1e-4) - { - comps[i] = pca.component(i); - } + comps[i] = pca.component(i); } + } - ////////////////// - xCombobox.setSelectedIndex(0); - yCombobox.setSelectedIndex(1); - zCombobox.setSelectedIndex(2); - - top = pca.getM().rows - 1; - - Vector points = new Vector(); - float[][] scores = pca.getComponents(top - 1, top - 2, top - 3, 100); - - for (int i = 0; i < pca.getM().rows; i++) - { - SequencePoint sp = new SequencePoint(seqs[i], scores[i]); - points.addElement(sp); - } + ////////////////// + xCombobox.setSelectedIndex(0); + yCombobox.setSelectedIndex(1); + zCombobox.setSelectedIndex(2); - rc.setPoints(points, pca.getM().rows); - rc.repaint(); + top = pca.getM().rows - 1; - addKeyListener(rc); + Vector points = new Vector(); + float[][] scores = pca.getComponents(top - 1, top - 2, top - 3, 100); + for (int i = 0; i < pca.getM().rows; i++) + { + SequencePoint sp = new SequencePoint(seqs[i], scores[i]); + points.addElement(sp); } - catch(OutOfMemoryError er) - { JOptionPane.showInternalMessageDialog(Desktop.desktop, - "Out of memory calculating PCA!!" - + - "\nSee help files for increasing Java Virtual Machine memory." - , "Out of memory", - JOptionPane.WARNING_MESSAGE); - System.out.println("PCAPanel: "+er); - System.gc(); - } + rc.setPoints(points, pca.getM().rows); + rc.repaint(); - } + addKeyListener(rc); - /** - * DOCUMENT ME! - */ - void doDimensionChange() + } + catch (OutOfMemoryError er) { - if (top == 0) - { - return; - } - - int dim1 = top - xCombobox.getSelectedIndex(); - int dim2 = top - yCombobox.getSelectedIndex(); - int dim3 = top - zCombobox.getSelectedIndex(); + JOptionPane.showInternalMessageDialog(Desktop.desktop, + "Out of memory calculating PCA!!" + + + "\nSee help files for increasing Java Virtual Machine memory." + , "Out of memory", + JOptionPane.WARNING_MESSAGE); + System.out.println("PCAPanel: " + er); + System.gc(); - float[][] scores = pca.getComponents(dim1, dim2, dim3, 100); - - for (int i = 0; i < pca.getM().rows; i++) - { - ((SequencePoint) rc.points.elementAt(i)).coord = scores[i]; - } - - rc.img = null; - rc.rotmat.setIdentity(); - rc.initAxes(); - rc.paint(rc.getGraphics()); } - /** - * DOCUMENT ME! - * - * @param e DOCUMENT ME! - */ - protected void xCombobox_actionPerformed(ActionEvent e) - { - doDimensionChange(); - } + } - /** - * DOCUMENT ME! - * - * @param e DOCUMENT ME! - */ - protected void yCombobox_actionPerformed(ActionEvent e) + /** + * DOCUMENT ME! + */ + void doDimensionChange() + { + if (top == 0) { - doDimensionChange(); + return; } - /** - * DOCUMENT ME! - * - * @param e DOCUMENT ME! - */ - protected void zCombobox_actionPerformed(ActionEvent e) + int dim1 = top - xCombobox.getSelectedIndex(); + int dim2 = top - yCombobox.getSelectedIndex(); + int dim3 = top - zCombobox.getSelectedIndex(); + + float[][] scores = pca.getComponents(dim1, dim2, dim3, 100); + + for (int i = 0; i < pca.getM().rows; i++) { - doDimensionChange(); + ( (SequencePoint) rc.points.elementAt(i)).coord = scores[i]; } + rc.img = null; + rc.rotmat.setIdentity(); + rc.initAxes(); + rc.paint(rc.getGraphics()); + } - public void outputValues_actionPerformed(ActionEvent e) - { - CutAndPasteTransfer cap = new CutAndPasteTransfer(); - Desktop.addInternalFrame(cap, "PCA details", 500, - 500); + /** + * DOCUMENT ME! + * + * @param e DOCUMENT ME! + */ + protected void xCombobox_actionPerformed(ActionEvent e) + { + doDimensionChange(); + } - cap.setText(pca.getDetails()); - } + /** + * DOCUMENT ME! + * + * @param e DOCUMENT ME! + */ + protected void yCombobox_actionPerformed(ActionEvent e) + { + doDimensionChange(); + } - public void showLabels_actionPerformed(ActionEvent e) - { - rc.showLabels(showLabels.getState()); - } + /** + * DOCUMENT ME! + * + * @param e DOCUMENT ME! + */ + protected void zCombobox_actionPerformed(ActionEvent e) + { + doDimensionChange(); + } - public void print_actionPerformed(ActionEvent e) + public void outputValues_actionPerformed(ActionEvent e) + { + CutAndPasteTransfer cap = new CutAndPasteTransfer(); + Desktop.addInternalFrame(cap, "PCA details", 500, + 500); + + cap.setText(pca.getDetails()); + } + + public void showLabels_actionPerformed(ActionEvent e) + { + rc.showLabels(showLabels.getState()); + } + + public void print_actionPerformed(ActionEvent e) + { + PCAPrinter printer = new PCAPrinter(); + printer.start(); + } + + public void originalSeqData_actionPerformed(ActionEvent e) + { + // this was cut'n'pasted from the equivalent TreePanel method - we should make this an abstract function of all jalview analysis windows + if (seqstrings == null) { - PCAPrinter printer = new PCAPrinter(); - printer.start(); + jalview.bin.Cache.log.info("Unexpected call to originalSeqData_actionPerformed - should have hidden this menu action."); + return; } - - public void originalSeqData_actionPerformed(ActionEvent e) + // decide if av alignment is sufficiently different to original data to warrant a new window to be created + // create new alignmnt window with hidden regions (unhiding hidden regions yields unaligned seqs) + // or create a selection box around columns in alignment view + // test Alignment(SeqCigar[]) + Object[] alAndColsel = seqstrings.getAlignmentAndColumnSelection(av. + getGapCharacter()); + + if (alAndColsel != null && alAndColsel[0] != null) { - // this was cut'n'pasted from the equivalent TreePanel method - we should make this an abstract function of all jalview analysis windows - if (seqstrings==null) + // AlignmentOrder origorder = new AlignmentOrder(alAndColsel[0]); + + Alignment al = new Alignment( (SequenceI[]) alAndColsel[0]); + Alignment dataset = av.getAlignment().getDataset(); + if (dataset != null) { - jalview.bin.Cache.log.info("Unexpected call to originalSeqData_actionPerformed - should have hidden this menu action."); - return; + al.setDataset(dataset); } - // decide if av alignment is sufficiently different to original data to warrant a new window to be created - // create new alignmnt window with hidden regions (unhiding hidden regions yields unaligned seqs) - // or create a selection box around columns in alignment view - // test Alignment(SeqCigar[]) - Object[] alAndColsel = seqstrings.getAlignmentAndColumnSelection(av. - getGapCharacter()); - - if (alAndColsel != null && alAndColsel[0]!=null) + if (true) { - // AlignmentOrder origorder = new AlignmentOrder(alAndColsel[0]); - - Alignment al = new Alignment((SequenceI[]) alAndColsel[0]); - Alignment dataset = av.getAlignment().getDataset(); - if (dataset != null) - { - al.setDataset(dataset); - } - - if (true) - { - // make a new frame! - AlignFrame af = new AlignFrame(al, (ColumnSelection) alAndColsel[1], - AlignFrame.DEFAULT_WIDTH, - AlignFrame.DEFAULT_HEIGHT -); - - //>>>This is a fix for the moment, until a better solution is found!!<<< - // af.getFeatureRenderer().transferSettings(alignFrame.getFeatureRenderer()); - - // af.addSortByOrderMenuItem(ServiceName + " Ordering", - // msaorder); - - Desktop.addInternalFrame(af, "Original Data for " + this.title, + // make a new frame! + AlignFrame af = new AlignFrame(al, (ColumnSelection) alAndColsel[1], AlignFrame.DEFAULT_WIDTH, - AlignFrame.DEFAULT_HEIGHT); - } - } - /* CutAndPasteTransfer cap = new CutAndPasteTransfer(); - for (int i = 0; i < seqs.length; i++) - { - cap.appendText(new jalview.util.Format("%-" + 15 + "s").form( - seqs[i].getName())); - cap.appendText(" " + seqstrings[i] + "\n"); - - } - - Desktop.addInternalFrame(cap, "Original Data", - 400, 400); - */ + AlignFrame.DEFAULT_HEIGHT + ); + + //>>>This is a fix for the moment, until a better solution is found!!<<< + // af.getFeatureRenderer().transferSettings(alignFrame.getFeatureRenderer()); + + // af.addSortByOrderMenuItem(ServiceName + " Ordering", + // msaorder); + + Desktop.addInternalFrame(af, "Original Data for " + this.title, + AlignFrame.DEFAULT_WIDTH, + AlignFrame.DEFAULT_HEIGHT); + } } + /* CutAndPasteTransfer cap = new CutAndPasteTransfer(); + for (int i = 0; i < seqs.length; i++) + { + cap.appendText(new jalview.util.Format("%-" + 15 + "s").form( + seqs[i].getName())); + cap.appendText(" " + seqstrings[i] + "\n"); + } + Desktop.addInternalFrame(cap, "Original Data", + 400, 400); + */ + } - class PCAPrinter extends Thread implements Printable + class PCAPrinter + extends Thread implements Printable + { + public void run() { - public void run() - { - PrinterJob printJob = PrinterJob.getPrinterJob(); - PageFormat pf = printJob.pageDialog(printJob.defaultPage()); + PrinterJob printJob = PrinterJob.getPrinterJob(); + PageFormat pf = printJob.pageDialog(printJob.defaultPage()); - printJob.setPrintable(this, pf); + printJob.setPrintable(this, pf); - if (printJob.printDialog()) + if (printJob.printDialog()) + { + try + { + printJob.print(); + } + catch (Exception PrintException) { - try - { - printJob.print(); - } - catch (Exception PrintException) - { - PrintException.printStackTrace(); - } + PrintException.printStackTrace(); } } + } - public int print(Graphics pg, PageFormat pf, int pi) - throws PrinterException - { - pg.translate( (int) pf.getImageableX(), (int) pf.getImageableY()); + public int print(Graphics pg, PageFormat pf, int pi) + throws PrinterException + { + pg.translate( (int) pf.getImageableX(), (int) pf.getImageableY()); - rc.drawBackground(pg, rc.bgColour); - rc.drawScene(pg); - if (rc.drawAxes == true) - { - rc.drawAxes(pg); - } + rc.drawBackground(pg, rc.bgColour); + rc.drawScene(pg); + if (rc.drawAxes == true) + { + rc.drawAxes(pg); + } - if (pi == 0) - return Printable.PAGE_EXISTS; - else - return Printable.NO_SUCH_PAGE; + if (pi == 0) + { + return Printable.PAGE_EXISTS; + } + else + { + return Printable.NO_SUCH_PAGE; } } + } + /** + * DOCUMENT ME! + * + * @param e DOCUMENT ME! + */ + public void eps_actionPerformed(ActionEvent e) + { + makePCAImage(jalview.util.ImageMaker.EPS); + } + /** + * DOCUMENT ME! + * + * @param e DOCUMENT ME! + */ + public void png_actionPerformed(ActionEvent e) + { + makePCAImage(jalview.util.ImageMaker.PNG); + } - /** - * DOCUMENT ME! - * - * @param e DOCUMENT ME! - */ - public void eps_actionPerformed(ActionEvent e) - { - makePCAImage(jalview.util.ImageMaker.EPS); - } + void makePCAImage(int type) + { + int width = rc.getWidth(); + int height = rc.getHeight(); - /** - * DOCUMENT ME! - * - * @param e DOCUMENT ME! - */ - public void png_actionPerformed(ActionEvent e) - { - makePCAImage(jalview.util.ImageMaker.PNG); - } + jalview.util.ImageMaker im; - void makePCAImage(int type) + if (type == jalview.util.ImageMaker.PNG) { - int width = rc.getWidth(); - int height = rc.getHeight(); - - jalview.util.ImageMaker im; - - if(type == jalview.util.ImageMaker.PNG) - im = new jalview.util.ImageMaker(this, - jalview.util.ImageMaker.PNG, - "Make PNG image from PCA", - width, height, - null, null); - else - im = new jalview.util.ImageMaker(this, - jalview.util.ImageMaker.EPS, - "Make EPS file from PCA", - width, height, - null, this.getTitle()); - - if(im.getGraphics()!=null) - { - rc.drawBackground(im.getGraphics(), Color.black); - rc.drawScene(im.getGraphics()); - if (rc.drawAxes == true) - { - rc.drawAxes(im.getGraphics()); - } - im.writeImage(); - } + im = new jalview.util.ImageMaker(this, + jalview.util.ImageMaker.PNG, + "Make PNG image from PCA", + width, height, + null, null); } - - - public void viewMenu_menuSelected() + else { - buildAssociatedViewMenu(); + im = new jalview.util.ImageMaker(this, + jalview.util.ImageMaker.EPS, + "Make EPS file from PCA", + width, height, + null, this.getTitle()); } - - void buildAssociatedViewMenu() + if (im.getGraphics() != null) { - AlignmentPanel [] aps = PaintRefresher.getAssociatedPanels(av.getSequenceSetId()); - if(aps.length==1 && rc.av ==aps[0].av) + rc.drawBackground(im.getGraphics(), Color.black); + rc.drawScene(im.getGraphics()); + if (rc.drawAxes == true) { - associateViewsMenu.setVisible(false); - return; - } - - associateViewsMenu.setVisible(true); - - if( (viewMenu.getItem(viewMenu.getItemCount()-2) instanceof JMenuItem)) - { - viewMenu.insertSeparator(viewMenu.getItemCount()-1); + rc.drawAxes(im.getGraphics()); } + im.writeImage(); + } + } - associateViewsMenu.removeAll(); + public void viewMenu_menuSelected() + { + buildAssociatedViewMenu(); + } + void buildAssociatedViewMenu() + { + AlignmentPanel[] aps = PaintRefresher.getAssociatedPanels(av. + getSequenceSetId()); + if (aps.length == 1 && rc.av == aps[0].av) + { + associateViewsMenu.setVisible(false); + return; + } - JRadioButtonMenuItem item; - ButtonGroup buttonGroup = new ButtonGroup(); - int i, iSize = aps.length; - final PCAPanel thisPCAPanel = this; - for(i=0; i