2011-02-03 4 views
8

를 기반으로하지만 내 데이터 파일은 기본적으로 내가 위도 긴 값을 보간하고 zvalue에 따라 색상을 그릴 필요이은 R의 히트 맵이 된 우편 번호 나 R의 히트 맵을 그릴하려는 경우에만

Lat, Long, Zip, Zvalue 

같다으로 한 도면

어떻게 R

에 나는 결국이

enter image description here

같은 것을 얻고 싶은 것을 할 수있는
+1

'raster' 패키지를 사용해보십시오. –

+0

위도와 경도 데이터가있는 경우 관련성이없는 우편 번호입니까? –

+0

예, 우편 번호 또는 위도/경도 데이터를 사용할 수 있습니다. – Mark

답변

3

The pa ckage spatstat은 (는) 친구입니다! 데이터가 이벤트 집합을 경우

http://www.spatstat.org/spatstat/

, 다음 생성하기 위해 커널 밀도 추정을 사용할 수 있습니다 (같은은 X, Y)이 (많은 "범죄는 (x, y)에 발생"참조) 히트 맵. 당신은 여기에서 예를 볼 수 있습니다

https://github.com/drewconway/WikiLeaks_Analysis/blob/master/density.r

를 특별히 선 다음 사용하는 spatstat의 smooth 기능을 사용할 수 있습니다

당신의 Z 값은 (아마도 귀하의 질문에 읽기의 경우) 실제 값은 경우 (72) 보간을 수행하고 데이터를 보간하여 생성 된 픽셀 집합을 되 찾는 가우스 커널.

또는 akima 패키지 (spatstat 작성자가 권장)를 사용하여 지정된 위치에서 보간 할 수 있습니다. 선형 또는 스플라인 기반 보간 중 하나를 사용하며 매우 직관적 인 것처럼 보입니다 (경험이 없지만!).

더 일반적으로 당신이하려고하는 것은 "Kriging"이라고 불리는 것이므로, Google에 해당하는 경우 많은 결과가 발생합니다. 예 : gstatgeoR 패키지

마지막으로 (및 FTW) Gaussian Process를 사용하여 동일한 종류의 작업을 수행 할 수 있습니다. 이것은 실제로 데이터에 주어진 보간에 대한 분포를 제공합니다. R 패키지 kernlab에는 구현 방법이 있지만 사용 방법은 모르겠습니다.

3

제프리 브린은 ggplot2와 우편 번호 패키지를 사용하여 플롯 이러한 유형의 수행했습니다 http://jeffreybreen.wordpress.com/2011/01/05/cran-zipcode/

당신이

을 시작하는 것을
+0

하지만 히트 맵조차도 아니며 최고로 Choropleth지도입니다. 자동으로 보간을 할 수있는 함수/패키지가 있는지 궁금합니다. – Mark

+0

@Mark : 아마도'heat.colors'? – Aaron

2
당신은 ggplot 케이스의 코드를 볼 수 있습니다

는 연구 경쟁의 승자 here

geom_tile을 사용하여 heatmap을 만들었습니다. 그래프의 제목은 입니다. Houston Houston의 폭력 범죄 날씨지도