球形地球上二维矢量场的散度 (Python/Basemap)
Divergence of a 2D vector field on a spherical Earth (Python/Basemap)
我在几个点上有风的 u（纬向）和 v（经向）分量的二维数组，它们的位置给出了相应的二维纬度和经度数组。
我想计算这个风场的散度。我知道分歧只是 du/dx + du/dy，但是因为我的 x 和 y 坐标是纬度和经度，所以结果不太正确。例如，地球上任何地方的均匀经向风（v = 到处都是常数）实际上在民意调查中会有 divergence/convergence，即使 v 数组在每个索引处只是一个常数。
我最终想在底图实例上绘制散度图。
我应该怎么做？
散度公式取决于您发现的坐标系。计算公式（使用坐标变化公式和链式法则）是一项值得练习的练习，但是 wikipedia 有一个有用的 table 必要公式。
您需要球坐标的公式，并且由于您在球体的表面上工作，所以您需要保留 r。这具有删除第一项并缩放其他两项的效果。
我在几个点上有风的 u（纬向）和 v（经向）分量的二维数组，它们的位置给出了相应的二维纬度和经度数组。
我想计算这个风场的散度。我知道分歧只是 du/dx + du/dy，但是因为我的 x 和 y 坐标是纬度和经度，所以结果不太正确。例如，地球上任何地方的均匀经向风（v = 到处都是常数）实际上在民意调查中会有 divergence/convergence，即使 v 数组在每个索引处只是一个常数。
我最终想在底图实例上绘制散度图。
我应该怎么做？
散度公式取决于您发现的坐标系。计算公式（使用坐标变化公式和链式法则）是一项值得练习的练习，但是 wikipedia 有一个有用的 table 必要公式。
您需要球坐标的公式，并且由于您在球体的表面上工作，所以您需要保留 r。这具有删除第一项并缩放其他两项的效果。