GDALを使用して1440x720データセットから4096x4096 pngを作成する方法


8

かなり簡単に思えるようなことをやろうとしています。私はftp://eclipse.ncdc.noaa.gov/pub/OI-daily-v2/NetCDF/2011/AVHRR-AMSR/からnetcdfファイルを取得し、4096x4096 pngに変換してタイル表示できるようにしていますグーグルマップで使用するためにそれを。私は次のコマンドを使用してきましたが、東と西の境界に常に空白の領域が作成されるようです。

/usr/local/bin/gdal_translate -a_srs EPSG:4326 -outsize 4096 4096 -a_ullr 0 90 360 -90 NETCDF:"input.nc":sst -of netCDF output.nc
/usr/local/bin/gdaldem color-relief output.nc -b 1 -alpha colorscale.txt -of PNG output.png
/usr/bin/convert output.png -transparent #FFFFFF output.png

私はこれに間違った方法でアプローチしていますか?間違ったツールを使用していますか?

FWIW-gdal_translateはv1.8.0です

前もって感謝します。

更新: Dija(私の同僚)からの以下の "Answer"は、実際には質問自体に対するコメントへの応答です。

更新 コメントではなく誤って回答を追加しました。ここに移動します。

私はここで体重を測ると思いました。私はこの問題にしばらく取り組んできました。私はプロセスの多くの反復を試みてきました。これらの画像を実際に4096x4096にスケーリングしました。つまり、画像をネイティブサイズのままにすると、(このデータセットの)子午線と他のデータセットの180番目の子午線では、1440x720のギャップが大きくなります。出力画像を256の倍数にスケーリングして(Googleにうまくフィットさせるため)、さまざまな結果を得ました。

1536x768にスケーリングすると、ギャップは0子午線の右側まで表示されます。1280x640に拡大すると、左側になります。そのままにしておけば、かなり分割されてしまいます。

私はgdaldemが吐き出す基本画像を撮り、それを半分にカットし、エッジを一致させました...それらは完全に整列しています。これにより、gdal2tiles.pyで奇妙なことが起こっていると思います

私は自分のタイラーを書いて、ギャップが残っているかどうかを確認するために、画像を緩める準備をしています。

スクリーンショットをいくつか組み合わせて、話している内容とギャップの方法を示し、タイラーの機能についてレポートします。

ご意見をお寄せいただきありがとうございます。ありがたいです。

更新2オリジナルのラスターに対して作成したマップティラーを実行しましたが、ギャップがなくても問題なく動作しました。gdal2tiles.pyが問題のようです。180番目の子午線に接する最後のタイルに透明なエッジを導入します。


アスペクト比が1:2の画像をアスペクト比が1:1:の画像に変換すると、クリッピングが発生するか、空白の境界線が空のスペースを埋めるか、2:1のストレッチが必要になります。
whuber

私はあなたが提供したFTPサイトからのファイルの1つで同じプロセスを実行しました、そして私はその境界に適切に適合する画像を持っています。得る効果のスクリーンショットを提供できますか?効果はgdaldemステップまたは変換ステップの後に発生しますか?
MerseyViking 2011

ええ、そうです、whuberが言うように、テキストで提供したパラメーターを使用する場合、アスペクト比は間違っていますが、私が使用したものであるコマンドラインパラメーターは正しいアスペクト比を持っています。
マージーバイキング、2011

そうは言っても、4096 x 4096として作成しようとしたところ、QGISは問題なくロードしました。さらに、両方のバージョンをGIMPにロードすると、左端と右端まで右に行くマップが表示されます。
MerseyViking 2011

@whuper-もちろん、私は完全に同意します。1:2のデータセットを1:1の画像にしようとすると、ほとんどデータが無効になることに気づきました。そうは言っても、私の期待は、空のスペースではなく、基本的に2:1ストレッチの出力を取得することでした。私は自分の前に設定されたビジネス要件のなすがままになっています:(ああ、きれいな絵のためにやることです。フィードバックをありがとうございました
11

回答:


2

gdal_gridには、出力ラスターの範囲座標を指定できるパラメーターがあります。

gdal_grid -txe [xcoord1] [xcoord2] -tye [ycoord1] [ycoord2]

左下の座標は(xcoord1、ycoord1)、右上は(xcoord2、ycoord2)

GDALは、出力範囲に従ってデータを補間します

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.