[GRASS-user] issues with netCDF import

Damian Maddalena damianmaddalena at gmail.com
Wed Mar 19 10:47:32 PDT 2014


I am having some difficulty with the import of a netCDF dataset into 
GRASS 7. The dataset I am attempting to import is a 30 second global LAI 
dataset that can be found here: 
http://globalchange.bnu.edu.cn/research/lai/

I have successfully imported other netCDF files into this build 
(worldCLIM and the soils data from the same source as this LAI data, for 
example), but this particular data is giving me issues. I have tried 
several files from several of the years available of the global LAI 
product. I get the same results.

I am running the following commands in a global lat/lon location:

g.region n=90 s=-90 e=-180 w=180 res=0.008333333767951  (I have also 
tried res=00:00:30.)
r.in.gdal input=NETCDF:"global_30s_2009_153.nc":lai output=testNETCDF -o

The file imports with some warnings.

Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
WARNING: Over-riding projection check
WARNING: Fixing subtle input data rounding error of south boundary
          (-9.38773e-06>4.62963e-08)

It appears to be null when plotted, though r.info reports min of 0 and 
max of 70.  r.info hints at an issue with the e-w resolution:

  |            N:        90N    S:        90S   Res: 0:00:30
  |            E: 179:59:59.932408W    W:       180W   Res: 0:00:00.000002

r.stats seems to indicate that all data is null?

r.stats testNETCDF -cp
100%
* 900072000 100.00%

d.histogram plots nothing.

Any hints would be greatly appreciated.

Thanks!!

-Damian






More information about the grass-user mailing list