data.datasets.power_plants.wind_offshore.map_ONEP_areas() - Code Metrics - Inspection of "Features/#960 add vom for ror and reservoir" - openego/eGon-data - Measure and Improve Code Quality continuously with Scrutinizer

Passed

Pull Request — dev (#1255)

unknown

created 2025-06-11 08:09 UTC

map_ONEP_areas() A

↳ Parent: data.datasets.power_plants.wind_offshore

Complexity

Conditions

Size

Total Lines	32
Code Lines	32

Duplication

Lines	0
Ratio	0 %

Importance

Changes

Metric	Value
eloc	32
dl	0
loc	32
rs	9.112
c	0
b	0
f	0
cc	1
nop	0

from pathlib import Path
import logging

from shapely.geometry import Point
import geopandas as gpd
import pandas as pd

from egon.data import db
import egon.data.config


def map_id_bus(scenario):
    # Import manually generated list of wind offshore farms with their
    # connection points (OSM_id)
    osm_year = egon.data.config.datasets()["openstreetmap"]["original_data"][
        "source"
    ]["url"]

    if scenario in ["eGon2035", "eGon100RE"]:
        id_bus = {
            "Büttel": "136034396",
            "Suchraum Gemeinden Ibbenbüren/Mettingen/Westerkappeln": "114319248",
            "Suchraum Zensenbusch": "76185022",
            "Rommerskirchen": "24839976",
            "Oberzier": "26593929",
            "Garrel/Ost": "23837631",
            "Diele": "177829920",
            "Dörpen/West": "142487746",
            "Emden/Borßum": "34835258",
            "Hagermarsch": "79316833",
            "Hanekenfähr": "61918154",
            "Inhausen": "29420322",
            "Unterweser": "32076853",
            "Wehrendorf": "33411203",
            "Rastede": "23837631",
            "Bentwisch": "32063539",
            "Lubmin": "460134233",
            "Suchraum Gemeinde Papendorf": "32063539",
            "Suchraum Gemeinden Brünzow/Kemnitz": "460134233",
            "inhausen": "29420322",
            "Cloppenburg": "50643382",
        }
        if "200101" in osm_year:
            id_bus2 = {
                "Heide/West": "289836713",
                "Emden/Ost": "177829920",
                "Wilhelmshaven 2": "23837631",
            }
        elif ("220101" in osm_year) | ("240101" in osm_year):
            id_bus2 = {
                "Heide/West": "603661085",
                "Emden/Ost": "34835258",
                "Wilhelmshaven 2": "23744346",
            }
        else:
            raise Exception(
                """The OSM year used is not yet compatible with
                            this function"""
            )
        id_bus = {**id_bus, **id_bus2}

    elif "status" in scenario:
        year = int(scenario[-4:])

        id_bus = {
            "UW Inhausen": "29420322",
            "UW Bentwisch": "32063539",
            "UW Emden / Borssum": "34835258",
            "UW Emden Ost": "34835258",
            "UW Cloppenburg": "50643382",
            "UW Hagermarsch": "79316833",
            "UW Büttel": "136034396",
            "UW Dörpen West": "142487746",
            "UW Diele": "177829920",
            "UW Lubmin": "460134233",
        }

        if year >= 2023:
            # No update needed as no new stations used for offshore wind
            # between 2019 and 2023
            pass

        # TODO: If necessary add new stations when generating status quo > 2023

    else:
        id_bus = {}

    return id_bus


def assign_ONEP_areas():
    return {
        "Büttel": "NOR-4-1",
        "Heide/West": "NOR-10-2",
        "Suchraum Gemeinden Ibbenbüren/Mettingen/Westerkappeln": "NOR-9-2",
        "Suchraum Zensenbusch": "NOR-7-1",
        "Rommerskirchen": "NOR-7-1",
        "Oberzier": "NOR-7-1",
        "Garrel/Ost": "NOR-7-1",
        "Diele": "NOR-6-1",
        "Dörpen/West": "NOR-2-2",
        "Emden/Borßum": "NOR-0-1",
        "Emden/Ost": "NOR-3-3",
        "Hagermarsch": "NOR-2-1",
        "Hanekenfähr": "NOR-6-3",
        "Inhausen": "NOR-0-2",
        "Unterweser": "NOR-9-1",
        "Wehrendorf": "NOR-7-1",
        "Wilhelmshaven 2": "NOR-11-1",
        "Rastede": "NOR-7-1",
        "Bentwisch": "OST-3-1",
        "Lubmin": "OST-1-1",
        "Suchraum Gemeinde Papendorf": "OST-7-1",
        "Suchraum Gemeinden Brünzow/Kemnitz": "OST-1-4",
        "inhausen": "NOR-0-2",
        "Cloppenburg": "NOR-4-1",
    }


def map_ONEP_areas():
    return {
        "NOR-0-1": Point(6.5, 53.6),
        "NOR-0-2": Point(8.07, 53.76),
        "NOR-1": Point(6.21, 54.06),
        "NOR-1-1": Point(6.21, 54.06),
        "NOR-2-1": Point(6.54, 53.99),
        "NOR-2-2": Point(6.54, 53.99),
        "NOR-2-3": Point(6.54, 53.99),
        "NOR-3-1": Point(6.95, 54.02),
        "NOR-3-3": Point(6.95, 54.02),
        "NOR-4-1": Point(7.70, 54.44),
        "NOR-4-2": Point(7.70, 54.44),
        "NOR-5-1": Point(7.21, 55.14),
        "NOR-6-1": Point(5.92, 54.30),
        "NOR-6-2": Point(5.92, 54.30),
        "NOR-6-3": Point(5.92, 54.30),
        "NOR-7": Point(6.22, 54.32),
        "NOR-7-1": Point(6.22, 54.32),
        "NOR-8-1": Point(6.35, 54.48),
        "NOR-9-1": Point(5.75, 54.5),
        "NOR-9-2": Point(5.75, 54.5),
        "NOR-10-2": Point(6, 54.75),
        "NOR-11-1": Point(6.5, 54.75),
        "OST-1-1": Point(14.09, 54.82),
        "OST-1-2": Point(14.09, 54.82),
        "OST-1-3": Point(14.09, 54.82),
        "OST-1-4": Point(14.09, 54.82),
        "OST-2": Point(13.86, 54.83),
        "OST-3-1": Point(13.16, 54.98),
        "OST-3-2": Point(13.16, 54.98),
        "OST-7-1": Point(12.25, 54.5),
    }


def insert():
    """
    Include the offshore wind parks in egon-data.

    Parameters
    ----------
    *No parameters required
    """
    # Read file with all required input/output tables' names
    cfg = egon.data.config.datasets()["power_plants"]

    scenarios = egon.data.config.settings()["egon-data"]["--scenarios"]

    for scenario in scenarios:
        # Delete previous generators
        db.execute_sql(
            f"""
            DELETE FROM {cfg['target']['schema']}.{cfg['target']['table']}
            WHERE carrier = 'wind_offshore'
            AND scenario = '{scenario}'
            """
        )

        # load file
        if scenario == "eGon2035":
            offshore_path = (
                Path(".")
                / "data_bundle_egon_data"
                / "nep2035_version2021"
                / cfg["sources"]["nep_2035"]
            )

            offshore = pd.read_excel(
                offshore_path,
                sheet_name="WInd_Offshore_NEP",
                usecols=[
                    "Netzverknuepfungspunkt",
                    "Spannungsebene in kV",
                    "C 2035",
                ],
            )
            offshore.dropna(subset=["Netzverknuepfungspunkt"], inplace=True)
            offshore.rename(columns={"C 2035": "el_capacity"}, inplace=True)
            offshore = offshore[offshore["el_capacity"] > 0]

        elif scenario == "eGon100RE":
            offshore_path = (
                Path(".")
                / "data_bundle_egon_data"
                / "nep2035_version2021"
                / cfg["sources"]["nep_2035"]
            )

            offshore = pd.read_excel(
                offshore_path,
                sheet_name="WInd_Offshore_NEP",
                usecols=[
                    "Netzverknuepfungspunkt",
                    "Spannungsebene in kV",
                    "B 2040 ",
                ],
            )
            offshore.dropna(subset=["Netzverknuepfungspunkt"], inplace=True)
            offshore.rename(columns={"B 2040 ": "el_capacity"}, inplace=True)
            offshore = offshore[offshore["el_capacity"] > 0]

        elif "status" in scenario:
            year = int(scenario[-4:])

            offshore_path = (
                Path(".")
                / "data_bundle_egon_data"
                / "wind_offshore_status2019"
                / cfg["sources"]["wind_offshore_status2019"]
            )
            offshore = pd.read_excel(
                offshore_path,
                sheet_name="wind_offshore",
                usecols=[
                    "Name ONEP/NEP",
                    "NVP",
                    "Spannung [kV]",
                    "Inbetriebnahme",
                    "Kapazität Gesamtsystem [MW]",
                ],
            )
            offshore.dropna(subset=["Name ONEP/NEP"], inplace=True)
            offshore.rename(
                columns={
                    "NVP": "Netzverknuepfungspunkt",
                    "Spannung [kV]": "Spannungsebene in kV",
                    "Kapazität Gesamtsystem [MW]": "el_capacity",
                },
                inplace=True,
            )
            offshore = offshore[offshore["Inbetriebnahme"] <= year]

        else:
            raise ValueError(f"{scenario=} is not valid.")

        id_bus = map_id_bus(scenario)

        # Match wind offshore table with the corresponding OSM_id
        offshore["osm_id"] = offshore["Netzverknuepfungspunkt"].map(id_bus)

        buses = db.select_geodataframe(
            f"""
                SELECT bus_i as bus_id, base_kv, geom as point, CAST(osm_substation_id AS text)
                as osm_id FROM {cfg["sources"]["buses_data"]}
                """,
            epsg=4326,
            geom_col="point",
        )

        # Drop NANs in column osm_id
        buses.dropna(subset=["osm_id"], inplace=True)

        # Create columns for bus_id and geometry in the offshore df
        offshore["bus_id"] = pd.NA
        offshore["geom"] = Point(0, 0)

        # Match bus_id
        for index, wind_park in offshore.iterrows():
            if not buses[
                (buses["osm_id"] == wind_park["osm_id"])
                & (buses["base_kv"] == wind_park["Spannungsebene in kV"])
            ].empty:
                bus_ind = buses[buses["osm_id"] == wind_park["osm_id"]].index[
                    0
                ]
                offshore.at[index, "bus_id"] = buses.at[bus_ind, "bus_id"]
            else:
                print(f'Wind offshore farm not found: {wind_park["osm_id"]}')

        offshore.dropna(subset=["bus_id"], inplace=True)

        # Overwrite geom for status2019 parks
        if scenario in ["eGon2035", "eGon100RE"]:
            offshore["Name ONEP/NEP"] = offshore["Netzverknuepfungspunkt"].map(
                assign_ONEP_areas()
            )

        offshore["geom"] = offshore["Name ONEP/NEP"].map(map_ONEP_areas())
        offshore["weather_cell_id"] = pd.NA

        offshore.drop(["Name ONEP/NEP"], axis=1, inplace=True)

        if "status" in scenario:
            offshore.drop(["Inbetriebnahme"], axis=1, inplace=True)

        # Scale capacities for eGon100RE
        if scenario == "eGon100RE":
            # Import capacity targets for wind_offshore per scenario
            cap_100RE = db.select_dataframe(
                f"""
                    SELECT SUM(capacity)
                    FROM {cfg["sources"]["capacities"]}
                    WHERE scenario_name = 'eGon100RE' AND
                    carrier = 'wind_offshore'
                    """
            ).iloc[0, 0]

            # Scale capacities to match  target
            scale_factor = cap_100RE / offshore.el_capacity.sum()
            offshore["el_capacity"] *= scale_factor

        # Assign voltage levels to wind offshore parks
        offshore["voltage_level"] = 0
        offshore.loc[
            offshore[offshore["Spannungsebene in kV"] == 110].index,
            "voltage_level",
        ] = 3
        offshore.loc[
            offshore[offshore["Spannungsebene in kV"] > 110].index,
            "voltage_level",
        ] = 1

        # Delete unnecessary columns
        offshore.drop(
            [
                "Netzverknuepfungspunkt",
                "Spannungsebene in kV",
                "osm_id",
            ],
            axis=1,
            inplace=True,
        )

        # Set static columns
        offshore["carrier"] = "wind_offshore"
        offshore["scenario"] = scenario

        offshore = gpd.GeoDataFrame(offshore, geometry="geom", crs=4326)

        # Look for the maximum id in the table egon_power_plants
        next_id = db.select_dataframe(
            "SELECT MAX(id) FROM "
            + cfg["target"]["schema"]
            + "."
            + cfg["target"]["table"]
        ).iloc[0, 0]

        if next_id:
            next_id += 1
        else:
            next_id = 1

        # Reset index
        offshore.index = pd.RangeIndex(
            start=next_id, stop=next_id + len(offshore), name="id"
        )

        # Insert into database
        offshore.reset_index().to_postgis(
            cfg["target"]["table"],
            schema=cfg["target"]["schema"],
            con=db.engine(),
            if_exists="append",
        )

        logging.info(
            f"""
              {len(offshore)} wind_offshore generators with a total installed capacity of
              {offshore['el_capacity'].sum()}MW were inserted into the db
              """
        )


1			from pathlib import Path
2			import logging
3
4			from shapely.geometry import Point
5			import geopandas as gpd
6			import pandas as pd
7
8			from egon.data import db
9			import egon.data.config
10
11
12			def map_id_bus(scenario):
13			# Import manually generated list of wind offshore farms with their
14			# connection points (OSM_id)
15			osm_year = egon.data.config.datasets()["openstreetmap"]["original_data"][
16			"source"
17			]["url"]
18
19			if scenario in ["eGon2035", "eGon100RE"]:
20			id_bus = {
21			"Büttel": "136034396",
22			"Suchraum Gemeinden Ibbenbüren/Mettingen/Westerkappeln": "114319248",
23			"Suchraum Zensenbusch": "76185022",
24			"Rommerskirchen": "24839976",
25			"Oberzier": "26593929",
26			"Garrel/Ost": "23837631",
27			"Diele": "177829920",
28			"Dörpen/West": "142487746",
29			"Emden/Borßum": "34835258",
30			"Hagermarsch": "79316833",
31			"Hanekenfähr": "61918154",
32			"Inhausen": "29420322",
33			"Unterweser": "32076853",
34			"Wehrendorf": "33411203",
35			"Rastede": "23837631",
36			"Bentwisch": "32063539",
37			"Lubmin": "460134233",
38			"Suchraum Gemeinde Papendorf": "32063539",
39			"Suchraum Gemeinden Brünzow/Kemnitz": "460134233",
40			"inhausen": "29420322",
41			"Cloppenburg": "50643382",
42			}
43			if "200101" in osm_year:
44			id_bus2 = {
45			"Heide/West": "289836713",
46			"Emden/Ost": "177829920",
47			"Wilhelmshaven 2": "23837631",
48			}
49			elif ("220101" in osm_year) \| ("240101" in osm_year):
50			id_bus2 = {
51			"Heide/West": "603661085",
52			"Emden/Ost": "34835258",
53			"Wilhelmshaven 2": "23744346",
54			}
55			else:
56			raise Exception(
57			"""The OSM year used is not yet compatible with
58			this function"""
59			)
60			id_bus = {id_bus, id_bus2}
61
62			elif "status" in scenario:
63			year = int(scenario[-4:])
64
65			id_bus = {
66			"UW Inhausen": "29420322",
67			"UW Bentwisch": "32063539",
68			"UW Emden / Borssum": "34835258",
69			"UW Emden Ost": "34835258",
70			"UW Cloppenburg": "50643382",
71			"UW Hagermarsch": "79316833",
72			"UW Büttel": "136034396",
73			"UW Dörpen West": "142487746",
74			"UW Diele": "177829920",
75			"UW Lubmin": "460134233",
76			}
77
78			if year >= 2023:
79			# No update needed as no new stations used for offshore wind
80			# between 2019 and 2023
81			pass
82
83			# TODO: If necessary add new stations when generating status quo > 2023
84
85			else:
86			id_bus = {}
87
88			return id_bus
89
90
91			def assign_ONEP_areas():
92			return {
93			"Büttel": "NOR-4-1",
94			"Heide/West": "NOR-10-2",
95			"Suchraum Gemeinden Ibbenbüren/Mettingen/Westerkappeln": "NOR-9-2",
96			"Suchraum Zensenbusch": "NOR-7-1",
97			"Rommerskirchen": "NOR-7-1",
98			"Oberzier": "NOR-7-1",
99			"Garrel/Ost": "NOR-7-1",
100			"Diele": "NOR-6-1",
101			"Dörpen/West": "NOR-2-2",
102			"Emden/Borßum": "NOR-0-1",
103			"Emden/Ost": "NOR-3-3",
104			"Hagermarsch": "NOR-2-1",
105			"Hanekenfähr": "NOR-6-3",
106			"Inhausen": "NOR-0-2",
107			"Unterweser": "NOR-9-1",
108			"Wehrendorf": "NOR-7-1",
109			"Wilhelmshaven 2": "NOR-11-1",
110			"Rastede": "NOR-7-1",
111			"Bentwisch": "OST-3-1",
112			"Lubmin": "OST-1-1",
113			"Suchraum Gemeinde Papendorf": "OST-7-1",
114			"Suchraum Gemeinden Brünzow/Kemnitz": "OST-1-4",
115			"inhausen": "NOR-0-2",
116			"Cloppenburg": "NOR-4-1",
117			}
118
119
120			def map_ONEP_areas():
121			return {
122			"NOR-0-1": Point(6.5, 53.6),
123			"NOR-0-2": Point(8.07, 53.76),
124			"NOR-1": Point(6.21, 54.06),
125			"NOR-1-1": Point(6.21, 54.06),
126			"NOR-2-1": Point(6.54, 53.99),
127			"NOR-2-2": Point(6.54, 53.99),
128			"NOR-2-3": Point(6.54, 53.99),
129			"NOR-3-1": Point(6.95, 54.02),
130			"NOR-3-3": Point(6.95, 54.02),
131			"NOR-4-1": Point(7.70, 54.44),
132			"NOR-4-2": Point(7.70, 54.44),
133			"NOR-5-1": Point(7.21, 55.14),
134			"NOR-6-1": Point(5.92, 54.30),
135			"NOR-6-2": Point(5.92, 54.30),
136			"NOR-6-3": Point(5.92, 54.30),
137			"NOR-7": Point(6.22, 54.32),
138			"NOR-7-1": Point(6.22, 54.32),
139			"NOR-8-1": Point(6.35, 54.48),
140			"NOR-9-1": Point(5.75, 54.5),
141			"NOR-9-2": Point(5.75, 54.5),
142			"NOR-10-2": Point(6, 54.75),
143			"NOR-11-1": Point(6.5, 54.75),
144			"OST-1-1": Point(14.09, 54.82),
145			"OST-1-2": Point(14.09, 54.82),
146			"OST-1-3": Point(14.09, 54.82),
147			"OST-1-4": Point(14.09, 54.82),
148			"OST-2": Point(13.86, 54.83),
149			"OST-3-1": Point(13.16, 54.98),
150			"OST-3-2": Point(13.16, 54.98),
151			"OST-7-1": Point(12.25, 54.5),
152			}
153
154
155			def insert():
156			"""
157			Include the offshore wind parks in egon-data.
158
159			Parameters
160			----------
161			*No parameters required
162			"""
163			# Read file with all required input/output tables' names
164			cfg = egon.data.config.datasets()["power_plants"]
165
166			scenarios = egon.data.config.settings()["egon-data"]["--scenarios"]
167
168			for scenario in scenarios:
169			# Delete previous generators
170			db.execute_sql(
171			f"""
172			DELETE FROM {cfg['target']['schema']}.{cfg['target']['table']}
173			WHERE carrier = 'wind_offshore'
174			AND scenario = '{scenario}'
175			"""
176			)
177
178			# load file
179			if scenario == "eGon2035":
180			offshore_path = (
181			Path(".")
182			/ "data_bundle_egon_data"
183			/ "nep2035_version2021"
184			/ cfg["sources"]["nep_2035"]
185			)
186
187			offshore = pd.read_excel(
188			offshore_path,
189			sheet_name="WInd_Offshore_NEP",
190			usecols=[
191			"Netzverknuepfungspunkt",
192			"Spannungsebene in kV",
193			"C 2035",
194			],
195			)
196			offshore.dropna(subset=["Netzverknuepfungspunkt"], inplace=True)
197			offshore.rename(columns={"C 2035": "el_capacity"}, inplace=True)
198			offshore = offshore[offshore["el_capacity"] > 0]
199
200			elif scenario == "eGon100RE":
201			offshore_path = (
202			Path(".")
203			/ "data_bundle_egon_data"
204			/ "nep2035_version2021"
205			/ cfg["sources"]["nep_2035"]
206			)
207
208			offshore = pd.read_excel(
209			offshore_path,
210			sheet_name="WInd_Offshore_NEP",
211			usecols=[
212			"Netzverknuepfungspunkt",
213			"Spannungsebene in kV",
214			"B 2040 ",
215			],
216			)
217			offshore.dropna(subset=["Netzverknuepfungspunkt"], inplace=True)
218			offshore.rename(columns={"B 2040 ": "el_capacity"}, inplace=True)
219			offshore = offshore[offshore["el_capacity"] > 0]
220
221			elif "status" in scenario:
222			year = int(scenario[-4:])
223
224			offshore_path = (
225			Path(".")
226			/ "data_bundle_egon_data"
227			/ "wind_offshore_status2019"
228			/ cfg["sources"]["wind_offshore_status2019"]
229			)
230			offshore = pd.read_excel(
231			offshore_path,
232			sheet_name="wind_offshore",
233			usecols=[
234			"Name ONEP/NEP",
235			"NVP",
236			"Spannung [kV]",
237			"Inbetriebnahme",
238			"Kapazität Gesamtsystem [MW]",
239			],
240			)
241			offshore.dropna(subset=["Name ONEP/NEP"], inplace=True)
242			offshore.rename(
243			columns={
244			"NVP": "Netzverknuepfungspunkt",
245			"Spannung [kV]": "Spannungsebene in kV",
246			"Kapazität Gesamtsystem [MW]": "el_capacity",
247			},
248			inplace=True,
249			)
250			offshore = offshore[offshore["Inbetriebnahme"] <= year]
251
252			else:
253			raise ValueError(f"{scenario=} is not valid.")
254
255			id_bus = map_id_bus(scenario)
256
257			# Match wind offshore table with the corresponding OSM_id
258			offshore["osm_id"] = offshore["Netzverknuepfungspunkt"].map(id_bus)
259
260			buses = db.select_geodataframe(
261			f"""
262			SELECT bus_i as bus_id, base_kv, geom as point, CAST(osm_substation_id AS text)
263			as osm_id FROM {cfg["sources"]["buses_data"]}
264			""",
265			epsg=4326,
266			geom_col="point",
267			)
268
269			# Drop NANs in column osm_id
270			buses.dropna(subset=["osm_id"], inplace=True)
271
272			# Create columns for bus_id and geometry in the offshore df
273			offshore["bus_id"] = pd.NA
274			offshore["geom"] = Point(0, 0)
275
276			# Match bus_id
277			for index, wind_park in offshore.iterrows():
278			if not buses[
279			(buses["osm_id"] == wind_park["osm_id"])
280			& (buses["base_kv"] == wind_park["Spannungsebene in kV"])
281			].empty:
282			bus_ind = buses[buses["osm_id"] == wind_park["osm_id"]].index[
283			0
284			]
285			offshore.at[index, "bus_id"] = buses.at[bus_ind, "bus_id"]
286			else:
287			print(f'Wind offshore farm not found: {wind_park["osm_id"]}')
288
289			offshore.dropna(subset=["bus_id"], inplace=True)
290
291			# Overwrite geom for status2019 parks
292			if scenario in ["eGon2035", "eGon100RE"]:
293			offshore["Name ONEP/NEP"] = offshore["Netzverknuepfungspunkt"].map(
294			assign_ONEP_areas()
295			)
296
297			offshore["geom"] = offshore["Name ONEP/NEP"].map(map_ONEP_areas())
298			offshore["weather_cell_id"] = pd.NA
299
300			offshore.drop(["Name ONEP/NEP"], axis=1, inplace=True)
301
302			if "status" in scenario:
303			offshore.drop(["Inbetriebnahme"], axis=1, inplace=True)
304
305			# Scale capacities for eGon100RE
306			if scenario == "eGon100RE":
307			# Import capacity targets for wind_offshore per scenario
308			cap_100RE = db.select_dataframe(
309			f"""
310			SELECT SUM(capacity)
311			FROM {cfg["sources"]["capacities"]}
312			WHERE scenario_name = 'eGon100RE' AND
313			carrier = 'wind_offshore'
314			"""
315			).iloc[0, 0]
316
317			# Scale capacities to match target
318			scale_factor = cap_100RE / offshore.el_capacity.sum()
319			offshore["el_capacity"] *= scale_factor
320
321			# Assign voltage levels to wind offshore parks
322			offshore["voltage_level"] = 0
323			offshore.loc[
324			offshore[offshore["Spannungsebene in kV"] == 110].index,
325			"voltage_level",
326			] = 3
327			offshore.loc[
328			offshore[offshore["Spannungsebene in kV"] > 110].index,
329			"voltage_level",
330			] = 1
331
332			# Delete unnecessary columns
333			offshore.drop(
334			[
335			"Netzverknuepfungspunkt",
336			"Spannungsebene in kV",
337			"osm_id",
338			],
339			axis=1,
340			inplace=True,
341			)
342
343			# Set static columns
344			offshore["carrier"] = "wind_offshore"
345			offshore["scenario"] = scenario
346
347			offshore = gpd.GeoDataFrame(offshore, geometry="geom", crs=4326)
348
349			# Look for the maximum id in the table egon_power_plants
350			next_id = db.select_dataframe(
351			"SELECT MAX(id) FROM "
352			+ cfg["target"]["schema"]
353			+ "."
354			+ cfg["target"]["table"]
355			).iloc[0, 0]
356
357			if next_id:
358			next_id += 1
359			else:
360			next_id = 1
361
362			# Reset index
363			offshore.index = pd.RangeIndex(
364			start=next_id, stop=next_id + len(offshore), name="id"
365			)
366
367			# Insert into database
368			offshore.reset_index().to_postgis(
369			cfg["target"]["table"],
370			schema=cfg["target"]["schema"],
371			con=db.engine(),
372			if_exists="append",
373			)
374
375			logging.info(
376			f"""
377			{len(offshore)} wind_offshore generators with a total installed capacity of
378			{offshore['el_capacity'].sum()}MW were inserted into the db
379			"""
380			)
381

openego / eGon-data

Pull Request — dev (#1255)

map_ONEP_areas() A

Complexity

Size

Duplication

Importance

Duplication Side-by-Side

Filter issues like