global cover_object; global watermark_en; global CWI; % 将载体图像进行DCT变换 dct_carrier = dct2(cover_object); % 对水印图像进行二值化处理 bw_watermark = imbinarize(watermark_en); % 将二值化的水印图像进行重复，以达到和载体图像相同的大小 bw_watermark = repmat(bw_watermark,[size(carrier,1)/size(bw_watermark,1),size(carrier,2)/size(bw_watermark,2)]); % 将二值化后的水印图像转换为1和-1的矩阵 bw_watermark(bw_watermark==0) = -1; % 获取DCT变换后的载体图像的大小 dct_size = size(dct_carrier); % 循环遍历DCT变换后的载体图像的每一个8x8的块 for i=1:8:dct_size(1)-8 for j=1:8:dct_size(2)-8 % 获取当前8x8块的DCT系数 dct_block = dct_carrier(i:i+7,j:j+7); % 获取当前8x8块的最低频率DCT系数，并将其用于嵌入水印 lowest_dct = dct_block(1,1); if bw_watermark((i-1)/8+1,(j-1)/8+1) == 1 dct_block(1,1) = lowest_dct + 10; else dct_block(1,1) = lowest_dct - 10; end % 将修改后的8x8块的DCT系数放回载体图像 dct_carrier(i:i+7,j:j+7) = dct_block; end end % 将修改后的DCT变换后的载体图像进行反DCT变换 watermarked_carrier = idct2(dct_carrier); CWI=uint8(watermarked_carrier); axes(handles.axes2); imshow(CWI); title('含水印信息的载体图像'); global CWI;

这段代码是一个简单的数字水印嵌入算法。它通过将水印图像二值化和重复，使其与载体图像具有相同的大小。然后将二值化后的水印图像转换为1和-1的矩阵。接下来，该算法循环遍历DCT变换后的载体图像的每一个8x8的块，并获取当前8x8块的DCT系数。然后，它获取当前8x8块的最低频率DCT系数，并将其用于嵌入水印。如果当前像素为1，则将最低频率DCT系数加上10，否则将其减去10。最后，该算法将修改后的DCT变换后的载体图像进行反DCT变换，以得到嵌入水印后的载体图像。代码中的变量 cover_object 和 watermark_en 分别代表载体图像和水印图像，而 CWI 则是嵌入了水印的载体图像。