球形地球上的2D矢量场的分歧(Python /底图)

时间:2016-10-25 21:42:10

标签: python vector matplotlib-basemap

我在几个点上有风的u(纬向)和v(子午)成分的二维阵列,其位置给出了相应的2D纬度和经度阵列。

我想计算这个风场的分歧。我知道分歧只是du / dx + du / dy,但由于我的x和y坐标是纬度和经度,所以它并不是很正确。例如,地球上到处都是均匀的子午线风(v =常数到处)实际上在民意调查中会有分歧/收敛,即使v阵列在每个指数上只是一个常数。

我最终想在Basemap实例上绘制分歧。

我该怎么做呢?

1 个答案:

答案 0 :(得分:2)

发散的公式取决于您发现的坐标系。计算公式(使用坐标公式和链规则的变化)是一项值得的练习,但wikipedia有一个有用的公式表。

您需要球面坐标的公式,并且由于您正在球体的表面上工作,因此您需要按住 r 。这具有删除第一个术语和缩放其他两个术语的效果。