我正在使用包含代表国家/地区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)
但这还没有把我带到任何地方......任何想法/建议欢迎。在此先感谢：）