QGIS에서 물리적 장벽 / 경계로 커널 밀도 추정을 수행하는 방법은 무엇입니까?


18

누구나 불가능한 영역을 마스킹하는 대신 경계를 계산하는 커널 링을 허용하는 프로그램을 알고 있습니까?

지금까지 내가 발견했습니다 : Arc10 (Hawthorne Beyer) 용 GME 애드온-경계 모양 파일을 지정할 때마다 오류가 발생합니다. 다양한 shapefile 유형과 경계 복잡도 등을 많이 시도했습니다. 경계를 지정하지 않으면 정상적으로 작동합니다.

R의 AdehabitatHR 패키지 (Calange 2011)-이 방법은 잘 작동하지만 지정한 경계는 매우 단순해야합니다. 즉, 3 배의 커널 대역폭 길이와 너무 비 틀리지 않아야합니다. 내 데이터의 경우 이것은 지나치게 단순화 된 것입니다.

따라서 다른 소프트웨어가 GRASS 또는 QGIS와 같이 할 수 있는지 궁금합니다.

감사


3
잠재적 관심 R 코드, 경계 편향 및 가중치 커널이 포함 된 최근 논문 . 모든 코드는 경계 바이어스를 수정하기 위해 관측 값 에 가중치를 부여 하는 방법 일뿐입니다 . 가중치를 계산할 수있는 경우 가중치를 제외하고 kde를 추정하기 위해 모든 프로그램을 사용할 수 있습니다.
Andy W

1
정확하게 경계를 어떻게 설명 하고 싶습니까? 마스킹에서 확산 차단, 경계 효과 보정에 이르기까지 다양한 방법이 있습니다.
whuber

1
@AndyW와 whuber에게 감사합니다. 커널이 확산되지 않도록 경계를 설명하고 싶습니다. 데이터는 바다 거북 위치이므로 활용 분포가 육지로 퍼져서는 안된다는 것을 알고 있지만 많은 위치가 해안 근처에 있으므로 커널 부분을 단순히 가리는 것만으로 이러한 위치에 치우치고 싶지 않습니다. 땅에 퍼졌다. 논문에 대한 링크를 주셔서 감사합니다-이것은 내 데이터에 비슷한 코드를 사용하려고 시도하는 것이 좋습니다.
KimS

답변:


2

커널 밀도 추정의 핵심은 거리 개념입니다. 내가 아는 가장 좋은 솔루션은 경계와 다양한 여행 비용을 설명하는 더 나은 거리 측정법을 사용하는 것입니다. 해결하려는 문제에 맞는 거리 측정법을 선택하는 것이 가장 좋습니다. 육상 지형 마찰은 하이킹에는 좋지만 에어로졸 분산에는 적합하지 않습니다. 범선 추적에는 바람 흐름이 필수적이지만 운전 방향에는 관련이 없습니다.
적절한 거리 측정법의 개념이 충분히 동기 부여되기를 원하므로 이제 비용 표면을 좋은 범용 거리 측정법으로 추천 할 수 있습니다. ArcGIS에서 R, JavaScript에 이르기까지 모든 것이 가능하며 구성이 매우 간단합니다. 예를 들어 qGIS에서 래스터 마찰 표면을 구성하고이를 사용하여 경로를 계산할 수 있습니다. 경계를 고려하여 마찰 표면을 사용자 정의하면 지점 주변의 커널 질량이 장애물 주위에 깔끔하게 퍼져있는 것을 볼 수 있습니다.


+1 이것은 잠재적으로 훌륭한 접근법입니다. 많은 문제에 적용한 결과, 여기서 가장 어려운 점은 여기서 언급되지 않은 것입니다. 정확히이 메트릭에 따라 커널을 확산시키는 방법은 무엇입니까? 즉시 사용 가능한 솔루션으로는이 작업을 수행 할 수 없습니다.
whuber

@whuber True, 즉시 사용 가능한 시스템은 일반적으로 비용 표면과 커널 밀도를 동시에 설명 할 수있는 기능이 없습니다. 내 제안은 마찰면을 사용하여 점대 점 거리를 계산 한 다음 가중치에 커널 함수를 직접 적용하는 것입니다. qGIS 또는 ArcGIS에서는 파이썬에서 정규 분포 함수를 사용하여 작성할 수 있으며 R에는 dnorm()함수가 있습니다.
JasonRDalton

어떤 경우에는 정확히 그렇게했지만 프로세스를 생각하고 결과를 조사하면 몇 가지 고유 한 문제가 드러납니다. 가장 어려운 방법은 질량 보존을 보장하는 방법 인 것 같습니다. 분산 된 점의 총 결과 질량이 원래 질량과 같지 않기 때문에 고정 가우스 함수에 따라 값을 분산시킬 수는 없습니다. 이 계산을 수행하는 효율적인 방법은 없습니다.
whuber



-2

ESRI ArcGIS의 확장 인 Geostatistical Analyst를 사용하여 성공적으로 수행했습니다. 선 피쳐 데이터 세트를 장벽으로로드 할 수 있으며 결과는 꽤 괜찮습니다. 도구를 보정하는 방법에 대한 명확한 아이디어를 얻기 위해 함수의 매개 변수를 변경하는 동안 테스트 실행을 수행합니다. 결과의 유효성을 검사하려면 데이터 집합이 충분히 큰 경우 점 샘플을 제거하고 해당 점없이 밀도 표면을 생성 한 다음 제거 된 점의 값과 표면의 값 간의 차이를 비교하는 것이 좋습니다. 제거 된 지점의 위치


2
보간 방법 을 구현하는 GA 가 커널 밀도 추정값 을 계산할 수 있는 방법을 설명하기 위해이 답변을 약간 증폭시킬 수 있습니까? 솔루션의 정확성을 다시 확인하면서 GA 그리드의 적분이 모든 입력 데이터의 합과 거의 같은지 확인 했습니까?
whuber

1
10.1에서는 장벽으로 커널 보간을 수행 할 수 있지만 밀도 추정은되지 않습니다. 비록 원래 게시물이 참조하는 방법을 고려할 때 밀도 추정값이 아닌 가우스 커널 추정값에 관심이 없는지 궁금합니다.
Jeffrey Evans
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.