[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