def load_pv_systems(metadata_filename: str=METADATA_FILENAME, stats_filename: str=PV_STATS_FILENAME, timeseries_filename: str=TIMESERIES_FILENAME) -> xr.Dataset:
    'Load metadata about PV systems'
    pv_metadata = pd.read_csv(metadata_filename, index_col='system_id')
    pv_stats = pd.read_csv(stats_filename, index_col='system_id', parse_dates=['actual_date_from', 'actual_date_to', 'record_efficiency_date'])
    pv_systems = pv_metadata.join(pv_stats[['actual_date_from', 'actual_date_to', 'outputs']], how='left')
    pv_systems_filtered = pv_systems.query('status_interval_minutes <= 60 and outputs > 100')
    pv_systems_filtered = pv_systems_filtered.dropna(subset=['latitude', 'longitude'])
    system_ids = _get_system_ids_dataframe_from_timeseries(timeseries_filename)
    pv_systems_filtered = pv_systems_filtered.join(system_ids, how='inner')
    pv_systems_filtered = pv_systems_filtered[['system_name', 'latitude', 'longitude']]
    ds = xr.Dataset.from_dataframe(pv_systems_filtered)
    ds = _transform_pv_systems(ds)
    return ds
Load metadata about PV systems
def load_pv_systems(metadata_filename: str=METADATA_FILENAME, stats_filename: str=PV_STATS_FILENAME, timeseries_filename: str=TIMESERIES_FILENAME) -> xr.Dataset:
    
    pv_metadata = pd.read_csv(metadata_filename, index_col='system_id')
    pv_stats = pd.read_csv(stats_filename, index_col='system_id', parse_dates=['actual_date_from', 'actual_date_to', 'record_efficiency_date'])
    pv_systems = pv_metadata.join(pv_stats[['actual_date_from', 'actual_date_to', 'outputs']], how='left')
    pv_systems_filtered = pv_systems.query('status_interval_minutes <= 60 and outputs > 100')
    pv_systems_filtered = pv_systems_filtered.dropna(subset=['latitude', 'longitude'])
    system_ids = _get_system_ids_dataframe_from_timeseries(timeseries_filename)
    pv_systems_filtered = pv_systems_filtered.join(system_ids, how='inner')
    pv_systems_filtered = pv_systems_filtered[['system_name', 'latitude', 'longitude']]
    ds = xr.Dataset.from_dataframe(pv_systems_filtered)
    ds = _transform_pv_systems(ds)
    return ds