[GRASSLIST:5812] Problem with raster

David Orme d.orme at imperial.ac.uk
Thu Mar 20 10:38:26 EST 2003


Hello again,

Many thanks for all the suggestions about shell scripts [5801].

I've run into another problem with a floating point raster map in a 
Lat/Long location (1 degree resolution). The map plots fine provided I 
don't use a region which crosses the 180E meridian; if the region does 
cross 180E then I get an error - there's a transcript showing this 
below. A similar integer raster map (which the floating point map is 
calculated from) does not show a problem with any of these regions. The 
reference in the error message ("ERROR: cell_values_d: xdr_double 
failed for index 14.")  to an index seems to vary depending on which 
column lies over 180E.

Any suggestions as to what might be going wrong?

David

> Mapset <dorme> in Location <global>
> GRASS 5.0.0 > g.region n=90 s=-90 w=-180 e=-150; d.erase
>
> Mapset <dorme> in Location <global>
> GRASS 5.0.0 > d.rast TBWrich.corr
>  100%
>
> Mapset <dorme> in Location <global>
> GRASS 5.0.0 > g.region n=90 s=-90 w=150 e=180; d.erase
>
> Mapset <dorme> in Location <global>
> GRASS 5.0.0 > d.rast TBWrich.corr
>  100%
>
> Mapset <dorme> in Location <global>
> GRASS 5.0.0 > g.region n=90 s=-90 w=165 e=-165; d.erase
>
> Mapset <dorme> in Location <global>
> GRASS 5.0.0 > d.rast TBWrich.corr
> ERROR: cell_values_d: xdr_double failed for index 14.
>
> Mapset <dorme> in Location <global>
> GRASS 5.0.0 > g.region n=90 s=-90 w=151 e=-179; d.erase
>
> Mapset <dorme> in Location <global>
> GRASS 5.0.0 > d.rast TBWrich.corr
> ERROR: cell_values_d: xdr_double failed for index 0.
>
> Mapset <dorme> in Location <global>
> GRASS 5.0.0 > g.region n=90 s=-90 w=179 e=-151; d.erase
>
> Mapset <dorme> in Location <global>
> GRASS 5.0.0 > d.rast TBWrich.corr
> ERROR: cell_values_d: xdr_double failed for index 28.




More information about the grass-user mailing list