pic = 'dar.jpg';
figure(1)
numbins=10;
H = imagehist3D(pic,numbins);
[I,J,K]= ind2sub([numbins,numbins,numbins],find(H>200));
[I,J,K]
colourfeature('dar.jpg',10,[9,8,2;9,7,2;9,7,1;9,6,1]);
pause
figure(1)
clf
figure(2)
clf