显示要素周围的附加区域，get_elev_raster()（来自 elevatr）答案

【问题标题】：show additional area around feature, get_elev_raster() (from elevatr)显示要素周围的附加区域，get_elev_raster()（来自 elevatr）
【发布时间】：2021-06-15 11:26:03
【问题描述】：

在回答此问题"Create topographic map in R" 之后，我想扩展边界框以显示特征周围的其他区域。 documentation 提到了“expand”参数，但我不确定如何应用它，而且我的谷歌搜索没有返回我能够成功使用的代码示例。

我正在测试一个数据框，目前只有一组坐标 (Adam's peak, from Wikipedia)，如下所示：

dataframe image

代码同here，如下，除了第一行“#生成一个经纬度坐标的数据框”

library(elevatr)
library(raster)

#input dataframe
ex.df <- mtSlice

# Specify projection.
prj_dd <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"

# Use elevatr package to get elevation data for each point.
elev <- get_elev_raster(ex.df, prj = prj_dd, z = 10, clip = "bbox")

raster::contour(elev)

谁能解释如何展开以显示该功能周围的更多区域？

【问题讨论】：

标签： r maps elevatr

【解决方案1】：

get_elev_raster 上的 expand 参数将在您的输入位置以 CRS 为单位展开边界框。在这种情况下，这将是十进制度。因此，要在边界框上再增加 1 度：

elev_plus_one_deg <- get_elev_raster(ex.df, prj = prj_dd, z = 10, clip = "bbox", expand = 1)

注意：它应该会多出 1 度，但仔细观察似乎并没有。我需要深入研究！

【讨论】：

谢谢@jhollist，你的意思是它似乎在指定的度数上增加了一个不同的度数，还是根本没有增加一个度数？
它确实增加了该区域，只是不是精确的度数（如果以十进制度为单位）。 CRAN 上的版本使用缓冲区，非投影/地理数据的缓冲区并不理想。在这种情况下，它仍然会为您提供更大的区域，只是不会大一度（或您选择的任何值）。您需要进行实验。
@WithRegards 只是好奇您的问题是否已解决？还有其他问题吗？