[gdal-dev] Issue reprojecting geo tiff

Wed Aug 14 07:00:39 PDT 2013

Hi,

Gdalinfo says that the coordinate system is
["UTM Zone 17, Northern Hemisphere"

TFW file from another topo map I captured looks like the units are in meters as you can see

0.9995919650
-0.0204382296
-0.0204382296
-0.9995919650
339224.8572580147
4472388.0032269871

You were telling gdalwarp that the input is in EPSG:4269 which had degrees as units but it is not right. It may be EPSG:26917 but check it
http://spatialreference.org/ref/epsg/26917/

-Jukka Rahkonen-

Lähettäjä: Reaves, Timothy [mailto:treaves at silverfieldstech.com]
Lähetetty: 14. elokuuta 2013 16:23
Vastaanottaja: Rahkonen Jukka
Kopio: gdal-dev at lists.osgeo.org
Aihe: Re: [gdal-dev] Issue reprojecting geo tiff

For PA, it gives:
Driver: PDF/Geospatial PDF
Files: PA_Tionesta_20130621_TM_geo.pdf
Size is 13652, 17403
Coordinate System is:
PROJCS["UTM Zone 17, Northern Hemisphere",
    GEOGCS["unknown",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.257222101],
            TOWGS84[0.9738,-1.9453,-0.5486,0,0,0,0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-81],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]
GeoTransform =
  623457.1927497694, 1.015672572042695, 0.01831264025401483
  4596843.380742616, 0.01831264025401483, -1.015672572042695
Metadata:
  AUTHOR=USGS National Geospatial Technical Operations Center
  CREATION_DATE=D:20130621093456Z
  CREATOR=ESRI ArcSOC 10.0.2.3200
  DPI=600.0975256151248
  KEYWORDS=Topographic, Transportation, Hydrography, Orthoimage, U.S. National Grid, imageryBaseMapsEarthCover, Imagery and Base Maps, Geographic Names Information System
  NEATLINE=POLYGON ((635951.851063535665162 4581600.695550612173975,625386.314008590881713 4581410.198254038579762,625133.892755070934072 4595410.2190180150792,635699.429810017696582 4595600.716314590536058,635951.851063535665162 4581600.695550612173975))
  SUBJECT=This image map depicts geographic features on the surface of the earth.  It was created to provide a representation of accessible geospatial data which is readily available to enhance the capability of Federal, State, and local emergency responders for homeland security efforts.  This image map is generated from selected National Map data holdings and other cartographic data.
  TITLE=USGS 7.5-minute image map for Tionesta, Pennsylvania
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  623457.193, 4596843.381) ( 79d31'13.87"W, 41d30'49.37"N)
Lower Left  (  623775.888, 4579167.631) ( 79d31'13.16"W, 41d21'16.21"N)
Upper Right (  637323.155, 4597093.385) ( 79d21'15.69"W, 41d30'49.35"N)
Lower Right (  637641.850, 4579417.635) ( 79d21'16.44"W, 41d21'16.21"N)
Center      (  630549.521, 4588130.508) ( 79d26'14.79"W, 41d26' 2.90"N)
Band 1 Block=1024x1024 Type=Byte, ColorInterp=Red
  Mask Flags: PER_DATASET ALPHA
Band 2 Block=1024x1024 Type=Byte, ColorInterp=Green
  Mask Flags: PER_DATASET ALPHA
Band 3 Block=1024x1024 Type=Byte, ColorInterp=Blue
  Mask Flags: PER_DATASET ALPHA
Band 4 Block=1024x1024 Type=Byte, ColorInterp=Alpha

And for OH, it gives:
Driver: PDF/Geospatial PDF
Files: OH_Lodi_20100809_TM_geo.pdf
Size is 13869, 17680
Coordinate System is:
PROJCS["UTM Zone 17, Northern Hemisphere",
    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","9108"]],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-81],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]]
GeoTransform =
  403892.2876185254, 0.9998662713660824, -0.01216630733279582
  4554552.367552815, -0.01216630733279582, -0.9998662713660824
Metadata:
  AUTHOR=USGS National Geospatial Technical Operations Center
  CREATION_DATE=D:20100809011028-06'00'
  CREATOR=ESRI ArcSOC 9.2.0.1324
  DPI=609.6432144373267
  KEYWORDS=Topographic, Transportation, Hydrography, Orthoimage, U.S. National Grid, imageryBaseMapsEarthCover, Imagery and Base Maps, Geographic Names Information System
  NEATLINE=POLYGON ((416074.52674049010966 4553140.300224657170475,415905.009541373059619 4539208.831523580476642,405377.212486990378238 4539336.933068925514817,405546.729686107486486 4553268.40177000220865,416074.52674049010966 4553140.300224657170475))
  SUBJECT=This image map depicts geographic features on the surface of the earth.  It was created to provide a representation of accessible geospatial data which is readily available to enhance the capability of Federal, State, and local emergency responders for homeland security efforts.  This image map is generated from selected National Map data holdings and other cartographic data.
  TITLE=USGS 7.5-minute image map for Lodi OH
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  403892.288, 4554552.368) ( 82d 8'42.41"W, 41d 8'11.74"N)
Lower Left  (  403677.187, 4536874.732) ( 82d 8'41.69"W, 40d58'38.49"N)
Upper Right (  417759.433, 4554383.633) ( 81d58'47.61"W, 41d 8'11.75"N)
Lower Right (  417544.333, 4536705.997) ( 81d58'48.32"W, 40d58'38.49"N)
Center      (  410718.310, 4545629.182) ( 82d 3'45.01"W, 41d 3'25.23"N)
Band 1 Block=1024x1024 Type=Byte, ColorInterp=Red
Band 2 Block=1024x1024 Type=Byte, ColorInterp=Green
Band 3 Block=1024x1024 Type=Byte, ColorInterp=Blue

The PA may have issues, but, the OH should not have.

On Wed, Aug 14, 2013 at 9:20 AM, Jukka Rahkonen <jukka.rahkonen at mmmtike.fi<mailto:jukka.rahkonen at mmmtike.fi>> wrote:
Reaves, Timothy <treaves <at> silverfieldstech.com<http://silverfieldstech.com>> writes:

>
>
>
> I'm using GDAL 1.10.0, released 2013/04/24 .  I am playing with TileMill,
and want to add USGS topo maps.
<snip>

> Can anyone help me with this?  I'd appreciate it!

What info does does gdalinfo give with

a) gdalinfo PA_Tionesta_20130621_TM_geo.pdf
b) gdalinfo PA_Tionesta_20130621_TM_geo.tiff

-Jukka Rahkonen-

_______________________________________________
gdal-dev mailing list
gdal-dev at lists.osgeo.org<mailto: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/20130814/4c9003b4/attachment.html>