Menyelaraskan dua kisi-kisi raster resolusi-tidak-kebetulan di ArcGIS Desktop?

Saya punya dua raster: model elevasi digital (DEM) dan band Landsat. Resolusi grid keduanya adalah 30 m, tetapi grid bukan kebetulan (lihat gambar: raster gelap di atas adalah DEM; raster cahaya di bawahnya adalah band Landsat).

Saya ingin melakukan beberapa transformasi band Landsat sehingga gridnya sejajar sempurna dengan grid DEM. Saya sudah mencoba resampling band Landsat ke 30 m di ArcGIS menggunakan DEM sebagai Snap Raster dalam pengaturan Lingkungan, tetapi ini tidak memiliki efek apa pun. Saya juga telah mencoba memotong band Landsat ke tingkat vektor poligon (bahwa DEM juga terpotong), sambil mengatur DEM sebagai Snap Raster, tetapi tidak berhasil.

Saya menduga ini terjadi karena kedua raster sudah pada resolusi yang sama (hanya pada grid yang berbeda), tetapi saya benar-benar membutuhkan grid mereka untuk disejajarkan agar dapat melakukan analisis lebih lanjut.

Bisakah saya melakukan ini di ArcGIS Desktop?

Setelah penelitian lebih lanjut, saya menemukan bahwa pengaturan Snap Raster yang tidak berfungsi sebenarnya disebabkan oleh bug di ArcGIS 10. Solusinya adalah mematikan Background Geoprocessing ( Geoprocessing \ Geoprocessing Options ... lalu hapus centang Aktifkan di bagian Background Processing ), dan melakukan operasi kliping, dengan Snap Raster yang ditentukan, di ArcCatalog.

Ada diskusi di forum ArcGIS di sini: http://forums.arcgis.com/threads/50808-Another-Snap-to-Raster-problem

Rupanya bug diperbaiki di ArcGIS SP 4. Saya sudah menginstal paket layanan tetapi belum diuji.

Menariknya (atau lebih tepatnya, menyebalkan), jika saya melakukan klip ke poligon bidang-minat saya dengan Snap Raster yang ditentukan (Snap Raster adalah DEM) sesuai solusi di posting forum, di atas, Snap Raster berfungsi tetapi Landsat band memiliki satu kolom lebih banyak dan satu baris lebih banyak daripada raster DEM, yang juga dipotong di ArcGIS ke poligon area-of-interest. Ini menjengkelkan karena untuk analisis lebih lanjut saya perlu dua raster memiliki jumlah baris dan kolom yang persis sama. Solusi sederhana adalah dengan melakukan perhitungan dengan Con()fungsi yang membuat raster baru dengan jumlah baris dan kolom yang sama dengan DEM, tetapi yang memiliki nilai band Landsat.

Saya sebenarnya berjuang dengan masalah ini untuk waktu yang lama, dan saya menemukan apa yang saya anggap sebagai solusi yang sangat efektif dan cukup mudah. Pilih salah satu raster Anda, salah satu yang Anda yakini berada di tempat yang "benar", perbesar sejauh yang diperbolehkan ArcGIS, dan buat titik baru shapefile langsung di atas persimpangan empat sel raster. Satu trik di sini adalah memastikan bahwa Anda diperbesar sejauh mungkin, dan bahwa Anda membuat titik baru langsung di atas 'salib' yang dibentuk oleh empat sel raster. Kemudian, rekam koordinat titik yang baru saja Anda buat. Akhirnya, proyeksi ulang raster kedua Anda (tetapi biarkan ukuran sel proyeksi dan output tidak berubah), dan gunakan koordinat titik yang baru Anda buat sebagai "Titik Referensi". Ini akan dengan sempurna menyesuaikan dua raster Anda. Tidak ada kesalahan,

Anda bisa mencoba menggunakan alat masking untuk menutupi DEM dengan Landsat.

Anda mungkin harus membuat kotak area studi untuk menangkap piksel yang benar-benar jatuh di atas satu sama lain.

Ukur jarak antara satu raster dan yang lain, dan kemudian geser satu raster sesuai.

Anda dapat menggunakan alat shift dari kotak alat manajemen data.

Saya memiliki masalah yang sama. Saya mencoba solusi Louis Longchamps , yang bekerja sangat baik pada sampel uji kecil, tetapi dataset saya terlalu besar untuk menggunakan metode ini secara praktis. Setelah beberapa saat dan banyak cara untuk mencoba melakukan ini, saya akhirnya, pada kemauan, menggunakan "Ekstrak oleh Rectangle" setelah resample, yang berfungsi seperti fungsi snap raster pada clip tool seharusnya bekerja. Saya mengatur dan mencoba menggunakan pembuat model, tetapi itu memberi saya kesalahan, jadi akhirnya harus melakukan langkah ini satu per satu, tetapi akhirnya semua raster saya cocok.

Saya merasa tidak nyaman dengan ide menggeser raster karena, menurut saya, Anda mengubah informasinya. Untuk menghindarinya, saya menggunakan urutan konversi.

Berikut langkah-langkahnya:

1. Untuk raster yang ingin Anda gunakan sebagai templat, masuk Properties/Sourcedan tulis Cellsizedi selembar kertas dengan semua desimal tersedia.
2. Ubah raster ini menjadi poin menggunakan Conversion Tools/Raster to Point.
3. Ekstrak nilai dari semua raster lain yang Anda tertarik gunakan Spatial Analyst Tools/Extraction/Extract Multi Values to Points, tentukan fitur Point yang baru saja Anda buat (pada langkah 2) sebagai fitur titik input dan semua raster yang ingin Anda sampel dalam raster Input.
4. Buat raster dari titik-titik ini satu per satu menggunakan Conversion Tools/To Raster/Point to Raster, tentukan layer mana yang Anda buat di bidang Nilai dan menggunakan nama file yang sesuai DAN Cellsizenilai yang Anda tulis di selembar kertas.

Mereka harus menyelaraskan dengan sempurna dan Anda dapat mentransfer ini ke R misalnya dan memiliki matriks dan vektor dengan dimensi yang sama persis.

Kami sangat khawatir tentang gambar DEM kami, itu juga tidak selaras dengan gambar raster rahasia yang diawasi. Kami juga telah menggunakan metode "Matthew Long" .

kami telah menambahkan snapshot dari gambar yang tidak selaras dan selaras.

Setelah metode oleh Matthew Long gambar DEM kami benar-benar selaras dengan gambar raster rahasia yang diawasi.

Ini adalah gambar alat yang kami gunakan untuk menyelaraskan gambar DEM kami, dalam ARCGIS 10.2.1 dalam alat Manajemen Data:

Ini adalah gambar sebelum menyelaraskan, gambar hitam adalah DEM dan gambar berwarna diawasi gambar raster diklasifikasikan.

Ini adalah gambar setelah disejajarkan, gambar hitam adalah DEM dan gambar berwarna diawasi gambar raster diklasifikasikan.

Saya memecahkan masalah saya dengan metode "Matthew Long". Namun, untuk membuatnya lebih mudah dan akurasi 100%. Tanpa zooming apa pun, cukup gunakan koordinat Kiri atas dari raster yang selaras.

Untuk melakukannya: Klik kanan, properti, panel Sumber di bawah properti Extent

Juga, menggunakan koordinat batas kiri atas, dimungkinkan untuk membuat solusi python tanpa intervensi manual.

Saya memiliki masalah serupa memiliki dua file dari sumber yang berbeda. Pertama, file raster untuk tingkat referensi deforestasi di tingkat negara (Peru). Kedua, file konsolidasi pada platform GEE diunduh dan diproyeksikan ulang ke sistem elemen pertama. Sebagai kesimpulan, tumpang tindih keduanya menunjukkan kepada saya bahwa ada perpindahan antar sel (30 * 30 sel) tidak lebih dari 5 meter. Untuk menyelaraskan dan menyesuaikan file kedua dengan yang pertama, terpikir oleh saya untuk membuat * basis yang diperluas dari file pertama yang dapat mencakup kedua file * dan mempertahankan penyelarasan piksel dari file pertama dengan sistem proyeksi dan fitur-fitur lainnya. . Jadi, pada saat memproyeksi ulang file kedua saya atau produk sampingan dengan offset yang sama,