def project_tilt_to_pa(tilt, inclination):
"""
This function ...
:param tilt:
:param inclination:
:return:
"""
# Get the tilt angle and inclination in radians
tilt_radian = tilt.to("radian").value
i_radian = inclination.to("radian").value
denominator = math.sqrt(math.sin(tilt_radian)**2 * math.sin(i_radian)**2 + math.cos(tilt_radian)**2)
cos_pa = math.sin(tilt_radian) * math.sin(i_radian) / denominator
sin_pa = math.cos(tilt_radian) / denominator
pa_radian = math.atan2(sin_pa, cos_pa) * Unit("radian")
return pa_radian.to("deg")
# -----------------------------------------------------------------
评论列表
文章目录