1.一种基于QR码的瓦片数据数字水印生成嵌入方法，包括以下步骤：步骤一、水印信息生成

步骤1.1：读取一段有确定含义的文字信息；

步骤1.2：根据二维码生成规则将读取的文字信息转换成版本三的二维码，并进行二维码预处理；

步骤1.3：按照编码方式将二维码转化为二进制的二维数组W0＝{0,1}，转换得到的数组长度为L＝Length*Length；

步骤二、水印嵌入

步骤2.1：利用GDAL以文件夹形式读取瓦片遥感图像数据；

步骤2.2：嵌入标识信息

对瓦片遥感图像数据进行16*16分块，在瓦片数据左上角、左下角、右下角的16*16分块嵌入固定的水印标识信息；

步骤2.3：对其余各16*16分块进行DCT得到该块的低频值Dij，计算量化值λij＝round(Dij/ξ)

其中，round运算为舍入取整操作，ξ为量化步长；

步骤2.4:水印嵌入，利用水印信息按照下式调制当前嵌入块的低频值：式中：Wij是嵌入的水印信息；

步骤2.5：利用调制后的低频值进行DCT逆变换，所有16*16分块逆变换完成后，得到含水印的瓦片遥感数据。

2.根据权利要求1所述一种基于QR码的瓦片数据数字水印生成嵌入方法，其特征是：所述二维码预处理是：去除二维码空白区、位置探测图形、定位图形和校正图形部分的水印信息。

3.根据权利要求1或2所述一种基于QR码的瓦片数据数字水印生成嵌入方法，其特征是：步骤2.2所述在瓦片数据左上角、左下角、右下角的16*16分块嵌入固定的水印标识信息，其嵌入规则为：左上角、左下角、右下角的16*16分块的直流分量嵌入水印信号“1”，比较中频系数[2,3]和[4,2]的大小，使得每个分块的DCT中频系数中[2,3]大于[4,2]。

4.一种权利要求1所述一种基于QR码的瓦片数据数字水印生成嵌入方法生成嵌入的基于QR码的瓦片数据数字水印的提取方法，其步骤是：步骤A：利用GDAL以文件夹形式读取瓦片遥感图像；

步骤B：标识信息提取

选择第一个256*256分块，判断左上角、左下角、右下角提出的水印标识信息是否有误，即每个块的直流分量提取出的水印信息是否为“1”，及每个块中频系数[2,3]是否大于[4,

2]；

如果每个块的直流分量提取出的水印信息为“1”，则进入步骤C，否则，以遍历的方式计算下一个256*256分块；

步骤C：计算该256*256分块中各16*16分块的DCT变换后的低频值D′ij；

步骤D：水印提取，计算量化值λij＝floor(Dij/ξ)。其中，floor为取整操作；

式中：Wijˊ是提取的水印信息；

步骤E：处理完所有16*16分块，利用ThoughtWork.QRCode.dll处理得到的0-1数据生成二维码，解码该二维码；

若正确提取出水印信息则输出水印信息，

若无法解码信息则提取出的水印信息错误，转向步骤B继续以遍历的方式计算下一个

256*256分块。