<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 TRANSITIONAL//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=UTF-8">
<META NAME="GENERATOR" CONTENT="GtkHTML/3.16.1">
</HEAD>
<BODY>
Hello,<BR>
<BR>
I would like to create Web Coverage Service from one netCDF file in following <A HREF="http://mapserver.gis.umn.edu/docs/howto/WCSServerFormatHowTo">http://mapserver.gis.umn.edu/docs/howto/WCSServerFormatHowTo</A><BR>
<BR>
My netCDF store rainfall matrix every five minutes in the "pluie" variable. This the CDL :<BR>
<BR>
dimensions:
<PRE>
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";
<B> double pluie(time=496, y=50, x=50);</B>
<B> :units = "mm/h";</B>
<B> :long_name = "Intensité de pluie";</B>
<B> :coordinates = "lat lon";</B>
<B> :grid_mapping = "Lambert_Conformal";</B>
<B> :_CoordinateSystems = "ProjectionCoordinateSystem LatLonCoordinateSystem";</B>
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";
}
</PRE>
<BR>
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.<BR>
When I use :<BR>
<B>gdal_translate -ot Float64 -of gTiff -b 300 -sds netCDF:"pluie-4-20020908.nc":pluie test.tif</B><BR>
A tif file is created based on the 300th band of "pluie" variable<BR>
<BR>
When I use :<BR>
<B>gdaltindex test.shp NETCDF:pluie-3-20020908.nc:pluie</B><BR>
A shape file is created but without every band, only a record with the netcdf extent and his name<BR>
<BR>
I suppose i have to create a new record for each time's band of netcdf file in the tileindex after creating with gdaltindex ?!<BR>
<BR>
I would like to use the Python Script of Norman Barker <A HREF="http://www.nabble.com/Map-File-Problem--td5541099.html#a5544967">http://www.nabble.com/Map-File-Problem--td5541099.html#a5544967</A><BR>
<BR>
But it is possible or the only solution is to create a tif file for each band ?<BR>
<BR>
Thanks<BR>
<BR>
Vincent<BR>
<BR>
<BR>
<BR>
</BODY>
</HTML>