[gdal-dev] gdal_translate ERROR 1: nc_put_vara_double(IMAGE.grd): NetCDF: Start+count exceeds dimension bound

[Joaquim Luis]

Well,

It's a relatively big file (46800 x 55200) and since GMT internally uses  
floats, that converts to an array of ~9.8 Gb so you need to have available  
a memory chunk of that size to process it. Aside from grdcuting it and  
process it on pieces I don't see what else you can do.

> Hello Joaquim,
>
>
> I thought it is clever to use -of GMT because gmt grdmath is not working
> with the GeoTIFF.
>
>      grdmath: gdalread: failure to allocate enough memory
>      grdmath: ERROR reading file with gdalread.
>      grdmath (GMTAPI_Import_Grid): Could not create file [IMAGE.tif]
>      [Session gmt (0)]: Error returned from GMT API: GMT_GRID_READ_ERROR  
> (18)
>
> Now I tried:
>
>      gdal_translate -of netCDF IMAGE.tif IMAGE.nc
>
> but
>
>      Input file size is 46800, 55200
>      0...10...20...30...40...50Warning 1: Unsupported GDAL datatype (2),  
> treat as NC_FLOAT.
>      ...60...70...80...90...100 - done.
>
> and gmt grdmath can't handle the output file.
>
> Any suggestions?
>
>
> Thanks for your help!
>
> Sebastian.
>
>
>
>
> Am 12/09/14 15:08, schrieb Joaquim Luis:
>> Sebastian,
>>
>> Any particular reason why you would want to use that so old and
>> abandoned format (-of GMT)? Do not confuse it with the netCDF CF that
>> GMT uses since GMT 4.0 (released many years ago). Note also that in
>> GMT5 you can use the  GeoTiff directly, or convert it with grdreformat
>> or with gdal_translate but to netCDF (not to 'GMT')
>>
>> Joaquim
>>
>>> Hello,
>>>
>>> I have a problem with the following GeoTIFF:
>>>
>>> gdalinfo IMAGE.tif
>>> Driver: GTiff/GeoTIFF
>>> Files: IMAGEDATA/IMAGE.tif
>>> Size is 46800, 55200
>>> Coordinate System is:
>>> PROJCS["WGS 84 / UTM zone 30S",
>>>      GEOGCS["WGS 84",
>>>          DATUM["WGS_1984",
>>>              SPHEROID["WGS 84",6378137,298.257223563,
>>>                  AUTHORITY["EPSG","7030"]],
>>>              AUTHORITY["EPSG","6326"]],
>>>          PRIMEM["Greenwich",0],
>>>          UNIT["degree",0.0174532925199433],
>>>          AUTHORITY["EPSG","4326"]],
>>>      PROJECTION["Transverse_Mercator"],
>>>      PARAMETER["latitude_of_origin",0],
>>>      PARAMETER["central_meridian",-3],
>>>      PARAMETER["scale_factor",0.9996],
>>>      PARAMETER["false_easting",500000],
>>>      PARAMETER["false_northing",10000000],
>>>      UNIT["metre",1,
>>>          AUTHORITY["EPSG","9001"]],
>>>      AUTHORITY["EPSG","32730"]]
>>> Origin = (422999.375000000000000,2201000.625000000000000)
>>> Pixel Size = (1.250000000000000,-1.250000000000000)
>>> Metadata:
>>>    AREA_OR_POINT=Point
>>> Image Structure Metadata:
>>>    COMPRESSION=PACKBITS
>>>    INTERLEAVE=BAND
>>> Corner Coordinates:
>>> Upper Left  (  422999.375, 2201000.625) (  5d 2'44.07"W, 70d17' 7.59"S)
>>> Lower Left  (  422999.375, 2132000.625) (  5d 6'33.10"W, 70d54'13.43"S)
>>> Upper Right (  481499.375, 2201000.625) (  3d29'30.00"W, 70d17'46.95"S)
>>> Lower Right (  481499.375, 2132000.625) (  3d30'25.10"W, 70d54'54.17"S)
>>> Center      (  452249.375, 2166500.625) (  4d17'17.74"W, 70d36' 6.69"S)
>>> Band 1 Block=46800x1 Type=UInt16, ColorInterp=Gray
>>>
>>> I want to translate it into a GMT grid file using gdal_translate:
>>> gdal_translate -of GMT IMAGE.tif IMAGE.grd
>>>
>>> But this produces an error and the file is not useable afterwards.
>>> ERROR 1: nc_put_vara_double(IMAGE.grd): NetCDF: Start+count exceeds
>>> dimension bound
>>>
>>> gdalinfo  IMAGE.grd
>>> Warning 1: IMAGE.grd is a netCDF file, but not in GMT configuration.
>>> gdalinfo failed - unable to open 'IMAGE.grd'.
>>>
>>> Can anybody explain this and offer a solution?
>>>
>>> Thanks in advance!
>>> Sebastian.
>>>
>>> _______________________________________________
>>> gdal-dev mailing list
>>> gdal-dev at lists.osgeo.org
>>> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev