[GRASS-user] Problem importing netCDF

Hanlie Pretorius hanlie.pretorius at gmail.com
Fri Jun 4 06:13:57 EDT 2010


Hi,

I'm using Win XP and GRASS 6.4RC6 to import a TRMM netCDF file using
the command:
-----
r.in.gdal -o input =c:/data/3B42.000202.0.6.nc output=test band=precipitation
-----

and am getting the error:
-----
G_set_window(): Illegal latitude for North
-----

gdalinfo for the file outputs:
-----
Driver: netCDF/Network Common Data Format
Files: C:\data\3B42.000202.0.6.nc
       C:\data\3B42.000202.0.6.nc.aux.xml
Size is 512, 512
Coordinate System is `'
Metadata:
  NC_GLOBAL#Conventions=CF-1.0
Subdatasets:
  SUBDATASET_1_NAME=NETCDF:"C:\data\3B42.000202.0.6.nc":precipitation
  SUBDATASET_1_DESC=[400x1440] precipitation (32-bit floating-point)
  SUBDATASET_2_NAME=NETCDF:"C:\data\3B42.000202.0.6.nc":relativeError
  SUBDATASET_2_DESC=[400x1440] relativeError (32-bit floating-point)
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  512.0)
Upper Right (  512.0,    0.0)
Lower Right (  512.0,  512.0)
Center      (  256.0,  256.0)
-----

In Panoply (http://www.giss.nasa.gov/tools/panoply/), the information
for the file displays as:
-----
netcdf file:/C:/data/3B42.000202.0.6.nc {
 dimensions:
   lat = 400;
   lon = 1440;
 variables:
   float precipitation(lat=400, lon=1440);
     :long_name = "precipitation";
     :units = "mm hr-1";
     :_FillValue = -9999.9f; // float
   float relativeError(lat=400, lon=1440);
     :long_name = "relativeError";
     :units = "mm hr-1";
     :_FillValue = -9999.9f; // float
   float lat(lat=400);
     :standard_name = "latitude";
     :units = "degrees_north";
     :long_name = "latitude";
     :_CoordinateAxisType = "Lat";
   float lon(lon=1440);
     :standard_name = "longitude";
     :units = "degrees_east";
     :long_name = "longitude";
     :_CoordinateAxisType = "Lon";

 :Conventions = "CF-1.0";
}
-----

In the Panoply array display, the coordinates range from 49.875° to
-49.875° (lattitude) and  -179.875° to +179.875° longitude. These are
coordinates for the centres of 0.25° grid cells.

I have set my region according to the documentation for the dataset to:
-----
projection: 3 (Latitude-Longitude)
zone:       0
datum:      wgs84
ellipsoid:  wgs84
north:      50N
south:      50S
west:       180W
east:       180E
nsres:      0:15
ewres:      0:15
rows:       400
cols:       1440
cells:      576000
-----

Can someone perhaps spot my error?

Thanks
Hanlie


More information about the grass-user mailing list