function [HILT] = imgbandfind(img,minb,maxb)
%IMGBANDFIND inputs an image array and minimum and maximum value,
% then finds all values of the array within that range, then plots with
%values in range highlighted in red against a black background

img=double(img);
maxv=max(max(img));
sizm=size(img);
rows=sizm(1,1);
cols=sizm(1,2);
yesm=zeros(rows,cols);
r=1;
c=1;
while r<(rows+1)
    while c<(cols+1)
        yesm(r,c)=double(~((minb<img(r,c))&&(img(r,c)<maxb)));
        img(r,c)=double(~(yesm(r,c)))*maxv;
        c=c+1;
    end
    c=1;
    r=r+1;
end
yesm=double(yesm);
HILT=cat(3,img,yesm.*img,yesm.*img);
imshow(HILT);
title(['Intensities between ' num2str(minb) ' and ' num2str(maxb)]);

end