npの配列を持つPythonでgdal_gridを使用

Question

私はPython 2.7を使用しています。 gdal 2.1.2ライブラリがインストールされています。

私はこのようなNPアレイを有する第1及び第2の列は私の座標（具体的には、緯度及び経度）を表し、第三の一つの位置に関連した変数値を収集

[[ 2.34845824 48.84626174 15.20369 ] 
[ 2.34966283 48.84910128 13.88528 ] 
[ 2.35120647 48.85098931 15.76322 ] 
..., 
[ 2.3556567 48.8415611 15.42184 ] 
[ 2.34394085 48.84248672 15.51128 ] 
[ 2.34926763 48.85128456 15.0685 ]] 

。私はgdal_gridをPythonで直接使用して画像を取得したいと思います。ピクセルごとに逆距離重み付けなどの方法を使って配列の値を補間する方法があります。

私は私のデータをCSVファイルに収集された場合、gdal_gridでそれを行う方法を見出した（私は「のtest.CSV」を呼ぶこと）と私はこのような「test.vrt」ファイル作成する場合：

を

gdal.Grid('test_python.tiff', 'test.vrt', layers='test', algorithm='invdistnn:max_points=10:min_points=1:radius=0.001', format='GTiff', width=333, height=240) 

をし、それが動作します：

<OGRVRTDataSource> 
    <OGRVRTLayer name="test"> 
     <SrcDataSource>test.csv</SrcDataSource> 
     <GeometryType>wkbPoint</GeometryType> 
     <LayerSRS>WGS84</LayerSRS> 
     <GeometryField separator=" " encoding="PointFromColumns" x="field_1" y="field_2" z="field_3"/> 
    </OGRVRTLayer> 
</OGRVRTDataSource> 

私はコマンドでpythonでそれを実行しようとしました。問題は、この操作を何回か繰り返す必要があるため、np配列を.csvに変換しないようにしたいということです。 どうすればnp配列を直接使用できますか？

ご協力いただきありがとうございます。

Answer 1

できません。

への入力は、GDALで判読可能な形式（VRTが1）でなければなりません。または、numpy配列をTiff/GeoTiffファイルとして保存することもできます。

GDAL> = 2.1の別の方法では、コマンドラインユーティリティではなくpython bindingを使用し、numpy配列（少なくともディスクに保存しないでください）から読み込まれた 'MEM'データセットを使用します。