surfer制作冲淤变化图

1、冲淤变化图制图要点注意：在研究一个新的区域的冲淤变化时，最好首先去确定冲淤变化的研究范围，根据这个范围取点、白化等，这很重要1、原始数据的处理：水深点、背景底图目的是得到可以在surfer中处理的数据，主要包括水深数据和工程图（范围、建筑、标记，用作背景）原始数据：一般为各种cad工程图（如果有处理过的水深点数据就可以直接用），其中用到的主要数据是各时间段某个区域的水深数据处理方法：（1）检查cad图中水深点数据的属性（特性），有时有些数据点的属性中没有水深值，则必须提取出图中标记在点上的那个值，此时可用韩志远给的插件提取，得到TXT文件，如图1所示。（2）如果是用作背景的工程范围图，则可直接

2、另存为dxf格式的数据，然后在surfer中作为背景图直接使用。文件的编鼠E）悟式（0直看M新助的30,54309007.050499627.3600.00030.04309012,100499034.1500.00030.04309004.7004兆65,0800.00029.14309029.970499651.0100.0002904309033.790499644.5700.00029.34309037.770499037.3000.00031-54300932.670500345.1700.00。31.4430831,均0500328.9700.00091.9430831.09050

3、0313.2200.00031.54308930.540500297.4000.00031.5430B928.370500257.3100.00032.14306921.570500200.4100.00032.043069ZL430500184.8000.00032,24308921.%。500L75.4100.00032,14306923.划500160.3600.00032.24308923.960500152.3500.00032.34308924.160500144.2800.000图1用插件从dwg数据中提取的水深点数据2、用水深数据生成网格源数据：水深点数据，dat、txt等格式

4、处理方法：（1）打开surfer,点击Grid->Data,弹出后选择水深数据文件。（2）根据水深文件的具体情况及制图的需求对数据进行网格化，保存为.grd网格数据，如图2、图3所示。（3）对各个时期各个工程区域的水深数据都进行网格化注意：sufer计算冲淤变化时要求参与计算的两个grd文件有相同大小和数量的网格，即他们的xy坐标范围和节点数量（node）都相同，根据此设置图3中的参数。（可参考网格化文件夹中的参数截图）图2选择水深文件后注意文件内点数据的排序图3注意XYZ对应的数据列，根据点的密度选择合适的网格搜索半径，根据需要设置网格大小和XY的范围3、白化处理blankSurfer

5、对水深数据插值计算生成网格数据时是覆盖整个矩形范围的，无论每个网格区域内是否有数据点都会在该网格中进行插值计算。因此计算冲淤变化之前先要确定是否有需要白化的区域，如果有就要先制作bln白化文件源数据：网格化后的grd数据处理方法：打开网格化后的grd数据（以contourmap的形式），可看到有效工程范围内有等深线，右键点击地图，弹出菜单中选择Digitize,此时可在地图上点选一个多边形，将无效区域围起来（可能需要多个多边形才能将无效区域包围，不过尽可能用一个多边形包围），Digitized对话框会记录所有点的坐标，如图4所示。多边形绘制玩后，在Digitized对话框中选择File->

6、;Saveas,将其彳存为bln文件，完成后还要再用记事本打开改文件，把第一个坐标行复制一份到最后一行，使多边形闭合，第一行第一个数字表示点数，加1,第二个数字为标志位，选择0,保存，如图5所示。如果有多个多边形，则对相应的bln文件修改后可将其内容合并到一个bln文件中。最后用Grid->Blank对grd数据进行白化处理口DigitizedCoordinates-digitiz.gileEditQpti6rl:图4用数字化工具选择需要白化的多边形区域13,0498903.59210673,4989U,97522016,48843.23412452,496817.29688101149

7、130.51490524,499619.18241507,4199042.97738882,499649.86140617,499687.08514506,499308.10666059,499123.55052207,499101.49664456,498903.59210078,4309324.10599874309262.00692084309246.33706014309236.47085144309118-05&71314303921,.91290574308323-07303014309926.555S2744309020,37864414309179.5950605430

8、9248.07BL5574309252.72107744309324.1059987,图5bln文件格式4、用网格数据计算冲淤变化图源数据：白化后的grd数据（如不需要白化，则用网格化后的grd数据），一般分为不同时期不同区域处理方法：选才GGrid->Math,弹出的对话框中先后选择同一区域不同时间的grd数据，用时间小的数据减去时间大的数据（如13年-14年），保存计算结果，如图6所示。TnputGridY*ri.日工arJ?KutJ.AddGrids,.D:caobeilian.xy1303301406leftJx52013(13L-hBlanl<outputD-ii1013

9、03-21)1406leftKS2D1.EElmktutputS.»KStheformf(即瓦)where角aretintvariablesinthelistatravt.A-BGridJnfb,.OutoLitGric：Le:图6计算冲淤变化5、调整冲淤变化图在冲淤变化图（contourmap）的属性管理器中设置其填充、等深线级别和颜色卡等等。颜色卡可保存为clr文件供其他图使用。调整颜色条colorscale属性，使其位置、大小合适。用划线和文字的方式添加图例，指示所计算的冲淤图的2个时间点。6、添加2个时间的等深线以contourmap的形式分别添加两个白化后的网格数据，设置其

10、属性为不填充，并设置显示线的颜色，标记，间隔等等。如图7所示。GeneralLevelnmetfiodAdvarcedCortourlevel&EdtItIs.FilledComouriIHIIcontoursColorscaleCsntcurItvHs:Displaythecontourlevsedito图7等深线的属性设置，不填充，没有colorscale,线的颜色和level的级别7、添加背景图如图8所示。立口fontaurs4Ut_Lthl_bln,grd3Mmp刁iRightAxis7hmiLe+tAxis3TopAxis习nBottomAxis引回Po5t-uut-fi.d

11、at引DBase-Drawing2rdxfContours-KS-序|LevelsLayermngeleral?lmethoditoLrlevelsmdContours:ontoursDrscale日渡layerfromani*CutCtrl+X口CopvCtrl-CXDqlfteDfkURename口bjectF22oomSelectedCtrHLOrderObjects>Rotate.FreeRotateTransfoririib.LockPositionEditContourLabels.Addb*Measure个DigitsbreakApartLayer11WIEa?eLayer

12、.,EmptyBastLayer.,.CentaurL科w”Po竞Layer.ClassedPostLayer.,irnagcLsycr.MShadedReliefLayer.1GridVectorLayer.-GridVectorLayer,H'A'atershedLayer.3DWireframeLayer.30SurfaceLayer”.LeftAxbRightAxisFopAxisBottomAxisIAjcisScaleBarProfile图8为冲淤图添加工程背景图Thenewormodifiedlayerexceedsthecurrentmaplimits.以昨li

13、mitsare：X:49891S.19134368to501467.7134723SV：4?(JW47,Q215175to4309723.46D。youwanttoadjustthemaplimitstoincludethislayer?星CO香(N)图9如果背景图范围大于冲淤图，则会弹出该对话框，选择否，使范围以冲淤图为准8、添加冲淤厚度值（此处是我会的方法，应该有其他方法？）冲淤厚度值就是冲淤变化的数值，即把数字尽量均匀的添加到冲淤变化图上。首先要生成dat格式的厚度值数据文件，选才IGrid->Convert,将冲淤图的grd文件转换为dat文件，此时dat文件中所有白化区域的网格

14、会填充无效值，如图无效值的行去掉，只保留有效值。10所示。编程将这些含有文件时德箱式（3苣看凹ft4988344308923L7O141E+O3E49883543089231.7O141E+O3849斶3643039231.7O141E4O3849883743085231.TO141M38498S3843039231.7O141EM3849“3s43009231.70141540304988404308923L70141E-+O384988414303923L70141EO38498s424300923L7014山403849跪434908923L加141E的咫49盟4443089231.7O

15、141E4O3849S84543089231.TO141E+03849884643089231-7O141E+O3849884743089231.7O141E+O384908484308923L7O141E+O9849884943089231.70141E-KJ384988504308923l.TO141E-K)38图10冲淤厚度值，白化区域均填充为无效值，需要过滤掉postlayer,在其属性管理器中设置属Vcoordinatesf.1ColumnIBSymbolLColuimnCJAngleINoneDefauItSymbol处理完后，类似添加背景图一样，为冲淤图添加一个性，如图11、图1

16、2所示。-.：ri.Defaultangledegr'FirstrowLastrowFrequencySymbolSizeSizingrnethodSymbolsize97947800Fixed引工好：21AH图11设置数据列、显示密度、字体大小；Frequency最重要，控制显示密度，多试几次选最佳效果gLabelsWorksheetcollLimnColumnCAngle(degrees)0PlaneX-Vplane0FontProperties±1Labe-IFornnatd.dPositionrelativetoAbove叵L«ad«rLines图

17、12设置数据列和字体、格式等属性9、遮挡添加的背景图、冲淤厚度值等可能会超出有效区域，此时可画一个多边形并填充为白色，从而将其遮挡住。10、出图制图完成后(保存为srf文件，以备修改)，选才iFile->Export将地图输出为需要的格式如jpg,选择合适的分到?率。如图13所示。430950043094004309300430920D43091004309000430890042510.50.2,0山2-05T-2.5Y11/制匚1百。W4a。41砧P7梗201406201303bl.-014.1Jfll-心4-A1N口0!11>D.301-4.1,，-心宜.t0J5%-ai力134-QI<*aj+1u箝»41*345-4.34.1_,rr-，»41知M心01&*qn4*a.4，AN刘&li,“一4J.°3二二,引J'y1船114J”/工.包1.4.1.-j62211”-r.3.jjj