多边形和点层无法投影计算获取

问题描述

我想用 fetchR 计算海岸的风暴露。 我准备了两层；站点（点）和海岸线（多边形），并使用 sp 函数设置相同的 CRS。但是代码显示错误“它们必须被投影......”。这两层不是被投影了吗？

library(sp)
library(fetchR)
long <- c(133.2982,132.6715,133.2375,133.3048,133.2594,133.2165)
lat <- c(35.5716,35.3551,35.5504,35.5707,35.5680,35.5708)
lonlat <-data.frame(cbind(long,lat))

fetch_locs_proj = sp::SpatialPoints(lonlat,# points at coasts
                           proj4string = sp::CRS("+init=epsg:4326"))

nz_poly_proj = rgdal::readOGR("file.shp")  # coast line shapefile 
nz_poly_proj1 = spTransform(nz_poly_proj,sp::CRS("+init=epsg:4326"))

my_fetch_proj = fetch(polygon_layer = nz_poly_proj1,site_layer = fetch_locs_proj,max_dist = 3,n_directions = 9,site_names = fetch.df$name)

*Error: polygon_layer and/or site_layer must be projected to calculate fetch*

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）