Python+Surfer自动化

大老鹰

文件来自官网，仅供参考

rabin_xu

Surfer都支持python了啊，孤陋寡闻了

931404656

我也是第一次听说

rgyxtq2009

多谢分享。

lon91ong

收费不厚道

1. 
   
2. Surfer can be called from a Python script. Although writing scripts in python are not officially supported, we do have some examples to share.
   
3. 
   
4. 
   
5. 
   
6.     Here is an example that shows how to start Surfer:
   
7. 
   
8.     import win32com.client def main():
   
9.     app = win32com.client.Dispatch("Surfer.Application")
   
10.     plot = app.Documents.Add(1)
    
11.     app.Visible = True main()
    
12. 
    
13.     Here is an example that shows how to call Surfer:
    
14. 
    
15.     Surfer = win32com.client.Dispatch("Surfer.Application")
    
16.     Grid = Surfer.NewGrid()
    
17.     Grid.LoadFile("c:/test/test.grd", HeaderOnly=True)
    
18.     print "BlankValue" + str(Grid.BlankValue)
    
19.     print "NumCols" + str(Grid.NumCols)
    
20.     print "NumRows" + str(Grid.NumRows)
    
21. 
    
22.     Here is an example of how to create and edit a vector map:
    
23. 
    
24.     import constants, Dispatch, CastTo
    
25.     import os.path def rgb(r, g, b): return b << 16 / g < 0:
    
26.     level.ShowLabel = True
    
27.     mapframe = plot.Shapes.AddVectorMap(GridFileName1=grid2, GridFileName2=grid3,  CoordSys=constants.srfVecPolar, AngleSys=constants.srfVecAngle)
    
28.     vectormap = CastTo(mapframe.Overlays(1), 'IVectorMap')
    
29.     orange = rgb(255, 102, 0)
    
30.     vectormap.ColorMap.SetNodes(Positions=[0.0, 1.0], Colors=[orange]*2)
    
31.     vectormap.ColorScaleMethod = constants.srfVecMagnitude
    
32.     plot.Shapes.SelectAll()
    
33.     plot.Selection.OverlayMaps()
    
34.     plot.SaveAs(output)
    
35.     plot.Export(exportfile)
    
36. 
    
37.     Here is a customer example of how to grid data:
    
38. 
    
39.     import glob
    
40.     import datetime
    
41. 
    
42.     #call this script like this:
    
43.     #\programs\python25\python.exe krig_data.py
    
44. 
    
45.     #get an instance of the Surfer application  
    
46.     #Surfer = win32com.client.Dispatch('Surfer.Application')
    
47.     lFile = glob.glob('wl*.csv')
    
48.     for i in range(len(lFile)):
    
49. 
    
50.     #the min and max listed below are center coordinates while gdal_grid uses edge
    
51. 
    
52.     #coordinates
    
53.     sFilePrefix = lFile[i][:lFile[i].index('.')]
    
54.     Surfer.GridData(DataFile = lFile[i], xCol = 1, yCol = 2, zCol = 3, NumCols = 4133,  NumRows = 3017, xMin = 455355, xMax = 579315, yMin = 1954585, yMax = 2045065,  ShowReport = False, DupMethod = win32com.client.constants.srfDupAvg, OutFmt = win32com.client.constants.srfGridFmtBinary, OutGrid = '%s.grd' %sFilePrefix) print 'finished kriging %s at %s' % (lFile[i], datetim{source}
    
55. 
    
56. 
    
57.     Here is an example on how to call Surfer and grid your datafile to create a contour map.
    
58. 
    
59.     def main():
    
60.         app = win32com.client.gencache.EnsureDispatch('Surfer.Application')
    
61.         Plot = app.Documents.Add(1)
    
62.         app.Visible = True
    
63.         DataFile = "C:\Program Files\Golden Software\Surfer 12\Samples\demogrid.dat"
    
64.         OutFile = "C:\Program Files\Golden Software\Surfer 12\Samples\outgrid.grd"
    
65.         app.GridData (DataFile=DataFile, Algorithm = win32com.client.constants.srfMinCurvature, NumRows=150, NumCols=150, ShowReport=False, OutGrid= OutFile)
    
66. 
    
67.         #Creates a contour map and assigns the map frame to the variable "MapFrame"
    
68. 
    
69.         MapFrame = Plot.Shapes.AddContourMap(GridFileName=OutFile)
    
70. 
    
71.         #Changes the limits and scale of the map
    
72.         MapFrame.SetLimits (xMin=0.5, xMax=4.5, yMin=0.5, yMax=3.5)
    
73.         MapFrame.xLength=6
    
74.         MapFrame.yLength=4
    
75. 
    
76.         #Declares ContourMap as an Object and assigns the contour map to variable "ContourMap"
    
77.         ContourMap = MapFrame.Overlays(1)
    
78. 
    
79.     main()
    
80. 
    

复制代码

xyangtian

十分感谢楼主，楼主在那个官网找的。

xyangtian

lon91ong 发表于 2017-12-8 09:38
收费不厚道

我用的comtypes调用，这个GirdData方法报异常

xyangtian

lon91ong 发表于 2017-12-8 09:38
收费不厚道

执行这个GridData方法总是报错啊，

xyangtian

lon91ong 发表于 2017-12-8 09:38
收费不厚道

执行这个GridData方法总是报错啊，