起因是这样一段对话~~~~
经过一番研究,研究出大概的思路,先将有经纬度的表中的数据筛选出表并生成xy事件,接着利用核密度工具生成栅格,最后呢裁剪栅格通过mapping包出图。
多线程第一次用,一开始打算用数据库存放数据表和栅格,不过锁表比较难解决,于是选择了简单粗暴地在文件地理数据库(gdb)里处理,每个文件使用一个gdb,这样速度快也不用考虑会被锁表,处理完以后就直接删库(跑路)。 多线程中尽量不要有删除操作,容易冲突,把删除都放在单线程中这样不容易出错。 MyGPpool = multiprocessing.Pool(19) results = MyGPpool.map(batch, parameters) MyGPpool.close() MyGPpool.join() worker中使用别的脚本要使用importScripts导入,当需要传递参数可以transferable object(通讯时会比较快,参考worker使用),本次项目中使用Float64Array.看几个初始化范例 mxd = arcpy.mapping.MapDocument(mxd_file) df = arcpy.mapping.ListDataFrames(mxd, '图层')[0] lyr = arcpy.mapping.ListLayers(mxd)[0] datasetName = os.path.basename(data_raster) dataRasterWorkspace = os.path.dirname(data_raster) if lyr.supports("DATASOURCE"): lyr.replaceDataSource(dataRasterWorkspace,"FILEGDB_WORKSPACE", datasetName) df.extent = extent df.panToExtent(extent) arcpy.mapping.ExportToPNG(mxd, out_file, df, color_mode='24-BIT_TRUE_COLOR', world_file=True,background_color="No Color", transparent_color="No Color") else: arcpy.AddMessage("error") arcpy.RefreshActiveView() del mxd, df, ly |
|