2016-11-04 9 views
1

gdal_gridを使用して、一部の3Dサーフェスの標高モデルを作成しています。ポイントでgdal_gridを使用するには

私は、このコマンドでにGeoJSON・ファイルを使用してそれを行うことができます。

ds2 = gdal.Grid('outputfile.tif', 'inputfile.geojson', format = 'GTiff', algorithm = 'linear:radius=0') 

これは正常に動作しますが、私は個別に各機能のためにそれを行うことができるようにしたいと思います。

[[12.135253194446484, 55.590235278979236, 44.500800000000005], 
[12.136885609925141, 55.58968131535586, 44.500800000000005], 
[12.149742647277185, 55.59946751368944, 89.5008], 
[12.14443275453964, 55.601269628832526, 89.5008], 
[12.135253194446484, 55.590235278979236, 44.500800000000005]] 

私の質問は、したがって、次のとおりです:

  1. することができます私は、私は、例えばにGeoJSONファイルをループをすることができますし、それぞれの特徴を得るが、ただの点でgdal.Gridを使用する方法がありますgeojsonの代わりにポイントでgdal.Gridを使用しますか?
  2. どこでgdal.Gridに使用できる入力パラメータを正確に見ることができますか?

答えて

0

これは私がこの問題を解決した方法です。最も洗練されたソリューションではないかもしれませんが、うまくいくようです。 geojsonファイル(辞書として)からサーフェスを読み込み、最初のフィーチャを取得してjson-stringに変換します。

with open(surfaceFileName,'r') as file: 
    data = json.load(file) 
# the first feature: 
dataJson = json.dumps(data['features'][0]['geometry']) 
# this feature as geojson-string 
featureJson = """{"type":"FeatureCollection", 
       "features": [ 
       {"type": "Feature", 
       "geometry": """+dataJson+""", 
       "properties": {} 
       }]}""" 
# Using gdal_grid: 
ds2 = gdal.Grid('test10py.tif', featureJson, format = 'GTiff', algorithm = 'linear:radius=0') 
関連する問題