CN103591904A - 一种两步三维傅里叶变换测量物体内部三维变形场的方法 - Google Patents
一种两步三维傅里叶变换测量物体内部三维变形场的方法 Download PDFInfo
- Publication number
- CN103591904A CN103591904A CN201310426262.XA CN201310426262A CN103591904A CN 103591904 A CN103591904 A CN 103591904A CN 201310426262 A CN201310426262 A CN 201310426262A CN 103591904 A CN103591904 A CN 103591904A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- interior
- articles
- dimensional fourier
- sub
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000009466 transformation Effects 0.000 title abstract 4
- 238000006073 displacement reaction Methods 0.000 claims abstract description 36
- 238000001228 spectrum Methods 0.000 claims abstract description 23
- 239000013598 vector Substances 0.000 claims abstract description 12
- 238000002591 computed tomography Methods 0.000 claims abstract description 5
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000004587 chromatography analysis Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000010603 microCT Methods 0.000 claims description 4
- 238000001218 confocal laser scanning microscopy Methods 0.000 claims description 3
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 2
- 238000012886 linear function Methods 0.000 claims description 2
- 238000003709 image segmentation Methods 0.000 claims 2
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 238000003384 imaging method Methods 0.000 claims 1
- 238000005481 NMR spectroscopy Methods 0.000 abstract description 2
- 239000002131 composite material Substances 0.000 abstract 2
- 238000005259 measurement Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000011859 microparticle Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Landscapes
- Image Analysis (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供了一种利用两次三维傅里叶变换测量物体内部三维位移场及应变场的方法。利用微焦点计算机层析,微核磁共振或激光扫描共焦显微镜获得内部具有天然或人工制作的散斑结构的物体变形前后的三维图像。将变形前后的三维图像分别分割成若干相同尺寸的子体块,相对应的成对子块体逐一进行第一步三维傅里叶变换,在频谱域进行数值干涉得到合成谱,对合成谱进行第二步的三维傅里叶变换,得到含位移信息的扩展脉冲函数,通过求取峰值位置得到相对应子块体的位移矢量,由所有子块***移矢量获得物体内部三维变形场。通过插值获得位移的亚体素精度。由应变与变形间的关系方程,求出物体内部的应变场。
Description
技术领域
本发明涉及物体内部三维变形场的测量,属固体实验力学领域。
背景技术
数字散斑相关法或数字图像相关法(Digital Image Correlation,DIC)已广泛用于二维平面变形及应变的测量。数字体相关法(Digital Volume Correlation,DVC)是数字图像相关法向物体内部三维变形测量的拓展,该方法利用微焦点CT,微焦点核磁共振或激光扫描共焦显微镜获取得变形前后物体的三维立体图像,直接在变形后的体图像中搜索与变形前体图像中相对应的点,得出点的位移矢量,由于要在搜索区域中进行逐一搜索,致使计算时间消耗较大,尤其当物体产生大变形,搜索区域增大,耗费的时间将大为增加,计算效率降低。
为提高计算效率,本发明提供一种利用两步三维傅里叶变换测量物体内部三维变形场的方法。
发明内容
本发明的目的是提供一种高效地测量计算内部具有散斑结构物体的三维变形及应变场的方法。
本发明用于测量计算内部具有散斑结构物体的三维变形及应变场,物体可以是透明的或不透明的,散斑结构可以是物体内部的天然结构所形成的,如骨骼、木材、孔隙类岩石或含多种成份的岩石等,也可以是人为在待研究物体中加入某种在图像中易于辨识的微颗粒形成散斑结构。
本发明利用微焦点计算机层析(Micro focus computer tomography,Micro-CT),微核磁共振成像(Micronuclear magnetic resonance imaging,Micro-MRI),激光扫描共焦显微镜(Laser Scanning Confocal Microscope,LSCM)或其他可以获取反映物体内部结构的三维图像的方法获取具有散斑结构物体变形前后的三维立体图像,所得立体图像分别定义为参考体图像和变形体图像,将参考体图像与变形体图像分割成若干大小相等的子体块,参考体图像与变形体图像中相对应的一对子体块构成“斑对”,对构成“斑对”的两个子体块分别进行第一步三维傅里叶变换,得到相应的频谱函数,在频谱域中对所得到的频谱函数进行数值干涉,得到合成光谱函数,对该函数进行第二步傅里叶变换,得到一个峰值中心偏移的扩展脉冲函数,通过求取该扩展脉冲函数峰值的坐标,得到参考子体块的位移矢量。对参考体图像与变形体图像分割后所有相对应的“斑对”逐一进行上述两步三维傅里叶变换,就可以得到物体内部的变形场。通过应变与变形的关系,可以求解物体内部应变场。
本发明的具体步骤是:
步骤一:利用微焦点计算机层析(Micro focus computer tomography,Micro-CT),微核磁共振成像(Micronuclear magnetic resonance imaging,Micro-MRI),激光扫描共焦显微镜(Laser Scanning Confocal Microscope,LSCM)或其他可以获取反映物体内部结构的三维图像的方法获取具有散斑结构物体变形前后的三维立体图像,分别定义为参考体图像和变形体图像;
步骤二:将参考体图像与变形体图像分割成若干大小如32×32×32体素的子体块,相对应的参考子体块和变形子体块构成一个“斑对”;当子体块尺寸较小时,可以认为子体块的内各点的位移是相同的;
本步骤中子块体的尺寸可以是4×4×4体素,8×8×8体素,16×16×16体素,32×32×32体素,64×64×64体素,……,根据物体内部散斑结构的大小,适当选取子块体的大小,但每个子块体中散斑结构的数量应大于9。
本步骤中体图像的分割相邻子体块之间可以有二分之一子块体尺寸的重叠,提高计算的准确性。
步骤三:令h1(x,y,z)和h2(x,y,z)分别表示一个“斑对”中的参体子体块和变形子体块,h1(x,y,z)与h2(x,y,z)具有如下关系:
h2(x,y,z)=h1[x-u(x,y,z),y-v(x,y,z),z-w(x,y,z)] (1)
式中,u,v,w分别为变形后参考子体块在x,y,z三个方向上的位移,对参考子体块和变形子体块分别进行第一步三维傅里叶变换,则有:
式是Δ表示子体块区域,(fx,fy,fz)为第一次三维傅里叶变换后的频谱域坐标,H1(fx,fy,fz)为h1(x,y,z)傅里叶变换后的函数,|H1(fx,fy,fz)|和φ1(fx,fy,fz)分别为H1(fx,fy,fz)的幅值与相位,H2(fx,gy,fz)为h2(x,y,z)傅里叶变换后的函数,|H2(fx,fy,fz)|和φ2(fx,fy,fz)分别为H2(fx,fy,fz)的幅值与相位;
步骤四:对步骤三中得到了H1(fx,fy,fz)和H2(fx,fy,fz)通过下式进行数值干涉生成新的合成光谱
也可以写成
F(fx,fy,fz)=|H1(fx,fy,fz)H2(fx,fy,fz)|αexp{j[φ1(fx,fy,fz)-φ2(fx,fy,fz)]}
(4)
式中0≤α≤1。
步骤五:对步骤四中所得的合成光谱进行第二步傅里叶变换。
将式(2)代入式(3)中可得
F(fx,fy,fz)≈|H1(fx,fy,fz)|2αexp[j2π(ufx+vfy+wfz)] (5)
式中Δf表示子体块经第一步三维傅里叶变换后的频谱区域,(ξ,η,ζ)表示第二步三维傅里叶变换后的频谱域的坐标。
是一个峰值点在(u,v,w)处的扩展脉冲函数。
步骤六:搜索求取步骤五中扩展脉冲函数峰值点的坐标值(u,v,w),从而得到参考子体块的位移矢量;
步骤七:在数字图像中,步骤六中所求得的三个方向的位移均为整体素,为了提高测量精度,通过插值获得位移矢量的亚体素精度。在扩展脉冲函数中,围绕着步骤六中所求得的峰值点(其坐标为整体素),选取3×3×3或是5×5×5的矩阵,可采用已有的三维线性(trilinear),三维多义线(tri-spline)或三维三次多项式插值(tri-cubic)等插值方法来求取亚体素精度的位移;
步骤八:对参考体图像与变形体图像分割成的若干子体块形成的“斑对”,逐一进行上述步骤三至步骤七的操作,最终获得整个物体内部的变形场。
步骤九:计算物体内部的应变场,由步骤八得到的物体内部变形场,通过应变与变形的关系方程,可以求解物体内部的应变场。本发明中采用逐点最小二乘法来求解应变。由步骤八所得物体内部的变形场,以待求点为中心,围绕其选取(2N+1)×(2N+1)×(2N+1)立方区域,当所选区域足够小时,其内部各点的位移可以看作是坐标的线性函数,因此有:
u(x,y,z)=a0+a1x+a2y+a3z
v(x,y,z)=b0+b1x+b2y+b3z
w(x,y,z)=c0+c1x+c2y+c3z
式中,x,y,z取值范围为[-N,N],且为整数。u(x,y,z),v(x,y,z)和w(x,y,z)为步骤八中所求得的物体内部变形。ai=0,1,2,3,bi=0,1,2,3和ci=0,1,2,3是待定的多项式系数,可以通过已有的最小二乘法或多元拟合方法来求取。这些系数确定后,待求点的Cauchy应变值为:
附图说明
图1为本发明的原理示意图;
其中1.参考体图像,2.变形体图像,3.分割后的参考子体图像矩阵,4.分割后的变形子体图像矩阵,5.参考子体块,6.变形子体块,7.扩展脉冲函数的等值分布
图2为中砂岩变形后的内部三维位移分量,x方向的位移u场截面图;
图3为中砂岩变形后的内部三维位移分量,y方向的位移v场截面图;
图4为中砂岩变形后的内部三维位移分量,z方向的位移w场截面图;
具体实施方式
下面结合附图与实施例对本发明作进一步说明
步骤一:利用微焦点计算机层析(Micro focus computer tomography,Micro-CT)扫描获得中砂岩试件变形前后的三维立体图像,分别定义为参考体图像1和变形体图像2;
步骤二:将参考体图像1和变形体图像2分别分割成若干大小如32×32×32体素的子体块矩阵3和4,相对应的每个子块体如5和6构成一个“斑对”,令h1(x,y,z)和h2(x,y,z)分别表示变形前参考子体块5及变形后相应的子体块6,h1(x,y,z)与h2(x,y,z)具有如下关系:
h2(x,y,z)=h1[x-u(x,y,z),y-v(x,y,z),z-w(x,y,z)] (1)
式中,u,v,w分别为子块体5在x,y,z三个方向上的位移,
步骤三:将构成“斑对”的变形前后的子块体5和6分别进行第一步三维傅里叶变换,则有:
步骤四:在频谱域中,对上述两个三维散斑体频谱H1和H2进行数值干涉,得到合成谱如下:
式中*表示共轭,Φ1(fx,fy,fz)和Φ2(fx,fy,fz)分别为H1(fx,fy,fz)和H2(fx,fy,fz)的相位角,由相位差与光程差之间的关系则有:
φ1(fx,fy,fz)-φ2(fx,fy,fz)=2π(ufx+vfy+wfz) (4)
步骤五:对数值干涉得到的合成谱F(fx,fy,fz)进行第二次三维傅里叶变换,有:
上式表示中心峰值点位于(u,v,w)的扩展脉冲函数,其三维分布如图1中7所示。
步骤六:通过搜索步骤五所得扩展脉冲函数的峰值点坐标,得到参考子体块5的位移矢量;
步骤七:为了提高测量精度,通过插值获得位移矢量的亚体素精度。在扩展脉冲函数中,围绕着步骤六中所求得的峰值点坐标(为整体素),选取3×3×3或是5×5×5的矩阵,可采用已有的三维线性(tri-linear),三维多义线(tri-spline)或三维三次多项式插值(tri-cubic)等方法来求取亚体素精度的位移;
步骤八:对所有子块体形成的“斑对”进行上述步骤三到步骤七的运算,最终获得整个物体内部的变形场。图2为中砂岩变形后的内部三维位移分量,x方向的位移u场截面图;图3为中砂岩变形后的内部三维位移分量,y方向的位移v场截面图;图4为中砂岩变形后的内部三维位移分量,z方向的位移w场截面图;
步骤九:计算物体内部的应变场,由步骤八得到的物体内部的变形场,通过应变与变形的关系方程,可以求解物体内部的应变场。
Claims (10)
1.一种利用两步三维傅里叶变换测量物体内部三维变形的方法,其特征是:利用计算机层析技术,核磁共振成像技术或激光扫描共焦显微镜等方法获取具有散斑结构特征物体在变形前后的三维立体图像,将物体变形前后的三维立体图像分割成一系列子体块,相对应的子体块构成“斑对”,分别进行第一步三维傅里叶变换,在频谱域中进行数值干涉,对所得合成光谱进行第二次三维傅里叶变换,得到该子体块的位移矢量,对所有相对应的子体块进行上述两步三维傅里叶变换,就可以求出物体内各子体块的位移矢量。根据应变与变形的关系方程,可以获得物体内部应变场。
2.根据权利要求1所述的一种两步三维傅里叶变换测量物体内部三维变形场的方法,该物体可以是透明的或是非透明的,但其内部应具有天然或人为制作的结构作为散斑结构。
3.根据权利要求1所述的一种两步三维傅里叶变换测量物体内部三维变形场的方法,物体的三维立体图像可以由微焦点计算机层析(Micro focus computer tomography,Micro-CT),微核磁共振成像(Micro nuclear magnetic resonance imaging,Micro-MRI)或激光扫描共焦显微镜(Laser Scanning Confocal Microscope,LSCM)或其他可以获取反映物体内部结构的三维图像的方法获取。
4.根据权利要求1所述的一种两步三维傅里叶变换测量物体内部三维变形场的方法,在将参考体图像与变形体图像进行分割时,子块体的尺寸可以是4×4×4体素,8×8×8体素,16×16×16体素,32×32×32体素,64×64×64体素……,根据物体内部散斑结构的大小,适当选取子块体的大小,但每个子块体中散斑结构的数量应大于9。在将参考体图像与变形体图像进行分割时,相邻子体块之间可以有二分之一子块体尺寸的重叠,提高计算的准确性。
5.根据权利要求1所述的一种两步三维傅里叶变换测量物体内部三维变形场的方法,参考体图像与变形体图像分割后相对应的子体块构成“斑对”,分别进行第一步三维傅里叶变换。
令h1(x,y,z)和h2(x,y,z)分别表示变形前参考子体块及其变形后相应变形子体块,h1(x,y,z)与h2(x,y,z)具有如下关系:
h2(x,y,z)=h1[x-u(x,y,z),y-v(x,y,z),z-w(x,y,z)]
式中,u,v,w分别为参考子块体在x,y,z三个方向上的位移,
将变形前后的子块体进行第一步三维傅里叶变换,有:
7.根据权利要求1所述的一种两步三维傅里叶变换测量物体内部三维变形场的方法,对数值干涉得到的合成光谱F(fx,fy,fz)进行第二次三维傅里叶变换,有:
上式表示中心峰值点位于点(u,v,w)的扩展脉冲函数。通过搜索得到扩展脉冲函数的峰值点坐标,从而得到参考子体块的位移矢量。
9.根据权利要求1所述的一种两步三维傅里叶变换测量物体内部三维变形场的方法,对所有的子块体形成的“斑对”进行两步三维傅里叶变换的运算,最终获得整个物体内部的变形场。
10.根据权利要求1所述的一种两步三维傅里叶变换测量物体内部三维变形场的方法,通过应变与变形之间的关系方程,由得到的物体内部变形场可以求解物体内部的应变场。由所求得物体内部的变形场,以待求点为中心,围绕其选取(2N+1)×(2N+1)×(2N+1)立方区域,当所选区域足够小时,其内部各点的位移可以看作是坐标的线性函数,因此有:
u(x,y,z)=a0+a1x+a2y+a3z
v(x,y,z)=b0+b1x+b2y+b3z
w(x,y,z)=c0+c1x+c2y+c3z
式中,x,y,z取值范围为[-N,N],且为整数。u(x,y,z),v(x,y,z)和w(x,y,z)是由两步三维傅里叶变换所求得的物体内部变形。ai=0,1,2,3,bi=0,1,2,3和ci=0,1,2,3是待定的多项式系数,可以通过已有的最小二乘法或多元拟合方法来求取。这些系数确定后,待求点的Cauchy应变值为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310426262.XA CN103591904B (zh) | 2013-09-18 | 2013-09-18 | 一种两步三维傅里叶变换测量物体内部三维变形场的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310426262.XA CN103591904B (zh) | 2013-09-18 | 2013-09-18 | 一种两步三维傅里叶变换测量物体内部三维变形场的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103591904A true CN103591904A (zh) | 2014-02-19 |
CN103591904B CN103591904B (zh) | 2016-08-17 |
Family
ID=50082137
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310426262.XA Expired - Fee Related CN103591904B (zh) | 2013-09-18 | 2013-09-18 | 一种两步三维傅里叶变换测量物体内部三维变形场的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103591904B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103900506A (zh) * | 2014-03-28 | 2014-07-02 | 东南大学 | 一种基于位移梯度分解的异位数字体积相关方法 |
CN103927786A (zh) * | 2014-04-29 | 2014-07-16 | 盐城工学院 | 一种刀具刃口三维模型构建方法 |
CN104914118A (zh) * | 2015-06-09 | 2015-09-16 | 河海大学 | 一种基于ct重建材料构件内部位移场的测量*** |
CN111060527A (zh) * | 2019-12-30 | 2020-04-24 | 歌尔股份有限公司 | 一种字符缺陷检测方法及装置 |
CN113624358A (zh) * | 2021-07-14 | 2021-11-09 | 中国电子科技集团公司第十三研究所 | 用于光热反射显微热成像的三维位移补偿方法及控制装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007240465A (ja) * | 2006-03-10 | 2007-09-20 | Wakayama Univ | 3次元変位ひずみ計測方法及び装置 |
CN101980304A (zh) * | 2010-10-20 | 2011-02-23 | 北京大学 | 一种三维数字体积图像变形测量方法 |
-
2013
- 2013-09-18 CN CN201310426262.XA patent/CN103591904B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007240465A (ja) * | 2006-03-10 | 2007-09-20 | Wakayama Univ | 3次元変位ひずみ計測方法及び装置 |
CN101980304A (zh) * | 2010-10-20 | 2011-02-23 | 北京大学 | 一种三维数字体积图像变形测量方法 |
Non-Patent Citations (3)
Title |
---|
潘兵 等: "基于梯度的数字体图像相关方法测量物体内部变形", 《光学学报》, vol. 31, no. 6, 30 June 2011 (2011-06-30) * |
潘兵 等: "数字体图像相关方法中基于迭代最小二乘法的物体内部变形测量", 《实验力学》, vol. 26, no. 6, 31 December 2011 (2011-12-31), pages 665 - 673 * |
胡宇明 冉英华: "一种基于快速傅立叶变换的数字散斑***", 《大众科技》, no. 4, 31 December 2011 (2011-12-31), pages 38 - 39 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103900506A (zh) * | 2014-03-28 | 2014-07-02 | 东南大学 | 一种基于位移梯度分解的异位数字体积相关方法 |
CN103900506B (zh) * | 2014-03-28 | 2016-08-17 | 东南大学 | 一种基于位移梯度分解的异位数字体积相关方法 |
CN103927786A (zh) * | 2014-04-29 | 2014-07-16 | 盐城工学院 | 一种刀具刃口三维模型构建方法 |
CN104914118A (zh) * | 2015-06-09 | 2015-09-16 | 河海大学 | 一种基于ct重建材料构件内部位移场的测量*** |
CN111060527A (zh) * | 2019-12-30 | 2020-04-24 | 歌尔股份有限公司 | 一种字符缺陷检测方法及装置 |
CN111060527B (zh) * | 2019-12-30 | 2021-10-29 | 歌尔股份有限公司 | 一种字符缺陷检测方法及装置 |
US12002198B2 (en) | 2019-12-30 | 2024-06-04 | Goertek Inc. | Character defect detection method and device |
CN113624358A (zh) * | 2021-07-14 | 2021-11-09 | 中国电子科技集团公司第十三研究所 | 用于光热反射显微热成像的三维位移补偿方法及控制装置 |
CN113624358B (zh) * | 2021-07-14 | 2023-09-26 | 中国电子科技集团公司第十三研究所 | 用于光热反射显微热成像的三维位移补偿方法及控制装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103591904B (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pan | Digital image correlation for surface deformation measurement: historical developments, recent advances and future goals | |
CN109163672B (zh) | 一种基于白光干涉零光程差位置拾取算法的微观形貌测量方法 | |
White et al. | Transverse coronal loop oscillations seen in unprecedented detail by AIA/SDO | |
CN103591904A (zh) | 一种两步三维傅里叶变换测量物体内部三维变形场的方法 | |
CN102621549B (zh) | 多基线/多频段干涉相位解缠频域快速算法 | |
CN102788558B (zh) | 散斑相关和散斑干涉相结合的三维变形测量***及方法 | |
CN103471520B (zh) | 面结构光与锥光偏振全息组合的反光复杂曲面测量方法 | |
CN106289109A (zh) | 一种基于结构光的三维重建***及方法 | |
CN104864819A (zh) | 一种基于数字散斑的高速三维应变测量方法 | |
CN101655357B (zh) | 一种用于二维相位展开的相位梯度相关质量图获取方法 | |
CN101915559B (zh) | 利用电子散斑相移技术测量物体三维面形的方法及其*** | |
Brücker et al. | Single-view volumetric PIV via high-resolution scanning, isotropic voxel restructuring and 3D least-squares matching (3D-LSM) | |
CN106032976A (zh) | 基于波长选择的三条纹投影相位展开方法 | |
Cai et al. | Automatic target recognition based on alignments of three-dimensional interferometric ISAR images and CAD models | |
CN101980304A (zh) | 一种三维数字体积图像变形测量方法 | |
CN110161501A (zh) | 一种多时序sar图像的目标区域地表起伏信息提取方法 | |
CN105043301A (zh) | 一种用于三维测量的光栅条纹相位求解方法 | |
CN111947599A (zh) | 基于学习的条纹相位恢复和散斑相关的三维测量方法 | |
Li et al. | Digital speckle correlation method based on phase vortices | |
Muller et al. | A Statistical investigation of HI in the Magellanic Bridge | |
CN109738852A (zh) | 基于低秩矩阵重建的分布式源二维空间谱估计方法 | |
CN104121864B (zh) | 一种岩石错动裂缝面的间隙评价方法 | |
CN105043270A (zh) | 混合分形插值和频域数字散斑相关的mems面内位移测量方法 | |
Singh et al. | Nanoscale measurement of in-plane and out-of-plane displacements of microscopic object by sensor fusion | |
Metzner et al. | A comparison of two data analysis approaches for quantitative magnetic resonance imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
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: 20160817 Termination date: 20190918 |