def data_manipulation_sh (network):
from shapely.geometry import Point, LineString, MultiLineString
from geoalchemy2.shape import from_shape, to_shape
#add connection from Luebeck to Siems
new_bus = str(int(network.buses.index.max())+1)
new_trafo = str(int(network.transformers.index.max())+1)
new_line = str(int(network.lines.index.max())+1)
network.add("Bus", new_bus,carrier='AC', v_nom=220, x=10.760835, y=53.909745)
network.add("Transformer", new_trafo, bus0="25536", bus1=new_bus, x=1.29960, tap_ratio=1, s_nom=1600)
network.add("Line",new_line, bus0="26387",bus1=new_bus, x=0.0001, s_nom=1600)
network.lines.loc[new_line,'cables']=3.0
#bus geom
point_bus1 = Point(10.760835,53.909745)
network.buses.set_value(new_bus, 'geom', from_shape(point_bus1, 4326))
#line geom/topo
network.lines.set_value(new_line, 'geom', from_shape(MultiLineString([LineString([to_shape(network.buses.geom['26387']),point_bus1])]),4326))
network.lines.set_value(new_line, 'topo', from_shape(LineString([to_shape(network.buses.geom['26387']),point_bus1]),4326))
#trafo geom/topo
network.transformers.set_value(new_trafo, 'geom', from_shape(MultiLineString([LineString([to_shape(network.buses.geom['25536']),point_bus1])]),4326))
network.transformers.set_value(new_trafo, 'topo', from_shape(LineString([to_shape(network.buses.geom['25536']),point_bus1]),4326))
return
评论列表
文章目录