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:[email protected]]
Lähetetty: 14. elokuuta 2013 16:23
Vastaanottaja: Rahkonen Jukka
Kopio: [email protected]
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 
<[email protected]<mailto:[email protected]>> 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
[email protected]<mailto:[email protected]>
http://lists.osgeo.org/mailman/listinfo/gdal-dev

_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to