[gdal-dev] gdalwarp use of Michigan Oblique Mercator
Chaitanya kumar CH
chaitanya.ch at gmail.com
Mon Jan 4 09:46:29 EST 2010
Carolyn,
There seems to be a problem with EPSG:3078.
There was a similar report on this list (
http://lists.osgeo.org/pipermail/gdal-dev/2009-June/020981.html )
where the target and source projections were reversed.
Could you file a ticket at http://trac.osgeo.org/gdal/newticket with
some sample data in EPSG:3078 projection?
On Mon, Jan 4, 2010 at 2:24 PM, Carolyn Gacke <cgacke at usgs.gov> wrote:
>
> 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
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
--
Best regards,
Chaitanya kumar CH.
/tʃaɪθənjə/ /kʊmɑr/
+91-9848167848
17.241582N 80.142635E
More information about the gdal-dev
mailing list