[发明专利]一种基于DFT的土壤表面起伏特征分形数据生成算法

说明书

本发明公开了一种基于DFT的土壤表面起伏特征分形数据生成算法，以离散傅里叶变换为基础的分形数据生成算法，用于生成表征土壤表面起伏特征的高度数据。本发明的有益效果是：可实现真实土壤起伏表面建模，客观再现电磁波与真实土壤表面的相互作用过程，缩小仿真结果与实测结果之间的差异，实现电磁波无损探测结果的准确分析和解释。

技术领域

本发明涉及一种土壤表面起伏特征分形数据生成算法，具体为一种基于DFT的土壤表面起伏特征分形数据生成算法，属于地球物理勘探领域电磁波无损探测方面的数值仿真技术领域。

背景技术

土壤是指地球表面的一层疏松的物质，土壤颗粒在自然条件下是聚集在一起以土壤结构的形式表现出来。土壤表面起伏特征的复杂程度可用分形维数进行衡量，分形在这个世界是普遍存在的，分型可以归纳出两个特征。第一，分型具有自相似性及分型，自身可以看成是有许多与自己相似的大小不一的部分组成；第二，分形具有无穷多的层次，及无论在分型的哪一个层次，总能看到有更精细的下一个层次存在，或者说分形的图形有无限的细节，可以不断的放大，永远都有结构。分形维数主要描述分形最主要的参量，简称分维。分维反映了复杂形体占有空间的有效性，它是复杂形体不规则性的量度。

真实土壤表面呈现凸凹不平的起伏特征，通常会引起反射电磁波相当程度的抖动。目前，常见的电磁场仿真物理模型，通常采用标准平面模拟复杂起伏的土壤表面。由于表面结构过于简单、理想，难以再现电磁波与真实土壤表面的相互作用过程，仿真结果与实测结果偏差较大，严重影响电磁波无损探测结果的分析和解释。

发明内容

本发明的目的就在于为了解决上述问题而提供的通过构建土壤表面凸凹不平的起伏特征，形成拟真的土壤表面形态的一种基于DFT的土壤表面起伏特征分形数据生成算法。

本发明通过以下技术方案来实现上述目的：一种基于DFT的土壤表面起伏特征分形数据生成算法，包括以下步骤：

步骤1、生成大小为x×y的随机二维矩阵A；

步骤2、对矩阵A进行二维快速傅里叶变换，并将结果进行平移，使零频率处于二维矩阵的中心；

步骤3、设矩阵A的中心点坐标为(ic,jc)，对于每一个数据点A(i,j)，计算其与中心点的模L；即

步骤4、根据近似幂律关系对矩阵A进行变换；采用公式为：

步骤5、再次对矩阵A进行逆傅里叶变换，将所得结果归一化到[-k,k]的区间范围内，此时矩阵A即代表土壤表面的起伏程度，则该数据可直接用于土壤表面的建模，实现具有分形特征的逼真土壤表面形态模型。

优选的，所述步骤1中，生成二维矩阵A的x和y值皆为400，且生成的矩阵为随机矩阵。

优选的，所述步骤4中，对分形特征数据进行离散傅里叶变换后，其频率f与相应的频率系数mf近似呈幂律关系，即

其中，C为常数，不影响土壤表面起伏特征的衡量；D为分形维度，其值是介于1到3之间的实数。

优选的，所述步骤4中，频域数据经过变换公式处理后能量集中到低频区域。

优选的，所述步骤5中，对矩阵A进行逆傅里叶变换之前需要再次对A进行平移操作。

优选的，所述步骤5中，生成的土壤起伏数据不包含海平面信息，且实际建模时，需要与模型中土壤的基准高度相加，实际的土壤表面高度坐标为：