擅长:python、mysql、java
<p>这里,在这个代码中,alpha是0.025</p>
<pre><code>import numpy as np
vp_ref = 30
surf_lay = np.ones(([1000,10000]), dtype=np.float32);
gvp_ref = vp_ref * surf_lay
alpha = 0.025
substration_array = np.array([[alpha*i]*gvp_ref.shape[1] for i in range(gvp_ref.shape[0])])
gvp_ref.shape
substration_array.shape
output = np.subtract(gvp_ref, substration_array)
</code></pre>