擅长:python、mysql、java
<p>该实用程序可以用一个调用执行多个坐标。例如,使用坐标对准备一个简单的<code>coords.txt</code>文本文件:</p>
<pre><code>1.0 2.0
3.0 4.0
5.0 6.0
</code></pre>
<p>然后从一个OSGeo4W外壳将其导入并输出:</p>
^{pr2}$
<p>它将生成一个<code>values.txt</code>文件,其中包含每个坐标的值。可以使用管道stdin和stdout参数对Popen执行相同的操作。在</p>
<pre><code>coords = [(1.0, 2.0), (3.0, 4.0), (5.0, 6.0)]
coord_txt = ''.join(['{0} {1}\n'.format(*c) for c in coords])
p = subprocess.Popen(
['gdallocationinfo', r'C:\path\to\raster.tif', '-wgs84', '-valonly'],
stdin=subprocess.PIPE, stdout=subprocess.PIPE, universal_newlines=True, shell=True)
values_txt, err = p.communicate(coord_txt)
values = values_txt.splitlines()
</code></pre>
<p><code>values</code>将是一个值列表,长度与<code>coords</code>相同。在</p>