<div dir="ltr"><div>BINGO !<br></div><div><br></div><div>Il faut retourner les images de 180° lorsque orientation + ST_Azimuth(du premier au 3ème point du polygone d'emprise) sont proches de 180 et pas de 360...<br><br></div><div>En postgis ça donne un truc comme:<br><br>(orientation::numeric + round(degrees(st_azimuth(st_pointn(st_exteriorring(p.geom),1)::geography,st_pointn(st_exteriorring(p.geom),3)::geography))::numeric,0))%360 < 270 as inverser<br><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">Le 6 mars 2016 à 20:32, Christian Quest <span dir="ltr"><<a href="mailto:cquest@openstreetmap.fr" target="_blank">cquest@openstreetmap.fr</a>></span> a écrit :<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>J'ai récupéré les photos de 1921 prises au dessus de ma commune, il y en a une vingtaine.<br><br></div>En gros c'est plutôt correct sauf que pour 3 clichés on est à 180° :(<br><br>En rouge sur la copie d'écran QGis: <a href="https://framapic.org/MMcURbVTytm9/iZoQxP0IIY5Y.png" target="_blank">https://framapic.org/MMcURbVTytm9/iZoQxP0IIY5Y.png</a><br><br></div><div>Malheureusement, entre le polygone d'emprise et l'orientation, je n'ai pas trouvé de critère pour les détecter.<br><br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">Le 6 mars 2016 à 18:17, Christian Quest <span dir="ltr"><<a href="mailto:cquest@openstreetmap.fr" target="_blank">cquest@openstreetmap.fr</a>></span> a écrit :<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Pour moi c'est ok en procédant ainsi:<br><br>{<br> "type": "FeatureCollection",<br> "features": [<span><br> {"type":"Feature",<br> "geometry":{<br> "type":"Polygon",<br></span> "coordinates":[[[2.5126989,48.8199908],<br> [2.5213251,48.8149892],<br> [2.5298119,48.8101662],<br> [2.5203523,48.8023041],<br> [2.5102159,48.7943484],<br> [2.5013197,48.7993482],<br> [2.4924227,48.8042573],<br> [2.5025584,48.8122145],<br> [2.5126989,48.8199908]]]},<br> "properties":{"idcliche":"CCF00B-361_1921_CAF_B-36_0144","mission":"CF00B-361","numcli":"144","idta":"1921_CAF_B-36_P_10000","date":"1921-06-11","res":"10010","support":"Ag","type":"P","surface":3.87479004551169,"lon":2.51125156305732,"lat":48.8071519050115,"orientation":"229","url":"IGNF_PVA_1-0__1921-06-11__CCF00B-361_1921_CAF_B-36_0144.jp2"}}]<br>}<br><br><br>w=6774<br>h=5050<br>gdal_translate -of GTiff IGNF_PVA_1-0__1921-06-11__CCF00B-361_1921_CAF_B-36_0144.jp2 temp.tif \<br> -gcp $w 0 2.5126989 48.8199908 \<br> -gcp $w $h 2.5298119 48.8101662 \<br> -gcp 0 $h 2.5102159 48.7943484 \<br> -gcp 0 0 2.4924227 48.8042573<br>gdalwarp -s_srs EPSG:4326 -r cubic -co COMPRESS=LZW -dstalpha -overwrite temp.tif IGNF_PVA_1-0__1921-06-11__CCF00B-361_1921_CAF_B-36_0144.geotiff<br><br></div>Il y a 9 coordonnées, les coins sont en positions impaires (en comptant de 1) et ça démarre en haut à droite...<br><br><br></div><div>Le calage n'est pas nickel, mais c'est une bonne première approximation !<br><br></div><div>Exemple ici: <a href="https://framapic.org/dLIvBVa9WTJz/RXxKlfWqAwrk.png" target="_blank">https://framapic.org/dLIvBVa9WTJz/RXxKlfWqAwrk.png</a><br><br></div><div><div><div><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>Le 6 mars 2016 à 18:02, <span dir="ltr"><<a href="mailto:osm.sanspourriel@spamgourmet.com" target="_blank">osm.sanspourriel@spamgourmet.com</a>></span> a écrit :<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>
<div bgcolor="transparent" text="#000000">
On y arrive, je vais laisser les Normands remettre la photo à
l'endroit.<br>
J'ai utilisé Qgis en affichant les coordonnées des points de
contrôle, càd étiquettes : name || ' (' || $x || ' ' || $y ||
')'.<br>
<img src="cid:part1.03070705.01080500@gmx.net" alt=""><br>
<br>
Ce qui équivaut à :<br>
gdal_translate -of GTiff -gcp 0 0 0.154072 49.4906 -gcp 6804 0
0.162634 49.4931 -gcp 6804 5090 0.167743 49.4856 -gcp 0 5090
0.159181 49.4831
"/Temp/IGNF_PVA_1-0__1919__CCF000C-71_1919_CAF_C-7_0006.jp2"
"C:/Users/landrac/AppData/Local/Temp/IGNF_PVA_1-0__1919__CCF000C-71_1919_CAF_C-7_0006.jp2"
<p style="margin-top:0px;margin-bottom:0px;margin-left:0px;margin-right:0px;text-indent:0px">gdalwarp
-r lanczos -order 1 -co COMPRESS=LZW -dstalpha
"/Temp/IGNF_PVA_1-0__1919__CCF000C-71_1919_CAF_C-7_0006.jp2"
"/Temp/IGNF_PVA_1-0__1919__CCF000C-71_1919_CAF_C-7_0006_georef.tif"</p>
</div>
<br></div></div><span>_______________________________________________<br>
Talk-fr mailing list<br>
<a href="mailto:Talk-fr@openstreetmap.org" target="_blank">Talk-fr@openstreetmap.org</a><br>
<a href="https://lists.openstreetmap.org/listinfo/talk-fr" rel="noreferrer" target="_blank">https://lists.openstreetmap.org/listinfo/talk-fr</a><br>
<br></span></blockquote></div><span><br><br clear="all"><br>-- <br><div><div dir="ltr">Christian Quest - OpenStreetMap France</div></div>
</span></div>
</blockquote></div><br><br clear="all"><br>-- <br><div><div dir="ltr">Christian Quest - OpenStreetMap France</div></div>
</div>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr">Christian Quest - OpenStreetMap France</div></div>
</div>