최근에 비슷한지도를 작성해야했습니다. 내가 생각해 낸 솔루션은 ggplot2 (RapidVis 패키지) 대신 rasterVis 패키지를 사용합니다.
제 경우에는 시간 경과에 따른 트렌드 맵 (수치 변수이기도 함)과 트렌드의 중요성을 보여주는 맵 (통계 테스트에서 얻은)이있었습니다.
필자의 경우 신뢰도 수준이 더 높은 곳을 보여주기 위해 숫자 맵 위에 유의성 맵을 플롯하여 "비틀림"패턴과 비슷했습니다.
이것은 대략지도를 만드는 데 사용한 코드입니다.
require(raster)
require(rasterVis)
# Scratch raster objects
data(volcano)
r1 <- raster(volcano)
over <- ifelse(volcano >=160 & volcano <=180, 1, NA) # This is the "mask" raster
r2 <- raster(over)
# And this is the key step:
# To convert the "mask" raster to spatial points
r.mask <- rasterToPoints(r2, spatial=TRUE)
# Plot
levelplot(r1, margin=F) +
layer(sp.points(r.mask, pch=20, cex=0.3, alpha=0.8))
귀하의 경우 다음과 같이 코드를 변경할 수 있습니다.
1) "이진"래스터를 생성합니다. 예를 들어 0은 비가오고 1은 관개되며, 내가 한 것처럼 또는 플롯 한 것처럼 플롯합니다.
2) 두 개의 서로 다른 마스크를 만듭니다 (하나는 비를위한 것과 다른 하나는 관개를위한 것). 그리고 둘 다를 서로 다른 물체 (예 : 비가 오는 원과 비가 오는 곳의 십자가)로 플로팅합니다.
두 번째 대안은 맵이 너무 시각적으로 "오염 된"것처럼 보이게한다고 생각하지만 sp
패키지 설명서 를 참고하여 포인트의 매개 변수를 항상 제어 할 수 있습니다 ?sp.points
.
도움이 되길 바랍니다.