擅长:python、mysql、java
<p>假设数据在一个<em>n</em>x3numpy数组中,首先将elevation列重新解释为矩阵(表示统一网格):</p>
<pre><code>m=data[:,2].reshape(ny,nx)
</code></pre>
<p>然后进行几次切片和减法,得到细胞中心的导数:</p>
^{pr2}$
<p>该系数校正单位(否则为米/单元格</em>,而不是每米),并将总和转换为平均值。(如果每个维度中的间距不同,您将分别将参数缩放为<code>hypot</code>。)请注意,结果数组在每个维度上都比输入小一个;如果大小需要相同,可以使用更复杂的差分方案。<code>numpy.gradient</code>实现了其中的一些,允许</p>
<pre><code>mag=numpy.hypot(*numpy.gradient(m,abs(data[1][1]-data[1][0])))
</code></pre>