Menentukan hubungan statistik antara raster menggunakan R vs ArcGIS Desktop?

8

Saya mencoba untuk menganalisis bagaimana set raster berhubungan satu sama lain menggunakan beberapa teknik statistik. Karena, saya tidak punya banyak pengalaman menggunakan alat statistik spasial di ArcGIS saya mengekspor raster saya sebagai file Ascii , dan menganalisisnya menggunakan R (khusus maptoolspaket, dan readAsciiGrid()). Ini telah berfungsi dengan baik (tetapi karena dataset memiliki 90.000 poin, lambat untuk menjalankan analisis), tetapi saya tidak tahu apakah saya membuat ulang di R , fungsionalitas yang ada di ArcGIS .

Sebagai contoh, saya ingin melakukan regresi antara masing-masing raster ini menggunakan beberapa transformasi yang berbeda (logaritmik, eksponensial, dll). Bisakah ini dilakukan dalam ArcGIS ? Pertanyaan kedua yang lebih luas adalah apakah ada metode statistik standar untuk memeriksa data jenis ini?

Setiap pasangan raster memiliki nilai data yang cocok / tidak ada data dan semua parameter identik, selain dari nilai gridcell.

djq
sumber

Jawaban:

9

Saya akan tetap berpegang pada R. Jika kecepatan benar-benar masalah (saya ragu jadi 90.000 bukan angka yang besar), Anda dapat mencoba menemukan hubungan antara subset data Anda. Sebenarnya hal pertama yang akan saya lakukan adalah membuat plot untuk mencari hubungan yang jelas.

Bahkan jika arcgis berisi alat untuk membandingkan raster, R akan selalu memberi Anda lebih banyak alat statistik.

Misalnya:

library(rgdal)
map1<-readGDAL('file.asc')
map2<-readGDAL('file2.asc')
samplenr<-sample(length(map1$band1), 1000)
smallset<-data.frame(map1=map1$band1[samplenr],map2=map2$band1[samplenr])
plot(smallset)
lm(map2~map1, smallset)
...

Saya harus benar-benar menambahkan bahwa sering kali lebih tepat untuk bekerja dengan subset data Anda kemudian dengan dataset lengkap Anda. Dalam banyak kasus sel-sel kisi tidak independen dari sel-sel data di sekitarnya, yang akan menghasilkan nilai p yang terlalu optimis untuk misalnya kecocokan regresi (Anda akan menemukan lebih banyak info jika Anda mencari pada declustering).

johanvdw
sumber
2
+1. Baik ArcGIS maupun, dalam hal ini, GIS apa pun tidak akan pernah menyediakan kemampuan seperti R atau paket statistik berfitur lengkap lainnya. Akan bodoh bagi mereka untuk mencoba. Apa yang dapat kita harapkan - dan ArcGIS tampaknya berada pada tahap awal - adalah untuk platform GIS untuk menyediakan integrasi yang efisien dengan aplikasi lain untuk analisis statistik dan visualisasi.
whuber
Itu benar - kecepatan sebenarnya bukan masalah. Dibutuhkan sekitar 10 detik untuk memuat dua raster dan 10-15 detik per operasi sejauh ini. Jika saya menjaga hal-hal dalam perspektif, saya sudah menunggu di ArcGIS selama hal-hal terjadi! Saya berencana untuk menggunakan sampling lebih banyak, meskipun yang akan menghapus sedikit keterlambatan.
djq
4

Lihatlah ke dalam paket Raster , itu dirancang dengan jenis masalah tertentu dalam pikiran. Ia mencoba untuk menyimpan sedikit dari raster dalam memori, dan melakukan sejumlah operasi spasial dasar - melalui GDAL ia dapat menangani file ESRI Grid asli, di antara banyak format lainnya. Sketsa memberikan beberapa contoh bagus penggunaannya.

scw
sumber
1

Saya baru saja membaca tentang paket ArcGIS 10 yang merupakan plugin R. Saya belum memiliki kesempatan untuk menjelajahinya sepenuhnya, tetapi dimungkinkan untuk memodifikasi ini untuk melakukan apa yang saya jelaskan di atas.

djq
sumber