3

pyresample モジュールを使用して、ASCAT 海風ベクトルとその wvc 品質フラグの完全なスワス軌道をプロットしようとしています。netcdf ファイルの ASCAT へのリンクは次の場所にあります : ftp://podaac-ftp.jpl.nasa.gov/allData/ascat/preview/L2/metop_a/12km/形状領域の定義を満たすためにファイル内の情報を見つける方法については、あまり詳しく説明されていません。衛星データのスワスをプロットするコードの例を以下に示します

from netCDF4 import Dataset     
import numpy as np
from pyresample import image, geometry
import pyresample as pr

netcdf ファイルから lons、lats、および wvc_quality_flag を抽出します

area_id = 'ease_sh'
name = 'Antarctic EASE grid'
proj_id = 'ease_sh'
proj4_args = 'proj=laea, lat_0=-90, lon_0=0, a=6371228.0, units=m'
x_size = 425
y_size = 425
area_extent = (-5326849.0625,-5326849.0625,5326849.0625,5326849.0625)
proj_dict = {'a': '6371228.0', 'units': 'm', 'lon_0': '0',
          'proj': 'laea', 'lat_0': '-90'}
area_def = geometry.AreaDefinition(area_id, name, proj_id, proj_dict, x_size,y_size, area_extent)
swath_def = geometry.SwathDefinition(lons=lon, lats=lat)
result = pr.kd_tree.resample_nearest(swath_def, wvc_quality_flag, area_def, radius_of_influence=20000, fill_value=None)
pr.plot.save_quicklook('/tmp/tb37v_pc.png', area_def, result, num_meridians=0, num_parallels=0, label='Flags')

AttributeError: 'module' object has no attribute 'plot'

最初に、pyresample モジュールにはその属性 plot がなく、ASCAT netcdf ファイルのどこにも「area_def」が定義されていないというエラーが表示されます。pyresample は、これらのタイプのファイルには妥当ではありませんか、それとも ASCAT ファイルのメタデータ内でこれらの定義の適切な場所を探していませんか? このモジュールに関するいくつかの説明は本当に役に立ちます! 再度、感謝します！

4

1 に答える 1

0

ここで実際に 2 つの質問があります。まずはplotサブモジュール。

plotはサブモジュールであり、他のサブモジュールと一緒にここにインポートできます。

from pyresample import image, geometry, plot

次に、ターゲット領域の定義は、このブロックでオンザフライで定義されます。

area_id = 'ease_sh'
name = 'Antarctic EASE grid'
proj_id = 'ease_sh'
proj4_args = 'proj=laea, lat_0=-90, lon_0=0, a=6371228.0, units=m'
x_size = 425
y_size = 425
area_extent = (-5326849.0625,-5326849.0625,5326849.0625,5326849.0625)
proj_dict = {'a': '6371228.0',
             'units': 'm',
             'lon_0': '0',
             'proj': 'laea',
             'lat_0': '-90'}

area_def = geometry.AreaDefinition(area_id,
                                   name,
                                   proj_id,
                                   proj_dict,
                                   x_size,
                                   y_size,
                                   area_extent)

ソース スワス定義は次のように定義されます。

swath_def = geometry.SwathDefinition(lons=lon, lats=lat)

lonそして、リサンプリングの前にとlat配列が取得されていると仮定します。

PS。はい、pyresample はまさにこの種のタスクを対象としています。swath データを含むファイルから余分なメタ データは必要ありません。経度と緯度の配列で十分です。

于 2016-11-29T14:45:01.027 に答える