如何根据离散点自动绘制等值线
- 格式:docx
- 大小:16.01 KB
- 文档页数:2
生成等值线图在surfer中确定等值线图的色标打开mapgis 在空间分析--DTM分析--Grd模型--离散数据网格化,打开原始数据,选择成图数据行列,然后确定(生成网格化文件默认保存在原始数据所在文件位置)。
Grd模型--平面等值线图绘制,打开生成的网格化文件,在设置选项中画勾,制图幅面--选择“原始数据范围”点击--等值层值--删除当前分段,起始Z和终止Z输入在surfer中确定的色标分段和颜色,然后点击添加分段,然后同样操作输入第二个分段,(一个分段的终止和下一个分段的起始不能重合.比如分段0-288-934-1571-2086-2400,第一个分段为0-288,第二个分段可以设为288.1-934,第三个分段为934.1-1571 ……),最后再点确定,然后要删除288.1 934.11571.1…这几层。
然后确定。
生成图文件--另存数据于—保存点,线,区文件到事先新建的文件夹,然后关掉工作区。
打开mapgis—图形处理—输入编辑—新建工程—在左边工程文件区右键添加项目打开生成的点,线,面文件.按照wp,wl,wt顺序排列在工作区,然后点击1:1或者刷新,然后成图.1.生成图框在上图中,将鼠标放在图的四个角读出四个角的坐标,确定x轴和y轴的最大值和最小值,记下来.重新打开mapgis—实用服务—投影变换—系列标准图框—鼠标生成矩形图框轮廓参数—输入x,y的起始值,例如:x的值15690200—15693200,y的值5191750—5194300 为了好看(框比图大),在这里输入的起值要比x,y的最小值小,比最大值大.比如可以输入x:15690150—15693250,y:5191700—5194350图廓内网线参数—为xy的最小值,间距dx,dy由图的xy方向的长度来定,网线类型--绘制实线坐标线坐标系—图形坐标其他的可以不管.点击确定文件—另存到事先建好的文件夹—关闭工程编辑中添加项目,把刚才生成的图框文件添加进去这不就看到框框了吗?3添加点位图打开mapgis—实用服务—投影变换—p投影转换—u用户投影转换--打开文件—原始数据设置分隔符选择xy所在行列不需要投影—生成数据—确定将生成的文件保存到事先建好的文件夹在工程编辑中添加刚才生成的点位文件添加后这不就看到点了吗?4添加标注选择需要标注的对象,标注的位置通过xy来移动(自己试验)确定后可以更改标注的大小确定这不就有标注了吗?5裁剪在工程编辑中添加一个线文件和一个区文件(注意该路径为事先建好的文件夹),如下图所示:下面造剪切框,首先选择剪切线和剪切区使他们处于编辑状态,接下来是先用用线把你想裁剪的区域画下来如下图(注意所造的区域一定要是封闭的线,因为我们接下来要用这个剪切线生成区):先编辑—输入线—输入线E—折线—确定,圈出所需要的区域R区编辑—区编辑…R—图形造区,鼠标点击勾出来的线,右键点中区域输入颜色—保存所有项—保存工程—删掉剪切线点文件和区文件—T其他—工程裁剪—保存到另一个文件夹—选择要裁剪的文件(等值线文件)载入剪切框—添加刚才生成的区文件开始裁剪—保存到新建的另外一个文件夹—关闭所有打开mapgis—输入编辑—新建工程—添加刚生成的裁剪文件—再把以前生成的点位文件,图框文件添加到这个文件夹—画出最终图谢谢!!。
•用surfer软件绘制等值线图详细说明共分以下几步完成:第一步:建立数据文件每个人需要建立两个文件便可完成所有等值线图的绘制。
(小曼与景立由于没有边界可不建立边界文件)1)绘图文件数据:一次在excell中建立完成(当然也可在surfer工作表中建立)。
格式如下:钻孔号(A栏) X坐标(B栏)Y坐标(C栏)绘图数据1(D栏)绘图数据2(E栏)。
延52 389640.0 197683.2 15.6 0.82 。
保存格式后缀为x.txt (如1.txt)2)边界数据:由于边界数据没有几个,最好在surfer 工作表中建立,这样不需数据格式转换。
步骤:打开surfer软件,点击左上脚空白处(第二行第一个位置键),选worksheet,则可建立边界数据文件。
格式如下:X坐标(A栏)Y坐标(B栏)5 0389640.0 1976830.2392300.4 1985301.4392100.8 1988023.4392235.6 1982003.1389640.0 1976830.2第一行5代表边界坐标的个数,注意后一行坐标数据必须重复一下第一行坐标数据,这样才能形成一个闭合边界,0代表边界内的等高线要绘制。
存盘格式为x.bln (如1.bln)第二步:数据网格化建立好两个数据文件后,打开surfer软件。
在Grid(网格)中——点data(数据)——把绘图数据文件x.txt打开——在图标左边对X,Y,Z坐标进行选项,由于先画等值线,因此,必须选绘图数据文件上B、C、D栏(以后做其它等值线可将D栏该选为E栏或F栏,等等),即X选columnB,Y选columnC,Z选columnD(必须切记)——其它项目默认不变,点击ok——若出现是否替代存在的文件,再点击yes(若出现将改动保存到report1?,点击N o,surfer8.0有此要求,surfer7.0也许没有)——自动生成x.grid文件。
第三步:绘等值线图在map(地图)中点contour map(等值线图)再点new contour map(新建等值线图),再点网格化文件x.grid——打开——等值线图形便绘出。
1、输入数据
打开surfer12,新建工作表,要求输入三列数据,分别坐标点(X、Y)和每个坐标点的数值(顺序可以任意)。
然后,打开“数据-指定XYZ各列”,设置X、Y、Z对应的各列标号,并保存文件。
2、生成网格
点击“网格-数据”,并选择第一步生成的文件,在弹出的对话框中设置网格化参数,网格化算法选择克里金插值法,再点击“确定”进行网格化操作并保存文件。
3、生成等高线图
点击工具栏左侧图标建立新场景。
选择“图形--新建—等值线图”选项,并打开第二步生成的文件,即可自动建立等值线图。
4、调整显示选项。
为了使得生成的图像更容易观察,需要适当调整。
先按照下图进行勾选,并在”顶边轴”的属性管理中选中”标注-显示“,以显示顶轴刻度。
在等值线图的属性管理器中,更改”层次-等值线着色-填充等值线“,并更改等值线颜色。
可选中“比色刻度尺“选项,以显示比色刻度尺,便于观察。
最终生成的等值线图如下图所示:。
Mapgis作等值线图方法用Mapgis作等值线图,方法很简单,用Mapgis中的空间分析、DTM分析子块中的Grd模型进行:1、数据准备 Mapgis6.7能识别的数据的格式是DAT格式,不能用excel 或者TXT格式。
这后两种格式别的地方可以用,唯此处用DAT格式。
将检查无误码的包含了平面位置数据和分析数据的TXT文件,通过别的软件(如SURFER等)转换即可。
在SURFER中,直接在文件菜单中打开TXT格式数据,再另存DAT格式。
2、数据网格化打开“DTM分析“子模块,出现“三角剖分显示窗口1“,单击“Grd模型”工具中的“离散数据网格化”,选择要处理的DAT格式数据文件。
出现“离散数据网格化“对话框。
在该对话框中,要有几个选项:1)在“数据列D“中,将横坐标、纵坐标和分析数据,分别对应着X、Y、和Z。
2)在“网格参数设置“中,将你的图面处理范围设置好。
即上、下、左、右侧坐标输入进去。
3)在“网格化方法“中,选进行处理的方法,一般选用”泛克里格“法比较好,直观一些。
4)网格化后的数据用“GRD“作为后缀。
默认存于你存放TXT数据和位置。
文件名是可自拟的。
默认的文件名“TmpGrid“。
5)点击确定后,网格化数据处理完毕。
对话框消失。
注意:此时的GRD格式数据,SURFER系统是可以使用的。
但SURFER系统生成的GRD文件,Mapgis不能识别。
怪事情!!3、作等值线在“三角剖分显示窗口1“下,单击“Grd模型”工具中的“平面等值线绘制“,出现”设置等值线参数“对话框,对等值线是否套区、是否绘制色阶、是否保留边框、是否光滑等进行选择。
对等值线层是否显示,在”注记参数“正下方勾选”YES或NO。
单击“注记参数”还可以对图面上的标注文字进行大小、字体、颜色进行设置。
线参数和区参数亦可分别设置。
视情况,还可以增加阶层可删除阶层。
注意“制图幅面”框下的选择,一般选择“原始数据范围”。
设置好、确定后,系统会提示是否保留文件。
从事煤炭资源勘查工作中,需要做很多等值线图,周围同事多是用传统的三角网插值方法,一个钻孔一个钻孔地量取之间的距离,然后在Excel表格中计算分点距离,再在CAD 中描点连线……虽然这是一种严密的插值法,但是感觉过程复杂,与传统手工作图无异,需耗费很大精力。
同学推荐用Surfer软件,我用了一天时间,基本掌握了这种作图方法。
这里,我要感谢水文07-1班的邓同学以及其他同学朋友给我提供的帮助。
Surfer做等值线图主要用到两种数据,一种是钻孔数据,包括钻孔坐标和需要描述的目标数据(如煤厚、灰分、硫分等),另一种是边界数据,主要是矿区边界,即拐点数据。
这两种数据分别对应着Surfer里的dat文件(也可以是其它可识别格式,如txt)和bln文件。
第一步:制作数据文件1、制作钻孔数据文件这一步可以在Excel、WPS、记事本中完成,也可以在Suerfer工作表中建立。
其格式如下图:数据录入完成后,保存为txt格式。
注意:我用的是WPS,另存为的时候须选择制表符分割的txt文件。
2、制作bln文件(边界数据)一般情况下,边界数据都是比较少的,可用Surfer工作表直接建立,保存为bln文件即可。
其格式如下图:第一行31代表边界拐点坐标的个数,后面的0值表示后面的等值线绘制在边界内部,如果写成1,表示等值线绘制在边界外部。
后面的数值为钻孔坐标,第一栏为横轴(Y值)。
保存为bln格式。
注意:我们工作区域只有30个拐点,写31的原因是最后一个钻孔坐标与第一个相同,形成一个封闭区域。
Ps:有时候拐点数较多,又没有拐点坐标文件,需要我们直接从图上量取,比较方便的方法是用工具展点。
比如,我用的CAD作图,作图比例为1:1000(保证了展点出来的坐标不变),矿区边界是多段线圈成,拐点较多,我就用到了Dxf2xyz2.0这个工具来展点。
把矿区边界带基点复制到空白文档中,并校正。
保存为dxf格式,据说dwg格式也行。
然后打开Dxf2xyz2.0,如下图:点击Output Format:选择输出格式为Blank Outside(0),即可输出为标准的bln文件。
1.将mopinfo中的散点(点)和边界(是面而非线)图层转出,生成.mif文件;(注意:散点图层必须使用相应的excel2003文件生成,2007版本生成的.mf文件会导致topmap无法识别。
)2.打开topmap,工具——MIF文件导入工具,将之前生成的散点和边界的mif文件导入,首先选择topmap图层保存目录,然后添加要导入的图层文件,最后点击开始导入,显示导入成功后退出;3.新建地图,地图——新建;添加刚才导入的散点和边界图层,图层——添加图层,选择数据源目录,在地图数据一栏中选中要添加的两个图层文件,单击确定;4.分析——等值分析;参数选择:等值线;计算方法,多面函数法;散点图层,对应的点图层文件;数值字段,勾绘等值线依据的数据字段名;边界图层:对应的边界图层文件;目标图层:生成等值线的文件名;数值列表:选择起始值、步进值和终止值后单击步进分布,然后双击选择首数值和末数值对应的颜色,单击颜色渐进,然后确定即可;(注意确定后是否有等值线生成)5. 图层——保存当前图层(即等值线图层);到此即完成了等值线的绘制,下面需要把等值线导入mapinfo 中进行编辑,有两种导入方法:1. 直接打开法。
在mapinfo 中打开刚才保存的等值线图层,打开,选择*.shp 格式,找到等值线图层,单击打开,填写保存的文件名称,单击保存;出现shap 文件信息,单击确定;现在完成了文件的导入,但导入的文件不可编辑,需要将导入的等值线文件另存副本,然后打开刚才另存的文件,分别选中每个值的线,分解方可编辑。
(右键——编辑对象——分解)2.利用导入工具导入:填写源文件格式(Format:ESRI Shape),选择文件;选择投影体系为Longitude/Latitude (WGS84) [EPSG:4326],注意投影前面的Use projection setting in sou 不勾选,否则导入后不能在原地图窗口显示;填写目标文件的格式(MapInfo TAB),选择目标文件存放位置,单击OK即可。
DTM空间分析(绘制等值线图)MAPGIS--DTM空间分析1、GRD模型电子表格(.XLS)→另存为CSV(逗号分隔)格式文件/逗号分隔的文本文件.TXT(无表头)→用于DTM空间分析→GRD模型→离散数据网格化(注意:输出网格数据文件“*.Grd”换名保存)→平面等值线图绘制→存数据于---点、线、区数据文件。
(1)数据准备:形成CSV(逗号分隔)(*.CSV)格式文件。
或:生成逗号分隔的文本文件(*.TXT)(无表头)。
(2)打开MAPGIS——空间分析——DTM分析,选择“GRD 模型”—“离散数据网格化…”。
按下图步骤分别设定好参数后,确定。
(3)选择“GRD模型”—“平面等值线图绘制”选择网格化数据文件(*.GRD)打开后,按下图步骤进行等值线图绘制参数设置。
(4)各参数设置好后,按确定,系统自动生成等值线图。
在确认无误后,将其点、线、区数据文件分别另名保存下来。
(5)套标准图框新建一个工程文件,工程地图参数“从文件导入…”标准图框的投影参数。
先将标准图框添加进来,然后把等值线图也添加进来,有时系统提示地图参数不配,不理会直接按确定添加进来,1:1刷新显示。
另名保存工程。
2、TIN模型(1)电子表格(.XLS)→另存为文本文件(制表符分隔)→添加表头“NOTGRID”→将“空格”替换成“逗号”,成为逗号分隔的文本文件→改后缀名保存为“.DET”格式文件(保存类型:文本文档;编码:ANSI)→用于DTM空间分析→TIN模型。
(2)选择“文件”→“打开/装入三角剖分文件…”,打开“.DET”格式数据文件;选择“Tin模型”→“快速生成三角剖分网”/{“生成高程初始三角剖分”(三角网构造Tin判别系数-0.001);选择“Tin 模型”→“优化高程初始三角剖分”};选择“Tin模型”→“整理三角剖分网…”,设置相关参数。
(3)选择“文件”→“另存剖分于…”另名保存三角剖分文件(*.TIN);选择“Tin模型”→“追踪剖分等值线…”,进行等值线图绘制参数设置。
MAPGIS制作等值线图【等值线图】1.做好Excel原始数据表:第一行为属性行(Y X Au);第二行以下为野外采集数据及Au含量化验数据。
2.将Excel原始数据(文件菜单)另存为:文本文档(制表符分隔).txt。
3.利用MAPGIS投影变换模块生成点数据文件 (﹡WT)注意参数设置!①指定数据起始位置→选第2行。
②设置用户文件选项→选按指定分隔符→设置分隔符→分隔符号:选Tab 键和空格。
③属性名称所在行→选第一行(Y X Au)。
④数据类型→选5(双精度型)。
4.MAPGIS主菜单→空间分析模块→DTM分析→文件菜单→打开数据文件→点数据文件(﹡WT)→处理点线菜单→点数据高程提取(高程属性字段选Au)→Grd模型菜单→离散数据网格化,见对话框。
设置对应参数(主要是网格化方法的选择:距离幂函数反比加权网格化,Kring泛克立格法网格化),选择网格化方法和输出网格文件(.grd)路径→Grd模型菜单→平面等值线绘制→选择上一步生成的网格文件(默认名为:TmpGrid.GRD)→设置等值线参数(可双击等值层值,线参数,进行修改)→确定即可。
注意:制图幅面一般选择原始数据范围,以方便生成的等值线与原始地图复合叠加。
或:① MAPGIS主菜单→空间分析模块→DTM分析→Grd模型→离散数据网格化→打开.TXT文件,选择数据列(X、Y、Z)—网格化方法(按不同数据选择不同网格化)—输出网格文件名F(存储数据的位置)—确定。
②文件→打开三角部分文件→选择上面保存的.grd文件打开。
③ grd模型→绘制平面等值线图(设置相关项)→确定即可!④文件→存数据于—点数据文件(.wt)、线数据文件(.wl)、区数据文件(.wp)。
5、菜单:Grd模型→平面等值线绘制。
选择上一步生成的网格文件(默认名为:TmpGrid.GRD)弹出设置等值线参数对话框,可双击等值层值,线参数,进行修改。
注记参数设置相应参数。
注意:1)制图幅面一般选择原始数据范围,以方便生成的等值线与原始地图复合叠加。
如何用Section作平面等值线图
在我们平时的工作中,常常会遇到需作等高线或等值线图的情况,如根据煤层底板作煤层底板等高线图,根据GPS点高层数据作地形等高线图等情况。
以前我也为这事烦恼,有了Section这件事就变得简单了,下面我就根据自己的实际操作步骤写出来,和大家一起分享。
1、平时我们可以将各地质点的X、Y、Z作好一个Excel表,等
要作图时将其另存为文本文件。
2、打开Section—1辅助工具—绘制等值线---高程数据网格化----
找到要作等值线的那个文本文件,打开。
此时弹出如下界面:
3、在数据列中,选择好X、Y、Z所对应的列和网格化方法,确
定。
4、Section—1辅助工具—绘制等值线---高程数据网格化—平面
等值线图,找到刚才网格化得到的文件,打开。
5、设置等值线参数:
⑴设置是否等值线套区,保留边界;
⑵等值线是否光滑处理;
⑶图幅范围一般选择原始数据范围;
⑷等值层值可以根据需要设置,只要双击就可以设置;设置线参数;
⑸点击注记参数,可以设置标注字体、格式、起始层、频度等;
⑹单击区中和颜色可以修改不同等值范围的颜色;
⑺全部设置完成后,点击确定。
6、然后点击1:1,即可看到等值线了。
7、将这个文件的点、线、文件保存,然后添加到相应的图中,就形成了底板等高线或等值线图了。
等值线图的自动勾绘DFSurface模块是一个具有带断层网格化功能的程序,它处理大数据量及断层(包括逆掩断层)的能力强,并且平滑效果合理,是进行速度研究及自动形成各类等值线的理想工具。
3.13.2菜单介绍主要介绍“建立”菜单中的三个命令:设置、分逆掩断层数据、获取轮廓曲线。
1、设置设置软件的加密方式及勾绘结果的显示方式等,如图3-1-78所示。
加密有加密在唯一一台电脑上、加密在唯一一只加密锁上两种方式。
图3-1-78 设置过程参数定义2、分逆掩断层数据输入文件中包括散点数据及断层数据,指定断层所在的层名,点击“确定”按钮后自动划分逆掩断层数据,散点数据及断层数据分别输出在两个文件中。
如图3-1-79所示:图3-1-79 分离带有逆掩断层的坐标文件和断层文件3、获取轮廓曲线自动获取指定数据的轮廓曲线,第一步指定输入文件,第二步指定输出的轮廓曲线文件及轮廓曲线采样点步长,然后点击“确定”按钮结束。
勾绘完成后,可以利用此曲线对输出结果进行剪切。
如图3-1-80、图3-1-81所示:图3-1-80图3-1-813.3操作过程启动DFSurface程序,主界面见图3-1-82所示:图3-1-82 DFSurface界面一、点击“如果您是个新手…”按钮后见图3-1-83所示,点击“xyz 文件”按钮指定散点文件,格式列表中列出了您如何获得数据,缺省为双狐散点格式,其中还有自定义(固定符号、固定宽度)、从数据库中获得数据三种格式。
图3-1-83 输入文件双狐缺省格式中X 、Y 表示坐标;Z 代表对应点的值,可以是高程、重力值、T0值、元素含量、地层埋深、某一种气象信息等数据等,S 代表对(X ,Y ,Z )这个空间数据点的可信度,S 越小表示可信度越高,S 越大表示可信度越低,当S 为0时表示此点完全可信,对S 值的控制可以实现对局部数据参与运算能力的控制;F 表示断层号,一般在对逆掩断层勾绘时用到此参数。
这种数据为散点数据,它是平面上分布的一系列坐标数值。
使用测绘技术绘制地下水位等值线图的方法随着城市化的快速发展和人口的不断增加,对地下水资源的保护和合理利用成为了当今社会亟待解决的问题之一。
而绘制地下水位等值线图是了解地下水分布和水文特征的重要途径之一。
本文将介绍使用测绘技术绘制地下水位等值线图的方法。
一、地下水位等值线图的意义地下水位等值线图是通过对地下水位进行测量和分析,绘制而成的地下水位分布图。
它可以直观地展示地下水位的高低、斜降、变化趋势以及地下水流向等信息,非常有助于我们了解地下水资源的现状和变化情况。
通过定期绘制地下水位等值线图,我们可以及时发现地下水资源的利用状况,为科学合理地利用和保护地下水提供参考。
二、测绘技术在绘制地下水位等值线图中的应用1. 高精度水位计的使用高精度水位计是绘制地下水位等值线图的重要工具之一。
它能够准确测量地下水位,提供可靠的数据基础。
在使用高精度水位计时,我们需要选择合适的测点,尽量避免人为和环境因素对测量结果的影响。
同时,还需要使用多个测点进行测量,以获取更加准确的数据,确保地下水位等值线图的可靠性和准确度。
2. 地下水位数据的处理和分析在获得地下水位数据之后,我们需要对其进行处理和分析。
首先,对数据进行筛选和整理,去除异常值和干扰因素。
然后,使用插值算法将离散的数据点插值为连续的等值线,以形成地下水位等值线图。
在插值时,通常会选择合适的插值方法,如逆距离加权法和克里金插值法等,以充分利用已有数据的信息。
3. GIS软件的应用在实际绘制地下水位等值线图的过程中,我们可以借助地理信息系统(GIS)软件来辅助完成。
GIS软件具有强大的数据处理和地图绘制功能,能够实现数据的可视化展示和分析。
通过将测量数据导入到GIS软件中,我们可以进行更加灵活和高效的数据处理和图像绘制,从而得到更加准确和美观的地下水位等值线图。
三、地下水位等值线图的应用和展望地下水位等值线图的应用范围非常广泛。
它可以被用于水资源管理、环境保护、灌溉农业等领域。
地基微波辐射计资料的等值线自动绘制方法地基微波辐射计是一种用于测量太阳辐射和地表反射辐射的仪器,可以获取到地表各种要素的辐射数据,从而可以用于地表反照率等参数的计算和分析。
然而,这些数据的处理和分析需要使用到等值线自动绘制方法,以便能够更好地展现出数据的空间特征和变化规律。
目前,地基微波辐射计所采集到的数据以及其它多源数据的处理与分析越来越受到关注,因为这种方法十分适用于对地表要素进行反演和监测。
等值线自动绘制方法是一种常用的数据处理方法之一,能够自动根据给定的数据绘制出等值线图,并可以对数据进行统计和分析。
本文将介绍一种基于MATLAB的地基微波辐射计资料的等值线自动绘制方法。
1. 数据预处理在进行等值线自动绘制前,需要对地基微波辐射计采集的数据进行预处理,包括数据清洗、数据插值、数据平滑和数据标准化等步骤。
首先,需要将采集到的数据进行清洗,即去除缺失值和异常值。
然后,采用插值方法对缺失点进行补充,以便能够建立起完整的数据集。
其次,对数据进行平滑处理,以减少数据中的噪声对等值线绘制的影响。
最后,对数据进行标准化,以便能够更好地展示数据的空间特征和变化规律。
2. 等值线自动绘制原理等值线是将多个等值点(即相同值的点)用线段连接而成的线,用于展示数据的空间分布和变化规律。
等值线自动绘制方法是将给定的数据表面用等值线表示出来的过程。
等值线绘制的过程是将数据表面分割成等间距平面,根据一定的间隔依次搜索每个点的四个方向,如果发现下一个点的值与当前值相同,那么就将这两个点之间的线段划入等值线集;否则跳过这个点,继续搜索下一个点。
3. 等值线自动绘制步骤(1)制作等值线网格将地基微波辐射计采集到的数据做成网格形式,采用Meshgrid函数将数据集转化成网格形式,以便进行等值线的计算和绘制。
(2)计算等值线集采用Contour函数,用等值线绘制的方法计算等值线集。
初步绘制时,可以通过选定不同等值线的数目和间隔来进行初步绘制,以便确认等值线的走向和线的数目是否适宜。
如何进行地形高程数据的等值线绘制地形高程数据的等值线绘制是制作地图中非常重要的一环。
通过等值线图,人们可以直观地了解到一个地区的地势起伏情况,以及地形特征。
本文将介绍如何进行地形高程数据的等值线绘制。
一、收集高程数据进行等值线绘制的第一步是收集高程数据。
高程数据是指地表上各点的海拔高度或地面起伏情况。
通常可通过航空激光雷达(LiDAR)或卫星覆盖数据来获取高程数据。
此外,还可以在实地进行测量,并通过GPS获取各点的具体坐标和高程信息。
二、数据预处理在进行等值线绘制之前,需要对收集到的高程数据进行预处理,以消除错误和噪声。
首先,需要检查数据的准确性,并删除可能存在的异常值。
然后,通过平滑等方法,对数据进行滤波处理,以去除不规则的噪声。
同时,还可以采用插值算法,填充数据中的缺失值,以保证绘制等值线的准确性和连续性。
三、等值线的选取和绘制等值线绘制是根据高程数据中的等值线进行描绘的。
等值线是连接具有相同高程值的点的曲线,相邻等值线之间的高程差是一个固定值,通常被称为等值线间隔。
选择合适的等值线间隔是绘制等值线的重要步骤。
如果等值线间隔设置过大,会导致绘制出来的等值线数量较少,地形特征不够准确;如果等值线间隔设置过小,会导致绘制出来的等值线数量过多,不易观察。
在选择等值线间隔之后,可以使用计算机软件或绘图工具进行等值线绘制。
一般来说,绘制等值线的流程是:首先选取一个起始点,然后根据等值线间隔,在该起始点的周围找到一条等高线,并标注上相应的高程值。
接着,从该等高线的一个未标注点出发,继续找到下一条等高线,直到将整个地图区域都绘制上等值线。
四、等值线的样式和标注在绘制等值线时,除了要选取适当的等值线间隔外,还需要考虑等值线的样式和标注。
一般来说,等值线可以使用不同的颜色或线型来表示不同的高程区间。
低海拔的区域可以使用浅色表示,高海拔的区域可以使用深色表示。
此外,还可以根据等值线的密度,调整线型的宽窄,以突出地区的地形特征。
基于不规则三角网(TIN)的追踪等值线算法基于TIN绘制等高线直接利用原始观测数据,避免了DTM内插的精度损失,因而等高线精度较高;对高程注记点附近的较短封闭等高线也能绘制;绘制的等高线分布在采样区域内而并不要求采样区域有规则四边形边界。
而同一高程的等高线只穿过一个三角形最多一次,因而程序设计也较简单。
但是,由于TIN的存贮结构不同,等高线的具体跟踪算法跟踪也有所不同[1],下面首先介绍TIN模型的存贮结构。
1不规则三角网TIN存贮结构:(1) 三角形的存贮结构struct triangle{long ID;//三角形的标识;BOOL IsBorder;//标识是否边界三角形;BOOL Flag;//标识三角形是否使用了;Triangle* la,lb,lc;//指向三角形三个边指针;line* ta,tb,tc;//指向la,lb,lc三边相邻三角形指针;}(1) 三角形边的存贮结构struct line{long ID;BOOL IsBorder;//标识是否边界;BOOL Flag;//标识三角形边是否使用了;Point* pb,pe;//指向三角形边的两个顶点;triangle * tr,tl;//指向与此三边相邻左(left)右(right)三角形;}(1) 顶点存贮结构struct point{long ID;//边的标识;BOOL Flag;//标识顶点是否使用了;float x,y,z;//顶点的三维坐标;}2 在建立TIN模型的基础上,自动生成等值线的步骤如下:(1) 设定等值特征值Z,清除所有三角形标志flag为false,逐个搜索三角形,已经处理过则把标志设为true,直到高程值改变。
(2) 先搜索边界三角形的边界边,查找是否有等值线经过,然后搜索内部三角形;判断是有等值线通过:根据三角形三个顶点的特征值分别进行以下分析:a.对于给定等值线特征值,三个值全部大于或小于当前特征值Z,则当前三角形内不存在等值点;(4)对下一个等值点的追踪,与步(3)相同对另两边进行判断和插值处理,得到等值线在此三角形的出点,把此三角形flag设置为true;(5)把上步得到的上一三角形的出点作为共用此边的邻三角形的入点,重复(4),直到得到等值点在没有其他邻接三角形的边上,或者得到下一等值点所在边已经处理过的三角形上(形成封闭等值线)。
Mapgis画等值线的方法1.在Excel中按照x、y、z的格式输入数据,x、y表示平面直角坐标数据,z代表该点数值(高程、某种元素含量等数据),不要标题行。
注意Y坐标带不带18或其他的分度带的带号。
另外MAPGIS中X是指18等分度带开头的坐标,Y是指255开头的坐标,注意别弄反。
2.将Excel表格另存为“文本文件*.txt”,注意最后一空白行要删除,光标删到最后一个字符后面。
3.打开mapgis,“空间分析”模块-DTM;在菜单栏的“Grd模型”中找到“离散数据网格化”窗口,打开由Excel表格另存的.txt文件。
在“网格化方法”下拉列表中选择“Kring泛克立格法网格化”后,设置最大最小值,再点击“文件换名”按钮保存你的.Grd格式的文件(设置保存到的文件路径)-确定。
网格间距太小则会内存不足,最好在200以下,300-400也都可以。
4. DTM对话框下-在菜单栏的“Grd模型”中找到“平面等值线图绘制”项,将已保存好的.Grd 格式的文件打开,进行相应的设置(按使用者对等值线间隔等参数的设置)后即可看到绘制成功的平面等值线图。
分别在设置、光滑(中程度)前全部打∨1、不选″等值线套区″→不选″绘制色阶″→不选″保留边界″→选″等值线光滑处理″→在″光滑度″选项中选择″高程度″;2、点击″等值层值″→对等值线层分段参数中的″起始Z″、″终止Z″、″步长增量1或2(沙漠区),20(沙漠外围)″进行设置→点击″更新当前分段″→确定;制图幅面(原始数值范围)注记参数-频度改为1,-注记间最小允许距离为10000,其余不改。
线参数改为光滑曲线注记参数一般都显示-选yes。
确定,处理完成后保存。
保留边界-是指生成的等值线图文件有个黑框。
等值线套区-是指各个等值线之间为一个区。
绘制色阶-是指给各等值线间的区赋一个颜色。
注记参数中频度是指多少根线标一根的值。
注记间最小允许距离-是指标记的等值线值的间距,越大每根线标记的越稀,数字越小标记的越密。
如何根据离散点自动绘制等值线(等高线)之三角形法
2001-04-06 16:35 6683人阅读评论(1) 收藏举报算法网格图形c
自动绘制等值线的方法从技术方向上看可以分为两大类,插值和曲线拟合.其中曲线拟合总的来说效果不如插值算法经典和应用广泛,效果也较逊色.这里着重介绍插值算法.
其中插值算法中,按照方式不同分为离散点客观化和三角网方式.两者区别在于三角网计算主要在生成三角网过程,省去了插值到格点的过程.而客观分析过程则是将离散点分析到格点后再内插到细网格,然后大多利用追踪法生成等值线,也有在这里再使用曲线拟合.
1、三角形算法
a、首先生成delaunay三角形,这一点在我的帖子"delaunay triangulation之丰衣足食“内有源程序,大家可以参考。
b、随后需要在三角形的边上插补等值点。
要确定某个三角形的边上是否有等值点,需要进行判断和处理。
注意:如果某原始数据点和等值线值相同,将该点改变一个微量。
如果一个三角形三顶点的值相同则各边无等值点。
如果一个三角形的任意边两端点(A、B〕的Z值(Za、Zb)满足满足(Zd-Za)*(Zd-Zb)<0,其中Zd代表等值线的值,则该边必有等值点,其平面位置是
Xd=Xa+(Xb-Xa)*(Zd-Za)/(Zd-Za) , Yd=Ya+(Yb-Ya)*(Zd-Za)/(Zb-Za)。
每个三角形上不可能三边都有同值的等值点,另一边上必定有同值的等值点。
c、等值点的追踪。
为了能将内插的等值点顺序追踪排列,绘出等值线,还必须找出相互重叠的环形网内所计算的等值点间的平面位置关系。
因每个环形网都是由多个三角形组成的,我们先简单分析一下单个三角形中存在等值点的情况。
由于不必考虑等值线穿过端点,如果一个三角形的边上存在等值点的话,只可能在某两条边上存在等值点,而不可能三条边上同时都有。
也就是说,只要三角形一边上存在等值点,则其余的两条边中必有一边存在等值点。
根据上面的约定,我们再研究等值线穿过任一环形网中两条及两条以上相邻的径边时,可能出现的几种情形:
① 等值线不通过环形网的界边。
在这个环形网中,必然所有的径边上都存在等值点,如果这个环形网由非边界点联结而成,内插的等值点就可顺序连接为一条闭合曲线(图5(a));若此环形网由边界点联成,那么这些等值点则连成一条开口曲线(图5(b))。
② 等值线通过环形网的界边,且次数不超过两次。
这是最常见的一种情形,如图5(c)、5(d)所示。
相邻径边上内插的等值点顺序排列,点数至少为两个,其起点为环形网的入口点,终点为环形网的出口点。
③ 等值线四次通过环形网的界边。
环形网中内插的等值点分为两部分顺序排列,每个部分都包括一个入口点和一个出口点,这个网所在的位置应该是地形的鞍部
由于离散点环形联网是沿同一方向(逆时针)进行的,环形网中相邻径边上内插的等值点所排列的顺序,也相应围绕中心点位逆时针旋转。
从对图形的分析中,我们还注意到,如果等值点不是位于边界上的话,那么一个环形网的入口点,必然是另一个环形网的入口点;一个环形网的出口点,也必然是另一个环形网的出口点;而内插入口点(或出口点)的径边的两个端点,就是联结这两个环形网的中心点。
利用这个原理,我们就可以成功地设计出等值线的追踪方案,且在追踪等值线时,只需将各环形网中内插的等值点进行单向比较,即入口点对入口点比较,出口点对出口点比较。
联网结束后,凡是没能联成闭合环形网的离散点,即为绘图区域的边界点,而在两个边界点连接的边上内插的等值点,就是开曲线的线头。
找到线头后,根据上述原理,就可顺序追踪出各条开曲线的全部等值点。
对于闭曲线来说,任一环形网中内插等值点中的起点都可作为线头,按上述方法追踪,直至又回到该点为止。
等值线追踪完成后,即可进行曲线的光滑输出。
目前,常用的是样条插值。
这方面的例子很多,不再详述。
另一种离散点分析等值线的方法是将离散的资料分析到网格点上然后进行格点的曲线拟合或者等值线追踪。
常用的插值方案很多,一般来说比较精确的是逐步订正法和最优插值。
最优插值着重于各插值点之间的物理或者其他的关系,在所有的插值中最为准确,并且物理意义明确。
但是有的时候并不需要这么复杂,或者你所提供的资料仅仅有一项,没有相互关联的几组资料。
下面介绍逐步订正法。
主要包括几步1、选取预备场2、进行客观分析3、输出资料的平滑
其中心思想顾名思义,就是用实际资料与预备场之差去改变和订正预备场或初值场,得到一个新场,在由新场去求出和实际值之差,去订正上一次的场,直到订正场逼近实际资料为止。
在资料离待分析的格点太远时,初值将不改变。
对于在网格点以内的区域可以使用二次全区间插值即可,大家应该比较熟悉,不再详述。
对于在所选取的网格区域以外的点,则使用cressman插值。
cressman 插值权重系数为Wij=(i**2-d**2)/(Ri**2+d**2) 当d<=Ri 当d>=Ri 时Wij=0 其中Ri为影响半径,d为离散点到格点的2位或三维距离。
经过每次比较插值点和实际资料的误差,生成订正值,直到订正值和实际资料的误差达到你所满意的精度时,可以进行下一步的计算。
一般订正10足够。
然后对得到的格点资料进行细网格嵌套。
等值线的追踪和三角形的稍有不同,追踪等值点和三角形的基本相同,主要在于追踪等值线的位置。
每当进入一个矩形时,也就是矩形一边有等值点,那么较三角形复杂,有几种情况,需要考虑,也就是等值线的出口有3条需要注意高程点的记录和跟踪,同时对起始点的判断也是需要注意。