Current Code

current code. 

clc
close all
clear all

i = imread('Image1.tif');

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

y = fspecial('laplacian', 0.9);
d = fspecial('gaussian', [3 3], 500);
f = imfilter(x,y);

r = imfilter(f,d);
figure, imshow(r*4);
clearvars f;

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);