findpeaks関​数を使って暗い画像の​中に存在する広がりの​ある輝点を検出したい​のですが、方法はあり​ますか？

今直面している問題は、例えば列方向にfindpeaksをかけると、複数の列にわたって同じ輝点が列ごとに検出されることと、それが同じ輝点であることが判別しにくい点です。

Ad-onとして2Dピーク検出のコードがあることは知っていますが、できればfindpeaksで対応したいです。

[Akira Agata]

findpeaks 関数を使いたい理由が何かあるということでしょうか？

輝点を検出したいということであれば、ある輝度値を閾値として画像を２値化するのがいちばん単純な方法になるかと思います。

[ヨナ エマヤ]

ご回答ありがとうございます。

質問の本文に書くべきでしたが、２値化で対応するプログラムは既に試していて、現在は違うアプローチを模索しています。

他のところでIslocalmaxを行方向と列方向にそれぞれ走らせて論理和をとればうまくいくのでは、という方法を教えていただいたのですが、結果は良さそうです。同じように論理和で対応する方法をfindpeaksでも試してみたら、意外とfindpeaksの方が速かったです。

しかし、これだけで輝点が抜けたり、２重検出されたりするところはまだ残っているので、現在はどのように改善できるかを考えています。

Answer 1

[Akira Agata]

試しに簡単なプログラムを作成してみました。行方向と列方向でピークを検出してプロットすると、検出されたピーク点は輝点の上に十字型に表示されます（右図）。ただ、輝点内の輝度分布によっては十字の縦と横の線がきれいなラインにならない場合や、中央でクロスしない場合などがあり、ここから輝点の中心を検出するのは難しい印象です。

もし @ヨナ エマヤ 様が想定されている対象画像を１，２枚共有頂けると、別の方法もご提案可能かと思いますがいかがでしょうか？

% 輝点のある画像を作成
I = imread('coins.png');
I = imgaussfilt(I,10);

% findpeaksで行,列それぞれの方向のピークを検出
J = false(size(I));
for kk = 1:size(I,2)
  [~, pt] = findpeaks(double(I(:,kk)),...
    'MinPeakProminence',50);
  J(pt,kk) = true;
end
for kk = 1:size(I,1)
  [~, pt] = findpeaks(double(I(kk,:)),...
    'MinPeakProminence',50);
  J(kk,pt) = true;
end

% 元画像と検出した行,列各方向のピーク位置を表示
imshowpair(I,J,'montage')

[ヨナ エマヤ]

サンプルプログラムを書いていただきありがとうございます。

こちらのサンプルデータですが、使っている画像については研究の都合上、提供できません。せっかくご対応いただいたのに、申し訳ございません。

@Akira Agata様のコードで作成した輝点はこちらで使っている画像に類似しているので、試しに使わせていただきました。以下のようなプログラムを掛けてみたところ、いい結果が求まったかと思います。

I = imread('coins.png');
I = imgaussfilt(I,10);

TF = islocalmax(I,1) & islocalmax(I,2);

J = zeros(size(I));
J(TF) = 1;

imshowpair(I,J,'montage')