我正试图通过ArcGIS REST API从澳大利亚政府获得一个光栅层(5米数字高尺模型(DEM))
的url: https://services.ga.gov.au/site_9/rest/services/DEM_LiDAR_5m/MapServer/
我可以对边界和足迹层执行查询，结果很好
边界查询:
https://services.ga.gov.au/site_9/rest/services/DEM_LiDAR_5m/MapServer/1/query?where=&text=&objectIds=&time=&geometry=xmin:133.226100657682, ymin:-23.6165145646468, xmax:133.235899342318，ymax: -23.6074854353532, geometryType = esriGeometryEnvelope& inSR = 4283, spatialRel = esriSpatialRelEnvelopeIntersects& relationParam =和外场=,returnGeometry = true& maxAllowableOffset =, geometryPrecision =, outSR =, returnIdsOnly = false& returnCountOnly = false& orderByFields =, groupByFieldsForStatistics =, outStatistics =, returnZ = false& returnM = false& gdbVersion =, returnDistinctValues = false& f = json"
响应(以R表示):
通过jsonlite响应绘制R图
足迹查询:
https://services.ga.gov.au/site_9/rest/services/DEM_LiDAR_5m/MapServer/2/query?where=&text=&objectIds=&time=&geometry=xmin:133.226100657682, ymin:-23.6165145646468, xmax:133.235899342318，ymax: -23.6074854353532, geometryType = esriGeometryEnvelope& inSR = 4283, spatialRel = esriSpatialRelEnvelopeIntersects& relationParam =和外场=,returnGeometry = true& maxAllowableOffset =, geometryPrecision =, outSR =, returnIdsOnly = false& returnCountOnly = false& orderByFields =, groupByFieldsForStatistics =, outStatistics =, returnZ = false& returnM = false& gdbVersion =, returnDistinctValues = false& f = json"
反应:
jsonlite response in R
我得到了我正在寻找的TIF文件的名称，尽管如此，我对栅格层的所有查询:https://services.ga.gov.au/site_9/rest/services/DEM_LiDAR_5m/MapServer/3没有给我任何结果
我试过rcran包" arcpull& quot;同样，但是没有结果
帮忙吗?
library(arcpullr)
library(sf)
library(jsonlite)
loc_lng <- 133.231
loc_lat <- -23.612
study_area <- 1000
bbox_area <- osmar::center_bbox(loc_lng, loc_lat, study_area, study_area)
loc_bbox_pol <- tmaptools::bb_poly(bbox_area)
loc_bbox_pol <- sf::st_as_sf(loc_bbox_pol)
sf::st_crs(loc_bbox_pol) = 4326
loc_bbox_pol <- sf::st_transform(loc_bbox_pol, 4283)
DEM_general <- arcpullr::get_layer_by_poly("https://services.ga.gov.au/site_9/rest/services/DEM_LiDAR_5m/MapServer/0/", loc_bbox_pol, sp_rel = "esriSpatialRelIndexIntersects")
DEM_Boundary <- arcpullr::get_layer_by_poly("https://services.ga.gov.au/site_9/rest/services/DEM_LiDAR_5m/MapServer/1/", loc_bbox_pol, sp_rel = "esriSpatialRelIndexIntersects")
DEM_Info <- arcpullr::get_layer_by_poly("https://services.ga.gov.au/site_9/rest/services/DEM_LiDAR_5m/MapServer/2/", loc_bbox_pol, sp_rel = "esriSpatialRelIndexIntersects")
arcpullr::get_service_type("https://services.ga.gov.au/site_9/rest/services/DEM_LiDAR_5m/MapServer/3/")
所以很可能数据提供者已经在该地图服务的第三层禁用了查询操作。此外，在第3层上没有启用字段，因此没有任何可查询或通过查询操作返回的内容。
所以只有导出操作可以通过提供一个bbox(您也可以通过使用占用层的结果来使用它)来处理它。
https://services.ga.gov.au/site_9/rest/services/dem_lidar_5m/mapserver/export?bbox= - 5611089.372356728% - 2 c - 5642303.594604798% - 2 - c5611089.372356728%2c - 3900762.342155805, bboxSR = % 7 b % 22 wkt % 3 22% % 22 projcs % 5 b % 5 c % 22 wgs_1984_web_mercator_auxiliary_sphere % 5 c % 22% 2 cgeogcs % 5 b % 5 c % 22 gcs_wgs_1984 % 5 c % 22% 2 cdatum % 5 b % 5 c % 22 d_wgs_1984 % 5 c % 22% 2 cspheroid % 5 b % 5 c % 22 wgs_1984 % 5 c % 22% 2 c6378137.0 % 2 c298.257223563 % 5 d % d % 2 cprimem % 5 b % 5 c % 22格林威治% 5 c % 22% 2 c0.0 % 5 d % 2 cunit % 5 b % 5 c % 22度% 5 c % 22% 2 c0.0174532925199433 cprojecti % 5 d % 5 d % 2% 5 b % 5 c % 22 mercator_auxiliary_sphere % 5 c % 22% 5 d % 2 cparameter % 5 b % 5 c % 22 false_easting % 5 c % 22% 2 c0.0 % 5 d % 2 cparameter % 5 b % 5 c % 22 false_northing % 5 c % 22% 2 c0.0 % 5 d % 2 cparameter % 5 b % 5 c % 22 central_meridian % 5 c % 22% 2 c149.430392074675 % 5 d % 2 cparameter % 5 b % 5 c % 22 standard_parallel_1 % 5 c % 22% 2 c0.0 % 5 d % 2 cparameter % 5 b % 5 c % 22 auxiliary_sphere_type % 5 c % 22% 2 c0.0 % 5 d % 2 cunit % 5 b % 5 c % 22计% 5 c % 22% 2 c1.0 % 5 d % 5 d % 22% 7 d& imageSR = % 7 b % 22 wkt % 3 22% % 22 projcs % 5 b % 5 c % 22 wgs_1984_web_mercator_auxiliary_sphere % 5 c % 22% 2 cgeogcs % 5 b % 5 c % 22 gcs_wgs_1984 % 5 c % 22% 2 cdatum% 5 b % 5 c % 22 d_wgs_1984 % 5 c % 22% 2 cspheroid % 5 b % 5 c % 22 wgs_1984 % 5 c % 22% 2 c6378137.0 % 2 c298.257223563 % 5 d % d % 2 cprimem % 5 b % 5 c % 22格林威治% 5 c % 22% 2 c0.0 % 5 d % 2 cunit % 5 b % 5 c % 22度% 5 c % 22% 2 c0.0174532925199433 % 5 d % 5 d % 2 cprojection % 5 b % 5 c % 22 mercator_auxiliary_sphere % 5 c % 22% 5 d % 2 cparameter % 5 b % 5 c % 22 false_easting % 5 c % 22% 2 c0.0 % 5 d % 2 cparameter % 5 b % 5 c % 22 false_northing % 5 c % 22% 2 c0.0 % 5 d % 2 cparameter % 5 b % 5 c % 22 central_meridian % 5 c % 22% 2 c149.430392074675 % 5 d % 2 cparameter % 5 b % 5 c % 22 standard_parallel_1 % 5 c % 22% 2 c0.0 % 5 d % 2 cparameter % 5 b % 5 c % 22 auxiliary_sphere_Type % 5 c % 22% 2 c0.0 % 5 d % 2 cunit % 5 b % 5 c % 22计% 5 c % 22% 2 c1.0 % 5 d % 5 d % 22% 7 d&大小= 3441% 2 c534& dpi = 144,格式= png32&透明= true&层=显示% 3 a3& f = html