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

[Gustavo]

I've also asked in the GDAL list but I had no solution. My system was Ubuntu 
9.10 and yesterday I managed to upgrade it to 10.04 and the problem got away.

Hamish,
regarding the modis data, I am using the daily 4km data. What I want to do is 
to relate some VMS (fishing vessel position) with catch and some environmental 
data. I am using a postgis database combined with some processing with GRASS.

Gustavo

On Vie 03 Dic 2010 18:04:35 usted escribió:
> Gustavo wrote:
> > I am using the SMI format. I've tried the following with
> > the same result.
> 
> that should be right, it's the L3m SMI/HDF download, not the
> BIN/HDF one.
> 
> >  gdal_translate -a_srs "+init=epsg:4326" -a_nodata 65535
> > 
> > -a_ullr -180 90 180
> > -90 -co "COMPRESS=PACKBITS"
> > HDF4_SDS:UNKNOWN:"A20100252010032.L3m_8D_CHL_chlor_a_9km":0
> > prueba.tif
> > Input file size is 4320, 2160
> > 0ERROR 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
> 
> I'm working with this data reasonably often and I've never seen
> that error, sorry. As the error is coming from GDAL maybe ask
> on their mailing list? maybe someone there has a better idea.
> I'm going to process some new chlorophyll-a modis data in the
> next week, I'll let you know if I run into the same problem.
> 
> Eduardo Klein escribió:
> > > 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)
> 
> good luck,
> Hamish
> 
> ps- it defaults to 9km, but on the ocean color website you can
> tell it to download the 4km resolution too :) (but I usually get
> the much smaller 9km data first and check for cloud cover, etc)