入力:
tws = ds1.lwe_thickness
print(tws.size)
print(tws.shape)
# tws
# vectorization
stacked1 = tws.stack(space =("latitude", "longitude")) #
stacked
print(stacked1.shape)
# standardization
Y = (stacked1-stacked1.mean())/(stacked1.std())
print(Y)
print(Y.shape)
print(Xt.shape) # precipitation dataset
print(Y.shape) # total water storage dataset
# co-variance between precipitation and total water storage
Cxy = (Xt*Y)/(72)
print(Cxy.shape)
print(Cxy)
print(Cxy.dims)
出力:
(64800, 72) (72, 64800) (64800, 0) <xarray.DataArray (スペース: 64800, 時間: 0)> array([], shape=(64800, 0), dtype=float64) 座標: * 時間(時刻) datetime64[ns] * スペース (スペース) マルチインデックス
- 緯度 (スペース) float64 -90.0 -90.0 -90.0 -90.0 ... 89.0 89.0 89.0 89.0
- 経度 (スペース) float64 0.0 1.0 2.0 3.0 4.0 ... 356.0 357.0 358.0 359.0 ('スペース', '時間')
最初の Xt マトリックスのサイズは 64800 72 で、Y は 72 64800です。これら 2 つのマトリックスの乗算後、サイズは 64800 64800 になるはずですが、64800 0 になっています。なぜこれが起こっているのでしょうか? この問題を解決するのを手伝ってください。基本的に、ここでは xarray を使用しています。
画像はここに添付されています