Thanks, Klokan. I'll check this out. <br><br><div class="gmail_quote">On Wed, Apr 1, 2009 at 2:52 PM, Klokan Petr Přidal <span dir="ltr"><<a href="mailto:klokan@klokan.cz">klokan@klokan.cz</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;">
Hi,<br>
<br>
I am not familiar with OziExplorer Map files, but I have seen a<br>
project trying to implement what you need, maybe you can reuse it (and<br>
probably improve it and add support of different projections):<br>
<br>
<a href="http://ozi2geotiff.svn.sourceforge.net/viewvc/ozi2geotiff/trunk/ozi2geotiff.py?view=markup" target="_blank">http://ozi2geotiff.svn.sourceforge.net/viewvc/ozi2geotiff/trunk/ozi2geotiff.py?view=markup</a><br>
<br>
I guess you would like to have affine transformation georeference from<br>
your GCPs points in your example with gdal_translate.<br>
Then you should either warp your raster with gdalwarp utility or use a<br>
tool like gcps2wld.py to calculate a best fitting World file and keep<br>
raster as it is.<br>
<br>
Klokan<br>
<br>
P.S.<br>
A standalone application for World file generation then you can use<br>
independently on GDAL is:<br>
<a href="http://blog.oldmapsonline.org/2008/06/georeferencing-images-by-control-points.html" target="_blank">http://blog.oldmapsonline.org/2008/06/georeferencing-images-by-control-points.html</a><br>
or very simple python script doing the calculation as well:<br>
<a href="http://code.google.com/p/oldmapsonline/source/browse/trunk/gcps2wld/python-no-dep/gcps2wld.py" target="_blank">http://code.google.com/p/oldmapsonline/source/browse/trunk/gcps2wld/python-no-dep/gcps2wld.py</a><br>
<br>
2009/3/30 erik <<a href="mailto:whatevar89@gmail.com">whatevar89@gmail.com</a>>:<br>
<div><div></div><div class="h5">> Hello,<br>
><br>
> I'm trying to convert a map from OziExplorer Map file fromat to Geotiff.<br>
> FOr this i'm using the calibration points in the Map file as Ground Control<br>
> Points.<br>
><br>
> I do the following:<br>
><br>
> gdal_translate \<br>
> -gcp 203 52 529000 4191000 \<br>
> -gcp 5174 3096 544000 4182000 \<br>
> -gcp 5201 99 544000 4191000 \<br>
> -gcp 176 3051 529000 4182000 \<br>
> -gcp 1861 1067 534000 4188000 \<br>
> -gcp 3517 2081 539000 4185000 \<br>
> -gcp 3526 1083 539000 4188000 \<br>
> -gcp 1851 2066 534000 4185000 \<br>
> -gcp 2844 2742 537000 4183000 \<br>
> -a_srs "+proj=utm +ellps=intl +k=0.999600 +x_0=500000.0 +y_0=0.0<br>
> +zone=30 +units=m<br>
> +towgs84=-125.098545,-76.000054,-156.198703,0.0,0.0,-1.129,8.30463103<br>
> +no_defs" \<br>
> map.bmp \<br>
> map.tif<br>
><br>
> However doing this still gives the maps upper left corner start at 0-0.<br>
> I was wondering if there was a way to make gdal infer the correct corners of<br>
> the map from the gcp's without specifying them directly.<br>
> Or is there another way to achieve this without specifying the coord<br>
> corners? (The calibration points are not corner-based)<br>
> One solution would be to calculate the corners myself programmatically (from<br>
> the pixels), but this would be the least good solution.<br>
><br>
> Many thanks,<br>
> Erik<br>
><br>
</div></div>> _______________________________________________<br>
> gdal-dev mailing list<br>
> <a href="mailto:gdal-dev@lists.osgeo.org">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>
><br>
<font color="#888888"><br>
<br>
<br>
--<br>
<a href="http://blog.klokan.cz/" target="_blank">http://blog.klokan.cz/</a><br>
<a href="http://www.maptiler.org/" target="_blank">http://www.maptiler.org/</a><br>
<a href="http://www.oldmapsonline.org/" target="_blank">http://www.oldmapsonline.org/</a><br>
</font></blockquote></div><br>