/* * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$) * Copyright (C) $$Year-Rel$$ The Jalview Authors * * This file is part of Jalview. * * Jalview is free software: you can redistribute it and/or * modify it under the terms of the GNU General Public License * as published by the Free Software Foundation, either version 3 * of the License, or (at your option) any later version. * * Jalview is distributed in the hope that it will be useful, but * WITHOUT ANY WARRANTY; without even the implied warranty * of MERCHANTABILITY or FITNESS FOR A PARTICULAR * PURPOSE. See the GNU General Public License for more details. * * You should have received a copy of the GNU General Public License * along with Jalview. If not, see . * The Jalview Authors are detailed in the 'AUTHORS' file. */ package jalview.gui; import jalview.analysis.scoremodels.ScoreModels; import jalview.analysis.scoremodels.SimilarityParams; import jalview.api.analysis.ScoreModelI; import jalview.api.analysis.SimilarityParamsI; import jalview.datamodel.Alignment; import jalview.datamodel.AlignmentI; import jalview.datamodel.AlignmentView; import jalview.datamodel.HiddenColumns; import jalview.datamodel.SequenceI; import jalview.jbgui.GPCAPanel; import jalview.util.MessageManager; import jalview.viewmodel.AlignmentViewport; import jalview.viewmodel.PCAModel; import java.awt.BorderLayout; import java.awt.Color; import java.awt.Dimension; import java.awt.Graphics; import java.awt.event.ActionEvent; import java.awt.event.ActionListener; import java.awt.print.PageFormat; import java.awt.print.Printable; import java.awt.print.PrinterException; import java.awt.print.PrinterJob; import javax.swing.ButtonGroup; import javax.swing.JCheckBoxMenuItem; import javax.swing.JColorChooser; import javax.swing.JMenuItem; import javax.swing.JRadioButtonMenuItem; import javax.swing.event.InternalFrameAdapter; import javax.swing.event.InternalFrameEvent; /** * DOCUMENT ME! * * @author $author$ * @version $Revision$ */ public class PCAPanel extends GPCAPanel implements Runnable, IProgressIndicator { private IProgressIndicator progressBar; RotatableCanvas rc; AlignmentPanel ap; AlignmentViewport av; PCAModel pcaModel; private static final int MIN_WIDTH = 470; private static final int MIN_HEIGHT = 250; int top = 0; private boolean working; /** * Creates a new PCAPanel object using default score model and parameters * * @param alignPanel */ public PCAPanel(AlignmentPanel alignPanel) { this(alignPanel, ScoreModels.getInstance() .getDefaultModel( !alignPanel.av.getAlignment().isNucleotide()) .getName(), SimilarityParams.SeqSpace); } /** * Constructor given sequence data, a similarity (or distance) score model * name, and score calculation parameters * * @param alignPanel * @param modelName * @param params */ public PCAPanel(AlignmentPanel alignPanel, String modelName, SimilarityParamsI params) { super(); this.av = alignPanel.av; this.ap = alignPanel; boolean nucleotide = av.getAlignment().isNucleotide(); progressBar = new ProgressBar(statusPanel, statusBar); addInternalFrameListener(new InternalFrameAdapter() { @Override public void internalFrameClosed(InternalFrameEvent e) { close_actionPerformed(); } }); boolean selected = av.getSelectionGroup() != null && av.getSelectionGroup().getSize() > 0; AlignmentView seqstrings = av.getAlignmentView(selected); SequenceI[] seqs; if (!selected) { seqs = av.getAlignment().getSequencesArray(); } else { seqs = av.getSelectionGroup().getSequencesInOrder(av.getAlignment()); } ScoreModelI scoreModel = ScoreModels.getInstance() .getScoreModel(modelName, ap); pcaModel = new PCAModel(seqstrings, seqs, nucleotide, scoreModel, params); PaintRefresher.Register(this, av.getSequenceSetId()); rc = new RotatableCanvas(alignPanel); this.getContentPane().add(rc, BorderLayout.CENTER); Thread worker = new Thread(this); worker.start(); } /** * Ensure references to potentially very large objects (the PCA matrices) are * nulled when the frame is closed */ protected void close_actionPerformed() { pcaModel = null; } /** * Repopulate the options and actions under the score model menu when it is * selected. Options will depend on whether 'nucleotide' or 'peptide' * modelling is selected (and also possibly on whether any additional score * models have been added). */ @Override protected void scoreModel_menuSelected() { scoreModelMenu.removeAll(); for (final ScoreModelI sm : ScoreModels.getInstance().getModels()) { final String name = sm.getName(); JCheckBoxMenuItem jm = new JCheckBoxMenuItem(name); /* * if the score model doesn't provide a description, try to look one * up in the text bundle, falling back on its name */ String tooltip = sm.getDescription(); if (tooltip == null) { tooltip = MessageManager.getStringOrReturn("label.score_model_", name); } jm.setToolTipText(tooltip); jm.setSelected(pcaModel.getScoreModelName().equals(name)); if ((pcaModel.isNucleotide() && sm.isDNA()) || (!pcaModel.isNucleotide() && sm.isProtein())) { jm.addActionListener(new ActionListener() { @Override public void actionPerformed(ActionEvent e) { if (!pcaModel.getScoreModelName().equals(name)) { ScoreModelI sm2 = ScoreModels.getInstance() .getScoreModel(name, ap); pcaModel.setScoreModel(sm2); Thread worker = new Thread(PCAPanel.this); worker.start(); } } }); scoreModelMenu.add(jm); } } } @Override public void bgcolour_actionPerformed(ActionEvent e) { Color col = JColorChooser.showDialog(this, MessageManager.getString("label.select_background_colour"), rc.bgColour); if (col != null) { rc.bgColour = col; } rc.repaint(); } /** * DOCUMENT ME! */ @Override public void run() { long progId = System.currentTimeMillis(); IProgressIndicator progress = this; String message = MessageManager.getString("label.pca_recalculating"); if (getParent() == null) { progress = ap.alignFrame; message = MessageManager.getString("label.pca_calculating"); } progress.setProgressBar(message, progId); working = true; try { calcSettings.setEnabled(false); pcaModel.run(); // //////////////// xCombobox.setSelectedIndex(0); yCombobox.setSelectedIndex(1); zCombobox.setSelectedIndex(2); pcaModel.updateRc(rc); // rc.invalidate(); nuclSetting.setSelected(pcaModel.isNucleotide()); protSetting.setSelected(!pcaModel.isNucleotide()); top = pcaModel.getTop(); } catch (OutOfMemoryError er) { new OOMWarning("calculating PCA", er); working = false; return; } finally { progress.setProgressBar("", progId); } calcSettings.setEnabled(true); repaint(); if (getParent() == null) { addKeyListener(rc); Desktop.addInternalFrame(this, MessageManager .getString("label.principal_component_analysis"), 475, 450); this.setMinimumSize(new Dimension(MIN_WIDTH, MIN_HEIGHT)); } working = false; } @Override protected void nuclSetting_actionPerfomed(ActionEvent arg0) { if (!pcaModel.isNucleotide()) { pcaModel.setNucleotide(true); pcaModel.setScoreModel( ScoreModels.getInstance().getDefaultModel(false)); Thread worker = new Thread(this); worker.start(); } } @Override protected void protSetting_actionPerfomed(ActionEvent arg0) { if (pcaModel.isNucleotide()) { pcaModel.setNucleotide(false); pcaModel.setScoreModel( ScoreModels.getInstance().getDefaultModel(true)); Thread worker = new Thread(this); worker.start(); } } /** * DOCUMENT ME! */ void doDimensionChange() { if (top == 0) { return; } int dim1 = top - xCombobox.getSelectedIndex(); int dim2 = top - yCombobox.getSelectedIndex(); int dim3 = top - zCombobox.getSelectedIndex(); pcaModel.updateRcView(dim1, dim2, dim3); rc.img = null; rc.rotmat.setIdentity(); rc.initAxes(); rc.paint(rc.getGraphics()); } /** * DOCUMENT ME! * * @param e * DOCUMENT ME! */ @Override protected void xCombobox_actionPerformed(ActionEvent e) { doDimensionChange(); } /** * DOCUMENT ME! * * @param e * DOCUMENT ME! */ @Override protected void yCombobox_actionPerformed(ActionEvent e) { doDimensionChange(); } /** * DOCUMENT ME! * * @param e * DOCUMENT ME! */ @Override protected void zCombobox_actionPerformed(ActionEvent e) { doDimensionChange(); } @Override public void outputValues_actionPerformed(ActionEvent e) { CutAndPasteTransfer cap = new CutAndPasteTransfer(); try { cap.setText(pcaModel.getDetails()); Desktop.addInternalFrame(cap, MessageManager.getString("label.pca_details"), 500, 500); } catch (OutOfMemoryError oom) { new OOMWarning("opening PCA details", oom); cap.dispose(); } } @Override public void showLabels_actionPerformed(ActionEvent e) { rc.showLabels(showLabels.getState()); } @Override public void print_actionPerformed(ActionEvent e) { PCAPrinter printer = new PCAPrinter(); printer.start(); } @Override 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 (pcaModel.getSeqtrings() == null) { jalview.bin.Cache.log.info( "Unexpected call to originalSeqData_actionPerformed - should have hidden this menu action."); return; } // 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[]) char gc = '-'; try { // we try to get the associated view's gap character // but this may fail if the view was closed... gc = av.getGapCharacter(); } catch (Exception ex) { } ; Object[] alAndColsel = pcaModel.getSeqtrings() .getAlignmentAndHiddenColumns(gc); if (alAndColsel != null && alAndColsel[0] != null) { // AlignmentOrder origorder = new AlignmentOrder(alAndColsel[0]); AlignmentI al = new Alignment((SequenceI[]) alAndColsel[0]); AlignmentI dataset = (av != null && av.getAlignment() != null) ? av.getAlignment().getDataset() : null; if (dataset != null) { al.setDataset(dataset); } if (true) { // make a new frame! AlignFrame af = new AlignFrame(al, (HiddenColumns) 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, MessageManager.formatMessage( "label.original_data_for_params", new String[] { 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 { @Override public void run() { PrinterJob printJob = PrinterJob.getPrinterJob(); PageFormat pf = printJob.pageDialog(printJob.defaultPage()); printJob.setPrintable(this, pf); if (printJob.printDialog()) { try { printJob.print(); } catch (Exception PrintException) { PrintException.printStackTrace(); } } } @Override 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); } if (pi == 0) { return Printable.PAGE_EXISTS; } else { return Printable.NO_SUCH_PAGE; } } } /** * DOCUMENT ME! * * @param e * DOCUMENT ME! */ @Override public void eps_actionPerformed(ActionEvent e) { makePCAImage(jalview.util.ImageMaker.TYPE.EPS); } /** * DOCUMENT ME! * * @param e * DOCUMENT ME! */ @Override public void png_actionPerformed(ActionEvent e) { makePCAImage(jalview.util.ImageMaker.TYPE.PNG); } void makePCAImage(jalview.util.ImageMaker.TYPE type) { int width = rc.getWidth(); int height = rc.getHeight(); jalview.util.ImageMaker im; if (type == jalview.util.ImageMaker.TYPE.PNG) { im = new jalview.util.ImageMaker(this, jalview.util.ImageMaker.TYPE.PNG, "Make PNG image from PCA", width, height, null, null, null, 0, false); } else if (type == jalview.util.ImageMaker.TYPE.EPS) { im = new jalview.util.ImageMaker(this, jalview.util.ImageMaker.TYPE.EPS, "Make EPS file from PCA", width, height, null, this.getTitle(), null, 0, false); } else { im = new jalview.util.ImageMaker(this, jalview.util.ImageMaker.TYPE.SVG, "Make SVG file from PCA", width, height, null, this.getTitle(), null, 0, false); } if (im.getGraphics() != null) { rc.drawBackground(im.getGraphics(), Color.black); rc.drawScene(im.getGraphics()); if (rc.drawAxes == true) { rc.drawAxes(im.getGraphics()); } im.writeImage(); } } @Override 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; } associateViewsMenu.setVisible(true); if ((viewMenu .getItem(viewMenu.getItemCount() - 2) instanceof JMenuItem)) { viewMenu.insertSeparator(viewMenu.getItemCount() - 1); } associateViewsMenu.removeAll(); JRadioButtonMenuItem item; ButtonGroup buttonGroup = new ButtonGroup(); int i, iSize = aps.length; final PCAPanel thisPCAPanel = this; for (i = 0; i < iSize; i++) { final AlignmentPanel ap = aps[i]; item = new JRadioButtonMenuItem(ap.av.viewName, ap.av == rc.av); buttonGroup.add(item); item.addActionListener(new ActionListener() { @Override public void actionPerformed(ActionEvent evt) { rc.applyToAllViews = false; rc.av = ap.av; rc.ap = ap; PaintRefresher.Register(thisPCAPanel, ap.av.getSequenceSetId()); } }); associateViewsMenu.add(item); } final JRadioButtonMenuItem itemf = new JRadioButtonMenuItem( "All Views"); buttonGroup.add(itemf); itemf.setSelected(rc.applyToAllViews); itemf.addActionListener(new ActionListener() { @Override public void actionPerformed(ActionEvent evt) { rc.applyToAllViews = itemf.isSelected(); } }); associateViewsMenu.add(itemf); } /* * (non-Javadoc) * * @see * jalview.jbgui.GPCAPanel#outputPoints_actionPerformed(java.awt.event.ActionEvent * ) */ @Override protected void outputPoints_actionPerformed(ActionEvent e) { CutAndPasteTransfer cap = new CutAndPasteTransfer(); try { cap.setText(pcaModel.getPointsasCsv(false, xCombobox.getSelectedIndex(), yCombobox.getSelectedIndex(), zCombobox.getSelectedIndex())); Desktop.addInternalFrame(cap, MessageManager .formatMessage("label.points_for_params", new String[] { this.getTitle() }), 500, 500); } catch (OutOfMemoryError oom) { new OOMWarning("exporting PCA points", oom); cap.dispose(); } } /* * (non-Javadoc) * * @see * jalview.jbgui.GPCAPanel#outputProjPoints_actionPerformed(java.awt.event * .ActionEvent) */ @Override protected void outputProjPoints_actionPerformed(ActionEvent e) { CutAndPasteTransfer cap = new CutAndPasteTransfer(); try { cap.setText(pcaModel.getPointsasCsv(true, xCombobox.getSelectedIndex(), yCombobox.getSelectedIndex(), zCombobox.getSelectedIndex())); Desktop.addInternalFrame(cap, MessageManager.formatMessage( "label.transformed_points_for_params", new String[] { this.getTitle() }), 500, 500); } catch (OutOfMemoryError oom) { new OOMWarning("exporting transformed PCA points", oom); cap.dispose(); } } /* * (non-Javadoc) * * @see jalview.gui.IProgressIndicator#setProgressBar(java.lang.String, long) */ @Override public void setProgressBar(String message, long id) { progressBar.setProgressBar(message, id); // if (progressBars == null) // { // progressBars = new Hashtable(); // progressBarHandlers = new Hashtable(); // } // // JPanel progressPanel; // Long lId = new Long(id); // GridLayout layout = (GridLayout) statusPanel.getLayout(); // if (progressBars.get(lId) != null) // { // progressPanel = (JPanel) progressBars.get(new Long(id)); // statusPanel.remove(progressPanel); // progressBars.remove(lId); // progressPanel = null; // if (message != null) // { // statusBar.setText(message); // } // if (progressBarHandlers.contains(lId)) // { // progressBarHandlers.remove(lId); // } // layout.setRows(layout.getRows() - 1); // } // else // { // progressPanel = new JPanel(new BorderLayout(10, 5)); // // JProgressBar progressBar = new JProgressBar(); // progressBar.setIndeterminate(true); // // progressPanel.add(new JLabel(message), BorderLayout.WEST); // progressPanel.add(progressBar, BorderLayout.CENTER); // // layout.setRows(layout.getRows() + 1); // statusPanel.add(progressPanel); // // progressBars.put(lId, progressPanel); // } // // update GUI // // setMenusForViewport(); // validate(); } @Override public void registerHandler(final long id, final IProgressIndicatorHandler handler) { progressBar.registerHandler(id, handler); // if (progressBarHandlers == null || !progressBars.contains(new Long(id))) // { // throw new // Error(MessageManager.getString("error.call_setprogressbar_before_registering_handler")); // } // progressBarHandlers.put(new Long(id), handler); // final JPanel progressPanel = (JPanel) progressBars.get(new Long(id)); // if (handler.canCancel()) // { // JButton cancel = new JButton( // MessageManager.getString("action.cancel")); // final IProgressIndicator us = this; // cancel.addActionListener(new ActionListener() // { // // @Override // public void actionPerformed(ActionEvent e) // { // handler.cancelActivity(id); // us.setProgressBar(MessageManager.formatMessage("label.cancelled_params", // new String[]{((JLabel) progressPanel.getComponent(0)).getText()}), id); // } // }); // progressPanel.add(cancel, BorderLayout.EAST); // } } /** * * @return true if any progress bars are still active */ @Override public boolean operationInProgress() { return progressBar.operationInProgress(); } @Override protected void resetButton_actionPerformed(ActionEvent e) { int t = top; top = 0; // ugly - prevents dimensionChanged events from being processed xCombobox.setSelectedIndex(0); yCombobox.setSelectedIndex(1); top = t; zCombobox.setSelectedIndex(2); } /** * Answers true if PCA calculation is in progress, else false * * @return */ public boolean isWorking() { return working; } }