擅长:python、mysql、java
<blockquote>
<p>I have arrays associated with both A and B which store the (x,y) positions of the gridpoints for each.</p>
</blockquote>
<p>在这种情况下,答案应该很简单。。。在</p>
<p>这两个网格是否严格在同一网格方案上?假设他们是,你可以做一些类似的事情:</p>
<pre><code>np.argwhere((Ax == Bx.min()) & (Ay == By.min()))
</code></pre>
<p>假设两个网格的世界坐标以与网格的指示符相同的方向增加,这将给出子网格的左下角。(如果它们不在同一方向上增加(即负数<code>dx</code>或<code>dy</code>),则只给出另一个角点)</p>
<p>在下面的例子中,我们显然可以从<code>ix = (Bxmin - Axmin) / dx</code>等计算适当的指标,但是假设你有一个更复杂的网格系统,这仍然可以工作。然而,这是假设两个网格<em>在同一网格方案</em>上!如果他们不是。。。在</p>
^{pr2}$