calculate_daylight

[gregory-halverson-jpl]

Calculate the number of daylight hours for a given day and location.

This function is flexible in its inputs:
- If SHA_deg (sunrise hour angle in degrees) is provided, it is used directly.
- Otherwise, SHA_deg is computed from the day of year (DOY) and latitude (lat).
- If lat is not provided but a geometry object is, latitude is extracted from geometry.
- If DOY is not provided but a datetime_UTC is, DOY is computed from the datetime.

Parameters
----------
DOY : Union[Raster, np.ndarray, int], optional
    Day of year (1-366). Can be a Raster, numpy array, or integer.
lat : Union[Raster, np.ndarray, float], optional
    Latitude in degrees. Can be a Raster, numpy array, or float.
SHA_deg : Union[Raster, np.ndarray, float], optional
    Sunrise hour angle in degrees. If not provided, it will be calculated.
datetime_UTC : datetime, optional
    Datetime in UTC. Used to determine DOY if DOY is not provided.
geometry : SpatialGeometry, optional
    Geometry object containing latitude information.

Returns
-------
Union[Raster, np.ndarray]
    Daylight hours for the given inputs.