matlab - 使用 Matlab 进行肿瘤分割的模糊 C 均值

Question

我有一个分段的肝脏。我需要分割里面的肿瘤。我使用了 FCM 方法。这是一个 3 级 FCM 阈值。当我将它应用于图像时，我需要单独对肿瘤区域(比其余部分更暗的区域)进行分割。但我得到了相反的结果。肿瘤周围的所有区域都被分割。请帮帮我。该程序有两个文件，testfcmthresh.m 和一个函数fcmthresh.m

输入“分割肝脏(使用区域生长)”和 FCM 输出图像:

from here

我尝试对使用 imcomplement() 获得的图像进行补充，但我得到的整个背景也是白色的，因为背景最初是黑暗的。请帮帮我。

function [bw,level]=fcmthresh(IM,sw)
%FCMTHRESH Thresholding by 3-class fuzzy c-means clustering
%  [bw,level]=fcmthresh(IM,sw) outputs the binary image bw and threshold level of
%  image IM using a 3-class fuzzy c-means clustering. It often works better
%  than Otsu's methold which outputs larger or smaller threshold on
%  fluorescence images.
%  sw is 0 or 1, a switch of cut-off position.
%  sw=0, cut between the small and middle class
%  sw=1, cut between the middle and large class
%
%  Contributed by Guanglei Xiong (xgl99@mails.tsinghua.edu.cn)
%  at Tsinghua University, Beijing, China.

% check the parameters
if (nargin<1)
    error('You must provide an image.');
elseif (nargin==1)
    sw=0;
elseif (sw~=0 && sw~=1)
    error('sw must be 0 or 1.');
end

data=reshape(IM,[],1);
[center,member]=fcm(data,3);
[center,cidx]=sort(center);
member=member';
member=member(:,cidx);
[maxmember,label]=max(member,[],2);
if sw==0
    level=(max(data(label==1))+min(data(label==2)))/2;
else
    level=(max(data(label==2))+min(data(label==3)))/2;
end
bw=im2bw(IM,level);

---

%testfcmthresh.m

clear;clc;
im=imread('mliver3.jpg');
fim=mat2gray(im);
level=graythresh(fim);
bwfim=im2bw(fim,0.1);
[bwfim0,level0]=fcmthresh(fim,0);
[bwfim1,level1]=fcmthresh(fim,1);
subplot(2,2,1);
imshow(fim);title('Original');
subplot(2,2,2);
imshow(bwfim);title(sprintf('Otsu,level=%f',level));
subplot(2,2,3);
imshow(bwfim0);title(sprintf('FCM0,level=%f',level0));
subplot(2,2,4);
imshow(bwfim1);title(sprintf('FCM1,level=%f',level1));
% imwrite(bwfim1,'fliver6.jpg');

Answer 1

Ghaul 在我之前的问题中已经回答了我的问题 'Extracting image region within boundary' .如果有人需要引用，请仔细阅读 Ghaul 的评论。谢谢。