我尝试在现有的 shapefile 中进行扩展(如果没有,我创建它然后扩展它)。我在扩展和附加方面遇到了问题,比如知道我是否应该做“data.append(list(pointlist))”,我总是得到一个
AttributeError:“多边形”对象没有属性“几何”
if os.path.isfile([name ".shp"])==False:
save_shp(name,poly)
else:
open_shp(name,poly)
if os.path.isfile([name ".shp"])==False:
save_shp(name,pointlist,area)
else:
open_shp(name,pointlist,area)
#def open_shp(name,data1,data2=None):
# if data2 is None:
# nc = gpd.read_file([name +".shp"])
# data=nc['geometry'].append(data1)
# shp = gpd.GeoDataFrame(geometry=data)
# shp.set_crs(epsg=4326, inplace=True)
# shp.to_file([name +".shp"])
# else:
# nc = gpd.read_file([name +".shp"])
# data_output=nc['area'].append(data2)
# data_geometry = nc['geometry'].append(data1)
# df = pd.DataFrame({'area':data_output})
# df["Coordinates"] = data_geometry
# shp = gpd.GeoDataFrame(geometry="Coordinates")
# shp.set_crs(epsg=4326, inplace=True)
# shp.to_file([name +".shp"])
# return
#def save_shp(name,data1,data2=None):
# if data2 is None:
# shp = gpd.GeoDataFrame(geometry=data1)
# shp.set_crs(epsg=4326, inplace=True)
# shp.to_file([name +".shp"])
# else:
# df = pd.DataFrame({'area':data2})
# df["Coordinates"] = data1
# shp = gpd.GeoDataFrame(df,geometry="Coordinates")
# shp.to_file([name +".shp"])
# return