def METHOD_NAME(self):
    for band in range( self.res_ds.RasterCount ):
        band += 1
        if band:
            exp_band = self.exp_ds.GetRasterBand(band)
            res_band = self.res_ds.GetRasterBand(band)
            array1 = np.array(exp_band.ReadAsArray()).flatten()
            array2 = np.array(res_band.ReadAsArray()).flatten()
            regress_result = linregress(array1,array2)
            self.assertGreaterEqual(regress_result.rvalue, 0.9