Případ chybějících předpon: ArcGIS ...Geometries

JoshuaBixby

<\/HEAD>

V \/blogs\/tilting\/2014\/07\/25\/semantic-overloading<\/A> se dotýkám jednoho z hlavních motivů pro psaní blogu Tilting at Globes<\/A>:<\/P>

Věřím, že sémantika je důležitá ve všech aspektech života.  Ať už v právu, medicíně, vědě, podnikání, informačních technologiích nebo v jakémkoli jiném oboru; mít společný jazyk není příliš užitečné, pokud neexistuje společné porozumění slovům, která jazyk tvoří.  Jelikož se jazyky vyvíjejí, udržování společného porozumění slovům v průběhu času je neustálou výzvou.<\/P><\/BLOCKQUOTE>

<\/P>

Ideálně by společné porozumění znamenalo, že slovo nebo termín má jediný význam a tento jediný význam je znám a pochopen všemi, kdo slovo nebo termín používají.  Bohužel realita není vždy ideální a někdy kontext slova nebo termínu hraje velkou roli v jeho významu.  Například způsob, jakým je termín "large-scale" aplikován na mapy, může být zdánlivě nelogický ve srovnání s tím, jak je aplikován na akce, události a běžné objekty.  Když člověk pochopí, že měřítko mapy se vztahuje k reprezentaci dat na mapě, tj. poměr vzdálenosti na mapě k vzdálenosti na zemi, pomáhá to vysvětlit, proč "large-scale" mapy pokrývají malé geografické oblasti místo velkých geografických oblastí.  Jistě, large-scale mapa by mohla pokrýt velkou geografickou oblast, pokud by byla vytištěna na obrovském médiu, ale předpokládám typické velikosti tisku.<\/EM><\/P><\/P>Přejdeme-li od map ke geometriím, zmiňuji případ prostorových operátorů v

\/blogs\/tilting\/2015\/05\/14\/whats-within-when-esri-clementini<\/A>.  V tomto blogovém příspěvku poukazuji na to, jak odpověď na otázku, zda je jedna geometrie "within" jiné geometrie, může záviset na tom, koho se zeptáte, nebo přesněji řečeno, čí definice "within" se použije k zodpovězení otázky.  Často jsou kvalifikátory buď neúplné nebo úplně vynechané z odpovědí a nakonec spoléháme na kontext k doplnění mezer.<\/P>

<\/P>

Kromě toho, jak geometrie spolu souvisejí, existují dokonce různé kontexty pro pochopení struktury geometrií nebo typů geometrií.  Současný generace Esri ArcGIS

Geometry Object Model<\/A> existuje od vydání ArcMap 8.0.  Ačkoli se model časem trochu vyvíjel, není zásadně odlišný od doby svého vydání v roce 1999.  Nejsem si jistý historií Esri ST_Geometry storage type<\/A>, ale vím, že existuje alespoň od dob ArcGIS 9.x.  I při práci v kontextu Esri typů geometrií může vynechání kvalifikátorů vést k nejednoznačnosti.  Lze tvrdit, že Esri REST API představuje třetí Esri geometrický model, ale přidání tohoto modelu do tohoto blogového příspěvku nemění pozorování a závěry.<\/EM><\/P><\/P>Místo ponoření se do diagramů objektového modelu a referenční dokumentace si ukážeme příklad toho, jak se typy ArcGIS a ST_Geometry liší co do struktury.   Níže je obrázek a kód pro vytvoření čtyř geometrií v SQLite databázi pomocí Esri ST_Geometry typu.  Každá geometrie je označena jak svým ST_Geometry tak ArcGIS typem geometrie.<\/P><\/P><\/P><\/P><SPAN class="keyword token">import<\/SPAN> arcpy <SPAN class="keyword token">import<\/SPAN> os <SPAN class="keyword token">import<\/SPAN> sqlite3 <SPAN class="keyword token">from<\/SPAN> itertools <SPAN class="keyword token">import<\/SPAN> chain <SPAN class="comment token"># Definujte vzorky geometrií<\/SPAN> PCSCode <SPAN class="operator token">=<\/> <SPAN class="number token">3857<\/></SPAN> geoms <SPAN class="operator token">=<\/> [ &nbsp;&nbsp;&nbsp; (<SPAN class="string token">"poly"<\/></SPAN>, <SPAN class="string token">"polygon"<\/></SPAN>, <SPAN class="string token">"POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))"<\/></SPAN>) , &nbsp;&nbsp;&nbsp; (<SPAN class="string token">"multipoly"<\/></SPAN>, <SPAN class="string token">"multipolygon"<\/></SPAN>, (<SPAN class="string token">"MULTIPOLYGON(((15 0, 25 0, 25 10, 15 10, 15 0)),"<\/></SPAN> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="string token">"((18 13, 18 18, 24 18, 24 13, 18 13)))"<\/></SPAN>) ) , &nbsp;&nbsp;&nbsp; (<SPAN class="string token">"line"<\/></SPAN>, <SPAN class="string token">"linestring"<\/></SPAN>, <SPAN class="string token">"LINESTRING(3 13, 3 18)"<\/></SPAN>) , &nbsp;&nbsp;&nbsp; (<SPAN class="string token">"multiline"<\/></SPAN>, <SPAN class="string token">"multilinestring"<\/></SPAN>, <SPAN class="string token">"MULTILINESTRING((6 13, 9 13),(9 14, 9 17))"<\/></SPAN>) ] <SPAN class="comment token"># Definujte parametry inicializace SQLite databáze<\/></SPAN> sqlitedb <SPAN class="operator token">=<\/> r<SPAN class="string token">"D:\tmp\test.sqlite"<\/></SPAN> <SPAN class="comment token"># cesta k SQLite DB<\/></SPAN> st_geometry_dll <SPAN class="operator token">=<\/> r<SPAN class="string token">"DatabaseSupport\\SQLite\\Windows32\\stgeometry_sqlite.dll"<\/></SPAN> st_geometry_dll <SPAN class="operator token">=<\/> os.path.join(arcpy.GetInstallInfo()[<SPAN class="string token">"InstallDir"<\/>], &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; st_geometry_dll) sql_stmts <SPAN class="operator token">=<\/> [[ &nbsp;&nbsp;&nbsp; <SPAN class="string token">"SELECT load_extension('%s','SDE_SQL_funcs_init');"</SPANTOKEN> % st_geometry_dll &nbsp;&nbsp; ],[ &nbsp;&nbsp; (<spanclass=stringtoken> "CREATE TABLE %s "</spanclass=stringtoken> &nbsp; <spanclass=stringtoken> "(id integer primary key autoincrement not null, "</spanclass=stringtoken> &nbsp; <spanclass=stringtoken> "stgeom_type text, arcgis_type text);"</spanclass=stringtoken> % table) &nbsp; for table, geom_type, WKT in geoms &nbsp; ],[ &nbsp; (<spanclass=stringtoken> "SELECT AddGeometryColumn "</spanclass=stringtoken> &nbsp; <spanclass=stringtoken> "(null,'%s','geom',%s,'%s',2,'null');"</spanclass=stringtoken> % (table, PCSCode, geom_type) &nbsp; for table, geom_type, WKT in geoms &nbsp; ],[ &nbsp; <SPAN class="punctuation token">(<\/SPAN><SPAN class="string token">"INSERT INTO %s "<\/SPAN> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="string token">"(geom) VALUES (st_geometry('%s', %s));"<\/SPAN> <SPAN class="operator token">%<\/SPAN> <SPAN class="punctuation token">(<\/SPAN>table<SPAN class="punctuation token">,<\/SPAN> WKT<SPAN class="punctuation token">,<\/SPAN> PCSCode<SPAN class="punctuation token">)<\/SPAN><SPAN class="punctuation token">)<\/SPAN> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">pro<\/SPAN> table<SPAN class="punctuation token">,<\/SPAN> geom_type<SPAN class="punctuation token">,<\/SPAN> WKT <SPAN class="keyword token">v<\/SPAN> geoms &nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">]<\/SPAN><SPAN class="punctuation token">,<\/SPAN><SPAN class="punctuation token">[<\/SPAN> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">(<\/SPAN><SPAN class="string token">"UPDATE %s "<\/SPAN> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="string token">"SET stgeom_type = st_geometrytype(geom);"<\/SPAN> <SPAN class="operator token">%<\/SPAN> table<SPAN class="punctuation token">)<\/SPAN> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">pro<\/SPAN> table<SPAN class="punctuation token">,<\/SPAN> geom_type<SPAN class="punctuation token">,<\/SPAN> WKT <SPAN class="keyword token">v<\/SPAN> geoms &nbsp;&nbsp;&nbsp; <SPAN class="punctuation token">]<\/SPAN> <SPAN class="punctuation token">]<\/SPAN> <SPAN class="comment token"># Vytvoření SQLite databáze a nastavení tabulek<\/SPAN> arcpy<SPAN class="punctuation token">.<\/SPAN>CreateSQLiteDatabase_management<SPAN class="punctuation token">(<\/SPAN>sqlitedb<SPAN class="punctuation token">,<\/SPAN> <SPAN class="string token">"ST_GEOMETRY"<\/SPAN><SPAN class="punctuation token">)<\/SPAN> conn <SPAN class="operator token">=<\/SPAN> sqlite3<SPAN class="punctuation token">.<\/SPA N>Connection<SP ANclass= "punctu ation toke n ">(</SP AN >sqlitedb<S PANclass= "p un c t u a t i o n to k e n ">)</SP AN > conn<S PANclass= "p u n c t u a t i o n to k e n ">.< / SPAN >enable_load_extension<S PANclass= "p u n c t u a t i o n to k e n ">(</SP AN >True<S PANclass= "p u n c t u a t i o n to k e n ">)</SP AN > cur <SP ANclass= "p u n c t u a t i o n to k e n ">=< / SPAN > conn<S PANclass= "p u n c t u a t i o n to k e n ">.< / SPAN >cursor<S PANclass= "p u n c t u a t i o n to k e n ">(</SP AN >)<SP ANclass= "p u n c t u a t i o n to k e n "> </SP AN > <K EYWORDtoken for </KEYWORDtoken > sql <KEYWORDtoken in </KEYWORDtoken > chain<S PANclass= "p u n c t u a t i o n to k e n ">(</SP AN >*sql_stmts<S PANclass= "p u n c t u a t i o n to k e n ">)< / SPAN >:<SP ANclass= "p u n c t u a t i o n to k e n "> &nbsp; & nbsp ; cur<S PANclass= "p u n c t u a t i o n to k e n ">.< / SPAN >execute<S PANclass= "p u n c t u a t i o n to k e n ">(</SP AN >sql<S PANclass= "p u n c t u a t i o n to k e n ">)< / SPAN > conn<S PANclass= "p u n c t u a t i o n to k e n ">.< / SPAN >commit<S PANclass= "p u n c t u a t i o n to k e n ">(</SP AN >)<SP ANclass= "p u n c t u a t i o n to k e n "> </SP AN > <K EYWORDtoken del </KEYWORDtoken > cur<S PANclass= "p u n c t u a t i o n to k e n ">,</SP AN > conn <SP ANclass= "c omment tok en "># Naplnění typu geometrie ArcGISu</SP AN > <K EYWORDtoken for </KEYWORDtoken > table<S PANclass= "p un c tu ation tok en ">,</SP AN > geom_type<S PANclass= "p un c tu ation tok en ">,</SP AN > WKT <KEYWORDtoken in </KEYWORDtoken > geoms:<SP ANclass= "p un c tu ation tok en "> &nbsp; & nbsp ; <KEYWORDtoken with </KEYWORDtoken > arcpy<S PANclass= "p un c tu ation tok en ">.< / SPAN >da<S PANclass= "p un c tu ation tok en ">.< / SPAN >UpdateCursor<S PANclass= "p un c tu ation tok en ">(</SP AN > &nbsp; & nbsp ; & nbsp ; & nbsp ; os<S PANclass= "p un c tu ation tok en ">.< / SPAN >path<S PANclass= "p un c tu ation tok en ">.< / SPAN >join<S PANclass= "p un c tu ation tok en ">(</SP AN >sqlitedb<S PANclass= "p un c tu ation tok en ">,</SP AN > table<S PANclass= "p un c tu ation tok en ">)< / SPAN ><SP ANclass= "p un c tu ation tok en ">,</SP AN > &nbsp; & nbsp ; & nbsp ; & nbsp ; ["shape@", "arcgis_type"] &nbsp; & nbsp ; ) <KEYWORDtoken as </KEYWORDtoken > cur:<SP ANclass= "p un c tu ation tok en "> &nbsp; & nbsp ; & nbsp ; & nbsp ; <KEYWORDtoken for </KEYWORDtoken > shape<S PANclass= "p un c tu ation tok en ">,</SP AN > _ <KEYWORDtoken in </KEYWORDtoken > cur:<SP ANclass= "p un c tu ation tok en "> &nbsp; & nbsp ; & nbsp ; & nbsp ; & nbsp ; & nbsp ; cur<S PANclass= "p un c tu ation tok en ">.< / SPAN >updateRow<S PANclass= "p un c tu ation tok en ">(</SP AN >[shape, shape.type]<SP ANclass= "p un c tu ation tok en ">)< / SPAN > <K EYWORDtoken del </KEYWORDtoken > cur‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SP ANclass= "line-numbers-rows"><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍< / SPAN ><SP AN>‍</span>

Několik pozorování:

• Systém ArcGIS má jediný typ geometrie polyline, ST_Geometry má typy ST_LINESTRINGa ST_MULTILINESTRING.
• Systém ArcGIS má jediný typ geometrie polygon, ST_Geometry má typy ST_POLYGONa ST_MULTIPOLYGON.
• Zatímco ArcGIS má typy pointa multipoint, při práci s typy polylinea polygon, předpona multi- není použita.

Když typy ST_Geometry vypadají povědomě těm, kdo pracují s otevřenými standardy, není to náhoda. Dokumentace na stránce Jak je implementováno ST_Geometry? — Nápověda | ArcGIS Desktop, uvádí na několika místech, že ST_Geometry „je vysoce výkonný typ úložiště zahrnující ISO- a OGC-kompatibilní SQL přístup ke prostorovým datům.“ Existuje několik ISO- a OGC standardů týkajících se geometrií a část věty „SQL přístup ke prostorovým datům“ je ve skutečnosti velmi důležitá jako upřesnění.

Soulad s ISO a OGC je poměrně široký napříč různými modely geometrií nebo typy úložišť, např. Microsoft STGeometry, Oracle SDO_GEOM, PostGIS/PostgreSQL ST_Geometry a další. Model objektů geometrie OGC/OpenGIS je popsán v dokumentu OpenGIS4 Implementační standard pro geografické informace - Jednoduchý přístup k prvkům - Část 1: Společná architektura. Místo vložení diagramu hierarchie tříd geometrie nebo výčtu všech typů geometrií sdílím zde jsou typy LineString a MultiLineString stejně jako typy Polygon a MultiPolygon.<\/P>

<\/P>

To vše vede k jednomu bodu, a ano, mám bod, že ArcGIS polyline a polygon typy jsou v podstatě multi-typy, jen s chybějícím prefixem multi-. Ještě horší je, že většina dokumentace ArcGIS uživatele ještě více abstrahuje od struktury ArcGIS Geometry Object Model. Místo aby polyliny měly cesty a polygony měly prstence, což mají, všechno má jen "parts." Dotýkám se tohoto problému "parts" v /blogs/tilting/2016/02/20/the-single-multipart-polygons-with-interior-boundaries:<\/P>

Problém s třídami ArcPy Geometry, alespoň s jednou z nich, je ten, že Esri nahradilo více specifických geometrických komponent jedinou generickou, "part." Abstrahováním geometrií a prstenců pomocí parts nejenže odklonili od geospatialních standardů a norem, ale vytvořili i cirkusovou multi-part single-part polygon.<\/SPAN><\/BLOCKQUOTE><\/P>Ať už chybí prefixy nebo je všudypřítomný part, pochopení nuancí implementace a dokumentace ArcGIS Geometry Object Model může pomoci pochopit, jak software ArcGIS včetně ArcPy interaguje s jinými typy geometrií.
<\/SPAN><\/P><\/BODY><\/HTML>