1

非常によく似た質問が同じ方法で解決されました:matplotlib.pyplot.imshowで「範囲」を使用する方法


地理的な軌跡を説明する地理座標 (「トラックログ」) のリストがあります。また、トラックログのカバレッジにまたがる画像を取得する手段があり、画像の隅の「地理的座標」がわかっています。

現在、私のプロットは次のようになっています (目盛りに注意してください - x=経度、y=緯度、UTM では WGS84):

ここに画像の説明を入力

次に、次の画像 (または青いトラックのないバージョン) のコーナー座標を知っていて、プロットの座標系に合うようにプロットしたいとします。

ここに画像の説明を入力

どうすればいいですか?

(ちなみに、気になる場合はタイルを使用する予定です)


ジョー・キングトンのコメント (私が受け入れることができるように彼の実際の回答を待っている) によると、次のコードは期待どおりに機能し、トラックログをプロットできるパン可能およびズーム可能な固定アスペクトの「ジオリファレンス」タイルを提供します。

import matplotlib.pyplot as plt
import Image
import numpy

imarray = numpy.asarray(Image.open('map.jpg'))

plt.plot([0,1], [0,1], 'o', c='red', ms=20)  ## some reference circles for debugging
plt.imshow(imarray, extent=[0,1,0,1])   ## some random map whose corners have known coordinates
plt.axis('equal')
plt.show()

ここに画像の説明を入力

4

1 に答える 1

0

ここでは実際にはあまり答えはありませんが、matplotlibを使用していて、geos-tuffを使用している場合は、matplotlib.basemapを参照してください
デフォルトでは、すべての操作はUTMマップで実行されますが、独自の投影法を選択できます。たとえば、 http://www.geophysique.beにある優れたチュートリアルのリストも参照してください。

于 2012-09-11T06:01:09.300 に答える