Editace vektorových dat v ArcPy

Programování v GIS 3

Jan Caha

2025-10-07

Objekty pro práci s vektorovými daty

• SpatialReference - souřadnicový systém
• Geometry - geometrie

SpatialReference

• souřadnicový systém s řadou parametrů
• lze vytvořit nový nebo získat z existujících dat
• identifikace názvem nebo lépe EPSG kódem

crs_1 = arcpy.SpatialReference("Sinusoidal (world)")
crs_2 = arcpy.SpatialReference(5514)
crs_3 = arcpy.Describe("data.shp").spatialReference

Geometry

• velice komplexní objekt reprezentující geometrii
• řada metod a vlastností
• aktuálně nás zajímají funkce pro tvorbu geometrií - arcpy.Point(), arcpy.Multipooint(), acpy.Polyline(), acpy.Polygon()

Funkce

• arcpy.management.CreateFeatureclass() - vytvoření nové vrstvy
• arcpy.management.Copy() - kopírování dat - vrstvy

CreateFeatureclass()

• vytvoření nové vrstvy dle zadaných parametrů

arcpy.management.CreateFeatureclass(
    out_path, 
    out_name, 
    {geometry_type}, {template}, 
    {has_m}, {has_z}, {spatial_reference}, 
    {config_keyword}, {spatial_grid_1}, 
    {spatial_grid_2}, {spatial_grid_3}, 
    {out_alias}, {oid_type}
)

Copy()

• kopie vrstvy

arcpy.management.Copy(
    in_data, 
    out_data, 
    {data_type}, {associated_data}
)

Přístup k vektorovým datům

• pomocí tzv. cursor (kurzor) objektů (databázová terminologie)
• v Pythonu fungují jako kontextové manažery

with cursor(parmeter_1, parameter_2, ...) as cursor:
    for row in cursor:
        pass

• součást modulu arcpy.da (např. arcpy.da.UpdateCursor)

Pozor

Verze volané přímo z arcpy nepoužívat (zastaralé - deprecated).

Takže arcpy.UpdateCursor ne!

Mají i jiné parametry.

Typy kurzorů

• SearchCursor - pro čtení a vyhledávání dat
• InsertCursor - pro vkládání nových dat
• UpdateCursor - pro aktualizaci dat

SearchCursor

SearchCursor (
    in_table, 
    field_names, 
    {where_clause}, 
    {spatial_reference}, 
    {explode_to_points}, 
    {sql_clause}, 
    {datum_transformation}, 
    {spatial_filter}, 
    {spatial_relationship}, 
    {search_order}
)

SearchCursor - parametry

• in_table - cesta k datové vrstvě (soubor nebo i v databázi)
• field_names - seznam polí, které mají být součástí výsledku
• where_clause - část SQL dotazu (označované jako query expression viz dokumentace)
• spatial_reference - souřadnicový systém v němž bude výsledek (provede se trasformace, je-li třeba)
• explode_to_points - rozbití geometrie na body (pokud je True) vzniká záznam (Row) pro každý bod
• sql_clause - SQL prefix a postfix (default je python (None, None))
• spatial_filter - prostorový filtr arcpy.Geometry objekt
• spatial_relationship - prostorový vztah (default je INTERSECT)
• search_order - pořadí filtrů - buď "ATTRIBUTEFIRST" nebo "SPATIALFIRST" - vliv na dobu trvání dotazu

Parametr fields_names

• lze zadat hodnotu jako * pro všechna pole, ale z důvodu výkonu a přehlednosti je lepší zadat jen potřebná pole
• speciální pole - mají označením @
• id prvku získáme jako OID@
• geometrie prvku jako SHAPE@ - vrací objekt Geometry, ale náročné na výkon
  • SHAPE@XY - vrací tuple s X a Y souřadnicemi centroidu geometrie
  • SHAPE@XYZ- viz výše, ale s výškou Z
  • SHAPE@JSON - Esri JSON reprezentace geometrie
  • SHAPE@WKT - WKT reprezentace geometrie
  • SHAPE@WKB - WKB reprezentace geometrie
  • SHAPE@AREA - plocha geometrie
  • SHAPE@LENGTH - délka geometrie

Objekt Row

• objekt reprezentující řádek v tabulce (dokumentace)
• 4 funkce pro přístup k datům
  • getValue(field_name)
  • setValue(field_name, value)
  • isNull(field_name)
  • setNull(field_name)

Přístup k hodnotám Row

• dle specifikace v field_names
• pomocí indexů

row[0]  # první pole
row[1]  # druhé pole

UpdateCursor

UpdateCursor (
    in_table, 
    field_names, 
    {where_clause}, 
    {spatial_reference}, 
    {explode_to_points}, 
    {sql_clause}, 
    {datum_transformation}, 
    {explicit}, 
    {spatial_filter}, 
    {spatial_relationship}, 
    {search_order}
)

• po úpravě je nutno změnu uložit

row[0] = "Nová hodnota"
cursor.updateRow(row)

UpdateCursor - parametry

• datum_transformation - transformace souřadnicových systémů (zřídka potřeba nastavovat manuálně)
• explicit - přepis NULL na defaultní hodnoty (pokud je vrstva má nastavené)

Mazání záznamů

• v podstatě update
• logiku výběru prvků je vhodné řešit v nastavení cursoru
• cursor maže aktivní řádek (určený podmínkou for row in cursor:)
• teoreticky lze mazat i v nějaké podmínce, ale je to neefektivní

cursor.deleteRow()

InsertCursor

InsertCursor (
    in_table, 
    field_names, 
    {datum_transformation}, 
    {explicit}
)

• vložení do tabulky - jako python Tuple (hodnoty v kulatých závorkách)

cursor.insertRow(("hodnota 1", 23, arcpy.Point(1, 2)))

Dotazy?