<div dir="ltr">Cool, thanks. I'll add some print statements to the script to see what's being written, will flip the y-axis coords in the output, and will verify that the origin values and pixel dimensions are properly expressed in the output.<br>
<br>It's neat that you can do this sort of thing with gdal.<br>--<br><br><div class="gmail_quote">On Sun, Aug 10, 2008 at 10:21 AM, Jamie Adams <span dir="ltr"><<a href="mailto:jaadfoo@gmail.com">jaadfoo@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div dir="ltr">It looks like Matt's script doesn't modify the geotransform, just reverses the rows in the image. You probably need to calculate a new ulx and make sure the pixel dimensions are expressed correctly (+/-).<br>
<br>From the sound of it, your original image would actually be right if the viewer software read the corner coordinates and displayed the image based on those. :-)<br><br>-Jamie<br><br><div class="gmail_quote"><div><div>
</div><div class="Wj3C7c">On Fri, Aug 8, 2008 at 4:51 PM, Roger André <span dir="ltr"><<a href="mailto:randre@gmail.com" target="_blank">randre@gmail.com</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div><div></div><div class="Wj3C7c"><div dir="ltr">My apologies to Matt, the flip_raster.py does flip something, as when I compare the images in OpenEV, I can see that they are reversed. But when I compare their corner coords with gdalinfo, they appear identical.<br>
<br>$ gdalinfo test.tif<br><snip><br>Corner Coordinates:<br>Upper Left ( -71.3600000, -23.5000000)<br>Lower Left ( -71.3600000, -8.5000000)<br>Upper Right ( -55.7000000, -23.5000000)<br>Lower Right ( -55.7000000, -8.5000000)<br>
<br>$ gdalinfo bolivia_fullres.tif<br><snip><br>Corner Coordinates:<br>Upper Left ( -71.3600000, -23.5000000)<br>Lower Left ( -71.3600000, -8.5000000)<br>Upper Right ( -55.7000000, -23.5000000)<br>Lower Right ( -55.7000000, -8.5000000)<br>
<br>My head hurts.<br><font color="#888888">--</font><div><div></div><div><br><br><div class="gmail_quote">On Fri, Aug 8, 2008 at 4:24 PM, Roger André <span dir="ltr"><<a href="mailto:randre@gmail.com" target="_blank">randre@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div dir="ltr">Ok, this is driving me nuts, I can't handle "South is up" anymore. I have NetCDF derived GMT grids which have a flipped y-axis. This makes for georeferencing that looks like this:<br>Corner Coordinates:<br>
Upper Left ( -90.0000000, -90.0000000)<br>Lower Left ( -90.0000000, 0.0000000)<br>Upper Right ( 0.0000000, -90.0000000)<br>Lower Right ( 0.0000000, 0.0000000)<br><br>I've seen some old discussions about this, but haven't really found a good solution. I've tried a "flip_raster.py" utility that Matt Perry wrote, but it doesn't seem to do anything with my data. (probably due to my crappy data, not Matt's tool) I've also read that there have been some rewrites of the NetCDF driver that are supposed to be able to deal with this, but I'm not sure about that.<br>
<br>How have others out there dealt with this issue? I'd like to be able to convert the grids into GeoTIFF's and have the TIFF's be "right" side up.<br><br>Thanks,<br><br>Flipped-off Roger<br>--<br>
</div>
</blockquote></div><br></div></div></div>
<br></div></div>_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br></blockquote></div><br></div>
</blockquote></div><br></div>