Xác định mối quan hệ thống kê giữa các raster bằng R vs ArcGIS Desktop?

Tôi đang cố gắng phân tích làm thế nào các bộ raster liên quan với nhau bằng cách sử dụng một số kỹ thuật thống kê. Vì, tôi không có nhiều kinh nghiệm khi sử dụng các công cụ thống kê không gian trong ArcGIS. Tôi đã xuất các trình quét của mình dưới dạng tệp Ascii và phân tích chúng bằng R (cụ thể là maptoolsgói và readAsciiGrid()). Điều này đã hoạt động ổn (nhưng vì các bộ dữ liệu có 90.000 điểm nên chạy phân tích chậm), nhưng tôi không biết liệu tôi có tái tạo trong R , chức năng hiện có trong ArcGIS hay không .

Ví dụ, tôi muốn thực hiện hồi quy giữa mỗi trình quét này bằng một vài phép biến đổi khác nhau (logarit, hàm mũ, v.v.). Điều này có thể được thực hiện trong ArcGIS không? Một câu hỏi rộng hơn thứ hai là nếu có các phương pháp thống kê tiêu chuẩn để kiểm tra loại dữ liệu này?

Mỗi cặp raster có các giá trị dữ liệu / không có dữ liệu phù hợp và tất cả các tham số đều giống hệt nhau, ngoài giá trị Gridcell.

Tôi sẽ bám lấy R. Nếu tốc độ thực sự là một vấn đề (tôi nghi ngờ vì vậy 90.000 không phải là một con số lớn như vậy) bạn có thể thử tìm mối quan hệ giữa một tập hợp con dữ liệu của bạn. Trên thực tế, điều đầu tiên tôi sẽ làm là lập một âm mưu để tìm kiếm các mối quan hệ rõ ràng.

Ngay cả khi arcgis chứa các công cụ để so sánh các trình quét, R sẽ luôn cung cấp cho bạn nhiều công cụ thống kê hơn.

Ví dụ:

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)
...

Tôi thực sự nên thêm rằng thường thì sẽ chính xác hơn khi làm việc với một tập hợp con dữ liệu của bạn sau đó với bộ dữ liệu đầy đủ của bạn. Trong nhiều trường hợp, các ô lưới không độc lập với các ô dữ liệu xung quanh, điều này sẽ dẫn đến các giá trị p quá lạc quan, ví dụ như hồi quy phù hợp (bạn sẽ tìm thấy thêm thông tin nếu bạn tìm kiếm trên bộ giải mã).

Nhìn vào gói R raster , nó được thiết kế với loại vấn đề cụ thể này. Nó cố gắng giữ ít raster trong bộ nhớ và thực hiện một số thao tác không gian cơ bản - thông qua GDAL, nó có thể xử lý các tệp ESRI Grid gốc, trong số nhiều định dạng khác. Các họa tiết cung cấp một số ví dụ tốt đẹp về việc sử dụng nó.

Tôi vừa đọc về gói ArcGIS 10 là plugin R. Tôi chưa có cơ hội khám phá nó đầy đủ, nhưng có thể sửa đổi điều này để làm những gì tôi mô tả ở trên.