QGIS'te ağırlıklı noktalardan yoğunluk haritası oluşturmanın en iyi yolu nedir?

ArcGIS Çekirdek Yoğunluğu aracının işlevselliğini tekrarlamak istiyorum - bunu nasıl yapabilirim?

Yüklemesi kolay ücretsiz ve açık kaynaklı yazılımlara verilen yanıtları sınırlamak istiyorum - yani QGIS, tüm platformlara GRASS ile kolayca yüklenir, böylece QGIS + GRASS iyi olur, ancak SAGA bunu yapmaz (ne yazık ki, ideal bir araç olabilir).

QGIS'te korunan alanlarda vahşi yaşam yoğunluk haritaları oluşturmaya çalışıyorum. Aşağıdaki örnekte yoğunluk haritası, Çekirdek Yoğunluğu (SpatialAnalyst) aracındaki her ızgara hücresinin yoğunluğunu ağırlıklandırmak için kullanılan grup büyüklüğü (vektör dosyasındaki bir alan) ham sayılarıyla ArcGIS'de vahşi yaşamın nokta gözlemlerinden üretildi. seçilen bir arama yarıçapı ve ızgara hücre boyutu:

Yoğunluk ile ilgili önceki bir soruda , ArcGIS Çekirdek Yoğunluğu aracını taklit etmek için GRASS v.kernel aracını kullanmanız önerildi, ancak v.kernel aynı işi yapmıyor . El kitabına baktıktan ve (başarıyla) bir yoğunluk haritası ürettikten sonra, v.kernel işlevinin sadece nokta yoğunluğu ile çalıştığı ve her noktayı (ham sayımlar gibi) her noktayı ağırlıklandırmak için bir değişken verme fırsatı olmadığı anlaşılıyor. .

GÜNCELLEME

*.surf.*GRASS'ta yoğunluk yüzeyi oluşturmaya yardımcı olabilecek çeşitli araçlar var gibi görünüyor ve bunlar ağırlık sütununu veya z değerini kabul ediyor veya rasterler üzerinde gerçekleştiriliyor. @underdark, v.surf.rst'yi önerdi - ve 'zcolumn' ağırlıklandırma (sayım) değişkenim olacaktı - ancak araçtan belirli bir ızgara boyutu oluşturmasını veya belirli bir yarıçapı nasıl kullanmasını isteyeceğim.

V.surf.rst veya başka bir yöntem nasıl kullanılır?

Örnek veri

x,y,count
431250,8707500,0
418750,8707500,5
413750,8707500,3
411250,8707500,1
408750,8685000,0
411250,8685000,0
416250,8685000,0
416250,8682500,6
411250,8682500,3
418750,8680000,0
433750,8677500,3
421250,8677500,0
423750,8675000,1
431250,8672500,0
428750,8672500,2
426250,8672500,2
423750,8670000,0

Kılavuz sayfasına göre, GRASS komutu r.resamp.filter , ArcGIS'in nokta katmanları için tam olarak ne yapacağını gösteren nokta verilerini temsil eden rasterler için yapacaktır : "basit" bir raster seçeneğini kullanın ve diğer ArcGIS çekirdeği seçeneğini kullanın . Boşları yaymaktan kaçınmak için bayrağı kullanın .filter=boxfilter=gauss-n

O Not çekirdek yoğunluk tahminleri (aka "ısı haritaları") vardır değil veri interpolasyonları. Bir yerde bir KDE değeri x değeri "miktarını tahmin Z yakın birim alan başına" x . (Yarıçap veya "bant genişliği", "yakın" öğesinin ne anlama geldiğini belirler.) Z değerlerinin , haritadaki her olası konumda tanımlanması gerekmez. Örneğin Z , kişi gibi bir şeyin varlığını temsil edebilir, bu durumda KDE nüfus yoğunluğu verir . Z değerlerinin harita üzerinde sürekli değişmesi de gerekmez. Enterpolasyon için, bir farz olduğunu Ztüm konumlarda tanımlanır ve verilerin Z'deki değerlerin belirtilen noktalardaki gözlemleridir . İnterpolatör , diğer tüm noktalarda Z'nin gözlemlenmemiş değerlerini tahmin etmeye çalışır . Bu, Z , örneğin bir sıcaklık veya basınç olduğunda anlamlıdır , ancak Z bir şeyin varlığını kaydettiğinde veya veriler tam bir sayım olduğunda genellikle saçmadır . (İkinci durumda, bir bölge için bir yol yoğunluğu haritasının ne anlama gelebileceğini ve yol dışı alanlarda yolların "enterpolasyonlu" yollarını nasıl anlamlı hale getirebileceğini düşünün.)

SAGA'nın 'Çekirdek Yoğunluğu Tahmini' modülü aradığınız şeydir.

SAGA Modülü arayüzünü QGIS'e kurun (menü: Eklentiler -> Python eklentilerini getir ..) ve modülü kullanın. İyi şanslar!

ArcMap'ta Çekirdek Yoğunluğundan Nokta Yoğunluğuna daha yakın olan GRASS GIS ile gerçekten basit bir yöntem:

1. r.in.xyzBelirtilen raster çözünürlüğünde (ile ayarlanır g.region) method = sum komutunu kullanarak noktaları bir raster haritasına alın .

2. r.neighborsHaritayı yöntem = ortalama (varsayılan) ile pürüzsüzleştirmek için kullanın ve arama yarıçapını ayarlamak için seçenek boyutunu kullanın.

(Şu anda GRASS'a erişimim yok, bu yüzden bunu gerçekten denemedim!)

V.surf.rst hakkında daha fazla rehberlik talep ettiğiniz gibi, işte girdilerim

İlk olarak, ızgara boyutu hakkında - Eklenti -> GRASS -> Mevcut GRASS bölgesini düzenleyebilir ve çıkış çözünürlüğünü ayarlayabilirsiniz. V.surf.rst dosyasından aldığınız çıktı bu çözünürlüğe sahip olacaktır.

Yarıçap için 'gerginlik' parametre gibi görünüyor. Bu algoritma konusunda uzman değilim, ancak kılavuzdan okumak, bu ilgili bit gibi görünüyor

"... Yüksek gerilim" noktalar arasındaki mesafeleri arttırır "ve her noktanın etki aralığını azaltır, düşük gerilim" mesafeyi azaltır "ve noktalar daha uzun aralıkta birbirini etkiler) .."

Böylece, radius parametresini kullandığınız gibi gerilim parametresini kabaca kullanabilirsiniz.

Örnek verilerinizden, v.surf.rst sonucu aşağıdaki gibi görünür ve sayıları enterpolasyon için ağırlık olarak kullandığından makul görünür

Test etmeme rağmen, QGIS Katkıda Bulunanları deposunda 'R ile ev aralığı tahmini' adlı bir eklenti var. Bu, Çekirdek (yoğunluk?) Hesaplamalarını içerir. Eğer işe yararsa, bunun en iyi yöntem olacağını tahmin ediyorum. R çekirdek yoğunluğunu hesaplamak için gerçek istatistiksel yöntemi yapacaktır.

R'yi yüklediyseniz, eklentiyi yükleyip bir şans vermeniz gerekir.

Qgis dışında biraz java programlama yapmayı kabul ederseniz, bu yoğunluk haritası oluşturucu kitaplığını kullanabilirsiniz. .

Yapıcıyı kullanarak, HeatMapBuilder(int w, int h, int[][] pts, int[] weights)ihtiyacınız olan her nokta için bir ağırlık vermek mümkündür. Çıktı görüntüsü yöntemle alınabilir getImage()ve birImageIO.write("mymap.png") .

İşte çıktı örneği:

Düzeltme gücünü ve renk paletini değiştirmek mümkündür.