如何从 ARCGIS 地图服务查询图块集以提取栅格图像？

问题描述

我正在尝试构建一个 url 查询来访问 ARCGIS 图块集（map here 和 services directory info here）并将其转换为 R 中的光栅图像。我看到“允许导出图块”条件是在服务目录上设置为 false，因此我尝试通过 "Tile Map" 选项进行查询。在文档中，我看到正确的格式是 https://<mapservice-url>/tilemap/<level>/<row>/<column>/<width>/<height>，但我不确定如何构建此 url 以便稍后提取栅格。

总的来说，我想查询裁剪到人口普查块组边界框的图块集并提取光栅图像。以下是人口普查块组多边形和获取边界框点的示例：

library(sf)
library(tigris)
library(tidyverse)
library(httr)
library(raster)
library(exactextractr)

ex_cbg <-
  block_groups("06","081",cb = T) %>%
  st_transform(2227) %>% 
  .[1,]

Box_list <-
    st_bBox(ex_cbg$geometry[1]) %>% unlist()

  element1 <- Box_list[[1]]
  element2 <- Box_list[[2]]
  element3 <- Box_list[[3]]
  element4 <- Box_list[[4]]

在其他非 tileset 查询中，我从以下代码中获取了栅格，其中包括边界框约束：

query_result <- GET(paste0("https://gis.smcgov.org/image/rest/services/SanMateoCounty_DEM2017/ImageServer/exportimage?bBox=",element1,",element2,element3,element4,"&size=2500,2500&format=tiff&f=json"))

result <- raster(content(query_result)$href)

解决方法

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

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

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