方位角投影代码python

下面是一个简单的Python代码示例，用于计算和绘制方位角投影。

import matplotlib.pyplot as plt
import numpy as np

def azimuthal_projection(lon, lat, center_lon, center_lat):
    # Convert degrees to radians
    lon, lat, center_lon, center_lat = np.radians([lon, lat, center_lon, center_lat])
    
    # Calculate the distance and azimuth from the center point
    delta_lon = lon - center_lon
    cos_lat = np.cos(lat)
    sin_lat = np.sin(lat)
    cos_center_lat = np.cos(center_lat)
    sin_center_lat = np.sin(center_lat)
    d = np.arccos(sin_center_lat * sin_lat + cos_center_lat * cos_lat * np.cos(delta_lon))
    az = np.arctan2(cos_lat * np.sin(delta_lon), cos_center_lat * sin_lat - sin_center_lat * cos_lat * np.cos(delta_lon))
    
    # Convert back to degrees
    d = np.degrees(d)
    az = np.degrees(az)
    
    # Return the distance and azimuth
    return d, az

# Example usage
lon = np.arange(-180, 180, 5)
lat = np.arange(-90, 90, 5)
center_lon, center_lat = 0, 0

fig = plt.figure(figsize=(8, 8))
ax = fig.add_subplot(111, projection='polar')
ax.set_theta_zero_location('N')
ax.set_theta_direction(-1)

for lon_i in lon:
    for lat_i in lat:
        d, az = azimuthal_projection(lon_i, lat_i, center_lon, center_lat)
        ax.plot(np.radians(az), d, marker='o', markersize=5, color='r')

plt.show()

此代码使用matplotlib库绘制了极坐标图，其中点的位置代表了经纬度坐标的方位角投影。在此示例中，我们使用经度和纬度每隔5度的网格来绘制所有点的方位角投影。