[gdal-dev] Problem with GDAL data folder in 1.9

[Ivan Lucena]

Hi There,

I am getting some weird results from a simple geotiff file and it seems to be a problem with the data folder.

If I run gdalinfo using GDAL 1.9 or 1.8 with the data folder installed from 1.8 the WKT comes with no EPSG Authority code:

D:\> set GDAL_DATA=D:\GDAL-1.9\data
D:\> gdalinfo sf1.gtif
PROJCS["NAD_1983_StatePlane_California_III_FIPS_0403",
    GEOGCS["GCS_North_American_1983",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS_1980",6378137,298.257222101]],
        PRIMEM["Greenwich",0],
        UNIT["Degree",0.0174532925199432955]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["False_Easting",2000000],
    PARAMETER["False_Northing",500000],
    PARAMETER["Central_Meridian",-120.5],
    PARAMETER["Standard_Parallel_1",37.06666666666667],
    PARAMETER["Standard_Parallel_2",38.43333333333333],
    PARAMETER["Latitude_Of_Origin",36.5],
    UNIT["Meter",1]]

Even tough listgeo shows:

PCS = 26943 (NAD83 / California zone 3)
Projection = 10433 (SPCS83 California zone 3 (meters))
Projection Method: CT_LambertConfConic_2SP
   ProjFalseOriginLatGeoKey: 36.500000 ( 36d30' 0.00"N)
   ProjFalseOriginLongGeoKey: -120.500000 (120d30' 0.00"W)
   ProjStdParallel1GeoKey: 38.433333 ( 38d26' 0.00"N)
   ProjStdParallel2GeoKey: 37.066667 ( 37d 4' 0.00"N)
   ProjFalseEastingGeoKey: 2000000.000000 m
   ProjFalseNorthingGeoKey: 500000.000000 m
GCS: 4269/NAD83
Datum: 6269/North American Datum 1983
Ellipsoid: 7019/GRS 1980 (6378137.00,6356752.31)
Prime Meridian: 8901/Greenwich (0.000000/  0d 0' 0.00"E)
Projection Linear Units: 9001/metre (1.000000m)

But if I use gdalinfo 1.8 or 1.9 with the data folder from 1.8 the WKT appears to be better.

D:\> set GDAL_DATA=D:\GDAL-1.8\data
D:\> gdalinfo sf1.gtif
PROJCS["NAD83 / California zone 3",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",38.43333333333333],
    PARAMETER["standard_parallel_2",37.06666666666667],
    PARAMETER["latitude_of_origin",36.5],
    PARAMETER["central_meridian",-120.5],
    PARAMETER["false_easting",2000000],
    PARAMETER["false_northing",500000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","26943"]]

The authority code are correct and standard_parallel_1, standard_parallel_2 are identical to listgeo information.

I got this files from the SVN branch 1.8 and 1.9 and I compiled and install it myself but I also tried that same steps with the OSGeo4W environment and got the same results. 

Does it means that there is something wrong with the data folder in 1.9 or is that part of a new approach of dealing with that particular kind of projection?

Regards,

Ivan