<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
<br>
<br>
Jukka Rahkonen wrote:
<blockquote cite="mid:loom.20090208T163124-633@post.gmane.org"
type="cite">
<pre wrap="">Jan Hartmann <j.l.h.hartmann <at> uva.nl> writes:
</pre>
<blockquote type="cite">
<pre wrap="">No, this doesn't seem to work. Although the coordinate system *is*
retained when I add -a_srs to the gcps, the boundaries of the input file
are discarded: they are reset to the pixel dimensions of the file. Is
this a bug?
Jan
Frank Warmerdam wrote:
</pre>
<blockquote type="cite">
<pre wrap="">Jan Hartmann wrote:
</pre>
<blockquote type="cite">
<pre wrap="">Hi,
I am trying to georeference old maps in two steps: first
georeferencing them from the old projection to a modern one, and then
ameliorating the georeferenced image by rubber-sheeting it on the
base of control points, mostly triangulation points for which the
coordinates in the old and new maps are known.
The second step causes trouble: if I add control points to the
georeferenced image with gdal_translate -gcp , the projection
information gets lost.
</pre>
</blockquote>
</blockquote>
</blockquote>
<pre wrap=""><!---->
Hi,
I am not sure if I understood how you trying to do the job, but hopefully my
task and work flow are close enough to be useful for you.
1. Let's say that the original image is in epsg:2393 and target is epsg:3067.
2. I insert some ground control points in epsg:3067 with gdal_translate and
-gcp. Image format of the temporal output file can be GeoTIFF but also virtual
raster(vrt).
3. I warp the temporary image, now containing gpc's in epsg:3067 projection, to
final product using gdalwarp as
gdalwarp -s_srs epsg:3067 -t_srs epsg:3067.
-Jukka Rahkonen-
_______________________________________________
gdal-dev mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="http://lists.osgeo.org/mailman/listinfo/gdal-dev">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a>
</pre>
</blockquote>
Hi Jukka,<br>
<br>
You have described exactly what I want to do, but the problem is that
in step 2, inserting gcps, the georeferenced boundaries of the file are
lost. Eventually you get a somewhat displaced raw raster file.<br>
<br>
Jan <br>
</body>
</html>