ラスタイメージの領域を別のイメージの値で置き換えます

Question

同じ領域の2つのラスタイメージとnumpy配列のx、y次元を持っています。画像1は土地利用分類（例えばクラス0~5）、画像2は雲影マスク（値は0 =雲霧、255 =雲/陰影領域）である。
これらの画像を結合したいと思います。画像2の255個の値をすべて取ってクリップし、画像1にモザイクします。または、画像2のすべての0の値を画像1の特定のピクセル位置にある値に置き換えます。

2d配列1dと0の値を置き換えますが、2dに正しく変換できませんでした。

このラスタ計算を完全にオープンソースのpythonで行うには、どのような方法が一番簡単ですか？

Answer 1

これは、numpyのブール型インデックス機能で実行できます。

img1 = np.array([[0, 1, 0, 1],[1, 0, 1, 0]]) 
img2 = np.array([[1, 2, 3, 4],[5, 6, 7, 8]]) 
bool_arr = img1 == 0 
img1[bool_arr] = img2[bool_arr] 
print(img1) 
# outputs: [[1 1 3 1] 
#   [1 6 1 8]]