imageData = imread('~/Documents/Large Melt Pond Images/Cropped/largempimage_2_0730_8.tif');
        BW = im2bw(imageData, .7);
        
        BW = BW(1:1000,1000:2000);
        
        %disp(BW);
        figure, imshow(BW);
        
        BW = imcomplement(BW); % water = 1, ice = 0
        
        [sz1,sz2] = size(BW);
        pxlcnt = sz1*sz2; %calculates # of pixels with dim of matrix
        propWater = sum(BW(:)==1)/length(BW(:));
        
        lk = bwlabel(BW,4); %labels groups of nearest neighbors

        lk(lk==0) = NaN;
        
        disp('propWater ='); disp(propWater);