Làm cách nào để ước tính mật độ hạt nhân với hàng rào / ranh giới vật lý trong QGIS?


18

Có ai biết về một chương trình sẽ cho phép kernelling tính đến một ranh giới trong tính toán của nó chứ không chỉ đơn giản là che giấu các khu vực không thể?

Cho đến nay tôi đã tìm thấy: GME add on Arc10 (Hawthorne Beyer) - Tôi gặp lỗi mỗi khi tôi chỉ định một shapefile ranh giới. Tôi đã thử rất nhiều loại shapefile khác nhau và độ phức tạp ranh giới, v.v ... Nó hoạt động tốt khi tôi không chỉ định một ranh giới.

Gói AdehabitatHR trong R (Calange 2011) - điều này hoạt động tốt nhưng ranh giới bạn chỉ định phải rất đơn giản - băng thông hạt nhân 3x có độ dài và không quá quanh co. Đối với dữ liệu của tôi, đây là một sự đơn giản hóa lớn.

Vì vậy, tôi tự hỏi liệu có phần mềm nào khác có thể làm điều đó không, ví dụ GRASS hoặc QGIS.

Cảm ơn


3
Bài báo gần đây với mã R quan tâm tiềm năng, Biên giới biên và hạt nhân có trọng số . Tất cả các mã chỉ là cách để cân nhắc các quan sát để điều chỉnh độ lệch biên (nếu bạn có thể tìm ra các trọng số, bạn có thể sử dụng bất kỳ chương trình nào để ước tính kde miễn là nó sẽ trừ trọng số).
Andy W

1
Chính xác làm thế nào để bạn muốn nó "tài khoản" cho một ranh giới? Có nhiều cách có thể, từ che giấu đến chặn lây lan đến sửa lỗi cho các hiệu ứng biên.
whuber

1
Cảm ơn bạn @AndyW và whuber đã trả lời. Tôi muốn giải thích cho ranh giới bằng cách ngăn chặn kernel lan rộng ra nó. Dữ liệu là các vị trí rùa biển, vì vậy tôi biết rằng phân phối sử dụng không nên lan rộng trên đất liền, nhưng nhiều vị trí rất gần bờ biển nên tôi không muốn thiên vị với các vị trí này bằng cách che giấu các phần nhân lan vào đất liền. Cảm ơn bạn đã liên kết đến bài báo - điều này có vẻ tuyệt vời tôi sẽ cố gắng sử dụng mã tương tự cho dữ liệu của tôi.
KimS

Câu trả lời:


2

Tại cốt lõi của ước tính mật độ hạt nhân là khái niệm khoảng cách. Giải pháp tốt nhất mà tôi biết là sử dụng thước đo khoảng cách tốt hơn chiếm các ranh giới và chi phí đi lại khác nhau. Tốt nhất là chọn một thước đo khoảng cách phù hợp với vấn đề bạn đang cố gắng giải quyết. Ma sát địa hình trên đất liền là tuyệt vời để đi bộ đường dài, nhưng không phải để phân tán khí dung. Các luồng gió rất cần thiết cho việc theo dõi thuyền buồm, nhưng không liên quan đến các hướng lái xe.
Bây giờ, khái niệm về số liệu khoảng cách phù hợp hy vọng đủ động lực, tôi có thể đề xuất các bề mặt chi phí như một thước đo khoảng cách mục đích chung tốt. Chúng có sẵn trong mọi thứ, từ ArcGIS, đến R, đến JavaScript và khá đơn giản để xây dựng. Ví dụ, trong qGIS, bạn có thể xây dựng bề mặt ma sát raster và sử dụng bề mặt ma sát đó. Tùy chỉnh bề mặt ma sát để giải thích cho ranh giới của bạn và bạn sẽ thấy khối lượng từ các hạt nhân xung quanh các điểm của bạn nằm gọn gàng xung quanh các chướng ngại vật.


+1 Đây có khả năng là một cách tiếp cận tuyệt vời. Đã áp dụng nó trong nhiều vấn đề, tôi thấy rằng khó khăn chính là một vấn đề không được đề cập ở đây: chính xác làm thế nào để bạn đề xuất việc truyền bá hạt nhân theo số liệu này? Các giải pháp vượt trội không thể làm được điều này.
whuber

@whuber Đúng, các hệ thống bên ngoài thường không có khả năng tính toán cho cả bề mặt chi phí và mật độ hạt nhân cùng một lúc. Đề nghị của tôi sẽ là sử dụng bề mặt ma sát để tính khoảng cách điểm-điểm, sau đó áp dụng trực tiếp một hàm nhân cho trọng lượng. Trong qGIS hoặc ArcGIS, bạn có thể viết sử dụng hàm phân phối bình thường trong python và trong R có dnorm()hàm này.
JasonRDalton

Tôi đã thực hiện chính xác điều đó trong một số trường hợp, nhưng suy nghĩ về quá trình và kiểm tra kết quả cho thấy một số vấn đề cố hữu. Cái khó nhất dường như là làm thế nào để đảm bảo khối lượng đó được bảo toàn: bạn không thể chỉ trải các giá trị theo hàm Gaussian cố định, bởi vì tổng khối lượng kết quả của một điểm trải rộng sẽ không bằng khối lượng ban đầu. Không có cách hiệu quả để thực hiện tính toán này.
whuber



-2

Tôi đã thực hiện thành công việc này bằng cách sử dụng Phân tích địa lý, một phần mở rộng của ESRI ArcGIS. Bạn có thể tải tập dữ liệu tính năng dòng làm rào cản và kết quả khá tốt. Tôi thực hiện chạy thử trong khi thay đổi các tham số của các chức năng để có ý tưởng rõ ràng về cách hiệu chỉnh công cụ. Để xác thực kết quả, điều tôi khuyên là nếu bạn có bộ dữ liệu đủ lớn, hãy xóa một mẫu điểm, tạo bề mặt mật độ không có các điểm đó, sau đó so sánh sự khác biệt giữa các giá trị của các điểm bị loại bỏ với các giá trị của bề mặt tại vị trí của các điểm bị loại bỏ.


2
Bạn có thể khuếch đại câu trả lời này một chút để giải thích cách GA - thực hiện các phương pháp nội suy - có thể tính toán ước tính mật độ hạt nhân không? Khi kiểm tra kỹ tính chính xác của giải pháp, bạn có xác minh rằng tích phân của lưới GA xấp xỉ bằng tổng của tất cả dữ liệu đầu vào không?
whuber

1
Trong 10.1 có thể thực hiện phép nội suy hạt nhân w / rào cản nhưng điều này không dẫn đến ước tính mật độ. Mặc dù, được đưa ra các phương thức mà tài liệu tham khảo bài gốc ban đầu, tôi tự hỏi nếu chúng không quan tâm đến ước tính nhân Gaussian và không phải là ước tính mật độ.
Jeffrey Evans
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.