发明名称 一种用于对岩屑荧光图像进行快速分析方法
摘要 本发明涉及一种用于对岩屑荧光图像进行快速分析方法,其包括以下步骤:1)基于RGB三维颜色空间进行油性组分识别训练,获得聚类文件;2)根据步骤1)得到的聚类文件对待测的岩屑荧光图像进行油性组分分析,具体包括以下步骤:2.1)对输入的待测的岩屑荧光图像进行颜色特征提取,并生成特征向量数组;2.2)对特征向量数组中的每一个特征向量进行分析;2.3)直到完成向量数组中所有特征向量的分析工作,计算出岩屑荧光图像中所含油性组分以及各组分所占比重。本发明的基于非均匀颜色空间色差计算的岩屑荧光分层聚类方法在油性组分识别和岩屑荧光图像分析结果精确、速度快,能够满足海上平台地质人员智能化录井需要,因此可以广泛应用在地质勘探开发过程中。
申请公布号 CN103927547B 申请公布日期 2017.04.19
申请号 CN201410125638.8 申请日期 2014.03.31
申请人 中海油能源发展股份有限公司;中海油信息科技有限公司 发明人 符耀庆;刘耀华;王亚楠;苑舒斌;麦文;苏文辉;赵伦;刘换来;晏菲;陈金定;杨宝伟;沈雪峰
分类号 G06K9/62(2006.01)I;G06T7/11(2017.01)I;E21B49/00(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 北京纪凯知识产权代理有限公司 11245 代理人 徐宁;关畅
主权项 一种用于对岩屑荧光图像进行快速分析方法,其包括以下步骤:1)基于RGB三维颜色空间进行油性组分识别训练,获得聚类文件,具体过程为:1.1)将RGB三维颜色空间转换成256个RGB二维平面,具体过程为:计算出RGB三维颜色空间中每一个颜色的饱和度值,饱和度值的取值范围为0~255之间的整数,然后根据不同颜色的饱和度值将RGB三维颜色空间分成256个RGB二维平面,并对每一个RGB二维平面进行编号,此编号就是每一个RGB二维平面所对应的饱和度值;1.2)基于RGB二维平面进行相近颜色区域分割,即对每一个RGB二维平面进行聚类分析,并根据每一个RGB二维平面的聚类分析结果,在相应RGB二维平面上绘制区域边界,具体过程为:1.2.1)读取某个RGB二维平面,并提取该RGB二维平面中每个像素的特征向量,其中,像素集S={X<sub>1</sub>,X<sub>2</sub>,…,X<sub>m</sub>},即此RGB二维平面中所有像素的集合;1.2.2)采用区域生长方法对上述RGB二维平面进行聚类分析,具体过程为:A)定义聚类半径R;B)对像素集S中的所有颜色逐一进行聚类处理:对某一个二维平面进行聚类分析时,首先选择像素集S中第一个像素X<sub>1</sub>作为第一个待分割区域的种子点,由于此时聚类中心集为空,所以将该种子点记作聚类中心O<sub>1</sub>,并将X<sub>1</sub>添加到O<sub>1</sub>对应的聚类集C<sub>1</sub>中;然后按顺序从像素集中取出未分类像素,记作X<sub>a</sub>,并采用非均匀颜色空间色差计算方法计算该像素到聚类中心集中所有聚类中心的距离:<maths num="0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>D</mi><mo>=</mo><mo>{</mo><msub><mi>D</mi><mn>1</mn></msub><mo>,</mo><msub><mi>D</mi><mn>2</mn></msub><mo>,</mo><mn>...</mn><mo>,</mo><msub><mi>D</mi><mi>n</mi></msub><mo>}</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>D</mi><mi>i</mi></msub><mo>=</mo><msqrt><mrow><mo>(</mo><msub><mi>k</mi><mi>r</mi></msub><mo>*</mo><msup><mrow><mo>(</mo><mrow><msub><mi>r</mi><mn>1</mn></msub><mo>-</mo><msub><mi>r</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>k</mi><mi>g</mi></msub><mo>*</mo><msup><mrow><mo>(</mo><mrow><msub><mi>g</mi><mn>1</mn></msub><mo>-</mo><msub><mi>g</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>k</mi><mi>b</mi></msub><mo>*</mo><msup><mrow><mo>(</mo><mrow><msub><mi>b</mi><mn>1</mn></msub><mo>-</mo><msub><mi>b</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo><mo>/</mo><mn>4</mn><mo>*</mo><msup><mn>255</mn><mn>2</mn></msup><mo>+</mo><mi>&theta;</mi></mrow></msqrt></mrow></mtd></mtr></mtable><mo>,</mo></mrow>]]></math><img file="FDA0001204133180000011.GIF" wi="1711" he="183" /></maths>式中,(r<sub>1</sub>,g<sub>1</sub>,b<sub>1</sub>)表示颜色样本X<sub>a</sub>颜色特征,(r<sub>2</sub>,g<sub>2</sub>,b<sub>2</sub>)表示聚类中心O<sub>i</sub>颜色特征,(k<sub>r</sub>,k<sub>g</sub>,k<sub>b</sub>)表示像素与聚类中心空间距离的加权系数,θ表示修正后的像素与聚类中心的角度距离;C)取D中的最小值D<sub>i</sub>,并比较D<sub>i</sub>与聚类半径R的大小,如果D<sub>i</sub>≤R,则返回步骤B);否则进入步骤F);D)将像素X<sub>a</sub>添加到聚类中心O<sub>i</sub>对应的聚类集C<sub>i</sub>中;E)更新聚类中心O<sub>i</sub>,使其始终处于聚类集C<sub>i</sub>的中心位置,更新完成后跳转到步骤B):<maths num="0002"><math><![CDATA[<mrow><msub><mi>O</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>K</mi><mi>i</mi></msub><mo>*</mo><msub><mi>O</mi><mi>i</mi></msub><mo>+</mo><msub><mi>X</mi><mi>a</mi></msub></mrow><mrow><msub><mi>K</mi><mi>i</mi></msub><mo>+</mo><mn>1</mn></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0001204133180000012.GIF" wi="316" he="102" /></maths>式中,K<sub>i</sub>表示聚类中心O<sub>i</sub>对应的聚类集C<sub>i</sub>中像素的个数;F)生成新的聚类,并将X<sub>a</sub>作为新的聚类中心O<sub>new</sub>添加到聚类中心集O中,同时创建该聚类中心所对应的聚类集C<sub>new</sub>,并将X<sub>a</sub>添加进去,进入步骤D);1.2.3)检查是否已对像素集中所有的像素进行了聚类处理,如果处理未完成返回步骤B)继续进行处理,否则进入步骤1.2.4);1.2.4)根据不同聚类中心所对应的聚类集,寻找区域边界,并将边界绘制在RGB二维平面上;1.3)当所有RGB二维平面分割完成后,为每一个分割完成的RGB二维平面的所有的颜色区域标定油性组分;1.4)对每一个标定完成的RGB二维平面创建相应聚类文件,根据每一RGB二维平面的饱和度值为其命名,同时将每一RGB二维平面的聚类结果以及相应分析结果保存到相应聚类文件中;2)根据步骤1)得到的聚类文件对待测的岩屑荧光图像进行油性组分分析,具体包括以下步骤:2.1)对输入的待测的岩屑荧光图像进行颜色特征提取,并生成特征向量数组;2.2)对特征向量数组中的每一个特征向量进行分析;2.3)完成向量数组中所有特征向量的分析工作,计算出岩屑荧光图像中所含油性组分以及各组分所占比重。
地址 100027 北京市东城区东直门外小街6号