[Gdal-dev] projection with ogr2ogr

[Lijian Shi]

Hi,

 

I have a shp file of point with longitude, latitude and other attribution.
And I have a geotiff file with the followed information:

 

Driver: GTiff/GeoTIFF

Size is 256, 256

Coordinate System is:

PROJCS[" Projection Name = Lambert Azimuthal Equal-ar",

    GEOGCS["WGS 84",

        DATUM["WGS_1984",

            SPHEROID["WGS 84",6378137,298.2572235629972,

                AUTHORITY["EPSG","7030"]],

            AUTHORITY["EPSG","6326"]],

        PRIMEM["Greenwich",0],

        UNIT["degree",0.0174532925199433],

        AUTHORITY["EPSG","4326"]],

    PROJECTION["Lambert_Azimuthal_Equal_Area"],

    PARAMETER["latitude_of_center",45],

    PARAMETER["longitude_of_center",-100],

    PARAMETER["false_easting",0],

    PARAMETER["false_northing",0],

    UNIT["metre",1,

        AUTHORITY["EPSG","9001"]]]

 

Now I want to use ogr2ogr to project the shp file in order that the shp file
has the same geoinformation(same projection, same datum) with the geotiff
file. I use the command like: 

ogr2ogr -s_srs WGS84 -t_srs EPSG:code utm2.shp point.shp

 

Am I right? If right, how can I get the code? If wrong, how should I use the
"ogr2ogr"?

 

Thank you very much!

 

Lijian

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20070306/158ffc19/attachment.html