src / pysoilmap / features.py

sun_exposure

sun_exposure(grad_dir, slope, latitude, declination=0.0)[source]

Return the cosine of the angle between surface normal and the sun at midday.

Parameters:

  • grad_dir – direction of the gradient in XY plane in radians in counter-clockwise direction with respect to east

  • slope – inclination of the gradient in radians

  • latitude – point latitude in radians

  • declination – sun declination angle in radians. Can be used to adjust for season. A positive value corresponds to a northward deviation. Zero means the sun is in the zenith over the equator.

All input arguments must be of the same shape or scalar.

This is very similar to rad_angle but not exactly the same.