5ea29551b9
* add REUSE.software license information * add release note on licensing * Update doc/index.rst Co-authored-by: euronion <42553970+euronion@users.noreply.github.com> * Update doc/index.rst Co-authored-by: euronion <42553970+euronion@users.noreply.github.com> * Update doc/release_notes.rst [skip travis] Co-authored-by: euronion <42553970+euronion@users.noreply.github.com>
71 lines
2.2 KiB
Python
71 lines
2.2 KiB
Python
#!/usr/bin/env python
|
||
|
||
# SPDX-FileCopyrightText: : 2017-2020 The PyPSA-Eur Authors
|
||
#
|
||
# SPDX-License-Identifier: GPL-3.0-or-later
|
||
|
||
"""
|
||
Extracts capacities of HVDC links from `Wikipedia <https://en.wikipedia.org/wiki/List_of_HVDC_projects>`_.
|
||
|
||
Relevant Settings
|
||
-----------------
|
||
|
||
.. code:: yaml
|
||
|
||
enable:
|
||
prepare_links_p_nom:
|
||
|
||
.. seealso::
|
||
Documentation of the configuration file ``config.yaml`` at
|
||
:ref:`toplevel_cf`
|
||
|
||
Inputs
|
||
------
|
||
|
||
*None*
|
||
|
||
Outputs
|
||
-------
|
||
|
||
- ``data/links_p_nom.csv``: A plain download of https://en.wikipedia.org/wiki/List_of_HVDC_projects#Europe plus extracted coordinates.
|
||
|
||
Description
|
||
-----------
|
||
|
||
*None*
|
||
|
||
"""
|
||
|
||
import logging
|
||
logger = logging.getLogger(__name__)
|
||
from _helpers import configure_logging
|
||
|
||
import pandas as pd
|
||
|
||
if __name__ == "__main__":
|
||
if 'snakemake' not in globals():
|
||
from _helpers import mock_snakemake #rule must be enabled in config
|
||
snakemake = mock_snakemake('prepare_links_p_nom', simpl='', network='elec')
|
||
configure_logging(snakemake)
|
||
|
||
links_p_nom = pd.read_html('https://en.wikipedia.org/wiki/List_of_HVDC_projects', header=0, match="SwePol")[0]
|
||
|
||
def extract_coordinates(s):
|
||
regex = (r"(\d{1,2})°(\d{1,2})′(\d{1,2})″(N|S) "
|
||
r"(\d{1,2})°(\d{1,2})′(\d{1,2})″(E|W)")
|
||
e = s.str.extract(regex, expand=True)
|
||
lat = (e[0].astype(float) + (e[1].astype(float) + e[2].astype(float)/60.)/60.)*e[3].map({'N': +1., 'S': -1.})
|
||
lon = (e[4].astype(float) + (e[5].astype(float) + e[6].astype(float)/60.)/60.)*e[7].map({'E': +1., 'W': -1.})
|
||
return lon, lat
|
||
|
||
m_b = links_p_nom["Power (MW)"].str.contains('x').fillna(False)
|
||
def multiply(s): return s.str[0].astype(float) * s.str[1].astype(float)
|
||
|
||
links_p_nom.loc[m_b, "Power (MW)"] = links_p_nom.loc[m_b, "Power (MW)"].str.split('x').pipe(multiply)
|
||
links_p_nom["Power (MW)"] = links_p_nom["Power (MW)"].str.extract("[-/]?([\d.]+)", expand=False).astype(float)
|
||
|
||
links_p_nom['x1'], links_p_nom['y1'] = extract_coordinates(links_p_nom['Converterstation 1'])
|
||
links_p_nom['x2'], links_p_nom['y2'] = extract_coordinates(links_p_nom['Converterstation 2'])
|
||
|
||
links_p_nom.dropna(subset=['x1', 'y1', 'x2', 'y2']).to_csv(snakemake.output[0], index=False)
|