gdalを使用してWebメルカトルタイルを正しくジオリファレンスする方法は?
例として、次のタイルhttp://a.tile.openstreetmap.org/3/4/2.pngを取得し、「4_2.png」として保存します。 このタイルのWGS84座標を計算するか、読み取ることができますが、対応するタイルをクリックして: 0 66.51326044311185 45 40.97989806962013 (West North East South) タイルを正しくジオリファレンスする方法(gdalを使用してジオティフまたは別のジオリファレンス形式を生成する): ビットマップを引き伸ばす必要はありません(=ジオティフのピクセルは元のビットマップのピクセルとまったく同じです) 結果の画像は、GISビューア/エディタの適切な場所で開かれますか(TatukGIS Free Viewerなど)? (2011年9月19日に編集して、質問を明確にし、結論を含めました) 私の結論: 私は最初に、3番目のアイデア(以下を参照)が正しいものであると考えました。GISビューアーでジオティフを開き、表示された座標を予想と比較しました。2番目のアイデアのジオティフは、2ピクセル北にシフトしているようです。だからこそ、私はアイデア3(または4)を正しいものと考えました。 しかし、はるかに高いズームレベルでタイルを試してみると、アイデア3からのジオティフは決定的に南にシフトします。ズームレベル3のタイル上の座標を比較するのはばかげていました。そのようなズームレベルでの国境は、比較が良い結果をもたらさないように単純化されています。 Dan S.は正しかった、タイル画像はすでにEPSG:3857にあります。2番目のアイデアは正しいアイデアです(高ズームレベルでも良い結果が得られます) 最初のアイデア:EPSG:4326 WGS84座標のEPSGコードはEPSG:4326です。そのため、単にgdal_translateを使用してWGS84座標を使用して、タイルをgeotifとしてジオリファレンスします。 gdal_translate -of Gtiff -co tfw=yes -a_ullr 0 66.51326044311185 45 40.97989806962013 -a_srs EPSG:4326 4_2.png t4326.tif 結果のマップは適切な場所に表示されますが、投影が正しくなく、タイルの中央にシフトがある可能性があります。gdalwarpを使用してマップを再投影して確認するために長い時間を試した後、Global Mapperのデモバージョンをダウンロードしました。EPSG:4326座標を使用できるように、画像をストレッチする必要があります。 2番目のアイデア:EPSG:3857 このタイルは、「web mercator」プロジェクション(エイリアスgoogleマッププロジェクション)を使用します。EPSG コードはEPSG:3857(エイリアスEPSG:900913)になりました。私は単にgdaltransformを使用して座標を変換します: gdaltransform -s_srs EPSG:4326 -t_srs EPSG:3857 0 66.51326044311185 …