git.gag.com Git - debian/openrocket/blob - src/net/sf/openrocket/gui/dialogs/optimization/OptimizationPlotDialog.java

   1 package net.sf.openrocket.gui.dialogs.optimization;
   2
   3 import java.awt.Color;
   4 import java.awt.Paint;
   5 import java.awt.Window;
   6 import java.awt.event.ActionEvent;
   7 import java.awt.event.ActionListener;
   8 import java.util.ArrayList;
   9 import java.util.Iterator;
  10 import java.util.List;
  11 import java.util.Map;
  12
  13 import javax.swing.BorderFactory;
  14 import javax.swing.JButton;
  15 import javax.swing.JDialog;
  16 import javax.swing.JLabel;
  17 import javax.swing.JPanel;
  18
  19 import net.miginfocom.swing.MigLayout;
  20 import net.sf.openrocket.gui.components.StyledLabel;
  21 import net.sf.openrocket.l10n.Translator;
  22 import net.sf.openrocket.logging.LogHelper;
  23 import net.sf.openrocket.optimization.general.Point;
  24 import net.sf.openrocket.optimization.rocketoptimization.OptimizableParameter;
  25 import net.sf.openrocket.optimization.rocketoptimization.SimulationModifier;
  26 import net.sf.openrocket.startup.Application;
  27 import net.sf.openrocket.unit.Unit;
  28 import net.sf.openrocket.unit.UnitGroup;
  29 import net.sf.openrocket.unit.Value;
  30 import net.sf.openrocket.util.GUIUtil;
  31 import net.sf.openrocket.util.LinearInterpolator;
  32 import net.sf.openrocket.util.MathUtil;
  33
  34 import org.jfree.chart.ChartFactory;
  35 import org.jfree.chart.ChartPanel;
  36 import org.jfree.chart.JFreeChart;
  37 import org.jfree.chart.annotations.XYBoxAnnotation;
  38 import org.jfree.chart.annotations.XYLineAnnotation;
  39 import org.jfree.chart.annotations.XYPointerAnnotation;
  40 import org.jfree.chart.axis.AxisLocation;
  41 import org.jfree.chart.axis.NumberAxis;
  42 import org.jfree.chart.labels.CustomXYToolTipGenerator;
  43 import org.jfree.chart.plot.PlotOrientation;
  44 import org.jfree.chart.plot.XYPlot;
  45 import org.jfree.chart.renderer.PaintScale;
  46 import org.jfree.chart.renderer.xy.XYLineAndShapeRenderer;
  47 import org.jfree.chart.renderer.xy.XYShapeRenderer;
  48 import org.jfree.chart.title.PaintScaleLegend;
  49 import org.jfree.data.xy.DefaultXYZDataset;
  50 import org.jfree.data.xy.XYSeries;
  51 import org.jfree.data.xy.XYSeriesCollection;
  52 import org.jfree.ui.RectangleEdge;
  53 import org.jfree.ui.TextAnchor;
  54
  55 /**
  56  * A class that plots the path of an optimization.
  57  *
  58  * @author Sampo Niskanen <sampo.niskanen@iki.fi>
  59  */
  60 public class OptimizationPlotDialog extends JDialog {
  61         private static final LogHelper log = Application.getLogger();
  62         private static final Translator trans = Application.getTranslator();
  63
  64
  65         private static final LinearInterpolator RED = new LinearInterpolator(
  66                         new double[] { 0.0, 1.0 }, new double[] { 0.0, 1.0 }
  67                         );
  68         private static final LinearInterpolator GREEN = new LinearInterpolator(
  69                         new double[] { 0.0, 1.0 }, new double[] { 0.0, 0.0 }
  70                         );
  71         private static final LinearInterpolator BLUE = new LinearInterpolator(
  72                         new double[] { 0.0, 1.0 }, new double[] { 1.0, 0.0 }
  73                         );
  74
  75         private static final Color OUT_OF_DOMAIN_COLOR = Color.BLACK;
  76
  77         private static final Color PATH_COLOR = new Color(220, 0, 0);
  78
  79
  80         public OptimizationPlotDialog(List<Point> path, Map<Point, FunctionEvaluationData> evaluations,
  81                         List<SimulationModifier> modifiers, OptimizableParameter parameter, UnitGroup stabilityUnit, Window parent) {
  82                 super(parent, trans.get("title"), ModalityType.APPLICATION_MODAL);
  83
  84
  85                 JPanel panel = new JPanel(new MigLayout("fill"));
  86
  87                 ChartPanel chart;
  88                 if (modifiers.size() == 1) {
  89                         chart = create1DPlot(path, evaluations, modifiers, parameter, stabilityUnit);
  90                 } else if (modifiers.size() == 2) {
  91                         chart = create2DPlot(path, evaluations, modifiers, parameter, stabilityUnit);
  92                 } else {
  93                         throw new IllegalArgumentException("Invalid dimensionality, dim=" + modifiers.size());
  94                 }
  95                 chart.setBorder(BorderFactory.createLineBorder(Color.BLACK));
  96                 panel.add(chart, "span, grow, wrap para");
  97
  98
  99                 JLabel label = new StyledLabel(trans.get("lbl.zoomInstructions"), -2);
 100                 panel.add(label, "");
 101
 102
 103                 JButton close = new JButton(trans.get("button.close"));
 104                 close.addActionListener(new ActionListener() {
 105                         @Override
 106                         public void actionPerformed(ActionEvent e) {
 107                                 OptimizationPlotDialog.this.setVisible(false);
 108                         }
 109                 });
 110                 panel.add(close, "right");
 111
 112
 113                 this.add(panel);
 114
 115                 GUIUtil.setDisposableDialogOptions(this, close);
 116         }
 117
 118
 119
 120         /**
 121          * Create a 1D plot of the optimization path.
 122          */
 123         private ChartPanel create1DPlot(List<Point> path, Map<Point, FunctionEvaluationData> evaluations,
 124                         List<SimulationModifier> modifiers, OptimizableParameter parameter, UnitGroup stabilityUnit) {
 125
 126                 SimulationModifier modX = modifiers.get(0);
 127                 Unit xUnit = modX.getUnitGroup().getDefaultUnit();
 128                 Unit yUnit = parameter.getUnitGroup().getDefaultUnit();
 129
 130                 // Create the optimization path (with autosort)
 131                 XYSeries series = new XYSeries(trans.get("plot1d.series"), true, true);
 132                 List<String> tooltips = new ArrayList<String>();
 133                 for (Point p : evaluations.keySet()) {
 134                         FunctionEvaluationData data = evaluations.get(p);
 135                         if (data != null) {
 136                                 if (data.getParameterValue() != null) {
 137                                         Value[] state = data.getState();
 138                                         series.add(xUnit.toUnit(state[0].getValue()), yUnit.toUnit(data.getParameterValue().getValue()));
 139                                         tooltips.add(getTooltip(data, parameter));
 140                                 }
 141                         } else {
 142                                 log.error("Could not find evaluation data for point " + p);
 143                         }
 144                 }
 145
 146
 147                 String xLabel = modX.getRelatedObject().toString() + ": " + modX.getName() + " / " + xUnit.getUnit();
 148                 String yLabel = parameter.getName() + " / " + yUnit.getUnit();
 149
 150                 JFreeChart chart = ChartFactory.createXYLineChart(
 151                                 trans.get("plot1d.title"),
 152                                 xLabel,
 153                                 yLabel,
 154                                 null,
 155                                 PlotOrientation.VERTICAL,
 156                                 false, // Legend
 157                                 true, // Tooltips
 158                                 false); // Urls
 159
 160
 161                 // Set the scale of the plot to the limits
 162                 double x1 = xUnit.toUnit(modX.getMinValue());
 163                 double x2 = xUnit.toUnit(modX.getMaxValue());
 164
 165                 chart.getXYPlot().getDomainAxis().setRange(x1, x2);
 166
 167                 // Add lines to show optimization limits
 168                 XYLineAnnotation line = new XYLineAnnotation(x1, -1e19, x1, 1e19);
 169                 chart.getXYPlot().addAnnotation(line);
 170                 line = new XYLineAnnotation(x2, -1e19, x2, 1e19);
 171                 chart.getXYPlot().addAnnotation(line);
 172
 173                 // Mark the optimum point
 174                 Point optimum = path.get(path.size() - 1);
 175                 FunctionEvaluationData data = evaluations.get(optimum);
 176                 if (data != null) {
 177                         if (data.getParameterValue() != null) {
 178                                 Value[] state = data.getState();
 179                                 double x = xUnit.toUnit(state[0].getValue());
 180                                 double y = yUnit.toUnit(data.getParameterValue().getValue());
 181
 182                                 XYPointerAnnotation text = new XYPointerAnnotation(trans.get("plot.label.optimum"),
 183                                                 x, y, Math.PI / 2);
 184                                 text.setTextAnchor(TextAnchor.TOP_LEFT);
 185                                 chart.getXYPlot().addAnnotation(text);
 186                         }
 187                 } else {
 188                         log.error("Could not find evaluation data for point " + optimum);
 189                 }
 190
 191
 192                 XYLineAndShapeRenderer lineRenderer = new XYLineAndShapeRenderer(true, true);
 193                 lineRenderer.setBaseShapesVisible(true);
 194                 lineRenderer.setSeriesShapesFilled(0, false);
 195                 //lineRenderer.setSeriesShape(0, shapeRenderer.getBaseShape());
 196                 lineRenderer.setSeriesOutlinePaint(0, PATH_COLOR);
 197                 lineRenderer.setSeriesPaint(0, PATH_COLOR);
 198                 lineRenderer.setUseOutlinePaint(true);
 199                 CustomXYToolTipGenerator tooltipGenerator = new CustomXYToolTipGenerator();
 200                 tooltipGenerator.addToolTipSeries(tooltips);
 201                 lineRenderer.setBaseToolTipGenerator(tooltipGenerator);
 202
 203                 XYPlot plot = chart.getXYPlot();
 204
 205                 plot.setDataset(0, new XYSeriesCollection(series));
 206                 plot.setRenderer(lineRenderer);
 207
 208
 209
 210                 return new ChartPanel(chart);
 211         }
 212
 213         /**
 214          * Create a 2D plot of the optimization path.
 215          */
 216         private ChartPanel create2DPlot(List<Point> path, Map<Point, FunctionEvaluationData> evaluations,
 217                         List<SimulationModifier> modifiers, OptimizableParameter parameter, UnitGroup stabilityUnit) {
 218
 219                 Unit parameterUnit = parameter.getUnitGroup().getDefaultUnit();
 220
 221                 SimulationModifier modX = modifiers.get(0);
 222                 SimulationModifier modY = modifiers.get(1);
 223
 224                 Unit xUnit = modX.getUnitGroup().getDefaultUnit();
 225                 Unit yUnit = modY.getUnitGroup().getDefaultUnit();
 226
 227                 // Create the optimization path dataset
 228                 XYSeries pathSeries = new XYSeries(trans.get("plot2d.path"), false, true);
 229                 List<String> pathTooltips = new ArrayList<String>();
 230                 for (Point p : path) {
 231                         FunctionEvaluationData data = evaluations.get(p);
 232                         if (data != null) {
 233                                 Value[] state = data.getState();
 234                                 pathSeries.add(xUnit.toUnit(state[0].getValue()), yUnit.toUnit(state[1].getValue()));
 235                                 pathTooltips.add(getTooltip(data, parameter));
 236                         } else {
 237                                 log.error("Could not find evaluation data for point " + p);
 238                         }
 239                 }
 240
 241
 242                 // Create evaluations dataset
 243                 double min = Double.POSITIVE_INFINITY;
 244                 double max = Double.NEGATIVE_INFINITY;
 245                 double[][] evals = new double[3][evaluations.size()];
 246                 List<String> evalTooltips = new ArrayList<String>();
 247
 248                 Iterator<FunctionEvaluationData> iterator = evaluations.values().iterator();
 249                 for (int i = 0; i < evaluations.size(); i++) {
 250                         FunctionEvaluationData data = iterator.next();
 251                         Value param = data.getParameterValue();
 252                         double value;
 253                         if (param != null) {
 254                                 value = parameterUnit.toUnit(data.getParameterValue().getValue());
 255                         } else {
 256                                 value = Double.NaN;
 257                         }
 258
 259                         Value[] state = data.getState();
 260                         evals[0][i] = xUnit.toUnit(state[0].getValue());
 261                         evals[1][i] = yUnit.toUnit(state[1].getValue());
 262                         evals[2][i] = value;
 263
 264                         if (value < min) {
 265                                 min = value;
 266                         }
 267                         if (value > max) {
 268                                 max = value;
 269                         }
 270
 271                         evalTooltips.add(getTooltip(data, parameter));
 272                 }
 273                 DefaultXYZDataset evalDataset = new DefaultXYZDataset();
 274                 evalDataset.addSeries(trans.get("plot2d.evals"), evals);
 275
 276
 277
 278                 String xLabel = modX.getRelatedObject().toString() + ": " + modX.getName() + " / " + xUnit.getUnit();
 279                 String yLabel = modY.getRelatedObject().toString() + ": " + modY.getName() + " / " + yUnit.getUnit();
 280
 281                 JFreeChart chart = ChartFactory.createXYLineChart(
 282                                 trans.get("plot2d.title"),
 283                                 xLabel,
 284                                 yLabel,
 285                                 null,
 286                                 //evalDataset,
 287                                 PlotOrientation.VERTICAL,
 288                                 true, // Legend
 289                                 true, // Tooltips
 290                                 false); // Urls
 291
 292
 293                 // Set the scale of the plot to the limits
 294                 double x1 = xUnit.toUnit(modX.getMinValue());
 295                 double x2 = xUnit.toUnit(modX.getMaxValue());
 296                 double y1 = yUnit.toUnit(modY.getMinValue());
 297                 double y2 = yUnit.toUnit(modY.getMaxValue());
 298
 299                 chart.getXYPlot().getDomainAxis().setRange(x1, x2);
 300                 chart.getXYPlot().getRangeAxis().setRange(y1, y2);
 301
 302                 XYBoxAnnotation box = new XYBoxAnnotation(x1, y1, x2, y2);
 303                 chart.getXYPlot().addAnnotation(box);
 304
 305                 int n = pathSeries.getItemCount();
 306                 XYPointerAnnotation text = new XYPointerAnnotation(trans.get("plot.label.optimum"),
 307                                 (Double) pathSeries.getX(n - 1), (Double) pathSeries.getY(n - 1), -Math.PI / 5);
 308                 text.setTextAnchor(TextAnchor.BASELINE_LEFT);
 309                 chart.getXYPlot().addAnnotation(text);
 310
 311                 PaintScale paintScale = new GradientScale(min, max);
 312
 313                 XYShapeRenderer shapeRenderer = new XYShapeRenderer();
 314                 shapeRenderer.setPaintScale(paintScale);
 315                 shapeRenderer.setUseFillPaint(true);
 316                 CustomXYToolTipGenerator tooltipGenerator = new CustomXYToolTipGenerator();
 317                 tooltipGenerator.addToolTipSeries(evalTooltips);
 318                 shapeRenderer.setBaseToolTipGenerator(tooltipGenerator);
 319
 320
 321                 shapeRenderer.getLegendItem(0, 0);
 322
 323
 324                 XYLineAndShapeRenderer lineRenderer = new XYLineAndShapeRenderer(true, true);
 325                 lineRenderer.setBaseShapesVisible(true);
 326                 lineRenderer.setSeriesShapesFilled(0, false);
 327                 lineRenderer.setSeriesShape(0, shapeRenderer.getBaseShape());
 328                 lineRenderer.setSeriesOutlinePaint(0, PATH_COLOR);
 329                 lineRenderer.setSeriesPaint(0, PATH_COLOR);
 330                 lineRenderer.setUseOutlinePaint(true);
 331                 tooltipGenerator = new CustomXYToolTipGenerator();
 332                 tooltipGenerator.addToolTipSeries(pathTooltips);
 333                 lineRenderer.setBaseToolTipGenerator(tooltipGenerator);
 334
 335
 336                 XYPlot plot = chart.getXYPlot();
 337
 338                 plot.setDataset(0, new XYSeriesCollection(pathSeries));
 339                 plot.setRenderer(lineRenderer);
 340
 341                 plot.setDataset(1, evalDataset);
 342                 plot.setRenderer(1, shapeRenderer);
 343
 344
 345                 // Add value scale
 346                 NumberAxis numberAxis = new NumberAxis(parameter.getName() + " / " + parameterUnit.getUnit());
 347                 PaintScaleLegend scale = new PaintScaleLegend(paintScale, numberAxis);
 348                 scale.setPosition(RectangleEdge.RIGHT);
 349                 scale.setMargin(4.0D, 4.0D, 40.0D, 4.0D);
 350                 scale.setAxisLocation(AxisLocation.BOTTOM_OR_RIGHT);
 351                 chart.addSubtitle(scale);
 352
 353
 354                 return new ChartPanel(chart);
 355         }
 356
 357
 358
 359         private String getTooltip(FunctionEvaluationData data, OptimizableParameter parameter) {
 360                 String ttip = "<html>";
 361                 if (data.getParameterValue() != null) {
 362                         ttip += parameter.getName() + ": " +
 363                                         parameter.getUnitGroup().getDefaultUnit().toStringUnit(data.getParameterValue().getValue());
 364                         ttip += "<br>";
 365                 }
 366                 if (data.getDomainReference() != null) {
 367                         ttip += trans.get("plot.ttip.stability") + " " + data.getDomainReference();
 368                 }
 369                 return ttip;
 370         }
 371
 372         private class GradientScale implements PaintScale {
 373
 374                 private final double min;
 375                 private final double max;
 376
 377                 public GradientScale(double min, double max) {
 378                         this.min = min;
 379                         this.max = max;
 380                 }
 381
 382                 @Override
 383                 public Paint getPaint(double value) {
 384                         if (Double.isNaN(value)) {
 385                                 return OUT_OF_DOMAIN_COLOR;
 386                         }
 387
 388                         value = MathUtil.map(value, min, max, 0.0, 1.0);
 389                         value = MathUtil.clamp(value, 0.0, 1.0);
 390
 391                         float r = (float) RED.getValue(value);
 392                         float g = (float) GREEN.getValue(value);
 393                         float b = (float) BLUE.getValue(value);
 394
 395                         return new Color(r, g, b);
 396                 }
 397
 398                 @Override
 399                 public double getLowerBound() {
 400                         return min;
 401                 }
 402
 403                 @Override
 404                 public double getUpperBound() {
 405                         return max;
 406                 }
 407         }
 408
 409 }