使用python调用surfer绘制等值线图的方法属性问题。

qianz98

几年前用VB搞过surfer自动绘制等值线图。最近尝试用PYTHON来实现类似功能，可以实现插值、白化、合并图层。工作环境是python3.6 + surfer16 。但是在调整等值线属性时候遇到问题，编写如下代码：
import win32com.client
import os

def changedate():
    app = win32com.client.gencache.EnsureDispatch('Surfer.Application')
    Plot = app.Documents.Add(1)
    app.Visible = True
    xpath=os.path.dirname(os.path.realpath(__file__))
    datafile=os.path.join(xpath,'temp.dat')
    gridfile=os.path.join(xpath,'temp.grd')
    gridoutf=os.path.join(xpath,'temp-out.dat')
    blankgrid=os.path.join(xpath,'temp_blank.grd')
    outlinefile=os.path.join(xpath,'ZJ_outline.bln')
    LS=os.path.join(xpath,'windlevel.lvl')
    colormap=os.path.join(xpath,'windcolor.clr')
    basefile=os.path.join(xpath,'ZjBack2.bln') #基底图
    postfile=os.path.join(xpath,'StaPotst.dat') #张贴图
    export_img=os.path.join(xpath,'temp.jpg') #输出图像
    #进行插值
    app.GridData(DataFile=datafile,
                           xCol=1,
                           yCol=2,
                           zCol=3,
                           ShowReport=False,
                           Algorithm=2,
                           xMin=117.5,
                           xMax=123.5,
                           yMin=27,
                           yMax=31.5,
                           NumCols="1201",
                           NumRows="901",
                           OutGrid=gridoutf,
                           OutFmt=4)
   
    #白化网格文件
    app.GridBlank(InGrid=gridfile,
            BlankFile=outlinefile,
            OutGrid=blankgrid)
   
   
    ###增加等值线图
    MapFrame1 = Plot.Shapes.AddContourMap(GridFileName=blankgrid)
    #改变等值线图边界范围
    MapFrame1.SetLimits(xMin=118, xMax=123, yMin=27, yMax=31.4)
    MapFrame1.xLength=6
    MapFrame1.yLength=5

    ContourLayer = MapFrame1.Overlays(1)
    ContourLayer.Name = "等值线图"   #修改等值线图名称

    #ContourLayer.Levels.Set(LS)

    ###增加基底图
    MapFrame2 = Plot.Shapes.AddBaseMap(ImportFileName=basefile)
    BaseLayer1 = MapFrame2.Overlays(1)
    BaseLayer1.Name = "基底图"   #修改等值线图名称

    ###增加散点图
    MapFrame3 = Plot.Shapes.AddPostMap(postfile,LabCol=3)
    PostMap = MapFrame3.Overlays(1)
    PostMap.Name = "张贴图"

    ###合并图层生成一个新图层
    MapFrame1.Selected=True
    MapFrame2.Selected=True
    MapFrame3.Selected=True
    NewMapFrame = Plot.Selection.OverlayMaps()
    Plot.Shapes.SelectAll()
    Plot.Export(export_img) #导出图片
   
    #关闭调用程序
    Plot.Close()
    app.Quit()

代码可以执行，能够生成等值线图。
当我打算修改ContourLayer.Levels.Set(LS)
ContourLayer的的Levels属性的时候发现抱错。但是ContourLayer.Name 这个属性是可以修改的。
结论是：
1、PYTHON不是SURFER并不是surfer推荐的开发语言，发现有些语法和对象属性在自带的脚本编辑器中是可用的，但是换成python就不行。网上可参考的PYTHON代码较少。
2、感觉目前用PYTHON还是没法完美的调用surfer绘制等值线图，也用matplotlib库绘制过等值线图，但是插值的方法和效果和克里金法差距还是蛮大。


如果有其他人有更好的办法可以在此交流。

chengxf

要调用 Surfer 的绘图功能，用它的脚本语言，我的经验是：
1、目标计算机安装 Surfer13。没有为什么，上面说了是我的经验。
2、一定要参考 Surfer14 或以上版本的帮助文档，Surfer13 及其以下版本的帮助文档不能用。
3、目标计算机上，Surfer13 一定要能够运行，不能出现要求输入注册码的对话框。
第 1、3 两条，是保证你的程序能够运行的必要条件，第 2 条，那是对你编程的指导。你目前提出的问题，也可以在 Surfer14 的帮助文档里找到答案。
当然，你的程序不能有Bug。

尽头的尽头

都用python了，何必调用surfer呢

qianz98

感谢程老师的建议。

qianz98

尽头的尽头 发表于 2019-6-10 18:05
都用python了，何必调用surfer呢

因为我感觉surfer的克里金插值算法比较好，python 调用surfer的方法走不通，所以我现在用matlotlib库的插值算法，它的方法就比较少了，一共就3种，并且感觉都没有克里金效果好。

holz

楼主这个问题只要好好看帮助文件就可以了。
levels.set() 方法要求的参数是一个数组。
楼主其实想载入一个 .lvl 文件，这时候应该用 levels.loadfile() 方法。

PS：对于用了 python 就无需用 surfer 的说法，我本人持有异议，每个人需求不一样，而且 surfer 自有它的优势所在。虽然我自己用的自动化代码多用 surfer 自带的 scripter，但也有很多时候，VB、C、python 等语言写 surfer 自动化也很必要。

PPS：程老师的意见非常好，surfer 的帮助文件对于其二次开发简直不可或缺。不过我补充两点，一是要用对应版本的帮助，因为每个版本都有可能出现方法、属性、参数的变化；二是用其他语言对 surfer 做二次开发，可以先在 scripter 测试一下自动化代码，对于不熟悉 surfer automation 的人，这可以事半功倍。

qianz98

holz 发表于 2019-6-14 09:31
楼主这个问题只要好好看帮助文件就可以了。
levels.set() 方法要求的参数是一个数组。
楼主其实想载入一 ...

感谢回复，我纠结的不是levels的各种方法的使用，而是python根本不支持levels这个对象，它不仅不支持levels，语句MapFrame.FillContours=1 这样的FillContours对象也不支持。ContourLayer.LevelMethod里的LevelMethod也不支持。另外问一下，如果 scripter 编辑器里面都调试好的文件。不用程序，能否命令行直接执行？

chengxf

qianz98 发表于 2019-6-18 15:15
感谢回复，我纠结的不是levels的各种方法的使用，而是python根本不支持levels这个对象，它不仅不支持leve ...

一般用 C 语言，我用 Pascal，也不错。这个 Python，不懂。帮不到你。
还有，你说的：“另外问一下，如果 scripter 编辑器里面都调试好的文件。不用程序，能否命令行直接执行？”
请将“能否”两字，还有问号，都去掉，这就是回答。

qianz98

chengxf 发表于 2019-6-18 16:06
一般用 C 语言，我用 Pascal，也不错。这个 Python，不懂。帮不到你。
还有，你说的：“另外问一下，如 ...

谢谢

1. "C:\Program Files\Golden Software\Surfer 16\Scripter.exe" -x BaseMap.bas

复制代码

holz

qianz98 发表于 2019-6-18 15:15
感谢回复，我纠结的不是levels的各种方法的使用，而是python根本不支持levels这个对象，它不仅不支持leve ...

刚刚有点时间，试验了一下，一切正常，并没有不支持 levels 对象、方法的现象出现。
我的测试代码：

1. '''
   
2. python 3.7 + win32com + surfer 16.5
   
3. '''
   
4. import win32com.client as w32c
   
5. 
   
6. def main():
   
7.     # 获取已经运行的 Surfer 实例
   
8.     # gsApp = w32c.GetObject(Class = 'Surfer.Application')
   
9. 
   
10.     # 创建新的 Surfer 实例
    
11.     gsApp = w32c.Dispatch('Surfer.Application')
    
12.     gsApp.visible = True
    
13.     Plot = gsApp.documents.add(1)
    
14.     MapFrame = Plot.Shapes.AddContourMap(gsApp.Path + '\Samples\demogrid.grd')
    
15.     ContourLayer = MapFrame.Overlays(1)
    
16.     ContourLayer.Name = '等值线'
    
17. 
    
18.     # 首先要设置等值线填充
    
19.     ContourLayer.FillContours = True
    
20.     # 然后设置的带填充等级属性才能看到
    
21.     ContourLayer.Levels.LoadFile(gsApp.Path + '\Samples\Rainbow.lvl')
    
22. 
    
23. if __name__ == '__main__':
    
24.     main()
    

复制代码