histograms

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

    cameras = [tiff0, tiff1, tiff2, tiff3]

    diff = []

    i = 0

    for cam in cameras:

    fig, axs = plt.subplots(nrows=2, ncols=2, figsize=(15,12))

    for cam, ax in zip(cameras, axs.ravel()):

        diff = []

        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

                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)

        ax.hist(diff)

        i += 1

    plt.tight_layout()

    plt.show()