R 收集高程数据-从geotiff文件中提取地理位置
我正在尝试使用R 收集高程数据-从geotiff文件中提取地理位置,r,maps,raster,geotiff,R,Maps,Raster,Geotiff,我正在尝试使用rayshader包将高程数据添加到绘图中。我可以使用传单软件包绘制要查找高程数据的区域 library(whitebox) library(leaflet) library(rayshader) library(rayrender) library(raster) # define bounding box with longitude/latitude coordinates bbox <- list( p1 = list(long = -3.6525599, la
rayshader
包将高程数据添加到绘图中。我可以使用传单
软件包绘制要查找高程数据的区域
library(whitebox)
library(leaflet)
library(rayshader)
library(rayrender)
library(raster)
# define bounding box with longitude/latitude coordinates
bbox <- list(
p1 = list(long = -3.6525599, lat = 40.4065001),
p2 = list(long = -3.7525599, lat = 40.4965001)
)
leaflet() %>%
addTiles() %>%
addRectangles(
lng1 = bbox$p1$long, lat1 = bbox$p1$lat,
lng2 = bbox$p2$long, lat2 = bbox$p2$lat,
fillColor = "transparent"
) %>%
fitBounds(
lng1 = bbox$p1$long, lat1 = bbox$p1$lat,
lng2 = bbox$p2$long, lat2 = bbox$p2$lat,
)
我用我的坐标替换了链接中的lat长坐标:
https://elevation.nationalmap.gov/arcgis/rest/services/3DEPElevation/ImageServer/exportImage?bbox=-3.6525599%2C40.4065001%2C-3.7525599%2C40.4965001&bboxSR=4326&size=600%2C480&imageSR=4326&time=&format=jpgpng&pixelType=F32&noData=&noDataInterpretation=esriNoDataMatchAny&interpolation=+RSP_BilinearInterpolation&compression=&compressionQuality=&bandIds=&mosaicRule=&renderingRule=&f=html
它返回了一个空图像-这并不奇怪,因为USGS
是美国的数据提供商。因此,我下载了以下文件:
DTM西班牙大陆2000万
该文件大小为2.4 GB,是整个西班牙的一个文件,但我只需要lat&long点边界的高程数据
rgdal::GDALinfo("DTM Spain_Mainland (2019) 20m.tif")
rows 48394
columns 58187
bands 1
lower left origin.x -23380
lower left origin.y 3901190
res.x 20
res.y 20
ysign -1
oblique.x 0
oblique.y 0
driver GTiff
projection +proj=utm +zone=30 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m
+no_defs
file DTM Spain_Mainland (2019) 20m.tif
apparent band summary:
GDType hasNoDataValue NoDataValue blockSize1 blockSize2
1 Float32 TRUE -32767 1 58187
apparent band statistics:
Bmin Bmax Bmean Bsd
1 -4294967295 4294967295 NA NA
Metadata:
AREA_OR_POINT=Point
我的问题是,我如何在
bbox
中获取长和横向位置的geotiff,或者将当前DTM西班牙大陆(2019)20m.tif文件过滤到长、横向坐标?您可以查看elevatr
软件包,该软件包将为您提供所需的光栅:
library(elevatr)
library(raster)
bbox2 <- data.frame(x = c(-3.6525599, -3.7525599), y = c(40.4065001, 40.4965001))
elev <- get_elev_raster(bbox2, z = 13, clip = "bbox",
prj = "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")
库(elevatr)
图书馆(光栅)
bbox2
library(elevatr)
library(raster)
bbox2 <- data.frame(x = c(-3.6525599, -3.7525599), y = c(40.4065001, 40.4965001))
elev <- get_elev_raster(bbox2, z = 13, clip = "bbox",
prj = "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")