Code updated, K-Means

clc
close all
clear all

i = imread('Image1.tif');

j = i(:, :, 1:3);
red = i(:,:,1);
clearvars i;
x = rgb2gray(j);
%m = x;

%k = fspecial('sobel');
%g = imfilter(m,k);
y = fspecial('laplacian', 0.9);
d = fspecial('gaussian', [3 3], 500);
f = imfilter(x,y);
%m1 = reshape(f, 8702*6002, 1);
r = imfilter(f,d);
%m2 = reshape(r, 8702*6002, 1);
clearvars f;
%figure, imshow(r*3);
%BW = im2bw(r*3, 0.11);
b = imresize(r, 0.5);
r2 = imresize(red,0.5);
clearvars r;

m3 = reshape(b, size(b,1)*size(b,2),1);
m4 = reshape(r2, size(r2,1)*size(r2,2),1);
m = [m3,m4];
m = double(m);
%figure,imshow(BW);
clearvars BW;
clearvars d;
clearvars j;
clearvars x;
clearvars y;
km = kmeans(m, 3);
plot(km);
%t = kmeans(z, 3);
%figure,imshow(t);
%z = imadjust(x,[0.0 0.3],[0 1]);
%figure, imshow(x);
%figure, imshow(z);

%bw1 = edge(z,'sobel');
%bw2 = edge(z, 'canny');
%figure,imshow(bw1);
%figure,imshow(bw2);