155GIS2 / 9. cvičení

← předchozí cvičení • stránky předmětu • další cvičení →

Úvod do skriptování v jazyku Python (arcpy)

Náplň

• Přednáška GIT2
• Tutorial Programování v Pythonu pro ArcGIS 10.0 pro začátečníky

Poznámky k systému ArcGIS

Python Window

View • Python Window

Ukázka interaktivní Python konzole v ArcGIS Pro

arcpy.env.workspace = r'S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb'
obce = arcpy.Describe('ObcePolygony')
len(obce.fields)
pocet = arcpy.GetCount_management('ObcePolygony')
pocet.getOutput(0)

Notebook

Založeno na open source projektu Jupyter.

Insert • New Notebook

Užitečné odkazy

• Python in ArcGIS Pro
• A quick tour of creating tools with Python
• What is ArcPy?
• ArcGIS API for Python
• Overview of the ArcGIS API for Python
• ArcGIS Pro Python reference

PyCharm

Cesta k Python interpretu: C:\Program Files\ArcGIS\Pro\bin\Python\envs\arcgispro-py3\python.exe

Úlohy

V následujících úlohách budeme, pokud nebude uvedeno jinak, používat data z datové sady ArcCR500. Tato datová sada je poskytovaná firmou ArcData zdarma.

V učebně B-870 jdou data nainstalována v adresáři

S:\K155\Public\data\ArcGIS\ArcCR500 3.3

Přístupová práva k tomu adresáři jsou omezena pouze na čtení! Nelze tedy datové vrstvy z ArcCR500 upravovat (včetně např. vytváření pyramid).

Další ukázky skriptů zde (zip).

1

Výpis názvů atributů

• Describe
  • FeatureClass properties
  • Table properties

import arcpy
arcpy.env.workspace = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb"
 
try:
   data = arcpy.Describe("ObcePolygony")   
   for i in range(len(data.fields)):
      print(data.fields[i].name)
except Exception as e:
   print(e)

Tip:

1.

for field in data.fields:
    print(field.name)

2. Pokud budete skript spouštět mimo notebook, tak použijte

arcpy.AddMessage()

2

Výpis typu geoprvků

import arcpy

inputFC = arcpy.GetParameterAsText(0)

shapeType = arcpy.Describe(inputFC).shapeType
arcpy.AddMessage("Shape Type: {0}".format(shapeType))

3

Výpis obcí, které mají počet obyvatel vyšší než Liberec

import arcpy
arcpy.env.workspace = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb"

try:
  rows = arcpy.SearchCursor("ObcePolygony", "NAZ_OBEC = 'Liberec'")   #vytvori cursor
  for row in rows:
      pocet = row.POCET_OBYV # pocet obce Liberec
      break
  rows = arcpy.SearchCursor("ObcePolygony")
  for row in rows:
      if row.POCET_OBYV is not None and row.POCET_OBYV > pocet:
         print("{0} - {1}".format(row.NAZ_OBEC, row.POCET_OBYV))
except Exception as e:
   print(e)

4

Vytvoření souborů Shapefile pro jednotlivé kraje

import os
import arcpy
arcpy.env.overwriteOutput = True
arcpy.env.workspace = r"C:\Users\landamar\Documents\u04"
os.makedirs(arcpy.env.workspace)

kraje = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb\KrajePolygony"

rows = arcpy.SearchCursor(kraje)   #vytvori cursor
for row in rows:
      nazev = row.KOD_CZNUTS3
      print("{0}...".format(nazev))
      kraje_Select_shp = os.path.join(arcpy.env.workspace, nazev + ".shp")
      arcpy.Select_analysis(kraje, kraje_Select_shp, "\"KOD_CZNUTS3\" = '{}'".format(nazev))

5

Výpis souřadnic lomových bodů uzemí obce 'Lhota'

import arcpy

arcpy.env.workspace = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb"
 
obce = "ObcePolygony"

try:
   data = arcpy.Describe(obce)
   shapename = data.ShapeFieldName
   rows = arcpy.SearchCursor(obce, '"NAZ_OBEC" = \'Lhota\'')
   for row in rows:
       prvek = row.getValue(shapename)
       arcpy.AddMessage("Prvek {0}: ".format(row.NAZ_OBEC))
       partnum = 0
       for part in prvek:
           arcpy.AddMessage("Part {0}: ".format(partnum))
           for bod in prvek.getPart(partnum):
               arcpy.AddMessage("{0} {1}".format(bod.X, bod.Y))
           partnum += 1  
except Exception as e:
   print(e)

---

Další úlohy najdete v cvičení k předmětu 153GIS2.