first try at histograms

            elif s[-6] == '3':

                tiff3.append(os.path.join(scene, s))

    cameras = [tiff0, tiff1, tiff2, tiff3]

    diff = []

    for ii in range(len(tiff1)):

        image = Image.open(tiff1[ii])    #Open the image and read it as an Image object

    i = 0

    for cam in cameras:

        for ii in range(len(cam)):

            image = Image.open(cam[ii])    #Open the image and read it as an Image object

            image = np.array(image)[1:,:]    #Convert to an array, leaving out the first row because the first row is just housekeeping data

            ar1, ar2 = image.shape

                else:

                    surrounding.append(image[ind1+i, ind1+j])               #Add the other 8 pixels to the list

            diff.append(np.max(surrounding)-np.min(surrounding))

        plt.subplot(2,2,i+1)

        plt.hist(diff)

        i += 1

    plt.tight_layout()

    plt.show()