BINGO ! 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...
En postgis ça donne un truc comme: (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 Le 6 mars 2016 à 20:32, Christian Quest <cqu...@openstreetmap.fr> a écrit : > J'ai récupéré les photos de 1921 prises au dessus de ma commune, il y en a > une vingtaine. > > En gros c'est plutôt correct sauf que pour 3 clichés on est à 180° :( > > En rouge sur la copie d'écran QGis: > https://framapic.org/MMcURbVTytm9/iZoQxP0IIY5Y.png > > Malheureusement, entre le polygone d'emprise et l'orientation, je n'ai pas > trouvé de critère pour les détecter. > > > Le 6 mars 2016 à 18:17, Christian Quest <cqu...@openstreetmap.fr> a écrit > : > >> Pour moi c'est ok en procédant ainsi: >> >> { >> "type": "FeatureCollection", >> "features": [ >> {"type":"Feature", >> "geometry":{ >> "type":"Polygon", >> "coordinates":[[[2.5126989,48.8199908], >> [2.5213251,48.8149892], >> [2.5298119,48.8101662], >> [2.5203523,48.8023041], >> [2.5102159,48.7943484], >> [2.5013197,48.7993482], >> [2.4924227,48.8042573], >> [2.5025584,48.8122145], >> [2.5126989,48.8199908]]]}, >> >> "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"}}] >> } >> >> >> w=6774 >> h=5050 >> gdal_translate -of GTiff >> IGNF_PVA_1-0__1921-06-11__CCF00B-361_1921_CAF_B-36_0144.jp2 temp.tif \ >> -gcp $w 0 2.5126989 48.8199908 \ >> -gcp $w $h 2.5298119 48.8101662 \ >> -gcp 0 $h 2.5102159 48.7943484 \ >> -gcp 0 0 2.4924227 48.8042573 >> 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 >> >> Il y a 9 coordonnées, les coins sont en positions impaires (en comptant >> de 1) et ça démarre en haut à droite... >> >> >> Le calage n'est pas nickel, mais c'est une bonne première approximation ! >> >> Exemple ici: https://framapic.org/dLIvBVa9WTJz/RXxKlfWqAwrk.png >> >> >> >> Le 6 mars 2016 à 18:02, <osm.sanspourr...@spamgourmet.com> a écrit : >> >>> On y arrive, je vais laisser les Normands remettre la photo à l'endroit. >>> J'ai utilisé Qgis en affichant les coordonnées des points de contrôle, >>> càd étiquettes : name || ' (' || $x || ' ' || $y || ')'. >>> >>> >>> Ce qui équivaut à : >>> 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" >>> >>> >>> 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" >>> >>> _______________________________________________ >>> Talk-fr mailing list >>> Talk-fr@openstreetmap.org >>> https://lists.openstreetmap.org/listinfo/talk-fr >>> >>> >> >> >> -- >> Christian Quest - OpenStreetMap France >> > > > > -- > Christian Quest - OpenStreetMap France > -- Christian Quest - OpenStreetMap France
_______________________________________________ Talk-fr mailing list Talk-fr@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-fr
