<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br>
<br><div><div>On Feb 11, 2011, at 12:14 PM, Stuart Edwards wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">I'm not familiar with the Kyngesburye plugin and haven't used gdalwarp, but here:<div><br></div><div><blockquote type="cite"><div><blockquote type="cite"><blockquote type="cite">I want to import it into a UTM zone 17 project. So I did the following:<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">gdalwarp -t_srs '+proj=utm +zone=17 +datum=WGS84' input.sid output_utm.sid</blockquote></blockquote></div></blockquote><div><br></div>it looks like you are just trying to create a new .sid file in zone 17 (output_utm.sid). From the meta file at the MI source, I think the file is already projected in utm - probably in zone 17, and with units of 'meters'. So what you maybe intended was to just create a .tiffg from the .sid so that it can be imported into GRASS.</div><div><br></div><div>In mrsiddecode (the app that does all the work inside the plugin - I assume) I would just issue the command</div><div><br></div><div><div>./mrsiddecode -i Hancock.sid -o Hancock_south.tif -of tifg -ulxy 0 0 -lrxy 71145 29160 </div></div><div><br></div><div>where</div><div><br></div><div>-i - input</div><div>-o - output</div><div>-of - output file type</div><div>ulxy - upper left x,y</div><div>lrxy - lower right x,y</div><div><br></div><div>and you get a nice geotiff that r.in.gdal will recognize - already in utm zone 17 most likely. Actually, in my case, the OH mrsids are projected in state plane coordinates but you get the picture.</div><div><br></div><div>To get the mrsiddecode app you must register at </div><div><br></div><div><span class="Apple-style-span" style="font-family: verdana; font-size: 11px; color: rgb(102, 102, 102); line-height: 15px; "><span style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; padding-top: 0px; padding-right: 0px; padding-bottom: 0px; padding-left: 0px; font-weight: bold; "><a href="http://developer.lizardtech.com/">http://developer.lizardtech.com/</a></span> </span></div><div><span class="Apple-style-span" style="font-family: verdana; font-size: 11px; color: rgb(102, 102, 102); line-height: 15px; "><br></span></div><div>so you can access to the Mac version of the decoder and mrsidinfo download page. And be slightly comfortable on the command line.</div><div><br></div><div>Stu</div><div><br></div><div><br></div></div></blockquote></div><br><div>Nope, that's not doing it. I downloaded mrsiddecode and ran it. It nicely generates a geotiff. But, if I try to import that I still get an error because the project is UTM, while the image is Oblique Mercator with omerc projection. I tried to use gdalwarp to go from the geotiff to utm, but I get an error:</div><div><br></div><div><div> gdalwarp -t_srs '+proj=utm +zone=17 +datum=WGS84' oak_grove_sw.tifg oak_grove_sw_utm.tifg</div><div>ERROR 1: latitude or longitude exceeded limits</div><div>ERROR 1: Too many points (441 out of 441) failed to transform,</div><div>unable to compute output bounds.</div></div><div><br></div><div>If I run gdalinfo on the geotiff I get it looks to me like the projection has lat long near 0. So, I don't think that the problem is a MrSID problem. But instead there is just something wrong with how I am doing the projection, or something is wrong with the file itself.</div><div><br></div><div>Here is the output of gdalinfo oak_grove_sw.tifg:</div><div><br></div><div><br></div><div><div>Driver: GTiff/GeoTIFF</div><div>Files: oak_grove_sw.tifg</div><div>Size is 6415, 8061</div><div>Coordinate System is:</div><div>PROJCS["Projection = Oblique Mercator (Hotine)",</div><div> GEOGCS["NAD83",</div><div> DATUM["North_American_Datum_1983",</div><div> SPHEROID["GRS 1980",6378137,298.2572221010002,</div><div> AUTHORITY["EPSG","7019"]],</div><div> AUTHORITY["EPSG","6269"]],</div><div> PRIMEM["Greenwich",0],</div><div> UNIT["degree",0.0174532925199433],</div><div> AUTHORITY["EPSG","4269"]],</div><div> PROJECTION["Hotine_Oblique_Mercator"],</div><div> PARAMETER["latitude_of_center",0.9996],</div><div> PARAMETER["longitude_of_center",0.9996],</div><div> PARAMETER["azimuth",0.9996],</div><div> PARAMETER["rectified_grid_angle",90],</div><div> PARAMETER["scale_factor",0.9996],</div><div> PARAMETER["false_easting",0.9996],</div><div> PARAMETER["false_northing",0.9996],</div><div> UNIT["metre",1,</div><div> AUTHORITY["EPSG","9001"]]]</div><div>Origin = (663144.500000000000000,240035.500000000000000)</div><div>Pixel Size = (1.000000000000000,-1.000000000000000)</div><div>Metadata:</div><div> AREA_OR_POINT=Area</div><div>Image Structure Metadata:</div><div> INTERLEAVE=PIXEL</div><div>Corner Coordinates:</div><div>Upper Left ( 663144.500, 240035.500) ( 6d57'27.03"E, 3d 9'14.09"N)</div><div>Lower Left ( 663144.500, 231974.500) ( 6d57'25.53"E, 3d 4'52.96"N)</div><div>Upper Right ( 669559.500, 240035.500) ( 7d 0'53.76"E, 3d 9'12.88"N)</div><div>Lower Right ( 669559.500, 231974.500) ( 7d 0'52.24"E, 3d 4'51.77"N)</div><div>Center ( 666352.000, 236005.000) ( 6d59'9.64"E, 3d 7'2.93"N)</div><div>Band 1 Block=6415x1 Type=Byte, ColorInterp=Red</div><div>Band 2 Block=6415x1 Type=Byte, ColorInterp=Green</div><div>Band 3 Block=6415x1 Type=Byte, ColorInterp=Blue</div></div><div><br></div></body></html>