3Dポリゴンから標高モデルを作成するにはどうすればよいですか？

Question

私はgeojsonファイルから3次元で多数のポリゴンを持っています。標高モデルを作りたいと思います。つまり、すべてのピクセルがこの位置にあるポリゴンの高さであるラスタが必要です。

私はgdal_rasterizeを見てみましたが、説明は、今のよう

を言う、唯一の点や線が3Dで描画されます。

gdal_rasterize

Answer 1

私はgriddataと呼ばれるscipy.interpolat-機能を使用して終了。これは、グリッド内の座標を取得するためにmeshgridを使用しています。そのため、meshgridのメモリ制限のためにそれをタイルアップする必要がありました。

import scipy.interpolate as il #for griddata 
# meshgrid of coords in this tile 
gridX, gridY = np.meshgrid(xi[c*tcols:(c+1)*tcols], yi[r*trows:(r+1)*trows][::-1]) 

## Creating the DEM in this tile 
zi = il.griddata((coordsT[0], coordsT[1]), coordsT[2], (gridX, gridY),method='linear',fill_value = nodata) # fill_value to prevent NaN at polygon outline 

直線補間は、私が望むものとまったく同じようです。説明を参照してくださいhttps://docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.griddata.html