GTiffラスターの時系列スタックを単一のNetCDFに変換する
gdal-devメーリングリストからの移行: 2013年9月2日月曜日午後7時9分、David Sheanは次のように書いています。 こんにちはリスト、私は配布用の単一のNetCDFファイルとして同一の投影/範囲/解像度を持つGTiffラスターの時系列をパッケージ化しようとしています。過去1時間、オンラインドキュメントを参照し、gdal_translate、gdalbuildvrt、gdalwarpで遊んでみましたが、成功しませんでした。 既存のgdalコマンドラインユーティリティを使用してこれを行う簡単な方法はありますか?NetCDF Python APIを使用したカスタムソリューションに頼る前に、私は尋ねたいと思いました。 ありがとう。-デビッド 2013年9月3日火曜日、午前10時15分、エティエンヌトゥーリニーは次のように書いています。 あなたが望むのはおそらくgdalの範囲外です。gdal_translateがそれらを単一のファイルに入れるためには、巧妙なメタデータ管理が必要になります... gdal_translateを使用してそれらをすべてnetcdfに変換してから、python-netcdf4(numpy / scipyのものではない)を使用して時間ディメンションにスタックすることをお勧めします。 2013年9月3日火曜日午前7時55分、「Signell、Richard」は次のように書いています。 David、質問をGIS stackexchangeグループ/gis// に投稿する場合、 参考になるサンプルコードを提供します。 -リッチ ===================== 更新日9/3/13 17:04 PDT 入力データセットの1つに対するgdalinfoの出力を次に示します。 gdalinfo 20120901T2024_align_x+22.19_y+3.68_z+14.97_warp.tif Driver: GTiff/GeoTIFF Files: 20120901T2024_align_x+22.19_y+3.68_z+14.97_warp.tif Size is 10666, 13387 Coordinate System is: PROJCS["unnamed", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]], PROJECTION["Polar_Stereographic"], PARAMETER["latitude_of_origin",70], …