2016-12-20 6 views
1

나는 셰이프 파일의 점을 맨 위에 래스터 레이어에 플롯하려고합니다. 이것에 대한 이전 답변을 확인했지만 여전히 문제가 있습니다. 포인트 셰이프 파일과 래스터 레이어를 문제없이 별도로 플로팅 할 수 있지만 함께 사용할 수는 없습니다.래스터 레이어에 모양 파일을 그릴 때 R

나는 그들이 볼 수있는 한 동일한 투영과 위치에 있어야합니다.

require(maptools) 

myproj <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0" 
shape <- readShapeSpatial("directory/mypoints.shp", proj4string = CRS(myproj)) 

plot(r <- raster(listVI[200])) 
plot(shape) 

답변

1

나는 같은 문제가 발생할 수있는 다른 사람들을 위해 대답을 찾았습니다.

솔루션

은 간단하다 : (래스터 및 Shape 파일이 동일한 CRS에있는만큼)

plot(r) 

plot(shape, add = TRUE) 
+0

도 rasterVis''와'layer'를 사용할 수 있습니다. '? rasterVis :: levelplot'의 예제를보십시오. – jbaums