Export třídy prvků do WKT souboru

[MartinaPetrovičová]

Dobrý den,

prosím o radu, jak vyexportuji třídu prvků typu polygon do formátu WKT?

Chtěla bych vytvořit polygon pro definování územní části DTI v ArcGIS Pro 3.2.2 a načíst ho při definování územní části DTI v portálu DMVS. Je možné tam načíst soubor ve formátu WKT pro definování územní části DTI. 

Nebo je na to vytvořena nějaká funkce v rámci produktu DTM Connect?

Děkuju za radu.

Petrovičová

[JaroslavŠkrobák · Accepted solution]

Dobrý den, ano, taky jsem pro definování rozsahu do DMVS pomohl tímto skriptem v arcpy.

## cesta ke zdrojovému souboru, kde je ta hranice obce nebo požadovaný polygon

obec = "\\\\MMJ3-GIS2\\gisdata\\ruian\\obec_p.shp"

# # tady si v části where_clause definju, aby mi skript vygeneroval do WKT jen obec Jihlava

obec_p_Layer = "obec_p_Layer"
arcpy.MakeFeatureLayer_management(in_features=obec, out_layer="obec_p_Layer", where_clause=""""nazev" = 'Jihlava'""")

## tady už je samotný výpis, který to ukáže, to si zkopírujete, vložíte do poznámkového bloku  uložíte jako wkt soubor

for row in arcpy.da.SearchCursor(obec_p_Layer, ["SHAPE@WKT"]):
    print row[0]

A takto mi to do DMVS prošlo.

[LukasHampl1]

"""
Script documentation
- Tool parameters are accessed using arcpy.GetParameter() or 
                                     arcpy.GetParameterAsText()
- Update derived parameter values using arcpy.SetParameter() or
                                        arcpy.SetParameterAsText()
"""
import arcpy
def script_tool(param0, param1):
    counter = 0
    for row in arcpy.da.SearchCursor(param0, ["OID@", "SHAPE@WKT"], param1):
        if (counter == 0):
            counter = counter + 1
            return(row[1])
    
if __name__ == "__main__":
    param0 = arcpy.GetParameterAsText(0)
    param1 = arcpy.GetParameterAsText(1)
    param2 = arcpy.GetParameterAsText(2)
    wkt = script_tool(param0, param1)
    f = open(param2, "w")
    f.write(wkt)
    f.close()