[GRASS-user] modis ocean and r.in.gdal

[Eduardo Klein]

Gustavo, are you working with the SMI/HDF product? I've just run the 
script in grass 6.4.0 and it worked perfect. Check with gdalinfo the 
name of the subset to process (usually something like 
HDF4_SDS:UNKNOWN:"A20101292010136.L3m_8D_SST_9":0)

Saludos, EKS


On 12/02/2010 06:19 PM, Gustavo wrote:
> Hola Eduardo,
>
> although the error message corresponded to the grass6.5 r.in.gdal, when I
> tried method 1 with gdal_translate, the error was the same.
>
> Gustavo
>
> On Jue 02 Dic 2010 19:24:04 Eduardo Klein escribió:
>> Hola Gustavo,
>>
>> It seems to be a problem with the projection and the boundaries. The
>> previous gdal_translate command ran OK?
>>
>> Saludos, EKS
>>
>> On 12/02/2010 05:34 PM, Gustavo wrote:
>>> Hi everyone,
>>>
>>> I am trying to import some L3 MODIS satellite images into GRASS 6.4 by
>>> following the instructions from:
>>>
>>> http://grass.osgeo.org/wiki/MODIS#Import
>>>
>>>    I've tried method 1 and 2 but I got some error messages.
>>> I've installed the grass6.5 snapshot from 27/11 and when run the
>>> following command:
>>>
>>> r.in.gdal -l -o in=HDF4_SDS:UNKNOWN:A2010021.L3m_DAY_NSST_4:0
>>> out=A2010021.L3m_DAY_NSST_4
>>>
>>> I got some errors:
>>>
>>> WARNING: Map bounds have been constrained to geographic coordinates. You
>>>            will almost certainly want to check map bounds and resolution
>>> with r.info and reset them with r.region before going any further.
>>> WARNING: Over-riding projection check
>>> ERROR 1: SDreaddata() failed for block.
>>> ERROR 1: IReadBlock failed at X offset 0, Y offset 0
>>> ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0
>>> ERROR 1: SDreaddata() failed for block.
>>> ERROR 1: IReadBlock failed at X offset 0, Y offset 1
>>> ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 1
>>> ERROR 1: SDreaddata() failed for block.
>>> .
>>> .
>>> .
>>> ERROR 1: IReadBlock failed at X offset 0, Y offset 332
>>> ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 332
>>> ERROR 1: SDreaddata() failed for block.
>>> More than 1000 errors or warnings have been reported. No more will be
>>> reported from now.
>>>    100%
>>> r.in.gdal complete. Raster map<A2010021.L3m_DAY_NSST_4>   created.
>>>
>>> When I loaded the raster and tried r.what, I got 0 everywhere.
>>>
>>> What can I do?
>>> Thanks,
>>>
>>> Gustavo Martínez
>>> _______________________________________________
>>> grass-user mailing list
>>> grass-user at lists.osgeo.org
>>> http://lists.osgeo.org/mailman/listinfo/grass-user
>>
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user

-- 
Eduardo Klein
Lab Sensores Remotos e
Instituto de Tecnología y Ciencias Marinas
Universidad Simón Bolívar
Caracas, Venezuela
ph (58) (212) 906-3111 ext 6700
fax (58) (212( 906-3111 ext 6701