Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/76.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
R 计算多边形的坡度和坡向_R_Geospatial_Raster_R Raster_Terrain - Fatal编程技术网

R 计算多边形的坡度和坡向

R 计算多边形的坡度和坡向,r,geospatial,raster,r-raster,terrain,R,Geospatial,Raster,R Raster,Terrain,我有一个由125个唯一多边形组成的多边形sf对象(代表分散在美国各州的不同区域)。所有多边形都不共享边界。我还有一个30米的山形数字高程模型光栅图像,它覆盖了多边形区域(还有一些)。我想使用光栅软件包中的地形函数来计算每个多边形的唯一坡度和坡向。我的最终产品是让polygon sf对象具有两个新的列,slope和aspect,因此每个多边形都有它的slope和aspect 我在原始全州光栅图像和仅多边形区域的遮罩光栅图像上使用了地形函数,该函数返回整个光栅图像的最小和最大坡度和坡向 除此之外,我

我有一个由125个唯一多边形组成的多边形sf对象(代表分散在美国各州的不同区域)。所有多边形都不共享边界。我还有一个30米的山形数字高程模型光栅图像,它覆盖了多边形区域(还有一些)。我想使用光栅软件包中的地形函数来计算每个多边形的唯一坡度和坡向。我的最终产品是让polygon sf对象具有两个新的列,slope和aspect,因此每个多边形都有它的slope和aspect

我在原始全州光栅图像和仅多边形区域的遮罩光栅图像上使用了地形函数,该函数返回整个光栅图像的最小和最大坡度和坡向

除此之外,我对如何计算每个唯一多边形的面积的坡度和坡向,然后将数据输入到多边形sf对象中感到非常困惑。如果有人有什么建议,我很乐意听

谢谢大家!

示例数据

library(raster)
elevation <- getData('alt', country='CHE')
x <- terrain(elevation, opt=c('slope', 'aspect'), unit='degrees')
库(光栅)
提出一个人们可以帮助解决的R问题。这包括一个数据样本、所有必要的代码,以及对您正在尝试做什么和没有做什么的清晰解释。
sw <- getData("GADM", country="CHE", level=1)
e <- extract(x, sw, fun=mean, na.rm=TRUE)
x <- data.frame(sw$NAME_1, e)
head(x)
#               sw.NAME_1    slope   aspect
#1                 Aargau 3.065123 170.2449
#2 Appenzell Ausserrhoden 4.905536 187.5591
#3  Appenzell Innerrhoden 7.380191 170.1497
#4       Basel-Landschaft 3.742552 178.2771
#5            Basel-Stadt 1.002250 125.0466
#6                   Bern 9.025566 195.5117
swe <- cbind(sw, e)