CN104572580A - 一个表达矿物颗粒光散射空间分布特征的函数的构建方法 - Google Patents
一个表达矿物颗粒光散射空间分布特征的函数的构建方法 Download PDFInfo
- Publication number
- CN104572580A CN104572580A CN201410799627.8A CN201410799627A CN104572580A CN 104572580 A CN104572580 A CN 104572580A CN 201410799627 A CN201410799627 A CN 201410799627A CN 104572580 A CN104572580 A CN 104572580A
- Authority
- CN
- China
- Prior art keywords
- phase function
- function
- expression
- integral
- light scattering
- 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
Links
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Image Generation (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一个表达矿物颗粒光散射空间分布特征的函数的构建方法,它有4个步骤:步骤一、相函数的选定;步骤二、相函数积分表达式的推导;步骤三、相函数积分表达式的分析;步骤四、新相函数的构建。与目前采用的相函数最大的不同是,它构建的相函数满足了散射相函数必须满足的一个基本物理条件,即对散射相函数在0~2π范围内积分时,积分值恒定为2π。本发明属于遥感地质领域,它适用于岩石与矿物反射与发射光谱的数值模拟与仿真,由于满足了颗粒光散射基本物理条件,有助于提高岩矿光谱的数值模拟精度。
Description
技术领域
本发明涉及一个表达矿物颗粒光散射空间分布特征的函数的构建方法,构建的函数是一个定量表达散射能量空间分布特征的函数,属于遥感地质领域,它适用于岩石与矿物反射与发射光谱的数值模拟与仿真。
背景技术
定量表达光线被矿物颗粒散射后散射能量空间分布特征的函数,是岩石与矿物反射与发射光谱的数值模拟与仿真模型的重要的数学基础。目前主要采用数学多项式来表达,应用最广泛的是勒让德多项式。该表达式也称为相函数,是相位角的函数。相位角是入射光线与散射光线的夹角,变化范围为0~2π,即0-180度。
严格来讲,散射相函数必须满足的一个光散射特性的基本物理条件是,对散射相函数在0~2π范围内积分时,积分值恒定为2π。而事实上,有些相函数无法满足这一基本条件,直接用于岩矿光谱的数值模拟会影响模拟精度。因此,开发满足这一基本条件的相函数的构建方法有助于提高岩矿光谱数值模拟精度。
发明内容
1、目的:本发明的目的是提供一个表达矿物颗粒光散射空间分布特征的函数的构建方法,构建的相函数满足相函数的基本物理条件,即对散射相函数在0~2π范围内积分时积分值恒定为2π。
2、需解决的技术问题
通过对采用的相函数在0~2π范围内的积分表达式进行推导,分析积分表达式是否恒定为2π。如不为2π,且相对于2π的变化量与相函数变化无关,则可通过在原有相函数的表达式中减去一个变化量,来使新的相函数满足光散射特性的基本物理条件。
3、技术方案
本发明针对上述需要解决的技术问题,提出了相应的解决方案。整体解决方案见附图1。本发明是一个表达矿物颗粒光散射空间分布特征的函数的构建方法,该方法具体步骤如下:
步骤一:相函数的选定
选定采用什么数学表达式来作为相函数表达式。常用的相函数表达式有1阶勒让德多项式,2阶勒让德多项式。
步骤二:相函数的积分表达式的推导
推导相位角在0~2π范围内函数的积分表达式。积分表达式推导的公式为,
其中,A为相函数在0~2π范围内函数的积分表达式,P为相函数,g为相位角。
步骤三:相函数积分表达式的分析
分析相函数积分表达式是否恒定为2π,如果是,则表明相函数满足光散射特性的基本物理条件,不需构建新的相函数。否则需构建,应分析积分值相对于2π的变化量是否与相函数变化无关,如无关,则可构建新的相函数。
步骤四:新相函数的构建
上一步计算的相函数积分值减去2π,得到一个变化量。该变化量除以2π,得到一个相函数的修正系数。在步骤一选定的相函数表达式基础上减去这一修正系数,得到新的相函数。
4、优点与功效
本发明一个表达矿物颗粒光散射空间分布特征的新函数的构建方法,其优点是:构建的相函数满足满足相函数的基本物理条件,即对散射相函数在0~2π范围内积分时积分值恒定为2π。
附图说明
图1为本发明一个表达矿物颗粒光散射空间分布特征的函数的构建方法的实现流程示意图
图2为二阶勒让德多项式表达的相函数随相位角的关系图
图3为二阶勒让德多项式表达的相函数与新函数的对比图(b=1,c=1)
图4为二阶勒让德多项式表达的相函数与新函数的对比图(b=-1,c=1)
具体实施方式
见图1,为了更好的说明本发明一个表达矿物颗粒光散射空间分布特征的函数的构建方法,构建了一个新的相函数。具体步骤如下:
步骤一:相函数的选定
选定2阶勒然德多项式作为相函数表达式,其表达式为,
其中P为散射相函数,g为相位角,是入射光线与散射光线的夹角,变化范围为0~2π,即0-180度,b、c为控制散射光空间分布特征的系数,b、c的变化范围为分别-1.0~1.0、0~1.0,二者全部为0时,散射为各向同性,即散射相函数不随相位角变化,恒定为1。b为负值时为前向散射,为正值时为后向散射,随着c增加,侧向散射(即相位角为90度时)能量逐渐减小,见图2。
步骤二:相函数的积分表达式的推导
步骤三:相函数积分表达式的分析
分析上一步骤推导的表达式可知,相函数在0~2π范围内的积分表达式不恒定为2π,且变化量cp/2与相位角无关。
步骤四:新相函数的构建
上一步计算的变化量cp/2除以2p,得到相函数修正系数c/4,在步骤一选定的2阶勒然德多项式的基础上减去修正系数c/4,得到新的相函数,其表达式为,
其中P为散射相函数,g为相位角,是入射光线与散射光线的夹角,变化范围为0~2π,即0-180度,b、c为控制散射光空间分布特征的系数。b、c的变化范围为分别-1.0~1.0、0~1.0,二者全部为0时,散射为各向同性,即散射相函数不随相位角变化,恒定为1。b为负值时为前向散射,为正值时为后向散射,随着c增加,侧向散射(即相位角为90度时)能量逐渐减小。图3、图4仅表示b为1或-1,c为1时的情况。
新的相函数,在0~2π范围内的积分表达式恒定为2π,推导如下:
Claims (1)
1.一个表达矿物颗粒光散射空间分布特征的函数的构建方法,其特征在于:该方法具体步骤如下:
步骤一:相函数的选定
选定采用什么数学表达式来作为相函数表达式,常用的相函数表达式有1阶勒让德多项式,2阶勒让德多项式;
步骤二:相函数的积分表达式的推导
推导相位角在0~2π范围内函数的积分表达式,积分表达式推导的公式为,
其中,A为相函数在0~2π范围内函数的积分表达式,P为相函数,g为相位角;
步骤三:相函数积分表达式的分析
分析相函数积分表达式是否恒定为2π,如果是,则表明相函数满足光散射特性的基本物理条件,不需构建新的相函数;否则需构建,应分析积分值相对于2π的变化量是否与相函数变化无关,如无关,则构建新的相函数;
步骤四:新相函数的构建
上一步计算的相函数积分值减去2π,得到一个变化量,该变化量除以2π,得到一个相函数的修正系数,在步骤一选定的相函数表达式基础上减去这一修正系数,得到新的相函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410799627.8A CN104572580B (zh) | 2014-12-19 | 2014-12-19 | 一个表达矿物颗粒光散射空间分布特征的函数的构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410799627.8A CN104572580B (zh) | 2014-12-19 | 2014-12-19 | 一个表达矿物颗粒光散射空间分布特征的函数的构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104572580A true CN104572580A (zh) | 2015-04-29 |
CN104572580B CN104572580B (zh) | 2018-06-19 |
Family
ID=53088686
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410799627.8A Expired - Fee Related CN104572580B (zh) | 2014-12-19 | 2014-12-19 | 一个表达矿物颗粒光散射空间分布特征的函数的构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104572580B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1595110A (zh) * | 2004-06-29 | 2005-03-16 | 中国国土资源航空物探遥感中心 | 高光谱多端元线性拟合矿物识别与矿化指示方法 |
CN1595109A (zh) * | 2004-06-29 | 2005-03-16 | 中国国土资源航空物探遥感中心 | 基于光谱混合组成的高光谱矿物极大相关识别方法 |
CN1595203A (zh) * | 2004-06-29 | 2005-03-16 | 中国国土资源航空物探遥感中心 | 高光谱矿物分层谱系识别方法 |
CN101551471A (zh) * | 2009-05-19 | 2009-10-07 | 中国国土资源航空物探遥感中心 | 一种高光谱遥感矿物含量定量反演方法 |
CN101782516A (zh) * | 2010-03-18 | 2010-07-21 | 华中科技大学 | 一种计算密实性颗粒介质光散射特性的方法 |
-
2014
- 2014-12-19 CN CN201410799627.8A patent/CN104572580B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1595110A (zh) * | 2004-06-29 | 2005-03-16 | 中国国土资源航空物探遥感中心 | 高光谱多端元线性拟合矿物识别与矿化指示方法 |
CN1595109A (zh) * | 2004-06-29 | 2005-03-16 | 中国国土资源航空物探遥感中心 | 基于光谱混合组成的高光谱矿物极大相关识别方法 |
CN1595203A (zh) * | 2004-06-29 | 2005-03-16 | 中国国土资源航空物探遥感中心 | 高光谱矿物分层谱系识别方法 |
CN101551471A (zh) * | 2009-05-19 | 2009-10-07 | 中国国土资源航空物探遥感中心 | 一种高光谱遥感矿物含量定量反演方法 |
CN101782516A (zh) * | 2010-03-18 | 2010-07-21 | 华中科技大学 | 一种计算密实性颗粒介质光散射特性的方法 |
Non-Patent Citations (2)
Title |
---|
NADIA T. ZAKHAROVA等: "Scattering properties of needlelike and platelike ice spheroids with moderate size parameters", 《OSA PUBLISHING》 * |
秦继民等编著: "《普通物理学专题研究(上)》", 31 July 1993, 广西师范大学出版社 * |
Also Published As
Publication number | Publication date |
---|---|
CN104572580B (zh) | 2018-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sato et al. | First Monte Carlo analysis of fragmentation functions from single-inclusive e+ e-annihilation | |
WO2015030990A3 (en) | Three-dimensional reservoir pressure determination using real time pressure data from downhole gauges | |
WO2014097124A3 (en) | Quantitative imaging | |
WO2016063235A3 (en) | Visualization of imaging uncertainty | |
CN102346808A (zh) | 一种从hj-1卫星数据反演lai的方法 | |
CN104237867B (zh) | 引信接收机的有效多普勒带宽的确定方法和装置 | |
CN104572580A (zh) | 一个表达矿物颗粒光散射空间分布特征的函数的构建方法 | |
CN105093265A (zh) | 一种模拟地震波在ti介质中传播规律的方法 | |
Omerbashich | Earth-model discrimination method | |
Santiago Pullarello et al. | Two Matlab routines for rock material propagation modelling: PLANOS, a volume assessment by infinite planar discontinuities, and PrePostDAN3D, a data treatment tool for before and after DAN3D processing | |
Somogyvári et al. | Assessing the convergence of transdimensional inversion used for subsurface geometry reconstruction | |
Krivonos et al. | Comparative analysis of methods of evaluating the lower ionosphere parameters by tweek atmospherics | |
Alcorn et al. | The 1D and 2D Hα Kinematics of Galaxies in ZFIRE at z~ 2 | |
Carrano et al. | New Opportunities for Remote Sensing Ionospheric Irregularities by Fitting Scintillation Spectra | |
Blake et al. | Dynamics of basaltic plumbing systems-a theoretical model of eruptive output and timescales | |
Aldaya et al. | Determination of the top quark mass from the cross section measured by CMS at√ s= 7 TeV | |
Minguez et al. | Magnetic Susceptibility Cyclostratigraphy of the Ediacaran Wonoka Formation, South Australia, from a Portable Susceptibility Meter | |
Hassanzadeh et al. | Reduced-Order Modeling of 3D Rayleigh-Benard Turbulent Convection | |
Hoffmann | Beyond Cloud Microphysics: Representing Subgrid-Scale Processes in Lagrangian Cloud Models | |
Ullmann de Brito et al. | Well control optimization using a two-step surrogate treatment | |
Shemer et al. | Measurements of wind-waves under transient wind conditions. | |
Sheibani et al. | Optimizing Shear Stresses at the Tip of a Hydraulic Fracture-What Is the Ideal Orientation of Natural Fractures with respect to Hydraulic Fracture? | |
Garth et al. | Global and local waveform simulations using the VERCE platform | |
Masson et al. | Accounting for propagation outside of the model boundaries in regional full waveform inversion based on adjoint methods | |
Bock et al. | A statistical approach for validating eSOTER and digital soil maps in front of traditional soil maps |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180619 Termination date: 20191219 |