生成dem.doc

方法一：由DWG地形图生成DEM由DWG地形图生成DEM1.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吧 31.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站 3“工程”-“宏”-“加载工程”-加载“提取文字信息”VBA代码-找到宏中的VBA编辑器运行就可以了（也许需要添加引用EXCEL），结果保存在C盘下的EXCEL文件xyz中（他的X和Y和CAD图上是反的）。输出距离小于100米内的点个数，这样避免下面这样的情况，两个都可以。设置上限100米，这样可以分析是不是找对了.235125.1.2根据高程点数据生成DEM1.2.1生成点SHP文件ArcMap的工具添加XY数据至少可以打开ACCESS数据表和TXT文件,ACCESS文件数据表直接打开没问题，TXT文件的格式为：第一行为属性字段，如（ID,X,Y,Z），然后从第二行开始就是数据了，如（1,12.45,23.4,234.5n2，45.23，45.5，236.5.)。具体步骤为：通过ACEESS创建数据表，新建数据库-附加数据库-打开CEXCEL高程xyz数据文件-根据提示可以生成mdf数据表-运行ArcMap-“工具”-“添加XY数据”-选择刚生成的数据表-X,Y字段选择位置-OK-然后保存为shp文件，在图层列表选择刚打开的点层-数据-导出为SHP格式-到此点SHP文件创建完毕OK。1.2.2生成TIN数据和栅格有的机器上装的ArcMap中的3D分析可以用，但是我这台机器上装的不能用，不过，我发现我的机器上ArcScene是好使的，幸亏能用，呵呵。那么具步骤就是：中国3S吧 3运行ArcScene-调出3D分析工具（“工具”-“扩展”中可，也可右键工具栏打开3D分析工具，当然还可“工具”-“定制”）-3D分析工具中创建/修改TIM-从要素生成TIN弹出对话框-打开高程点SHP层-不用修改了其他的（高度源就是高程字段）-确定就OK了。转换到栅格的方法：3D分析-转换-TIN转换到栅格OK了。从网上上看到的“利用等高线生成DEM,最好还要有高程点数据,生成方法最好不要用构TIN的方法,要用ARCTOOLBOX-空间分析-内插-TOPO TO RASTER 命令,效果要好得多。此命令是专门用于DEM生成的内插方法”，做了第一个插值的，效果好像是好点，山顶没有构成TIN的图尖锐。1.2.3DEM数据的显示设置如果觉得高程差别太小或者太大，起伏不明显或玄虚太大，设置基高也许会帮助你解决视觉上的难题，按以下步骤做：右键图层名-属性-基表面高度（倒数第2个）-选中“从表面为图层获得高度（这个栅格图像起始没选中，tin是选中的）”-Z单位转换自定义后面的数字就是转换因子，数字越大，高程差别越明显。栅格起始是黑白来区分高度，tin是用起伏来表示高低，可以设置由高到低用不同的颜色表示，两种数据设置有点不一样的。具体步骤是：专业的3S站 3tin的设置：右键打开属性-符号-显示下面的“添加”添加渲染-可以选第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的下拉菜单中选择被用来进行插值的离散点数据；（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、最后一步，分别设置建筑物等的高程值，方法是在图层上右键，打开