Dear Devs,
i found a strange behavior in r.out.gdal while updating the GRASS test suite.
Im using the grass6.4svn snapshot: grass-6.4.svn_src_snapshot_2009_08_01.tar.gz.
on openSuse 10.3 Linux 2.6.22.19-0.2-bigsmp
gdal version 1.6.1
The test i am implementing tries to export a DCELL raster map (named elevation ) to a DTED file with r.out.gdal.
r.out.gdal reports plenty of errors (all of the same type) an exits with return code 0.
A dted file is created, but im not sure if it is valid.
Exporting the same file with r.out.gdal in ENVI and ELAS formats seems to work correctly.
These steps are performed:
GRASS 6.4.0svn > g.region -p
projection: 1 (UTM)
zone: 32
datum: wgs84
ellipsoid: wgs84
north: 5582000
south: 5580000
west: 400000
east: 402000
nsres: 20
ewres: 20
rows: 100
cols: 100
cells: 10000
GRASS 6.4.0svn > r.info elevation
±---------------------------------------------------------------------------+
| Layer: elevation Date: Sun Apr 30 12:48:58 2006 |
| Mapset: PERMANENT Login of Creator: soeren |
| Location: TestLocation |
| DataBase: /1/gebbert/src/GRASS_Testsuite |
| Title: ( elevation ) |
Timestamp: none |
---|
Type of Map: raster Number of Categories: 255 |
Data Type: DCELL |
Rows: 100 |
Columns: 100 |
Total Cells: 10000 |
Projection: UTM (zone 32) |
N: 5582000 S: 5580000 Res: 20 |
E: 402000 W: 400000 Res: 20 |
Range of data: min = 68.246322 max = 363.529491 |
Data Description: |
generated by r.mapcalc |
Comments: |
sin(col() * 5) * col() + cos(row() * 5) * row() + 200 |
±---------------------------------------------------------------------------+ |
GRASS 6.4.0svn > r.out.gdal --verbose format=DTED type=Float32 input=elevation output=/tmp/test.dted
…
ERROR 5: Access window out of range in RasterIO(). Requested
(0,994) of size 201x1 on raster of 100x100.
ERROR 5: Access window out of range in RasterIO(). Requested
(0,995) of size 201x1 on raster of 100x100.
ERROR 5: Access window out of range in RasterIO(). Requested
(0,996) of size 201x1 on raster of 100x100.
More than 1000 errors or warnings have been reported. No more will be reported from now.
r.out.gdal complete.
GRASS 6.4.0svn > gdalinfo -nomd -noct -stats /tmp/test.dted
Driver: DTED/DTED Elevation Raster
Files: /tmp/test.dted
/tmp/test.dted.aux.xml
Size is 201, 1201
Coordinate System is:
GEOGCS[“WGS 84”,
DATUM[“WGS_1984”,
SPHEROID[“WGS 84”,6378137,298.257223563]],
PRIMEM[“Greenwich”,0],
UNIT[“degree”,0.0174532925199433],
AUTHORITY[“EPSG”,“4326”]]
Origin = (399.997500000000002,559.000416666666638)
Pixel Size = (0.005000000000000,-0.000833333333333)
Corner Coordinates:
Upper Left ( 399.998, 559.000) (399d59’51.00"E,559d 0’1.50"N)
Lower Left ( 399.998, 558.000) (399d59’51.00"E,557d59’58.50"N)
Upper Right ( 401.002, 559.000) (401d 0’9.00"E,559d 0’1.50"N)
Lower Right ( 401.002, 558.000) (401d 0’9.00"E,557d59’58.50"N)
Center ( 400.500, 558.500) (400d30’0.00"E,558d30’0.00"N)
Band 1 Block=1x1201 Type=Int16, ColorInterp=Palette
Min=0.000 Max=0.000
Minimum=0.000, Maximum=0.000, Mean=0.000, StdDev=0.000
NoData Value=-32767
Unit Type: m
Color Table (RGB with 365 entries)
Any suggestions what i did wrong?
Maybe the region settings are incompatible with the dted output coordinate system?
Or is this a correct behavior of r.out.gdal?
Thanks in advance
best regards
Soeren