Создание рыболовной сети из шаблонного класса объектов с использованием ArcPy?

9

Я не могу использовать инструмент arcpy.CreateFishnet_management, поскольку определяя параметр «templateExtent» с помощью шейп-файла, он не заполняет автоматически параметры «originCoordinate» и «yAxisCoordinate».

import arcpy
from arcpy import env
env.overwriteOutput = True
env.workspace = r"D:\Users\julia\erste_aufg"

#Process: Create Fishnet
outFeatureClass = r"D:\Users\julia\erste_aufg\at001l_wien\at001l_wien\wien.shp"
cellSizeWidth = '200'
cellSizeHeight = '200'
templateExtent = r"D:\Users\julia\erste_aufg\at001l_wien\at001l_wien\at001l_wien.shp"

arcpy.CreateFishnet_management(outFeatureClass, "", "", cellSizeWidth, cellSizeHeight, '0', '0', "", "NO_LABELS", templateExtent, "POLYGON")

введите описание изображения здесь

Он работает в ModelBulider, поэтому что-то работает в фоновом режиме ModelBulider, чтобы он мог создавать параметры «originCoordinate» и «yAxisCoordinate», когда у него есть «templateExtent». Как я могу запустить этот инструмент в ArcPy, используя только параметр «templateExtent»?

Я был бы очень рад, если у кого-то есть решение, потому что я нуждаюсь в Fishnet в скриптовом инструменте и не могу обойтись без него, потому что в конце есть цикл, поэтому значения экстента всегда разные. первая часть всего сценария

Нора
источник
Кто-то знает, почему мы добавляем 10 к части решения выше? arcpy.CreateFishnet_management (fc [: - 4] + "_ c200.shp", str (desc.extent.lowerLeft), str (desc.extent.XMin) + "" + str (desc.extent.YMax + 10), " 200 "," 200 "," 0 "," 0 ", str (desc.extent.upperRight)," NO_LABELS "," # "," POLYGON ")
user5956986
Это не дает ответа на вопрос. Как только у вас будет достаточно репутации, вы сможете комментировать любой пост ; вместо этого предоставьте ответы, которые не требуют разъяснений от автора . - Из обзора
Дэн С

Ответы:

14

вот пример. Вам необходимо извлечь ограничивающую рамку из объекта описания.

desc = arcpy.Describe(fc)
arcpy.CreateFishnet_management(fc[:-4]+"_c200.shp",str(desc.extent.lowerLeft),str(desc.extent.XMin) + " " + str(desc.extent.YMax + 10),"200","200","0","0",str(desc.extent.upperRight),"NO_LABELS","#","POLYGON")
radouxju
источник
@@ radouxju, что цель + 10в str(desc.extent.YMax + 10)?
maycca
хороший вопрос. На самом деле не нужно в этом случае. У меня есть привычка добавлять произвольное значение Ymin для построения вертикальной оси, но здесь я использовал Ymax, поэтому это излишне.
Radouxju
4

Вот альтернативный подход, который я использовал для создания нескольких рыболовных сетей в пределах каждого объекта в классе объектов. Переменная search_extents определяет путь к этому классу пространственных объектов, определяя экстенты каждой рыболовной сети, которую я хотел создать. Вращения рыболовной сети не было.

search_extents = "path to extents" 
rows = arcpy.SearchCursor(search_extents)
shapeName = arcpy.Describe(search_extents).shapeFieldName
for row in rows:
    print("Starting Extent" + row.getValue("Extent_Num"))
    feat = row.getValue(shapeName)
    extent = feat.extent
    arcpy.CreateFishnet_management(arcpy.env.workspace + "/extents/extentgrid" + row.getValue("Extent_Num"),str(extent.lowerLeft), str(extent.upperLeft),"0","0","200","200",str(extent.upperRight),"NO_LABELS","#","POLYGON")
    print("Finishing Extent" + row.getValue("Extent_Num"))
GForce
источник
1

Вот код, который я наконец-то успешно начал работать (с помощью примеров выше), чтобы решить проблему, описанную здесь:

    env.workspace = "C:/Holly/Work/Projects/NavigationStudy2019/Data"

    # Fetch each feature from the cursor and examine the extent properties
    for row in arcpy.da.SearchCursor(feature_class, ['SHAPE@', 'id']):
        extent = row[0].extent
        print('Extent of home range {}:'.format(row[1]))
        print('XMin: {}, YMin: {}'.format(extent.XMin, extent.YMin))
        print('XMax: {}, YMax: {}'.format(extent.XMax, extent.YMax))
        arcpy.CreateFishnet_management("fishnet_temp.shp",
                                       str(extent.XMin) + " " + str(extent.YMax),
                                       str(extent.XMin) + " " + str(extent.YMax + 10),
                                       "100",
                                       "100",
                                       "",
                                       "",
                                       "",
                                       "LABELS",
                                       feature_class,
                                       "POLYGON")
Холли Коупленд
источник