Mon parcours: programmeur expérimenté, historien amateur, débutant complet en SIG (donc j'utilise QGIS). J'ai cependant appris à utiliser l'outil de géoréférencement QGIS.
L'entrée: supposons que j'ai une image comme celle-ci Wikipedia: File: VoyagesOfRabbanBarSauma.jpg . Supposons qu'il s'agit d'un fichier JPG "nu" sans métadonnées géographiques et que l'auteur ne peut pas fournir de telles informations. L'image d'entrée est sous licence GNU FDL, sous laquelle je suis heureux de publier mes travaux dérivés.
Le but: importer le chemin indiqué sur la carte Wikipedia dans QGIS pour pouvoir modifier la projection, superposer l'altitude, ajouter des étiquettes pour la géographie physique, etc.
La mise en garde: une solution de contournement pourrait simplement être d'entrer le lat / long des emplacements sur la carte dans un nouveau projet SIG. Je ne veux pas emprunter cette voie parce que j'aimerais avoir un enregistrement de l'estimation de cet auteur particulier des informations géographiques.
Ce que j'ai essayé: j'ai appliqué le QGIS Georeferencer au fichier JPG raster avec un tas de points sur les bords des continents. L'alignement des points à l'intérieur du continent est moyen.
Ce que j'aimerais savoir: les distorsions de cette carte particulière sont assez visibles autour de la bordure nord de l'Eurasie, et il semble que ce soit une projection assez standard. Il me semble qu'il devrait y avoir un moyen algorithmique d'identifier la projection et les paramètres de projection qui ont généré cette carte, et de récupérer une conversion de non-interpolation précise entre les pixels de l'image et les lat / long?
Je vois deux questions à ce sujet:
- Identifier le système de coordonnées du fichier de formes lorsqu'il est inconnu?
- Système de coordonnées inconnu sur l'ancien dessin
Ceux-ci semblent parler d'une méthode de force brute impliquant GDAL ou QGIS, et pour un débutant SIG, semblent quelque peu laborieux (essentiellement géoréférencer et les adapter aux projections connues) mais je suis prêt à les essayer et voir si je peux les automatiser.
Les instructions sur ces deux GIS.SE sont un peu plus élevées que mon expérience ne le permet, et tout indice sur la façon de faire exactement cette méthode de force brute serait très apprécié. Par exemple, la question 2 ci-dessus montre comment utiliser gdaltransform pour convertir une liste de coordonnées en lat / long. Les coordonnées d'entrée sont-elles des emplacements de pixels xy?
la source
La carte utilise la latitude et la longitude comme s'il s'agissait de coordonnées cartésiennes. Il s'agit d'une version de la projection équirectangulaire ou "Plate Carree". Ainsi, tout ce que vous avez à faire est de décaler et de redimensionner les coordonnées x et y pour obtenir respectivement la longitude et la latitude.
L'outil de géoréférencement d'ArcGIS 10.0 permet de créer rapidement des liens pour estimer la transformation. Bien qu'il ne signale pas la transformation, il est possible d'exporter la table de liens
R
et d'y estimer les coefficients. En plaçant le pixel central dans le coin supérieur gauche de la version haute résolution de l'image à (0,0), le résultat estLes contours du pays noir de la couche de référence (fournis par ESRI avec les versions antérieures d'ArcView) montrent à quel point les couches coïncident après l'application d'une telle transformation. Les erreurs quadratiques moyennes sont de 6 à 7 km le long de chaque axe: environ un pixel dans l'image, presque aussi bon que l'on est en droit d'attendre.
la source
La carte est susceptible d'être schématique et l'emplacement des villes peut ne pas être parfait, ce qui les rend risquées à utiliser lors du géoréférencement.
La meilleure approche pourrait être de capturer la carte à l'œil en utilisant une base appropriée dans QGIS. Il ne correspondra pas parfaitement à l'original, mais les itinéraires sont susceptibles d'être indicatifs, il est donc peu probable que vous perdiez toute précision.
la source