• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 14
    • 下载费用:30 金币  

    重庆时时彩二星杀号技巧: 一种数字图像隐写分析特征提取方法及提取系统.pdf

    关 键 词:
    一种 数字图像 分析 特征 提取 方法 系统
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201611264561.8

    申请日:

    2016.12.30

    公开号:

    CN106683031A

    公开日:

    2017.05.17

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06T 1/00申请日:20161230|||公开
    IPC分类号: G06T1/00 主分类号: G06T1/00
    申请人: 深圳大学
    发明人: 李斌; 李忠蓬; 黄继武
    地址: 518060 广东省深圳市南山区南海大道3688号
    优先权:
    专利代理机构: 深圳市君胜知识产权代理事务所(普通合伙) 44268 代理人: 王永文;刘文求
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201611264561.8

    授权公告号:

    |||

    法律状态公告日:

    2017.06.09|||2017.05.17

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明公开了一种数字图像隐写分析特征提取方法及提取系统,方法包括:将待提取隐写分析特征图像通过10个滤波器进行滤波,得到残差图像;其中,10个滤波器包括4个无方向性滤波器和6个有方向性滤波器;根据LBPriu2模型和/或DLBP模型对残差图像进行编码后得到编码图,并根据编码图得到4个方向的二阶共生矩阵,对其中6个有方向性滤波器,相加合并不同方向计算得到的共生矩阵;对计算得到的共生矩阵,合并不同方向共生矩阵并降维后进行LOG函数映射得到共生矩阵特征。本发明与现有的主流通用隐写分析方法SRM具有相近的性能,分类效果较SRM具有更高的准确率。

    权利要求书

    1.一种数字图像隐写分析特征提取方法,其特征在于,所述方法包括以下步骤:
    A、将待提取隐写分析特征图像通过10个滤波器进行滤波,得到残差图像;其中,10个滤
    波器包括4个无方向性滤波器和6个有方向性滤波器;
    B、根据LBPriu2模型和/或DLBP模型对残差图像进行编码后得到编码图,并根据编码图
    得到4个方向的二阶共生矩阵,对其中6个有方向性滤波器,相加合并不同方向计算得到的
    共生矩阵;
    C、对计算得到的共生矩阵,合并不同方向共生矩阵并降维后进行LOG函数映射得到共
    生矩阵特征。
    2.根据权利要求1所述数字图像隐写分析特征提取方法,其特征在于,所述步骤A中10
    个滤波器所包括的4个基本滤波器分别记为D12、D13、D14和D15,其中:
    [1,0]-[0,1]=[1,-1] (D12)
    [1,-1,0]-[0,1,-1]=[1,-2,1] (D13)
    [1,-2,1,0]-[0,1,-2,1]=[1,-3,3,-1] (D14)
    [1,-3,3,-1,0]-[0,1,-3,3,-1]=[1,-4,6,-4,1] (D15);
    10个滤波器所包括的6个由基本滤波器推导出的滤波器分别记为D22、D33、D44、D55、
    D23和D35,其中:
    <mrow> <mi>D</mi> <mn>22</mn> <mo>=</mo> <mi>D</mi> <msup> <mn>12</mn> <mo>&prime;</mo> </msup> <mo>&times;</mo> <mi>D</mi> <mn>12</mn> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> <mi>D</mi> <mn>33</mn> <mo>=</mo> <mi>D</mi> <msup> <mn>13</mn> <mo>&prime;</mo> </msup> <mo>&times;</mo> <mi>D</mi> <mn>13</mn> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </mtd> <mtd> <mn>4</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
    <mrow> <mi>D</mi> <mn>44</mn> <mo>=</mo> <mi>D</mi> <msup> <mn>14</mn> <mo>&prime;</mo> </msup> <mo>&times;</mo> <mi>D</mi> <mn>14</mn> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> <mtd> <mn>3</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>3</mn> </mrow> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mn>3</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>9</mn> </mrow> </mtd> <mtd> <mn>9</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>3</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>3</mn> </mrow> </mtd> <mtd> <mn>9</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>9</mn> </mrow> </mtd> <mtd> <mn>3</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>3</mn> </mrow> </mtd> <mtd> <mn>3</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> <mi>D</mi> <mn>55</mn> <mi>=</mi> <mi>D</mi> <msup> <mn>15</mn> <mo>&prime;</mo> </msup> <mo>&times;</mo> <mi>D</mi> <mn>15</mn> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>4</mn> </mrow> </mtd> <mtd> <mn>6</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>4</mn> </mrow> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>4</mn> </mrow> </mtd> <mtd> <mn>16</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>24</mn> </mrow> </mtd> <mtd> <mn>16</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>4</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>6</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>24</mn> </mrow> </mtd> <mtd> <mn>36</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>24</mn> </mrow> </mtd> <mtd> <mn>6</mn> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mn>4</mn> </mrow> </mtd> <mtd> <mn>16</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>24</mn> </mrow> </mtd> <mtd> <mn>26</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>5</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>4</mn> </mrow> </mtd> <mtd> <mn>6</mn> </mtd> <mtd> <mrow> <mo>-</mo> <mn>4</mn> </mrow> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>


    其中,D22、D33、D44、D55是4个无方向性滤波器,D12、D13、D14、D15、D23、D35是6个有方
    向性滤波器,对有方向性滤波器,分别计算各个方向上的余量以及计算其中各个方向余量
    的最大值余量和最小值余量。
    3.根据权利要求2所述数字图像隐写分析特征提取方法,其特征在于,所述步骤A中通
    过10个滤波器进行滤波得到残差图像的计算式为:
    <mrow> <mi>R</mi> <mo>=</mo> <mi>I</mi> <mo>&CircleTimes;</mo> <mi>f</mi> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mo>-</mo> <mi>M</mi> </mrow> <mi>M</mi> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mo>-</mo> <mi>N</mi> </mrow> <mi>N</mi> </munderover> <mi>f</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mi>I</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mi>m</mi> <mo>,</mo> <mi>j</mi> <mo>-</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow>

    其中,I为待提取隐写分析特征图像,f是滤波器,是卷积符号,Row(f)和Colum(f)分
    别为求f的行和列数,M、N分别是编码图的长和宽。
    4.根据权利要求1所述数字图像隐写分析特征提取方法,其特征在于,所述步骤B中
    LBPriu2模型的计算式为:

    其中,gp是相邻像素值,gc是中间像素值,P、R是LBP编码中定义的相邻像素的个数和位
    置,取P=8和R=1。
    5.根据权利要求1所述数字图像隐写分析特征提取方法,其特征在于,所述步骤B中
    DLBP模型的计算式为:

    其中,其中gp是相邻像素值,gc是中间像素值,P、R是LBP编码中定义的相邻像素的个数
    和选取半径,取P=8和R=1,α根据经验设置成0.3和0.5,fi,j是滤波器系数。
    6.根据权利要求1所述数字图像隐写分析特征提取方法,其特征在于,所述步骤B中4个
    方向的二阶共生矩阵的生成式为:
    C0(i,j)=#{(x1,y1),(x1,y1+1)∈M×N|RLBP(x1,y1)=i,RLBP(x1,y1+1)=j}
    C45(i,j)=#{(x1,y1),(x1-1,y1+1)∈M×N|RLBP(x1,y1)=i,RLBP(x1-1,y1+1)=j}
    C90(i,j)=#{(x1,y1),(x1-1,y1)∈M×N|RLBP(x1,y1)=i,RLBP(x1-1,y1)=j}
    C135(i,j)=#{(x1,y1),(x1-1,y1+1)∈M×N|RLBP(x1,y1)=i,RLBP(x1-1,y1+1)=j};
    其中,C0为0°方向上二阶共生矩阵、C45为45°方向上二阶共生矩阵、C90为90°方向上二阶
    共生矩阵、C135为135°方向上二阶共生矩阵,M、N分别是编码图的长和宽,RLBP是LBP或者DLBP
    编码图。
    7.根据权利要求1所述数字图像隐写分析特征提取方法,其特征在于,所述步骤C中具
    体包括:
    C1、对6个有方向性滤波器进行滤波得到不同方向的残差图像,计算共生矩阵并相加合
    并同一个滤波器不同方向的共生矩阵,对6个有方向性滤波器最大值残差图像和最小值残
    差图像、以及4个无方向性滤波器残差图像得到的共生矩阵不合并;
    C2、将共生矩阵沿左上到右下对角线对折相加,并拉伸成为0°方向上的一维向量Cf0、
    45°方向上的一维向量Cf45、90°方向上的一维向量Cf90、及135°方向上的一维向量Cf135;
    C3、将0°方向上的一维向量Cf0与90°方向上的一维向量Cf90相加得到Cf0+90,将45°方向
    上的一维向量Cf45与135°方向上的一维向量Cf135相加得到Cf45+135;
    C4、分别计算f1=log10(Cf0+90+1)和f2=log10(Cf45+135+1),将f1和f2串联后得到共生矩
    阵特征。
    8.一种数字图像隐写分析特征提取系统,其特征在于,包括:
    滤波???,用于将待提取隐写分析特征图像通过10个滤波器进行滤波,得到残差图像;
    其中,10个滤波器包括4个无方向性滤波器和6个有方向性滤波器;
    编码???,用于根据LBPriu2模型和/或DLBP模型对残差图像进行编码后得到编码图,
    并根据编码图得到4个方向的二阶共生矩阵,对其中6个有方向性滤波器,相加合并不同方
    向计算得到的共生矩阵;
    共生矩阵特征提取???,用于对计算得到的共生矩阵,合并不同方向共生矩阵并降维
    后进行LOG函数映射得到共生矩阵特征。
    9.根据权利要求8所述数字图像隐写分析特征提取方法,其特征在于,所述编码??橹?br />LBPriu2模型的计算式为:

    其中,gp是相邻像素值,gc是中间像素值,P、R是LBP编码中定义的相邻像素的个数和位
    置,取P=8和R=1。
    10.根据权利要求8所述数字图像隐写分析特征提取方法,其特征在于,所述编码???br />中DLBP模型的计算式为:

    其中,其中gp是相邻像素值,gc是中间像素值,P、R是LBP编码中定义的相邻像素的个数
    和选取半径,取P=8和R=1,α根据经验设置成0.3和0.5,fi,j是滤波器系数。

    说明书

    一种数字图像隐写分析特征提取方法及提取系统

    技术领域

    本发明涉及数字图像隐写分析技术领域,尤其涉及的是一种数字图像隐写分析特
    征提取方法及提取系统。

    背景技术

    数字图像隐写分析有主动隐写分析和被动隐写分析两种类型,其中被动隐写分析
    主要判断载体是否隐藏秘密信息,而主动隐写分析在判别出载体有秘密信息的同时还要通
    过估计隐写算法的参数进一步提取出秘密信息。根据使用场景,隐写分析又分为专用隐写
    分析和通用隐写分析。通用隐写方法需要对多种隐写算法进行检测,其操作过程主要分为
    特征提取和分类器分类两个部分,其中特征设计部分是研究的主要内容也是本专利设计的
    算法作用环节,而分类操作现在主流的方法是使用集成分类器。

    图像空域通用隐写分析特征中,最有代表性的是空域富模型(SRM),它首先通过许
    多线性和非线性滤波器得到图像不同的余量参差,之后对余量进行量化截断,再通过将不
    同量化截断滤波余量图生成四阶共生矩阵,最后用不同的合并方式降维生成特征。富模型
    特征在隐写分析上对内容自适应隐写取得了不错的效果,SRM在滤波器的选择上较充分考
    虑但是其滤波器规律不足这方面值得改进。

    随着内容自适应隐写方法HUGO的出现,以最小失真框架为基础的隐写算法陆续出
    现,这些算法通过定义一种失真函数来反映嵌入带来的代价,通过最小化整体代价值来将
    信息嵌入在纹理和边缘区域,以此更难被检测。针对隐写算法将秘密信息嵌入载体图像过
    程中,会造成图像纹理的改变,选择对纹理区域更敏感的检测可以对分析隐写有很大帮助。
    局部二值模型(LBP)是一种将结构和统计相结合的算子,它的直方图特征在图像纹理分类
    中得到很好的应用,利用局部二值模型直方图特征在空域隐写分析上已经取得一定效果,
    但是其性能还是劣于SRM,直方图特征对于图像像素间相关性上表达不足,难以捕捉隐写带
    来图像统计特性的改变,另外LBP的定义和隐写分析问题不是直接匹配也是其不足之处。

    因此,现有技术还有待于改进和发展。

    发明内容

    鉴于上述现有技术的不足,本发明的目的在于提供一种数字图像隐写分析特征提
    取方法及提取系统,旨在解决现有技术中利用局部二值模型直方图特征在空域隐写分析上
    直方图特征对于图像像素间相关性上表达不足,难以捕捉隐写带来图像统计特性的改变,
    LBP的定义和隐写分析问题不是直接匹配的缺陷等问题。

    本发明的技术方案如下:

    一种数字图像隐写分析特征提取方法,其中,所述方法包括以下步骤:

    A、将待提取隐写分析特征图像通过10个滤波器进行滤波,得到残差图像;其中,10
    个滤波器包括4个无方向性滤波器和6个有方向性滤波器;

    B、根据LBPriu2模型和/或DLBP模型对残差图像进行编码后得到编码图,并根据编
    码图得到4个方向的二阶共生矩阵,对其中6个有方向性滤波器,相加合并不同方向计算得
    到的共生矩阵;

    C、对计算得到的共生矩阵,合并不同方向共生矩阵并降维后进行LOG函数映射得
    到共生矩阵特征。

    所述数字图像隐写分析特征提取方法,其中,所述步骤A中10个滤波器所包括的4
    个基本滤波器分别记为D12、D13、D14和D15,其中:

    [1,0]-[0,1]=[1,-1] (D12)

    [1,-1,0]-[0,1,-1]=[1,-2,1] (D13)

    [1,-2,1,0]-[0,1,-2,1]=[1,-3,3,-1] (D14)

    [1,-3,3,-1,0]-[0,1,-3,3,-1]=[1,-4,6,-4,1] (D15);

    10个滤波器所包括的6个由基本滤波器推导出的滤波器分别记为D22、D33、D44、
    D55、D23和D35,其中:





    其中,D22、D33、D44、D55是4个无方向性滤波器,D12、D13、D14、D15、D23、D35是6个
    有方向性滤波器,对有方向性滤波器,分别计算各个方向上的余量以及计算其中各个方向
    余量的最大值余量和最小值余量。

    所述数字图像隐写分析特征提取方法,其中,所述步骤A中通过10个滤波器进行滤
    波得到残差图像的计算式为:



    其中,I为待提取隐写分析特征图像,f是滤波器,是卷积符号,Row(f)和Colum
    (f)分别为求f的行和列数,M、N分别是编码图的长和宽。

    所述数字图像隐写分析特征提取方法,其中,所述步骤B中LBPriu2模型的计算式
    为:


    其中


    其中

    其中,gp是相邻像素值,gc是中间像素值,P、R是LBP编码中定义的相邻像素的个数
    和位置,取P=8和R=1。

    所述数字图像隐写分析特征提取方法,其中,所述步骤B中DLBP模型的计算式为:


    其中


    其中,其中gp是相邻像素值,gc是中间像素值,P、R是LBP编码中定义的相邻像素的
    个数和选取半径,取P=8和R=1,α根据经验设置成0.3和0.5,fi,j是滤波器系数。

    所述数字图像隐写分析特征提取方法,其中,所述步骤B中4个方向的二阶共生矩
    阵的生成式为:

    C0(i,j)=#{(x1,y1),(x1,y1+1)∈M×N|RLBP(x1,y1)=i,RLBP(x1,y1+1)=j}

    C45(i,j)=#{(x1,y1),(x1-1,y1+1)∈M×N|RLBP(x1,y1)=i,RLBP(x1-1,y1+1)=j}

    C90(i,j)=#{(x1,y1),(x1-1,y1)∈M×N|RLBP(x1,y1)=i,RLBP(x1-1,y1)=j}

    C135(i,j)=#{(x1,y1),(x1-1,y1+1)∈M×N|RLBP(x1,y1)=i,RLBP(x1-1,y1+1)=
    j};

    其中,C0为0°方向上二阶共生矩阵、C45为45°方向上二阶共生矩阵、C90为90°方向上
    二阶共生矩阵、C135为135°方向上二阶共生矩阵,M、N是编码图的长和宽,RLBP是LBP或者DLBP
    编码图。

    所述数字图像隐写分析特征提取方法,其中,所述步骤C中具体包括:

    C1、对6个有方向性滤波器进行滤波得到不同方向的残差图像,计算共生矩阵并相
    加合并同一个滤波器不同方向的共生矩阵,对6个有方向性滤波器最大值残差图像和最小
    值残差图像、以及4个无方向性滤波器残差图像得到的共生矩阵不合并;

    C2、将共生矩阵沿左上到右下对角线对折相加,并拉伸成为0°方向上的一维向量
    Cf0、45°方向上的一维向量Cf45、90°方向上的一维向量Cf90、及135°方向上的一维向量Cf135;

    C3、将0°方向上的一维向量Cf0与90°方向上的一维向量Cf90相加得到Cf0+90,将45°
    方向上的一维向量Cf45与135°方向上的一维向量Cf135相加得到Cf45+135;

    C4、分别计算f1=log10(Cf0+90+1)和f2=log10(Cf45+135+1),将f1和f2串联后得到共
    生矩阵特征。

    一种数字图像隐写分析特征提取系统,其中,包括:

    滤波???,用于将待提取隐写分析特征图像通过10个滤波器进行滤波,得到残差
    图像;其中,10个滤波器包括4个有方向性滤波器和6个无方向性滤波器;

    编码???,用于根据LBPriu2模型和/或DLBP模型对残差图像进行编码后得到编码
    图,并根据编码图得到4个方向的二阶共生矩阵,对其中6个有方向性滤波器,相加合并不同
    方向计算得到的共生矩阵;

    共生矩阵特征提取???,用于对计算得到的共生矩阵,合并不同方向共生矩阵并
    降维后进行LOG函数映射得到共生矩阵特征。

    所述数字图像隐写分析特征提取方法,其中,所述编码??橹蠰BPriu2模型的计算
    式为:


    其中


    其中

    其中,gp是相邻像素值,gc是中间像素值,P、R是LBP编码中定义的相邻像素的个数
    和位置,取P=8和R=1。

    所述数字图像隐写分析特征提取方法,其中,所述编码??橹蠨LBP模型的计算式
    为:


    其中


    其中,其中gp是相邻像素值,gc是中间像素值,P、R是LBP编码中定义的相邻像素的
    个数和选取半径,取P=8和R=1,α根据经验设置成0.3和0.5,fi,j是滤波器系数。

    本发明所提供的数字图像隐写分析特征提取方法及提取系统,方法包括:将待提
    取隐写分析特征图像通过10个滤波器进行滤波,得到残差图像;其中,10个滤波器包括4个
    无方向性滤波器和6个非有方向性滤波器;根据LBPriu2模型和/或DLBP模型对残差图像进
    行编码后得到编码图,并根据编码图得到4个方向的二阶共生矩阵;对计算得到的共生矩
    阵,合并不同方向共生矩阵并降维后进行LOG函数映射得到共生矩阵特征。本发明与现有的
    主流通用隐写分析方法SRM具有相近的性能,分类效果较SRM具有更高的准确率。

    附图说明

    图1为本发明所述数字图像隐写分析特征提取方法较佳实施例的流程图。

    图2为6个相邻点选取范围尺度的示意图。

    图3为本发明所述数字图像隐写分析特征提取系统较佳实施例的功能结构框图。

    具体实施方式

    本发明提供一种数字图像隐写分析特征提取方法及提取系统,为使本发明的目
    的、技术方案及效果更加清楚、明确,以下对本发明进一步详细说明。应当理解,此处所描述
    的具体实施例仅仅用以解释本发明,并不用于限定本发明。

    如图1所示,为本发明所述一种数字图像隐写分析特征提取方法实施例的流程图,
    所述方法包括以下步骤:

    步骤S100、将待提取隐写分析特征图像通过10个滤波器进行滤波,得到残差图像;
    其中,10个滤波器包括4个无方向性滤波器和6个有方向性滤波器。

    对于4个无方向性滤波器,每个滤波器可以得到一个残差图像,对于6个有方向性
    滤波器,每个滤波器分别可以得到4个方向残差图像、1个最小值残差图像和1个最大值残差
    图像。

    步骤S200、根据LBPriu2模型和/或DLBP模型对残差图像进行编码后得到编码图,
    并根据编码图得到4个方向的二阶共生矩阵,对其中6个有方向性滤波器,相加合并不同方
    向计算得到的共生矩阵。

    步骤S300、对计算得到的共生矩阵,合并不同方向共生矩阵并降维后进行LOG函数
    映射得到共生矩阵特征。

    本发明的实施例中,对于待提取隐写分析特征图像(即给定的一副数字灰度图像
    I),首先需要通过不同的滤波器得到许多的余量R。对于滤波器的选择上,这里通过导数的
    方式得到四个基本滤波器,四个基本滤波器的设计较为规律易懂而且多样性适中。

    SRM特征充分考虑了图像像素间相关性,但是对于现在发展较为迅速的自适应隐
    写分析方法重点隐写区域纹理噪声区域没有特别的策略。而传统LBP模型对于图像纹理较
    为敏感,能够充分识别图像纹理特征,但是隐写对图像的改变十分细微,原始定义利用相对
    关系得出的模型对隐写并不是十分适用,因此本专利中公开了根据相对距离的DLBP模型,
    将变化幅度考虑进来,增加对隐写的敏感程度。

    对于共生矩阵,由于其是对图像给定方向相邻像素间关系的统计,它有一定稀疏
    性以及最大最小值相差几个数量级,这不利于集成分类器训练分类,因此本发明中将共生
    矩阵特征共LOG函数进行映射,使特征值能在同一个数量级里。

    优选的,在所述数字图像隐写分析特征提取方法中,所述步骤S100中10个滤波器
    所包括的4个基本滤波器分别记为D12、D13、D14和D15,其中:

    [1,0]-[0,1]=[1,-1] (D12)

    [1,-1,0]-[0,1,-1]=[1,-2,1] (D13)

    [1,-2,1,0]-[0,1,-2,1]=[1,-3,3,-1] (D14)

    [1,-3,3,-1,0]-[0,1,-3,3,-1]=[1,-4,6,-4,1] (D15);

    10个滤波器所包括的6个由基本滤波器推导出的滤波器分别记为D22、D33、D44、
    D55、D23和D35,其中:





    对于D23和D35两个滤波器,需要将其用零填充成3*3和5*5正方形矩阵,在D23下部
    填充1*3的零矩阵,在D35下部填充2*5的零矩阵。

    D22、D33、D44、D55是的4个无方向性滤波器,D12、D13、D14、D15、D23、D35是6个有方
    向性滤波器,对有方向性滤波器,分别计算他们各个方向上的余量以及计算其中各个方向
    余量的最大和最小值余量。

    所述数字图像隐写分析特征提取方法,其中,所述步骤S100中通过10个滤波器进
    行滤波得到残差图像的计算式为:



    其中,I为待提取隐写分析特征图像,f是滤波器,是卷积符号,Row(f)和Colum
    (f)分别为求f的行和列数。

    对于以上10个滤波器,其中因为D22、D33、D44、D55已经具有中心对称性,不需要旋
    转,所以各自可以得到一个残差图像;对于D12、D13、D14、D15、D23、D35六个滤波器,将他们
    以90度旋转4次得到4个残差图像,再由各自4个方向的残差图像取最大值和最小值可以得
    到两个非线性残差图像(取名为min和max),因此分别可以得到6个残差图像。因此通过以上
    10个滤波器总共可以40个残差图像。

    优选的,在所述数字图像隐写分析特征提取方法中,所述步骤S200中LBPriu2模型
    的计算式为:


    其中


    其中

    其中,gp是相邻像素值,gc是中间像素值,P、R是LBP编码中定义的相邻像素的个数
    和位置,取P=8和R=1。

    这种模型根据旋转不变特性和LBP编码零和一变换次数两个原则将原始LBP模型
    从256种降低到10种,这有利于后续共生矩阵转换维数控制。由以上几个操作,首先将残差
    图像量化截断,之后进行LBPriu2模型编码可以得到具有10个取值范围的LBPriu2编码图
    RLBP。

    优选的,在所述数字图像隐写分析特征提取方法中,所述步骤S200中DLBP模型的
    计算式为:


    其中


    其中,其中gp是相邻像素值,gc是中间像素值,P、R是LBP编码中定义的相邻像素的
    个数和选取半径,取P=8和R=1,α根据经验设置成0.3和0.5,fi,j是滤波器系数。

    DLBP这里通过编码之后具有256种模式,之后线性将相邻32种模式归为一种模式,
    以此将256种降为8种。对于DLBP模型编码在选取相邻像素点的选择上,本发明中选择选取
    相邻八个像素的点作为LBP编码参考点,另外采取一个已有的多种范围尺度取相邻点的方
    法,即同样取中心像素附近八个相邻点作为模型编码参考点,但是不只是选取相邻的像素,
    如图2所示,由此可以得到6个选取范围尺度,跟LBP编码模型进行相同操作。根据以上操作,
    进行DLBP模型编码可以得到具有8个取值范围的DLBP编码图RDLBP。

    优选的,在所述数字图像隐写分析特征提取方法中,所述步骤S200中4个方向的二
    阶共生矩阵的生成式为:

    C0(i,j)=#{(x1,y1),(x1,y1+1)∈M×N|RLBP(x1,y1)=i,RLBP(x1,y1+1)=j}

    C45(i,j)=#{(x1,y1),(x1-1,y1+1)∈M×N|RLBP(x1,y1)=i,RLBP(x1-1,y1+1)=j}

    C90(i,j)=#{(x1,y1),(x1-1,y1)∈M×N|RLBP(x1,y1)=i,RLBP(x1-1,y1)=j}

    C135(i,j)=#{(x1,y1),(x1-1,y1+1)∈M×N|RLBP(x1,y1)=i,RLBP(x1-1,y1+1)=
    j};

    其中,C0为0°方向上二阶共生矩阵、C45为45°方向上二阶共生矩阵、C90为90°方向上
    二阶共生矩阵、C135为135°方向上二阶共生矩阵,M、N是编码图的长和宽,RLBP是LBP或者DLBP
    编码图。

    具体的,在所述数字图像隐写分析特征提取方法中,所述步骤S300中具体包括:

    步骤S301、对滤波器进行滤波得到的残差图像,计算共生矩阵;

    步骤S302、将共生矩阵沿左上到右下对角线对折相加,并拉伸成为0°方向上的一
    维向量Cf0、45°方向上的一维向量Cf45、90°方向上的一维向量Cf90、及135°方向上的一维向
    量Cf135;

    步骤S303、将0°方向上的一维向量Cf0与90°方向上的一维向量Cf90相加得到
    Cf0+90,将45°方向上的一维向量Cf45与135°方向上的一维向量Cf135相加得到Cf45+135;

    步骤S304、分别计算f1=log10(Cf0+90+1)和f2=log10(Cf45+135+1),将f1和f2串联后
    得到共生矩阵特征。

    对于x角度方向上的共生矩阵Cx,本专利将其沿左上到右下对角线对折相加,对角
    线上像素不做处理,之后拉伸成一维向量Cfx。对一个余量图生成的4个方向一维向量,将
    Cf0和Cf90相加、Cf45和Cf135相加得到Cf0+90和Cf45+135,之后分别求f1=log10(Cf0+90+1)和f2
    =log10(Cf45+135+1),将f1和f2串联起来即可以得到该余量图在某个范围尺度下的共生矩
    阵变变换特征。由以上步骤一个残差图像,可以得到的LBP共生矩阵特征维数是660*1维,得
    到的DLBP共生矩阵特征维数是432*1维。

    对于D12、D13、D14、D15、D23、D35这六个具有方向性的滤波器,可以分别得到四个
    方向以及最小值和最大值共六个残差图像,本方案中对其中的四个方向余量图在分别得到
    的共生矩阵之后,进行相加合并成一个共生矩阵,以此再进行余下操作,通过这种方法进行
    降维合并,而最小值和最大值余量图则分别得到各自的共生矩阵特征不做另外的合并操
    作。最小值和最大值余量的DLBP共生矩阵特征性能不理想,因此只用LBP共生矩阵特征。因
    此这六个滤波器D12、D13、D14、D15、D23、D35总共有(660*2+432*2)*6+660*2*2*6=28944
    维;对于D22、D33、D44、D55这四个滤波器,总共有(660*2+432*2)*4=8736维。因此一幅数字
    图像总共可以得到37680维特征。

    基于上述方法实施例,本发明还提供一种数字图像隐写分析特征提取系统。如图3
    所示,所述数字图像隐写分析特征提取系统包括:

    滤波???00,用于将待提取隐写分析特征图像通过10个滤波器进行滤波,得到残
    差图像;其中,10个滤波器包括4个无方向性滤波器和6个有方向性滤波器;

    编码???00,用于根据LBPriu2模型和/或DLBP模型对残差图像进行编码后得到
    编码图,并根据编码图得到4个方向的二阶共生矩阵,对其中6个有方向性滤波器,相加合并
    不同方向计算得到的共生矩阵;

    共生矩阵特征提取???00,用于对计算得到的共生矩阵,合并不同方向共生矩阵
    并降维后进行LOG函数映射得到共生矩阵特征。

    优选的,在所述数字图像隐写分析特征提取方法中,所述编码???00中LBPriu2
    模型的计算式为:


    其中


    其中

    其中,gp是相邻像素值,gc是中间像素值,P、R是LBP编码中定义的相邻像素的个数
    和位置,取P=8和R=1。

    优选的,在所述数字图像隐写分析特征提取方法中,所述编码???00中DLBP模型
    的计算式为:


    其中


    其中,其中gp是相邻像素值,gc是中间像素值,P、R是LBP编码中定义的相邻像素的
    个数和选取半径,取P=8和R=1,α根据经验设置成0.3和0.5,fi,j是滤波器系数。

    综上所述,本发明所提供的数字图像隐写分析特征提取方法及提取系统,方法包
    括:将待提取隐写分析特征图像通过10个滤波器进行滤波,得到残差图像;其中,10个滤波
    器包括4个无方向性滤波器和6个有方向性滤波器;根据LBPriu2模型和/或DLBP模型对残差
    图像进行编码后得到编码图,并根据编码图得到4个方向的二阶共生矩阵;对滤波得到的残
    差图像,计算共生矩阵,合并后进行LOG函数映射得到共生矩阵特征。本发明与现有的主流
    通用隐写分析方法SRM具有相近的性能,分类效果较SRM具有更高的准确率。

    应当理解的是,本发明的应用不限于上述的举例,对本领域普通技术人员来说,可
    以根据上述说明加以改进或变换,所有这些改进和变换都应属于本发明所附权利要求的保
    护范围。

    关于本文
    本文标题:一种数字图像隐写分析特征提取方法及提取系统.pdf
    链接地址://www.4mum.com.cn/p-6079461.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 大乐透走势图表 即开彩官网 北京pk10靠谱刷水方案 飞艇计划人工在线app 一分快三怎么玩稳赚 极速28稳赚 12个码复式二中二多少组 二分彩技巧稳赚保本玩法 7八卦网正版一肖包平特 福彩3d组6复试 3d5码组六遗漏统计表 北京pk10稳赚技巧方案 重庆时时开奖官方同步 辽宁12选5技巧稳赚 秒速时时开奖规律 必赢客北京pk拾