Opened 14 years ago
Closed 13 years ago
#3757 closed defect (fixed)
Problem with sphere representation in the RST driver
Reported by: | sjurk | Owned by: | warmerdam |
---|---|---|---|
Priority: | normal | Milestone: | 1.8.1 |
Component: | GDAL_Raster | Version: | 1.7.0 |
Severity: | normal | Keywords: | RST sphere |
Cc: |
Description
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]],
PROJECTIONLambert_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]],
PROJECTIONLambert_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
Change History (1)
comment:1 by , 13 years ago
Milestone: | → 1.8.1 |
---|---|
Resolution: | → fixed |
Status: | new → closed |
Fixed in trunk (r22445) and in branches/1.8 (r22446)