Swmmio面文件不能正确显示

面文件不能与点线正确投影

• Swmmio version 0.6.7
• 在ArcGIS中打开，面文件没有范围显示
• 原因是当前版本的Swmmio在write_shapefile时并未给添加面文件添加geometry，这个bug在2023/7/24号下午刚提交上去
  

解决方案

• 需要在swmmio的安装目录下的utils文件中对spatial.py文件进行修改(line 218)
• if geomtype== 'polygon': w.poly([row.coords])
• 或者等Swmmio更新版本后安装更新的版本

导出时的投影问题

• Swmmio时根据Pyshp提供的api导出shp文件，可以需要根据投影文件(crs文件)自定义投影
• 而目前的GDAL中是利用数据库文件（proj.db）存储投影信息，这导致无法直接调用crs文件去定义Swmmio的输出信息
• 可以通过下述代码输出crs文件解决

from osgeo import ogr, osr

spatialRef = osr.SpatialReference()
spatialRef.ImportFromEPSG(32650)

spatialRef.MorphToESRI()
file = open('G:\\Desktop\\yourshpfile.prj', 'w')
file.write(spatialRef.ExportToWkt())
file.close() ```