问题描述
我在合并两个具有不同分辨率和范围且具有不规则点数据的栅格数据集时遇到问题。以下是每个栅格的信息。我需要根据网格点合并这些数据集,以便我可以运行火势蔓延模型(需要风和坡度数据)。
我尝试在合并之前转换为普通数据帧(使用 rasterToPoints),但差异导致大量网格点丢失。
我曾尝试将栅格与项目栅格和栅格化对齐,但我没有设法让它工作。如果有人有可以提供帮助的想法,我将非常感谢您的回答!
'''
来自https://globalwindatlas.info/downloads/gis-files的数据
> wind <- raster("ZAF_wind-speed_10m.tif"); wind
class : RasterLayer
dimensions : 11271,11804,133042884 (nrow,ncol,ncell)
resolution : 0.0025,0.0025 (x,y)
extent : 13.33407,42.84407,-50.31423,-22.13673 (xmin,xmax,ymin,ymax)
crs : +proj=longlat +datum=wgs84 +no_defs
source : C:/Users/s2000128/Documents/Flammability modelling/ZAF/ZAF_wind-speed_10m.tif
names : ZAF_wind.speed_10m
来自https://datacatalog.worldbank.org/dataset/world-slope-model的数据
> slope <- raster("ZAF1_msk_alt.Grd"); slope
class : RasterLayer
dimensions : 1548,1992,3083616 (nrow,ncell)
resolution : 0.008333333,0.008333333 (x,y)
extent : 16.4,33,-34.9,-22 (xmin,ymax)
crs : +proj=longlat +ellps=wgs84 +no_defs
source : C:/Users/s2000128/Documents/Flammability modelling/slope_deg_0/slope_deg/ZAF1_msk_alt.Grd
names : ZAF1_msk_alt
values : -26,3264 (min,max)
以下是我尝试过的一些示例:
1
windresampled <- projectRaster(wind,slope,method = 'ngb'); windresampled
2
wind_points <- rasterToPoints(wind);
coordinates(wind_points) = ~x+y;
proj4string(wind_points) = CRS("+init=epsg:4326");
gridded(wind_points) = TRUE;
g_wind <- raster(wind_points); g_wind;
extent(wind) <- c(16.4,-22);
res(wind) <- c(0.0025,0.0025);
r_wind <- rasterize(wind_points,g_wind,field = wind_points$z,fun = mean,na.rm = TRUE); r_wind
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)