6

静止衛星の 3712x3712 ピクセル サイズの画像があります。画像が次のように見えるように、地球の周りに黒い部分があります。

eumetsat 画像ではない画像の例

地球の各ピクセルについて、緯度と経度を取得したいと思います。私は pyproj があることを知っており、次のようにプロジェクションをインスタンス化できました:

sat = pyproj.Proj('+proj=geos +lon_0 +h=035785831.0 +x_0=0 +y_0=0')

しかし、プロジェクションは私の画像の寸法（3712x3712）を認識していないため、ピクセルのラトロンを取得すること（画像内のピクセルの座標であるsat(x,y,inverse=True)wherexとare を使用）は明らかに不可能です。y

私は何が欠けていますか？

4

2 に答える 2

2

正しい投影ライブラリと設定を使用していると思います。

典型的なピクセル解像度 (1 ピクセルあたりの km) は、ユーメサットの Web サイトで報告されています。ここ では、1 ピクセルあたり約 3 km と報告されています。

子午線上の x、y への経度/緯度変換を実行し、これをピクセル数で割ることで確認できます (-81 度、81 度が最大範囲です。参照については、euetsat サイトを参照してください。 euetsat.int/ ):

import pyproj
sat = pyproj.Proj('+proj=geos +lon_0 +h=035785831.0 +x_0=0 +y_0=0')
x,y =  sat( 81.299, 0, radians = False, errcheck = True)
print (x * 2.0 / 3712.0 ) / 1000.0

eumetsat によって提供された情報に適合する 2.927 という値が得られます。

次に、マップのよく知られたポイントのセット (たとえば、海岸の特徴) を定義し、それらの x/y 位置を決定し、オンラインで緯度/経度座標を検索することにより、さらにキャリブレーションを行うことができます。さまざまなピクセル解像度を試して、どれが最適かを確認するか、より複雑なルーチンを使用することができます。

解像度は、赤道にどれだけ近いかによって異なる場合があります。こちらを参照してください)。したがって、いくつかの緯度で上記のルーチンを使用する必要がある場合があります。

于 2012-07-06T15:11:53.440 に答える
0

まず、これは EUMETSAT の画像ではありません。画像はオーストラリアを示しており、ゼロ度サービス (ZDS) とインド洋データ カバレッジ (IODC) の静止 EUMETSAT 衛星は (およそ) (0°N, 0°E) と (0°N, 41.5°E) の上にあります。 ) サービス、それぞれ。どちらの視点からもオーストラリアを見ることはできません。あなたの衛星画像はおそらく、2015 年に HIMAWARI に置き換えられるまで日本の気象庁 (JMA) によって運用されていた運輸多機能衛星 (MTSAT)イメージャーから来ています。

Meteosat SEVIRI、MTSAT、HIMAWARI AHI、またはその他の地球観測イメージャーから画像内のピクセルの位置情報を取得するには、ライブラリを使用してデータを読み取ることも検討できます。

from satpy import Scene
from glob import glob
filenames = glob('/path/to/your/files/*')  # for SEVIRI, can be HRIT, Native, or NetCDF
sc = Scene(filenames=filenames)
sc.load(["IR_108"], upper_right_corner="NE")
print(sc["IR_108"].attrs["area"].get_lonlat_from_array_coordinates(1000, 2000))

与える：

(-24.259016114751823, -3.9842388279769465)

静止画像の場合、これは最速の方法ではありません。答えが一定であるため、実際にデータを読み取る必要はないからです。ただし、十分にテストされ、広く使用されているライブラリを使用すると、エラーのリスクが軽減されます。少なくとも、参考にはなるでしょう。多くの衛星イメージャーでは、チャネルが異なれば空間解像度も異なるため、答えはチャネルに依存する可能性があることに注意してください。レベル 1 データは、マッピング規則に反して、画像の上部に南が表示され、画像の下部に北が表示されるように、南が北の前に来るように保存することもできます。したがって、結果が正確であることを確認する価値があります。

于 2022-01-14T09:08:39.640 に答える