&&版权声明:本文为博主原创文章,未经博主允许不得转载。
一、绪论
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,自动绘制接边线,如下图所示,自动裁剪掉TM边缘“锯齿”。
(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内插算法引起的误差;