&&版权声明：本文为博主原创文章，未经博主允许不得转载。
一、绪论
1.1DEM误差分类
1.2DEM误差传播过程
二、软件使用
1、ArcMap 10.2
2、ENVI 5.1（此版本及更高版本添加了无缝拼接工具）或更高版本
三、主要流程及结果
3.1数据准备
完成 图像拼接与投影坐标系统一
3.1.1打开四幅影像
在ENVI中打开四幅STRM影像文件-->打开文件选项
3.1.2在ENVI的Toolbox中选择工具
3.1.2选择需要拼接的四幅影像
3.1.3打开拼接预览
勾选Show Preview，可以预览镶嵌效果。
3.1.4匀色处理
提供的匀色方法是直方图匹配。
（1） 在Corlor Correction选项中，勾选Histogram Matching：
Overlap Area Only：重叠区直方图匹配
Entire Scene：整景影像直方图匹配
（2）在Main选项中，Color Matching Action中右键设置参考（Reference）和校正（Adjust），根据预览效果确定参考图像。
3.1.5接边线与羽化
接边线包括自动和手动绘制两种方法，也可以结合起来使用。
（1） 选择下拉菜单Seamlines->Auto Generate Seamlines，自动绘制接边线，如下图所示，自动裁剪掉ＴＭ边缘“锯齿”。
（2） 自动生成的接边线比较规整，可以明显看到由于颜色不同而显露的接边线。下拉菜单Seamlines-> Start editing seamlines，编辑接边线，可以在接边处绘制多边形，之后自动将绘制的多边形作为新的接边线。
3.1.6输出结果
切换到Export选项，这是输出文件名、路径、格式、波段、背景值、重采样方法等信息。
3.1.7统一投影坐标系（ArcMap）
Ø 这一步我们选择在ArcMap中操作，在ArcCatalog右击影像在属性中查看投影坐标信息。
① SRTM拼接前的四幅原始影像投影坐标系（如 下图）基准：GCS-WGS-84（geographic coordinate systems，简称GCS），拼接后依旧是这个投影坐标系；
而我们所给的Landsat影像投影坐标系（如 下图）；
投影坐标系不相同，因此需要统一投影坐标系。
② 我们将四幅拼接而成的以整幅影像投影坐标系与Lansat影像统一
其中，点击按钮 选择要转换的投影坐标系时可以采用引入已有文件坐标系方式：
点击按钮下的引入坐标系->搜索到Landsat影像，将其坐标系引入
③ 现在查看转换投影坐标系结果：已经和Lansat相同
3.2在ENVI中生成DEM
3.2.1加载ASTER影像（3N和3B）
① 打开影像； 打开方式：Fileà打开外来文件—>thermalàASTER-->搜索文件路径
② 打开ASTER中的3N影像，新建窗口打开3B，这样可以对比显示。
3.2.2打开dem提取向导
3.2.3加载左右两幅图像
加载左影像：3N 右影像：3B，然后点击下面的Save按钮进行保存第一步骤；
3.2.4提取GCPs（Ground Control Points）类型
① 交互式手动添加。
② 裁剪好的SRTM影像与Landsat影像（投影坐标系已经统一）叠加寻找控制点，同时在ENVI中寻找寻找同一个同名点（控制点）
需注意：我们叠加SRTM与Landsat影像目的更好寻找同名点，因为直接在SRTM（精度较差，分辨率较低）上面找同名点不好找，而LANSAT影像（分表率较高）上面找同名点比较好找，所以地接起来比较容易确定同名点；
③ 进入下一步后，首先选择坐标系统（如 下图）；
④ 在ASTER 3N和3B影像上找到与裁剪好的SRTM影像与Landsat影像（投影坐标系已经统一）叠加所确定的并对应的第一个同名点，输入高程，最后点击”Add”来添加第一个同名点（控制点）；
⑤ 同样继续天机几个控制点，注意分布要均匀。
⑥ 最终添加了六个控制点，点击Show talble,可以修改误差较大的点；
⑦ 同样保存第3步的页面
3.2.5自动添加GTPs（连接点）
① 选择25个链接点，自动生成方式；
② 依旧保存页面。然后下一步
③ 进入下一步即第5步：将删除偏差较大的连接点，并添加接插位置的连接点，使得值尽可能小，一般实习20左右就可以提取DEM，项目要求可能达到0.5~0.8
3.2.6生成核线影像
同样保存上面页面。进入下一步，保存核线影像：
3.2.7设置DEM输出参数
3.2.8在ArcMap中加载生成的DEM
3.3根据生成DEM范围裁剪SRTM
① 首先查看两幅影像，大小范围如下；
② 新建面文件，用于确定裁剪范围，来掩磨SRTM影像
③ 打开编辑器，对mask文件新建面要素；
④ 选择AacMap下ToolBox中掩磨栅格数据的工具；
⑤ 参数选择：
⑥ 掩磨结果：
3.4采用格网点定量分析DEM
① 使用重采样工具（可以使用三种重采样方法：双线性内插、最邻近法内插、三次卷积内插，都做一遍比较那个效果更好），生成栅格大小为600x600m的栅格数据
② 转换为点图层，使用如下转换工具。
输入栅格数据集中的每个像元都会在输出中转换为点。每个新生成的点都被定位在它所表示的像元的中心。NoData像元不会转换为点。
③ 转换为点图层结果
④ 在ArcToolbox中，找到Spatial Analyst Tools->Extraction->Extract Values to Points工具，然后将DEM以及要对比的srtm_mask作为栅格数据，将以上采集的shapefile文件作为矢量文件checkpoints.shp，输出获得高程值的shapefile文件。
⑤ 如图，输出结果已经添加高程信息
⑥ 同样对我们生成的DEM做同样操作
⑦ 导出两个shp文件为.txt格式，以便在EXCEL中添加高程信息，分析生成DEM的精度
⑧ 将checkpoints_masksrtm.shp导出结果命名为参考高程.txt，checkpoints_dem.shp命名为生成DEM高程.txt
⑨ 在EXCEL中导入两个文本数据。点击菜单栏数据-->导入数据，注意要逗号隔开；添加进去注意修改两个字段
⑩ 最后呈现折线图（三种内插方式），可以看到走势起伏情况基本相同，但是不能完全重合表示存在误差。
3.5在SPSS软件中分析结果
① 进行双变量相关分析
② 分析结果如下：
③ 结果显示在0.01水平（双侧）上显著相关
四、问题总结
① 生成的DEM中存在负值；
② 在ENVI中生成DEM有一步，添加控制点时，出现需要添加坐标信息的错误，如下图：
五、问题解决
① 我们在生成DEM时需要添加控制点，即统一到真实DEM的高程的上面（否则出现下面图 1折线图情况，总体差同一个高程值），
这一步跟GPS中GPS卫星坐标计算相似
u 1.计算真近点角
u 2.计算摄动改正项
u 3.计算升交距角、卫星的地心距离及轨道倾角
u 4.计算卫星在轨道坐标系中的坐标（x, y, z）(这里的X轴指向了升交点）
u 5. 计算升交点的经度
u 6.计算在协议天球坐标系中的空间直角坐标
u 7.计算在协议地球坐标系中的空间直角坐标
u 8.考虑极移的影响，最后得到在协议地球坐标系中的空间直角坐标
图 1
就是需要将我们要生成的DEM拉到真实高度，否则生成的DEM只是一个相对的高程信息集合。
同时也跟摄影测量学中物空间坐标系确定相似，由像平面坐标系-->像空间坐标系-->像空间辅助坐标系-->摄影测量坐标系-->物空间坐标系 一步步转换相似。
② 第二个问题，需要在ENVI中添加控制点高程值时需要同时添加同名点的位置信息。如下图：
六、实习总结
通过这次学习，学到好多东西。平时操作中没有注意的地方，进行了学习与补充。
首先，要统一坐标系，这是非常重要的，因为不同坐标系之下没有可比性，另外不同坐标系下数据肯定存在位移，统一坐标系后，同一个位置才会重合；
其次，以前我都是用的ENVI 5.0，里面没有无缝拼接工具（但是有完成拼接的工具），这个学期，装得ENVI 5.1升级版，里面有无缝拼接工具，我在图像拼接方面学习并使用了无缝拼接工具，效果还不错；
然后，DEM生成过程，出现了问题，就是生成的DEM中存在负值（数据本身不应该有负值，因为采用高原影像），而且每个高程点总是和参考SRTM影像中高程点高程相差几乎同一个高程值。经过查阅资料以及和同学交流经验，发现需要在生成DEM过程中不仅要添加连接点（同名点的链接），还要添加GCPs即控制点，这样生成的DEM才不会有负值，而且跟SRTM在同一水平面上；
接着，在EXCEL生成折线图，查看了生成DEM与参考影像SRTM的高程偏差，同时，本学期学习了计量地理学课程，因此使用SPSS经行相关分析，得到精度分析后结果。
最后我分析了误差生成原因：
① 首先，数据本身的原因。得到更高精度的结果需要分辨率再高些的影像数据；
② 其次，在找GTPs（连接点）前需要添加较多GCPs（控制点），而且要分布均匀，不能只采集一个区域，忽略其他地方。因为完成作业并非项目要求，所以我采集了六个控制点，如果可能的话可以在添加更多的控制点；
③ 然后，对实际地形表面采用引起的误差，即原始数据的误差；
④ 接着，重采样引起的误差，即有DEM内插算法引起的误差；
3842
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
