Mein Hintergrund: erfahrener Programmierer, Amateurhistoriker, absoluter Anfänger in GIS (also benutze ich QGIS). Ich habe jedoch gelernt, das QGIS-Georeferencer-Tool zu verwenden.
Die Eingabe: Angenommen, ich habe ein Bild wie dieses Wikipedia: Datei: VoyagesOfRabbanBarSauma.jpg . Angenommen, es handelt sich um eine "nackte" JPG-Datei ohne geografische Metadaten, und der Autor kann keine solchen Informationen bereitstellen. Das Eingabebild ist unter GNU FDL lizenziert, unter dem ich gerne meine abgeleiteten Werke veröffentlichen kann.
Das Ziel: Importieren Sie den in der Wikipedia-Karte angezeigten Pfad in QGIS, damit ich die Projektion ändern, die Höhe überlagern, Beschriftungen für die physische Geografie hinzufügen usw. kann.
Die Einschränkung: Eine Problemumgehung besteht möglicherweise darin, die Breite / Länge der Positionen auf der Karte in ein neues GIS-Projekt einzugeben. Ich möchte diesen Weg nicht gehen, weil ich die Schätzung des jeweiligen Autors der geografischen Informationen aufzeichnen möchte.
Was ich versucht habe: Ich habe den QGIS Georeferencer auf die Raster-JPG-Datei mit einer Reihe von Punkten an den Rändern der Kontinente angewendet. Die Ausrichtung der Punkte im Inneren des Kontinents ist mittelmäßig.
Was ich gerne wissen würde: Die Verzerrungen dieser speziellen Karte sind am nördlichen Rand Eurasiens ziemlich auffällig, und es scheint, als wäre es eine ziemlich normale Projektion. Es scheint mir, dass es einen algorithmischen Weg geben sollte, um die Projektions- und Projektionsparameter zu identifizieren, die diese Karte erzeugt haben, und eine genaue Nichtinterpolationskonvertierung zwischen Bildpixeln und Lat / Long wiederherzustellen?
Ich sehe zwei Fragen dazu:
- Identifizieren des Koordinatensystems von Shapefile, wenn unbekannt?
- Unbekanntes Koordinatensystem auf alter Zeichnung
Diese scheinen über eine Brute-Force-Methode mit GDAL oder QGIS zu sprechen, und für einen GIS-Anfänger scheinen sie etwas arbeitsintensiv zu sein (im Grunde genommen Georeferenzierung und Anpassung an bekannte Projektionen), aber ich bin bereit, sie auszuprobieren und zu sehen, ob ich kann automatisieren sie.
Die Anweisungen auf diesen beiden GIS.SE sind etwas höher als es meine Erfahrung zulässt, und alle Hinweise, wie genau diese Brute-Force-Methode durchzuführen ist, wären sehr willkommen. Beispiel 2 oben zeigt, wie mit gdaltransform eine Liste von Koordinaten in lat / long konvertiert wird. Sind die Eingabekoordinaten xy Pixelpositionen?
quelle
Die Karte verwendet Breiten- und Längengrade, als wären sie kartesische Koordinaten. Dies ist eine Version der gleichwinkligen oder "Plate Carree" -Projektion. Sie müssen also nur die x- und y-Koordinaten verschieben und neu skalieren, um den Längen- bzw. Breitengrad zu erhalten.
Mit dem Georeferenzierungswerkzeug in ArcGIS 10.0 können Sie in kurzer Zeit Links erstellen, um die Transformation abzuschätzen. Obwohl die Transformation nicht gemeldet wird, kann die Verknüpfungstabelle in
Rdie Koeffizienten exportiert und dort geschätzt werden. Wenn Sie das mittlere Pixel in der oberen linken Ecke der hochauflösenden Version des Bildes bei (0,0) platzieren, erhalten Sie das ErgebnisDie schwarzen Länderumrisse der Referenzebene (von ESRI mit früheren Versionen von ArcView bereitgestellt) zeigen, wie gut die Ebenen nach dem Anwenden einer solchen Transformation übereinstimmen. Die mittleren quadratischen Fehler betragen 6 bis 7 km entlang jeder Achse: etwa ein Pixel im Bild, fast so gut, wie man es erwarten kann.
quelle
Die Karte ist wahrscheinlich schematisch und die Platzierung der Städte ist möglicherweise nicht perfekt, was die Verwendung bei der Georeferenzierung riskant macht.
Der beste Ansatz könnte darin bestehen, die Karte mit einer geeigneten Basis in QGIS per Auge zu erfassen. Es stimmt nicht perfekt mit dem Original überein, aber die Routen sind wahrscheinlich indikativ, sodass Sie wahrscheinlich keine Genauigkeit verlieren.
quelle