ArcMap中使用ArcPy实现Geometry与WKT的相互转换

最新推荐文章于 2025-05-28 20:46:52 发布

原创最新推荐文章于 2025-05-28 20:46:52 发布 · 2.4w 阅读

20 ·

CC 4.0 BY-SA版权

用户典型问题同时被 3 个专栏收录

99 篇文章

订阅专栏

ArcGIS Desktop

57 篇文章

订阅专栏

ArcGIS10.2

49 篇文章

订阅专栏

本文介绍了如何在ArcGIS中使用Geometry与WKT之间的相互转换，包括将Geometry转换为WKT格式，以及如何通过WKT转换回Geometry。详细步骤包括在属性表中使用Field Calculator进行转换，以及如何利用ArcPy类实现转换操作。此外，还展示了如何将WKT用于查询和更新数据的方法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

在Web GIS迅猛发展的今天，使用浏览器来进行交互以其方便性、快捷性被广大用户所推崇，那么在传输格式方面，都已比较简单的JSON或者WKT来解决网络带宽带来的数据压力。

在ArcGIS10.2版本，已经提供了Geometry与JSON的相互转换，而且在现实业务中，Geometry与WKT的相互转换也比较多见。例如在ArcSDE环境中，用户可以使用Esri或者原生提供的函数来进行转换。

SQL> select sde.st_astext(shape) from city where objectid<8;

SDE.ST_ASTEXT(SHAPE)
---------------------------------------------------------------------

POINT Z ( 116.38094000 39.92361500 0.00000000)
POINT Z ( 117.20350000 39.13111900 0.00000000)
POINT Z ( 114.48978000 38.04512800 0.00000000)
POINT Z ( 112.56935000 37.87111300 0.00000000)
POINT Z ( 111.66330000 40.82094200 0.00000000)
POINT Z ( 123.41168000 41.79661600 0.00000000)
POINT Z ( 125.31543000 43.89256300 0.00000000)

已选择7行。

如果数据格式并非存储在SDE中，而是存储在文件格式中，怎么获得相关的WKT呢？

----------------------------------------------------------------------------------

建议看到转载，请直接访问正版链接获得最新的ArcGIS技术文章

Blog: http://blog.youkuaiyun.com/linghe301

----------------------------------------------------------------------------------

一：Geometry2WKT

1：比如下面存储在FGDB的三个多边形

2：创建一个TEXT类型的字段来存储WKT串，注意这个字段的长度，如果你的要素包含点比较多的话，长度可能要求比较高。我创建的字段名称为wkt，字段长度为4000.

----------------------------------------------------------------------------------

建议看到转载，请直接访问正版链接获得最新的ArcGIS技术文章

Blog: http://blog.youkuaiyun.com/linghe301

----------------------------------------------------------------------------------

3：创建好字段之后，打开要素类的属性表，然后选择wkt字段，右键选择Field Calculator选项

4：打开字段计算器，选择Python，输入wkt=！SHAPE.wkt！

注意：!SHAPE!是Python的固有写法，而且ArcPy也提供了其他函数，比如json、wkb、area等

http://resources.arcgis.com/en/help/main/10.2/index.html#/Geometry/018z00000070000000/

----------------------------------------------------------------------------------

建议看到转载，请直接访问正版链接获得最新的ArcGIS技术文章

Blog: http://blog.youkuaiyun.com/linghe301

----------------------------------------------------------------------------------

5：执行之后，我们可以看到wkt字段已经存储了这些Geometry转换后的wkt串文本信息。而且在转换过程中，并不需要ArcMap开启编辑。

用途：其实使用这种方法不仅仅可以获得相关图形的WKT，也可以变相的为文件存储的数据保存历史。

意思很明白，在编辑数据之前，我可以把原来的数据以WKT存储在某个字段中，如果编辑后发现编辑的有问题，可以使用接下来介绍的WKT2Geometry方法把图形还原回来。

----------------------------------------------------------------------------------

建议看到转载，请直接访问正版链接获得最新的ArcGIS技术文章

Blog: http://blog.youkuaiyun.com/linghe301

----------------------------------------------------------------------------------

二：WKT2Geometry

在进行ST_Geometry的使用或者测试过程中经常会碰到这种情况

select ssdl, t.shape.minx x, t.shape.miny y
  from GX_FM_PT t
 where sde.st_within(t.shape,sde.st_polyfromtext('polygon((119.91009745117188 31.763774721679688,119.91198954101563 31.763774721679688,119.91198954101563 31.765422670898438,119.91009745117188 31.763774721679688,119.91009745117188 31.763774721679688))',9)) = 1

用户使用了一个WKT串来生成一个Geometry对象来进行空间查询，首先我需要验证这个WKT是什么样子的，是否正确，是否可能有查询结果，我就需要将WKT串转换成Geometry看一下，那么该功能就需要用户掌握了。

例如上面的截图，我希望将ObjectID=2的三角形更新为ObjectID=1的四边形，那我获得了ObjectID=1的WKT串就可以实现上述操作。

（该例子仅仅说明转换方法）

----------------------------------------------------------------------------------

建议看到转载，请直接访问正版链接获得最新的ArcGIS技术文章

Blog: http://blog.youkuaiyun.com/linghe301

----------------------------------------------------------------------------------

1：打开属性表，选中ObjectID=2的要素，然后选择SHAPE字段右键