[GRASS5] Import NOAA AVHRR .dat files

GRASS GIS GRASS GIS developer
1

Hello.

Im trying to import a .dat file with r.in.bin.

r.in.bin input=/home/neb_15.dat output=neb_exp bytes=1 north=44.80
south=35.01 east=-1.14 west=-12.20 rows=89 cols=79 -b

However, I always get an error message like this:

Bytes do not match File size

File size 171262... Total Bytes 7031

Try bytes=24 or adjusting input parameters

The file header looks like this

#LLMATRIX,

NORTH = 44.80,

WEST = 347.80,

GRID = 0.14/0.11,

NLAT = 89,

NLON = 79

#DATA

Does anyone have a solution?

Thanks in advance,

Ricardo

DISCLAIMER: This message may contain confidential information or privileged material and is intended only for the individual(s) named. If you are not a named addressee and mistakenly received this message you should not copy or otherwise disseminate it: please delete this e-mail from your system and notify the sender immediately. E-mail transmissions are not guaranteed to be secure or error-free as information could be intercepted, corrupted, lost, destroyed, arrive late or incomplete or contain viruses. Therefore, the sender does not accept liability for any errors or omissions in the contents of this message that arise as a result of e-mail transmissions. Please request a hard copy version if verification is required. Critical Software.

2

On Wed, May 04, 2005 at 05:08:50PM +0100, Ricardo Armas wrote:

Hello.
Im trying to import a .dat file with r.in.bin.

r.in.bin input=/home/neb_15.dat output=neb_exp bytes=1 north=44.80
south=35.01 east=-1.14 west=-12.20 rows=89 cols=79 -b

However, I always get an error message like this:

Bytes do not match File size

File size 171262... Total Bytes 7031

Try bytes=24 or adjusting input parameters

The file header looks like this

#LLMATRIX,
NORTH = 44.80,
WEST = 347.80,
GRID = 0.14/0.11,
NLAT = 89,
NLON = 79
#DATA

Does anyone have a solution?

Calculating:

7031/(89 * 79)

[1] 1

suggests 1 byte. But the boundary box above doesn't seem to
match this:

g.region n=44.80 s=35.01 e=-1.14 w=-12.20 -p res=0:08
projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: a=6378137 es=0.00669438
north: 44:48N
south: 35:00:36N
west: 12:12W
east: 1:08:24W
nsres: 0:08:02.794521
ewres: 0:07:59.710843
rows: 73
cols: 83

I get close to 89x79 playing with 'res=' but don't get it
precisely. Are you sure about the n,s,e,w values?

Also: what's the projection of these data (I tried latLong/WGS84)?
Often it's Goode or something else.

Markus