<div dir="ltr">Hi Jukka,<div><br></div><div>I finally found the time to produce a test set.</div><div><br></div><div>You can download it here: <a href="https://drive.google.com/file/d/1Y08Q-tIm5dxyKFKNdVqilvAO3H7FFFbx/view?usp=sharing" target="_blank">https://drive.google.com/file/d/1Y08Q-tIm5dxyKFKNdVqilvAO3H7FFFbx/view?usp=sharing</a></div><div><br></div><div>I started from an unreferenced raster file (raster2tps.tif) with gcp's (gcp4tps.gcp) and transformed it with tps (tpsraster.tif).</div><div>Then polygonized the unreferenced raster file (vector2tps.shp) and transformed the result with the same gcp's (gcp4tps.gcp) and with tps (tpsvector.shp).</div><div><br></div><div>The vector2tps.shp polygons are "flipped" because of the different Y-origin for rasters and vectors but this way both datasets can use the exact same gcp's. </div><div><br></div><div>When you lay the tpsvector-result on top of the tpsraster-result (in QGis for instance) you'll see the differences in how both are transformed.</div><div><br></div><div>Kind regards,</div><div><br></div><div>Stijn</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Op wo 16 aug 2023 om 13:16 schreef Stijn Tallir <<a href="mailto:stijn%2Bgdal-dev@strict.be">stijn+gdal-dev@strict.be</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Yes, I checked them visually for both raster and vector.<div><br></div><div>I compared the results also visually. The rasters are transformed in a way that the end ponts of the gcp's align exactly with the result so that is why I referred to it as "right". The vector data result is in the neighbourhood of the end points (sometimes a rather significant distance).</div><div><br></div><div>The result is different from order 1-3 transformations so I presume the tps option isn't ignored.</div><div><br></div><div>Stijn</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Op wo 16 aug 2023 om 11:52 schreef Rahkonen Jukka <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi" target="_blank">jukka.rahkonen@maanmittauslaitos.fi</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>
<div lang="FI">
<div>
<p class="MsoNormal"><span>Hi,<u></u><u></u></span></p>
<p class="MsoNormal"><span><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Did you check the ground control points? What is your reference when you say that one result is right, and another wrong? Have you used some other software for comparison? Or do you
only know that the results are different?<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">-Jukka-<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(225,225,225);padding:3pt 0cm 0cm">
<p class="MsoNormal"><b>Lähettäjä:</b> Stijn Tallir <<a href="mailto:stijn@strict.be" target="_blank">stijn@strict.be</a>> <br>
<b>Lähetetty:</b> keskiviikko 16. elokuuta 2023 12.27<br>
<b>Vastaanottaja:</b> Rahkonen Jukka <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi" target="_blank">jukka.rahkonen@maanmittauslaitos.fi</a>><br>
<b>Kopio:</b> <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<b>Aihe:</b> Re: [gdal-dev] tps - gdalwarp vs ogr2ogr<u></u><u></u></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">Hi Jukka,<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I thought of the density as an option for the "error" as you suggested and I made a point-file with a point for every pixel in my original image and used this as a source for the ogr2ogr transformation. So you could say the desnity for
both sources raster and vector) are then alike.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">The results were still the same (and wrong) ... <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Stijn<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<p class="MsoNormal">Op wo 16 aug 2023 om 10:22 schreef Rahkonen Jukka <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi" target="_blank">jukka.rahkonen@maanmittauslaitos.fi</a>>:<u></u><u></u></p>
</div>
<blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin-left:4.8pt;margin-right:0cm">
<div>
<div>
<div>
<p class="MsoNormal">Hi,<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US">Without test data it is very hard to say much. I believe that the promise of tps is that the ground control points stay where they are set. The intermediate points
follow the least tension surfaces and I do not know how exactly those spline algorithms are defined. Raster data is full of points to warp but probably in the vector data the transformation is done vertex by vertex. I would first check if the GCPs are in the
same place in both outputs. Then I would make a test by densifying the vector data before georeferencing to have much more vertices and see if it has any effect on the result.</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US">-Jukka Rahkonen-</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(225,225,225);padding:3pt 0cm 0cm">
<p class="MsoNormal"><b>Lähettäjä:</b> gdal-dev <<a href="mailto:gdal-dev-bounces@lists.osgeo.org" target="_blank">gdal-dev-bounces@lists.osgeo.org</a>>
<b>Puolesta </b>Stijn Tallir<br>
<b>Lähetetty:</b> keskiviikko 16. elokuuta 2023 10.29<br>
<b>Vastaanottaja:</b> <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<b>Aihe:</b> [gdal-dev] tps - gdalwarp vs ogr2ogr<u></u><u></u></p>
</div>
<p class="MsoNormal"> <u></u><u></u></p>
<div>
<p class="MsoNormal">Hi,<u></u><u></u></p>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">According to the documentation gdal and ogr use the same algorithm for the tps-transformation but I don't seem to get the same results using the same set of gcp's for images and
vectors.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">I have images that are unreferenced and vector data digitised on these images (in pixel coordinates).<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">The images are then georeferenced with +100 gcp's and warped with gdalwarp using the "tps" option.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">When I use the same gcp's (with adjusted y-origin to lower left corner) to georeference the vector data with ogr2ogr and the "tps" option I get different results. The vector-result
is similar to the image-result but never exactly the same and differences can be substantial.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Any thoughts?<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Stijn<u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<p class="MsoNormal"><br clear="all">
<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal"><span>-- </span><u></u><u></u></p>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">Stijn Tallir - StrICT BV<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Wijnveld 8<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">9112 Sinaai-Waas<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">GSM: 0486 750220<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">E-mail: <a href="mailto:Info@strict.be" target="_blank">Info@strict.be</a><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Web: <a href="http://www.strict.be/" target="_blank">www.strict.be</a><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">BTW: BE 0567.559.668<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><img border="0" width="289" height="162" style="width: 3.0104in; height: 1.6875in;" id="m_5154515844185920288m_8416825920031188548_x0000_i1025" src="http://www.strict.be/strict_web2.jpg"><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div></blockquote></div>
</blockquote></div>