[GRASS-user] error while displaying netCDF file

[Veronica Andreo]

Sorry for the noise :)

I solved it.
Just forgot "-l" in r.in.gdal and then set the bounds of the map properly
with r.region.

Cheers,
Vero

2015-12-16 18:01 GMT-03:00 Veronica Andreo <veroandreo at gmail.com>:

> Hi list,
>
> I'm trying to import a netCDF file I downloaded from here:
> http://oceandata.sci.gsfc.nasa.gov/MODISA/Mapped/8Day/4km/chl/
>
> I've tried two different ways of importing:
>
> r.in.gdal input=NETCDF:"A20030012003008.L3m_8D_CHL_chl_ocx_4km.nc":chl_ocx
> output=test_chlor -o
>
> and
>
> gdal_translate -a_srs EPSG:4326 -a_ullr -180 90 180 -90 NETCDF:"
> A20030012003008.L3m_8D_CHL_chl_ocx_4km.nc":chl_ocx test_chlor_b.tif
> r.in.gdal input=test_chlor_b.tif output=test_chlor_b
>
> In both cases I get this error when trying to display the map (both with
> d.mon - d.rast and in the GUI):
>
> *GRASS_INFO_ERROR(18597,1): Syntax error in cell header*
> *GRASS_INFO_END(18597,1)*
>
> though, r.info test_chlor_b, seems fine
>
> [...]
>
>  |   Type of Map:  raster               Number of Categories: 0
>
>  |   Data Type:    FCELL
>
>  |   Rows:         4320
>
>  |   Columns:      8640
>
>  |   Total Cells:  37324800
>
>  |        Projection: Latitude-Longitude
>
>  |            N:        90N    S:        90S   Res: 0:02:30
>
>  |            E:       180E    W:       180W   Res: 0:02:30
>
>  |   Range of data:    min = 0.001  max = 99.96646
>
> [...]
>
> While transforming with gdal_translate and also when unsing r.in.gdal with
> the original netCDF file, I noted this:
>
> gdal_translate -a_srs EPSG:4326 -a_ullr -180 90 180 -90 NETCDF:"
> A20030012003008.L3m_8D_CHL_chl_ocx_4km.nc":chl_ocx test_chlor_b.tif
> *Warning 1: dimension #1 (lon) is not a Longitude/X dimension.*
> *Warning 1: dimension #0 (lat) is not a Latitude/Y dimension.*
> Input file size is 8640, 4320
> 0...10...20...30...40...50...60...70...80...90...100 - done.
>
> I guess it is maybe related to the error? What does it mean? Can someone
> help me, please? QGIS displays it fine...
>
> I'm using grass71 r67160, under fedora 22. GDAL 1.11.3
>
> Thanks a lot! And sorry for my ignorance :/
> Vero
>
> ps: Here, the output of gdalinfo from original netCDF and the one
> converted with gdal_translate: http://pastebin.com/4c7Hk7cW
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20151216/75e04eda/attachment.html>