基于Matlab的模糊C均值肿瘤分割

基于Matlab的模糊C均值肿瘤分割,matlab,image-processing,Matlab,Image Processing,我有一个分段的肝脏。我需要把肿瘤分割开来。我使用FCM方法。这是一个3级FCM阈值。当我将其应用于图像时,我需要仅对肿瘤区域(比其余部分暗的区域)进行分割。但我的情况正好相反。肿瘤周围的所有区域都被分割。请帮帮我。该程序有两个文件,testfcmthresh.m和一个函数fcmthresh.m 输入“分割肝脏(使用区域生长)”和FCM输出图像: 我尝试使用imcomplete()来补充图像,但我得到的整个背景也是白色的,因为背景原本是黑色的。请帮帮我 Ghaul在我的上一个问题中已经给出

我有一个分段的肝脏。我需要把肿瘤分割开来。我使用FCM方法。这是一个3级FCM阈值。当我将其应用于图像时,我需要仅对肿瘤区域(比其余部分暗的区域)进行分割。但我的情况正好相反。肿瘤周围的所有区域都被分割。请帮帮我。该程序有两个文件,
testfcmthresh.m
和一个函数
fcmthresh.m

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

我尝试使用
imcomplete()
来补充图像,但我得到的整个背景也是白色的,因为背景原本是黑色的。请帮帮我


Ghaul在我的上一个问题中已经给出了我问题的答案。如果有人需要参考,请仔细阅读Ghaul在这方面的评论。谢谢。

我们可以有肿瘤和肝脏的图片吗?我是新用户。所以我不能在这里张贴图片。如果你给我你的身份证,我会把它寄给那个。提前谢谢。另外,我的项目正在使用Matlab将它们上传到某个地方,并将链接放入问题中。我将在中编辑图像。我已编辑了我的问题(添加了图像链接)。将有2个图像。输入的“分割肝脏图像”和最终的“fcm输出”图像。请检查一下。反转图像怎么样?
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');