surfer制作沖淤變化圖

1、沖淤變化圖制圖要點(diǎn)注意：在研究一個(gè)新的區(qū)域的沖淤變化時(shí)，最好首先去確定沖淤變化的研究范圍，根據(jù)這個(gè)范圍取點(diǎn)、白化等，這很重要1、原始數(shù)據(jù)的處理：水深點(diǎn)、背景底圖目的是得到可以在surfer中處理的數(shù)據(jù)，主要包括水深數(shù)據(jù)和工程圖（范圍、建筑、標(biāo)記，用作背景）原始數(shù)據(jù)：一般為各種cad工程圖（如果有處理過的水深點(diǎn)數(shù)據(jù)就可以直接用），其中用到的主要數(shù)據(jù)是各時(shí)間段某個(gè)區(qū)域的水深數(shù)據(jù)處理方法：（1）檢查cad圖中水深點(diǎn)數(shù)據(jù)的屬性（特性），有時(shí)有些數(shù)據(jù)點(diǎn)的屬性中沒有水深值，則必須提取出圖中標(biāo)記在點(diǎn)上的那個(gè)值，此時(shí)可用韓志遠(yuǎn)給的插件提取，得到TXT文件，如圖1所示。（2）如果是用作背景的工程范圍圖，則可直接

2、另存為dxf格式的數(shù)據(jù)，然后在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數(shù)據(jù)中提取的水深點(diǎn)數(shù)據(jù)2、用水深數(shù)據(jù)生成網(wǎng)格源數(shù)據(jù)：水深點(diǎn)數(shù)據(jù)，dat、txt等格式

4、處理方法：（1）打開surfer,點(diǎn)擊Grid->Data,彈出后選擇水深數(shù)據(jù)文件。（2）根據(jù)水深文件的具體情況及制圖的需求對(duì)數(shù)據(jù)進(jìn)行網(wǎng)格化，保存為.grd網(wǎng)格數(shù)據(jù)，如圖2、圖3所示。（3）對(duì)各個(gè)時(shí)期各個(gè)工程區(qū)域的水深數(shù)據(jù)都進(jìn)行網(wǎng)格化注意：sufer計(jì)算沖淤變化時(shí)要求參與計(jì)算的兩個(gè)grd文件有相同大小和數(shù)量的網(wǎng)格，即他們的xy坐標(biāo)范圍和節(jié)點(diǎn)數(shù)量（node）都相同，根據(jù)此設(shè)置圖3中的參數(shù)。（可參考網(wǎng)格化文件夾中的參數(shù)截圖）圖2選擇水深文件后注意文件內(nèi)點(diǎn)數(shù)據(jù)的排序圖3注意XYZ對(duì)應(yīng)的數(shù)據(jù)列，根據(jù)點(diǎn)的密度選擇合適的網(wǎng)格搜索半徑，根據(jù)需要設(shè)置網(wǎng)格大小和XY的范圍3、白化處理blankSurfer

5、對(duì)水深數(shù)據(jù)插值計(jì)算生成網(wǎng)格數(shù)據(jù)時(shí)是覆蓋整個(gè)矩形范圍的，無(wú)論每個(gè)網(wǎng)格區(qū)域內(nèi)是否有數(shù)據(jù)點(diǎn)都會(huì)在該網(wǎng)格中進(jìn)行插值計(jì)算。因此計(jì)算沖淤變化之前先要確定是否有需要白化的區(qū)域，如果有就要先制作bln白化文件源數(shù)據(jù)：網(wǎng)格化后的grd數(shù)據(jù)處理方法：打開網(wǎng)格化后的grd數(shù)據(jù)（以contourmap的形式），可看到有效工程范圍內(nèi)有等深線，右鍵點(diǎn)擊地圖，彈出菜單中選擇Digitize,此時(shí)可在地圖上點(diǎn)選一個(gè)多邊形，將無(wú)效區(qū)域圍起來(lái)（可能需要多個(gè)多邊形才能將無(wú)效區(qū)域包圍，不過盡可能用一個(gè)多邊形包圍），Digitized對(duì)話框會(huì)記錄所有點(diǎn)的坐標(biāo)，如圖4所示。多邊形繪制玩后，在Digitized對(duì)話框中選擇File->

6、;Saveas,將其彳存為bln文件，完成后還要再用記事本打開改文件，把第一個(gè)坐標(biāo)行復(fù)制一份到最后一行，使多邊形閉合，第一行第一個(gè)數(shù)字表示點(diǎn)數(shù)，加1,第二個(gè)數(shù)字為標(biāo)志位，選擇0,保存，如圖5所示。如果有多個(gè)多邊形，則對(duì)相應(yīng)的bln文件修改后可將其內(nèi)容合并到一個(gè)bln文件中。最后用Grid->Blank對(duì)grd數(shù)據(jù)進(jìn)行白化處理口DigitizedCoordinates-digitiz.gileEditQpti6rl:圖4用數(shù)字化工具選擇需要白化的多邊形區(qū)域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、用網(wǎng)格數(shù)據(jù)計(jì)算沖淤變化圖源數(shù)據(jù)：白化后的grd數(shù)據(jù)（如不需要白化，則用網(wǎng)格化后的grd數(shù)據(jù)），一般分為不同時(shí)期不同區(qū)域處理方法：選才GGrid->Math,彈出的對(duì)話框中先后選擇同一區(qū)域不同時(shí)間的grd數(shù)據(jù)，用時(shí)間小的數(shù)據(jù)減去時(shí)間大的數(shù)據(jù)（如13年-14年），保存計(jì)算結(jié)果，如圖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計(jì)算沖淤變化5、調(diào)整沖淤變化圖在沖淤變化圖（contourmap）的屬性管理器中設(shè)置其填充、等深線級(jí)別和顏色卡等等。顏色卡可保存為clr文件供其他圖使用。調(diào)整顏色條colorscale屬性，使其位置、大小合適。用劃線和文字的方式添加圖例，指示所計(jì)算的沖淤圖的2個(gè)時(shí)間點(diǎn)。6、添加2個(gè)時(shí)間的等深線以contourmap的形式分別添加兩個(gè)白化后的網(wǎng)格數(shù)據(jù)，設(shè)置其

10、屬性為不填充，并設(shè)置顯示線的顏色，標(biāo)記，間隔等等。如圖7所示。GeneralLevelnmetfiodAdvarcedCortourlevel&EdtItIs.FilledComouriIHIIcontoursColorscaleCsntcurItvHs:Displaythecontourlevsedito圖7等深線的屬性設(shè)置，不填充，沒有colorscale,線的顏色和level的級(jí)別7、添加背景圖如圖8所示。立口fontaurs4Ut_Lthl_bln,grd3Mmp刁iRightAxis7hmiLe+tAxis3TopAxis習(xí)nBottomAxis引回Po5t-uut-fi.d

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

12、.,EmptyBastLayer.,.CentaurL科w”Po競(jìng)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如果背景圖范圍大于沖淤圖，則會(huì)彈出該對(duì)話框，選擇否，使范圍以沖淤圖為準(zhǔn)8、添加沖淤厚度值（此處是我會(huì)的方法，應(yīng)該有其他方法？）沖淤厚度值就是沖淤變化的數(shù)值，即把數(shù)字盡量均勻的添加到?jīng)_淤變化圖上。首先要生成dat格式的厚度值數(shù)據(jù)文件，選才IGrid->Convert,將沖淤圖的grd文件轉(zhuǎn)換為dat文件，此時(shí)dat文件中所有白化區(qū)域的網(wǎng)格

14、會(huì)填充無(wú)效值，如圖無(wú)效值的行去掉，只保留有效值。10所示。編程將這些含有文件時(shí)德箱式（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沖淤厚度值，白化區(qū)域均填充為無(wú)效值，需要過濾掉postlayer,在其屬性管理器中設(shè)置屬Vcoordinatesf.1ColumnIBSymbolLColuimnCJAngleINoneDefauItSymbol處理完后，類似添加背景圖一樣，為沖淤圖添加一個(gè)性，如圖11、圖1

16、2所示。-.：ri.Defaultangledegr'FirstrowLastrowFrequencySymbolSizeSizingrnethodSymbolsize97947800Fixed引工好：21AH圖11設(shè)置數(shù)據(jù)列、顯示密度、字體大??；Frequency最重要，控制顯示密度，多試幾次選最佳效果gLabelsWorksheetcollLimnColumnCAngle(degrees)0PlaneX-Vplane0FontProperties±1Labe-IFornnatd.dPositionrelativetoAbove叵L«ad«rLines圖

17、12設(shè)置數(shù)據(jù)列和字體、格式等屬性9、遮擋添加的背景圖、沖淤厚度值等可能會(huì)超出有效區(qū)域，此時(shí)可畫一個(gè)多邊形并填充為白色，從而將其遮擋住。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