QGIS'de fiziksel bir bariyer / sınır ile çekirdek yoğunluğu tahmini nasıl yapılır?

18

Muhtemelen imkansız olan alanları maskelemek yerine, hesaplamasında bir sınırı dikkate alan çekirdeklendirmeye izin verecek bir program biliyor mu?

Şimdiye kadar buldum: Arc10 (Hawthorne Beyer) için GME eklentisi - Her bir sınır şekil dosyası belirlediğimde bir hata alıyorum. Bir çok farklı şekil dosyası türleri ve sınır karmaşıklıkları vb denedim. Bir sınır belirtmediğimde iyi çalışıyor.

AdehabitatHR paketi R (Calange 2011) - bu iyi çalışıyor, ancak belirttiğiniz sınır çok basit olmalı - 3x çekirdek bant genişliği uzunluğunda ve çok kıvrımlı değil. Verilerim için bu büyük bir aşırı basitleştirme.

Başka bir yazılımın bunu yapıp yapamayacağını merak ediyorum, örneğin GRASS veya QGIS.

Teşekkürler

qgis grass

— Kims
kaynak

3

Potansiyel ilgi R kodu içeren son makale, Sınır sapması ve ağırlıklı çekirdekler . Tüm kod, sınır sapmasını düzeltmek için gözlemleri ağırlıklandırmanın sadece yoludur (ağırlıkları anlayabiliyorsanız, ağırlıkları hariç olduğu sürece kde'yi tahmin etmek için hangi programı kullanabilirsiniz).

— Andy W

1

Tam olarak bir sınırı nasıl "açıklamasını" istersiniz? Maskelemeden yayılmayı engellemeye kadar sınır etkileri düzeltmeye kadar birçok olası yol vardır.

— whuber

1

Cevaplar için @AndyW ve whuber'a teşekkür ederim. Çekirdeğin yayılmasını engelleyerek sınırı hesaba katmak istiyorum. Veriler deniz kaplumbağası yerleridir, bu yüzden kullanım dağılımının toprağa yayılmaması gerektiğini biliyorum, ancak yerlerin çoğu sahile çok yakın olduğundan, sadece çekirdek kısımlarını maskeleyerek bu konumlara karşı önyargılı olmak istemiyorum. karaya yayıldı. Makaleye bağlantı için teşekkür ederim - bu harika görünüyor, verilerim için benzer kodu kullanmaya çalışacağım.

— KimS

2

Çekirdek yoğunluğu tahmininin merkezinde mesafe kavramı vardır. Bildiğim en iyi çözüm, sınırları ve değişen seyahat maliyetlerini açıklayan daha iyi bir mesafe metriği kullanmaktır. Çözmeye çalıştığınız soruna uyan bir mesafe metriği seçmek en iyisidir. Arazi arazisi sürtünmesi yürüyüş için mükemmeldir, ancak aerosol dağılımı için değildir. Rüzgar akıntıları yelkenli takibi için gereklidir, ancak yol tarifleri için önemsizdir.
Artık uygun mesafe ölçütleri kavramının yeterince motive edildiğine göre, maliyet yüzeylerini iyi bir genel amaçlı mesafe ölçütü olarak önerebilirim. ArcGIS'ten R'ye, JavaScript'e kadar her şeyde mevcuttur ve yapılandırılması oldukça kolaydır. Örneğin qGIS'de, bir raster sürtünme yüzeyi oluşturabilir ve bunu rotaları hesaplamak için kullanabilirsiniz. Sınırlarınızı hesaba katmak için sürtünme yüzeyini özelleştirin ve noktalarınızın etrafındaki çekirdeklerden gelen kütlenin engellerin etrafına düzgün bir şekilde yayıldığını göreceksiniz.

— JasonRDalton
kaynak

+1 Bu potansiyel olarak harika bir yaklaşımdır. Birçok problemde uyguladıktan sonra, asıl zorluğun burada belirtilmediğini gördüm: tam olarak çekirdekleri bu metriğe göre yaymayı nasıl öneriyorsunuz? Hazır çözümler bunu yapamaz.

— whuber

@whuber Doğru, kullanıma hazır sistemler genellikle hem maliyet yüzeylerini hem de çekirdek yoğunluğunu aynı anda hesaba katamaz. Benim önerim, noktadan noktaya mesafeyi hesaplamak için sürtünme yüzeyinin kullanılması, daha sonra doğrudan ağırlık için bir çekirdek fonksiyonu uygulanması olacaktır. QGIS veya ArcGIS'de python'daki normal dağıtım işlevini kullanarak yazabilirsiniz ve R'de bu dnorm()işlev vardır.

— JasonRDalton

Bazı durumlarda tam olarak bunu yaptım, ancak süreci düşünmek ve sonuçları incelemek bazı doğal sorunları ortaya çıkarıyor. En zoru, kütlenin korunmasını nasıl garanti edeceğidir: sabit bir Gauss fonksiyonuna göre değerleri sadece yayamazsınız, çünkü bir yayılma noktasının toplam kütle kütlesi orijinal kütleye eşit olmaz. Bu hesaplamayı gerçekleştirmenin etkili bir yolu yoktur.

— whuber

0

R için spatstat paketini denemenizi öneririm owin. Çalışma alanının sınırlarını belirlemek için bir nesne ayarlayabilirsiniz . Ayrıca bu paket için harika bir öğretici var .

— SS_Rebelious
kaynak

0

Aslında çekirdek yoğunluğu tahmini için R bağımlılığından kurtulmayı amaçlayan QGIS için AniMove eklentisine yardım ediyorum. Buraya bir bak .

EDIT: Eklenti aslında resmi QGIS eklenti deposunda deneysel olarak kullanılabilir

Ve herhangi bir öneriniz varsa , listede sormaktan çekinmeyin

— Jorge Arévalo
kaynak

-2

ESRI ArcGIS'in bir uzantısı olan Geostatistik Analisti kullanarak bunu başarıyla yaptım. Çizgi özellikleri veri kümesini bariyerler olarak yükleyebilirsiniz ve sonuçlar oldukça iyi. Aracın nasıl kalibre edileceğine dair net bir fikir edinmek için fonksiyonların parametrelerini değiştirirken test çalışmaları yapıyorum. Sonuçları doğrulamak için, önerdiğim şey, yeterince büyük bir veri kümeniz varsa, bir nokta örneğini kaldırmak, bu noktalar olmadan yoğunluk yüzeyini oluşturmak, daha sonra kaldırılan noktaların değerleri ile yüzeydeki yüzey değerleri arasındaki farkı karşılaştırmaktır. kaldırılan noktaların konumu.

— Petya
kaynak

2

İnterpolasyon yöntemlerini uygulayan GA'nın çekirdek yoğunluğu tahminini nasıl hesaplayabildiğini açıklamak için bu cevabı biraz güçlendirebilir misiniz ? Çözümünüzün doğruluğunun iki kez kontrol edilmesi olarak, GA ızgarasının integralinin yaklaşık olarak tüm giriş verilerinin toplamına eşit olduğunu doğruladınız mı?

— whuber

1

10.1'de bariyerlerle çekirdek enterpolasyonu yapmak mümkündür, ancak bu bir yoğunluk tahmini ile sonuçlanmaz. Her ne kadar, orijinal yazının referans aldığı yöntemler göz önüne alındığında, bir Gauss çekirdek tahmini ile ilgilenip ilgilenmediklerini merak ediyorum.

— Jeffrey Evans