CN109947394B - 一种基于三角形覆盖的图像二维分形维数计算方法 - Google Patents
一种基于三角形覆盖的图像二维分形维数计算方法 Download PDFInfo
- Publication number
- CN109947394B CN109947394B CN201910180449.3A CN201910180449A CN109947394B CN 109947394 B CN109947394 B CN 109947394B CN 201910180449 A CN201910180449 A CN 201910180449A CN 109947394 B CN109947394 B CN 109947394B
- Authority
- CN
- China
- Prior art keywords
- image
- triangular
- module
- square
- coverage
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title abstract description 10
- 238000000034 method Methods 0.000 claims description 13
- 238000004220 aggregation Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Abstract
本发明公开了一种基于三角形覆盖的图像二维分形维数计算方法,该方法用三角形模块对边长为2的整数次幂的正方形图像进行覆盖,计算其二维分形维数。本发明证明了三角形模块覆盖形式能够充分并精确地包含图像的细节信息,三角形模块覆盖形式计算出的图像维数相较常用的正方形模块覆盖形式计算出的图像维数更接近图像的理论hasudorff维数。因此三角形模块覆盖形式可以使得图像的分维分析更加准确。
Description
技术领域
本发明属于图像识别技术领域,具体涉及一种基于三角形覆盖的图像二维分形维数的计算方法,适用于以图像的分形维数作为特征值的图像识别研究。
背景技术
盒子法是一种经典的二维分形维数计算方法,它的特征是在空间域内用不同边长R的正方形模块覆盖边长为2的整数次幂的二维正方形图像,得到一系列不同尺度R下的测度N,将这一系列尺度R与测度N带入相关公式中进行运算,即可得出图像的二维分形维数,从而实现对图像二维分形特性-分形维数的计算。
常用的盒维法计算图像分形维数时,一般以正方形模块的形式覆盖图像。使用正方形模块使分形维数计算简便,但是该方法容易形成过度覆盖,使得图像目标区域的轮廓占据的正方形模块区域不够充分。
发明内容
为了解决上述技术问题,本发明提供了一种基于三角形覆盖的图像二维分形维数的计算方法。
本发明所采用的技术方案是:一种基于三角形覆盖的图像二维分形维数计算方法,其特征在于:用三角形模块代替正方形模块对边长为2的整数次幂的正方形图像进行覆盖,并计算其二维分形维数。
本发明提供的三角形模块覆盖形式可以充分并精确地包含图像的细节信息,使得图像的分维分析更加准确。
附图说明
图1为本发明实施例的基于三角形覆盖的结构示意图,其中(a)、(b)、(c)、(d)分别表示的为模块1、2、3、4,模块1是边长为2的整数次幂正方形模块右上方处的一个三角模块,占正方形模块面积的1/4。将模块1分别顺时针旋转90°,180°与270°即可得模块2、3与4。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明,下面结合附图及实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
针对盒子法中用正方形模块覆盖边长为2的整数次幂图像时可能存在过度覆盖的问题,本发明提出了一个技术方案是:用三角形模块(▲)代替正方形模块对边长为2的整数次幂的正方形图像进行覆盖,并计算其二维分形维数。
请见图1,本发明提供的一种基于三角形覆盖的图像二维分形维数计算方法,具体实现包括以下步骤:
步骤1:在一幅大小为M*M的图像中,可将任意一个边长为r的正方形模块分成4个尺寸相同方向不同的三角形模块,如图1(a)所示,包含“1”的阴影区域即为覆盖图像的三角形模块,其中h为三角形模块的高度,h等于r/2。
步骤2:使用这四种三角形模块对图像进行覆盖,若图像中模块阴影(包含1值)的三角形区域内灰度值存在不为0的像素点,则计为1个模块,统计四种模块的高度h与覆盖整幅图像所需的模块数即测度Nh;
步骤3:改变四个三角形模块的高h(h为2的整数次幂实数,取值范围为2至M/2),重复步骤2,获得一组不同尺度h下的测度Nh,对1/h和Nh取对数得到log(1/h)和log(Nh);
步骤4:计算基于三角形二维盒维数D;
以下对三角形覆盖的二维分形维数计算方法的有效性进行证明。
豪斯多夫测度的定义:U为n维欧式几何空间的Rn非空子集,U的直径定义为:
式(1)中x与y分别为U集合内的任意两点,sup为上确界(supremum)的缩写。因此,|U|即U集合内任意两点间距离的最大值。
如果{Ui}为有限个直径不超过δ的集构成的覆盖A(在n维欧式空间Rn)的集类,即:
且对任意的i值有0<|Ui|≤δ,则称n维欧氏空间Rn的{Ui}是A的一个δ-覆盖。
Hausdorff测度定义如下:假设s为非负数,对任意集合A为Rn中的任意子集,δ>0,令
式(3)中inf为下确界(infimum)的缩写。表示对A中所有的δ-覆盖取下确界。考察所有直径不超过δ时A的覆盖,并使这些直径的s次幂的和达到最小值,注式(3)中单调非减,且当δ→0时趋于一个极限值,记为。
对于任意Rn中的子集A这个极限都存在。称Hs(A)为A的s-维Hausdorff测度。此外,E、F为Rn中δ覆盖的任意子集,若E、F间的距离d(E,F)>δ(E、F为Rn的子集),它隐含则在E∪F的δ覆盖中没有一个集合能同时与E和F相交,故即有
Hs(E∪F)=Hs(E)+Hs(F) (5)
由于集合所以,集合B可以表示为B=A∪Δ1,且A∩Δ1=φ,其中,Δ1为▲覆盖集合A时过度覆盖的集合;集合C可表示为C=B∪Δ2,其中,Δ2为覆盖集合B时过度覆盖的集合。即C=A∪Δ1∪Δ2,且B∩Δ2=φ,A∩Δ1∩Δ2=φ。
所以,用▲覆盖集合A可以表示为:
B=A∪Δ1 (6)
根据式(5)可得:
即
C=B∪Δ2=A∪Δ1∪Δ2 (9)
根据式(5)可得:
即:
由此可以证明:▲方法计算A的维数Dt出更接近A的理论hasudorff维数。
故,当分别使用三角形模块(▲)和正方形模块对二维图像进行覆盖时(假设,是由4个▲组成)可推导出以下计算定理:假设集合A∈U,用▲和分别对A进行覆盖,可分别计算出集合A的分形维数Dt和Ds,▲方法计算的A的分形维数Dt相较于方法计算的A的分形维数Ds更接近A的理论Hasudorff维数Df。
本实施例证明了三角形模块覆盖形式能够充分并精确地包含图像的细节信息,三角形模块覆盖形式计算出的图像维数Dt相较正方形模块覆盖形式计算出的图像维数Ds更接近图像的理论hasudorff维数Df。因此三角形模块覆盖形式可以使得图像的分维分析更加准确。
应当理解的是,本说明书未详细阐述的部分均属于现有技术。
应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。
Claims (1)
1.一种基于三角形覆盖的图像二维分形维数计算方法,其特征在于:用三角形模块代替正方形模块对边长为2的整数次幂的正方形图像进行覆盖,并计算其二维分形维数;
具体实现包括以下步骤:
步骤1:在一幅大小为M*M的图像中,将任意一个边长为r的正方形模块分成4个尺寸相同方向不同的三角形模块,其中h为三角形模块的高度,h等于r/2;
步骤2:使用这四种三角形模块对图像进行覆盖,若图像中模块阴影的三角形区域存在灰度值不为0的像素点,则计为1个模块,统计四种模块的高度h与覆盖整幅图像所需的模块数即测度Nh;
步骤3:改变四种三角形模块的高度h,h为2的整数次幂实数,取值范围为2至M/2;重复步骤2,获得一组不同尺度h下的测度Nh,对1/h和Nh取对数得到log(1/h)和log(Nh);
步骤4:计算基于三角形二维盒维数D;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910180449.3A CN109947394B (zh) | 2019-03-11 | 2019-03-11 | 一种基于三角形覆盖的图像二维分形维数计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910180449.3A CN109947394B (zh) | 2019-03-11 | 2019-03-11 | 一种基于三角形覆盖的图像二维分形维数计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109947394A CN109947394A (zh) | 2019-06-28 |
CN109947394B true CN109947394B (zh) | 2023-03-10 |
Family
ID=67009543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910180449.3A Active CN109947394B (zh) | 2019-03-11 | 2019-03-11 | 一种基于三角形覆盖的图像二维分形维数计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109947394B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110751657B (zh) * | 2019-09-26 | 2023-05-02 | 湖北工业大学 | 一种基于三角形覆盖的图像三维分形维数计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101504727A (zh) * | 2008-01-31 | 2009-08-12 | 株式会社东芝 | 图像处理设备及其方法 |
CN103208003A (zh) * | 2013-04-12 | 2013-07-17 | 大连理工大学 | 一种基于几何图形特征点形状描述子的方法 |
CN108564609A (zh) * | 2018-04-23 | 2018-09-21 | 大连理工大学 | 一种基于盒维数法的计算分形维数的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2905637C (en) * | 2013-03-13 | 2022-04-05 | Fdna Inc. | Systems, methods, and computer-readable media for identifying when a subject is likely to be affected by a medical condition |
-
2019
- 2019-03-11 CN CN201910180449.3A patent/CN109947394B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101504727A (zh) * | 2008-01-31 | 2009-08-12 | 株式会社东芝 | 图像处理设备及其方法 |
CN103208003A (zh) * | 2013-04-12 | 2013-07-17 | 大连理工大学 | 一种基于几何图形特征点形状描述子的方法 |
CN108564609A (zh) * | 2018-04-23 | 2018-09-21 | 大连理工大学 | 一种基于盒维数法的计算分形维数的方法 |
Non-Patent Citations (4)
Title |
---|
一类Sierpinski垫的Hausdorff测度;王明华;《数学研究与评论》;20071115(第04期);全文 * |
二维数字图像分形维数的计算方法;彭瑞东等;《中国矿业大学学报》;20040130(第01期);全文 * |
李业学等.粗糙面分形计算理论研究进展.《数学进展》.2012,(第04期), * |
粗糙面分形计算理论研究进展;李业学等;《数学进展》;20120815(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109947394A (zh) | 2019-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109255811B (zh) | 一种基于可信度图视差优化的立体匹配方法 | |
WO2022077863A1 (zh) | 视觉定位方法及相关模型的训练方法及相关装置、设备 | |
US10467777B2 (en) | Texture modeling of image data | |
CN105184786B (zh) | 一种浮点型三角形特征描述方法 | |
CN109947394B (zh) | 一种基于三角形覆盖的图像二维分形维数计算方法 | |
CN109785221A (zh) | 一种数字图像隐写方法及秘密信息提取方法 | |
CN116129037B (zh) | 视触觉传感器及其三维重建方法、系统、设备及存储介质 | |
CN111553296B (zh) | 一种基于fpga实现的二值神经网络立体视觉匹配方法 | |
CN103093226B (zh) | 一种用于图像特征处理的ratmic描述子的构造方法 | |
CN115564926A (zh) | 基于影像建筑物结构学习的三维面片模型构建方法 | |
CN105469392A (zh) | 基于区域光谱梯度特征对比的高光谱图像显著性检测方法 | |
CN113409332B (zh) | 一种基于三维点云的建筑物平面分割方法 | |
JP4661879B2 (ja) | 変換テーブル作成装置、変換テーブル作成プログラム及び変換テーブル作成方法 | |
CN116505953B (zh) | 基于bim和gis的海量地图数据优化压缩处理方法 | |
US11348261B2 (en) | Method for processing three-dimensional point cloud data | |
CN117392237A (zh) | 一种鲁棒的激光雷达-相机自标定方法 | |
CN104112278A (zh) | 基于协方差的多光谱图像实时配准方法 | |
US20080129727A1 (en) | System and method for generating curvature adapted isosurface based on delaunay triangulation | |
CN116958434A (zh) | 多视图三维重建方法、测量方法及系统 | |
EP3742398A1 (en) | Determining one or more scanner positions in a point cloud | |
CN103489197A (zh) | 一种城市地区航空影像角特征匹配方法 | |
CN113592744B (zh) | 一种适用于高分辨率遥感影像的几何精校正方法 | |
CN110751657B (zh) | 一种基于三角形覆盖的图像三维分形维数计算方法 | |
CN111127667B (zh) | 基于区域曲率二进制描述符的点云初始配准方法 | |
CN115330619A (zh) | 适用于高分辨率遥感影像的局部几何精校正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20190628 Assignee: Hubei Huasheng Xinxin plate making Co.,Ltd. Assignor: HUBEI University OF TECHNOLOGY Contract record no.: X2023980053043 Denomination of invention: A Method for Calculating the 2D Fractal Dimension of Images Based on Triangle Coverage Granted publication date: 20230310 License type: Common License Record date: 20231220 |