地理＆＃x2B; MySQL- GIS数据 /无效的几何形状

Question

我有一个geopandas我插入MySQL数据库中的dataframe。当我尝试在MySQL表中设置几何时，我遇到了一个错误。我在Python中检查有效的几何形状，然后插入SQL。这是我的代码和示例数据：

示例geo-data： https://www.dropbox.com/s/14NI2MFPPT5DN7X/GDF%20%281%281%299999999.DL = 0

import geopandas as gpd
from shapely.geometry import Point, Polygon 
from shapely import wkt

# mysql connection
import pymysql
from sqlalchemy import create_engine
user = 'user'
pwd = 'pwd'
host =  'host1'
port = 3306
database = 'db'

engine = create_engine("mysql+pymysql://{}:{}@{}/{}".format(user,pwd,host,database))

gdf = gpd.read_file()

# To GeoPandas
gdf['geometry'] = gdf['zip_code_geom'].apply(wkt.loads)

# Fix Bad Geometries
def valid_geom(geom):
    try:
       return wkt.loads(geom)
    except:
       return np.nan


gdf['geometry'] = gdf.zip_code_geom.apply(valid_geom)

# Set Geometry
gdf = gdf.GeoDataFrame(gdf, geometry='geometry')

# MultiPolygon to Polygon
gdf = gdf.explode(column='geometry', ignore_index=True, index_parts=False)

# Insert into SQL
con = engine.connect()
gdf.to_sql(name="gdf_sql", con=con, if_exists = 'replace')
con.close()

现在，geterry列被存储为<代码>文本在mySQL中，我们需要将其转换为polygon。

我尝试了以下操作，但是继续获取​​错误代码：3037。提供给功能的无效GIS数据ST_POLYGONFROMTEXT。

UPDATE db1.gdf_sql SET geometry = ST_PolygonFromText(geometry) WHERE geometry!='' AND geometry IS NOT NULL;

ALTER TABLE db1.gdf_sql
ADD COLUMN geom_poly Polygon;

# insert column into table
UPDATE db1.gdf_sql 
SET geometry = (
    SELECT ST_GeomFromTxt(t1.geometry)
    FROM geom as t1
    WHERE t1.geometry = db1.gdf_sql.geometry
)

我需要在mySQL中保存几何>几何>，以便能够执行空间功能在我的查询中。

Answer 1

对于一个 - nan不是有效的几何形状。因此，您可以在发送到SQL之前丢弃NAN，也可以正确处理st_geomfromtxt正确处理它们。为此，您可以使用一个空的geomemycollection，例如shapely.geometry.geometrycollection（）。

第二 - 转换为任何字符串格式时，geopandas运行 to_wkt 在引擎盖下。依次调用 pygeos.pygeos.pygeos.to_wkt ，所有关键字参数（和默认值）都传递给其他功能（继承）。默认情况下，在发送到WKT（然后是SQL）之前，几何形状的精度降低到6个小数位，这可能会导致有效的形状变得无效，尤其是当您以程度谈论地理空间数据时。摘录摘自 pygeos docs ：

to_wkt （几何，roughing_precision = 6，trim = true，output_dimension = 3，old_3d = false，** kwargs）

转换为几何的众所周知的文本（WKT）表示。

参数

在

编写WKT字符串时的舍入精度。设置为-1的值以指示完整的精度。

为了解决这个问题，您可以在发送到SQL之前自己进行转换：

df = pd.DataFrame(gdf).assign(
    geometry=gdf.geometry.to_wkt(rounding_precision=-1),
)

最后，有一个“有效的几何形状”概念，与“通过wkt.loads可解行”不同，因此您的测试可能无法按预期工作。而可以使用共线内环和其他“和其他”来定义多边形无效的“形状在整齐的形状中，空间操作将根据语言和操作的不同（或可能无法按预期起作用）。我不确定MySQL的几何引擎如何处理此问题 - 它们可能只是丢失一个错误，因此这是一个可能会渗入错误的地方。您可以使用 gdf.geometry.is_valid.all（）带有gdf ​​[〜gd.geometry.is_valid]的无效。