Current Code 5/17/2013

up to date code below:

clc
close all
clear all

i = imread('Image1.tif');

j = i(:, :, 1:3);
%red = i(:,:,1);
clearvars i;
x = rgb2gray(j);
grad = edge(x, 'canny');
gradient = grad;

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

%fun = @(block_struct)block_struct.data;


r = imfilter(f,d);
%figure, imshow(r*4);
clearvars f;
clearvars y;
clearvars d;
clearvars x;
clearvars j;
clearvars grad;
b = imresize(4*r,0.5); %numbers for the filtered image
c = imresize(gradient,0.5); %numbers for gradient.

rs1 = reshape(b,size(b,1)*size(b,2),1); %creates an m by n matrix of size
rs2 = reshape(c, size(c,1)*size(c,2),1);

list = [rs1,max(rs1)*rs2];
list = double(list);

km = kmeans(list,2,'emptyaction', 'singleton');
SizeKM = size(km);
count = 1;
for i = 1:SizeKM
    if km(i) ~= 2
        storage(1,count) = i;
        count = count + 1;
    end
end