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

Répondre à