3.输入“&r srtmgrid.aml”,表示运行我们刚才复制粘贴生成的 AML 程序。可以看到这个命
令的用法。这里提示我们在运行该程序时要输入两个参数:输入文件名和输出文件;
4.明白了第 3 步以后,让我们正式开始执行该程序。输入“&r srtmgrid.aml N39E114.hgt
N39E114”。注意该空格的地方要空格。该命令表示运行(run)srtmgrid.aml 程序,输入文
件为 N39E114.hgt,输出文件为 N39E114 文件。由文件名我们也能看到:该文件的经纬度方
格的范围是北纬 39 东经 114.点击回车,即可执行。
5.再次输入 quit,就从 GRID 环境退回到 arc 环境中,再次重复输入第 4、5 步的命令,只需
将输入输出文件名改一下,就能将其他 hgt 文件转成 GRID 文件了,直到把所有的 hgt 文件
全部转换完成为止。
6.输入“grid”,进入 GRID 环境。我们要将第 4、5 步转换的 GRID 文件拼接(mosaic)成一
个文件。输入“mosaic”,回车,就能看到 mosaic 命令的用法。我们要指定一个输出文件
名 , 再将 第 4 、 5 步 生 成 的 GRID 文件 名 依 次 输 进 去 , 就 行 了 。 我 们 输 入 “ bjdem =
mosaic(N39E114,N39E115,N39E116,N39E117,N40E114,N40E115,N40E116,N40E117,N41E114,N
41E115,N41E116,N41E117,N42E114,N42E115,N42E116,N42E117)”,回车,就把这些 GRID 文
件全都拼接成一张图了。
7.有了第 6 步的成果,数据就算转换完成了,就可以用来进行自己的工作了。如果还想转成
tiff 格式或 shp 格式等,在 ArcMap 里面操作即可,也不难的。这里就不再赘述了。
srtmgrid.aml 代码如下:
/*--------------------------------------------------------------------------
/* Program: SRTMGRID.AML
/* Purpose: Convert raw prelminary SRTM elevation data files
/* (.hgt extension) to ESRI GRID format. Many thanks to
/* Arthur Tarr, USGS, (atarr@usgs.gov) for figuring out how
/* to convert SRTM image files to ArcInfo GRID format.
/*
/* The georeferencing and grid-cell resolution
/* for the output grid is determined from the file
/* name and the size of the file.
/*--------------------------------------------------------------------------
/* Usage: &r srtmgrid <dat_file> <out_grid>
/*
/* Arguments:
评论0
最新资源