2010-11-19 1 views
0

내가 가진이R의 좌표를 사용하여 정규 데이터 집합을 어떻게 매핑 할 수 있습니까?

위도 경도 mxgYield

38.96152878 -123.5578995 7.925971605

39.24150085 -123.6392975 19.6720449

39.52162933 -123.7214966 19.777747

39.80192947 -123.8044968 12.23218451

같은 데이터 세트

40.36297989 -123.9729004 -124.058403 1.443492322

38.74369812 -123.1184998 1.938840925

40.64371872 15.00058428 나는 컬러 그라데이션으로 mxgYield을지도하는 것을 시도하고있다. R에서 어떻게해야합니까?

감사

+0

중복 가능 : http://stackoverflow.com/questions/1896419/plotting-a-3d-surface-plot-contour-map-overlay-using-r –

답변

1

의 데이터 볼륨에 조금 약한 것 같습니다하지만 아키 마의 보정시 기능의 조합과 기본 윤곽 플로팅 기능을 사용하여 성공 :

yield <- read.table(textConnection("Lat Lon mxgYield 
38.96152878 -123.5578995 7.925971605 
39.24150085 -123.6392975 19.6720449 
39.52162933 -123.7214966 19.777747 
39.80192947 -123.8044968 12.23218451 
40.36297989 -123.9729004 15.00058428 
40.64371872 -124.058403 1.443492322 
38.74369812 -123.1184998 1.938840925"), header=TRUE) 
require(akima) 
ak.yld <- with(yield, interp(x=Lat, y=Lon, z=mxgYield)) 
contour(ak.yld) 

alt text

가 (에 예뻐 보이는을 내 화면이지만 데이터가있는 지역의 등고선 플롯을 계속 볼 수 있습니다.)

+0

감사합니다. 하지만 내가 원하는 것은 수확량이 다른지도의 그라디언트로 표시된지도입니다. 등고선 플롯이 아닙니다. – didimichael

+0

이제 MAPS 패키지에서 이미지를 사용할 수 있다고 생각합니다. 감사합니다. – didimichael

+0

이미지 (기본 함수)도 interp 객체와 함께 작동하지만 멋지게 보지 않습니다. contour (ak.yld, col = 1 : 10) 또는 image : ak.yld, col = terrain.colors (8)를 사용하여 등고선을 색칠 할 수도 있습니다. –

관련 문제