John Caron wrote:
I have remade this test dataset
ftp://ftp.unidata.ucar.edu/pub/galeon/striped.nc
in order to make the latitude coordinates evenly spaced ("rectified").
If you have already downloaded, please get it again so we all have the
same test datasets.
All 3 test datasets have evenly spaced coordinates in X,Y or Lat, Lon.
An interesting test later might be to deal with non-rectified X,Y.
I am realizing that striped.nc also has irregularly spaced vertical
coordinates. I was thinking that didnt matter, but now im realizing that
i cant use a RectifiedGrid like the one below. Can anyone give me an
example of how it should be changed to accomodate an irregular vertical
coordinate?
<CoverageDescription xmlns="http://www.opengis.net/wcs"; xmlns:xlink="http://www.w3.org/1999/xlink";
xmlns:gml="http://www.opengis.net/gml"; xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"; version="1.0.0"
xsi:schemaLocation="http://www.opengis.net/wcs http://schemas.opengis.net/wcs/1.0.0/describeCoverage.xsd";>
<CoverageOffering>
<name>u_wind</name>
<label>u wind @ isobaric</label>
<lonLatEnvelope srsName="WGS84(DD)">
<gml:pos>-133.70633026238878 11.747958797813363</gml:pos>
<gml:pos>-48.67045857136577 57.47373334674363</gml:pos>
</lonLatEnvelope>
<domainSet>
<spatialDomain>
<gml:Envelope srsName="WGS84(DD)">
<gml:pos>-133.70633026238878 11.747958797813363</gml:pos>
<gml:pos>-48.67045857136577 57.47373334674363</gml:pos>
</gml:Envelope>
<gml:RectifiedGrid dimension="3">
<gml:limits>
<gml:GridEnvelope>
<gml:low>0 0 0 </gml:low>
<gml:high>92 64 18 </gml:high>
</gml:GridEnvelope>
</gml:limits>
<gml:axisName>x</gml:axisName>
<gml:axisName>y</gml:axisName>
<gml:axisName>z</gml:axisName>
<gml:origin>
<gml:pos>-4226.1069969154705 -832.6982610175619 100.0 </gml:pos>
</gml:origin>
<gml:offsetVector>81.27099999999973 0.0 0.0 </gml:offsetVector>
<gml:offsetVector>0.0 81.27099999999996 0.0 </gml:offsetVector>
<gml:offsetVector>0.0 0.0 50.0 </gml:offsetVector>
</gml:RectifiedGrid>
</spatialDomain>
<temporalDomain>
<gml:timePosition>2002-12-02T22:00:00Z</gml:timePosition>
<gml:timePosition>2002-12-03T01:00:00Z</gml:timePosition>
</temporalDomain>
</domainSet>
<rangeSet>
<RangeSet>
<name>RangeSetName</name>
<label>RangeSetLabel</label>
<axisDescription>
<AxisDescription>
<name>Vertical</name>
<label>isobaric</label>
<values>
<singleValue> 100.0</singleValue>
<singleValue> 150.0</singleValue>
<singleValue> 200.0</singleValue>
<singleValue> 250.0</singleValue>
<singleValue> 300.0</singleValue>
<singleValue> 350.0</singleValue>
<singleValue> 400.0</singleValue>
<singleValue> 450.0</singleValue>
<singleValue> 500.0</singleValue>
<singleValue> 550.0</singleValue>
<singleValue> 600.0</singleValue>
<singleValue> 650.0</singleValue>
<singleValue> 700.0</singleValue>
<singleValue> 750.0</singleValue>
<singleValue> 800.0</singleValue>
<singleValue> 850.0</singleValue>
<singleValue> 900.0</singleValue>
<singleValue> 950.0</singleValue>
<singleValue> 1000.0</singleValue>
</values>
</AxisDescription>
</axisDescription>
<nullValues>
<singleValue>NaN</singleValue>
</nullValues>
</RangeSet>
</rangeSet>
<supportedCRSs>
<requestCRSs>EPSG:4326</requestCRSs>
<responseCRSs>EPSG:4326</responseCRSs>
</supportedCRSs>
<supportedFormats>
<formats>GeoTIFF</formats>
<formats>GeoTIFFfloat</formats>
<formats>NetCDF3</formats>
</supportedFormats>
<supportedInterpolations>
<interpolationMethod>none</interpolationMethod>
</supportedInterpolations>
</CoverageOffering>
</CoverageDescription>