Site info refactor

examples/CSP_PV_Battery_Analysis/simulation_init.py

@@ -82,8 +82,7 @@ def init_hybrid_plant(techs_in_sim: list, is_test: bool = False, ud_techs: dict
         "lon": -116.7830,
         "elev": 561,
         "tz": 1,
-        "no_wind": True
-        }
+    }
     solar_file = example_root + "02_weather_data/daggett_ca_34.865371_-116.783023_psmv3_60_tmy.csv"
     prices_file = example_root + "03_cost_load_price_data/constant_norm_prices.csv"
     desired_schedule_file = example_root + "03_cost_load_price_data/desired_schedule_normalized.csv"
@@ -98,7 +97,8 @@ def init_hybrid_plant(techs_in_sim: list, is_test: bool = False, ud_techs: dict
     site = SiteInfo(site_data, 
                     solar_resource_file=solar_file, 
                     grid_resource_file=prices_file,
-                    desired_schedule=desired_schedule
+                    desired_schedule=desired_schedule,
+                    wind=False
                     )
 
     # Load in system costs

examples/offshore-hybrid/wind-h2.py

@@ -160,14 +160,14 @@ def setup_hopp(plant_config, turbine_config, wind_resource, orbit_project, flori
     hopp_site_input_data["lat"] = plant_config["project_location"]["lat"]
     hopp_site_input_data["lon"] = plant_config["project_location"]["lon"]
     hopp_site_input_data["year"] = plant_config["wind_resource_year"]
-    hopp_site_input_data["no_solar"] = not plant_config["project_parameters"]["solar"]
+    solar = plant_config["project_parameters"]["solar"]
 
     # set desired schedule based on electrolyzer capacity
     desired_schedule = [plant_config["electrolyzer"]["rating"]]*8760
     desired_schedule = []
 
     # generate HOPP SiteInfo class instance
-    hopp_site = SiteInfo(hopp_site_input_data, hub_height=turbine_config["hub_height"], desired_schedule=desired_schedule)
+    hopp_site = SiteInfo(hopp_site_input_data, hub_height=turbine_config["hub_height"], desired_schedule=desired_schedule, solar=solar)
 
     # replace wind data with previously downloaded and adjusted wind data
     hopp_site.wind_resource = wind_resource

hopp/eco/hopp_mgmt.py

@@ -21,8 +21,8 @@ def setup_hopp(
     hopp_site_input_data["lat"] = plant_config["project_location"]["lat"]
     hopp_site_input_data["lon"] = plant_config["project_location"]["lon"]
     hopp_site_input_data["year"] = plant_config["wind_resource_year"]
-    hopp_site_input_data["no_wind"] = not plant_config["project_parameters"]["wind"]
-    hopp_site_input_data["no_solar"] = not plant_config["project_parameters"]["solar"]
+    wind = plant_config["project_parameters"]["wind"]
+    solar = plant_config["project_parameters"]["solar"]
 
     # set desired schedule based on electrolyzer capacity
     desired_schedule = [plant_config["electrolyzer"]["rating"]] * 8760
@@ -33,6 +33,8 @@ def setup_hopp(
         hopp_site_input_data,
         hub_height=turbine_config["hub_height"],
         desired_schedule=desired_schedule,
+        wind=wind,
+        solar=solar
     )
 
     # replace wind data with previously downloaded and adjusted wind data

hopp/simulation/technologies/sites/site_info.py

@@ -1,19 +1,32 @@
+from typing import Optional, Union
+from pathlib import Path
+
+from attrs import define, field
 import matplotlib.pyplot as plt
-from shapely.geometry import *
-from shapely.geometry.base import *
-from shapely.validation import make_valid
-from fastkml import kml
+import numpy as np
+from numpy.typing import NDArray
+from shapely.geometry import Polygon, MultiPolygon, Point 
+from shapely.geometry.base import BaseGeometry
 from shapely.ops import transform
+from shapely.validation import make_valid
+from fastkml import kml, KML
 import pyproj
 import utm
 
-from hopp.simulation.technologies.resource.solar_resource import SolarResource
-from hopp.simulation.technologies.resource.wind_resource import WindResource
-from hopp.simulation.technologies.resource.elec_prices import ElectricityPrices
+from hopp.simulation.technologies.resource import (
+    SolarResource,
+    WindResource,
+    ElectricityPrices
+)
 from hopp.simulation.technologies.layout.plot_tools import plot_shape
 from hopp.utilities.log import hybrid_logger as logger
 from hopp.utilities.keys import set_nrel_key_dot_env
-
+from hopp.type_dec import (
+    hopp_array_converter as converter, NDArrayFloat, resource_file_converter,
+    hopp_float_type
+)
+from hopp.simulation.base import BaseClass
+from hopp.utilities.log import hybrid_logger as logger
 
 def plot_site(verts, plt_style, labels):
     for i in range(len(verts)):
@@ -25,140 +38,123 @@ def plot_site(verts, plt_style, labels):
 
     plt.grid()
 
-
-class SiteInfo:
-    """
-    Site specific information
-
-    Attributes
-    ----------
-    data : dict 
-        dictionary of initialization data
-    lat : float
-        site latitude [decimal degrees]
-    long : float
-        site longitude [decimal degrees]
-    vertices : np.array
-        site boundary vertices [m]
-    polygon : shapely.geometry.polygon
-        site polygon
-    valid_region : shapely.geometry.polygon
-        `tidy` site polygon
-    solar_resource : :class:`hybrid.resource.SolarResource`
-        class containing solar resource data
-    wind_resource : :class:`hybrid.resource.WindResource`
-        class containing wind resource data
-    elec_prices : :class:`hybrid.resource.ElectricityPrices`
-        Class containing electricity prices
-    n_timesteps : int
-        Number of timesteps in resource data
-    n_periods_per_day : int
-        Number of time periods per day
-    interval : int
-        Number of minutes per time interval 
-    urdb_label : string
-        `Link Utility Rate DataBase <https://openei.org/wiki/Utility_Rate_Database>`_ label for REopt runs
-    capacity_hours : list
-        Boolean list where ``True`` if the hour counts for capacity payments, ``False`` otherwise
-    desired_schedule : list
-        Absolute desired load profile [MWe]
-    follow_desired_schedule : boolean
-        ``True`` if a desired schedule was provided, ``False`` otherwise
+@define
+class SiteInfo(BaseClass):
     """
+    Represents site-specific information needed by the hybrid simulation class and layout optimization.
 
-    def __init__(self, data,
-                 solar_resource_file="",
-                 wind_resource_file="", 
-                 grid_resource_file="",
-                 hub_height=97,
-                 capacity_hours=[],
-                 desired_schedule=[]):
-        """
-        Site specific information required by the hybrid simulation class and layout optimization.
-
-        :param data: dict, containing the following keys:
-
-            #. ``lat``: float, latitude [decimal degrees]
-            #. ``lon``: float, longitude [decimal degrees]
-            #. ``year``: int, year used to pull solar and/or wind resource data. If not provided, default is 2012 [-]
-            #. ``elev``: float (optional), elevation (metadata purposes only) [m] 
-            #. ``tz``: int (optional), timezone code (metadata purposes only) [-]
-            #. ``no_solar``: bool (optional), if ``True`` solar data download for site is skipped, otherwise solar resource is downloaded from NSRDB
-            #. ``no_wind``: bool (optional), if ``True`` wind data download for site is skipped, otherwise wind resource is downloaded from wind-toolkit
-            #. ``site_boundaries``: dict (optional), with the following keys:
+    Attributes:
+        data (dict): Dictionary of initialization data.
+        lat (numpy.float64): Site latitude in decimal degrees.
+        lon (numpy.float64): Site longitude in decimal degrees.
+        tz (int, optional): Timezone code for metadata purposes only. Defaults to None.
+        vertices (:obj:`NDArray`): Site boundary vertices in meters.
+        polygon (:obj:`shapely.geometry.polygon.Polygon`): Site polygon.
+        valid_region (:obj:`shapely.geometry.polygon.Polygon`): Tidy site polygon.
+        solar_resource (:obj:`hopp.simulation.technologies.resource.SolarResource`): Class containing solar resource data.
+        wind_resource (:obj:`hopp.simulation.technologies.resource.WindResource`): Class containing wind resource data.
+        elec_prices (:obj:`hopp.simulation.technologies.resource.ElectricityPrices`): Class containing electricity prices.
+        n_timesteps (int): Number of timesteps in resource data.
+        n_periods_per_day (int): Number of time periods per day.
+        interval (int): Number of minutes per time interval.
+        urdb_label (str): Link to `Utility Rate DataBase <https://openei.org/wiki/Utility_Rate_Database>`_ label for REopt runs.
+        capacity_hours (:obj:`NDArray`): Boolean list indicating hours that count for capacity payments.
+        desired_schedule (:obj:`NDArray`): Absolute desired load profile in MWe.
+        follow_desired_schedule (bool): Indicates if a desired schedule was provided. Defaults to False.
 
-                * ``verts``: list of list [x,y], site boundary vertices [m]
-                * ``verts_simple``: list of list [x,y], simple site boundary vertices [m]
+    Args:
+        data (dict): Dictionary containing site-specific information.
+        solar_resource_file (Union[Path, str], optional): Path to solar resource file. Defaults to "".
+        wind_resource_file (Union[Path, str], optional): Path to wind resource file. Defaults to "".
+        grid_resource_file (Union[Path, str], optional): Path to grid pricing data file. Defaults to "".
+        hub_height (float, optional): Turbine hub height for resource download in meters. Defaults to 97.0.
+        capacity_hours (:obj:`NDArray`, optional): Boolean list indicating hours for capacity payments. Defaults to [].
+        desired_schedule (:obj:`NDArray`, optional): Absolute desired load profile in MWe. Defaults to [].
+        solar (bool, optional): Whether to set solar data for this site. Defaults to True.
+        wind (bool, optional): Whether to set wind data for this site. Defaults to True.
+    """
+    # User provided
+    data: dict
+    solar_resource_file: Union[Path, str] = field(default="", converter=resource_file_converter)
+    wind_resource_file: Union[Path, str] = field(default="", converter=resource_file_converter)
+    grid_resource_file: Union[Path, str] = field(default="", converter=resource_file_converter)
+    hub_height: hopp_float_type = hopp_float_type(97.)
+    capacity_hours: NDArray = field(default=[], converter=converter(bool))
+    desired_schedule: NDArrayFloat = field(default=[], converter=converter())
+    solar: bool = True
+    wind: bool = True
 
-            #. ``kml_file``: string (optional), filepath to KML with "Boundary" and "Exclusion" Placemarks
-            #. ``urdb_label``: string (optional), `Link Utility Rate DataBase <https://openei.org/wiki/Utility_Rate_Database>`_ label for REopt runs
+    # Set in post init hook
+    n_timesteps: int = field(init=False, default=None)
+    lat: hopp_float_type = field(init=False)
+    lon: hopp_float_type = field(init=False)
+    year: int = field(init=False, default=2012)
+    tz: Optional[int] = field(init=False, default=None)
+    solar_resource: Optional[SolarResource] = field(init=False, default=None)
+    wind_resource: Optional[WindResource] = field(init=False, default=None)
+    elec_prices: Optional[ElectricityPrices] = field(init=False, default=None)
+    n_periods_per_day: int = field(init=False)
+    interval: int = field(init=False)
+    follow_desired_schedule: bool = field(init=False)
+    polygon: Union[Polygon, BaseGeometry] = field(init=False)
+    vertices: NDArrayFloat = field(init=False)
+    kml_data: Optional[KML] = field(init=False, default=None)
 
-            .. TODO: Can we get rid of verts_simple and simplify site_boundaries
+    # .. TODO: Can we get rid of verts_simple and simplify site_boundaries
 
-        :param solar_resource_file: string, location (path) and filename of solar resource file (if not downloading from NSRDB)
-        :param wind_resource_file: string, location (path) and filename of wind resource file (if not downloading from wind-toolkit)
-        :param grid_resource_file: string, location (path) and filename of grid pricing data 
-        :param hub_height: int (default = 97), turbine hub height for resource download [m]
-        :param capacity_hours: list of booleans, (8760 length) ``True`` if the hour counts for capacity payments, ``False`` otherwise
-        :param desired_schedule: list of floats, (8760 length) absolute desired load profile [MWe]
-        """
+    def __attrs_post_init__(self):
         set_nrel_key_dot_env()
-        self.data = data
+
+        data = self.data
         if 'site_boundaries' in data:
             self.vertices = np.array([np.array(v) for v in data['site_boundaries']['verts']])
-            self.polygon: Polygon = Polygon(self.vertices)
+            self.polygon = Polygon(self.vertices)
             self.polygon = self.polygon.buffer(1e-8)
         if 'kml_file' in data:
             self.kml_data, self.polygon, data['lat'], data['lon'] = self.kml_read(data['kml_file'])
             self.polygon = self.polygon.buffer(1e-8)
+
         if 'lat' not in data or 'lon' not in data:
             raise ValueError("SiteInfo requires lat and lon")
         self.lat = data['lat']
         self.lon = data['lon']
-        self.n_timesteps = None
+
         if 'year' not in data:
             data['year'] = 2012
+        if 'tz' in data:
+            self.tz = data['tz']
 
-        if 'no_solar' not in data:
-            data['no_solar'] = False
-
-        if not data['no_solar']:
-            self.solar_resource = SolarResource(data['lat'], data['lon'], data['year'], filepath=solar_resource_file)
+        if self.solar:
+            self.solar_resource = SolarResource(data['lat'], data['lon'], data['year'], filepath=self.solar_resource_file)
             self.n_timesteps = len(self.solar_resource.data['gh']) // 8760 * 8760
 
-        if 'no_wind' not in data:
-            data['no_wind'] = False
-
-        if not data['no_wind']:
+        if self.wind:
             # TODO: allow hub height to be used as an optimization variable
-            self.wind_resource = WindResource(data['lat'], data['lon'], data['year'], wind_turbine_hub_ht=hub_height,
-                                            filepath=wind_resource_file)
+            self.wind_resource = WindResource(data['lat'], data['lon'], data['year'], wind_turbine_hub_ht=self.hub_height,
+                                            filepath=self.wind_resource_file)
             n_timesteps = len(self.wind_resource.data['data']) // 8760 * 8760
             if self.n_timesteps is None:
                 self.n_timesteps = n_timesteps
             elif self.n_timesteps != n_timesteps:
                 raise ValueError(f"Wind resource timesteps of {n_timesteps} different than other resource timesteps of {self.n_timesteps}")
 
-        self.elec_prices = ElectricityPrices(data['lat'], data['lon'], data['year'], filepath=grid_resource_file)
+        self.elec_prices = ElectricityPrices(data['lat'], data['lon'], data['year'], filepath=self.grid_resource_file)
         self.n_periods_per_day = self.n_timesteps // 365  # TODO: Does not handle leap years well
         self.interval = int((60*24)/self.n_periods_per_day)
         self.urdb_label = data['urdb_label'] if 'urdb_label' in data.keys() else None
 
-        if len(capacity_hours) == self.n_timesteps:
-            self.capacity_hours = capacity_hours
-        else:
-            self.capacity_hours = [False] * self.n_timesteps
+        if len(self.capacity_hours) != self.n_timesteps:
+            self.capacity_hours = np.array([False] * self.n_timesteps)
 
         # Desired load schedule for the system to dispatch against
-        self.desired_schedule = desired_schedule
-        self.follow_desired_schedule = len(desired_schedule) == self.n_timesteps
-        if len(desired_schedule) > 0 and len(desired_schedule) != self.n_timesteps:
+        self.follow_desired_schedule = len(self.desired_schedule) == self.n_timesteps
+        if len(self.desired_schedule) > 0 and len(self.desired_schedule) != self.n_timesteps:
             raise ValueError('The provided desired schedule does not match length of the simulation horizon.')
-
             # FIXME: this a hack
-        if 'no_wind' in data and data["no_wind"]:
+
+        if not self.wind:
             logger.info("Set up SiteInfo with solar resource files: {}".format(self.solar_resource.filename))
-        elif 'no_solar' in data and data["no_solar"]:
+        elif not self.solar:
             logger.info("Set up SiteInfo with wind resource files: {}".format(self.wind_resource.filename))
         else:
             logger.info(
@@ -220,6 +216,9 @@ def plot(self,
         return figure, axes
 
     def kml_write(self, filepath, turb_coords=None, solar_region=None, wind_radius=200):
+        if self.kml_data is None:
+            raise AttributeError("No KML data to write.")
+
         if turb_coords is not None:
             turb_coords = np.atleast_2d(turb_coords)
             for n, (x, y) in enumerate(turb_coords):

hopp/to_organize/hopp_tools_steel.py

@@ -114,11 +114,10 @@ def set_site_info(hopp_dict, site_df, sample_site):
     lon = float(lon)
     sample_site['lat'] = lat
     sample_site['lon'] = lon
-    sample_site['no_solar'] = False
     # if solar_size_mw>0:
-    #     sample_site['no_solar'] = False
+    #     sample_site['solar'] = True
     # else:
-    #     sample_site['no_solar'] = True
+    #     sample_site['solar'] = False
 
     hopp_dict.add('Configuration', {'sample_site': sample_site})
 

hopp/tools/dispatch/csp_pv_battery_plot.py

@@ -65,8 +65,7 @@ def init_hybrid_plant():
         "elev": 641,
         "year": 2012,
         "tz": -8,
-        "no_wind": True
-        }
+    }
 
     root = "C:/Users/WHamilt2/Documents/Projects/HOPP/CSP_PV_battery_dispatch_plots/"
     solar_file = root + "34.865371_-116.783023_psmv3_60_tmy.csv"
@@ -84,7 +83,8 @@ def init_hybrid_plant():
     site = SiteInfo(site_data,
                     solar_resource_file=solar_file,
                     grid_resource_file=prices_file,
-                    desired_schedule=desired_schedule
+                    desired_schedule=desired_schedule,
+                    wind=False
                     )
 
     technologies = {'tower': {

hopp/tools/hopp_tools.py

@@ -39,7 +39,6 @@ def set_site_info(site_df, sample_site):
     lon = float(lon)
     sample_site['lat'] = lat
     sample_site['lon'] = lon
-    sample_site['no_solar'] = False
 
     return site_df, sample_site