<html>
<body>
<blockquote type=cite class=cite cite>Rich,<br><br>
I appologise for not responding to your earlier message yet. I have
been<br>
bogged down in deliverables, and an addiction to an MMORPG
(<a href="http://www.atitd.com/" eudora="autourl">www.atitd.com</a>).<br><br>
I tried something similar to your earlier example with a UTM 11 file
I<br>
have here, and it seemed to work pretty well. I used this
command:<br><br>
gdalwarp<br>
-s_srs "+proj=utm +ellps=GRS80 +zone=11 +north +datum=NAD83
+units=m +no_defs"<br>
-t_srs "+proj=tmerc +ellps=GRS80 +units=ft
+lon_0=-110.0833333333 +lat_0=40.5<br>
+x_0=800000 +y_0=100000 +k=0.9999375
+no_defs"<br>
-order 3 openev/utm.tif nw.tif<br><br>
Unlike your experience, my output did seem to get the right linear
units:<br><br>
Coordinate System is:<br>
PROJCS["unnamed",<br>
GEOGCS["unnamed",<br>
DATUM["unknown",<br>
SPHEROID["unnamed",6378137,298.2572221010042]],<br>
PRIMEM["Greenwich",0],<br>
UNIT[,0.0174532925199433]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",40.5],<br>
PARAMETER["central_meridian",-110.0833333333],<br>
PARAMETER["scale_factor",0.9999375],<br>
PARAMETER["false_easting",2624671.916010499],<br>
PARAMETER["false_northing",328083.9895013123],<br>
UNIT["foot",0.3048,<br>
AUTHORITY["EPSG","9002"]]]<br><br>
Based on that, I now wonder what is different about your build and
mine.<br>
Are you building from a recent CVS update? I suggest we try to
isolate<br>
this issue off the list today if you have the time.</blockquote><br>
Frank, <br><br>
No apologies are necessary on your part, I apologize for imposing on your
time. You suggested that we go off list, but I do not seem to have
an "off line" address for you.<br><br>
I built from cvs on sunday, win 2k, VC 6.0.<br><br>
Here are three examples of how I am launching gdalwarp, the resulting
gdalinfo, and a note on the resulting image (rotated or not
rotated)<br><br>
=== example 1 ===<br><br>
gdalwarp <br>
-s_srs "+proj=utm +ellps=GRS80 +zone=12 +north +datum=NAD83
+units=m"<br>
-t_srs "+proj=tmerc +lat_0=40.5 +lon_0=-110.0833333333333
+k=0.999938 +x_0=800000 +y_0=100000 +ellps=GRS80 +datum=NAD83
+to_meter=0.3048006096012192" <br>
nw.tif <br>
nw_spsc.tif<br><br>
gdalinfo<br>
Driver: GTiff/GeoTIFF<br>
Size is 5764, 7625<br>
Coordinate System is:<br>
PROJCS["unnamed",<br>
GEOGCS["NAD83",<br>
DATUM["North_American_Datum_1983",<br>
SPHEROID["GRS 1980",6378137,298.2572221010002,<br>
AUTHORITY["EPSG","7019"]],<br>
AUTHORITY["EPSG","6269"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4269"]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",40.5],<br>
PARAMETER["central_meridian",-110.0833333333333],<br>
PARAMETER["scale_factor",0.999938],<br>
PARAMETER["false_easting",800000],<br>
PARAMETER["false_northing",100000],<br>
UNIT["unknown",1]]<br>
Origin = (590671.941324,1423019.616154)<br>
Pixel Size = (3.282086,-3.282086)<br>
Corner Coordinates:<br>
Upper Left ( 590671.941, 1423019.616) (113d 9'24.33"W,
52d21'45.75"N)<br>
Lower Left ( 590671.941, 1397993.709) (113d 8'28.54"W,
52d 8'17.19"N)<br>
Upper Right ( 609589.886, 1423019.616) (112d52'45.75"W,
52d22'10.56"N)<br>
Lower Right ( 609589.886, 1397993.709) (112d51'54.98"W, 52d
8'41.81"N)<br>
Center ( 600130.914, 1410506.663)
(113d 0'38.36"W, 52d15'14.12"N)<br>
Band 1 Block=5764x1 Type=Byte, ColorInterp=Gray<br><br>
<b>(image is rotated, lon and lat are totally off)<br><br>
</b>=== example 2 ===<br>
(only difference from example 1 is unit definition in t_srs)<br><br>
gdalwarp <br>
-s_srs "+proj=utm +ellps=GRS80 +zone=12 +north +datum=NAD83
+units=m" <br>
-t_srs "+proj=tmerc +lat_0=40.5 +lon_0=-110.0833333333333
+k=0.999938 +x_0=800000 +y_0=100000 +ellps=GRS80 +datum=NAD83
+units=ft" <br>
nw.tif <br>
nw_spsc.tif<br><br>
gdalinfo<br>
Driver: GTiff/GeoTIFF<br>
Size is 5764, 7625<br>
Coordinate System is:<br>
PROJCS["unnamed",<br>
GEOGCS["NAD83",<br>
DATUM["North_American_Datum_1983",<br>
SPHEROID["GRS 1980",6378137,298.2572221010002,<br>
AUTHORITY["EPSG","7019"]],<br>
AUTHORITY["EPSG","6269"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4269"]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",40.5],<br>
PARAMETER["central_meridian",-110.0833333333333],<br>
PARAMETER["scale_factor",0.999938],<br>
PARAMETER["false_easting",800000],<br>
PARAMETER["false_northing",100000],<br>
UNIT["unknown",1]]<br>
Origin = (590671.522667,1423022.262198)<br>
Pixel Size = (3.282093,-3.282093)<br>
Corner Coordinates:<br>
Upper Left ( 590671.523, 1423022.262) (113d 9'24.36"W,
52d21'45.83"N)<br>
Lower Left ( 590671.523, 1397996.305) (113d 8'28.57"W,
52d 8'17.27"N)<br>
Upper Right ( 609589.505, 1423022.262) (112d52'45.77"W,
52d22'10.65"N)<br>
Lower Right ( 609589.505, 1397996.305) (112d51'55.00"W, 52d
8'41.89"N)<br>
Center ( 600130.514, 1410509.284)
(113d 0'38.39"W, 52d15'14.21"N)<br>
Band 1 Block=5764x1 Type=Byte, ColorInterp=Gray<br><br>
<b>(image is rotated, lon and lat are totally off)<br><br>
</b>=== example 3 ===<br><br>
gdalwarp <br>
-s_srs "+proj=utm +ellps=GRS80 +zone=12 +north +datum=NAD83
+units=m" <br>
-t_srs WY83-WF.prj <b>(see below for my prj file)<br>
</b>nw.tif <br>
nw_prj.tif<br><br>
gdalinfo<br>
Driver: GTiff/GeoTIFF<br>
Size is 5680, 7562<br>
Coordinate System is:<br>
PROJCS["unnamed",<br>
GEOGCS["NAD83",<br>
DATUM["North_American_Datum_1983",<br>
SPHEROID["GRS 1980",6378137,298.2572221010002,<br>
AUTHORITY["EPSG","7019"]],<br>
AUTHORITY["EPSG","6269"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4269"]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",0],<br>
PARAMETER["central_meridian",-111],<br>
PARAMETER["scale_factor",0.9996],<br>
PARAMETER["false_easting",500000],<br>
PARAMETER["false_northing",0],<br>
UNIT["metre",1,<br>
AUTHORITY["EPSG","9001"]],<br>
AUTHORITY["EPSG","26912"]]<br>
Origin = (509637.500000,4886074.500000)<br>
Pixel Size = (1.000000,-1.000000)<br>
Corner Coordinates:<br>
Upper Left ( 509637.500, 4886074.500) (110d52'46.32"W,
44d 7'40.08"N)<br>
Lower Left ( 509637.500, 4878512.500) (110d52'46.82"W,
44d 3'34.98"N)<br>
Upper Right ( 515317.500, 4886074.500) (110d48'30.73"W, 44d
7'39.73"N)<br>
Lower Right ( 515317.500, 4878512.500) (110d48'31.52"W, 44d
3'34.63"N)<br>
Center ( 512477.500, 4882293.500)
(110d50'38.85"W, 44d 5'37.38"N)<br>
Band 1 Block=5680x1 Type=Byte, ColorInterp=Gray<br><br>
<b>image is NOT rotated, <br>
X and Y are the UTM input origin +/- 0.5 m<br><br>
</b>My WY83-WF.prj file:<br>
PROJCS[<br>
<x-tab> </x-tab>"NAD_1983_StatePlane_Wyoming_West_FIPS_4904_Feet",<br>
<x-tab> </x-tab>GEOGCS["GCS_North_American_1983",<br>
<x-tab> </x-tab><x-tab> </x-tab>DATUM["D_North_American_1983",<br>
<x-tab> </x-tab><x-tab> </x-tab><x-tab> </x-tab>SPHEROID["GRS_1980",6378137,298.257222101]<br>
<x-tab> </x-tab><x-tab> </x-tab>],<br>
<x-tab> </x-tab><x-tab> </x-tab>PRIMEM["Greenwich",0],<br>
<x-tab> </x-tab><x-tab> </x-tab>UNIT["Degree",0.017453292519943295]<br>
<x-tab> </x-tab>],<br>
<x-tab> </x-tab>PROJECTION["Transverse_Mercator"],<br>
<x-tab> </x-tab>PARAMETER["False_Easting",2624666.666666666],<br>
<x-tab> </x-tab>PARAMETER["False_Northing",328083.3333333333],<br>
<x-tab> </x-tab>PARAMETER["Central_Meridian",-110.0833333333333],<br>
<x-tab> </x-tab>PARAMETER["Scale_Factor",0.9999375],<br>
<x-tab> </x-tab>PARAMETER["Latitude_Of_Origin",40.5],<br>
<x-tab> </x-tab>UNIT["Foot_US",0.30480060960121924]<br>
]<br>
<x-sigsep><p></x-sigsep>
Richard W. Greenwood, PLS<br>
Greenwood Mapping, Inc.<br>
Rich@GreenwoodMap.com<br>
(307) 733-0203<br>
<a href="http://www.greenwoodmap.com/" eudora="autourl">http://www.GreenwoodMap.com</a></body>
</html>