オブジェクトのトリミング
显示 更早的评论
添付した図のように、複数の細胞が存在するバイナリイメージ(二値画像)から、細胞を1つ1つトリミングしたいと考えています。
何かいい案を頂けると幸いです。
※以前まで、個々の細胞の重心を求め、そこを中心として必要なピクセルにトリミングしようと考えていたのですが、
全体のバイナリイメージから個々の細胞の重心を求める方法が分からず断念致しました。
2 个评论
Akira Agata
2021-9-16
バイナリーイメージ内で細胞が１つ１つ分離しているものについては、regionprops 関数で Bounding Box を計算することでトリミング (クロップ) 可能です。また、同じ関数により個々の細胞の重心も求めることができます。ただ、複数の細胞が融合して "かたまり" になっている部分では、そのような処理は困難です。
まずは細胞が１つ１つ分離しているものについてのみ計算する、というアプローチでも問題ないでしょうか？
HIGUCHI 樋口
2021-9-23
采纳的回答
各位のコメント内容を実際に動かしてみました。尚、image processing toolboxが必要です。
I = imread('saibou.png');
BW = imbinarize(rgb2gray(I)); % カラー⇒グレースケール⇒二値化
BW = imopen(BW,strel('disk',10)); % モルフォロジーオープニングで小さなドットを除く
% Bounding Box, 重心, 面積を得る
stats = regionprops('table', BW, {'Centroid', 'BoundingBox', 'Area'});
idx = stats.Area > 10000;
stats(idx,:) = []; % 面積が10000ピクセル以上の領域は以降の処理から除外する
BB = stats.BoundingBox + [-40 -40 80 80]; % Bounding Boxに余白を付ける
imshow(I); hold on; % 加工前画像。小さなドットや接触した細胞は認識されない
plot(stats.Centroid(:,1),stats.Centroid(:,2),'b*'); % 重心を表示
for i = 1:size(BB,1)
rectangle('Position', BB(i,:), 'EdgeColor', 'y'); % Bounding Boxを表示
end
hold off;
for i = 1:size(BB,1)
crop{i} = imcrop(BW, BB(i,:)); % 細胞を1つ1つトリミングする(一部はくっついたまま)
if any([crop{i}(1,:) crop{i}(end,:) crop{i}(:,1)' crop{i}(:,end)'])
crop{i} = ~crop{i}; % Bounding Boxの境界に他の細胞が写り込んだら反転表示
end
end
montage(crop); % トリミングした二値化画像を並べて表示
4 个评论
HIGUCHI 樋口
2021-9-30
ありがとうございます。参考にさせて頂きました。
加えてもう一点質問がございます。可能でしたらご回答いただけると幸いです。
細胞個々のトリミングの前に、bwareaopen関数を用いてノイズのような非常に小さなオブジェクトを削除しました。
更に、細胞1つ1つを抽出するために、図のように細胞が重なっている大きな塊(例えば、一定の面積より大きいオブジェクト)を除去したいと考えていますが、難航しております。
何か良いアイディアはありますでしょうか。
Akira Agata
2021-9-30
regionprops を適用する際にあわせて面積も計算しておいて、ある値以上のものはその後の計算対象から除外するという方法で対応可能です。
stats = regionprops('table', BW, {'Centroid', 'BoundingBox', 'Area'}); % Bounding Box, 重心, 面積を得る
idx = stats.Area > 10000;
stats(idx,:) = [];
Atsushi Ueno
2021-9-30
なるほど。早速回答に反映しました。
HIGUCHI 樋口
2021-10-10
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Image Processing Toolbox 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
