PeterWilson的最新贡献

我目前仍在使用ArcGIS 10.3，必须为我们的一个项目的洪水风险分析任务生成整个达累斯萨拉姆的1m分辨率DTM。第一个问题是，我们输入数据的来源是1m分辨率的轮廓形状文件。我最初从轮廓数据集中创建了一个5m分辨率的DTM，用于初步结果，运行了近5天。这是一个问题，因为我必须将分辨率降低到1m才能运行我的分析。我生成了一个重叠索引网格，允许20%的重叠，并将轮廓数据集切成500块。然后，我使用Python的Multiprocessing模块生成了一个Python模块，在我们的服务器上跨16个核处理瓷砖。整个研究区域的DTM瓷砖在12小时内完成。我在下面附上的Python模块，是我利用Python的多处理的第一个草稿。我现在正在重写它，因为我更好地理解了Multiprocessing模块是如何工作的。注意:我要感谢Duncan Hornby的帖子:创建一个使用多处理的脚本工具，帮助我创建了以下脚本，希望以下内容可以帮助其他试图生成大型DTM的达累斯萨拉姆:100万DTM瓷砖“'创建于2018年2月1日多线程拓扑到光栅(3D) @作者: PeterW ''' # import site-packages and modules import multiprocessing import argparse import arcpy from functools import partial # set environment settings arcpy.env.overwriteOutput = True # check-out extensions arcpy.CheckOutExtension('3D') def dtm_work(tiles_directory, cont): """ Multiprocess Topo to Rasster for multiple DTMs """ index_name = cont[0] inContours = cont[1] try: output_dtm = '{0}\\{1}'.format(tiles_directory, index_name) arcpy.TopoToRaster_3d(inContours, out_surface_raster=output_dtm, cell_size=1) return True except Exception as e: print e.message return False def dtm_multi(index_grid, tiles_directory): """ Create overlapping DTM tiles from contour tiles & index grid """ try: cont_param = [] with arcpy.da.SearchCursor(index_grid, ['SHAPE@', 'PageName', 'Path']) as scur: # @UndefinedVariable for row in scur: index_name = row[1] input_cont = '{0} {1} {2}'.format(row[2], 'ELEVATION', 'Contour') cont_param.append([index_name, input_cont]) func = partial(dtm_work, tiles_directory) arcpy.AddMessage('Sending to pool') # declare the number of cores to use, use 2 less than the max cpu_num = multiprocessing.cpu_count()-2 # create the pool object pool = multiprocessing.Pool(processes=cpu_num) # fire off list to worker function # res is a list that is created with what ever the worker function returns res = pool.map(func, cont_param) pool.close() pool.join() # if an error has occurred report it if False in res: arcpy.AddError('A process\thread failed!') arcpy.AddMessage('Finished Multiprocessing DTM Tiles') except arcpy.ExecuteError: # Geoprocessor through error arcpy.AddError(arcpy.GetMessages(2)) if __name__ == '__main__': parser = argparse.ArgumentParser(description='Create overlapping DTM tiles for mosaic') parser.add_argument('--index_grid', metavar='path', required=True, help='path to input overlapping index grid feature class') parser.add_argument('--tiles_directory', metavar='path', required=True, help='path to output DTM tiles directory') args = parser.parse_args() dtm_multi(index_grid=args.index_grid, tiles_directory=args.tiles_directory) …查看更多

我能够根据我的学习范围使用网格索引功能工具创建一个合理的索引网格。我想弄清楚的是如何创建一个索引网格，有一定比例的重叠。下面将使用Topo to Raster生成DTM。源是整个达累斯萨拉姆的等高线，以测量等高线为基准，分辨率为1m。当我最后一次尝试使用完整的轮廓数据集作为单个特征类运行它时，它以5m分辨率运行了7天。索引网格:重叠索引网格:10% DTM: 5m分辨率 …查看更多

我有一个道路边缘多边形，我想在每个十字路口分割。我偶然发现了下面的帖子:如何缓冲折线而不在垂直处重叠，Richard Fairhurst建议了一个工作流来实现以下目标。不幸的是，我无法跟随他的脚步，因为他在道路中心线的交叉点周围建立了一个缓冲区。我已经附上了我的道路中心线的样本，以及我的当前缓冲区，这是我基于道路中心线和代表道路宽度的字段创建的。任何建议，如何实现一个分裂多边形的基础上的结果，理查德菲尔赫斯特实现将不胜感激。从Richard Fairhurst回复样本数据集:(道路中心线;道路缓冲器(合并);道路交叉口) …查看更多

我试图确定如何计算两个光栅之间的骰子相似系数。其中一个栅格是路面面积的真实结果，第二个栅格是计算机视觉和机器学习(卷积神经网络)的结果。我试图了解骰子相似系数是如何工作的，以便我可以将其复制为Python\NumPy函数，以运行手动确定的所有样本结果和地面真相结果。注:我也想确定下面的栅格描绘的真正;假阳性和假阴性维基百科:Sorensen-Dice Cofficient公式:光栅结果: …查看更多