[gdal-dev] gdalinfo, map unit DD and epsg:9108

yrrman yrrman yrrman at googlemail.com
Thu Nov 11 11:09:12 EST 2010


Hi everyone,

I'm new to this list and would therefore like to say "hello and thanks
for gdal" in the first place :-)


I just came across a minor problem where I would like to ask for your advice -

the Global Landcover set GLC2000 is available in grid format from [0].
Associated with it is the .hdr file with the projection description
[2].

Apparently the map units are in decimal degree.

When I ask gdalinfo, I get the WKT description of the SRS, see [3]


Everything's fine with this, I only came across the UNIT-Field

   UNIT["degree",0.0174532925199433, AUTHORITY["EPSG","9108"]],

where the unit is degree, defined according to epsg:9108. When you
look up epsg:9108 in the epsg registry [1] you get as Unit of Measure:

Code: EPSG::9108 - Name: degree minute second hemisphere (i.e. not
decimal degree)

Shouldn't this rather be epsg:9102 then?
Code: EPSG::9102 Name: degree, 1 degree: 0.017453292519943278 radian



As I said, it's not very important, just out of curiousity.

I'm using
GDAL 1.7.0b2, FWTools 2.4.7, released 2010/01/19


Thanks a lot and kind regards,


Chris


-------------------------------------------------------------------------------------
[0] http://bioval.jrc.ec.europa.eu/products/glc2000/products.php
[1] http://www.epsg-registry.org/
[2] glc2000_v1_1_projinfo.hdr:

BANDS:      1
ROWS:   16353
COLS:   40320
INTERLEAVING:    NA
DATATYPE: U8
BYTE_ORDER: NA
UL_X_COORDINATE:           -180.000000
UL_Y_COORDINATE:             89.991071
LR_X_COORDINATE:            179.991070
LR_Y_COORDINATE:            -56.008928
PIXEL_WIDTH:               0.0089285714
PIXEL_HEIGHT:              0.0089285714
MAP_UNITS:    DD
PROJECTION_NAME: Geographic (Lat/Lon)
PROJECTION_ZONE:     0
PROJECTION_PARAMETERS:
                    0.000000
                    0.000000
                    0.000000
                    0.000000
                    0.000000
                    0.000000
                    0.000000
                    0.000000
                    0.000000
                    0.000000
                    0.000000
                    0.000000
                    0.000000
                    0.000000
                    0.000000
SPHEROID_NAME: WGS 84
SEMI_MAJOR_AXIS:        6378137.000000
SEMI_MINOR_AXIS:        6356752.314245
E_SQUARED:                    0.006694
RADIUS:                 6371007.180918

[3] gdalinfo output:
Driver: AIG/Arc/Info Binary Grid
Files: glc2000_v1_1
      glc2000_v1_1\dblbnd.adf
      glc2000_v1_1\hdr.adf
      glc2000_v1_1\log
      glc2000_v1_1\metadata.xml
      glc2000_v1_1\prj.adf
      glc2000_v1_1\sta.adf
      glc2000_v1_1\vat.adf
      glc2000_v1_1\w001001.adf
      glc2000_v1_1\w001001x.adf
Size is 40320, 16353
Coordinate System is:
GEOGCS["WGS 84",
   DATUM["WGS_1984",
       SPHEROID["WGS 84",6378137,298.257223563,
           AUTHORITY["EPSG","7030"]],
       TOWGS84[0,0,0,0,0,0,0],
       AUTHORITY["EPSG","6326"]],
   PRIMEM["Greenwich",0,
       AUTHORITY["EPSG","8901"]],
   UNIT["degree",0.0174532925199433,
       AUTHORITY["EPSG","9108"]],
   AUTHORITY["EPSG","4326"]]
Origin = (-180.004464285700010,89.995535666300043)
Pixel Size = (0.008928571400000,-0.008928571400000)
Corner Coordinates:
Upper Left  (-180.0044643,  89.9955357) (180d 0'16.07"W, 89d59'43.93"N)
Lower Left  (-180.0044643, -56.0133924) (180d 0'16.07"W, 56d 0'48.21"S)
Upper Right ( 179.9955346,  89.9955357) (179d59'43.92"E, 89d59'43.93"N)
Lower Right ( 179.9955346, -56.0133924) (179d59'43.92"E, 56d 0'48.21"S)
Center      (  -0.0044649,  16.9910716) (  0d 0'16.07"W, 16d59'27.86"N)
Band 1 Block=256x4 Type=Byte, ColorInterp=Palette
 Min=1.000 Max=23.000
 NoData Value=255
 Metadata:
   LAYER_TYPE=thematic
 Color Table (RGB with 256 entries)
   0: 0,0,0,255 [..]


More information about the gdal-dev mailing list