|
|
| Type : | script pour la boite à outils traitement |
| Licence : | GPL |
Suite à cette discussion sur le forumSIG http://www.forumsig.org/showthread.php/37661-Logiciel-Ajustement-spatial... , où je cherchais une méthode pour faire du géoréférencement vectoriel, voici le protocole mis en place pour y parvenir.
Attention, cette méthode nécessite que QGIS utilise gdal 1.10 ou supérieur.
Faire un fichier texte de points de correspondance (une ligne par point) formaté ainsi :
pointNonGeoref_x pointNonGeoref_y pointGeoref_x pointGeoref_y
Pour cela à l'aide de l'outil identification / dérivées / copier les valeurs d'attributs (on fait ça pour le fichier à géoréférencer et sur le référentiel à utiliser pour géoréférencer, le plus simple étant d'avoir 2 instances de qgis simultanées, l'une avec le fichier à géoréférencer sans reprojection à la volée et l'autre avec le référentiel permettant de récupérer les bonnes coordonnées)
Exemple
-869784.220806 273824.996155 475353.856328 267608.211198
-869786.957227 273804.951538 475351.619893 267590.139363
-869767.137134 273845.613669 475369.467973 267625.62811
-869758.781489 273830.170332 475377.513896 267610.888946
A noter qu'il est également possible d'utiliser le plugin copy_Coords de QGIS, plus simple d'utilisation, mais moins précis.
On lance ensuite le script joint (celui-ci ne fait que remplacer la commande d'ogr2ogr proposée par gene http://www.forumsig.org/showthread.php/37661-Logiciel-Ajustement-spatial...) mais permet d'éviter la ligne de commande.
Décompresser le zip dans le dossier script de l'outil traitements (voir Traitements / options / scripts / Scripts folder)
Dans la boite à outils de traitements / scripts, une entrée "[Géoréférencement] / ajustementSpatial" a du apparaitre, lancer le script et remplir comme précisé dans la documentation.
Si tout se passe bien, le fichier nouvellement géoréférencé doit s'ouvrir, placé correctement.
Le nombre de points de correspondance définit l'ordre du polynôme utilisé pour le géoréférencement.
Le script a été testé sous linux et sous windows avec l'installateur OSGEO.
Sous Mac, il est nécessaire de remplacer la ligne
ogr = 'ogr2ogr -a_srs EPSG:'+epsg+' '
par
ogr = '/Library/Frameworks/GDAL.framework/Versions/1.10/Programs/ogr2ogr -a_srs EPSG:'+epsg+' '
En pratique, dans le script téléchargé, commenter la première ligne en y insérant un # devant et décommenter la ligne pour Mac en supprimant le #.
En espérant que cela vous sera utile, n'hésitez pas à faire un retour si vous avez des soucis sur la discussion du forumSIG.
