ArcGIS Desktop'ta rastlantısal olmayan iki eşitlikli tarama ızgaraları hizalanıyor mu?

İki rasterim var: bir dijital yükseklik modeli (DEM) ve bir Landsat bandı. Her ikisinin ızgara çözünürlüğü 30 m'dir, ancak ızgaralar çakışmaz (resme bakın: üstteki karanlık tarama DEM'dir; altındaki ışık tarama Landsat bandıdır).

Şebekenin DEM şebekesi ile mükemmel bir şekilde hizalanması için Landsat bandında bir dönüşüm yapmak istiyorum. Ortam ayarlarında bir Snap Raster olarak DEM kullanarak ArcGIS 30 m Landsat band yeniden örnekleme denedim, ama bu herhangi bir etkisi yoktur. Ayrıca DEM Snap Raster olarak ayarlarken, boşuna değil, bir vektör çokgen (DEM de kırpılır) ölçüde Landsat bant kırpma denedim.

Bunun olduğundan şüpheleniyorum çünkü her iki raster zaten aynı çözünürlükte (sadece farklı ızgaralarda), ancak daha fazla analiz yapabilmek için ızgaralarının hizalanmasına gerçekten ihtiyacım var.

Bunu ArcGIS Desktop'ta yapabilir miyim?

Daha fazla araştırma sonra ben arızalı çekin Raster ayarı aslında çözüm kapatmak için gibi görünüyor ArcGIS 10. bir hata neden olduğunu keşfetti Geoprocessing (Arkaplan Geoprocessing \ Geoprocessing Seçenekler ... sonra işaretini Enable içinde Arkaplan İşleme bölümünde), ArcCatalog'da Snap Raster belirtilmiş olarak kırpma işlemini gerçekleştirin.

ArcGIS forumlarında burada bir tartışma var: http://forums.arcgis.com/threads/50808-Another-Snap-to-Raster-problem

Görünüşe göre hata ArcGIS SP 4'te düzeltildi. Hizmet paketini yükledim ancak henüz test etmedim.

İlginçtir (ya da daha doğrusu, can sıkıcı bir şekilde), klibi ilgi alanı poligonuma yukarıdaki forum gönderisindeki çözüme göre belirtilen Snap Raster (Snap Raster DEM'dir) ile uygularsam, Snap Raster çalışır ancak Landsat bandın, DEM rasterinden bir sütun ve bir satır daha vardır, bu da ArcGIS'de ilgi alanı çokgene kırpılmıştır. Bu sinir bozucu çünkü daha fazla analiz için iki rasterin aynı sayıda satıra ve sütuna sahip olmasını istiyorum. Basit bir çözüm, Con()DEM ile aynı sayıda satır ve sütuna sahip yeni bir raster oluşturan ancak Landsat bandının değerlerine sahip bir hesaplama yapmaktır .

Aslında bu sorunla uzun süre mücadele ettim ve çok etkili ve oldukça kolay bir çözüm olduğunu düşündüğüm şeyi buldum. "Doğru" yerde olduğuna inandığınız rasterlerinizden birini seçin, ArcGIS'in izin verdiği ölçüde yakınlaştırın ve doğrudan dört raster hücresinin kesişiminin üstünde yeni bir nokta şekil dosyası oluşturun. Buradaki bir numara, mümkün olduğunca yakınlaştırdığınızdan ve yeni noktayı doğrudan dört tarama hücresinin oluşturduğu 'haç'ın üzerinde oluşturduğunuzdan emin olmaktır. Sonra yeni oluşturduğunuz noktanın koordinatlarını kaydedin. Son olarak, ikinci tarama cihazınızı yeniden projeksiyon yapın (ancak projeksiyon ve çıkış hücre boyutunu değiştirmeden bırakın) ve yeni oluşturduğunuz noktanın koordinatlarını "Referans Noktası" olarak kullanın. Bu, iki rasterinizi mükemmel şekilde hizalayacaktır. Hata yok,

DEM'yi Landsat ile kaplamak için maskeleme aracını kullanmayı deneyebilirsiniz.

Tamamen üst üste düşen pikselleri yakalamak için bir çalışma alanı kutusu oluşturmanız gerekebilir.

Bir tarama ile diğeri arasındaki boşluğu ölçün ve ardından bir tarama işlemini uygun şekilde kaydırın.

Vardiya aracını veri yönetimi araç kutusundan kullanabilirsiniz.

Ben de aynı problemi yaşadım. Küçük bir test örneği üzerinde harika çalışan Louis Longchamps çözümünü denedim , ancak veri kümem bu yöntemi pratik olarak kullanmak için çok büyüktü. Bir süre sonra ve bunu yapmaya çalışmanın birçok yolu, nihayet, bir kapris üzerinde, bir klipten sonra "Dikdörtgen ile Ayıkla" kullandım, bu da klips aracındaki yapıştırma raster işlevinin çalışması gerektiği gibi çalıştı. Bir model oluşturucu kurdum ve kullanmaya çalıştım, ama bana bir hata verdi, bu yüzden bunu her seferinde bir adım yapmak zorunda kaldım, ama sonunda tüm rasterlerim eşleşti.

Rasteri değiştirme fikrinden rahatsız değildim çünkü bana göre bilgileri değiştiriyorsunuz. Bundan kaçınmak için bir dizi dönüşüm kullandım.

İşte adımlar:

1. Şablon olarak kullanmak istediğiniz tarama için, içeri girin Properties/Sourceve Cellsizetüm ondalık sayıları içeren bir kağıda yazın.
2. Bu raster kullanarak puanlara dönüştürün Conversion Tools/Raster to Point.
3. Spatial Analyst Tools/Extraction/Extract Multi Values to PointsGiriş noktası özelliği olarak az önce oluşturduğunuz Nokta özelliğini (adım 2'de) ve Giriş sıralayıcılarında örneklemek istediğiniz tüm rasterleri belirterek, kullanmak istediğiniz diğer tüm rasterlerden değerleri çıkarın .
4. Conversion Tools/To Raster/Point to RasterDeğer alanında oluşturduğunuz katmanı ve karşılık gelen bir dosya adını VE Cellsizekağıda yazdığınız değeri kullanarak bu noktalardaki rasterleri tek tek oluşturun .

Mükemmel bir şekilde hizalanmaları gerekir ve bunları örneğin R'ye aktarabilir ve tam olarak aynı boyutlarda matris ve vektörlere sahip olabilirsiniz.

DEM imajımızla ilgili çok endişeliydik, aynı zamanda denetimli sınıflandırılmış raster imajına da hizalanmadı. "Matthew Long" yöntemini de kullandık .

hizalanmamış ve hizalanmış görüntülerin anlık görüntülerini ekledik.

Matthew Long'un yönteminden sonra DEM imajımız, denetlenen sınıflandırılmış raster imajımızla tamamen hizalandı.

Bu, Veri Yönetimi araçlarındaki ARCGIS 10.2.1'de DEM imajımızı hizalamak için kullandığımız aracın görüntüsüdür:

Bu, hizalamadan önceki görüntü, siyah görüntü DEM ve renkli görüntü denetimli gizli tarama görüntüsüdür.

Bu, hizalamadan sonraki görüntüdür, siyah görüntü DEM'dir ve renkli görüntü denetimli sınıflandırılmış tarama görüntüsüdür.

Sorunumu "Matthew Long" yöntemiyle çözüyorum. Ancak, daha da kolay ve% 100 doğruluk için. Herhangi bir yakınlaştırma olmadan, iyi hizalı rasterın sol üst genişlik koordinatını kullanın.

Bunu yapmak için: Extent özelliği altında sağ tıklayın, özellik, Kaynak paneli

Ayrıca, sol üst taraf koordinatını kullanarak, manuel müdahaleler olmadan bir python çözümü oluşturmak mümkündür.

Farklı kaynaklardan iki dosyamda benzer bir sorun yaşadım. İlk olarak, ülke düzeyinde (Peru) ormansızlaşma referans seviyesi için bir raster dosyası. İkinci olarak, GEE platformunda birleştirilmiş bir dosya indirildi ve ilk eleman sistemine yeniden yansıtıldı. Sonuç olarak, her ikisinin çakışması, hücreler (30 x 30 hücre) arasında 5 metreden fazla olmayan bir yer değiştirme olduğunu gösterdi. İkinci dosyayı birinciyle hizalamak ve ayarlamak için, her iki dosyayı * kapsayabilen * genişletilmiş bir ilk dosya tabanı oluşturmak ve ilk dosyanın projeksiyon sistemi ve diğer özellikleri ile piksel hizalamasını korumak benim başıma geldi . Böylece, ikinci dosyamı veya benzer bir kayma ile herhangi bir yan ürünü yeniden projelendirirken,