Viac

Pridanie vertikálnej referencie v Pythone do množiny údajov o funkcii ArcGIS

Pridanie vertikálnej referencie v Pythone do množiny údajov o funkcii ArcGIS


Vytváram pythonový skript na generovanie terénu. V procese vytváram dve množiny údajov funkcií. Na úspešné vytvorenie tohto FD však potrebujem horizontálne aj vertikálne základy. Dokážem úspešne kódovať horizontálny nulový bod, ale nemôžem zadať vertikálny nulový bod. Aj keď sa kód nemýli, nevytvára požadovanú množinu údajov funkcií. Použil som ESRI WKID na porovnanie horizontálneho a vertikálneho vzťažného bodu, ako je uvedené nižšie.

Tu je riadok kódu pre priestorovú referenciu, ktorý som použil:

sr = arcpy.SpatialReference (horizontal_wkid)

Skúsil som to urobiťsr = arcpy.SpatialReference (horizontal_wkid, vertical_wkid)ale bezvýsledne. Používam ESRI ArcGIS 10.2 a kódujem v Pythone 2.7


Pretože dané priestorové referencie v ArcGIS majú svoje vlastné parametre pre zvislú časť, pokiaľ viem, neexistuje spôsob, ako zvislú časť zaviesť pomocou daných výrobných kódov alebo podobne.

sr = arcpy.SpatialReference (horizontal_wkid, vertical_wkid)

že by si niekto mohol myslieť, že by to malo fungovať, ale nie je. Môžete si však vytvoriť vlastnú priestorovú referenciu vypožičaním potrebných parametrov z danej sady. Ak navštívite túto stránku, existuje príklad podobného dopytu. V otázke uvidíte túto definíciu:

sp = arcpy.CreateSpatialReference_management ("GEOGCS ['GCS_WGS_1984', DATUM ['D_WGS_1984', SPHEROID ['WGS_1984', 6378137.0,298.257223563]], PRIMEM ['Greenwich', 0,0], UNIT VERTCS ['NAVD_1988', VDATUM ['North_American_Vertical_Datum_1988'], PARAMETER ['Vertical_Shift', 0,0], PARAMETER ['Direction', 1,0], UNIT ['Meter', 1,0]];-400 -400 1000000000; -100000 10000 ; -100000 10000; 8,98315284119521E-09; 0,001; 0,001; IsHighPrecision ","#","#","#","#","#"," 0 ")

Prvá časť tu definuje geografické súradnicové systémy ako:

GEOGCS ['GCS_WGS_1984', DATUM ['D_WGS_1984', SPHEROID ['WGS_1984', 6378137.0,298.257223563]], PRIMEM ['Greenwich', 0,0], UNIT ['Degree', 0,0174532925199433]]

za ktorými nasledujú vertikálne systémové parametre (vrátane názvu, jednotky, parametrov definície atď.) ako:

VERTCS ['NAVD_1988', VDATUM ['North_American_Vertical_Datum_1988'], PARAMETER ['Vertical_Shift', 0,0], PARAMETER ['Direction', 1,0], UNIT ['Meter', 1,0]]

Potom by ste mali byť schopní predstaviť tento vlastný priestorový odkaz ako argument pre arcpy.CreateFeatureDataset.