Keyfi meta veri içermeyen harita görüntüsünü QGIS projesine dönüştürün

Geçmişim: deneyimli programcı, amatör tarihçi, CBS'ye tam başlangıç ​​(bu yüzden QGIS kullanıyorum). Yine de QGIS georeferencer aracını kullanmayı öğrendim.

Girdi: varsayalım böyle bir görüntü var Wikipedia: Dosya: VoyagesOfRabbanBarSauma.jpg . Herhangi bir coğrafi meta verisi olmayan "açık" bir JPG dosyası olduğunu ve yazarın bu tür bilgileri sağlayamayacağını varsayalım. Giriş görüntüsü, altındaki türev çalışmalarımı yayınlamaktan mutluluk duyduğum GNU FDL altında lisanslanmıştır.

Amaç: Wikipedia haritasında gösterilen yolu QGIS'e aktarın, böylece projeksiyonu değiştirebilir, yüksekliği üst üste alabilir, fiziksel coğrafya için etiketler ekleyebilirim.

Uyarı: bir çözüm, haritadaki konumların enlem / boylamını yeni bir CBS projesine girmek olabilir. Bu rotaya gitmek istemiyorum çünkü o yazarın coğrafi bilgileri tahmin etmesinin bir kaydını almak istiyorum.

Ne denedim: QGIS Georeferencer raster JPG dosyasına kıtaların kenarlarında bir sürü nokta ile uyguladım. Kıtanın iç kısmındaki noktaların hizası so-so.

Bilmek istediğim şey: Bu belirli haritanın çarpıklıkları Avrasya'nın kuzey kenarı etrafında oldukça belirgindir ve oldukça standart bir projeksiyon gibi görünüyor. Bana öyle geliyor ki, bu haritayı oluşturan projeksiyon ve projeksiyon parametrelerini tanımlamak ve görüntü pikselleri ile enlem / boylam arasında doğru bir enterpolasyon olmayan dönüşümü kurtarmak için algoritmik bir yol olmalı mı?

Bununla ilgili iki soru görüyorum:

1. Bilinmiyorsa Shapefile'in Koordinat Sistemini Belirleme?
2. Eski çizim üzerinde bilinmeyen koordinat sistemi

Bunlar GDAL veya QGIS içeren bir kaba kuvvet yönteminden bahsediyor gibi görünüyor ve bir CBS yeni başlayan için, biraz emek yoğun görünüyor (temel olarak coğrafi referans ve bunları bilinen projeksiyonlara uydurmaya) ama bunları deneyip görmek isteyip istemediğimi görmek istiyorum onları otomatikleştirin.

Bu GIS.SE'nin her ikisi de benim deneyimlerimin izin verdiğinden biraz daha yüksek ve bu kaba kuvvet yönteminin tam olarak nasıl yapılacağına dair ipuçları çok takdir edilecektir. Örneğin, yukarıdaki 2. soru koordinatların bir listesini enlem / boyuta dönüştürmek için gdaltransform'un nasıl kullanılacağını gösterir. Giriş koordinatları xy piksel konumları mı?

Belirtildiği gibi, örnek haritanız QGIS'de WGS84 / lat lon, EPSG: 4326 olan Plate Carree'dir. Helmert Linear yöntemini kullanarak sadece iki nokta ile QGIS'de güzel referanslar. Benim durumumda görüntünün sol üst köşesine yakın bir nokta ve görüntünün sağ alt köşesine yakın başka bir nokta kullandım. Coğrafi referans için kullanılan lat / lon koordinatları Google Earth'ten alındı. Coğrafi referans yaparken Latitude'un "y" ve Boylamın "x" olduğunu unutmayın.

Sadece eğlence için, coğrafi referanslı rasterleri QGIS'e yükledim ve hızlı bir şekilde rotayı izledim ve bir çizgi şekil dosyası olarak kaydettim. Daha sonra Natural Earth'ten bir raster görüntüsü yükledim ve bunu rota ile kapladım (ortaya çıkan haritanın bir kısmı için ekteki resme bakın).

Doğal Earth raster görüntüleri, güzel ve ücretsiz olarak coğrafi referanslıdır (EPSG: 4326'ya), buradan alabilirsiniz:

Doğal Toprak

Tarihsel verileri çizmek için temel haritalar olarak kullanmak için idealdirler.

Daha sonra eklendi:

Dünya ölçeğinde haritalama hakkında hiçbir şey bilmiyorum ama orijinal WGS84 lat / lon (EPSG: 4326) haritasındaki bariz bozulma hakkındaki yorumlarınızı not ettim. QGIS, vektör ve raster verilerini yeniden üretmek için iyi bir araç setine sahiptir, bu yüzden (kısmen kendi kendi eğitimim için!) Haritayı daha "doğal" bir şekilde sunan bir CRS bulmaya çalışacağımı sanıyordum.

Ekli resim, özel bir Albers eşit alanlı projeksiyona dönüştürülmüş kırpılmış bir Doğal Dünya rasterini göstermektedir. Haritadaki vektör verileri, hareket halindeyken WGS84 lat / lon'dan dönüştürüldü. Neden Albers? Denemek için makul bir şey gibi görünüyordu. Açıkçası, bitmiş bir haritada vektör graticule (ayrıca Doğal Dünya'nın izniyle) görüntüye kırpılır. Yine de, haritadaki şekillerin makul derecede doğalcı olduğunu ve haritanın Bar Sauma'nın dikkat çekici yolculuğunun bir dünya üzerinde yapıldığı gerçeğini de vurguladığını düşünüyorum.

Özel CRS için parametreler aşağıdaki gibidir:

+proj=aea +lat_1=22.0 +lat_2=56.0 +lat_0=20.0 +lon_0=55.0 +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs

Bu, kullanımı makul bir tahmin miydi? Dürüst olmak gerekirse bilmiyorum.

Nick.

Harita, Kartezyen koordinatlarıymış gibi enlem ve boylam kullanıyor. Bu, dikdörtgen şeklinde veya "Plate Carree" projeksiyonunun bir versiyonudur . Böylece tek yapmanız gereken, boylam ve enlem elde etmek için sırasıyla x ve y koordinatlarını kaydırmak ve yeniden ölçeklendirmektir.

ArcGIS 10.0'daki coğrafi referans aracı, dönüşümü tahmin etmek için bağlantılar oluşturma konusunda kısa bir çalışma yapar. Dönüşümü rapor etmese de, bağlantı tablosunu dışa aktarabilir Rve oradaki katsayıları tahmin edebiliriz. Orta pikseli , görüntünün (0,0) yüksek çözünürlüklü sürümünün sol üst köşesine yerleştirdiğinizde sonuç,

latitude = 75.71992178 + 0.0427073 * row index
longitude = -18.86514764 + 0.04275506 * column index

Referans katmanının kara ülke ana hatları (ESRI tarafından ArcView'in önceki sürümleriyle birlikte sağlanır), böyle bir dönüşüm uygulandıktan sonra katmanların ne kadar iyi çakıştığını gösterir. Kök ortalama kare hataları, her eksen boyunca 6 - 7 km'dir: görüntüde yaklaşık bir piksel, neredeyse birinin beklemeye hakkı olduğu kadar iyidir.

Haritanın şematik olması muhtemeldir ve şehirlerin yerleşimi, coğrafi referanslar sırasında onları kullanmak için riskli hale gelmeyebilir.

En iyi yaklaşım, haritayı QGIS içinde uygun bir temel kullanarak gözle yakalamak olabilir. Orijinal ile mükemmel bir şekilde eşleşmeyecektir, ancak yolların gösterici olması muhtemeldir, bu nedenle herhangi bir doğruluğu kaybetmeniz olası değildir.