155GIS2 / 9. cvičení: Porovnání verzí
m (→Náplň) |
m (→Úlohy) |
||
Řádek 54: | Řádek 54: | ||
** [https://pro.arcgis.com/en/pro-app/arcpy/functions/table-properties.htm Table properties] | ** [https://pro.arcgis.com/en/pro-app/arcpy/functions/table-properties.htm Table properties] | ||
< | <syntaxhighlight lang="python"> | ||
arcpy.env.workspace = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb" | arcpy.env.workspace = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb" | ||
Řádek 63: | Řádek 63: | ||
except Exception as e: | except Exception as e: | ||
print(e) | print(e) | ||
</ | </syntaxhighlight> | ||
Tip: | Tip: | ||
< | <syntaxhighlight lang="python"> | ||
for field in data.fields: | for field in data.fields: | ||
print(field.name) | print(field.name) | ||
</ | </syntaxhighlight> | ||
=== 2 === | === 2 === | ||
Řádek 76: | Řádek 76: | ||
; Výpis typu geoprvků | ; Výpis typu geoprvků | ||
< | <syntaxhighlight lang=python> | ||
inputFC = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb\ObcePolygony" | inputFC = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb\ObcePolygony" | ||
shapeType = arcpy.Describe(inputFC).shapeType | shapeType = arcpy.Describe(inputFC).shapeType | ||
print("Shape Type: {}".format(shapeType)) | print("Shape Type: {}".format(shapeType)) | ||
</ | </syntaxhighlight> | ||
=== 3 === | === 3 === | ||
Řádek 86: | Řádek 86: | ||
; Výpis obcí, které mají počet obyvatel vyšší než Liberec | ; Výpis obcí, které mají počet obyvatel vyšší než Liberec | ||
< | <syntaxhighlight lang="python"> | ||
arcpy.env.workspace = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb" | arcpy.env.workspace = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb" | ||
Řádek 100: | Řádek 100: | ||
except Exception as e: | except Exception as e: | ||
print(e) | print(e) | ||
</ | </syntaxhighlight> | ||
=== 4 === | === 4 === | ||
Řádek 106: | Řádek 106: | ||
; Vytvoření souborů Shapefile pro jednotlivé kraje | ; Vytvoření souborů Shapefile pro jednotlivé kraje | ||
< | <syntaxhighlight lang="python"> | ||
import os | import os | ||
arcpy.env.overwriteOutput = True | arcpy.env.overwriteOutput = True | ||
Řádek 121: | Řádek 121: | ||
kraje_Select_shp = os.path.join(arcpy.env.workspace, nazev + ".shp") | kraje_Select_shp = os.path.join(arcpy.env.workspace, nazev + ".shp") | ||
arcpy.Select_analysis(kraje, kraje_Select_shp, "\"KOD_CZNUTS3\" = '{}'".format(nazev)) | arcpy.Select_analysis(kraje, kraje_Select_shp, "\"KOD_CZNUTS3\" = '{}'".format(nazev)) | ||
</ | </syntaxhighlight> | ||
=== 5 === | === 5 === | ||
Řádek 127: | Řádek 127: | ||
; Výpis souřadnic lomových bodů uzemí obce 'Lhota' | ; Výpis souřadnic lomových bodů uzemí obce 'Lhota' | ||
< | <syntaxhighlight lang="python"> | ||
arcpy.env.workspace = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb" | arcpy.env.workspace = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb" | ||
Řádek 147: | Řádek 147: | ||
except Exception as e: | except Exception as e: | ||
print(e) | print(e) | ||
</ | </syntaxhighlight> | ||
---- | ---- | ||
''Další úlohy najdete v cvičení k předmětu [[153GIS2 / 10. cvičení#Úkoly|153GIS2]].'' | ''Další úlohy najdete v cvičení k předmětu [[153GIS2 / 10. cvičení#Úkoly|153GIS2]].'' |
Verze z 14. 4. 2021, 08:36
Úvod do skriptování v jazyku Python (arcpy)
Náplň
- Přednášky 155GIT1: Úvod do jazyka Python, Algoritmizace, čtení a zápis souborů
- Tutorial Programování v Pythonu pro ArcGIS 10.0 pro začátečníky (studentský projekt)
Poznámky k systému ArcGIS
Python Window
View • Python Window
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
Ú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ů
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:
for field in data.fields:
print(field.name)
2
- Výpis typu geoprvků
inputFC = r"S:\K155\Public\data\ArcGIS\ArcCR500 3.3\AdministrativniCleneni_v13.gdb\ObcePolygony"
shapeType = arcpy.Describe(inputFC).shapeType
print("Shape Type: {}".format(shapeType))
3
- Výpis obcí, které mají počet obyvatel vyšší než Liberec
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
arcpy.env.overwriteOutput = True
arcpy.env.workspace = r"C:\Users\landamar\Documents\u04"
if not os.path.exists(arcpy.env.workspace):
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'
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:
print("Part {0}: ".format(partnum))
for bod in prvek.getPart(partnum):
print("{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.