基于同态加密保护的宝石图像特征识别方法与流程

本发明涉及的是一种宝石鉴定领域的技术，具体是一种基于同态加密保护的宝石图像特征识别方法。

背景技术：

宝石合成技术越来越先进，合成宝石甚至使用一些传统的简易检测手段都无法检出，这对整个宝石行业已经造成了严重的问题。

技术实现要素：

本发明针对现有技术存在的上述不足，提出一种基于同态加密保护的宝石图像特征识别方法，使用宝石的放大镜下的清晰照片，针对每颗宝石的独特物理特征，诸如包裹体、内含物等，形成每颗独一无二的身份编码，类似人类的指纹，因此可以通过这种方式识别宝石的身份。但是，如果造假者掌握了这些物理特征数据，也有可能加以伪造。因此，需要对宝石物理特征数据进行保护，使之既可以运算来实现宝石身份识别，同时不会泄露这些数据。

本发明是通过以下技术方案实现的：

本发明通过对宝石检测图像进行特征提取并与宝石内含物特征一并编码后加密生成识别库，再通过对待测宝石的特征编码进行比较实现识别。

所述的特征提取，基于宝石检测图像采用具有旋转不变性和灰度不变性的rilbp算子提取得到后，将rilbp指标替代像素灰度，进行图像缩减和旋转后，再使用haar-like指标提取图像特征。

所述的宝石内含物特征是指：通过内含物边缘曲线基于不规则图形的曲线边界的四顶点表达式、八顶点表达式、内含物的重心的坐标以及边界点最小夹表征。

所述的加密，包括对提取到的特征编码后使用完全同态加密技术将编码中的参数取值进行加密作为身份特征编码。

所述的比较，利用公钥将待识别宝石进行特征编码，并将编码结果与预存的身份特征编码进行比对以实现识别。

本发明涉及一种实现上述方法的系统，包括：特征提取单元、内含物特征单元、加密单元以及比较单元，其中：特征提取单元和内含物特征单元分别获取真实宝石检测图像和待测宝石检测图像并将分别提取到的对应haar-like特征向量和宝石内含物特征输出至加密单元，比较单元通过公开方式获得公钥并输出至加密单元，加密单元采用多项式完全同态加密算法分别采用来自比较单元的公钥和内置的私钥进行加密并得到真实宝石和待测宝石的特征编码并实现比较。

技术效果

与现有技术相比，本发明采用两种截然不同的特征定义和提取方法，使用计算机图形学和宝石学两种独立的知识体系来定义特征。这样可以在宝石识别时，根据具体场景要求灵活选择恰当的特征提取方法，获得稳定可靠的识别结果。同时，采用这种独立双特征的识别模式，在识别图像质量较低，或者存在变形的条件下，也同样可以有效的进行识别；进一步配合完全同态加密方法，特征编码可以在不泄露明文的情况下直接在加密状态下对待识别宝石特征进行比对运算。

附图说明

图1为宝石特征结构示意图；

图2为实施例宝石示意图；

图3为图2特征提取示意图；

图4a～图4e为几何特征提取过程示意图；

图5为宝石特征提取过程示意图；

图6为实施例流程示意图。

具体实施方式

本实施例包括：

步骤一、图片灰度纹理特征提取：根据宝石检测图像，采用具有旋转不变性和灰度不变性的rilbp算子来提取检测图像特征，将rilbp指标替代像素灰度，进行图像缩减和旋转后，再使用haar-like指标提取图像特征。

所述的rilbp(rotationinvariantlbp)算子是基于lbp(localbinarypattern，局部二值模式)的改进方案，实现了lbp特征的旋转不变性，实现在输入图像任意角度旋转下始终保持特征取值的稳定性，即对转换前像素点对应的原始lbp二进制编码循环移位，一个像素点对应8种二进制编码，然后将这8种二进制编码转换为十进制数值，找到其中最小数值作为像素点的特征值，具体包括：

i)针对宝石图像进行预处理，包括：拍摄设备及环境标准、特征平面选取、照片清晰度筛选、照片颜色矫正、滤波处理、边缘锐化处理、二值化处理等，确定轮廓，缩放识别区域到标准像素尺寸得到如图2和图3所示的特征照片。

ii)遍历识别区域图像，计算每个3*3像素区域的rilbp指标，并以中心点替代像素区域，缩减识别区域图像，总像素缩减为原图的1/9；

iii)对缩减后的识别区域图像进行旋转变换，让轮廓上拥有最大rilbp指标之和的3*3像素区域位于图像底部。

所述的haar-like指标包括：针对旋转后的图像提取haar-like系列指标，包含16个特征结构，如图1所示，将图1中每个haar-like特征的白色区域内的像素值之和减去黑色部分的像素和，得到haar-like特征值，使用特征向量[h1，h2，h3，h4，h5，h6，h7，h8，h9，h10，h11，h12，h13，h14，h15，h16]作为待标识区域的特征，记作h，即haar-like特征向量。

步骤二、使用不规则图形的曲线边界表示法表征宝石内含物特征，即通过内含物边缘曲线基于不规则图形的曲线边界的四顶点表达式、八顶点表达式、内含物的重心的坐标以及边界点最小夹表征，具体步骤包括：

2.1)选取适当平面坐标系统以确定宝石边缘曲线，计算得到重心坐标o＝(xo,yo)；

2.2)测得边缘曲线上离重心最远的边界点集合a，在a中找到一点a0，使得以oa0作为0轴的极坐标系统下，a中的所有点极坐标夹角总和最小

2.3)更新极坐标系统至极轴为oa0，如图4所示；使用不规则图形的曲线边界表示法，在极坐标系统oa0下写出宝石边缘曲线的四顶点表达式b0和八顶点表达式bn。

如图4所示，所述的不规则图形的曲线边界表示法是指：

一、通过宝石边缘曲线上距离最远的两个点的连线d1d2将其分为两部分，分别在两段宝石边缘曲线上获得距离直线d1d2最远的点d3、d4，并将宝石边缘曲线分成四段；

二、在曲线d1d3上找到一与外部边界的矩形面积最大的点，例如d5，采用类似方式得到另外三段曲线上的三个点，组成{d5，d6，d7，d8}；

三、以d4和d1为顶点做椭圆形所得该1/4部分椭圆面积与该部分宝石边缘曲线所围面积之比为r41；采用类似方式得到另外三段曲线上的面积比，组成四个曲线段的内切椭圆面积覆盖率{r13，r32，r24，r41}；

四、以{d1，d2，d3，d4，d5，d6，d7，d8，r13，r32，r24，r41}作为该不规则曲线边界的八顶点表达式，以{d1，d2，d3，d4，r13，r32，r24，r41}作为该不规则曲线边界的四顶点表达式。

步骤三、以基于系数映射的多项式完全同态加密算法对{h，p，b0，[{bn，on，θn}]}进行加密，得到受保护的特征编码ecode0，其中：h为haar-like特征向量；p为宝石主要物理特征参数测量；b0为内含物边缘曲线基于不规则图形的曲线边界的四顶点表达式；bn为内含物边缘曲线基于不规则图形的曲线边界的八顶点表达式；on为内含物的重心的坐标；θn为边界点最小夹角。

所述的主要物理特征参数测量包括：重量、外形尺寸、莫氏硬度、折射率；

所述的加密，具体为：根据n维明文向量p＝(v1,v2,...,vn)，加密密钥为k＝{f,z}，设则有a1＝{a1i,i∈i},an＝{ani,i∈i},x＝{xi,i∈i},i＝{1,2}，得到密文为c＝{a1,a2,...,an,x}，该密文表示法比针对特征向量的每一个分量进行加密得到的密文向量节省几乎一半的存储空间和计算时间，这种表达形式下，由于对所有分量，其密文的x均相同，因此向量密文的加减法可以非常快速的进行处理，同时可以节省大量的密文存储空间。

所述的加密，优选先将{h，p，b0，[{bn，on，θn}]}展开成一维向量，然后进行加密，得到加密特征编码ecode0。

所述的加密特征编码ecode0可以直接保持在云服务中，或者通过qrcode等形式打印在宝石随附的证书上。

步骤四、当需要对宝石进行识别时，重复步骤一和步骤二的特征提取步骤，即通过放大镜照片获取待识别宝石的haar-like特征和内含物轮廓几何特征，然后使用识别服务所提供的比较公钥，对特征值进行加密，得到特征编码ecode1，对ecode0和ecode1求差，并将所得到的差deltacode中每个分量di与对应的阈值σi进行比较，使用密文比较运算得到di-σi＜0的结果，当对任意i∈i，存在di-σi＜0的结果为false，则可以判断待识别宝石不是登记的宝石；否则可以判断待识别宝石就是登记的宝石。

上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整，本发明的保护范围以权利要求书为准且不由上述具体实施所限，在其范围内的各个实现方案均受本发明之约束。