[发明专利]基于ORC采样和QGPCE变换的弥散张量图像鲁棒零水印方法

权利要求书

1.一种基于ORC采样和QGPCE变换的弥散张量图像鲁棒零水印方法，其特征在于：包括，

根据DTI图像数据构建空间张量描述场；

利用光线投射采样策略对所述DTI图像进行正交投射采样，获得基于冠状面、矢状面和横切面的三张正交特征投射图像；

将所述三张正交特征投射图像作为四元数的三个虚部，计算其QGPCE变换系数以构建特征图像；

对混沌置乱的所述特征图像和水印图像进行异或操作获得零水印图像；

泛型极复指数GPCE变换，定度图像f_g(r，θ)的阶数为n，且n为自然数；则得到迭代次数为m，且m为自然数，s为泛型极复指数GPCE变换的参数，如下，

其中，R_ns(r)为径向基函数，*表示复数共轭，为虚数单位；

根据四元数理论，将三幅大小相同特征图像视为一个纯四元数阵列的三个虚部，则利用纯四元数阵列的形式表示三张采样特征图，如下，

f(r，θ)＝f₁(r，θ)i+f₂(r，θ)j+f₃(r，θ)k

其中，f₁、f₂和f₃表示在极坐标下三幅光线投射采样特征图像，通过f(r，θ)可以把三张特征图像视为一个整体来描述；

其中，Re(G)表示复数的实部，Im(G)表示复数的虚部，G_nms(f₁)、G_nms(f₂)和G_nms(f₃)分别为三张投射特征图的泛型极复指数GPCE变换系数/矩，为QGPCE变换系数。

2.根据权利要求1所述的基于ORC采样和QGPCE变换的弥散张量图像鲁棒零水印方法，其特征在于：构建所述空间张量描述场包括，

T＝{D′_xyz|x∈[1，N_x]，y∈[1，N_y]，z∈[1，N_z]}

其中，将所述DTI图像f_D中的每一个弥散张量体素D进行张量分解为D′，得到DTI图像特征值和特征向量形式T，D′_xyz是一个弥散张量体素D_xyz的特征值和特征向量形式，(x，y，z)表示了弥散张量体素D在f_D空间下的局部空间坐标，(N_x，N_y，N_z)表示DTI图像的三维空间最大尺寸。

3.根据权利要求1或2所述的基于ORC采样和QGPCE变换的弥散张量图像鲁棒零水印方法，其特征在于：所述光线投射采样策略包括，

将所述DTI图像表示的生理解剖结构特征进行降维表示；

采样解剖结构方向与投射光线方向相一致的弥散张量体素；

将采样值在二维图像空间进行合成，获得二维投射采样特征图；

其中，R为给定光线矩阵，r₁ r₂ r₃分别表示垂直于冠状面，矢状面和横切面三个方向的正交光线，根据光线矩阵R对DTI图像同时进行冠状面、矢状面和横切面的采样。