为insar处理创建上采样dem
sardem的Python项目详细描述
DEM创建者
用于干涉合成孔径雷达(in sar)处理的roipac数据格式(16位整数,小尾数)数字高程图(dem)制作工具
createdem
创建裁剪(也可能是上采样)的数字高程图:
usage: createdem left_lon top_lat dlon dlat [-h][--rate RATE=1][--output OUTPUT=elevation.dem][--data-source {NASA,AWS}]
设置和安装
pip install sardem
这将创建命令行可执行文件createdem
或者,您可以克隆以生成/安装:
git clone https://github.com/scottstanie/sardem
cd sardem
make
它将运行pip install --upgrade .
,并创建命令行脚本。
如果您使用virtualenv
# Optional for using virtualenv virtualenv ~/envs/sardem &&source ~/envs/sardem/bin/activate # Or wherever you store your virtual envs# Or if you have virtualenv wrapper: mkvirtualenv sardem pip install sardem
命令行界面
可以使用sardem/cli.py
找到命令行工具的完整选项
$ createdem --help
usage: createdem left_lon top_lat dlon dlat
[-h] [--rate RATE=1] [--output OUTPUT=elevation.dem]
[--data-source {NASA,AWS}]
Stiches SRTM .hgt files to make (upsampled) DEM
Pick a lat/lon bounding box for a DEM, and it will download
the necessary SRTM1 tiles, stitch together, then upsample.
Usage Examples:
createdem -156.0 20.2 1 2 --rate 2 # Makes a box 1 degree wide, 2 deg high
createdem -156.0 20.2 0.5 0.5 -r 10 --data-source NASA -o my_elevation.dem
Default out is elevation.dem for the final upsampled DEM.
Also creates elevation.dem.rsc with start lat/lon, stride, and other info.
positional arguments:
left_lon Left (western) most longitude of DEM box (degrees, west=negative)
top_lat Top (northern) most latitude of DEM box (degrees)
dlon Width of DEM box (degrees)
dlat Height of DEM box (degrees)
optional arguments:
-h, --help show this help message and exit
--rate RATE, -r RATE Rate at which to upsample DEM (default=1, no upsampling)
--output OUTPUT, -o OUTPUT
Name of output dem file (default=elevation.dem)
--data-source {NASA,AWS}, -d {NASA,AWS}
Source of SRTM data (default NASA). See README for more.
在upsampling例程中用于双线性插值的代码位于cython/upsample.c
中,并包装在cython中,以便于安装和从python调用函数。
安装是通过pip install
,或通过运行make build
来处理的。
使用数字高程图(DEM)的函数主要包含在Downloader
和Stitcher
类中。
nasa srtm数据访问
默认数据源是nasa的srtm version 3全局1度数据。 见https://lpdaac.usgs.gov/dataset_discovery/measures/measures_products_table/srtmgl3s_v003。 该数据在北极地区以外(-60到60度纬度)有效,在公海上空为零。
此数据需要用户名和密码: https://urs.earthdata.nasa.gov/users/new
当使用nasa数据运行时,将提示您输入用户名和密码。 它将保存到~/.netrc文件中以备将来使用,这意味着您以后不必输入用户名和密码。 条目将如下所示:
machine urs.earthdata.nasa.gov
login USERNAME
password PASSWORD
如果要完全避免这种情况,可以通过指定
sardem 156.0 20.0 .5 0.5 --data-source AWS
--data-source NASA
是默认值。
mapzen将srtm数据与其他源相结合,因此.hgt文件将略有不同(但通常不明显)
警告:mapzen注意到他们正在停止一些服务,这就是为什么nasa是默认的。