python绘制极坐标轮廓图contourf

最新推荐文章于 2024-08-09 15:34:32 发布

小朱小朱绝不服输

最新推荐文章于 2024-08-09 15:34:32 发布

阅读量4.4k

点赞数 4

分类专栏： Python 文章标签： python contourf

本文链接：https://blog.csdn.net/weixin_44052055/article/details/118758568

版权

该博客介绍如何使用Python的极坐标系统绘制风速的轮廓图。首先计算每个象限的风速、距离和角度，然后通过`ax.contourf`在极坐标下绘制。文章特别提到在角度计算和插值过程中可能遇到的问题，以及解决方法，最后展示成功绘制的极坐标风速轮廓图。

摘要由CSDN通过智能技术生成

任务：将每个象限的风速，距离，角度绘制成极坐标轮廓图。
使用极坐标画等直线图，可以用极地图的ax.contour或ax.contourf 。

1.变量计算

每个象限的风速，距离就不再说怎么画了，这里说下角度的计算。

两个经纬度之间的夹角（与正北方向的夹角）：
在这里插入图片描述

# 点1到点2方向沿逆时针方向转到正北方向的夹角
def LatLng2Degree(LatZero,LngZero,Lat,Lng):
    """
    Args:
        point p1(latA, lonA)
        point p2(latB, lonB)
    Returns:
        bearing between the two GPS points,
        default: the basis of heading direction is north
    """
    radLatA = math.radians(LatZero)
    radLonA = math.radians(LngZero)
    radLatB = math.radians(Lat)
    radLonB = math.radians(Lng)
    dLon = radLonB - radLonA
    y = math.sin(dLon) * math.cos(radLatB)
    x = math.cos(radLatA) * math.sin(radLatB) - math