二维数据的分形维数分析（C++实现）

最新推荐文章于 2024-06-01 17:50:41 发布

改不了名字的二狗咂

最新推荐文章于 2024-06-01 17:50:41 发布

阅读量488

点赞数

文章标签： c++

本文链接：https://blog.csdn.net/zzlhyqz/article/details/124446495

版权

分形维数

20 世纪 70 年代 Mandelbrot 首次提出分形方法，分形几何是非线性科学的交叉研究领域，广泛应用于流体力学、计算机图形图像处理、经济学、社会学、生物学、地质学等学科，分形的引入为科学家们提供了一种全新的视角去揭示自然界中的规律和特性。将分形方法引入地球科学领域，可以用于海岸线的描绘、海底地形图的成图、地震预测、地学图像及数字信号处理等方面。
本程序主要用来处理重力异常数据，以期得到图区的地质构造复杂程度。处理得到的分维数D在一定程度上反应了对应地区地质构造的复杂性和活动强度。D值越大，反应的复杂性越高，反之亦然。

算法

本程序主要基于乔计花等[1]人文章中描述的技术思路改写，初步窗口选择处采用陈歌[2]提到的数据分块处理方法，对选定好每个子矩阵内数据进行一下操作：（为保证计算的可行，取子矩阵窗口大小时长度需为奇数）
取自文献[1]
具体流程可描述为：（1）对观测数据扩编。（2）根据窗口大小将观测数据分块。（3）对窗口内数据进行傅里叶变换，求取对应等效半径（即半径波数Kj）的二维平均功率谱密度，两边取对数后最小二乘拟合得到beta和logR，再利用测度与分维数的关系得到这一窗口的分维数D。（4）将窗口向后移动，重复（3）的计算，直至所有数据计算完毕。
本程序不改变原始数据的大小，输入（MxN）的数据，计算后得到同样大小（MxN）的分维数D结果。