2025-12-25 来源:MiniSAR
1)生成128bit唯一版权标识(包含生产者ID、采集时间、数据编号);(3)相位域嵌入:
2)采用Logistic混沌映射(参数μ=3.99,初始值x₀=0.31415)对版权标识进行置乱,生成混沌序列;
3)通过AES-256算法对置乱后的水印加密,密钥由主密钥与数据采集时间戳生成;
1)对SAR复数影像的相位信息进行分块处理(块大小8×8);(4)数据重构:将嵌入水印后的相位信息与原始幅度信息重组,生成带水印的SAR复数影像。
2)采用最小相位扰动算法,在每个相位块的冗余位中嵌入1bit加密水印,嵌入强度根据相位敏感度动态调整(边缘区域嵌入强度≤0.02,平坦区域≤0.05);
import numpy as np
from cryptography.fernet import Fernet
import chaosmap as cm
# 1. 加载SAR复数影像(I/Q分量)
sar_data = np.load('mini_sar_complex.npy') # 形状:(H, W, 2),I/Q分量
I = sar_data[..., 0]
Q = sar_data[..., 1]
phase = np.arctan2(Q, I) # 计算相位信息
# 2. 生成版权水印并加密
copyright_info = b"Producer:MiniSAR-Lab, Time:20250618, ID:MS2025001"
# 混沌置乱
logistic_map = cm.LogisticMap(mu=3.99, x0=0.31415)
scrambled_watermark = logistic_map.scramble(copyright_info)
# AES-256加密
key = Fernet.generate_key() # 实际应用中由主密钥派生
cipher = Fernet(key)
encrypted_watermark = cipher.encrypt(scrambled_watermark)
# 3. 相位域水印嵌入
block_size = 8
H, W = phase.shape
watermark_idx = 0
watermark_len = len(encrypted_watermark) * 8 # 转换为比特数
for i in range(0, H, block_size):
for j in range(0, W, block_size):
if watermark_idx >= watermark_len:
break
# 提取相位块
phase_block = phase[i:i+block_size, j:j+block_size]
# 计算相位块的冗余位(最小相位扰动)
phase_mean = np.mean(phase_block)
phase_std = np.std(phase_block)
embed_strength = 0.02 if phase_std > 0.1 else 0.05 # 动态调整嵌入强度
# 嵌入1bit水印(0→相位微调+embed_strength,1→-embed_strength)
bit = (encrypted_watermark[watermark_idx//8] >> (7 - watermark_idx%8)) & 0x01
phase_block_embed = phase_block + embed_strength * (1 - 2*bit)
# 更新相位信息
phase[i:i+block_size, j:j+block_size] = phase_block_embed
watermark_idx += 1
# 4. 重构SAR复数影像
I_embed = np.abs(sar_data).astype(np.float32) * np.cos(phase)
Q_embed = np.abs(sar_data).astype(np.float32) * np.sin(phase)
sar_data_embed = np.stack([I_embed, Q_embed], axis=-1)
# 保存带水印数据
np.save('mini_sar_watermarked.npy', sar_data_embed)
1)数据源:微型SAR传感器(X波段,分辨率0.5m×0.5m)采集的复数影像(512×512像素,单精度浮点型);(2)评价指标
2)硬件平台:Intel Corei7-12700H CPU,16GB DDR5内存,NVIDIA RTX3060 GPU;
3)对比算法:传统空域水印(LSB算法)、小波域水印(DWT算法)、本文提出的相位域加密水印算法。
1)不可见性:峰值信噪比(PSNR)、相位误差(PE)、后向散射系数误差(σ⁰误差);
2)鲁棒性:水印提取率(ER),即经过各类处理后成功提取水印的比例;
3)实时性:水印嵌入延迟、提取延迟;
4)安全性:抗水印移除攻击、抗替换攻击的能力。
|
算法 |
PSNR(dB) |
相位误差(rad) |
σ⁰误差(dB) |
|
LSB空域水印 |
32.8 |
0.087 |
0.32 |
|
DWT小波域水印 |
36.5 |
0.042 |
0.18 |
|
本文相位域加密水印 |
39.2 |
0.021 |
0.09 |
1)常规处理:JPEG2000压缩(压缩比20:1)、高斯噪声(σ=0.01)、中值滤波(3×3)、几何校正(旋转10°);结果显示,本文算法在常规处理后的水印提取率均≥95%,其中抗压缩、抗噪声能力尤为突出;面对恶意攻击时,水印提取率仍保持在88%以上,而传统LSB算法在压缩比10:1时提取率已降至65%,DWT算法在几何旋转后提取率仅为72%,验证了本文算法的强鲁棒性。
2)恶意攻击:水印移除攻击(均值滤波5×5)、替换攻击(嵌入伪造水印)。
MiniSAR聚焦于微型合成孔径雷达(SAR)制造研发,为用户提供定制化机载SAR、轻型MiniSAR、无人机载MiniSAR、SAR数据采集服务、SAR飞行服务等。如您有相关业务需求,欢迎联系!