data.frame和get_elev_points.位置错误$海拔[idx_s[i]：idx_e[i]] <- unlist(resp$height)：替换长度为零

Question

我试图在数据帧中使用函数get_elevation_points，但始终得到相同的错误：

谢谢!

> #get the API key : et an API Key from https://mapzen.com/developers. Since elevatr will search the environment for the appropriate key, 
> #it is suggested that you store this in your .Renviron file. To do this:
> 
> require(sf) 
  require(sp)
> 
> cat("mapzen_key=mapzen-XXXXXXX\n", file =
> file.path(normalizePath("~/"), ".Renviron"), 
>     append = TRUE) mydata <- st_as_sf(x = Coordinates_alt, 
>                         coords = c("Longitude", "Latitude"),
>                         crs = "+proj=longlat +datum=WGS84")

> 

> coordinates(Coordinates_alt) <- cbind(Coordinates_alt$Longitude ,
> Coordinates_alt$Latitude) proj4string(Coordinates_alt) =
> CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0") 
> 
> require(elevatr)
> 
> prj_dd <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"
> 
> mydata1 <- SpatialPointsDataFrame(df, proj4string = CRS(prj_dd), data
> = df)
> 
> head(mydata1)
> 
> df_elev <- get_elev_point(mydata1, prj = prj_dd, src = "mapzen")

df_elev <- get_elev_point(df，prj = prj_dd，src = "mapzen")在位置$idx_si:idx_ei[idx_si:idx_ei] <- unlist(resp$height)中的错误： 替换的长度为零，此外:警告信息: src 'mapzen‘被废弃，并将停止运行2018-01-31后，因为关闭Mapzen；使用'epqs’而不是美国的地点。仍在寻找全球标高服务的替代品。

Answer 1

很抱歉没有及时回复。我对此并不十分积极，但考虑到elevatr的几个问题，我会这么做的！

您遇到的问题是，src = "mapzen"是一个旧的来源，因为所有的Mapzen服务都已经关闭。如果您的坐标在美国，您可以使用src = "epqs"，如果它们在美国以外，您可以使用src = "aws"，但取决于您的边界框的多少点和范围，这可能需要一些时间。另外，一定要读取elevatr::get_elev_point()的文档，因为其中还有一些关于为aws源代码设置z值的附加信息。像elevatr::get_elev_point(mydata1, prj = prj_dd, src = "aws", z = 5)这样的东西会很快(Ish)，但是源光栅数据的分辨率相当粗糙。您可以更改z(最多14)，但是更高的值将为下载增加时间，这同样取决于您的坐标范围。