ラベリングした画像の​ラベルごとの面積を求​めるにはどうしたらよ​いですか？

バイナリイメージの白の部分の面積を求め，一定の閾値より小さいものを削除するプログラムを作成するために，ラベリング処理をしてラベルごとの面積を求めようとしているのですが方法がわかりません．

https://jp.mathworks.com/help/images/labeling-and-measuring-objects-in-a-binary-image.html

上記のページを参考にしているのですが，このページだと例えばラベル1は6ピクセル，ラベル2は6ピクセル，ラベル3だけ5ピクセルなのでラベル3の部分を0に置き換えるようなプログラムを作ろうとしています．

Answer 1

[Tohru Kikawada]

ある面積以下の領域を0に置き換えたい場合は bwareaopen が使えそうです。

それとも特定の面積の領域を0にしたいのでしょうか。その場合には bwarea が使えそうですね。

領域解析を詳細にやりたい場合は regionprops も使えるかもしれません。ご参考まで。

[Yumi Iwakami]

2値画像である面積以下の領域を0に置き換えるとともに，置き換えた元の領域の元画像（元画像はグレースケール）のピクセルの輝度値の合計もとりたいです．

[Tohru Kikawada]

regionprops のヒントの項目に置き換える方法の記載がありますよ。また、 accumarray を活用すると輝度値の合計もシンプルに求められます。

%%元画像の読み込み(グレイスケール)
G = imread('coins.png');
figure, imshow(G);

%%2値化
BW = imbinarize(G);
BW2 = imfill(BW,'holes');
figure, imshow(BW2);

%%面積が2000ピクセル以下を抽出
CC = bwconncomp(BW2); 
stats = regionprops('table',CC,'Area','Centroid');
figure, histogram(stats.Area,10);
idx = stats.Area <= 2000

%%2000ピクセル以下の領域を0にする
% 2000ピクセルを超える領域だけを抽出する
BW3 = ismember(labelmatrix(CC), find(~idx));
figure, imshowpair(BW2,BW3,'montage');

%%2000ピクセル以下の領域の輝度値の和

% 領域ごとの輝度値の和を求める
L = labelmatrix(CC)+1; % ラベルを1はじまりにするため+1
A = accumarray(L(:),double(G(:))); % ラベルごとの累積値計算
A = A(2:end); % 0に相当する項目(黒)を除去

% 面積が2000ピクセル以下の累積値を可視化
Iout = G;
BW4 = ismember(labelmatrix(CC), find(idx));
Iout(~BW4) = 0; % 該当箇所以外すべて0にする
Iout = insertText(Iout,stats.Centroid(idx,:),cellstr(num2str(A(idx))));
figure, imshow(Iout);

[Yumi Iwakami]

ご丁寧な解説ありがとうございます． やってみます．