マップタイル作成プロパティについて


10

タイルベースマップを作成しようとしています。タイルを作成したい画像があり、GDALを使用しています。

タイルを作成できましたが、結果を表示すると、タイルマップが次のように見えます。

  1. ひっくり返されます
  2. 比率が正しくありません

これは地図のオーバーレイ画像です: 地図画像

これはタイリングの結果です:

タイル結果

オルソフォトの画像のコーナーを比較して、座標が良いことを何度も確認しました。

もう1つは、同じ座標のGoogle Earthでは、マップ画像がGoogleマップに正確にフィットすることです。

グーグルアースオーバーレイ

編集済み

gdalinfo map_overlay.png

この出力では:

ドライバー:PNG / Portable Network Graphicsファイル:map_overlay.pngサイズは3527、2494座標系は `'画像構造メタデータ:
INTERLEAVE = PIXELコーナー座標:左上(0.0、0.0)左下(0.0、2494.0)右上(3527.0 、0.0)右下(3527.0、2494.0)中央(1763.5、1247.0)バンド1ブロック= 3527x1タイプ=バイト、ColorInterp =赤マスクフラグ:PER_DATASET ALPHAバンド2ブロック= 3527x1タイプ=バイト、ColorInterp =緑マスクフラグ:PER_DATASET ALPHAバンド3ブロック= 3527x1タイプ=バイト、ColorInterp =ブルー
マスクフラグ:PER_DATASET ALPHAバンド4ブロック= 3527x1タイプ=バイト、ColorInterp =アルファ

次に翻訳します:

gdal_translate -of VRT -a_srs EPSG:4326 -gcp 0 0 31.7431761644 35.1680410195 -gcp 3527 0 31.7493769674 35.1784535489 -gcp 3527 2492 31.7431011291 35.1784951643 -gcp 0 2494 31.7431761644 35.1680410195 map_overlay.pngrt_laylay.png map_overlay.png

この出力では:

入力ファイルのサイズは3527、2494です

そして最後に:

gdal2tiles.py -z 14-21 map_overlay.vrt

オーバーレイは、特別なデータのない単純な.png画像です。

私はどんな助けにも感謝します。

シャニ


2
gdal_translateここに投稿したパラメーターに不整合があるようです:(lat、lon)=(31.7431761644、35.1680410195)を(0,0)と(0,2494)の2つの異なる点に関連付けています。これは実際に発行したコマンドですか、それともコピーアンドペーストエラーですか?
whuber

回答:


5

次のようにタイリングする前に、メトリック調整画像を使用してみてください。

   gdalwarp -of GTiff -t_srs EPSG:3857 input.tif output.tif 

更新

注:gcpを指定してもgdal_translate、tiffのコーナー座標は指定されません。

これに加えて、ファイルが1つしかないため、仮想データセット(.vrt)を作成する必要はありません...

次のコードを使用します。

  gdal_translate -of  GTiff  -gcp 0 0 31.7431761644 35.1680410195 -gcp 3527 0 
       31.7493769674 35.1784535489 -gcp 3527 2492 31.7431011291 35.1784951643
      -gcp 0 2494 31.7431761644 35.1680410195 map_overlay.png map_overlay.tif

その後

  gdalwarp -s_srs epsg:4326 -t_srs epsg:4326  map_overlay.tif warped_map_overlay.tif

そして今gdalinfoを使用してみてください...問題がなければ、次のコマンドを実行します。

  gdal2tiles.py -z 14-21 warped_map_overlay.tif

タイルが合わない場合は、更新前に上記のコードを実行してください。

あなたはここでより多くの情報を得ることができます

お役に立てば幸いです...


ありがとう。いつ私はそれをするべきですか?今のところ、私は1.gdalinfo map.png 2.gdal_translate -of VRT -a_srs EPSG:4326 -gcp ....を実行しています。3. gdal2tiles.py -z 10-16 map.vrt。翻訳後に行う必要がありますか?変換方法でEPSGを変更する必要がありますか?
shannoga

マップが調整されていない場合、gdalinfoで情報を取得しようとしても無駄です。
アラゴン

代わりにgdalwarpを使用する必要がありますか?
shannoga

あなたの地図は投影と調整されていますか?はい、gdalwarpを使用して再投影しようとした場合、それ以外の場合はgdal_transformと調整します...
Aragon

患者さん、ありがとうございます。あなたが正しいと理解すれば違います。これは、特別なデータのない地図に基づいて描かれた単純な正方形の画像です。使用したコマンドで質問を更新しました。助けてくれてありがとう。
shannoga
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.