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
_______________________________________________ grass-dev mailing list grass-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-dev