0

空間参照 WSG84 で gdal2tiles.py を使用していましたが、正常に動作していました。ただし、空間参照 EPSG:3440 に変更すると。次のエラーが発生しました。

まず、次のコマンドを実行します。

gdal_translate -of VRT -a_srs EPSG:3440 -gcp 0 0 58.068451479718924 23.65512391903488 -gcp 21816 0 58.126966134442846 23.65512391903488 -gcp 21816 14871 58.126966134442846 23.6185834507829  myinputimage.png myoutput.vrt

それより、次のコマンドを実行します。

python gdal2tiles.py  myoutput.vrt

次のエラーが発生しました。 ここに画像の説明を入力

このエラーを修正する方法、またはこのエラーのない同様の空間参照を使用する方法はありますか? このエラーは、WSG84 以外の空間参照を使用している場合にのみ発生することに注意してください。

編集:問題は、空間参照 EPSG:3440 が含まれていない可能性があると思いますが、gdal 参照データベースに追加する必要があるかどうかという問題が残ります。

4

1 に答える 1

1

EPSG:3440 は投影システムです。精度 0.5 メートルの正しい定義は、おそらく次のようになります。

+proj=utm +zone=40 +ellps=clrk80 +towgs84=-180.624,-225.516,173.919,-0.81,-1.898,8.336,16.71006 +units=m +no_defs

http://epsg.io/3440-1439を見ると、次のように表示されます。

Projected bounds:
-35283.71 1840661.51
789562.66 2942956.58

上記のコマンドでは、58.068451479718924 などの測地 (緯度経度番号) のように見える地上基準点を割り当てます。これは間違っている可能性が非常に高いです - 割り当てた数値が予測された境界範囲内にないためです。

gdal_translate の数値の -gcp パラメータで使用する必要があります。これは、既に投影されており、上記の範囲内にあります。

オンライン インターフェイス ( http://epsg.io/3440-1439/map )を介して WGS84 から緯度/ 経度を変換し、変換された数値を -gcp 引数に割り当てることができます。

http://www.maptiler.com/maptilerを参照するタグでクエリをマークしました。このツールは、VRT を作成する必要なく視覚的にも座標を割り当てることができるため、画像をクリックして、投影されたシステムで上記の座標をコピー & ペーストできます。MapTiler は、OpenLayers で EPSG:3440 タイルを表示するビューアを生成するか、これらを Google や OSM とのオーバーレイやモバイル アプリで使用するための球状メルカトル タイルに変換します。

ところで、私たちが開発したシステム gdal2tiles、maptiler、およびhttp://epsg.io/を使用していると聞いてうれしいです...

于 2015-09-16T17:06:14.427 に答える