r - 绘制特立尼达和多巴哥的平滑降雨图
问题描述
我使用来自世界气候的数据为特立尼达创建了一张降雨图。但是,地图并不平滑。我想知道如何获得平滑的地图。我附上了我的代码和地图。任何帮助,将不胜感激。
library(dplyr)
library(tidyr)
library(sp)
library(ggplot2)
library(rgeos)
library(maptools)
library(rgdal)
library(raster)
library(rasterVis) #visualization library for raster
TT=getData('GADM', country='Trinidad and Tobago', level=1)
plot(TT) # this can be a little slow
clim=getData('worldclim', var='bio', res=2.5)
clim
gain(clim)=0.1
plot(clim[[12]])
gplot(clim[[12]])+geom_raster(aes(fill=value))+
facet_wrap(~variable)+
scale_fill_gradientn(colours=c("brown", "red", "yellow", "darkgreen", "green"),trans="log10") +
coord_equal()
## crop to a latitude/longitude box
r1 <- crop(clim[[12]], extent(-62,-60,10,12))
## Crop using a Spatial polygon
r1 <- crop(clim[[12]], bbox(TT))
plot(r1)
解决方案
您的数据:
library(raster)
TT <- getData('GADM', country='Trinidad and Tobago', level=1)
clim <- getData('worldclim', var='bio', res=2.5)
r <- crop(clim[[12]], TT)
解决方案:
x <- disaggregate(r, 5, method="bilinear")
y <- mask(x, TT)
plot(y)
推荐阅读
- mysql - 从一个表中获取所有 ID 并根据它在另一个表中插入值
- python - 如何从 obj 文件在 pyOpenGL 中导入 3d 纹理
- swift - 简而言之,使用 SwiftyJSON 将对象映射到数组中
- oracle - 从 hash_value 中提取 sql_id
- loops - 如果然后循环在 AutoHotkey 之间切换桌面
- angular - 检查后表达式已更改。以前的值:'null:1561881140820'。当前值:'null:1561881140880'
- python - 如果字符出现在字典中,则“加密”字符串
- r - 使用“neuralnet”包为不同层设置不同的激活函数
- python - 在文本文件中查找字符串:添加 else 语句后继续返回 False
- react-native - store.getstate() 评估为未定义