Opened 16 years ago

Last modified 13 years ago

#2550 closed defect

gdal grib driver gets the origin and corner coords wrong — at Version 4

Reported by: winkey Owned by: warmerdam
Priority: normal Milestone:
Component: GDAL_Raster Version: unspecified
Severity: normal Keywords: grib
Cc: erick@…, Daniel Morissette

Description (last modified by warmerdam)

rush@winkey ~ $ gdalinfo gfs.t06z.pgrb2f00
Driver: GRIB/GRIdded Binary (.grb)
Files: gfs.t06z.pgrb2f00
Size is 720, 361
Coordinate System is:
GEOGCS["Coordinate System imported from GRIB file",
    DATUM["unknown",
        SPHEROID["Sphere",6371229,0]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]
Origin = (0.000000000000000,270.000000000000000)
Pixel Size = (0.500000000000000,-0.500000000000000)
Corner Coordinates:
Upper Left  (       0.000,     270.000) (  0d 0'0.01"E,270d 0'0.00"N)
Lower Left  (   0.0000000,  89.5000000) (  0d 0'0.01"E, 89d30'0.00"N)
Upper Right (     360.000,     270.000) (360d 0'0.00"E,270d 0'0.00"N)
Lower Right (     360.000,      89.500) (360d 0'0.00"E, 89d30'0.00"N)
Center      (     180.000,     179.750) (180d 0'0.00"E,179d45'0.00"N)

degrib gfs.t06z.pgrb2f00 -out stdout -C -Csv | less
GDS | Number of Points | 259920
GDS | Projection Type | 0 (Latitude/Longitude)
GDS | Shape of Earth | sphere
GDS | Radius | 6371.229000 (km)
GDS | Nx (Number of points on parallel) | 720
GDS | Ny (Number of points on meridian) | 361
GDS | Lat1 | 90.000000
GDS | Lon1 | 0.000000
GDS | u/v vectors relative to | easterly/northerly
GDS | Lat2 | -90.000000
GDS | Lon2 | 359.500000
GDS | Dx | 0.500000 (degrees)
GDS | Dy | 0.500000 (degrees)
GDS | Input GRIB2 grid, scan mode | 64 (0100)
GDS | Output grid, scan mode | 64 (0100)
GDS | (.flt file grid), scan mode | 0 (0000)
GDS | Output grid, scan i/x direction | positive
GDS | Output grid, scan j/y direction | positive
GDS | (.flt file grid), scan j/y direction | negative
GDS | Output grid, consecutive points in | i/x direction
GDS | Output grid, adjacent rows scan in | same direction

the dataset can be found at http://atmos.ucsd.edu/gfs.t06z.pgrb2f00
i did not attach it because of its size
-rw-r--r-- 1 rush users 42806506 May 19 02:28 gfs.t06z.pgrb2f00

Change History (4)

comment:2 by Even Rouault, 16 years ago

Reproduced with GDAL 1.6.0dev

comment:3 by warmerdam, 15 years ago

Component: defaultGDAL_Raster
Keywords: grib added
Status: newassigned

comment:4 by warmerdam, 15 years ago

Description: modified (diff)
Note: See TracTickets for help on using tickets.