CN110927190A - 在劳厄衍射图谱中识别同一晶带轴衍射峰的方法 - Google Patents
在劳厄衍射图谱中识别同一晶带轴衍射峰的方法 Download PDFInfo
- Publication number
- CN110927190A CN110927190A CN201911313367.8A CN201911313367A CN110927190A CN 110927190 A CN110927190 A CN 110927190A CN 201911313367 A CN201911313367 A CN 201911313367A CN 110927190 A CN110927190 A CN 110927190A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- diffraction
- angle
- circle
- circles
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
- G01N23/2055—Analysing diffraction patterns
Landscapes
- Chemical & Material Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种在劳厄衍射图谱中识别同一晶带轴衍射峰的方法,方法包括以下步骤:建立探测器坐标系,计算在该探测器坐标系下各衍射峰的衍射矢量,以探测器坐标系的原点为球心建立单位球壳S;求取垂直于该衍射峰的衍射矢量且过该探测器坐标系的原点的各平面Pi,求取平面Pi与单位球壳S相交的圆Ci,定义位角θ与方位角步长Δθ,俯仰角与俯仰角步长依据方位角步长与俯仰角步长,将球壳的上半球壳划分为m个单位区域Dj,在单位区域Dj中求各圆Ci之间的广义交点,每个广义交点对应着一个晶带轴,该晶带轴的方向为以该广义交点为终点,以探测器坐标系的原点为起点的单位矢量,相交于该广义交点的圆Ci对应的衍射峰属于该广义交点对应的晶带轴。
Description
技术领域
本发明属于X射线衍射技术领域,特别是一种在劳厄衍射图谱中识别同一晶带轴衍射峰的方法。
背景技术
作为一种高效的材料表征手段,X射线劳厄衍射被广泛用于多种材料的表征。基于劳厄衍射图谱,能够通过分析计算得到材料的取向、组成相和应变等各种信息。由于劳厄衍射原理简单,在进行衍射实验时无需转动样品、X射线光源或探测器的空间位置,从而能够实现劳厄图谱的快速采集,相比其他相似的表征手段,在速度上具有极大的优势。而在对劳厄衍射图谱进行标定后,通过对衍射峰位置的回归拟合,可以求解包括三维晶体取向,二阶应力应变张量等信息。使用诸如传统的X射线衍射技术、电子衍射技术等方法,都无法实现同时对这些材料信息的高精度表征。进而在材料表征领域,劳厄衍射有其不可替代的优势。
而在分析劳厄衍射图谱前,需要对劳厄衍射图谱进行标定。对于如矿物样品等有着较大晶胞的材料,使用劳厄衍射往往能采集到数十甚至数百个衍射峰,而对于这些衍射峰的标定,往往需要花费大量的时间,且由于基于劳厄衍射图谱的所有其他分析都需要使用到衍射峰的标定结果,故分析这样的劳厄衍射图谱必须花费的时间一般很长。对于一般的个人电脑,以现有的算法和商用化的软件,分析这样的衍射峰往往需要数分钟甚至十数分钟。而对于衍射峰较多的劳厄衍射图谱,属于各晶带轴的衍射峰都较多,通过对属于同一晶带轴的衍射峰的识别并获得对应晶带轴的方向可以极大的简化标定过程,减少所需要的时间,对于分析取向等简单的材料信息,甚至无需继续对衍射峰进行标定,仅需要对晶带轴进行标定。
而现有的劳厄衍射图谱分析算法往往无法实现对属于同一晶带轴的衍射峰进行识别。已有的方法中,往往都使用椭圆或双曲线拟合的方法直接对衍射峰在劳厄衍射图谱的位置进行拟合。这样对衍射峰的位置的拟合过程往往耗时较多,即无法达到简化标定分析过程,减少标定分析耗时的目的。
在背景技术部分中公开的上述信息仅仅用于增强对本发明背景的理解,因此可能包含不构成在本国中本领域普通技术人员公知的现有技术的信息。
发明内容
针对现有技术中存在的问题,本发明提出一种在劳厄衍射图谱中识别同一晶带轴衍射峰的方法。
本发明的目的是通过以下技术方案予以实现,一种在劳厄衍射图谱中识别同一晶带轴衍射峰的方法包括以下步骤:
第二步骤中,以探测器坐标系的原点为球心建立单位球壳S;
第四步骤中,求取平面Pi与单位球壳S相交的圆Ci,其中i=1,2,3…n;
第六步骤中,在单位区域Dj中求各圆Ci之间的广义交点,每个广义交点对应着一个晶带轴,该晶带轴的方向为以该广义交点为终点,以探测器坐标系的原点为起点的单位矢量,相交于该广义交点的圆Ci对应的衍射峰属于该广义交点对应的晶带轴。
所述的方法中,第六步骤中,在单位区域Dj中求各圆Ci之间的广义交点时,
第一步,进入单位区域Dj的圆Ci组成圆的集合Ci′,其中i′=1,2,3…n′,n′表示该集合中圆的总数,单位区域Dj的左侧边界为方位角θ=θ1,右侧边界为方位角θ=θ2,若n′<np,则在该单位区域Dj中不存在广义交点,np为预定参数;
第四步,将序列L1和序列L2中位于相同位置的数字相减并求绝对值以组成序列L3;若序列L3中不存在非零值,则认为该单位区域Dj的中心为广义交点,若序列L3中存在非零最大值,求其所在位置,如序列L1该位置的数字为t,则将编号为t的圆从集合Ci′删去;
第五步,若删去后集合Ci′中圆的总数大于np,则跳至第二步(S6-2);若删去后集合Ci′中圆的总数少于np,则在该单位区域Dj中不存在广义交点。
所述的方法中,参数np为形成一个广义交点需要的最少的相交的圆的数量,
所述的方法中,参数np为一般为5、6或7。
所述的方法中,探测器坐标系为三维直角坐标系。
相比现有技术,本发明中所介绍的方法计算简单,对劳厄衍射峰的质量要求较低,同时可以适用于衍射峰数目从十数个到数百个的劳厄衍射图谱。本方法具有简便、耗时少等特点。
附图说明
通过阅读下文优选的具体实施方式中的详细描述,本发明各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。说明书附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。显而易见地,下面描述的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。而且在整个附图中,用相同的附图标记表示相同的部件。
在附图中:
图1是根据本发明一个实施例的在劳厄衍射图谱中识别同一晶带轴衍射峰的方法的步骤示意图;
图2是根据本发明一个实施例的劳厄衍射图谱示意图。
以下结合附图和实施例对本发明作进一步的解释。
具体实施方式
下面将参照附图1至图2更详细地描述本发明的具体实施例。虽然附图中显示了本发明的具体实施例,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
需要说明的是,在说明书及权利要求当中使用了某些词汇来指称特定组件。本领域技术人员应可以理解,技术人员可能会用不同名词来称呼同一个组件。本说明书及权利要求并不以名词的差异来作为区分组件的方式,而是以组件在功能上的差异来作为区分的准则。如在通篇说明书及权利要求当中所提及的“包含”或“包括”为一开放式用语,故应解释成“包含但不限定于”。说明书后续描述为实施本发明的较佳实施方式,然该描述乃以说明书的一般原则为目的,并非用以限定本发明的范围。本发明的保护范围当视所附权利要求所界定者为准。
为便于对本发明实施例的理解,下面将结合附图以具体实施例为例做进一步的解释说明,且各个附图并不构成对本发明实施例的限定。
为了更好地理解,如图1所示,一种在劳厄衍射图谱中识别同一晶带轴衍射峰的方法包括以下步骤:
第二步骤S2中,以探测器坐标系的原点为球心建立单位球壳S;
第四步骤S4中,求取平面Pi与单位球壳S相交的圆Ci,其中i=1,2,3…n;
第六步骤S6中,在单位区域Dj中求各圆Ci之间的广义交点,每个广义交点对应着一个晶带轴,该晶带轴的方向为以该广义交点为终点,以探测器坐标系的原点为起点的单位矢量,相交于该广义交点的圆Ci对应的衍射峰属于该广义交点对应的晶带轴。
为使对本发明的叙述更为清晰明了,现使用单晶硅的劳厄衍射图谱为例具体描述本发明的具体实施步骤与技术细节。
本具体实施方法中使用的劳厄衍射图谱如图2所示。具体实施方法为:
步骤1:使用公知的方法在图2所示劳厄衍射图谱中寻衍射峰,共寻到衍射峰39个,其在探测器上的位置为(xi yi),其中i=1,2,3…39。建立探测器坐标系,使用公知的方法计算该坐标系下各衍射峰的衍射矢量并以单位矢量表示其中i=1,2,3…39。这里衍射矢量的具体计算方法为:读取探测器的空间转角、和。有
得到矩阵A后,读取该点样品上照射点到探测器平面的距离d,利用公式可得各衍射峰出射束kout,i在探测器坐标系下为
步骤2:以探测器坐标系的原点为球心建立单位球壳S。
步骤4:求39个平面Pi与单位球壳S相交的39各在球壳S上的圆Ci。
步骤5:以方位角步长Δθ=5°,俯仰角步长将俯仰角的上半球壳划分为1296个单位区域Dj,j=1,2,3…1296。这里方位角θ为球壳上的点在探测器坐标系x-y平面的投影与探测器坐标系原点连线和探测器坐标系x轴的夹角,俯仰角为球壳上的点与探测器原点连线和探测器坐标系z轴的夹角。
步骤6:定义参数np=5。在每个单位区域Dj中求各圆Ci之间的广义交点。每个广义交点对应着一个晶带轴,该晶带轴的方向为以该广义交点为终点,以探测器坐标系的原点为起点的单位矢量。相交于该广义交点的圆Ci对应的衍射峰属于该广义交点对应的晶带轴。求单位区域中的广义交点的方法为:
第一步,定义所有进入单位区域Dj的圆Ci组成圆的集合Ci′,其中i′=1,2,3…n′,n′表示该集合中圆的总数。单位区域Dj的左侧边界为方位角θ=θ1,右侧边界为方位角θ=θ2。若n′<np,则在该单位区域Dj中不存在广义交点,且不需要执行以下各步骤,这里np=5。
第二步,对所有集合Ci′中的所有圆,求其当方位角为θ1时的俯仰角其中i1=1,2,3…n′。依据俯仰角的值的降序为集合Ci′中的所有圆编号为1,2,3…n′,并将这些编号按此顺序组成序列L1,n′表示此时集合中圆的总数。
第三步,对所有集合Ci′中的所有圆,求其当方位角为θ2时俯仰角其中i2=1,2,3…n′,依据俯仰角的值的升序将为集合Ci′中的所有圆排序,并以各圆的编号按此顺序组成序列L2,n′表示此时集合Ci′中圆的总数。
第四步,将序列L1和序列L2中位于相同位置的数字相减并求绝对值以组成序列L3。若序列L3中不存在非零值,则认为该单位区域Dj的中心为广义交点,在该广义交点相交的圆为集合Ci′中的所有圆,且不需要执行后续步骤。若序列L3中存在非零最大值,求其所在位置。假设序列L1该位置的数字为t,则将编号为t的圆从集合Ci′删去。
第五步,若删去后集合Ci′中圆的总数大于np,则跳至第二步。若删去后集合Ci′中圆的总数少于np,则在该单位区域Dj中不存在广义交点,且不需要执行后续步骤。这里np=5。
尽管以上结合附图对本发明的实施方案进行了描述,但本发明并不局限于上述的具体实施方案和应用领域,上述的具体实施方案仅仅是示意性的、指导性的,而不是限制性的。本领域的普通技术人员在本说明书的启示下和在不脱离本发明权利要求所保护的范围的情况下,还可以做出很多种的形式,这些均属于本发明保护之列。
Claims (8)
1.一种在劳厄衍射图谱中识别同一晶带轴衍射峰的方法,该方法包括以下步骤:
第二步骤(S2)中,以探测器坐标系的原点为球心建立单位球壳S;
第四步骤(S4)中,求取平面Pi与单位球壳S相交的圆Ci,其中i=1,2,3…n;
第六步骤(S6)中,在单位区域Dj中求各圆Ci之间的广义交点,每个广义交点对应着一个晶带轴,该晶带轴的方向为以该广义交点为终点,以探测器坐标系的原点为起点的单位矢量,相交于该广义交点的圆Ci对应的衍射峰属于该广义交点对应的晶带轴,即识别出该广义交点的圆Ci对应的衍射峰属于同一晶带轴。
2.根据权利要求1所述的方法,其中,优选的,第六步骤(S6)中,在单位区域Dj中求各圆Ci之间的广义交点时,
第一步(S6-1),进入单位区域Dj的圆Ci组成圆的集合Ci′,其中i′=1,2,3…n′,n′表示该集合中圆的总数,单位区域Dj的左侧边界为方位角θ=θ1,右侧边界为方位角θ=θ2,若n′<np,则在该单位区域Dj中不存在广义交点,np为预定参数;
第四步(S6-4),将序列L1和序列L2中位于相同位置的数字相减并求绝对值以组成序列L3;若序列L3中不存在非零值,则认为该单位区域Dj的中心为广义交点,若序列L3中存在非零最大值,求其所在位置,如序列L1该位置的数字为t,则将编号为t的圆从集合Ci′删去;
第五步(S6-5),若删去后集合Ci′中圆的总数大于np,则跳至第二步(S6-2);若删去后集合Ci′中圆的总数少于np,则在该单位区域Dj中不存在广义交点。
4.根据权利要求2所述的方法,其中,参数np为形成一个广义交点需要的最少的相交的圆的数量。
5.根据权利要求2所述的方法,其中,参数np为5、6或7。
7.根据权利要求1所述的方法,其中,探测器坐标系为三维直角坐标系。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911313367.8A CN110927190B (zh) | 2019-12-18 | 2019-12-18 | 在劳厄衍射图谱中识别同一晶带轴衍射峰的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911313367.8A CN110927190B (zh) | 2019-12-18 | 2019-12-18 | 在劳厄衍射图谱中识别同一晶带轴衍射峰的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110927190A true CN110927190A (zh) | 2020-03-27 |
CN110927190B CN110927190B (zh) | 2020-11-06 |
Family
ID=69864400
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911313367.8A Active CN110927190B (zh) | 2019-12-18 | 2019-12-18 | 在劳厄衍射图谱中识别同一晶带轴衍射峰的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110927190B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113325015A (zh) * | 2021-05-12 | 2021-08-31 | 西安交通大学 | 劳厄衍射图谱中识别同一衍射峰所劈裂衍射亚峰的方法 |
CN113325016A (zh) * | 2021-05-12 | 2021-08-31 | 西安交通大学 | 存在衍射峰劈裂的劳厄衍射图谱中标定多个亚晶的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2569862B2 (ja) * | 1990-02-13 | 1997-01-08 | 三菱電機株式会社 | X線露光装置およびx線露光方法 |
CN101710085A (zh) * | 2009-12-23 | 2010-05-19 | 重庆大学 | 基于Matlab的电子衍射指数标定方法 |
CN106950237A (zh) * | 2017-03-10 | 2017-07-14 | 西安交通大学 | 一种基于峰间夹角比对的扫描式劳厄衍射图谱分析方法 |
CN106950234A (zh) * | 2017-03-10 | 2017-07-14 | 西安交通大学 | 一种劳厄衍射扫描数据的衍射峰空间分布比对分析法 |
CN106959312A (zh) * | 2017-03-10 | 2017-07-18 | 西安交通大学 | 一种基于峰位比对的扫描式劳厄衍射图谱分析方法 |
-
2019
- 2019-12-18 CN CN201911313367.8A patent/CN110927190B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2569862B2 (ja) * | 1990-02-13 | 1997-01-08 | 三菱電機株式会社 | X線露光装置およびx線露光方法 |
CN101710085A (zh) * | 2009-12-23 | 2010-05-19 | 重庆大学 | 基于Matlab的电子衍射指数标定方法 |
CN106950237A (zh) * | 2017-03-10 | 2017-07-14 | 西安交通大学 | 一种基于峰间夹角比对的扫描式劳厄衍射图谱分析方法 |
CN106950234A (zh) * | 2017-03-10 | 2017-07-14 | 西安交通大学 | 一种劳厄衍射扫描数据的衍射峰空间分布比对分析法 |
CN106959312A (zh) * | 2017-03-10 | 2017-07-18 | 西安交通大学 | 一种基于峰位比对的扫描式劳厄衍射图谱分析方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113325015A (zh) * | 2021-05-12 | 2021-08-31 | 西安交通大学 | 劳厄衍射图谱中识别同一衍射峰所劈裂衍射亚峰的方法 |
CN113325016A (zh) * | 2021-05-12 | 2021-08-31 | 西安交通大学 | 存在衍射峰劈裂的劳厄衍射图谱中标定多个亚晶的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110927190B (zh) | 2020-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Savitzky et al. | py4DSTEM: A software package for four-dimensional scanning transmission electron microscopy data analysis | |
Kohno et al. | Real-space visualization of intrinsic magnetic fields of an antiferromagnet | |
CN110927190B (zh) | 在劳厄衍射图谱中识别同一晶带轴衍射峰的方法 | |
JP7232911B2 (ja) | 電子ビーム画像向上のための完全自動semサンプリングシステム | |
US9401013B2 (en) | Method of design-based defect classification and system thereof | |
JP5745573B2 (ja) | 製造ツールのレシピを生成する方法及びそのシステム | |
Heilmann et al. | Computerized method to determine crystal orientations from Kikuchi patterns | |
CN106950234A (zh) | 一种劳厄衍射扫描数据的衍射峰空间分布比对分析法 | |
Munshi et al. | Disentangling multiple scattering with deep learning: application to strain mapping from electron diffraction patterns | |
CN110927191B (zh) | 劳厄衍射图谱的标定方法 | |
Dingley et al. | Determination of crystal phase from an electron backscatter diffraction pattern | |
WO2020178296A1 (en) | Methods and apparatus for electron backscatter diffraction sample characterisation | |
CN111710360B (zh) | 一种预测蛋白质序列的方法、系统、装置及介质 | |
US5235523A (en) | Apparatus and methods for identifying and comparing lattice structures and determining lattice structure symmetries | |
CN116380954A (zh) | 单晶材料的二维衍射数据处理方法及系统 | |
Marquardt et al. | Diffusion in yttrium aluminium garnet at the nanometer-scale: Insight into the effective grain boundary width | |
US5168457A (en) | Apparatus for identifying and comparing lattice structures and determining lattice structure symmetries | |
Wardini et al. | Probing Multiscale Disorder in Pyrochlore and Related Complex Oxides in the Transmission Electron Microscope: A Review | |
Cabral et al. | Scanning transmission electron microscopy for advanced characterization of ferroic materials | |
Godó et al. | Residue assignment in crystallographic protein electron density maps with 3d convolutional networks | |
Thronsen et al. | Scanning precession electron diffraction data analysis approaches for phase mapping of precipitates in aluminium alloys | |
CN113325016B (zh) | 存在衍射峰劈裂的劳厄衍射图谱中标定多个亚晶的方法 | |
CN116071329B (zh) | 可视化掠入射衍射数据处理方法、系统、介质及设备 | |
US20230394852A1 (en) | Automatic selection of structures-of-interest for lamella sample preparation | |
CN113325015B (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 |