[GRASS-user] issues with netCDF import

[Anna Petrášová]

Hi Damian,


On Wed, Mar 19, 2014 at 1:47 PM, Damian Maddalena <damianmaddalena at gmail.com
> wrote:

> 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 does not respect region


> 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.
>

I wanted to try it myself, but it's too big to download it now. Could you
post output of:

gdalinfo NETCDF:"global_30s_2009_153.nc":lai

I am not an expert on netcdf but others might have some idea based on the
gdalinfo output.

Regards,

Anna


> Thanks!!
>
> -Damian
>
>
>
>
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20140320/afe78073/attachment.html>