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
_______________________________________________ Talk-fr mailing list Talk-fr@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-fr