1m DEM mozaiklerini 10m DEM'e toplamak için en iyi yöntem

1m LiDar türevli bir DEM mozaiğe sahibim. 10m DEM olarak verilerin bir alt kümesini çıktı gerekiyor. Şu anda ARCGIS 10'daki toplama aracını her yeni 10m piksel için ortalama bir değer üretmek için kullanıyorum. Bunun görev için en iyi teknik olup olmadığı konusunda herhangi bir tavsiye var mı? Ortalama değer bu tip verilerle en iyi yaklaşım mıdır?

Yaygın bir hata (benim de yaptığım), bilinear enterpolasyonlu yeniden örnekleme aracını kullanarak bir rasterden aşağı örneklemektir. Bunun neden iyi olmadığını açıklamak için bu cevaba bakınız . Bir raster üç adımda aşağı örneklenebilir.

1. İlk adım gerekli olmayabilir. Rasteri hedef uzantılara yeniden projekte edin. Bilineer enterpolasyon kullanın ve çıkış hücresi boyutunu giriş çözünürlüğü ile aynı tutun (örn. 1 m). Raster köşelerini projeksiyona "tutturmak" için kayıt noktasını kullanın. Çıktı uzantıları "Ortamlar" da belirtilebilir ve uzantıların 10 m'nin katları (veya herhangi bir çözünürlükle) belirtilmesini öneririm. Bu kapsamlar, son tarama için istatistiklerin nerede belirlendiğini kontrol edecektir.

2. Blok İstatistikleri Gerçekleştirin (Mekansal Analiz Araçları> Mahalle'de bulunur). Hem yükseklik hem de genişlik için 10 hücreli bir dikdörtgen kullanın ve istatistik türü için "ORTALAMA" seçeneğini seçin. İsterseniz farklı şekil ve türleri deneyin. Hücre boyutu aşağı örnekleme oranıdır.

3. Blok istatistikleri raster çözünürlüğünü değiştirmediğinden, son adım Yeniden Örnekleme'dir (Veri Yönetim Araçları> Raster> Raster İşleme bölümünde bulunur). 10 m seçin ve hücrenin merkezinde blok istatistik seçmek için "EN YAKIN" kullanın.

Adım 2 ve 3, Curtvprice'in Toplama aracını kullanma önerisiyle değiştirilebilir , bu da aynı sonuçları dikdörtgen araçlar kullanarak etkili bir şekilde verecektir.

Tüm bu özet değer hücrelerini ve Yeniden Örneklemeyi çoğaltmaktan kaçınan Toplama aracını kullanarak işlemin bir adımını kesebilirsiniz .

Yukarıda belirtilen tüm hususları dikkate almak önemlidir ve bilinear yeniden örneklemenin oldukça sorunlu olduğunu tamamen kabul ediyorum. Yine de neden kimsenin kübik evrişimi tartışmadığını merak ediyorum? Bir blok fonksiyonunun kullanılmasındaki problem, dağılım normal olmayan veya multimodal olduğunda, lidardan türetilen bir DEM ile beklendiği gibi ortalamanın oldukça alakasız olmasıdır.

Orijinal lidar verilerine erişiminiz varsa, ArcGIS'teki "Topo'dan Raster'e" aracını kullanarak verileri istenen çözünürlüğe enterpolasyonlayın. Sadece 1m DEM rasterine erişiminiz varsa, en iyi, en az verimli olsa da, rasterleri noktalara dönüştürmek ve ince bir plaka veya Bi-kübik spline kullanmak gibi bir yöntem olacaktır. Bu, yeniden örnekleme mahallinin verilere doğrusal olmayan bir eğri sığmasını sağlar.

Alternatif olarak, 1m rasterini istediğiniz yeniden örnekleme çözünürlüğünün (10x10) boyutuna yaklaşarak bir Gauss çekirdeği kullanarak düzeltebilirsiniz ve daha sonra bir çift doğrusal yeniden örnekleme çok daha uygun olacaktır. Bu yaklaşım, yumuşatma parametresi üzerinde doğrudan kontrole sahip olmanızı sağlar ve ortalamanın bir eğilimin göstergesi olarak alakalı hale geldiği ve doğrusal bir uyumun desteklendiği bir "yerel" normal dağılımla sonuçlanır.