微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

组合2个具有不同范围、分辨率和点规律性的栅格数据集 来自https://globalwindatlas.info/downloads/gis-files的数据来自https://datacatalog.worldbank.org/dataset/world-slope-model的数据12

如何解决组合2个具有不同范围、分辨率和点规律性的栅格数据集 来自https://globalwindatlas.info/downloads/gis-files的数据来自https://datacatalog.worldbank.org/dataset/world-slope-model的数据12

我在合并两个具有不同分辨率和范围且具有不规则点数据的栅格数据集时遇到问题。以下是每个栅格的信息。我需要根据网格点合并这些数据集,以便我可以运行火势蔓延模型(需要风和坡度数据)。

我尝试在合并之前转换为普通数据帧(使用 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 举报,一经查实,本站将立刻删除。

相关推荐


Selenium Web驱动程序和Java。元素在(x,y)点处不可单击。其他元素将获得点击?
Python-如何使用点“。” 访问字典成员?
Java 字符串是不可变的。到底是什么意思?
Java中的“ final”关键字如何工作?(我仍然可以修改对象。)
“loop:”在Java代码中。这是什么,为什么要编译?
java.lang.ClassNotFoundException:sun.jdbc.odbc.JdbcOdbcDriver发生异常。为什么?
这是用Java进行XML解析的最佳库。
Java的PriorityQueue的内置迭代器不会以任何特定顺序遍历数据结构。为什么?
如何在Java中聆听按键时移动图像。
Java“Program to an interface”。这是什么意思?