WorldWind系列十五:如何切割影像和DEM数据及其在WW中的应用配置
卫星数据的数据裁剪的方法
卫星数据的数据裁剪的方法我折腾了好久卫星数据的数据裁剪啊,总算找到点门道。
我一开始真的是瞎摸索,像个没头苍蝇似的。
最开始我想着这能有多难,不就是把不需要的数据去掉嘛。
我就直接按照数据的边界值去划分,觉得这样就能裁剪出来想要的部分。
结果呢,全错了,获得的数据要么少了重要部分,要么还夹杂了好多不相关的内容。
然后我就学聪明了些。
我了解到有一种按地理区域裁剪的方法。
比如说你要裁剪的卫星数据是关于某个特定城市的,就像画地图一样,你得先确定这个城市的经纬度范围。
我记得我当时处理关于纽约的数据。
我找到了纽约大致的经度从多少到多少,纬度从多少到多少。
这时候就像是用一把剪刀,沿着这个框架去裁剪出关于纽约的数据部分。
但是这里面也有坑啊,经纬度一不小心就容易弄错一个小数点之类的,那结果就偏差很大,我就因为这个反复调整了好多次。
还有一种方法是按照卫星影像的形状或者特定目标来裁剪。
比如说卫星数据中包含一个很小的湖泊的信息,我只想留下这个湖泊的数据,我就要先精确找出湖泊的轮廓。
这就像在一幅很复杂的画上想要单独裁出一朵花一样困难。
我得先识别出这个湖泊各个边缘的坐标点,这个过程中有时候会把旁边的陆地也误当成湖泊的一部分算进去,后来发现错误了,我就又得重新核对坐标点,一点一点调整,就像绣花一样精细。
我不确定是不是还有其他更好的方法了。
不过我觉得干这个事儿,一定要有耐心。
在裁剪之前,一定要多确认几遍自己的参数,不管是地理范围的经纬度,还是图像形状的坐标点,一旦弄错,那前面的功夫可就白做了。
而且,要是有条件的话,最好先在一小部分数据上做个裁剪实验,看看结果对不对,再大规模去裁剪整批数据。
我也试过使用一些工具,有些专业的图像处理软件可以专门处理卫星数据的裁剪。
这些软件有的操作比较复杂,但是功能也很强大。
像是有些软件,你在输入解析数据,比如说裁剪的坐标之类的时候,它有特定的格式要求,你要是搞错了这个,那也是不行的。
我当初就没仔细看要求,结果一直报错,后来才发现是输入格式不对。
高分辨率遥感影像分割方法综述
定工作 . 先后 主持 、 参 与完成 了多项 国家和 自治 区重大测绘地理信 息项 目的检 查验 收工作 , 并 主持 完成 了内蒙古 自治 区G P S  ̄ -
遥感测绘 2 0 1 6 年 第 六期
小, 或者灰度范围重叠较大时, 这种方法的结果往往不够理 言 , 由于其丰富的地表细节信息 , 使得外界条件变化 , 包括 想。根本原 因, 就是 阈值法只考虑了像元灰度这一单一属 太 阳光照射 、 阴影等因素 , 对其成像过程产生较大影 响 , 加 性, 忽略了其他诸多信息 , 所以这种方法在对高分辨率遥感 之 “ 同物异 谱 ” 和“ 异 物 同谱 ” 现象 的存 在 , 使 得 获取 较 好 的 影像进行分割时 , 难 以扩展到多波段彩色空间, 因此较少被 影像分割结果变得异常困难 。因此 , 通常的基 于物理模 型 采用 。 的方法 , 都要求满足一定的约束和条件 , 因此在应用中受到
2 0 1 6 年第六期 遥感测绘
WE S T R N R E S O U 只 C E S
西 部鸳
高分辨率遥感影像分割方法综述
傅 刚
内蒙古 自治区测绘产品质量监督检验站
呼和浩特 0 1 0 0 2 0
摘要 : 随着遥感技术的发展 , 遥感影像的空间分辨率不断提升 , 民用商业卫星已达亚米级 , 这给传统的图像处理方法 带来 了挑战, 于是提出了面向对象的分析方法 , 这种分析方法不再基于像元 , 而是基于由若干像元组成的对象 , 因此 , 面向 对象方法的基础 , 是要求准确地获取对象 , 而这一点依赖于影像分割技术 。本文调研整理了目前主要的遥感影像分割方 法, 并针对高分辨率影像的特点 , 对这些方法的适用性进行了分析说明 , 并指出了分割技术未来的发展方向。
dem多重分形的计算方法
dem多重分形的计算方法
DEM(Digital Elevation Model,数字高程模型)是描述地表地形和地貌的数字化模型。
多重分形理论是一种用于描述和分析具有多尺度特征的复杂系统的方法,可以应用于DEM数据的分析和计算。
下面是一种常用的计算DEM多重分形的方法:
1. 数据预处理:首先,对DEM数据进行预处理,包括数据的插值、平滑和去噪等操作,以确保数据的准确性和可靠性。
2. 分割数据:将DEM数据按照不同的尺度进行分割,可以通过使用不同的窗口大小或滑动窗口的方式来实现。
较大的窗口尺寸适用于较粗的尺度分析,较小的窗口尺寸适用于较细的尺度分析。
3. 计算局部均值:对于每个分割窗口,计算窗口内所有像元的高程值的均值。
这将产生一个局部均值图像。
4. 计算局部方差:对于每个分割窗口,计算窗口内所有像元的高程值的方差。
这将产生一个局部方差图像。
5. 计算分形维度:使用局部均值图像和局部方差图像计算分形维度。
一种常用的计算方法是通过计算局部方差与局部均值的对数之间的线性拟合来得到分
形维度。
6. 分析结果:通过分析分形维度的分布、尺度关系和变化趋势,可以揭示DEM数据的多重分形特征。
需要注意的是,DEM多重分形的计算方法可能因具体的研究目的和需求而有所差异,上述方法仅为一种常见的计算方式。
在实际应用中,还可以使用其他分形分析方法和工具进行DEM数据的多重分形计算和分析。
如何arcgis下用SHP文件切DEM和影像图
如何arcgis下用SHP文件切DEM和影像图转自:/s/blog_6758f2cf0100n0e9.html ArcGIS9软件,ArcToolbox -->Spatial Analysis Tools-->Extraction--->Extract by Mask(记得:必须安装SPatial Analysis扩展模块;地形图shp文件需要是面状;空间参考得和你得DEM用一样得,如果不一样,要用PROJECT先转换~。
)---ok(下面情况未验证~!)栅格数据的剪切可以用Arcmap spatial analysis 工具中的Raster caculation 首先在spatial analysis中的option中选择要剪切的范围(可利用矢量图层作为掩膜),计算中进行+0运算,就可以得到所要范围的栅格图像。
-------------------------------------Raster:栅格数据DEM:数字高程模型(ASCII)Image:图片Grid:grid是ARCGIS专有一种栅格数据类型。
---------------------------------------ArcGIS中对矢量和栅格数据进行裁剪切割的方法1、对矢量数据进行裁减: Arctoolbox中,spatial tool > extract>clip在InputFeatures中选择被裁剪的图层,在ClipFeatures中选择裁剪形状//2、利用矩形对栅格数据裁剪:Arctoolbox中,data management>raster>clip InputRaster中选择被裁剪的栅格数据,设定好矩形四个顶点即可//3、在Spatial Analyst Tools中提供了多种对栅格数据的提取方法:Arctoolbox中,Spatial Analyst Tools>extraction包括提取值到点,根据属性提取,用圆提取,用多边形提取,用掩模提取,用点提取,用矩形提取等,其中的用掩模提取功能可以让我们通过不规则边界来获取需要的栅格数据。
到手的DEM不会用?教你6个常用强大功能
到手的DEM不会用?教你6个常用强大功能一、概述DEM是数字高程模型(Digital Elevation Model)的简称,接触GIS,规划,设计类的多多少少会接触到DEM,可是这个直接查看黑溜溜一片DEM到底可以用来做什么呢?二、背景在大多应用中DEM只是是原始数据,想要让它具体发挥价值的话,就必须要对其进行分析以及处理操作,很多人第一反应是使用ArcGIS 进行分析,然后上网查各类ArcGIS的相关教程,并且如果本地没有ArcGIS的话,还需要花一些功夫去找软件进行下载,安装以及激活......这些都会花费不少的时间,并且对于一些低配机子来说,含有众多工具的ArcGIS进行DEM数据处理的效率还很低。
这时候我们最想要的效果就是,用轻量化的软件以最少的步骤来最快的获取所需要的结果,对于DEM的处理来说,Global Mapper(以下简称GM)就是一个较好的软件。
那么GM到底如何处理DEM呢?网上的教程太少?官方的帮助文档是英文?没关系!本文就会对6个功能点结合实际应用进行详细的讲解呢。
全球5米12米30米DEM样例数据免费下载:http://suo.im/5YUCpb三、功能介绍1. 晕渲图网上下载或者项目发来的DEM,直接看上去是黑漆漆的一片,完全看不出地形的样子,这类数据到底如何查看呢?黑漆漆一片的原因是因为,tif默认识别是图片,用图片浏览器查看dem由于是按图片的编码进行展示的,所以看不到具体的东西,展示出来就是一片黑。
这里就需要使用GIS类软件对dem进行加载了,直接将tif拖入GM即可看到如下效果:这么一看是不是就清楚多了,地形的高程以不同颜色显示着,一眼就看出了地形的样子。
这个其实是dem的晕渲图效果,在工具栏中我们也可以根据不同的晕渲效果制作不同的专题图。
这样的话我们就可以根据不同的需求生成不同的晕渲图,并将其进行成图进行展示,如下我们选择边坡渲染效果,即可将坡度信息以可视化的形式直接显示出来。
使用Photoscan生成DEM与正射影像流程
使用Photoscan生成DEM与正射影像流程(使用像控点)1.参数预设使用工具菜单的工具-偏好设置打开PhotoScan Preferences对话框一般(General)选项卡上的参数设置下列值:立体模式:浮雕(如果你的图形卡支持四轴缓冲,使用硬件)视差:1.0将日志写入文件:指定Agisoft日志的目录GPU选项卡设置如下:勾选在对话框中PhotoScan检测到的任何GPU设备。
当使用少于两个GPU时,勾选“在执行GPU加速时使用CPU”高级选项卡参数设置下列值:保持深度图:启用存储绝对图像路径:禁用启用VBO支持:启用2.添加照片从工作流菜单中“添加照片”选择添加照片命令或单击工作区工具栏上的Add Photos按钮。
在添加照片对话框中浏览源文件夹并选择要处理的文件。
点击打开按钮。
3.装载相机POS文件生成的模型使用的坐标系统是由这个步骤中设置的相机POS坐标系统决定的。
如果相机位置未知,这一步可以跳过。
对齐照片这种情况下需要更多的时间。
打开视图菜单中的参考面板,在参考面板工具栏上单击“导入”按钮,并在打开的对话框中选择包含POS信息的文件。
最简单的方法是载入字符分隔的文本文件(每张照片的x-和y坐标和高度(相机方位数据,即俯仰、滚动和偏航值,也可以导入,但数据不是必须)。
然后单击参考窗格中的Settings按钮,在参考设置对话框中选择相应的坐标系统,并根据测量准确度设置照片POS精度及标记、连接点、精度,如果没有在导入POS时指定坐标系,也可以在这个面板中指定坐标系。
地面高程:在倾斜拍摄的情况下,应该指定对应坐标系统椭球面上的平均地面高度。
点击确认后,相机位置会标记在模型视图中,如果在POS数据正确的情况下无法看到相机位置,点击工具栏中的显示相机按扭,然后点击工具栏上的重置视图按钮。
4.检查相机校准打开菜单栏“工具”-“相机校准”窗口。
默认情况下,Photoscan将在对齐照片和优化的过程中通过照片的相应参数自动估算矫正参数,如果相机相关参数缺失可以手工输入。
切割技巧详解 Photoshop中的裁剪和修剪工具
切割技巧详解:Photoshop中的裁剪和修剪工具Photoshop是一款强大的图像处理软件,广泛应用于设计、美术和摄影领域。
在使用Photoshop时,切割是一项非常基础且常用的技巧。
本文将详细介绍Photoshop中的裁剪和修剪工具,帮助读者更好地掌握这一技巧。
一、裁剪工具裁剪工具是Photoshop中最常用的切割工具之一。
主要用于裁剪图片的尺寸和比例。
下面是使用裁剪工具的步骤:1. 打开你需要裁剪的图片。
2. 选择工具栏上的裁剪工具(快捷键C)。
3. 在图片上拖拽鼠标,选中你想要保留的区域。
可以根据需要调整选中区域的大小和位置。
4. 单击“确定”按钮,或按下回车键,即可完成裁剪。
在使用裁剪工具时,你还可以通过调整选中区域的边缘来改变裁剪的尺寸和形状。
例如,将鼠标移动到选中区域的边缘,会出现一个双向箭头的图标,表示可以拖拽该边缘来调整裁剪区域的大小。
另外,裁剪工具还提供一些选项供你选择,如保留裁剪区域内的隐藏信息或图层。
这些选项可以在裁剪工具的设置中进行调整。
二、修剪工具修剪工具是另一种常用的切割工具,用于去除图片中的多余部分或不需要的区域。
下面是使用修剪工具的步骤:1. 打开你需要修剪的图片。
2. 选择工具栏上的修剪工具(快捷键C)。
3. 在图片上拖拽鼠标,选中你想要去除的区域。
可以根据需要调整选中区域的大小和位置。
4. 单击“确定”按钮,或按下回车键,即可完成修剪。
与裁剪工具类似,修剪工具也提供了更多的选项供你调整,如保留修剪区域的隐藏信息或图层。
这些选项可以在修剪工具的设置中进行设置。
三、其他注意事项在使用裁剪和修剪工具时,你还需要注意一些细节,以便获得更好的效果:1. 备份原图:在对图片进行切割之前,最好备份一份原图。
这样,即使切割效果不理想,你还可以随时恢复到原图状态。
2. 保持比例:如果你需要保持图片的宽高比例不变,可以按住Shift 键,然后调整选中区域的大小和位置。
这样可以确保裁剪或修剪后的图片比例与原图一致。
dem数据使用教程
DEM高程数据(2013-11-12 15:06:40)转载▼标签:杂谈DEM高程数据包括两个部分:ASTER GDEM30米分辨率高程数据和SRTM90米分辨率高程数据。
ASTER GDEM数据来源于NASA,数据覆盖范围为北纬83°到南纬83°之间的所有陆地区域,时间范围为2000年前后;SRTM数据来源于CIAT,数据覆盖范围为北纬60°至南纬60°之间的所有陆地区域,时间范围为2000年前后。
ASTER GDEM 30米分辨率高程数据本数据集利用ASTER GDEM第一版本(V1)的数据进行加工得来,是全球空间分辨率为30米的数字高程数据产品。
由于云覆盖,边界堆叠产生的直线,坑,隆起,大坝或其他异常等的影响,ASTER GDEM第一版本原始数据局部地区存在异常,所以由ASTER GDEMV1加工的数字高程数据产品存在个别区域的数据异常现象,可以和全球90米分辨率数字高程数据产品互相补充使用。
ASTER GDEM数据采用UTM/WGS84投影,数据格式为IMG栅格影像,数据的值域范围为-152-8806米之间,比例尺为1:25万,其垂直精度20米,水平精度30米。
数据命名规则:ASTER GDEM基本的单元按1度X1度分片。
每个GDEM数据包有两个文件,一个数据高程文件和一个质量评估(QA)文件。
每个文件的命名是根据影像几何中心左下角的经纬度产生。
例如,ASTGTM_N29E091代表左下角坐标是北纬29度,东经91度。
ASTGTM_N29E091_dem和ASTGTM_N29E091_num对应的分别是高程数据和质量控制数据。
SRTM 90米分辨率高程数据SRTM(ShuttleRadarTopographyMission)90米分辨率高程数据由美国太空总署(NASA)和国防部国家测绘局(NIMA)联合测量。
2000年2月11日,美国发射的“奋进”号航天飞机上搭载SRTM系统,共计进行了222小时23分钟的数据采集工作,获取北纬60度至南纬60度之间总面积超过1.19亿平方公里的雷达影像数据,覆盖地球80%以上的陆地表面。
DEM遥感图像合并裁剪
如何在ERDAS、ARCgis中实现用矢量线界裁切遥感影像总结一下利用Erdas和Arcgis来随意图形分割影像图:影像图格式为tif随意图形格式随意(就当shape格式)一:合并影像图:由于影像图的分割需要,则要全部覆盖shape格式的边界。
分Erdas和Arcgis两种合并法:Erdas合并:①打开Erdas,再打开viewer窗口(注意:再打开图层时,要将raster option的no stretch 和background transparent前打勾,这样图层就不会失真)。
②打开需要合并的图:files of type选择,选择对应的的图层③在窗口viewer—raster—mosaic images,弹出mosaic tool窗口④在窗口mosaic tool—process—run mosaic,弹出窗口⑤在窗口output file nam选择files of type为tif格式,存放路径自己选择,点击ok完成Arcgis合并: ①打开Arcgis(Arcmap), 加入要合并的图层②ArcToolBox—datamanagementtools—raster—mosaic或mosic to new rastevr然后按需求选择input raster 和 output raster还有名称。
二:生成分割边界:①先用arcmap打开shape图,选中shape边界②ArcToolBox—conversion tools—to coverage—feature to class coverage然后按需求选择input feature classse和在output coverage填入存储路径,点击ok完成,生成coverage格式。
③然后在Erdas用窗口viewer打开coverage图和tif图。
④再矢量图上用鼠标点击你要裁切的边界线,再点击“AOI”菜单,选择 copy select toaoi; View—>Arrange Layers Viewer打开Arrange Layers Viewer对话框,在Vector图层上单击右键,选择Show Properties,打开Properties对话框,选中Polygon,点击Apply按钮,再选中矢量图中重叠的部分。
worldwind切片规则
WorldWind是一个三维地球模拟器,它包含遥感影像、DEM及矢量国界线三种图形图像文件。
下面分别介绍各类型文件的分片及存储目录规则。
1. 遥感影像:第0级块的大小为
2.25度。
每增加一级,大小减小为原来的1/2,总共5级。
存储结构为——层级/行数/行数_列数.dds/jpg,每张图片大小均为512*512。
2. DEM:第0级块的大小为18度。
每增加一级,大小减小为上级的1/2,总共12级。
存储结构为——层级/行数/行数_列数.bil,每块均为150*150的高程矩阵。
3. 矢量国界线:第0级大小为36度。
每增加一级,大小减小为上级的1/2,总共6级。
存储结构为——层级/行数/行数_列数.dds/png,每块图片大小为512*512。
以上信息仅供参考,如需了解更多信息,建议查阅WorldWind 的官方文档或咨询相关技术人员。
WorldWind系列十三:地形数据(DEM)加载和应用(入门篇)
WorldWind系列十三:地形数据(DEM)加载和应用(入门篇)来源:博客园作者:无痕客学习WW中DEM的加载和应用对于基于WW搞三维GIS插件开发很重要。
我一直也想关注和搜集相关学习资料,之前没有开始研究它,主要是自己没能在WW看到高程应用功能,感觉没法切入。
我昨天下午和今天上午的时间就花在找研究切入点啦!实质上就是找如何在WW 查看高程的功能。
(现在回想一下,感觉浪费时间了,还是因为我对WW的功能没能完全知道)我在开始学习DEM加载和应用前,我主要参看了CSDN上的博文——《庆贺,绘制出World Wind的LOD地形网格》(/paul_xj/archive/2007/08/29/1763585.aspx),相信研究WW 的人都应该看过这篇和他的三篇学习总结。
我也看过研究过好多遍,结合自己的学习,每遍拜读都有些收获。
他的文章关注度很高,绝对是WW初学者的起航灯塔。
但是,我对上面提到的《绘制出World Wind的LOD地形网格》几乎要全文否定。
我不希望网友以该文为标准来学习WW 的地形网格(即DEM加载应用),因为我原本也是很相信他说的,按他说的去做的,事实不是那样的,相信他也是刚入门时写下的,很多DirectX的东西理解也不是很深刻的。
所谓的绘制出LOD地形网格有问题!怎么能将DirectX三维渲染与GIS中DEM混为一团?!地形图包括地物和地貌。
而地貌主要体现在等高线上,在三维中体现为DEM应用(数字高程模型)!上面博客中讲到的构建地形网格只是将DirectX的面渲染改为线渲染,这还是平面级别的,根本没有体现高程,怎么就是地形网格啦?!从DirectX三维技术角度讲,也是有问题的:DX中面渲染和线渲染存储点的方式是不同的,因为WW中使用LOD模型各层要渲染的面很多且相邻,所以简单更改渲染方式没能看出影响。
但如果你绘制单个面然后更改为线性渲染,就会看出问题啦!这还是平面级别上面,面渲染改为线渲染会出现问题,如果是三维物体(立方体)如果简单改为线渲染,问题更明显!(你学过DirectX后,试一下就知道啦!)至于博客中说到的:“按下…Ctrl+W‟出来的网格,是什么空中网格”,只是简单地操作,放大后看到后面是“天空”,就认为是什么空中网格,这真是典型的形而上学啦!真正学习过三维开发的人都会知道:为了提高效率,三维渲染只是渲染可视的前一面,不被看到的面统统隐去(不渲染)。
影像处理和配置
Sgis平台中影像图层基本处理思路是通过arcgis软件分级切片生成的金字塔结构的类似:其中每个小图片的大小为512*512,将每个图片网格化,平均分成32*32个网格,共需要33*33个点,并求每个点对应的高程值,将高程值保存到一个文件内,因此每个影像文件对应一个高程数据文件,高程数据组织方式和上面的一致,下面是详细的处理步骤:一、首先利用arcgis工具对影像进行分级切割,要求:1.相邻level的比例尺必须是2倍关系2.level 0的切片必须是一张图。
3.切图参考点在地图的左上方,即xmin ,ymax4.生成的图片指定为png32位的,要有alpha通道信息5.生成的图片大小必须是512*512二、配置影像图层的元数据信息在用arcgis工具分级生成影像后,会生成conf.cdi conf.xml两个文件,这两个文件记录了切割的相关信息。
其中conf.cdi中类似<XMin>349025.62499989558</XMin><YMin>3590592.8750000009</YMin><XMax >439756.37499989558</XMax><YMax>3642389.1250000009</YMax>键对记录的是影像的范围,这个范围可能比实际存在影像数据的范围大的多,如下图:。
conf.xml中<TileOrigin xsi:type="typens:PointN"><X>349028</X><Y>3642387</Y></TileOrigin>记录的是切割参考点的坐标;<LevelID>0</LevelID><Scale>1024000</Scale><Resolution>270.93387520108377</Resolution>这个是每个level对应的比例尺信息。
测绘技术的影像处理方法详解
测绘技术的影像处理方法详解引言:测绘技术是现代社会中至关重要的一项技术,涉及到地理信息的收集、处理和分析等方面。
而影像处理则是测绘技术中重要的一环,通过对图像的处理和分析,可以提取出有关地理信息的关键内容,为地图制作、资源管理和环境监测等领域提供数据支持。
本文将详细介绍测绘技术中常用的影像处理方法。
一、图像预处理图像预处理是影像处理的第一步,目的是提高图像质量,为后续处理做准备。
其中,图像去噪、图像增强和图像分割是常见的预处理方法。
图像去噪是指去除图像中的噪声信号,提高图像的清晰度和细节。
常用的去噪方法包括中值滤波、高斯滤波和小波去噪等。
中值滤波是一种空间滤波方法,通过计算像素周围邻域内像素的中值来代替原像素值,从而消除孤立噪声对图像的影响。
高斯滤波则利用高斯核来平滑图像,使噪声信号得到衰减。
小波去噪则利用小波变换将图像分解成不同频率的子带,通过去除低频子带中的噪声信息来实现去噪效果。
图像增强是指通过图像处理方法来改善图像的视觉效果,使图像更易于人眼观察和分析。
常用的图像增强方法包括直方图均衡化、灰度变换和滤波增强等。
直方图均衡化通过对图像像素值的统计分析,将图像的灰度级均衡分布,增强图像对比度。
灰度变换则通过对图像的灰度级进行非线性变换,调整图像的亮度和对比度。
滤波增强则通过选择合适的滤波器,增强图像的边缘和细节。
图像分割是将图像分成若干互不重叠的区域,每个区域具有相似的像素值和纹理特征。
图像分割是测绘技术中一项重要的任务,可以提取出有关地物边界的信息。
常用的图像分割方法包括阈值分割、边缘检测和区域生长等。
阈值分割是指通过设置合适的像素值阈值,将图像分成两个区域,分别对应地物和背景。
边缘检测则通过检测图像中灰度级变化较大的区域,将图像分成不同的边缘区域。
区域生长是指从某个种子点出发,将与之相邻且具有相似像素值的像素逐步加入到同一个区域。
二、图像配准图像配准是指将不同影像之间的几何校正和坐标统一,使其之间具有空间对应关系。
ENVI中裁剪DEM
从数字高程图上导出高程数据步骤:1.提取高程数据使用的软件:ENVI 4.52.装入数字高程数据(以现有的dem90-utm.img为例)单击“File”选择“Open Image File”打开高程数据双击“Available Bands List”窗口中的“Layer1”便会出现如下界面:上图三个灰黑色的图框内即为全疆的数字高程数据,只是放大比例不同,在双击图中红色方框后(上方的图框)即可产生带有经纬度坐标的对话框,此时经纬度坐标会跟着鼠标在图上的移动而相应变动(经纬度的显示可能会不正常,请到网上寻找相应的补丁进行安装)。
我们依据经纬度坐标来划定我们想要提取数据的范围。
3.数据提取操作在“Basic Tool”中选择“Region Of Interest”中的“ROI Tool”,此时我们便可以使用鼠标在灰黑色图框中进行区域的划定。
首先在“Scroll”窗口中用鼠标单击红色方框,按住左键拖动方框选择大致的区域位置。
再在“ROI Tool”中点选“Image”选项。
此时便可根据区域的角店坐标在图框中找寻并确定;在确定出大致的区域后单击鼠标右键得到一个封闭的区域。
再次单击右键。
而后选择“ROI Tool”窗口下的“File”,并选择“Subset Data via ROIs ”我们即得到了大致的区域。
在新得到的区域中继续用鼠标选取较为精确的区域(请根据需要将窗口进行放大再选取),注意此时仍然需要用到“ROI Tool”。
如认为第二次选取的区域不够精确可再选取好的区域中间继续进行更加精确的区域选取和保存。
区域选取和保存的步骤在此不再重复。
4.数据输出其中XY表示公里网坐标,Z表示相应的高程。
文件以txt格式保存,可以直接由mapgis打开经行投影变换或经行离散数据网格化。
浅谈1:2000DEM的制作与裁切
介 绍 了制 作 过 程 中需要 注意 的 一 些 问题 。 关 键 词 : E Ma e; rM p D M; 据 转换 D M krA c a ; E 数
为推进城市地理空间框架建设 , 更好地将测绘工作与 地方人 民政府 的工作需求相结合 ,为 国民经济 建设 服务 , 推进数字中国地理空间框架建设进程 , 国家测绘局正式将 本溪 市列入 20 0 9年数 字城市 地理空 间框架 建设试 点城 市 。我单位有幸参与数 字本溪 12 0 D M、 O :0 0 E D M的制作 。
结束 语 我们通过从 D M kr A c p 到 v t z 一系列过程 , E Mae 到 rMa 再 iu o ro 完 大于两米 的宅基地要用特征线表示完整 , 这样表示 即可将地型完整 E 这样可 以通过这些软件使 D M E 表示 , 又不影响作业精度 。本次作业精度高 , 要求在 2米以内 , 因此 成了 D M 的制作检查和数据转换 , 这样更适合使用者的要求 。 此次 D M 编缉特别细质 ,对宽度大于两米 的路 、水要求用双线表 数据实现各种格式的转换 , E 参 考 文献 示, 高差大于两米的坎均应表பைடு நூலகம் , 以此次作业的数据量 特别 大 , 所 比 l 1 A c I M] 北 清华 大学出 以往 D M制作有一定的难度 。在一个作业 区内所有人按分幅 图制 f 吴秀芹.rG S9地理信 息 系统应 用与 实践[ . 京: E
表示完整 , 由于当地地形特殊 , 民地依山而建 , 并且 居 因此对于高差
20. 作完后 , 出 D F 在 C D下进拼接 , 导 X, A 方便进入 A c I rGS下进行裁图 版社 0 9 f 李志林 , 2 1 朱庆. 字高程模 型『 . 数 M1 第二版 . 汉: 武 武汉 大学 出版社 , 和数据转换 。 20 . 0 0 3在 A e p软件 中对 D M数据进 行数据转换和裁切 rMa E
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
WorldWind系列十五:如何切割影像和DEM数据及其在WW中的应用配置来源:博客园作者:无痕客WorldWind学习系列十四中我从代码上分析如何加载DEM数据,里面涉及了算法,有学习和借鉴意义。
但对于只求应用或者说是急于求成的网友来说,实用价值不是太大!我们分析代码是一种学习过程,不是目的,终究要落在如何应用自己的影像和DEM数据、如何基于自己的数据开发满足自己项目需求的新的WW插件。
这是一过程,需要循序渐进,今天只跟大家分享一下如何切割自己的影像或DEM,及如何设置相应的XML配置。
我的学习和实践,主要是参考/wiki/Dstile_howto,对于英语较好的网友,建议阅读原文,反正我英语不太好,结合实践,我差不多看了五六遍原文。
我下面是翻译部分内容及自己的实践步骤和心得(将采用红色字部分)。
这是一篇教你如何一步步将你自己的具有地理参考系的影像数据添加到WorldWind中的文档。
该帮助文档将告诉你如何将影像数据切割成瓦片,及怎样在WW中利用这些瓦片数据创建图层。
至此,你可以将该图层作为一个Add-on插件来发布。
这个图层可以完全放在你的电脑上。
如果你想将瓦片数据作为发布服务提供给他人,该文档可能对你帮助不大。
WW的瓦片系统的组织如下图:阅读下面的内容前,请确保你已经理解上图中WW瓦片的工作原理和XML的配置原理。
(注:上图是WW的精华)下面我将带你创建一个图层实例,里面的数据可以换做你自己的。
(注:第一次实践的话,建议你使用同样的数据同样的步骤,将整个过程走一遍)必备软件工具:FWTools:这是一个操作具有地理参考系影像的免费的工具包,是开源工具。
dstile切图工具需要使用FWTools,最新版本V2.4.6。
下载地址:/安装FWTools时,你可以选择不安装开发文件,这样可以节省一些安装空间。
dstile工具与2.X版本兼容,与1.X版FWTools不兼容。
Dstile:使用版本是2007年8月13日,下载地址::8080/tisham/dstile-0.2_win32_whatnickpatches.zip这个程序按照Making Layers文章的处理过程。
请将压缩文件解压后拷贝到FWTools的安装目录下(如:D:\Program Files\FWTools2.4.6)WorldWind(1.4版)前期准备:除了软件,你还将需要一个具有参考系的影像数据。
本文中,我将使用特拉华州北部的地质图。
(下载:/publications/digitaldata/geologicdata.aspx文中使用的是第六个TIFF数据) 影像格式必须是FWTools所支持的(here)。
Dstile目前不能处理索引影像,但是这样的影像可以先通过FWTools中pct2rgb功能转换成RGB影像。
简便起见,我会将影像数据放到FWTool文件夹下,但也可放在其他地方,使用时确保正确的影像文件路径。
Dstile 也可以切合地形数据,只需要在命令行中添加--dem参数。
注意:地形数据必须是16位的。
创建瓦片金字塔:Dstile是个命令行程序,在你的开始菜单中运行“FWTools shell”(注:安装FWTools 会在桌面上出现该快捷图标)。
将会出现带FWTools路径的命令行窗口。
你也可以在FWTools 文件夹下运行“setfw.bat”命令脚本。
创建瓦片的通用命令行格式如下:dstile.exe tile --lztsd tile_size --wwcache --overviews output_directory georeferenced_image_file(s)参数如下:tile_size - 0图层瓦片大小,十进制数,能被180正常(最好为180/(2^N))。
output_directory - 瓦片金字塔的输出路径。
(如何没有盘符只要文件夹名,是在FWTool安装目录下)georeferenced_image_file(s) - 影像数据的完整路径(如果影像数据在FWTool安装目录下的文件夹下,可以使用相对路径).一系列的多重影像可以自动组合成单一的瓦片金字塔。
( A list of multiple image files can be supplied and will be automatically combined into a single set of tiles.)命令说明:--tile :告诉dstile将要创建瓦片金字塔--wwcache :告诉dstile 按World Wind 的内部格式组织和命名创建的瓦片--overviews :告诉dstile 去创建所有可以创建的瓦片图层,忽略该命令,将会导致只创建最详细的瓦片图层--png:通常用来创建.png格式的瓦片(速度慢且没压缩),忽略该命令,将默认为jpg格式瓦片--float :被用来切割浮点型DEM数据,比Int16稍微更准确些。
--srcProj: 被用来指明缺少.tif/.tfw文件时的投影坐标。
使用epsg代码。
例如 epsg:32754 is UTM Zone 54S , epsg:4326 is WGS84.学会在FWTools Shell中查看命令:我示例中的地质图,我使用如下命令创建瓦片金字塔:dstile.exe tile --lztsd 1.0 --wwcache --overviews output_tiles geomap10.tif 将会在FWTools的安装目录下创建一个新的output_tiles 文件夹,里面存放着瓦片金字塔系统。
输入影像文件为geomap10.tif。
注意:这个输入文件必须是影像数据自身,而不是世界文件(示例中的geomap10.tfw )尽管.tfw文件是必须的且必须与影像数据在同一个文件夹下。
在这个例子中,output_tiles 文件夹下创建了六层数据(文件夹名字从0到5)。
Dstile运行中,会依次读取输入参数,并显示进度。
处理该111M的影像数据花费了一两分钟。
据说dstile在Free Earth Foundation网站的应用中,单机从21GB中录制1.8G的瓦片树花了70分钟。
这是我使用自己的试验数据生成的瓦片金字塔,为9层。
0层瓦片大小(Level zero tile size):这个参数控制着你的影像图层在哪一缩放界别或视角高度上显示。
0层瓦片大小越大说明你放大时越早出现,同理相反。
作为参照,the blue marble的lzts为36°,WorldWind中默认的Landsat图层的LZTS为2.25°。
Zoomit!功能中的高分辨率影像使用的LZTS为0.1~0.2度。
在地质图中,我使用1°的LZTS。
这使得该图层将会在海拔390千米左右出现。
你也可以按照自己的需要改变该参数。
LZTS的大小理想值为最好为180/(2^N),N为整数。
该例子中理想的LZTS为1.40625(180/(2^7))。
如果dstile的切割值不满足该表达式,可能会导致图层在WW中绘制是出现水平定位问题。
本实例中的地图没有表现出该问题是因为相对粗略的比例尺。
高精度的正射影像更可能会出现水平定位问题,所以最好是遵循该表达式180/(2^N)规则。
向WW中添加图层:为了在WW中显示该图层,需要将切好的瓦片数据放到WW的目录下,同时要为该新图层创建一个XML配置文件。
本实例中,我将瓦片数据放在\Data\Earth目录下。
同样,瓦片数据可以在/Add-ons/的任何位置。
我没有将数据放到\Cache\目录下,是因为我偶尔会删除缓存文件夹下所以数据。
我会在\Cache下存放可以重新下载的数据。
因此我将output_files文件夹下的所有内容放在WW 目录的 "\Data\Earth\N. Delaware geology map\"里。
这些瓦片文件在WW目录下的组织方式由你来定。
如果我有很多地质地图放在WW离,我可能会创建个文件夹然后把瓦片数据放到里面。
(例如:\Data\Earth\Geology maps\N. Delaware\ )只要你能让WW按照你的组织方式找到它们。
一旦瓦片数据放好位置,我们需要为该图层创建一个XML文件。
该XML文件存放在WW中的\Config\Earth\。
将XML命名为图层名是很好的方式。
我使用的XML内容如下:<?xml version="1.0" encoding="UTF-8"?><LayerSet Name="Geology maps" ShowOnlyOneLayer="false" ShowAtStartup= "false" xmlns:xsi="/2001/XMLSchema-instance"xsi:noNamespaceSchemaLocation="LayerSet.xsd"><QuadTileSet ShowAtStartup="true"><Name>N. Delaware</Name><DistanceAboveSurface>0</DistanceAboveSurface><BoundingBox><North><Value>40.0</Value></North><South><Value>39.0</Value></South><West><Value>-76.0</Value></West><East><Value>-75.0</Value></East></BoundingBox><TerrainMapped>true</TerrainMapped><ImageAccessor><LevelZeroTileSizeDegrees>1.0</LevelZeroTileSizeDegrees><NumberLevels>6</NumberLevels><TextureSizePixels>512</TextureSizePixels><ImageFileExtension>jpg</ImageFileExtension><PermanentDirectory>Data\Earth\N. Delaware geology map</P ermanentDirectory></ImageAccessor><TransparentColor><Red>0</Red><Green>0</Green><Blue>0</Blue>< /TransparentColor></QuadTileSet></LayerSet>" LevelZeroTileSizeDegrees"为你切割影像时的LZTS的值"PermanentDirectory" :是你存放瓦片的文件夹路径图层集合的名字定义为:"Geology maps",该图层的名字为 "N. Delaware".影像外部的部分空白的,对应的瓦片是黑色的。