CN107045435A - 基于梯度插值的三维相干随机数生成方法 - Google Patents

基于梯度插值的三维相干随机数生成方法 Download PDF

Info

Publication number
CN107045435A
CN107045435A CN201710236412.9A CN201710236412A CN107045435A CN 107045435 A CN107045435 A CN 107045435A CN 201710236412 A CN201710236412 A CN 201710236412A CN 107045435 A CN107045435 A CN 107045435A
Authority
CN
China
Prior art keywords
interpolation
point
gradient
random number
interpolation point
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.)
Granted
Application number
CN201710236412.9A
Other languages
English (en)
Other versions
CN107045435B (zh
Inventor
朱国庆
董纯柱
殷红成
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Environmental Features
Original Assignee
Beijing Institute of Environmental Features
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Environmental Features filed Critical Beijing Institute of Environmental Features
Priority to CN201710236412.9A priority Critical patent/CN107045435B/zh
Publication of CN107045435A publication Critical patent/CN107045435A/zh
Application granted granted Critical
Publication of CN107045435B publication Critical patent/CN107045435B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
    • G06F7/58—Random or pseudo-random number generators
    • G06F7/588—Random number generators, i.e. based on natural stochastic processes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40—Means for monitoring or calibrating
    • G01S7/4052—Means for monitoring or calibrating by simulation of echoes

Landscapes

  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Complex Calculations (AREA)
  • Image Generation (AREA)

Abstract

公开了基于梯度插值的三维相干随机数生成方法，包括：S1、获取待匹配点周围的N个插值点，为该插值点生成非相干随机数；基于所述非相干随机数查询预设的映射关系，确定对应插值点的梯度矢量；S2、获取每个插值点的坐标、以及所述N个插值点所在几何体的中心点坐标；根据插值点坐标、中心点坐标和每个插值点的梯度矢量，确定每个插值点的线性函数；S3、采用三线性插值法对每个插值点进行插值处理，得到待匹配点的相干随机数。本发明基于待匹配点周围的N个插值点进行三线性梯度插值处理，能够降低插值处理维数，能够快速生成待匹配点的相干随机数，易于实现，计算量小。

Description

基于梯度插值的三维相干随机数生成方法
技术领域
本发明涉及雷达回波仿真研究领域，尤其涉及基于梯度插值的三维相干随机数生成方法。
背景技术
以下对本发明的相关技术背景进行说明，但这些说明并不一定构成本发明的现有技术。
由于能很好地描述地面环境的起伏特性和雷达探测环境回波的起伏特性，相干随机数广泛应用于各种地貌模型构建和地杂波数值模拟。与产生非相干随机数相比，产生相干随机数具有以下几点不同特性：1)输入相同的参数总会输出相同的随机数值；2)输入较小改变量的参数总会输出一个较小改变量的随机数值；3)输入较大改变量的参数会输出一个随机变化量的数值。一个一维的相干随机数序列和非相干随机数序列如图1a和图1b所示。。
目前，相干随机数生成应用最多的有两种方法，一种是零记忆非线性变换方法(ZMNL)，另一种是球不变随机过程法(SIRP)。ZMNL先生成具有特定空间相关性的标准高斯分布随机数，然后再通过非线性变换将高斯分布随机数变换到特定分布的随机数，它易于实现，并且在相关Gaussian序列产生后仿真速度较快，缺点是难以找到高斯分布与特定分布的非线性变换关系。SIRP模型属于外生模型，能够独立控制随机数序列的概率密度函数和协方差矩阵，它的缺点是受所需仿真序列阶数和自相关函数的限制，且当所需的仿真序列较长时计算量很大，不易实现快速算法。。
发明内容
本发明的目的在于提出基于梯度插值的三维相干随机数生成方法，能够快速生成相干随机数，易于实现，计算量小。
本发明基于梯度插值的三维相干随机数生成方法，包括：
S1、获取待匹配点周围的N个插值点，为该插值点生成非相干随机数；基于所述非相干随机数查询预设的映射关系，确定对应插值点的梯度矢量；
S2、获取每个插值点的坐标、以及所述N个插值点所在几何体的中心点坐标；根据插值点坐标、中心点坐标和每个插值点的梯度矢量，确定每个插值点的线性函数；
S3、采用三线性插值法对每个插值点的线性函数进行插值处理，得到待匹配点的相干随机数；
其中，所述映射关系是指非相干随机数与预设的梯度矢量集中每个梯度矢量之间的一一对应关系。
优选地，采用如下方法构建梯度矢量集：
根据地貌起伏的坡度确定梯度取值范围；
在所述梯度取值范围的幅值内等间隔地取A个值，在所述梯度取值范围的方向范围内按预设角度间隔取B个值，得到包含A×B个梯度矢量的梯度矢量集。
优选地，A的取值为10，B的取值为10°。
优选地，采用如下方法确定所述映射关系：
将非相干随机数的取值范围按等间隔分为A×B份；
使每个非相干随机数分别对应梯度矢量集中的一个梯度矢量，得到非相干随机数与梯度矢量集中每个的梯度矢量之间的一一对应关系。
优选地，N的取值为8，获取待匹配点周围的N个插值点包括：
在待匹配点处建立一个棱边长度为单位长度的立方体包围盒；
以所述立方体包围盒的八个顶点作为插值点。
优选地，所述单位长度为随机数步长。
优选地，以(i,j,k)标记为插值点坐标，其中： x、y、z分别为所述立方体包围盒的中心点坐标，符号表示向下取整。
优选地，每个插值点的线性函数为：gi,j,k·(x-i,y-j,z-k)，gi,j,k为插值点的梯度矢量，(i,j,k)为插值点坐标，(x,y,z)为中心点坐标。
优选地，采用如下公式进行三线性插值：s(t)＝6t5-15t4+10t3；其中，t为插值点的线性函数；s(t)为待匹配点的相干随机数。。
本发明基于待匹配点周围的N个插值点进行三线性梯度插值处理，能够降低插值处理维数，能够快速生成待匹配点的相干随机数，易于实现，计算量小。
附图说明
通过以下参照附图而提供的具体实施方式部分，本发明的特征和优点将变得更加容易理解，在附图中：
图1a是示出相干随机序列的示意图；图1b是示出非相干随机序列的示意图
图2是本发明优选实施例中基于梯度插值的三维相干随机数生成方法流程图；
图3是本发明优选实施例中慢起伏地貌模型的示意图；
图4是本发明优选实施例中快起伏地貌模型的示意图；
图5是本发明优选实施例中立方体包围盒的示意图。
具体实施方式
下面参照附图对本发明的示例性实施方式进行详细描述。对示例性实施方式的描述仅仅是出于示范目的，而绝不是对本发明及其应用或用法的限制。
地面环境是一类在局部范围内具有连续性，而在大尺度范围内具有随机性的粗糙表面，这种特性正好与相干随机数的特性一致，因而可用三维相干随机数进行描述。梯度矢量是描述局部连续的有力工具，若局部范围内的梯度矢量变化具有连续性，则在局部范围内的表面起伏具有连续性。大尺度范围的随机性则可通过梯度矢量选择的随机性进行刻画。为了生成空间中某一点P处的相干随机数，可以在P点相邻各个位置点随机选择一组梯度矢量，然后对该矢量集进行三线性插值，即可得到P点处的随机数和梯度矢量，同时也保证了P点局域梯度矢量变化具有连续性。
如图2所示，本发明基于梯度插值的三维相干随机数生成方法包括：
S1、获取待匹配点周围的N个插值点，为该插值点生成非相干随机数；基于非相干随机数查询预设的映射关系，确定对应插值点的梯度矢量；
S2、获取每个插值点的坐标、以及N个插值点所在几何体的中心点坐标；根据插值点坐标、中心点坐标和每个插值点的梯度矢量，确定每个插值点的线性函数；
S3、采用三线性插值法对每个插值点进行插值处理，得到待匹配点的相干随机数；
其中，映射关系是指非相干随机数与预设的梯度矢量集中每个梯度矢量之间的一一对应关系。
不同地貌的起伏不同，表现为梯度变化范围不一样，例如平原地区的梯度幅度范围可取为[-0.1,0.1]，而陡峭的山地的梯度幅度范围可取为[-2,2]，这些梯度的取值是依据不同地貌起伏的坡度确定的。根据地貌的取值范围和分布特性确定梯度取值范围，能够产生更加符合实际地貌的高程起伏数据。
梯度矢量不仅具有幅度，还具有方向。在一些实施例中，可以采用如下方法构建梯度矢量集：根据地貌起伏的坡度确定梯度取值范围；在梯度取值范围的幅值内等间隔地取A个值，在梯度取值范围的方向范围内按预设角度间隔取B个值，得到包含A×B个梯度矢量的梯度矢量集。本领域技术人员可以根据实际情况选择A和B的取值，例如，在梯度取值范围的幅度范围内等间隔的取10个值，在梯度取值范围的方向范围内按10°间隔取36个值，则梯度矢量集G具有360个梯度矢量。A的取值越大或者B的值越大，生成的相干随机数越准确，但是插值处理的计算量也就越大。图2和图3分别给出了利用本专利方法产生的慢起伏地貌模型和快起伏地貌模型。其中，慢起伏模型的梯度取值范围为[-0.2,0.2]，梯度取值服从高斯分布；快起伏模型的梯度取值范围为[-2,2]，梯度取值服从均匀分布。
在一些实施例中，采用如下方法确定所述映射关系：将非相干随机数的取值范围按等间隔分为A×B份；使每个非相干随机数分别对应梯度矢量集中的一个梯度矢量，得到非相干随机数与梯度矢量集中每个的梯度矢量之间的一一对应关系。例如，当梯度矢量集G具有360个梯度矢量时，可以把非相干随机数R的取值范围按等间隔分为360份，分别与梯度矢量集G中的每个梯度矢量一一对应起来。
为了生成空间中待匹配点P(xp,yp,zp)处的随机数，需要首先获取待匹配点周围的N个插值点。在一些实施例中，在待匹配点处建立一个棱边长度为单位长度的立方体包围盒；以立方体包围盒的八个顶点作为插值点。图5是本发明优选实施例中立方体包围盒的示意图，其中，p1、p2、p3、p4、p5、p6、p7、p8分别代表立方体包围盒的八个顶点。实际应用过程中，本领域技术人员也可以根据实际情况选择其他的方式获取待匹配点周围的插值点，若采用本发明上述实施例的获取思路，也可以根据实际情况通过建立椭球体或者其他多边形体来获取待匹配点周围的插值点，本发明对此不做具体限定。
在上述立方体包围盒中，单位长度可以是随机数步长。单位长度的取值越小，插值处理的计算量越大，插值结果越精确。一般来说，空间模型的分辨率越高，单位长度越小。
在一些实施例中，以(i,j,k)标记插值点的坐标，其中： x、y、z分别为所述立方体包围盒的中心点坐标，符号表示向下取整。取整处理能够简化插值处理的计算量，提高计算速度。
每个插值点的线性函数可以表示为：gi,j,k·(x-i,y-j,z-k)，gi,j,k为插值点的梯度矢量，(i,j,k)为插值点坐标，(x,y,z)为中心点坐标。本领域技术人员可以根据实际情况选择合适的三线性插值方法，在本发明的优选实施例中，采用如下公式进行三线性插值：s(t)＝6t5-15t4+10t3；其中，t为插值点的线性函数；s(t)为待匹配点的相干随机数。
在图2示出的优选实施例中，获得当前待匹配点的相干随机数之后，判断是否遍历所有位置。若否，按照本发明的方法生成下一位置的相干随机数；若是，则结束三维相干随机数的生成步骤。至此，一个具有三维空间相干分布的点云数据集已经生成完毕，若要进一步生成如图2和图3的地貌起伏网格模型，可以进一步根据地貌起伏模型对点云数据进行Delaunay三角剖分处理，获取下一位置作为待匹配点。
虽然参照示例性实施方式对本发明进行了描述，但是应当理解，本发明并不局限于文中详细描述和示出的具体实施方式，在不偏离权利要求书所限定的范围的情况下，本领域技术人员可以对所述示例性实施方式做出各种改变。

Claims (9)

1.基于梯度插值的三维相干随机数生成方法，其特征在于包括：
S1、获取待匹配点周围的N个插值点，为该插值点生成非相干随机数；基于所述非相干随机数查询预设的映射关系，确定对应插值点的梯度矢量；
S2、获取每个插值点的坐标、以及所述N个插值点所在几何体的中心点坐标；根据插值点坐标、中心点坐标和每个插值点的梯度矢量，确定每个插值点的线性函数；
S3、采用三线性插值法对每个插值点的线性函数进行插值处理，得到待匹配点的相干随机数；
其中，所述映射关系是指非相干随机数与预设的梯度矢量集中每个梯度矢量之间的一一对应关系。
2.如权利要求1所述的方法，其特征在于，采用如下方法构建梯度矢量集：
根据地貌起伏的坡度确定梯度取值范围；
在所述梯度取值范围的幅值内等间隔地取A个值，在所述梯度取值范围的方向范围内按预设角度间隔取B个值，得到包含A×B个梯度矢量的梯度矢量集。
3.如权利要求2所述的方法，其特征在于，A的取值为10，B的取值为10°。
4.如权利要求2所述的方法，其特征在于，采用如下方法确定所述映射关系：
将非相干随机数的取值范围按等间隔分为A×B份；
使每个非相干随机数分别对应梯度矢量集中的一个梯度矢量，得到非相干随机数与梯度矢量集中每个的梯度矢量之间的一一对应关系。
5.如权利要求1所述的方法，其特征在于，N的取值为8，获取待匹配点周围的N个插值点包括：
在待匹配点处建立一个棱边长度为单位长度的立方体包围盒；
以所述立方体包围盒的八个顶点作为插值点。
6.如权利要求5所述的方法，其特征在于，所述单位长度为随机数步长。
7.如权利要求5所述的方法，其特征在于，以(i,j,k)标记为插值点坐标，其中：x、y、z分别为所述立方体包围盒的中心点坐标，符号表示向下取整。
8.如权利要求7所述的方法，其特征在于，每个插值点的线性函数为：gi,j,k·(x-i,y-j,z-k)，gi,j,k为插值点的梯度矢量，(i,j,k)为插值点坐标，(x,y,z)为中心点坐标。
9.如权利要求8所述的方法，其特征在于，采用如下公式进行三线性插值：s(t)＝6t5-15t4+10t3；其中，t为插值点的线性函数；s(t)为待匹配点的相干随机数。
CN201710236412.9A 2017-04-12 2017-04-12 基于梯度插值的三维相干随机数生成方法 Active CN107045435B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710236412.9A CN107045435B (zh) 2017-04-12 2017-04-12 基于梯度插值的三维相干随机数生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710236412.9A CN107045435B (zh) 2017-04-12 2017-04-12 基于梯度插值的三维相干随机数生成方法

Publications (2)

Publication Number Publication Date
CN107045435A true CN107045435A (zh) 2017-08-15
CN107045435B CN107045435B (zh) 2019-08-23

Family

ID=59545664

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710236412.9A Active CN107045435B (zh) 2017-04-12 2017-04-12 基于梯度插值的三维相干随机数生成方法

Country Status (1)

Country Link
CN (1) CN107045435B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5909520A (en) * 1997-08-25 1999-06-01 The United States Of America As Represented By The Secretary Of The Navy Noise coding processor
CN101881830A (zh) * 2010-03-15 2010-11-10 中国电子科技集团公司第十研究所 重建雷达扫描数据生成三维可视地形的方法
CN102741854A (zh) * 2009-10-23 2012-10-17 埃克森美孚上游研究公司 利用梯度信息进行优化的方法
CN103034471A (zh) * 2012-12-10 2013-04-10 杜海洋 一种随机数的生成方法及系统
CN105653785A (zh) * 2015-08-07 2016-06-08 西安电子科技大学 一种雷达杂波模拟方法、装置及系统
US20160328211A1 (en) * 2011-09-30 2016-11-10 Los Alamos National Security, Llc Quantum random number generators

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5909520A (en) * 1997-08-25 1999-06-01 The United States Of America As Represented By The Secretary Of The Navy Noise coding processor
CN102741854A (zh) * 2009-10-23 2012-10-17 埃克森美孚上游研究公司 利用梯度信息进行优化的方法
CN101881830A (zh) * 2010-03-15 2010-11-10 中国电子科技集团公司第十研究所 重建雷达扫描数据生成三维可视地形的方法
US20160328211A1 (en) * 2011-09-30 2016-11-10 Los Alamos National Security, Llc Quantum random number generators
CN103034471A (zh) * 2012-12-10 2013-04-10 杜海洋 一种随机数的生成方法及系统
CN105653785A (zh) * 2015-08-07 2016-06-08 西安电子科技大学 一种雷达杂波模拟方法、装置及系统

Also Published As

Publication number Publication date
CN107045435B (zh) 2019-08-23

Similar Documents

Publication Publication Date Title
US8014568B2 (en) Method for computer-aided identification of the child octants of a parent octant, which are intersected by a beam, in an octree data structure by means of look-up tables
CN109633570B (zh) 一种快速任意形状随机箔条云团建模方法
CN112819962B (zh) 数字图像相关中非均匀网格划分及局部网格疏密方法
Leeper et al. Point clouds can be represented as implicit surfaces for constraint-based haptic rendering
CN104182571B (zh) 基于Delaunay和GPU的Kriging插值方法
Sane et al. An Interpolation Scheme for VDVP Lagrangian Basis Flows.
CN109993801A (zh) 一种用于二维相机与三维传感器的标定装置及标定方法
CN111402422B (zh) 三维表面重建方法、装置和电子设备
CN107481319B (zh) 一种隐式曲面随机点云生成器
CN105469415A (zh) 多视角遥感图像融合方法
CN110363863B (zh) 一种神经网络的输入数据生成方法和系统
CN107045435A (zh) 基于梯度插值的三维相干随机数生成方法
Garcia et al. CPU-based real-time surface and solid voxelization for incomplete point cloud
CN113126944B (zh) 深度图的显示方法、显示装置、电子设备及存储介质
CN109740234B (zh) 一种家居布局3d碰撞检测方法及系统
CN103020474A (zh) 基于杂波统计分布模型和功率谱模型确定杂波模型的方法
CN113971694A (zh) 一种点云数据处理方法及装置
Nengem et al. A Study of the Effects of the Shape Parameter and Type of Data Points Locations on the Accuracy of the Hermite-based Symmetric Approach Using Positive Definite Radial Kernels
CA2877507A1 (en) Truncation diagram determination for a pluri-gaussian estimation
CN111965694B (zh) 地震物理点位置确定方法、装置及地震观测系统
CN114862994A (zh) 基于柏林噪声的静态云雾仿真生成方法
Chen et al. An efficient numerical method for determination of shapes, sizes and orientations of flaws for nondestructive evaluation
Ye et al. Study on self-intersection detection algorithm based on OBJ model
JP2692587B2 (ja) 画像処理方法
Zhang et al. Bathymetric Data Fitting based on Fourier Series of Linear Complementarity

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