Add support for new metadata item "coordinate system string" to ENVI driver

[dszlucha]

ENVI 4.7 adds a new metadata item called "coordinate system string" which is an ESRI .prj compatible WKT string to the ENVI .hdr format in envidataset.cpp. I would like to have this new code integrated into the GDAL trunk as well as the ESRI branch. Attached is our working copy of the new code which adds parsing to ENVIDataset::ProcessMapinfo() and writes the new string in ENVIDataset::WriteProjectionInfo().

[warmerdam]

David,

It seems there have been several changes in trunk since the revision you worked with. Could you re-prepare your changes as a patch, or at least apply them to a current trunk version of the source?

Also, I would appreciate a very small dataset demonstrating this new feature.

[dszlucha]

coordinate system string additions applied to trunk version of envidataset.cpp

[Mateusz Łoskot]

@dszlucha

It would be better if you can ​prepare regular patch that consists only your changes. It is as simple as issuing this command in gdal source root directory of your SVN trunk copy:

svn diff > ticket3312.patch

then attach the ticket3312.patch file.

[dszlucha]

svn diff with changes to envidataset.cpp

[warmerdam]

Integrated in 1.6-esri branch by Gao (r18702).

[warmerdam]

David,

I'd still appreciate some sort of sample data before I try integrating this.

[dszlucha]

Sample data with "coordinate system string" metadata

[gaopeng]

Frank, Could you integrate this patch into GDAL 1.8 and GDAL 1.8-esri? ITT-VIS is asking for help. It looks like the sample data has been added.

Thanks, Gao

[warmerdam]

Slightly modified patch applied in trunk (r22303), 1.8 branch (r22304) and 1.8-esri branch (r22305).

[Even Rouault]

I think this is related with this change in the 1.8 branch (works fine in trunk). I now get this failure :

  TEST: envi_4 ... Got: 
PROJCS["Lambert_Conformal_Conic",
    GEOGCS["GCS_North_American_1983",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS_1980",6378137,298.257222101]],
        PRIMEM["Greenwich",0],
        UNIT["Degree",0.017453292519943295]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",33.90363402775256],
    PARAMETER["standard_parallel_2",33.62529002776137],
    PARAMETER["latitude_of_origin",33.76446202775696],
    PARAMETER["central_meridian",-117.4745428888127],
    PARAMETER["false_easting",20000],
    PARAMETER["false_northing",30000],
    PARAMETER["scale_factor",1.0],
    UNIT["Meter",1]]
Expected:
PROJCS["unnamed",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.257222101],
            TOWGS84[0,0,0,0,0,0,0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",33.90363402775256],
    PARAMETER["standard_parallel_2",33.62529002776137],
    PARAMETER["latitude_of_origin",33.76446202775696],
    PARAMETER["central_meridian",-117.4745428888127],
    PARAMETER["false_easting",20000],
    PARAMETER["false_northing",30000]]
fail
    Did not get expected projection reference.