输入经纬度求太阳高度角和太阳方位角的代码

要计算太阳高度角和太阳方位角，你需要知道经纬度和日期/时间。以下是一个计算太阳高度角和太阳方位角的Python代码示例：

import math
import datetime

def calculate_sun_position(latitude, longitude, date_time):
    # 计算一天中的时间参数
    day_of_year = date_time.timetuple().tm_yday
    time_offset = date_time.utcoffset().total_seconds() / 3600

    # 计算太阳赤纬
    solar_declination = -23.45 * math.cos(math.radians(360 * (day_of_year + 10) / 365))

    # 计算时角
    solar_time = date_time.hour + (4 * longitude) - time_offset
    hour_angle = 15 * (solar_time - 12)

    # 计算太阳高度角
    sin_elevation = math.sin(math.radians(latitude)) * math.sin(math.radians(solar_declination)) + \
                    math.cos(math.radians(latitude)) * math.cos(math.radians(solar_declination)) * \
                    math.cos(math.radians(hour_angle))
    elevation_angle = math.degrees(math.asin(sin_elevation))

    # 计算太阳方位角
    cos_azimuth = (math.sin(math.radians(latitude)) * math.sin(math.radians(elevation_angle)) -
                   math.sin(math.radians(solar_declination))) / \
                  (math.cos(math.radians(latitude)) * math.cos(math.radians(elevation_angle)))

    if math.sin(math.radians(hour_angle)) > 0:
        azimuth_angle = math.degrees(math.acos(cos_azimuth))
    else:
        azimuth_angle = 360 - math.degrees(math.acos(cos_azimuth))

    return elevation_angle, azimuth_angle

# 示例用法
latitude = 37.7749  # 纬度
longitude = -122.4194  # 经度
date_time = datetime.datetime.now()  # 当前日期和时间

elevation, azimuth = calculate_sun_position(latitude, longitude, date_time)

print("太阳高度角：", elevation)
print("太阳方位角：", azimuth)

请注意，这段代码假设输入的日期和时间为当前日期和时间，你可以根据需要进行修改。另外，这段代码使用的经纬度是旧金山的经纬度，你可以将其替换为你所需的经纬度