DEM'den pens çıkarımı mı yapıyorsunuz?

DEM'in yapıldığı 1x1 metrelik DEM ve .las içinde orijinal LiDAR nokta bulutum var. Vektör özelliği (nokta, çoklu çizgi) için nehir penseyi (pensenin en yüksek noktaları) ayıklamak gerekir.

Algoritma veya mevcut bir araç hakkında herhangi bir fikriniz var mı?

Dikenler 1. görüntüde açık gri renktedir ve aşağıda, dikenlere sahip örnek alanı için nokta bulutu görüntüsüdür. Mavi nehir merkezinde.

DEM doğrusal özellikleri çıkarma için herhangi bir araç?

ENVI'nın Mekansal Özellik Çıkarma Modülü (sayfa 7) gibi bir şeye ihtiyacım olduğuna inanıyorum , ancak bunu karşılayamadığım için ENVI parçası olmadan :)

Her şey çizgiyi nereye çizdiğinize bağlıdır. Ne olursa olsun, bu problem , Spatial Analyst'ta mevcut olan morfolojik fonksiyonlar , özellikle bileşenleri tanımlamak ve çıkarmak için özellikle eşikleme ("<" ve ">" yerel işlemler ile gerçekleştirilir) ve "RegionGroup" kullanılarak kolayca çözülebilecek gibi görünüyor .

Göstermek için DEM'e erişimim olmamasına rağmen, ilk görüntü zaten işi yapacak kadar iyi. Örneğin, burada 0 bandından başlayarak ve 0,02'lik artışlarla soldan sağa, yukarıdan aşağıya doğru hareket eden kırmızı bandın bölge gruplu sürümlerinin (0.0 = siyahtan 1.0 = beyaza ölçeklendirilmiş) bir sırası vardır:

(Bu görüntüler burada çoğaltılmak üzere küçülmüştür: tüm analizler orijinal görüntünün çözünürlüğünde yapılmıştır.)

Pervazlar, ikinci sıranın başlangıcında siyah bölgenin sınırı etrafında ortaya çıkar (eşik = 0.08). Üçüncü sıranın başlangıcında (eşik = 0.16), pensler kendi bileşenlerini (koyu mavi olarak) oluşturur ve bu noktada ayrı ızgaralar veya çokgenler olarak kolayca çıkarılabilir (ve kenarları bir kenar algılama adımından sonra polienler olarak çıkarılabilir ). Son satırda yalnızca en geniş pensler kalır (eşik = 0,24 ve daha büyük). Tam olarak istediğinizi elde etmek için uygun bir eşik seçmeniz gerekecektir .

Orijinal DEM'de, yükselmeler yoğunluk rolünü oynar, bu nedenle bu prosedürler DEM'in kendisi ile benzer şekilde etkili olmalıdır. DEM'in büyük bir kısmı varsa (nehir veya dikenlerden uzak özellikler dahil), ilgili özellikler nehir özelliğinin bulunduğu bileşene bitişik olarak seçilebilir.

RegionGroup tarafından bulunan ince istenmeyen bölgeler, küçük bir erozyon (negatif tampon) ve ardından eşdeğer bir dilatasyon (pozitif tampon) uygulanarak ortadan kaldırılabilir. Küçük istenmeyen bölgeler bir boyut kriterinde (toplam hücre sayısı veya alan) hariç tutulabilir. Her bir dike üzerindeki en yüksek noktalar (gerçekten gerekli olan buysa), DEM'i bölgesel maksimum ızgaralarla karşılaştırarak bulunabilir (dikenleri bölge olarak kullanarak).

Bu arada, bu görüntü dizisi Mathematica 8'de üretildi . İşte bu seçeneği sürdürmek isteyenler için komutlar.

original = Import["http://i.stack.imgur.com/gV7Du.jpg"];
{r, g, b} = ColorSeparate[original];
frames = ParallelTable[Colorize[MorphologicalComponents[Binarize[r, t/100]], 
                       ColorFunction -> "ThermometerColors"], {t, 0, 30, 2}];
Rasterize @ TableForm[Partition[frames, 4]]

Toplam saat süresi (görüntüyü içe aktardıktan sonra) 0.94 saniyeydi, bu sürenin 16 görüntünün tümünü daraltmak ve dışa aktarmak için yarısı gerekliydi: morfolojik işlemler hızlı olma eğilimindedir (LIDAR DEM'ler çok büyük olabilir).

İş için çeşitli algoritmalar hakkında okudum (yani @Hornbydds bağlantısına göre ).

Ben birkaç appoches denedim ve benim durumumda en iyi sonuçlar SAGA Standart Arazi Analizi verir. İşte yaptığım ve neden:

Dikler genellikle nehir kanalının en yüksek özelliğidir, bu yüzden onları DEM'i çevirerek kanallara dönüştürdüm (MapAlgebra DEM * -1 veya -10 yamaçlarının abartılması için). Bu noktadan sonra herhangi bir hidrolojik alet seti (ArcHydro, HEC-GeoRas veya SAGA hidro aletleri) kullanabilirim. SAGA / Terrain Analysis -compound analizi / Standard Terrain Analysis'i çift tıklamayla River Network ürettiği için seçtim. Nehir Ağı, tersine çevrilmiş penslerde en yüksek noktayı veren su akışı için en iyi rotayı gösteren polinler ürettiği için ulaşmak istediğim şeydir.

Üretilen şekil dosyası biraz dağınık (çoğu daha küçük polylines), ancak bazı ince ayarlardan sonra tatmin edici. Bir başka sıkıntı, 1 lezbiyen 300 kısa polietilenden yapılmış, ancak bununla başa çıkmanın bir yolunu bulacağımı düşünüyorum.

İşte sonuçlar örneği:

Tabii ki bu oldukça kaba bir çözüm ve muhtemelen konuyu daha fazla araştıracağım, ancak paylaşmaya değer olduğunu düşündüm.

Aşağıdakilerin orijinal DEM ile bazı kombinasyonları, istediğiniz tarama bilgisini vurgulamaya ve geri kalanını atmaya yardımcı olmalıdır. Bunlar, etkiyi iyileştirmek için farklı ölçeklerde ve mahalle boyutlarında gerçekleştirilebilir.

• Yeniden örnekleyin, ardından Mahalle -> Filtre: Yüksek Geçiş
• Yüzey -> Eğrilik
• Mahalle -> Odak İstatistikleri: Standart Sapma
• "Göreli Yükseklik" = (Piksel) - (Semt -> Odak İstatistikleri: Medyan)
• Ters Çevir, Yeniden Örnekle, ardından Komşu -> Odak Akışı

Bunu yaptıktan sonra, yeniden sınıflandırma ve konturlar, lezbiyenlerin kendisinin oldukça iyi bir taslağını vermelidir. Penseler düz eğimli ise ve daha geniş bir kapsama alanı istiyorsanız, eğrilik ve eğim toplamı gibi bir miktar ölçeklendirilmiş bir şey deneyebilir, hatta eğim konturlarını ayrı ayrı yapabilir ve bunları bir araya getirebilirsiniz.

Burada, yüksek oranda uygulanabilir olması gereken Spatial Analyst'ı kullanarak bir raster DEM'den vektör ridgelinlerin çıkarılmasıyla ilgili bir öğretici:

• http://support.esri.com/es/knowledgebase/techarticles/detail/39093

Bunun sadece cevabın bir parçası olabileceğini düşünüyorum, ancak belki belirli aralıklarla / yükseklik değerlerine dayalı olarak konturlarınızı / çokgenlerinizi görüntülerinizden çıkarabilirsiniz. Direğin daha yüksek yükseklik değerlerine sahip olması gerektiğini düşünüyordum ve bu nedenle sizi DEM'i çokgenlere dönüştürdüğünüzde, ilgilendiğiniz alanların çokgen çokgenlerle temsil edildiğini görebilirsiniz. Ayrıca, yalnızca DEM'yi yeniden sınıflandırarak ve bunun peşinde olduğunuz sonuçları üretip üretmediğini görebilirsiniz. İlgilendiğiniz alanı kapsayan başka görüntü verileriniz varsa, görüntü sınıflandırmanıza eklemeyi deneyin ve yardımcı olup olmadığını görün. İyi şanslar!!