[mapserver-users] Netcdf / WCS

[vincent.thierion at ema.fr]

Hello,

I'm not expert on OWS, but serve netCDF files within WCS seems not  
totally achieve. The
case of Barker concerns a set of netCDF files.
The Galeon project is trying to develop this issue :
http://galeon-wcs.jot.com/Implementation%20Progress%20Page.

For my case, i think the best way is to create one netCDF file per  
time step or to create
a script to generate TIF files and after serves them.

Best Regards

Vincent

Steve Lime <Steve.Lime at dnr.state.mn.us> a écrit :

> I need to think about this. For some reason I remember one of the test
> cases being
> just like this although with not so many bands. As the cobwebs clear I
> believe this
> was started but not completed. There's even a BANDSITEM member of the
> layer object
> defined in the source. This would allow definition of a column to map a
> temporal value
> to a particular band within a larger image. I'm thinking it's probably
> not far from functioning.
>
> Steve
>
>>>> On 2/19/2008 at 10:08 AM, in message <47BAFEE1.60508 at pobox.com>,
> Frank
> Warmerdam <warmerdam at pobox.com> wrote:
>> vincent wrote:
>>> Hello,
>>>
>>> I would like to create Web Coverage Service from one netCDF file in
>
>>> following
> http://mapserver.gis.umn.edu/docs/howto/WCSServerFormatHowTo
>>>
>>> My netCDF store rainfall matrix every five minutes in the "pluie"
>>> variable. This the CDL :
>>>
>>> dimensions:
>>>
>>>    time = UNLIMITED;   // (496 currently)   // (has coord.var)
>>>    x = 50;   // (has coord.var)
>>>    y = 50;   // (has coord.var)
>>> variables:
>>>    int Lambert_Conformal;
>>>      :grid_mapping_name = "lambert_conformal_conic";
>>>      :standard_parallel = 43.20317, 44.99683; // double
>>>      :longitude_of_central_meridian = 2.337229; // double
>>>      :latitude_of_projection_origin = 44.1; // double
>>>      :false_easting = 600.0; // double
>>>      :false_northing = 3200.0; // double
>>>      :_CoordinateTransformType = "Projection";
>>>      :_CoordinateSystems = "ProjectionCoordinateSystem";
>>>      :_CoordinateAxes = "y x";
>>>      :_CoordinateAxisTypes = "GeoX GeoY";
>>> *   double pluie(time=496, y=50, x=50);*
>>> *     :units = "mm/h";*
>>> *     :long_name = "Intensité de pluie";*
>>> *     :coordinates = "lat lon";*
>>> *     :grid_mapping = "Lambert_Conformal";*
>>> *     :_CoordinateSystems = "ProjectionCoordinateSystem
>> LatLonCoordinateSystem";*
>>
>> Vincent,
>>
>> Because GDAL treats this array as a single dataset having 496 bands,
> I don't
>> believe there is a way to properly create a tileindex with one entry
> per
>> time slice, referring back directly to that band though this is
> something
>> that
>> I think was discussed at one point.
>>
>> You *can* serve the pluie array, but it will be represented as having
> 496
>> bands along a dimension called Band instead of time.
>>
>> Hopefully Steve will correct me if I am wrong.  What would be nice
>> is a way to have WCS (and raster rendering in general) potentially
>> override the BANDS= processing option with information read from a
>> tile index, and then update the gdaltindex program accordingly to
>> generate such results.  This shouldn't - in theory - be too hard to
>> implement.
>>
>> Best regards,
>> --
>>
> ---------------------------------------+------------------------------------
>> --
>> I set the clouds in motion - turn up   | Frank Warmerdam,
>> warmerdam at pobox.com
>> light and sound - activate the windows | http://pobox.com/~warmerdam
>
>> and watch the world go round - Rush    | President OSGeo,
> http://osgeo.org
>>
>> _______________________________________________
>> mapserver-users mailing list
>> mapserver-users at lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/mapserver-users
>
>



----------------------------------------------------
Ce message a ete envoye par le serveur IMP de l'EMA.