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