#######################################################
## 使用矢量文件裁剪栅格数据
raster_data # 待裁剪的栅格数据
shapefile_data # 读入的shapefile文件
# 检查栅格数据和shapefile数据的投影是否一致
projection(raster_data)
projection(shapefile_data)

# 如果投影不一致,则需要将其中一个数据集投影到另一个的投影坐标系
shapefile_data_proj <- spTransform(shapefile_data, crs(raster_data))

# 获取shapefile数据的Extent
ext <- extent(shapefile_data_proj)

# 使用shapefile数据裁剪栅格数据
cropped_raster <- crop(raster_data, ext)

# 最后,使用mask()函数根据shapefile数据裁剪栅格数据
clipped_raster <- mask(cropped_raster, shapefile_data_proj)
clipped_raster
# 查看结果
plot(clipped_raster[[1]])


posted @ 2024-05-12 13:51 junjunang 阅读(97) 评论(0) 推荐(0) 编辑
摘要: ###########################################raster_data = target.raster## 确定范围大小clp.ras.ext = extent(raster_data)## 确定分辨率大小 1 minmy.res = 1 * (1/60)## 阅读全文
posted @ 2024-05-12 13:49 junjunang 阅读(28) 评论(0) 推荐(0) 编辑
摘要: 在使用R中的plot()函数绘制地图时,需要设置的参数。在R中绘制地图时,为了确保经纬度像素大小相同,需要设置asp参数。asp参数代表"宽高比"(Aspect Ratio),它控制着地图上经纬度的像素尺寸。当asp=1时,经纬度的像素大小是相等的,这可以确保地图上的距离和面积比例是正确的。 阅读全文
posted @ 2024-05-08 11:19 junjunang 阅读(16) 评论(0) 推荐(0) 编辑
摘要: bio2,3,4,6,7,9,12,13,14, 15,16,17,18,19 在格陵兰岛存在显著问题如下: 有明显的分割线。 阅读全文
posted @ 2024-05-07 06:45 junjunang 阅读(31) 评论(0) 推荐(0) 编辑
摘要: 如果直接使用extent(a)= extent(b) 则可能出现下面的状况 为了保证范围一致 可以使用 new.a = setExtent(a, ext = extent(b), keepres = T) 注意里面的keepres是关键 阅读全文
posted @ 2024-04-22 13:58 junjunang 阅读(102) 评论(0) 推荐(0) 编辑
摘要: 找到这个文件的输入信息栏,下面找到file attachments框,右键就可以clear了 阅读全文
posted @ 2024-04-02 15:50 junjunang 阅读(10) 评论(0) 推荐(0) 编辑
摘要: 在CMIP6(Coupled Model Intercomparison Project Phase 6)中,Variant Label是用于标识不同模型变体的标签。它由四个部分组成,分别是r1i1p1f1。 r1 表示实验重复次数(run number),即同一模型下运行的不同实验的编号。r1表示 阅读全文
posted @ 2023-05-31 14:17 junjunang 阅读(2211) 评论(0) 推荐(0) 编辑
摘要: image函数绘图,按照值的大小,默认赋予一定颜色。 如果希望按照自己的想法,设置离散颜色的话,需要注意: 首先,去掉NA值,然后对你的图像值取唯一值,然后排序,从小到大。 例如: testm = matrix(c(NA,5,5,2,2,2,3,3,3,4,4,4),nrow = 3,byrow = 阅读全文
posted @ 2022-08-13 19:45 junjunang 阅读(170) 评论(0) 推荐(0) 编辑
摘要: labels = bquote('('*.(letters[i*2])*')') *号无空格,~表示有一个空格 阅读全文
posted @ 2022-07-29 16:03 junjunang 阅读(27) 评论(0) 推荐(0) 编辑
摘要: ## clip if(projection(thistif) != projection(mo.range)){ thistif = projectRaster(from = thistif,to = mo.range) } r2 <- crop(thistif, extent(mo.range)) 阅读全文
posted @ 2022-07-16 10:15 junjunang 阅读(179) 评论(0) 推荐(0) 编辑
点击右上角即可分享
微信分享提示