我正在尝试使用 mapview 包来绘制栅格和 sf 特征。
streamsTransform <- st_transform(streams, 2958)
st_crs(streamsTransform)
dem17N <- projectRaster(dem, crs = "+proj=utm +zone=17 +ellps=GRS80
+towgs84=0,0,0,0,0,0,0 +units=m +no_defs")
course <- raster::aggregate(dem17N,8)
mapview(course) + mapview(streamsTransform)
Run Code Online (Sandbox Code Playgroud)
我已经使用 st_transform 和 projectRaster 确保它们在相同的投影坐标系中,但我仍然收到以下错误:
(!is.na(getProjection(lst[[i]]))) { :
argument is of length zero
Run Code Online (Sandbox Code Playgroud)
任何建议将不胜感激!