Changeset 191 for ether_megapoli/trunk/service/implementation/com/ether/EtherPlotContentServiceImpl.java

Timestamp:

09/12/11 18:11:46 (13 years ago)

Author:

vmipsl

Message:

Servlet _ Contour en cours

File:

• ether_megapoli/trunk/service/implementation/com/ether/EtherPlotContentServiceImpl.java (modified) (7 diffs)

ether_megapoli/trunk/service/implementation/com/ether/EtherPlotContentServiceImpl.java

@@ -5 +5 @@
 import gov.noaa.pmel.sgt.Axis;
 import gov.noaa.pmel.sgt.CartesianGraph;
+import gov.noaa.pmel.sgt.ColorMap;
+import gov.noaa.pmel.sgt.ContourLevels;
 import gov.noaa.pmel.sgt.Graph;
+import gov.noaa.pmel.sgt.GridAttribute;
+import gov.noaa.pmel.sgt.IndexedColorMap;
 import gov.noaa.pmel.sgt.JPane;
 import gov.noaa.pmel.sgt.Layer;
@@ -16 +20 @@
 import gov.noaa.pmel.sgt.StackedLayout;
 import gov.noaa.pmel.sgt.TimeAxis;
+import gov.noaa.pmel.sgt.demo.TestData;
+import gov.noaa.pmel.sgt.dm.SGTData;
 import gov.noaa.pmel.sgt.dm.SGTMetaData;
+import gov.noaa.pmel.sgt.dm.SimpleGrid;
 import gov.noaa.pmel.sgt.dm.SimpleLine;
+import gov.noaa.pmel.sgt.swing.JPlotLayout;
 import gov.noaa.pmel.util.Dimension2D;
 import gov.noaa.pmel.util.GeoDateArray;
@@ -30 +38 @@
 import java.awt.*;
 import java.util.ArrayList;
+import java.util.Collections;
 import java.util.Date;
 import java.util.List;
@@ -79 +88 @@
         final List<Pair> pIdPIdListForSecondAxis = listPairOfPIdPIdForTwoAxes.getSecondValue();
 
-        final List<SimpleLine> allLines = createLines( pIdPIdListForFirstAxis, megapoliPlot.getBeginDate(), megapoliPlot.getEndDate(), null, locale );
-        final List<SimpleLine> allLinesForSecondAxis = createLines( pIdPIdListForSecondAxis, megapoliPlot.getBeginDate(), megapoliPlot.getEndDate(), 10, locale );
+        final List<SimpleLine> allLines = createLines( pIdPIdListForFirstAxis, megapoliPlot.getBeginDate(), megapoliPlot.getEndDate(), 10, locale );
+        final List<SimpleLine> allLinesForSecondAxis = createLines( pIdPIdListForSecondAxis, megapoliPlot.getBeginDate(), megapoliPlot.getEndDate(), 50, locale );
         List<SimpleLine> lines = extractRealLines( allLines );
         List<SimpleLine> linesForSecondAxis = extractRealLines( allLinesForSecondAxis );
@@ -171 +180 @@
     }
 
+    @NotNull
+    public JPane create2DPlot( @NotNull final MegapoliPlot megapoliPlot, @NotNull final Integer plotWidth, @NotNull final Integer plotHeight, @Nullable final Locale locale )
+    {
+//        final JPane mainPane = new JPane( "2D Pane", new Dimension( plotWidth, plotHeight ) );
+
+        final ResourceBundle bundle = WebHelper.getBundle( locale );
+
+        final List<Double> parameterValues = (List<Double>) megapoliPlot.getData().getFirstArray();
+        final List<Double> latitudeValues = (List<Double>) megapoliPlot.getData().getSecondArray();
+        final List<Double> longitudeValues = (List<Double>) megapoliPlot.getData().getThirdArray();
+
+        SGTData newData;
+        TestData td;
+
+        /*
+       * Create a test grid with sinasoidal-ramp data.
+        */
+//        Range2D xr = new Range2D( 190.0f, 250.0f, 1.0f );
+//        Range2D yr = new Range2D( 0.0f, 45.0f, 1.0f );
+        Range2D xr = new Range2D( 190.0f, 200.0f, 1.0f );
+        Range2D yr = new Range2D( 0.0f, 4.0f, 1.0f );
+//        newData = testDataa( xr, yr, TestData.SINE_RAMP, 12.0f, 30.f, 5.0f, parameterArray, latitudeArray, longitudeValues );
+
+
+//            public SGTData testDataa( Range2D range1, Range2D range2, int type, float amp, float off, float per, final double[] parameterArray, final double[] latitudeArray, final double[] longitudeValues )
+        final SimpleGrid simpleGrid;
+        final SGTMetaData xMeta;
+        final SGTMetaData yMeta;
+        final SGTMetaData zMeta;
+
+        final int num1 = (int) ( ( xr.end - xr.start ) / xr.delta ) + 1;
+        final double[] axis1 = new double[num1];
+        final int num2 = (int) ( ( yr.end - yr.start ) / yr.delta ) + 1;
+        final double[] axis2 = new double[num2];
+
+        for( int count = 0; count < num1; count++ )
+            axis1[count] = xr.start + count * xr.delta;
+
+        for( int count = 0; count < num2; count++ )
+            axis2[count] = yr.start + count * yr.delta;
+
+//        final double[] valuesOrigPetites = getValues( axis1, num1, axis2, 2, TestData.SINE_RAMP, 12.0f, 30.f, 5.0f );
+//        final double[] valuesOrigPetites = getValues( longitudeValues, longitudeValues.length, latitudeValues, latitudeValues.length, TestData.SINE_RAMP, 12.0f, 30.f, 5.0f );
+
+//        final double[] valuesOrig = new double[num1 * num2];
+//        for( int count = 0; count < valuesOrigPetites.length; count++ )
+//        {
+//            valuesOrig[count] = valuesOrigPetites[count];
+//        }
+//        for( int count = valuesOrigPetites.length; count < num1 * num2; count++ )
+//        {
+//            valuesOrig[count] = Double.NaN;
+//        }
+
+        final double[] parameterValuesArray = EtherHelper.convertListDoubleToDoubleArray( parameterValues );
+        final double[] parameterValuesMap = new double[latitudeValues.size() * longitudeValues.size()];
+        for( int i = 0; i < parameterValuesArray.length; i++ )
+            parameterValuesMap[i] = parameterValuesArray[i];
+
+        for( int i = parameterValuesArray.length; i < latitudeValues.size() * longitudeValues.size(); i++ )
+            parameterValuesMap[i] = Double.NaN;
+
+        SGLabel keyLabel = new SGLabel( "Key Label", "", new Point2D.Double( 0.0, 0.0 ) );
+        keyLabel.setHeightP( 0.16 );
+        //
+        // create SimpleGrid
+        //
+        zMeta = new SGTMetaData( "ts", "m s-1" );
+        keyLabel.setText( "XY test grid" );
+        xMeta = new SGTMetaData( bundle.getString( "plot.longitude" ), bundle.getString( "plot.degres" ) );
+        yMeta = new SGTMetaData( bundle.getString( "plot.latitude" ), bundle.getString( "plot.degres" ) );
+//        final SimpleGrid sg = new SimpleGrid( valuesOrig, axis1, axis2, "Test Series" );
+//        simpleGrid = new SimpleGrid( parameterValues, longitudeValues, latitudeValues, "Test Series" );
+//        simpleGrid = new SimpleGrid( valuesOrig, longitudeValues, latitudeValues, "Test Series" );
+//        simpleGrid = new SimpleGrid( valuesOrig, axis1, axis2, "Test Series" );
+
+        final double[] longitudeValuesArray = EtherHelper.convertListDoubleToDoubleArray( longitudeValues );
+        final double[] latitudeValuesArray = EtherHelper.convertListDoubleToDoubleArray( latitudeValues );
+        simpleGrid = new SimpleGrid( parameterValuesMap, longitudeValuesArray, latitudeValuesArray, "Test Series" );
+
+        simpleGrid.setXMetaData( xMeta );
+        simpleGrid.setYMetaData( yMeta );
+        simpleGrid.setZMetaData( zMeta );
+        simpleGrid.setKeyTitle( keyLabel );
+        simpleGrid.setRealFullGrid( false );
+
+
+        // Create the contour levels and color map
+        final Double min = Collections.min( parameterValues );
+        final Double max = Collections.max( parameterValues );
+        final Double delta = ( max - min ) / 64;
+        final Range2D datar = new Range2D( min, max, delta );
+//        final Range2D datar = new Range2D( Collections.min( Arrays.asList( parameterValuesMap ) ), 35f, 0.5f );
+        final ContourLevels clevels = ContourLevels.getDefault( datar );
+        final ColorMap cmap = createColorMap( datar );
+
+//        final Range2D datarOrig = new Range2D( 5f, 10.0f, 1f );
+//        final ContourLevels clevelsOrig = ContourLevels.getDefault( datarOrig );
+//        final ColorMap cmapOrig = createColorMap( datarOrig );
+//        final GridAttribute gridAttributeOrig = new GridAttribute( clevelsOrig );
+//        gridAttributeOrig.setColorMap( cmapOrig );
+//        gridAttributeOrig.setStyle( GridAttribute.RASTER_CONTOUR );
+
+
+        // Create a grid for contour data
+        final GridAttribute gridAttribute = new GridAttribute( clevels );
+        gridAttribute.setColorMap( cmap );
+        gridAttribute.setStyle( GridAttribute.RASTER );
+
+        // Create the plot
+        final JPlotLayout jPlotLayout = new JPlotLayout( true, false, false, "JPlotLayout Pane", null, true );
+        jPlotLayout.setEditClasses( false );
+        jPlotLayout.addData( simpleGrid, gridAttribute, "First Data" );
+        jPlotLayout.setTitles( "", "", "" );
+        jPlotLayout.setSize( new Dimension( plotWidth, plotHeight ) );
+
+        /*
+       * Resize the key Pane, both the device size and the physical
+       * size. Set the size of the key in physical units and place
+       * the key pane at the "South" of the frame.
+        */
+//    rpl.setKeyLayerSizeP(new Dimension2D(6.0, 1.02));
+//    rpl.setKeyBoundsP(new Rectangle2D.Double(0.01, 1.01, 5.98, 1.0));
+
+        return jPlotLayout;
+//        return mainPane;
+    }
+
+    @NotNull
+    protected String formatTitle( @Nullable final String title )
+    {
+        if( null == title )
+            return "";
+        else
+        {
+            final String formatTitle = title.replaceAll( "\\n", "<br>" );
+            return "<html>" + formatTitle + "</html>";
+        }
+    }
+
     private List<SimpleLine> createEmptyLinesList( final List<SimpleLine> allLines, final List<SimpleLine> allLinesForSecondAxis, final List<SimpleLine> lines )
     {
@@ -211 +360 @@
 
         return lines;
-    }
-
-    @NotNull
-    public JPane create2DPlot( @NotNull final MegapoliPlot megapoliPlot, @NotNull final Integer plotWidth, @NotNull final Integer plotHeight, @Nullable final Locale locale )
-    {
-        final JPane mainPane = new JPane( "2D Pane", new Dimension( plotWidth, plotHeight ) );
-
-        final ResourceBundle bundle = WebHelper.getBundle( locale );
-
-        final Data realValuesFor2D = BouchonHelper.createRealValuesFor2D();
-//        final double[] parameterValues = (double[]) megapoliPlot.getData().getFirstArray();
-//        final double[] latitudeValues = (double[]) megapoliPlot.getData().getSecondArray();
-//        final double[] longitudeValues = (double[]) megapoliPlot.getData().getThirdArray();
-
-        final double[] parameterValues = (double[]) realValuesFor2D.getFirstArray();
-        final double[] latitudeValues = (double[]) realValuesFor2D.getSecondArray();
-        final double[] longitudeValues = (double[]) realValuesFor2D.getThirdArray();
-
-        return mainPane;
-    }
-
-    @NotNull
-    protected String formatTitle( @Nullable final String title )
-    {
-        if( null == title )
-            return "";
-        else
-        {
-            final String formatTitle = title.replaceAll( "\\n", "<br>" );
-            return "<html>" + formatTitle + "</html>";
-        }
     }
 
@@ -459 +577 @@
     }
 
+    private double[] getValues( double[] axis1, int num1, double[] axis2, int num2,
+                                int type, float amp, float off, float per )
+    {
+        double[] values;
+        int count1, count2, count;
+        values = new double[num1 * num2];
+        switch( type )
+        {
+            default:
+            case TestData.RANDOM:
+                for( count = 0; count < num1 * num2; count++ )
+                {
+                    values[count] = amp * Math.random() + off;
+                }
+                break;
+            case TestData.SINE:
+                count = 0;
+                for( count1 = 0; count1 < num1; count1++ )
+                {
+                    for( count2 = 0; count2 < num2; count2++ )
+                    {
+                        values[count] = amp * Math.sin( axis1[count1] / per ) *
+                                Math.sin( axis2[count2] / per ) + off;
+                        count++;
+                    }
+                }
+                break;
+            case TestData.SINE_RAMP:
+                count = 0;
+                double ax1factr = 0.08 * Math.abs( axis1[0] - axis1[num1 - 1] ) /
+                        Math.max( Math.abs( axis1[0] ), Math.abs( axis1[num1 - 1] ) );
+                double ax2factr = 0.08 * Math.abs( axis2[0] - axis2[num2 - 1] ) /
+                        Math.max( Math.abs( axis2[0] ), Math.abs( axis2[num2 - 1] ) );
+                for( count1 = 0; count1 < num1; count1++ )
+                {
+                    for( count2 = 0; count2 < num2; count2++ )
+                    {
+                        values[count] = amp * Math.sin( axis1[count1] / per ) *
+                                Math.sin( axis2[count2] / per ) + off +
+                                amp * ( ax1factr * count1 - ax2factr * count2 );
+                        count++;
+                    }
+                }
+        }
+        return values;
+    }
+
+
+    @NotNull
+    private ColorMap createColorMap( final Range2D datar )
+    {
+        final int[] red =
+                {0, 0, 0, 0, 0, 0, 0, 0,
+                        0, 0, 0, 0, 0, 0, 0, 0,
+                        0, 0, 0, 0, 0, 0, 0, 0,
+                        0, 7, 23, 39, 55, 71, 87, 103,
+                        119, 135, 151, 167, 183, 199, 215, 231,
+                        247, 255, 255, 255, 255, 255, 255, 255,
+                        255, 255, 255, 255, 255, 255, 255, 255,
+                        255, 246, 228, 211, 193, 175, 158, 140};
+        final int[] green =
+                {0, 0, 0, 0, 0, 0, 0, 0,
+                        0, 11, 27, 43, 59, 75, 91, 107,
+                        123, 139, 155, 171, 187, 203, 219, 235,
+                        251, 255, 255, 255, 255, 255, 255, 255,
+                        255, 255, 255, 255, 255, 255, 255, 255,
+                        255, 247, 231, 215, 199, 183, 167, 151,
+                        135, 119, 103, 87, 71, 55, 39, 23,
+                        7, 0, 0, 0, 0, 0, 0, 0};
+        final int[] blue =
+                {0, 143, 159, 175, 191, 207, 223, 239,
+                        255, 255, 255, 255, 255, 255, 255, 255,
+                        255, 255, 255, 255, 255, 255, 255, 255,
+                        255, 247, 231, 215, 199, 183, 167, 151,
+                        135, 119, 103, 87, 71, 55, 39, 23,
+                        7, 0, 0, 0, 0, 0, 0, 0,
+                        0, 0, 0, 0, 0, 0, 0, 0,
+                        0, 0, 0, 0, 0, 0, 0, 0};
+
+        final IndexedColorMap cmap = new IndexedColorMap( red, green, blue );
+        cmap.setTransform( new LinearTransform( 0.0, (double) red.length, datar.start, datar.end ) );
+        return cmap;
+    }
 
     public EtherService getEtherService()