土地被覆ごとのプロット数を数える

Question

前提

色分けした図と散布図を重ね合わせて、それぞれの色の上に存在するプロット数をカウントしようと考えています（どの土地被覆で一番プロットが多いかを調査）。
プロットの数ではなく、地図の画素の方を数えてしまっています(エラーなし)。
散布図を土地被覆に紐づけて区分する点で詰まっています。
2つのレイヤ間は同じ解像度で位置にずれはないと考えられます。
なにかお気づきの点がありましたらご教示いただけますと幸いです。

実現したいこと

土地被覆ごとのプロットの数を明らかにしたい。

試したこと

###画像の色分離
land_color = np.array(img)#画像の色をland_colorに格納
r=land_color[:,:,0]#R
g=land_color[:,:,1]#G
b=land_color[:,:,2]#B

land_color[np.where(r==194)]=A#grayを抽出してAとラベリング
land_color[np.where(r==140)]=B#brownを抽出してBとラベリング
land_color[np.where(r==61)]=C#greenを抽出してCとラベリング

###散布図
fig = plt.figure(figsize = (5,5))
ax_data = fig.add_axes((1, 1, 1, 1))
ax_data.scatter(x_list, y_list, 'o', color = 'black')

###カウント
sum_A=len(x_list[A])#(y座標でも検証済)
sum_B=len(x_list[B])
sum_C=len(x_list[C])

Answer 1

コードが断片的かつx_listなどの変数の型ももろもろ不明なので、基本的なコード例を示します。

測定点の位置リスト？から「測定点のみが1である行列」を作成します。
あとはこの行列と「特定の色である部分のみが1(True)となる行列」と掛け合わせると
その色上にある測定点のみ１となるのでそれを合計すればよいです。

import numpy as np

# 色分けされた図
land = np.array([
    [1,1,1,1,1],    # AAAAA
    [2,2,2,1,1],    # BBBAA
    [3,3,2,2,2],    # CCBBB
    [3,3,3,3,3]])   # CCCCC
A,B,C = [1,2,3]

# 測定点の位置リスト
x_list = [0,4,2,3,4,2]
y_list = [0,0,1,1,2,3]
# x___x
# __xx_
# ____x
# __x__

# 測定点のみが1である行列を作成
plot = np.zeros(shape=land.shape, dtype=int)
plot[y_list,x_list] = 1
print(plot)
#[[1 0 0 0 1]
# [0 0 1 1 0]
# [0 0 0 0 1]
# [0 0 1 0 0]]

# 「(land==A)」＝Aである部分のみが1(True)となる行列 と
# 測定点が1である行列を掛け合わせた結果
print( (land==A) * plot )
#[[1 0 0 0 1]
# [0 0 0 1 0]
# [0 0 0 0 0]
# [0 0 0 0 0]]

# 重ね合わせて残った部分の「1」を合計
print( np.sum( (land==A) * plot ) ) # 3
print( np.sum( (land==B) * plot ) ) # 2
print( np.sum( (land==C) * plot ) ) # 1

Comments

[reado]

ありがとうございます！
参考にさせていただきます。