[mapserver-users] Netcdf / WCS

[vincent]

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";
   char LatLonCoordinateSystem;
     :_CoordinateAxes = "time lat lon";
   char ProjectionCoordinateSystem;
     :_CoordinateAxes = "time y x";
     :_CoordinateTransforms = "LambertConformalProjection";
   double time(time=496);
     :long_name = "rainfall time";
     :standard_name = "time";
     :units = "seconds since 1970-01-01 01:00:00";
     :axis = "T";
     :_CoordinateAxisType = "Time";
     :bounds = "time_bnds";
   double x(x=50);
     :units = "km";
     :long_name = "x coordinate of projection";
     :standard_name = "projection_x_coordinate";
     :_CoordinateAxisType = "GeoX";
   double y(y=50);
     :units = "km";
     :long_name = "y coordinate of projection";
     :standard_name = "projection_y_coordinate";
     :_CoordinateAxisType = "GeoY";
   double lat(x=50, y=50);
     :units = "degrees_north";
     :long_name = "latitude coordinate";
     :standard_name = "latitude";
     :_CoordinateAxisType = "Lat";
   double lon(x=50, y=50);
     :units = "degrees_east";
     :long_name = "longitude";
     :standard_name = "longitude coordinate";
     :_CoordinateAxisType = "Lon";

 :Zone_Radar = "Gard-Ceze Amont";
 :history = "Tuesday 19 February 2008, 00:00";
 :Projection_cartographique = "Lambert III";
 :X0 = 700.7; // double
 :Y0 = 3189.1; // double
 :Conventions = "CF-1.0";
}


The HowTo of Mapserver Website seems adapt to a set of netcdf files but
not for one netCDF file including several bands corresponding to each
time step.
When I use :
gdal_translate -ot Float64 -of gTiff -b 300 -sds
netCDF:"pluie-4-20020908.nc":pluie test.tif
A tif file is created based on the 300th band of "pluie" variable

When I use :
gdaltindex test.shp NETCDF:pluie-3-20020908.nc:pluie
A shape file is created but without every band, only a record with the
netcdf extent and his name

I suppose i have to create a new record for each time's band of netcdf
file in the tileindex after creating with gdaltindex ?!

I would like to use the Python Script of Norman Barker
http://www.nabble.com/Map-File-Problem--td5541099.html#a5544967

But it is possible or the only solution is to create a tif file for each
band ?

Thanks

Vincent



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20080219/192e5db1/attachment.htm>