用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的评论。谢谢。