CN103812638A - 一种加密域surf图像特征提取方法 - Google Patents
一种加密域surf图像特征提取方法 Download PDFInfo
- Publication number
- CN103812638A CN103812638A CN201410031154.7A CN201410031154A CN103812638A CN 103812638 A CN103812638 A CN 103812638A CN 201410031154 A CN201410031154 A CN 201410031154A CN 103812638 A CN103812638 A CN 103812638A
- Authority
- CN
- China
- Prior art keywords
- partiald
- sigma
- image
- calculate
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
本发明涉及一种加密域的SURF特征提取方法,包括:构建Paillier加密系统和DGK加密系统,生成相应的公匙与私匙;用户端利用Paillier加密系统,以生成的公匙对图像进行加密,然后将加密后的图像发送给服务器端;服务器端对加密后的图像提取SURF特征点;服务器端对提取的SURF特征点进行校正;服务器端提取SURF特征描述子。本发明利用Pallier同态加密方法的同态特性,提出一种加密域的SURF特征提取方法。该方法无需解密即可对加密后的图像直接提取SURF特征,避免了图像信息的泄漏;而且取出的SURF特征点数和位置与明文域算法完全一致,描述子与明文域的误差也仅为0.0002932%。
Description
技术领域
本发明涉及一种加密域的SURF(Speeded Up Robust Features,加速健壮特征)特征提取方法,属于多媒体信息安全处理领域。
背景技术
近年来,随着云计算和物联网等应用对安全需求的不断增加,人们开始研究在加密域内进行信号处理的可能性,即对加密后的信号直接进行处理,可以在保证用户信息安全的同时又不牺牲信号处理的功能。近年来,随着云计算应用的不断推广,加密域信号处理已经成为信号处理领域的一个研究热点。
SURF特征是一种局部特征提取算法。该算法提取的特征具有尺度不变、旋转不变的性能,对光照变化、仿射变化、透视变换具有部分不变性。而且SURF算法在重复度、区分性、鲁棒性三个方面均超越或接近于以往提出的算法。相比于SIFT(Scale Invariant Feature Transform,尺度不变特征变换)算法,其特征提取速度具有明显优势,也因此常被人们看作是SIFT的快速算法,目前已被广泛应用于图像检索、图像配准、人脸识别等各种应用中。
同态加密允许人们对加密后的数据进行特定的代数运算,其得到的处理结果与对明文进行同样的运算再将结果加密一样。Paillier加密算法不仅具有加法同态特性,还可以实现数乘运算,是目前最常用的同态加密算法。相比于RSA加密算法,Paillier加密算法具有更高的运算效率。
申请号为200680019169.5的专利公开了一种用于在媒体信号中嵌入水印的方法和系统。当接收媒体信号的设备不被信任时,该方法和系统可以在媒体信号中嵌入水印,作为媒体信号的身份验证,但是随着媒体信号的近拷贝版本的变化,认证精度无法保证。不如直接检验媒体信号的内容,把媒体信号的内容当做其唯一身份验证信息。
发明内容
本发明基于Pallier同态加密算法,提出了一种加密域的SURF特征提取方法,该方法无需解密,可针对加密后的图像直接提取图像的SURF特征,避免了图像信息的泄漏。
一种加密域的SURF特征提取方法,包括以下步骤:
步骤1:构建Paillier加密系统和DGK加密系统,生成相应的公匙与私匙,用于图像的加密和解密。
步骤2:用户端利用Paillier加密系统,以步骤1生成的公匙对图像进行加密,然后将加 密后的图像发送给服务器端。
步骤3:服务器端对加密后的图像提取SURF特征点。
步骤3.1:计算图像每个点的海森矩阵。
步骤3.2:计算每个海森矩阵对应的行列式值。
(1)将图像的像素值放大到100倍。
(2)计算每个海森矩阵对应的行列式值。
利用方格滤波器与图像的卷积来代替高斯二阶导数与图像的卷积近似计算海森矩阵对应的行列式值。假设x、y和xy方向的滤波器与图像的卷积结果分别为Dxx、Dyy和Dxy,海森矩阵Happrox的行列式值的计算公式为:
det[Happrox]=DxxDyy-(0.9Dxy)2 (1)
步骤3.3:在所有行列式值中寻找局部极值,该局部极值对应的点即为特征点。
步骤4:服务器端对提取的SURF特征点进行校正,即曲线拟合。
设特征点坐标为X=(σ,y,x)T,校正后特征点的坐标X′=(σ′,y′,x′)T,则:
其中,
式中,σ为该点所在的尺度,y为该点的纵坐标,x为该点的横坐标,F表示该点对应的海森矩阵的行列式值。
步骤5:服务器端提取SURF特征描述子。
步骤5.1:以特征点为圆心,计算半径为6σ圆内x和y方向上的Haar小波响应系数HaarX和HaarY,公式如下:
HaarX=(E(x+σ,y+σ)+E(x,y-σ)-E(x,y+σ)-E(x+σ,y-σ))-
(3)
(E(x,y+σ)+E(x-σ,y-σ)-E(x-σ,y+σ)-E(x,y-σ))
HaarY=(E(x-σ,y)+E(x+σ,y+σ)-E(x+σ,y)-E(x-σ,y+σ))-
(4)
(E(x+σ,y)+E(x-σ,y-σ)-E(x-σ,y)-E(x+σ,y-σ))
式中,E(x,y)为积分图像中点(x,y)处的加密值。
步骤5.2:根据Haar小波响应系数计算主方向。
(1)旋转顶点在圆心、夹角为60°的扇形区域,旋转扇形区域的步长为0.15弧度,将整个区域分成42个部分。
(2)判断点(HaarX,HaarY)所属的扇形区域。
假设点(HaarX,HaarY)与原点的连线和横轴的夹角为β,某个扇形区域的两条边与横轴的夹角分别为α,α+60°。如果α<β<α+60°,则点(HaarX,HaarY)属于该扇形区域;否则,继续判断该点与该扇形区域相邻的扇形区域的关系,直到找出该点所属的扇形区域。
在加密域比较角度β与α大小的方法如下:
当与X_α同号时,如果 则β>α,如果 则β<α。
步骤5.3:将Haar响应旋转到主方向,即计算Haar响应值在主方向上的投影;
(1)计算E(sinθ)与E(cosθ),θ为主方向的角度。
(2)计算dx和dy在主方向上的投影d′x和d′y,具体公式如下:
式中,θ为主方向的角度。
所述公式(1)进行加密域乘法的实现方法如下:
若已知加密值c1和c2,其中c1=E(m1),c2=E(m2),ci=E(mi)代表明文mi加密后的值为ci,i=1,2,则加密域的乘法运算就是计算加密值c3,使得D(c3)=m1×m2;E代表加密过程,D代表解密过程。
(1)服务器端生成随机数r1和r2并加密,得到c4=E(r1)、c5=E(r2)和c6=E(-r1r2)。
(2)计算c4×c1modN2,c5×c2modN,并将结果发送给用户端;N为加密系统的参数。
(3)用户端将步骤(2)中的两个结果解密后,进行两个解密值的乘法运算,然后将相乘结果c7加密后发送回服务器端。
(4)服务器端计算 ,其中D(c3)=m1×m2。
所述公式(2)在加密域实现除法运算的方法如下:
假设加密值c1和c2,其中c1=E(m1),c2=E(m2),加密域除法计算就是得到加密值c3,使得 s=107。
(1)服务器端生成随机数r,然后对r加密,得到c5=E(r)。
(2)计算c5×c2=E(m2+r)=c6。
(3)服务器将c6发送给用户端,由于用户端不知道随机数r的大小,所以无法通过解密c6来获取m2的值;用户端解密c6,D(c6)=m3;然后计算并加密得到
(4)用户端将c4的值返回给服务器端。
所述公式(5)在加密域实现开方运算的方法如下:
(1)服务器端生成随机数r,然后对r2加密,得到c3=E(r2)。
(2)计算
(3)服务器将c4发送给用户端,由于用户端不知道随机数r的大小,所以无法通过解密c4来获取m的值;这时,需要用户端对c4解密,得到D(c4)=m2;然后计算并加密得到
(4)用户端将c5的值发送给服务器端。
(5)服务器计算r2的乘法逆元,其值为m3,然后计算 则
与现有技术相比,本发明具有以下优点:
本发明利用Pallier同态加密方法的同态特性,提出了一种加密域的SURF特征提取方法。 该方法无需解密,可针对加密后的图像直接提取图像的SURF特征,避免了图像信息的泄漏;而且取出的SURF特征点数量和位置与明文域算法提取的完全一致,描述子与明文域的误差也仅为0.0002932%。本发明为各种加密域图像处理应用提供了一种有用的工具。
附图说明:
图1为本发明提出的加密域SURF特征提取方法流程图;
图2为本发明提出的角度比较方法示意图;
图3为SURF特征描述子示意图;
图4为加密域SURF和明文域SURF特征点个数及位置的对比图。
具体实施方式
本发明的流程图如图1所示,具体包括以下步骤:
步骤1:构建Paillier加密系统和DGK加密系统,生成相应的公匙与私匙,用于图像的加密和解密。
步骤2:用户端利用Paillier加密系统,以步骤1生成的公匙对图像进行加密,然后将加密后的图像发送给服务器端。
步骤3:服务器端对加密后的图像提取SURF特征点,具体步骤包括:
步骤3.1:对图像的像素值进行100倍的放大。
步骤3.2:计算每个海森矩阵对应的行列式值。
步骤3.3:在所有行列式值中寻找局部极值作为特征点。
步骤4:服务器端对提取的SURF特征点进行校正,利用加密域乘法和加密域除法计算公式(2),获得特征点校正后的值。
步骤5:服务器端提取SURF特征描述子。SURF特征描述子示意图如图3所示。
步骤5.1:以特征点为圆心,计算半径为6σ圆内x和y方向上的Haar小波响应系数。
步骤5.2:根据Haar小波响应系数,计算主方向。
步骤5.3:将Haar响应旋转到主方向,即计算Haar响应值在主方向上的投影。
图4为加密域SURF特征与明文域SURF特征对比图。左边第一列为原始图像,第二列为明文域SURF特征提取结果,第三列为加密域SURF特征提取结果,其中方块代表特征点的位置。图4表明加密域SURF的特征点数量和位置与明文域完全一致。
表1是加密域与明文域提取出的SURF特征描述子的对比。从表1中可以看出,加密域描述子大小约是明文域的1021倍,1021是计算加密域除法和加密域开方运算时乘上的系数。如果将加密域描述子检测结果除以1021,可求得加密域和明文域检测结果的平均误差为 0.0002932%。这个误差来源于加密域近似开方运算和高斯加权值的近似。这一误差对后续的检索、匹配等应用不会造成太多的影响,是可以接受的。
表1加密域与明文域描述子对比
向量 | 明文域 | 加密域 |
第4维向量 | 0.0254506 | 25447978340820304956 |
第14维向量 | -0.144928 | 144913278128088816342 |
第26维向量 | -15.9526 | -15952583853243776628096 |
第32维向量 | 1.03341 | 1033403209007703418500 |
第54维向量 | -0.00993253 | -993051543396542800 |
64维向量总和 | -66.1773 | -66177105905710516176452 |
平均差值百分比 | 0.0002932% |
Claims (4)
1.一种加密域的SURF特征提取方法,其特征在于,无需解密即可针对加密后的图像直接提取图像的SURF特征,所述方法包括以下步骤:
步骤1:构建Paillier加密系统和DGK加密系统,生成相应的公匙与私匙,用于图像的加密和解密;
步骤2:用户端利用Paillier加密系统,以步骤1生成的公匙对图像进行加密,然后将加密后的图像发送给服务器端;
步骤3:服务器端对加密后的图像提取SURF特征点;
步骤3.1:计算图像每个点的海森矩阵;
步骤3.2:计算每个海森矩阵对应的行列式值;
(1)将图像的像素值放大到100倍;
(2)计算每个海森矩阵对应的行列式值;
利用方格滤波器与图像的卷积来代替高斯二阶导数与图像的卷积近似计算海森矩阵对应的行列式值;假设x、y和xy方向的滤波器与图像的卷积结果分别为Dxx、Dyy和Dxy,海森矩阵Happrox的行列式值的计算公式为:
det[Happrox]=DxxDyy-(0.9Dxy)2 (1)
步骤3.3:在所有行列式值中寻找局部极值,该局部极值对应的点即为特征点;
步骤4:服务器端对提取的SURF特征点进行校正,即曲线拟合;
设特征点坐标为X=(σ,y,x)T,校正后特征点的坐标X′=(σ′,y′,x′)T,则:
其中,
式中,σ为该点所在的尺度,y为该点的纵坐标,x为该点的横坐标,F表示该点对应的海森矩阵的行列式值;
步骤5:服务器端提取SURF特征描述子;
步骤5.1:以特征点为圆心,计算半径为6σ圆内x和y方向上的Haar小波响应系数HaarX和HaarY,公式如下:
HaarX=(E(x+σ,y+σ)+E(x,y-σ)-E(x,y+σ)-E(x+σ,y-σ))-
(3)
(E(x,y+σ)+E(x-σ,y-σ)-E(x-σ,y+σ)-E(x,y-σ))
HaarY=(E(x-σ,y)+E(x+σ,y+σ)-E(x+σ,y)-E(x-σ,y+σ))-
(4)
(E(x+σ,y)+E(x-σ,y-σ)-E(x-σ,y)-E(x+σ,y-σ))
式中,E(x,y)为积分图像中点(x,y)处的加密值;
步骤5.2:根据Haar小波响应系数计算主方向;
(1)旋转顶点在圆心、夹角为60°的扇形区域,旋转扇形区域的步长为0.15弧度,将整个区域分成42个部分;
(2)判断点(HaarX,HaarY)所属的扇形区域;
假设点(HaarX,HaarY)与原点的连线和横轴的夹角为β,某个扇形区域的两条边与横轴的夹角分别为α,α+60°;如果α<β<α+60°,则点(HaarX,HaarY)属于该扇形区域;否则,继续判断该点与该扇形区域相邻的扇形区域的关系,直到找出该点所属的扇形区域;
在加密域比较角度β与α大小的方法如下:
当与X_α同号时,如果 则β>α,如果 则β<α;
步骤5.3:将Haar响应旋转到主方向,即计算Haar响应值在主方向上的投影;
(1)计算E(sinθ)与E(cosθ),θ为主方向的角度;
(2)计算dx和dy在主方向上的投影d′x和d′y,具体公式如下:
式中,θ为主方向的角度;
2.根据权利要求1所述的一种加密域的SURF特征提取方法,其特征在于,所述公式(1)进行加密域乘法的实现方法如下:
若已知加密值c1和c2,其中c1=E(m1),c2=E(m2),ci=E(mi)代表明文mi加密后的值为ci,i=1,2,则加密域的乘法运算就是计算加密值c3,使得D(c3)=m1×m2;E代表加密过程,D代表解密过程;
(1)服务器端生成随机数r1和r2并加密,得到c4=E(r1)、c5=E(r2)和c6=E(-r1r2);
(2)计算c4×c1modN2,c5×c2modN,并将结果发送给用户端;N为加密系统的参数;
(3)用户端将步骤(2)中的两个结果解密后,进行两个解密值的乘法运算,然后将相乘结果c7加密后发送回服务器端;
(4)服务器端计算 ,其中D(c3)=m1×m2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410031154.7A CN103812638B (zh) | 2014-01-22 | 2014-01-22 | 一种加密域surf图像特征提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410031154.7A CN103812638B (zh) | 2014-01-22 | 2014-01-22 | 一种加密域surf图像特征提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103812638A true CN103812638A (zh) | 2014-05-21 |
CN103812638B CN103812638B (zh) | 2017-02-22 |
Family
ID=50708888
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410031154.7A Active CN103812638B (zh) | 2014-01-22 | 2014-01-22 | 一种加密域surf图像特征提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103812638B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104573058A (zh) * | 2015-01-22 | 2015-04-29 | 山东大学 | 一种基于HBase的海量图片存储方法及其实现系统 |
CN104883580A (zh) * | 2015-06-03 | 2015-09-02 | 合肥工业大学 | 一种基于同态加密的视频安全卷积计算系统及方法 |
CN106250898A (zh) * | 2016-07-28 | 2016-12-21 | 哈尔滨工业大学 | 一种基于尺度预测的图像局部区域特征提取方法 |
CN106529464A (zh) * | 2016-10-28 | 2017-03-22 | 信利光电股份有限公司 | 一种虹膜特征采集模组及电子设备 |
CN106952212A (zh) * | 2017-03-14 | 2017-07-14 | 电子科技大学 | 一种基于向量同态加密的hog 图像特征提取算法 |
CN108418995A (zh) * | 2018-03-20 | 2018-08-17 | 西安电子科技大学 | 基于同态加密算法的个人图像安全检索方法 |
WO2021068444A1 (zh) * | 2019-10-11 | 2021-04-15 | 云图技术有限公司 | 数据处理方法、装置、计算机设备和存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080212780A1 (en) * | 2005-06-03 | 2008-09-04 | Koninklijke Philips Electronics, N.V. | Homomorphic Encryption For Secure Watermarking |
CN101674299B (zh) * | 2009-10-16 | 2013-04-17 | 西安电子科技大学 | 密钥生成方法 |
CN103336970B (zh) * | 2013-06-19 | 2016-12-28 | 电子科技大学 | 基于多子树的分布式图像训练及检索方法 |
CN103442253B (zh) * | 2013-07-17 | 2016-08-10 | 宁波工程学院 | 一种加密域h.264/avc视频数据隐藏方法 |
-
2014
- 2014-01-22 CN CN201410031154.7A patent/CN103812638B/zh active Active
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104573058A (zh) * | 2015-01-22 | 2015-04-29 | 山东大学 | 一种基于HBase的海量图片存储方法及其实现系统 |
CN104573058B (zh) * | 2015-01-22 | 2018-11-13 | 山东大学 | 一种基于HBase的海量图片存储方法及其实现系统 |
CN104883580A (zh) * | 2015-06-03 | 2015-09-02 | 合肥工业大学 | 一种基于同态加密的视频安全卷积计算系统及方法 |
CN104883580B (zh) * | 2015-06-03 | 2020-12-11 | 合肥工业大学 | 一种基于同态加密的视频安全卷积计算系统及方法 |
CN106250898A (zh) * | 2016-07-28 | 2016-12-21 | 哈尔滨工业大学 | 一种基于尺度预测的图像局部区域特征提取方法 |
CN106250898B (zh) * | 2016-07-28 | 2019-04-26 | 哈尔滨工业大学 | 一种基于尺度预测的图像局部区域特征提取方法 |
CN106529464A (zh) * | 2016-10-28 | 2017-03-22 | 信利光电股份有限公司 | 一种虹膜特征采集模组及电子设备 |
CN106952212A (zh) * | 2017-03-14 | 2017-07-14 | 电子科技大学 | 一种基于向量同态加密的hog 图像特征提取算法 |
CN106952212B (zh) * | 2017-03-14 | 2019-10-15 | 电子科技大学 | 一种基于向量同态加密的hog图像特征提取算法 |
CN108418995A (zh) * | 2018-03-20 | 2018-08-17 | 西安电子科技大学 | 基于同态加密算法的个人图像安全检索方法 |
CN108418995B (zh) * | 2018-03-20 | 2019-10-11 | 西安电子科技大学 | 基于同态加密算法的个人图像安全检索方法 |
WO2021068444A1 (zh) * | 2019-10-11 | 2021-04-15 | 云图技术有限公司 | 数据处理方法、装置、计算机设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103812638B (zh) | 2017-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103812638A (zh) | 一种加密域surf图像特征提取方法 | |
CN103744976B (zh) | 一种基于同态加密的图像安全检索方法 | |
US9749128B2 (en) | Compact fuzzy private matching using a fully-homomorphic encryption scheme | |
Gao et al. | Bessel-Fourier moment-based robust image zero-watermarking | |
US9503266B2 (en) | Relational encryption | |
Barman et al. | Fingerprint-based crypto-biometric system for network security | |
Irawan et al. | Hiding and securing message on edge areas of image using LSB steganography and OTP encryption | |
CN100507936C (zh) | 一种基于Zernike矩的鲁棒哈希图像认证方法 | |
CN105046234B (zh) | 基于稀疏表示的云环境中人脸图像隐秘识别方法 | |
US9577831B2 (en) | Relational encryption | |
US20170142109A1 (en) | Relational encryption | |
CN103023633B (zh) | 基于混沌随机相位和相干叠加原理的数字图像隐藏方法 | |
CN109803142B (zh) | 一种基于邻域预测的加密域图像可逆数据隐藏方法 | |
JP2009129292A (ja) | 生体認証方法、装置およびシステム | |
CN104166955A (zh) | 基于保角变换图像哈希生成及图像篡改检测定位方法 | |
CN112948795B (zh) | 保护隐私的身份认证方法及装置 | |
Da et al. | A novel hybrid information security scheme for 2D vector map | |
CN105553980A (zh) | 基于云计算的安全指纹识别系统和方法 | |
Singh et al. | A robust image hashing based on discrete wavelet transform | |
Im et al. | Privacy-preserving palm print authentication using homomorphic encryption | |
Singh et al. | Images as graphical password: verification and analysis using non-regular low-density parity check coding | |
Hazarika et al. | A wavelet based partial image encryption using chaotic logistic map | |
Jin et al. | Privacy preserving face identification in the cloud through sparse representation | |
CN104717072A (zh) | 基于感知哈希与椭圆曲线的遥感影像认证方法 | |
CN103810407B (zh) | 一种gis矢量数据线面图层版权认证方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220817 Address after: 100044 1-14-2107, floor 2, No. 136, Xiwai street, Xicheng District, Beijing Patentee after: Shuanxin (Beijing) Technology Co.,Ltd. Address before: 100124 No. 100 Chaoyang District Ping Tian Park, Beijing Patentee before: Beijing University of Technology |
|
TR01 | Transfer of patent right |