surfer脚本编程作图

liumiaomiao

xyangtian 发表于 2020-5-18 20:50
就是Z值的滤波方式吧，。

你好，请问griddata()中设置z值的滤波方式，这么写exclusionfilter:="z<5"，是不是不对啊，你知道怎么写吗。例子里面也没有。

liumiaomiao

xyangtian 发表于 2020-5-18 20:50
就是Z值的滤波方式吧，。

你好，不好意思，请教您一下：griddata()中设置z值滤波范围的参数我写了后，感觉没有作用，成图后还有空白，负值。你遇到过这种情况没？grid z limits ()这个 参数是不是不能在函数中设置阿。菜单中设置网格z的极限值好像更起作用。

lmf6103

来苦点工分来苦点工分

xyangtian

liumiaomiao 发表于 2020-5-22 22:09
你好，请问griddata()中设置z值的滤波方式，这么写exclusionfilter:="z

用的什么编程语言

liumiaomiao

xyangtian 发表于 2020-5-30 15:26
用的什么编程语言

surfer自带的脚本编程

xyangtian

liumiaomiao 发表于 2020-6-2 10:14
surfer自带的脚本编程

1. '生成等值线图
   
2. Public Sub ContourMain()
   
3.     Set SurferApp = CreateObject("Surfer.Application")
   
4.     SurferApp.Visible = False
   
5.     Set Wks = SurferApp.Documents.Open(baseFile)
   
6.    
   
7.     maxValue = Wks.Columns(2).Statistics().Maximum '获取第二列最大值
   
8.     minValue = Wks.Columns(2).Statistics().Minimum
   
9.     maxHeight = Wks.Columns(4).Statistics().Maximum
   
10.     Wks.Close srfSaveChangesNo '不提示直接关闭
    
11.             
    
12.     '网格
    
13.     SurferApp.GridData DataFile:=baseFile, xCol:=xCombo.Text, yCol:=yCombo.Text, zCol:=zCombo.Text, NumCols:=(maxValue - minValue) / xText.Text + 1, NumRows:=(maxHeight - yMinText.Text) / yText.Text + 1, xMin:=minValue, xMax:=maxValue, yMin:=yMinText.Text, yMax:=maxHeight, _
    
14.     ExclusionFilter:=zFilterText.Text, Algorithm:=srfKriging, DupMethod:=srfDupFirst, ShowReport:=False, OutGrid:=baseName & ".grd"
    
15.             
    
16.     '滤波
    
17.     Call SurferApp.GridFilter(baseName & ".grd", srfFilterGaussian, srfFltEdgeReplicate, srfFltBlankLeave, OutGrid:=baseName & ".grd")
    
18.             
    
19.     '转化为dat
    
20.     'Call SurferApp.GridConvert(InGrid:=GridFile, OutGrid:=BaseName + ".dat")
    
21. 
    
22.     '平滑
    
23.     Call SurferApp.GridSplineSmooth(baseName & ".grd", nRow:=2, nCol:=2, Method:=srfSplineInsert, OutGrid:=baseName & ".grd")
    
24.     '加载网格图
    
25.     Set Plot = SurferApp.Documents.Add(srfDocPlot)
    
26.     Set MapFrame = Plot.Shapes.AddContourMap(baseName & ".grd")
    
27.    
    
28.     '叠加图形
    
29.     'Plot.Shapes.SelectAll
    
30.     'Plot.Selection.OverlayMaps
    
31.    
    
32.     Set plotMape = Plot.Selection.Item(1)
    
33.     Call plotMape.SetLimits(xMin:=minValue, xMax:=maxValue, yMin:=yMinText.Text, yMax:=maxHeight)
    
34.     With plotMape
    
35.         .xMapPerPU = scaleText.Text '设置比例
    
36.         .yMapPerPU = scaleText.Text '设置比例
    
37.     End With
    
38.    
    
39.     '填充网格图
    
40.     Set ContourMap = MapFrame.Overlays(1)
    
41.     ContourMap.FillContours = True
    
42.     'ContourMap.ShowColorScale = True
    
43.     ContourMap.SmoothContours = srfConSmoothHigh
    
44.     ContourMap.Levels.LoadFile (path + "1111.lvl")
    
45.    
    
46.     Plot.SaveAs (baseName + ".srf")
    
47.     SurferApp.Quit
    
48. 
    
49. End Sub
    
50. 
    
51. '加载地形
    
52. Public Sub BlankMian()
    
53.     Set SurferApp = CreateObject("Surfer.Application")
    
54.     SurferApp.Visible = False
    
55.     Set Wks = SurferApp.Documents.Open(baseFile + ".srf")
    
56.    
    
57.     '加载地形bln
    
58.     blankFile = path & "地表" & "" & fileName & ".bln"
    
59.     Set blankMap = Plot.Shapes.AddBaseMap(blankFile)
    
60.     Set blankMapSet = blankMap.Overlays(1)
    
61.     With blankMapSet
    
62.         .Fill.Pattern = "Solid"
    
63.         .Fill.ForeColor = srfColorWhite
    
64.     End With
    
65.    
    
66.     '叠加图形
    
67.     Plot.Shapes.SelectAll
    
68.     Plot.Selection.OverlayMaps
    
69.    
    
70.     Set plotMape = Plot.Selection.Item(1)
    
71.     Call plotMape.SetLimits(xMin:=minValue, xMax:=maxValue, yMin:=yMinText.Text, yMax:=maxHeight)
    
72.     With plotMape
    
73.         .xMapPerPU = scaleText.Text '设置比例
    
74.         .yMapPerPU = scaleText.Text '设置比例
    
75.     End With
    
76.    
    
77.     Plot.Save
    
78.     SurferApp.Quit
    
79. End Sub
    
80. 
    
81. '加载底板
    
82. Public Sub BaseMain()
    
83.     Set SurferApp = CreateObject("Surfer.Application")
    
84.     SurferApp.Visible = False
    
85.     Set Wks = SurferApp.Documents.Open(baseFile + ".srf")
    
86.    
    
87.     '加载底板bln
    
88.     baseMapFile = path & "底板" & "" & fileName & ".bln"
    
89.     Set baseMape = Plot.Shapes.AddBaseMap(baseMapFile)
    
90.    
    
91.     Set baseMapSet = baseMape.Overlays(1)
    
92.     baseMapSet.Line.Width = 0.02
    
93.    
    
94.     '叠加图形
    
95.     Plot.Shapes.SelectAll
    
96.     Plot.Selection.OverlayMaps
    
97.    
    
98.     Set plotMape = Plot.Selection.Item(1)
    
99.     Call plotMape.SetLimits(xMin:=minValue, xMax:=maxValue, yMin:=yMinText.Text, yMax:=maxHeight)
    
100.     With plotMape
     
101.         .xMapPerPU = scaleText.Text '设置比例
     
102.         .yMapPerPU = scaleText.Text '设置比例
     
103.     End With
     
104.    
     
105.     Plot.Save
     
106.     SurferApp.Quit
     
107. End Sub

复制代码