Eye2sky network

def get_eye2sky(station, start, end, file_type='cleaned', url=EYE2SKY_BASE_URL):
    """
    Retrieve ground-measured meteorological and solar irradiance data from the Eye2sky network.

    The Eye2Sky network [1]_ consists of 29 all-sky imager stations, of which 12 are equipped 
    with meteorological sensors in the northwest of Germany. Instrumentation at stations varies, 
    while all measure global horizontal irradiance, tilted irradiance, as well as temperature and humidity.
    Several stations measure diffuse horizontal and direct normal irradiance with a rotating shadowband irradiometer.
    Two Tier 1 stations have solar trackers measuring global, direct, and diffuse irradiance with thermopile sensors.

    Data is provided in two ways:
    1. Raw data, including quality flags from an in-house quality control procedure
    2. Cleaned data where data failing the quality control procedure has been removed

    Additionally, clear sky irradiance data and relative solar positions used for the QC are provided.

    Temporal resolution is 1 minute.
    
    Data is available from Zenodo [2]_.

    Parameters
    ----------
    station: str or list
        Station ID (5-digit ID). All stations can be requested by specifying
        station='all'.
    start: datetime-like
        First day of the requested period
    end: datetime-like
        Last day of the requested period
    file_type: str, default : 'cleaned'
        Data version: "flagged" (L0 - raw + flags) or "cleaned" (L1 - based on QC)

    Returns
    -------
    data: xarray Dataset
        Dataset of timeseries data from the Eye2Sky measurement network.

    Warns
    -----
    UserWarning
        If one or more requested files are missing, a UserWarning is returned.

    Examples
    --------
    >>> # Retrieve 2 weeks of data from Eye2Sky station OLWIN
    >>> data = get_eye2sky('OLWIN','2022-05-21','2022-07-21', type="cleaned")  # doctest: +SKIP

    References
    ----------
    .. [1] `Schmidt, Thomas; Stührenberg, Jonas; Blum, Niklas; Lezaca, Jorge; Hammer, Annette; Vogt, Thomas (2022)
            A network of all sky imagers (ASI) enabling accurate and high-resolution very short-term forecasts of solar irradiance**. 
            In: Wind & Solar Integration Workshop. IET Digital Library. 21st Wind & Solar Integration Workshop, 12.14. Okt. 2022, The Hague, Netherlands. 
            DOI: 
            <https://doi.org/10.1049/icp.2022.2778/>`_ 
    .. [2] `Zenodo
       <https://zenodo.org/records/12804613>`_ 
    """
    zipped_archive_url = os.path.join(EYE2SKY_BASE_URL, f"{station}.zip")

    response = requests.get(zipped_archive_url)

    zip_data = io.BytesIO(response.content)

    with zipfile.ZipFile(zip_data, "r") as zip_file:
        # List all files in the ZIP archive
        # print(zip_file.namelist())

        filename = f"data/{station}.{file_type}.nc"

        # Open a specific file inside the ZIP without extracting
        with zip_file.open(filename) as nc_file:
            nc_data = io.BytesIO(nc_file.read())  # Load NetCDF file into memory
            ds = xr.load_dataset(nc_data).sel(time=slice(start, end))

    return ds