R 如何放大光栅打印(将限制设置为光栅)

R 如何放大光栅打印(将限制设置为光栅),r,plot,zooming,raster,R,Plot,Zooming,Raster,我已经绘制了山体阴影和滑坡敏感性的叠加图。我想放大情节的某个特定区域 landslide_raster_pred <-predict(ta,landslide_model, type = "response") hillshade <- hillShade(slope_data, aspect_ratio, angle=45, 0) overlay_map <- overlay(landslide_raster_pred,hillshade, fun=function(x,

我已经绘制了山体阴影和滑坡敏感性的叠加图。我想放大情节的某个特定区域

landslide_raster_pred <-predict(ta,landslide_model, type = "response")

hillshade <- hillShade(slope_data, aspect_ratio, angle=45, 0)

overlay_map <- overlay(landslide_raster_pred,hillshade, fun=function(x,y){return(x*y)})    


dev_coordinate=data.frame(x=714717.7, y=9560497)

# Extract the susceptibility at the point of interest from the raster
dev_prediction <- extract(landslide_raster_pred, dev_coordinate)

# Plot the overlay map with the point of interest 
plot(overlay_map, main = "Landslide susceptibility and hillshade map", 
 xlab = "longitude",
 ylab ="latitude")
points(c(714717.7), c(9560497))

slademy\u graster\u pred请使用相关的示例数据(不是我们无法访问的您自己的数据)

还是像这样

s <- 500
e <- extent(pt[1]-s, pt[1]+s, pt[2]-s, pt[2]+s)
plot(r, ext=e)
zoom(r, e)
也可以通过在绘图上绘制矩形以交互方式执行此操作

plot(r)
zoom(r, drawExtent())

(请注意这里的“经度”和“纬度”标签明显错误)
plot(r)
zoom(r, drawExtent())
plot(r,  xlim=c(pt[1]-s, pt[1]+s), ylim=c(pt[2]-s, pt[2]+s))