方法一:由DWG地形图生成DEM
.由DWG地形图生成DEM
1.1从DWG中提取高程点数据
1.1.1切割DWG地形图
数据量太大,先切割再进行其他操作。具体步骤为:
用CAD2005把上、下两幅图转换成2000格式(CASS是CAD2002配套产品)-用CASS打开上、上两幅图(CAD中没有SAVET保存选择多边形内图形功能)-“插入”-“块”-名称中打开红线研究区-去掉“在屏幕上指定点”(X,Y,Z全是0)-确定后就可以显示红线研究区-用矩形圈出研究区-“SAVET命令”-输入比例尺(10 000)-多边形保存1-选中刚画的矩形-OK。
1.1.2合并上下两幅图
CAD中有一些命令,qselect可以选择满足条件的数据,就可以选择一层数据,wblock可以制作块保存选择的数据,具体步骤为:
打开裁剪后的图上-“插入”-“块”-打开裁剪后的图下-去掉“在屏幕上指定点”(X,Y,Z全是0)-选上左下角的“分解”(如果不分解,整个下图就是一块,选中一条线就把图下全部选中了,删除一条线就把整个删除了,当然现在不选,可以用CAD分解命令分解开)-确定后两幅图就拼接好了-然后打开红线-再次整体裁剪两幅合并的图-打开图层管理-只显示等高线和高程数据图层-另存为CAD图。
中国3S吧 3s8
1.1.3补充高程点数据
由于等高线质量太差了-断线或缺少线,没有高程属性等,不用等高线生成DEM,用高程点数据生成DEM)。具体步骤为:
设置文字样式通过“格式”-“文字样式”-设置和原来的高程文字相同样式-补点用TEXT命令-用鼠标确定文字位置-确定角度为0-输入高程数据-复制高程数据文字-沿着等高线粘贴该高程数据即可(以后用回车或空格完成粘贴)-换等高线时粘上错误高程后双击文字可改-然后再复制新文字
1.1.4获得高程点数据表
原先已有高程点是由“高程点和高程数据注记文字”组成的,高程点提供了准确的位置(X,Y)而没有Z属性,但文字注记提供了高程值而位置是不准的,有一个解决办法可以得到准确位置的准确高程值,先得到所有点的位置数据表(包含X,Y),再得到高程数据表(包含X,Y,H),再编程实现点和高程值的匹配,具体实现方法为:点的位置数据和高程数据分别保存在两个数组中,从第一个点开始在高程数据中距离与他小于一个定值的高程文字,这个文字的内容就是这个点的高程,到后马上去掉这个高程文字数据,减小以后的寻负担(在VC中可以用CUintArray作为数据数组,有删除函数,采用GetSize()得到要寻的数据个数;当然还有一种方法是,现在已经有EXCEL数据,转换成ACESS数据库,然后在VC中读取数据库,一个在VC中好实现读取ACESS数据库,再一个是不是速度比VC中读取EXCEL文件快呢?具体实现时在点数据表中新那一个字段,保存高程,在另一个高程表中读取XY值比较距离,打到高程就把高程数值更新到点数据表中的新字段中,当然到一个就把高程表那一条记录删除,当然到一条记录最好是再接着,要是到两个就说明那附近有问题,一个点和两个高程数据接近,或者说没有到任何一个点,是不是距离设置太小了。可以通过VC,VB访问数据库,当然也可以在ACESS的VBA中使用),开始没有做点和高程的
匹配,只是把高程数据文字的位置当作高程点的坐标了,在CAD图上看了一下,一般高程点和高程文字注记的距离为30多米。不过,后来把研究区分解成四部分在Excel中根据阀值和最小距离实现了坐标和调和的匹配,在测绘通报和其他测绘方面的期刊上有这样的论文。在这里提取文字信息也是在明经CAD论坛中到VBA的代码的。获得高程点数据的具体步骤为:
专业的3S站 3s8
“工程”-“宏”-“加载工程”-加载“提取文字信息”VBA代码-到宏中的VBA编辑器运行就可以了(也许需要添加引用EXCEL),结果保存在C盘下的EXCEL文件xyz中(他的X和Y和CAD图上是反的)。
输出距离小于100米内的点个数,这样避免下面这样的情况,两个都可以。设置上限100米,这样可以分析是不是对了
.235
125.
1.2.根据高程点数据生成DEM
1.2.1生成点SHP文件
ArcMap的工具添加XY数据至少可以打开ACCESS数据表和TXT文件,ACCESS文件数据表直接打开没问题,TXT文件的格式为:第一行为属性字段,如(ID,X,Y,Z),然后从第二行开始就是数据了,如(1,12.45,23.4,234.5\n2,45.23,45.5,)。具体步骤为:
通过ACEESS创建数据表,新建数据库-附加数据库-打开CEXCEL高程xyz数据文件-根据提示可以生成mdf数据表-运行ArcMap-“工具”-“添加XY数据”-选择刚生成的数据表-X,Y字段选择位置-OK-然后保存为shp文件,在图层列表选择刚打开的点层-数据-导出为SHP格式-到此点SHP文件创建完毕OK。
1.2.2生成TIN数据和栅格
有的机器上装的ArcMap中的3D分析可以用,但是我这台机器上装的不能用,不过,我发现我的机器上ArcScene是好使的,幸亏能用,呵呵。那么具步骤就是:
中国3S吧 3s8
运行ArcScene-调出3D分析工具(“工具”-“扩展”中可,也可右键工具栏打开3D分析工具,当然还可“工具”-“定制”)-3D分析工具中创建/修改TIM-从要素生成TIN弹出对话框-打开高程点SHP层-不用修改了其他的(高度源就是高程字段)-确定就OK了。
转换到栅格的方法:3D分析-转换-TIN转换到栅格OK了。
从网上gissky上看到的“利用等高线生成DEM,最好还要有高程点数据,生成方法最好不要用构TIN的方法,要用ARCTOOLBOX-空间分析-内插-TOPO TO RASTER 命令,效果要好得多。此命令是专门用于DEM生成的内插方法”,做了第一个插值的,效果好像是好点,山顶没有构成TIN的图尖锐。
1.2.3DEM数据的显示设置
如果觉得高程差别太小或者太大,起伏不明显或玄虚太大,设置基高也许会帮助你解决视觉上的难题,按以下步骤做:
右键图层名-属性-基表面高度(倒数第2个)-选中“从表面为图层获得高度(这个栅格图像起始没选中,tin是选中的)”-Z单位转换自定义后面的数字就是转换因子,数字越大,高程差别
越明显。
栅格起始是黑白来区分高度,tin是用起伏来表示高低,可以设置由高到低用不同的颜表示,两种数据设置有点不一样的。具体步骤是: 专业的3S站 3s8
tin的设置:右键打开属性-符号-显示下面的“添加”添加渲染-可以选第5个“面高程用颜梯度进行渲染”-添加,取消关闭对话框-然后去掉Faces前的勾-然后就可以在“颜梯度”中选择自己喜欢的颜梯度了,还可以在右边设置分级数。
栅格可以直接点击图层上的颜,设置梯度颜,也可以打开属性和tin一样在符号中设置。
1.2.4DEM与遥感影像的叠加显示
在ArcScene中打开DEM,再加载影像时,一定要注意影像边界不能超过DEM,而且影像必须是GRID栅格格式,一个办法是利用ArcToolbox中的栅格剪裁工具根据DEM边界裁剪,然后在属性中指定基调文件为DEM就可以了
方法二:
DEM可以通过点,等高线,TIN等通过插值生成.
以点为例:
(1) 在Spatial Analyst下拉菜单中选择Interpolate to Raster,在弹出的下一级菜单中点击Inverse Distance Weighted命令, 弹出IDW对话框。
(2) 在Input points的下拉菜单中选择被用来进行插值的离散点数据;
acess数据库(3) 在Z value field的下拉菜单中选择要加入的字段;
(4) 在Power栏中填入进行插值计算的幂值;幂值就是距离的指数。如幂指数为2时则进行反向距离平方插值。幂指数是一个正实数,其缺省值为2。 12
(5) 在Search radius type 栏中选择一种搜索半径设置类型;
1) Variable:当选择此项时,搜索半径由下面两Maximum distance。首先在Number of points中输入搜索的最近点的个数(缺省值为12),然后在Maximum distance中输入一个控制距离。如果最近点的个数超出控制距离,则将会以控制距离为限制来选取较少的点;
2) Fixed IDW对话框
来控制,Distance和Minimum number of points。首先在Distance中输入搜索半径距离(缺省值是输出栅格大小的五倍),然后在Minimum number of points中输入控制插值点个数的最小整数值。如果搜索半径距离内的点个数小于插值点个数的最小整数值,则搜索半径自动增大。
(6) Use barriers polyline为可选项,输入中断线文件。barriers是在插值中,如有某些地方出现异常,(如某些断裂带),而要求插值时考虑到这样的因素,所设置的选项。它是一个打断表面的线特征。这一线特征没有Z值。悬崖,峭壁,堤岸或某些障碍都是典型的barriers。barriers限制了插值计算,它使得计算只在线的两侧各自进行。而落在线上的点则会同时参与线两侧的计算。
(7) Output cell size:指定输
(8) Output raster:为输出结果指定目录及名称;
(9) 点击OK按钮。
基于GE的dem生成方法:
1.将采样点数据存为Excel格式.
2.Arcmap中,Tools--add xy data,将Excel加载进去,以经纬度为xy值,生成点状图层.
3.打开3D analysis工具,creat TIN,然后convert TIN to raster,生成DEM.
4.在ArcScene中,打开生成的DEM和照片,通过联合高程信息将照片覆盖在DEM数据层上(右击照片图层-属性-base heights-obtein heights for layer from surface 选择DEM).
上面两种方法是我在网上看的,具体操作我也不知道,下面我介绍一下我做dem的步骤。
我的数据是一些高程点,cad格式的平面图层。步骤如下:
1、在cass里面展高程点,生成三角网,再把三角网生成等高线,保存。
2、用arcgis里的arcscene的3d analyst菜单将三角网生成tin,再把tin转换成raster格式的一个gird文件,即tin to raster。
3、在gird文件右键打开属性表,选择base heights,height选项选择obtain heights for layer from surface,这里选择tin文件,在下面z值那里可以设置高程,复合后点击ok。
4、将cad格式的文件添加进来,用第三步的方法分别把各个图层复合上去,当然这些图层复合时z值要设置为同一个值。
5、最后一步,分别设置建筑物等的高程值,方法是在图层上右键,打开属性表,选择extrusion,有个选项的默认值是0,把它改为你需要的值。ok完成。当然你可以在做完以上步骤之后把所有的图层保存一下,生成一个文件类型为ESRI ArcScene Document的文件,以后每次用arcscene打开就行了。
相册里我上传了几张成果的截图,由于我也是半灌水,以上有错的地方还望多指教。在做完以上所有工作之后,其实还可以贴图,显得更真实。我一个师兄做的毕业设计就是三维可视化,他做的成果,可以进入房间查看里面的情况,可以自由的进入,非常逼真。
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论