我正在使用一个数据库，其中包含一个代表国家/地区 Admin1 边界的地理字段。我想绘制一张索马里 Admin1 边界的地图：
params = urllib.quote_plus(...)
engine = sqlalchemy.create_engine("mssql+pyodbc:///?odbc_connect=%s" % params)
query = ("SELECT SimplifiedCoords.STAsText() FROM Adm1GeoBoundaries WHERE (Adm0Code = '226')")
df = pd.read_sql(query,engine)
这将返回 18 个细分的列表，例如：
0 POLYGON ((44.661 3.2516, 44.611 3.2715, 44.410...
1 POLYGON ((49.0611 9.571, 48.561 9.541, 48.061 ...
2 POLYGON ((44.2362 10.4724, 44.2021 10.5066, 44...
我现在的问题是将这些转换为 shapefile，以便使用 arcpy 映射它们。到目前为止，我已经尝试过：
wkt = (df.loc[1]).to_string() # I only focus on one to test
g1 = loads(wkt)
但这会抛出：“ParseException：预期数字但遇到单词：'46.151 ...'”
我还尝试了以下方法：
sr = arcpy.SpatialReference(4326)
g1 = arcpy.FromWKT(wkt,sr)
但这还没有把我带到任何地方......欢迎任何想法/建议。提前致谢 ：）