海岸線に沿って Sentinel-3 シーンをプロットするためのPyresample の例に従っています。ただし、プロットしようとしている 3D (3 バンド) Sentinel シーンではなく、2D 配列をデータとして期待しているようpyresample.kd_tree.resample_nearest
に見えます。pyresample.save_quicklook
<xarray.DataArray 'reshape-fb54fb4cbc6c65d43122d7b5e1078a83' (bands: 3, y: 352, x: 507)>
dask.array<reshape, shape=(3, 352, 507), dtype=float64, chunksize=(3, 352, 507), chunktype=numpy.ndarray>
Coordinates:
* x (x) float64 -7.601e+04 -7.571e+04 ... 7.571e+04 7.601e+04
* y (y) float64 5.227e+04 5.197e+04 5.167e+04 ... -5.279e+04 -5.309e+04
* bands (bands) <U1 'R' 'G' 'B'
crs object PROJCRS["unknown",BASEGEOGCRS["unknown",DATUM["World Geod...
Attributes: (12/15)
resolution: 300
ancillary_variables: []
start_time: 2020-04-22 13:51:38
area: Area ID: area_id\nDescription: An Area of Interest\n...
end_time: 2020-04-22 13:54:38
sensor: olci
... ...
name: ocean_color
_satpy_id: DataID(name='ocean_color', resolution=300)
standard_name: ocean_color
prerequisites: [DataQuery(name='Oa08', modifiers=('effective_so...
optional_prerequisites: []
mode: RGB
疑似コードとエラーは次のとおりです。
lons, lats = scn_aoi[composite].area.get_lonlats()
composite_data = scn_aoi[composite].data.compute()
swath_def = pyresample.SwathDefinition(lons, lats)
roi = resolution - (resolution * 0.35)
# Below gives ValueError -> 'Mismatch between geometry and dataset'
composite_swath = (pyresample.kd_tree
.resample_nearest(swath_def, composite_data, aoi_def,
radius_of_influence=roi,
fill_value=None))
これらのデータはどのようにプロットできますか?