id	summary	reporter	owner	description	type	status	priority	milestone	component	version	severity	resolution	keywords	cc
3757	Problem with sphere representation in the RST driver	sjurk	warmerdam	"
Trying to convert a grib2 file to IDRISI, extracting band 1, 
the PROJCS information is distorted.

Windows 7, GDAL 1.7.0b2, FWTools 2.4.7, released 2010/01/19

I use the following sequence of commands:

gdalinfo har04_00_ext.grib2>gribinfo.txt

gdal_translate -of RST -b 1 har04_00_ext.grib2 har04_00_ext.rst
Input file size is 289, 489
0...10...20...30...40...50...60...70...80...90...100 - done.

gdalinfo har04_00_ext.rst>rstinfo.txt


Output of first GDALINFO:

Driver: GRIB/GRIdded Binary (.grb)
Files: har04_00_ext.grib2
Size is 289, 489
Coordinate System is:
PROJCS[""unnamed"",
    GEOGCS[""Coordinate System imported from GRIB file"",
        DATUM[""unknown"",
            SPHEROID[""Sphere"",6367470,0]],
        PRIMEM[""Greenwich"",0],
        UNIT[""degree"",0.0174532925199433]],
    PROJECTION[""Lambert_Conformal_Conic_2SP""],
    PARAMETER[""standard_parallel_1"",64.35],
    PARAMETER[""standard_parallel_2"",64.35],
    PARAMETER[""latitude_of_origin"",0],
    PARAMETER[""central_meridian"",337],
    PARAMETER[""false_easting"",0],
    PARAMETER[""false_northing"",0]]
Origin = (1192906.028257467600000,10107044.152202457000000)
Pixel Size = (4000.000000000000000,-4000.000000000000000)
Corner Coordinates:
Upper Left  ( 1192906.028,10107044.152) ( 19d40'29.23""E, 74d32'51.23""N)
Lower Left  ( 1192906.028, 8151044.152) (  1d52'57.71""W, 58d57'40.03""N)
Upper Right ( 2348906.028,10107044.152) ( 40d41'28.49""E, 66d46'40.32""N)
Lower Right ( 2348906.028, 8151044.152) ( 14d55'47.70""E, 54d18'2.03""N)
Center      ( 1770906.028, 9129044.152) ( 16d25'18.42""E, 64d26'49.23""N)
Band 1 Block=289x1 Type=Float64, ColorInterp=Undefined
Etc... (long file)


Output of second GDALINFO, note the “-1.#INF” instead of 0 under the SPHEROID label.

Driver: RST/Idrisi Raster A.1
Files: har04_00_ext.rst
       har04_00_ext.rst.aux.xml
       har04_00_ext.rdc
       har04_00_ext.ref
Size is 289, 489
Coordinate System is:
PROJCS[""unnamed"",
    GEOGCS[""har04_00_ext"",
        DATUM[""unknown"",
            SPHEROID[""Sphere"",6367470,""-1.#INF""],
            TOWGS84[0,0,0,0,0,0,0]],
        PRIMEM[""Greenwich"",0],
        UNIT[""degree"",0.0174532925199433]],
    PROJECTION[""Lambert_Conformal_Conic_2SP""],
    PARAMETER[""standard_parallel_1"",64.35],
    PARAMETER[""standard_parallel_2"",64.35],
    PARAMETER[""latitude_of_origin"",0],
    PARAMETER[""central_meridian"",337],
    PARAMETER[""false_easting"",0],
    PARAMETER[""false_northing"",0],
    UNIT[""degree"",0]]
Origin = (1192906.028257499900000,10107044.152202500000000)
Pixel Size = (4000.000000000000000,-4000.000000000000000)
Corner Coordinates:
Upper Left  ( 1192906.028,10107044.152) 
Lower Left  ( 1192906.028, 8151044.152) 
Upper Right ( 2348906.028,10107044.152) 
Lower Right ( 2348906.028, 8151044.152) 
Center      ( 1770906.028, 9129044.152) 
Band 1 Block=289x1 Type=Float32, ColorInterp=Gray


"	defect	closed	normal	1.8.1	GDAL_Raster	1.7.0	normal	fixed	RST sphere