[gdal-dev] Get EPSG code from WKT

[Gane R]

Here is a link to the GeoPDF

https://drive.google.com/file/d/0B4shTM1bPS5TTk00aGIwUkhfd2s/view?usp=sharing

gdalwarp.exe -t_srs EPSG:3857 "Landuse(P)_VNY.pdf" "Landuse(P)_VNY.tif"

Resulting TIF is offset from the actual location, how do I find the correct
EPSG code for this kind of PDF.

but when -s_srs EPSG:26945 is part of the command line it works fine, how
do I know the the correct EPSG code from the WKT as there are more than one
AUTHORITY code with EPSG value ?

gdalinfo Landuse(P)_VNY.pdf

Driver: PDF/Geospatial PDF
Files: Landuse(P)_VNY.pdf
Size is 1275, 1650
Coordinate System is:
PROJCS["NAD83 / California zone 5",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.257222101,
                AUTHORITY["EPSG","7019"]],
            TOWGS84[0,0,0,0,0,0,0],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",35.46666666666667],
    PARAMETER["standard_parallel_2",34.03333333333333],
    PARAMETER["latitude_of_origin",33.5],
    PARAMETER["central_meridian",-118],
    PARAMETER["false_easting",609601.2192024384],
    PARAMETER["false_northing",152400.3048006096],
    UNIT["Meter",1],
    AXIS["X",EAST],
    AXIS["Y",NORTH],
    AUTHORITY["EPSG","26945"]]
Origin = (1954721.377004246000000,580804.857233627580000)
Pixel Size = (6.774914195278887,-6.774914195278887)
Metadata:
  CREATION_DATE=D:20150428111946Z
  CREATOR=TerraGo Technologies Publisher for ArcGIS GeoPDF Export version:
6.5.1.12
  NEATLINE=POLYGON ((1963054.9640081285
571963.11863109597,1955025.6468122969 571963.11863109597,1955025.6468122969
580500.62045958859,1963054.9640081285 580500.62045958859,1963054.9640081285
571963.11863109597))
Corner Coordinates:
Upper Left  ( 1954721.377,  580804.857) (102d56'40.35"W, 36d27' 6.25"N)
Lower Left  ( 1954721.377,  569626.249) (102d57'47.23"W, 36d21' 7.77"N)
Upper Right ( 1963359.393,  580804.857) (102d50'57.53"W, 36d26'24.32"N)
Lower Right ( 1963359.393,  569626.249) (102d52' 4.82"W, 36d20'25.88"N)
Center      ( 1959040.385,  575215.553) (102d54'22.49"W, 36d23'46.09"N)
Band 1 Block=1275x1 Type=Byte, ColorInterp=Red
Band 2 Block=1275x1 Type=Byte, ColorInterp=Green
Band 3 Block=1275x1 Type=Byte, ColorInterp=Blue

On Thu, Feb 18, 2016 at 7:18 PM, Jukka Rahkonen <
jukka.rahkonen at maanmittauslaitos.fi> wrote:

> Gane R <gane.prog <at> gmail.com> writes:
>
>
> > When I try to convert to different projection using the gdalwarp without
> specifying the -s_srs option it reprojects to a different location, for
> some
> raster it projects to a correct location.
> >
> > But when try to run gdalinfo and I see a lot of EPSG code, when I use the
> last AUTHORITY EPSG code as -s_srs option the above mention problem is
> solved.
> >
> > Is it that I have use the last AUTHORITY EPSG code or any thing will help
> >
> > Expecting your suggestion
>
> Expecting some sample data that could be used for reproducing the problem.
>
> -Jukka Rahkonen-
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20160219/9ac55c16/attachment.html>