python环境下的遥感图像dct水印提取

插入在前面的文章 本文主要是提取算法，详细注释看插入算法

import gdal

import numpy

from skimage import data,img_as_float,img_as_ubyte

import cv2

filename = "f:/bs/band_11_255_0.0075.tif"#插入水印后的图片

dataset = gdal.Open(filename)

im_width = dataset.RasterXSize #获取长

im_height = dataset.RasterYSize #获取宽

im_bands = dataset.RasterCount #获取通道数

im_data = dataset.ReadAsArray(0,0,im_width,im_height)

BLOCK = numpy.zeros([32,32])

M = 3808

N = 238

K = 16

I = numpy.zeros([M,M])

I = im_data[2000:,2000:]

filename = "f:/bs/band_1.tif"

dataset = gdal.Open(filename)

im_width = da