[Gdal-dev] problem abotu reprojection with ogr2ogr

[Lijian Shi]

Hi, everyone,

 

I have a shp file and landsat data of the same field, but these data's
projections are different. The projection information of shp file is below:

 

GEOGCS["GCS_Assumed_Geographic_1",DATUM["D_North_American_1927",SPHEROID["Cl
arke_1866",6378206.4,294.9786982]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0
174532925199433]]

 

We can get this information from .prj file or use command 'ogrinfo' .

 

The projection information of landsat data is:

 

PROJCS["WGS 84 / UTM zone 14N",

    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["Transverse_Mercator"],

    PARAMETER["latitude_of_origin",0],

    PARAMETER["central_meridian",-99],

    PARAMETER["scale_factor",0.9996],

    PARAMETER["false_easting",500000],

    PARAMETER["false_northing",0],

    UNIT["metre",1,

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

AUTHORITY["EPSG","32614"]]

 

I copy this information into a empty file, then saved as 'landsat.prj'.

 

Now I use 'ogr2ogr' to reproject the shp file in order that the shp file has
the same projection information with landsat data. The command is as below:

 

ogr2ogr -s_srs input.prj -t_srs landsat.prj output.shp input.shp

 

When I put the output.shp and input.shp into ArcGIS together, there a shift
between two file. Dose anyone know the reason? Is there any problem with my
command?

 

Thanks in advance.

 

Lijian

 

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20070508/3516dae4/attachment.html