問題タブ [map-projections]

画像のコーナー座標は既知です：

• 左上: (70.20315、20.92749)
• 右上: (70.20315, 31.45240)
• 左下: (59.70571, 20.92749)
• 右下: (59.70571, 31.45240)

(基本的に、フィンランドとスウェーデン/ノルウェー/ロシアの一部をマップします)

できれば Python を使用して、画像投影を Google マップのオーバーレイ画像に変更する方法が見つかりませんでした。

pyprojとPILでそれを行うことは可能ですが、その方法を理解できませんでした。GDAL、より具体的にはgdalwarpも見つけましたが、正しい方法で投影を変更するためのオプションが見つかりませんでした。

緯度/経度で定義された位置をマップ (画像) の x、y 座標に投影できる関数を作成する必要があります。

マップ自体は、投影タイプとして「ベルギーのランバート」を使用しているようです。

この投影法とその構成に関する公式情報をいくつか見つけました: http://ign.be/FR/FR2-1-4.shtmおよびhttp://ign.be/Common/Lambert2008/Transformation_Geographic_Lambert_FR.pdf (フランス語で、私は必要に応じて翻訳します）。

基本的に、いくつかのライブラリを使用する方が簡単だという結論に達しました (ここでは必須ではありませんが)。いくつかの調査の結果、OpenMap がその作業を実行できるように思われます。

だから、これが私がこれまでに得たものです：

（私のgifが1000x1000であることを知っています）

LamberConformal コンストラクターの JavaDoc: http://openmap.bbn.com/doc/api/com/bbn/openmap/proj/LambertConformal.html#LambertConformal%28com.bbn.openmap.LatLonPoint,%20float,%20int,%20int, %20double,%20double,%20double,%20double,%20double,%20double,%20com.bbn.openmap.proj.Ellipsoid%29

私はそれを正しく構成していなかったと思います: マップ上にあるはずのポイント xxx で、次の結果が得られます。

パラメータとして「中心点」（？）を与える

(マップの真ん中、いい感じ)

これに精通している、またはリンクから構成を把握できる人はいますか?

編集：

メソッドの最後の行をこれに置き換えてみました：

しかし、良くはありません。

別の質問を作成しました（現在地が移動する地図（静止画像ファイル）を表示する方法）

しかし、それには 2 つの質問が含まれているため、分離する必要があります。

人を目的地まで案内できるプロトタイプを作ろうとしています。

• 場所は数階建ての広い建物です。
• 地図（静止画）を取得・取得できます。例: 現在地:1F 目的地:5F; 1階、2階…5階の静止画（5枚の画像ファイル）を取得できます。

シナリオ：

1. アプリケーションを開始する
2. 現在地を入力 (または現在地を使用して自動的に設定される場合があります) & 目的地
3. 検索ルートボタンをクリックして、使用する地図（静止画）を検索し、現在地と目的地をマークします
4. 目的地への移動/移動時に現在地を更新する

問題: WiFi/セル タワー/IP アドレスを介して現在の位置座標を取得できますが、それを静止画像に入れて現在の位置をマークする方法がわかりません。

概念/アイデアを共有するか、コード スニペットを含めますか。私の論文に大いに役立ちます。

正しい方向へのガイダンスをいただければ幸いです。

アップデート

実際の、期待される出力、テスト ケースを含む実際の例 (ここにあるコードの一部は、mapsforge から取得したものです)

メルカトル投影クラス.java

上記のコードの出力:

• PixelXYToLatitude テスト結果:PASSED
• LatLngToPixelXY テスト結果:PASSED

LatLng を Pixel に変換するための投影クラスが既にあります。私の問題は、上記のクラスを使用して、特定の LatLng の静止画像にマークを付ける方法です。

これが私の静止画像です（シカゴ、イリノイ州）：

マークを付けたい（これは例ですが、後でバルーンを小さなポインターに変更する必要があります）

関連する緯度と経度を持つ不規則にグリッド化されたデータセット (生の衛星データ) を、 で指定された規則的にグリッド化された緯度と経度のセットにマッピングしようとしていますbasemap.makegrid()。装着して使用matplotlib.mlab.griddataしていmpl_toolkits.natgridます。以下は、whosipython で出力として使用されている変数のリストと、変数に関するいくつかの統計です。

datalatとdatalonは元のデータ座標です

gridlatとgridlonは補間する座標です

var実際のデータが含まれています

これらの変数を使用すると、呼び出しgriddata(datalon, datalat, var, gridlon, gridlat)が完了するまでに 20 分もかかり、nan. データを見ると、緯度と経度は正しいように見え、元の座標は新しい領域の一部と重なり、いくつかのデータ ポイントは新しい領域の外にあります。誰か提案はありますか？ナン値は、私が何か愚かなことをしていることを示唆しています...

ユーザーがシェープファイルをアップロードできるシステムを構築しています。次に、これらのシェープファイルを を使用して PostGIS に変換しshp2pgsqlます。このコマンドには、EPSG コード形式の SRS ID が必要です。

*.prjしたがって、シェープファイルのファイル (WKT としてエンコードされた投影/空間参照システムを含む) を読み取って、対応する SRS ID を返すことができる ruby​​ gem が必要です。

ピクセル位置を画面から度に変換する際に問題があります（EPSG4326）。

http://proj4js.org/とグーグルマップによると、EPSG900913のlon/lat値は正しいです。

私が得ているいくつかの例の値は

変換後

値が必要な場合

私が間違っていることについて何か考えはありますか？

これが私のコードの関連部分です。

このコードで 3D クラウド ポイントのコレクションを表示しました。

このコードは、3D ポイント クラウドの配列とその色の値を表示します。2D ポイントの配列を取得し、これらのポイントを 2D で表示したい (3D ポイントの上部ビュー)

OpenLayers を使用して geoDjango で GIS アプリを構築し、投影 EPSG:3908 で PostGIS のデータを使用しています。ベースレイヤーはOSMです。ベクター データを追加すると、ボスニアではなくブラジルの近くでレンダリングされました。私はすべてを試しましたが、解決策を得ることができません。

誰でも助けてください。

これらのサイトで2つの異なるメルカトル式が議論されている理由についての説明を探しています。

私はこれが正しいメルカトル図法アルゴリズムであることを理解しています。

http://en.wikipedia.org/wiki/Mercator_projection

ただし、このサイトは完全に異なるものを参照しています： http ://wiki.openstreetmap.org/wiki/Mercator

何か案は？

タイトルですでに述べたように、ミラー図法で緯度/経度をx/y座標に変換するのに助けが必要です。私は現在、都市リストを入力として取得し、YahooPlacefinderから各都市の緯度/経度を取得するアプリ（Java）を作成しています。私は自分のコードでこれらの式を使用しました。これが私が得たものの例です。（画像は参照用であり、私が使用しているものではありません）。X位置がわずか数ピクセル（2〜3）ずれていることがわかるように、これは、このマップの本初子午線シフト（CENTRAL_MERIDIAN_OFFSET）の計算で問題になる可能性があります。しかし、主な問題はY座標が正しくないことです。

これが私のコードです（更新-赤道オフセットの34px補正）：

明らかに、Y座標の計算に問題があります。5.6が本当に正しい値であるかどうかはわかりませんが、ミラー図法の垂直範囲は、私が読んだ参考文献の1つで-2.3 .. + 2.3と言われていたので、それを使用しました。