[gdal-dev] gdalwarp use of Michigan Oblique Mercator

Carolyn Gacke cgacke at usgs.gov
Mon Jan 4 09:24:29 EST 2010


Hello-
Is anyone aware of the use of gdalwarp on GeoTIFF files; projecting from 
Michigan Oblique Mercator  (epsg 3078) to WGS84 (epsg 4326)?  ls
EXAMPLE: gdalwarp -t_srs EPSG:4326 085_r3c25.tif 085_r3c25warp.tif

We end up NE of Alaska, where as the original tifs overlay just fine in 
Michigan when using ArcMap 9.3.1 

These are tifs of high resolution orhto flown in 2006 over Lake County MI; 
the MI projection info appears correct in the tags.

We have quite a bit of high res ortho datasets; at times we have to modify 
tfws/tags, however this one has us stumped(!)

Here is the gdalinfo of the original GTIFF;  it also interprets the xys to 
lat/longs in Alaska:

TIF: gdalinfo 085_r3c25.tif
Driver: GTiff/GeoTIFF
Size is 1619, 1708
Coordinate System is:
PROJCS[" Projection Name = Oblique Mercator (Hotine) Units = meters 
GeoTIFF Units = meters",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.2572221010042,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Hotine_Oblique_Mercator"],
    PARAMETER["latitude_of_center",45.3091667],
    PARAMETER["longitude_of_center",-86],
    PARAMETER["azimuth",337.25556],
    PARAMETER["rectified_grid_angle",90],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",2546731.496],
    PARAMETER["false_northing",-4354009.816],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (519086.028692,400863.870183)
Pixel Size = (0.60960000,-0.60960000)
Metadata:
  AREA_OR_POINT=Area
Corner Coordinates:
Upper Left  (  519086.029,  400863.870) (173d36'53.89"W, 73d28'19.77"N)
Lower Left  (  519086.029,  399822.673) (173d34'57.28"W, 73d28'14.49"N)
Upper Right (  520072.971,  400863.870) (173d36'36.33"W, 73d28'51.24"N)
Lower Right (  520072.971,  399822.673) (173d34'39.66"W, 73d28'45.96"N)
Center      (  519579.500,  400343.272) (173d35'46.79"W, 73d28'32.87"N)
Band 1 Block=1619x5 Type=Byte, ColorInterp=Gray

TFW:  cgacke at edclxs84> more 085_r3c25.tfw
0.609600000000000030
0.0
0.0
-0.609600000000000030
519086.33349170600
400863.56538265612


Any advice or suggestions would be greatly appreciated!

Carolyn Gacke
SGT Inc.
Contractor to the USGS EROS
Sioux Falls SD
Email:     cgacke at usgs.gov
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20100104/25d68b63/attachment.html


More information about the gdal-dev mailing list