[GRASS-dev] Possible bug in r.out.gdal?

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

Soeren Gebbert wrote:

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

Using grass6.4svn with gdal1.6.0 here

Exporting elevation from nc_spm_08 with

r.out.gdal input=elevation@PERMANENT format=DTED type=Float32 output=test.dted

gives me

Driver <DTED> does not support direct writing. Using MEM driver for
intermediate dataset.
Exporting to GDAL data type: Float32
100%
Warning 1: The source does not appear to be a properly formatted cell.
Warning 1: The source projection coordinate system is PROJCS["Lambert Conformal Conic",GEOGCS["grs80",DATUM["North_American_Datum_1983",SPHEROID["Geodetic_Reference_System_1980",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic_2SP"],PARAMETER["standard_parallel_1",36.16666666666666],PARAMETER["standard_parallel_2",34.33333333333334],PARAMETER["latitude_of_origin",33.75],PARAMETER["central_meridian",-79],PARAMETER["false_easting",609601.22],PARAMETER["false_northing",0],UNIT["Meter",1]]. Only WGS 84 is supported.
The DTED driver will generate a file as if the source was WGS 84 projection coordinate system.
Warning 1: The corner coordinates of the source are not properly aligned on plain latitude/longitude boundaries.
Warning 1: The horizontal source size is not conformant with the one expected by DTED Level 1 at this latitude (1500 pixels found instead of 225).
r.out.gdal complete.

These steps are performed:

GRASS 6.4.0svn > g.region -p
projection: 1 (UTM)

At least in gdal1.6.0, only WGS84 is supported for DTED, see warning above. WGS72 is also supported, but the gdal driver issues a warning.

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.

These messages are generated by gdal, not r.out.gdal. See
http://www.gdal.org/frmt_dted.html
for more info on the DTED format, particularly the section Creation Issues and the link to the DTED specs:
http://www.nga.mil/ast/fm/acq/89020B.pdf
If you can confirm that the input data are exactly formatted as Level 0, 1 or 2 cell, it is a problem of r.out.gdal. Also note that for DTED, elevation data is returned as 16 bit signed integer (citation from gdal).

hope that helps,

Markus M

From the DTED specs:

A data file of DTED Level 0, DTED Level 1 or
DTED Level 2 is a 1° by 1° cell defined by whole degree latitude and
longitude
lines on WGS. A DTED file shall not cross whole degree latitude or longitude
lines.

Matrix intervals for DTED Level 0.
ZONE LATITUDE latitude x longitude INTERVAL
I 0° - 50° North-South 30 x 30 seconds
II 50° - 70° North-South 30 x 60 seconds
III 70° - 75° North-South 30 x 90 seconds
IV 75° - 80° North-South 30 x 120 seconds
V 80° - 90° North-South 30 x 180 seconds

Matrix intervals for DTED Level 1.
ZONE LATITUDE latitude x longitude INTERVAL
I 0° - 50° North-South 3 x 3 seconds
II 50° - 70° North-South 3 x 6 seconds
III 70° - 75° North-South 3 x 9 seconds
IV 75° - 80° North-South 3 x 12 seconds
V 80° - 90° North-South 3 x 18 seconds

Matrix intervals for DTED Level 2.
ZONE LATITUDE latitude x longitude INTERVAL
I 0° - 50° North-South 1 x 1 seconds
II 50° - 70° North-South 1 x 2 seconds
III 70° - 75° North-South 1 x 3 seconds
IV 75° - 80° North-South 1 x 4 seconds
V 80° - 90° North-South 1 x 6 seconds

No idea how strict the gdal driver for DTED is, but it seems that it is
very strict and won't work without reprojection, followed by setting
region extends and resolution according to the desired DTED level.

Markus M

Soeren Gebbert wrote:

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 <http://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

Hello Markus,
many thanks for your response.
I guess the behavior of r.out.gdal in case of a raster map in UTM projection and DTED export seems to be correct.
So i removed the DTED test from the GRASS testsuite, based on your suggestions and the gdal documentation.

soeren

2009/8/11 Markus Metz <markus.metz.giswork@googlemail.com>

From the DTED specs:

A data file of DTED Level 0, DTED Level 1 or
DTED Level 2 is a 1° by 1° cell defined by whole degree latitude and
longitude
lines on WGS. A DTED file shall not cross whole degree latitude or longitude
lines.

Matrix intervals for DTED Level 0.
ZONE LATITUDE latitude x longitude INTERVAL
I 0° - 50° North-South 30 x 30 seconds
II 50° - 70° North-South 30 x 60 seconds
III 70° - 75° North-South 30 x 90 seconds
IV 75° - 80° North-South 30 x 120 seconds
V 80° - 90° North-South 30 x 180 seconds

Matrix intervals for DTED Level 1.
ZONE LATITUDE latitude x longitude INTERVAL
I 0° - 50° North-South 3 x 3 seconds
II 50° - 70° North-South 3 x 6 seconds
III 70° - 75° North-South 3 x 9 seconds
IV 75° - 80° North-South 3 x 12 seconds
V 80° - 90° North-South 3 x 18 seconds

Matrix intervals for DTED Level 2.
ZONE LATITUDE latitude x longitude INTERVAL
I 0° - 50° North-South 1 x 1 seconds
II 50° - 70° North-South 1 x 2 seconds
III 70° - 75° North-South 1 x 3 seconds
IV 75° - 80° North-South 1 x 4 seconds
V 80° - 90° North-South 1 x 6 seconds

No idea how strict the gdal driver for DTED is, but it seems that it is
very strict and won’t work without reprojection, followed by setting
region extends and resolution according to the desired DTED level.

Markus M

Soeren Gebbert wrote:

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 <http://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

Soeren wrote:

| Data Description:
| generated by r.mapcalc

....

| sin(col() * 5) * col() + cos(row() * 5) * row() + 200

nice test pattern

H

Soeren Gebbert wrote:

Hello Markus,
many thanks for your response.
I guess the behavior of r.out.gdal in case of a raster map in UTM projection and DTED export seems to be correct.

To be precise, it's the behaviour of the gdal library. r.out.gdal leaves the format-specific tests to the gdal library, and it was the gdal library, not r.out.gdal, that issued the ERROR 5: messages. There is nothing that can be done in r.out.gdal to avoid this kind of error.

Markus M

So i removed the DTED test from the GRASS testsuite, based on your suggestions and the gdal documentation.

soeren

2009/8/11 Markus Metz <markus.metz.giswork@googlemail.com <mailto:markus.metz.giswork@googlemail.com>>

    From the DTED specs:

    A data file of DTED Level 0, DTED Level 1 or
    DTED Level 2 is a 1° by 1° cell defined by whole degree latitude and
    longitude
    lines on WGS. A DTED file shall not cross whole degree latitude or
    longitude
    lines.

    Matrix intervals for DTED Level 0.
    ZONE LATITUDE latitude x longitude INTERVAL
    I 0° - 50° North-South 30 x 30 seconds
    II 50° - 70° North-South 30 x 60 seconds
    III 70° - 75° North-South 30 x 90 seconds
    IV 75° - 80° North-South 30 x 120 seconds
    V 80° - 90° North-South 30 x 180 seconds

    Matrix intervals for DTED Level 1.
    ZONE LATITUDE latitude x longitude INTERVAL
    I 0° - 50° North-South 3 x 3 seconds
    II 50° - 70° North-South 3 x 6 seconds
    III 70° - 75° North-South 3 x 9 seconds
    IV 75° - 80° North-South 3 x 12 seconds
    V 80° - 90° North-South 3 x 18 seconds

    Matrix intervals for DTED Level 2.
    ZONE LATITUDE latitude x longitude INTERVAL
    I 0° - 50° North-South 1 x 1 seconds
    II 50° - 70° North-South 1 x 2 seconds
    III 70° - 75° North-South 1 x 3 seconds
    IV 75° - 80° North-South 1 x 4 seconds
    V 80° - 90° North-South 1 x 6 seconds

    No idea how strict the gdal driver for DTED is, but it seems that
    it is
    very strict and won't work without reprojection, followed by setting
    region extends and resolution according to the desired DTED level.

    Markus M

    Soeren Gebbert wrote:

        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 <http://r.info> <http://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 <mailto:grass-dev@lists.osgeo.org>
        http://lists.osgeo.org/mailman/listinfo/grass-dev