CN113405628A - 基于双层图形透射成像的液面变形和形貌测量方法和系统 - Google Patents
基于双层图形透射成像的液面变形和形貌测量方法和系统 Download PDFInfo
- Publication number
- CN113405628A CN113405628A CN202110691868.0A CN202110691868A CN113405628A CN 113405628 A CN113405628 A CN 113405628A CN 202110691868 A CN202110691868 A CN 202110691868A CN 113405628 A CN113405628 A CN 113405628A
- Authority
- CN
- China
- Prior art keywords
- image
- information carrier
- liquid level
- layer
- displacement vector
- 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
- 239000007788 liquid Substances 0.000 title claims abstract description 166
- 230000005540 biological transmission Effects 0.000 title claims abstract description 26
- 238000003384 imaging method Methods 0.000 title claims abstract description 21
- 238000000691 measurement method Methods 0.000 title claims abstract description 15
- 238000006073 displacement reaction Methods 0.000 claims abstract description 112
- 239000013598 vector Substances 0.000 claims abstract description 105
- 238000000034 method Methods 0.000 claims abstract description 65
- 238000005516 engineering process Methods 0.000 claims abstract description 28
- 230000003068 static effect Effects 0.000 claims abstract description 26
- 238000005259 measurement Methods 0.000 claims abstract description 19
- 239000010410 layer Substances 0.000 claims description 150
- 238000012545 processing Methods 0.000 claims description 44
- 238000004458 analytical method Methods 0.000 claims description 19
- 238000012876 topography Methods 0.000 claims description 14
- 238000000605 extraction Methods 0.000 claims description 10
- 239000011521 glass Substances 0.000 claims description 9
- 229920003023 plastic Polymers 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 239000002355 dual-layer Substances 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 13
- 230000003287 optical effect Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 239000000969 carrier Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 241001489698 Gerridae Species 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000010587 phase diagram Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000005290 field theory Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 239000003380 propellant Substances 0.000 description 1
- 238000001028 reflection method Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 230000009131 signaling function Effects 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F23/00—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
- G01F23/22—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water
- G01F23/28—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water by measuring the variations of parameters of electromagnetic or acoustic waves applied directly to the liquid or fluent solid material
- G01F23/284—Electromagnetic waves
- G01F23/292—Light, e.g. infrared or ultraviolet
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Thermal Sciences (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于双层图形透射成像的液面变形和形貌测量方法和系统，方法包括：提供图像采集设备；提供具有一定厚度的介质；液面静止时，利用图像获取方法通过图像采集设备对信息载体图形进行图像采集，得到液面静止时的信息载体图像；当液面发生形变时，利用图像获取方法对信息载体图形进行图像采集，得到液面形变后的信息载体畸变图像；求解上层信息载体图像各个时刻的全场位移矢量场与下层信息载体图像各个时刻的全场位移矢量场；求解液面各个时刻的液体高度数据。本发明通过记录液面静止时和液面变形后各时刻的双层信息载体图像，运用相关技术，可以实现对静态和动态液面形貌进行快速高精度地测量。
Description
技术领域
本发明涉及光学测量领域，更具体地，涉及一种基于双层图形透射成像的液面变形和形貌测量方法和系统。
背景技术
液体自由表面受到扰动产生的变形现象，既有重要的科学研究价值，也有实际的工程应用。其中以水黾的漂浮为代表的固体漂浮现象是重力、浮力与表面张力共同作用的结果，而水黾水上运动过程中水对水黾的作用力如何定量分析目前还没有得到明确的结论。不同激励下水波的传递与变形是流体力学中的重要研究方向，尤其是在航空航天领域，携带的大量液体推进剂在任务中会不可避免的发生晃动，容易导致航天器姿态失稳，也因此液体晃动成为重要的研究方向。由于液面的形貌与浮力、表面张力、水波与液体晃动等直接相关，因此完成对液面变形的表征能显著帮助相关研究方向的进展。
现有的液面形貌测量方法主要有两类方法：一类是透射方法，一类是反射方法。透射方法中，J.C.Wyant使用全息剪切干涉法(J.C.Wyant,Appl Opt,1973；12:2057-60)与K.D.Hinsch使用全息干涉法(K.D.Hinsch,Appl Opt,1978；17:3101-7)用于液面测量，变形测量灵敏度达到了纳米量级，但是测量量程较小，无法测量毫米及以上量级的液面变形；汤慧颖等人(Optics and Lasers in Engineering,2017,98:205-216)提出的双色点阵偏折法能测量微米到毫米量级的单层与双层液面变形，但是不适用于大斜率的液面变形；刘战伟等人(Optics and Lasers in Engineering,2012；51(2)：167-171)发展的透射条纹法，专利(申请号为CN201310355784.5)“基于数字散斑相关技术的透明液面微形貌测量方法和系统”与专利(申请号为CN201410429656.5)“贮箱液体推进剂液面形貌与剂量动态测量方法和系统”等方法由于均采用了迭代算法，无法避免误差的累积；Qian等人(Proceedingsof the IEEE conference on computer vision and pattern recognition.2017:1269-1278)提出的光线追踪法测量精度非常高，但需要多个相机，同时计算复杂，效率低。A.K.Asundi等通过反射法测量了动态液面的三维形貌(Huang L,A.K.Asundi等,OptExpress,2011；12809-14)，但是液面变形较大会导致反射光线超出相机视场。
因此，如何在保证测量灵敏度的同时，减少迭代误差，增大测量量程，减少装置制作成本，提高测量效率成为亟待解决的问题。
发明内容
有鉴于此，本发明提供了一种基于双层图形透射成像的液面变形和形貌测量方法与系统，以解决现有的液面形貌测量方法难以同时实现高精度，高效率，装置经济，操作简单的问题。
一方面，本发明公开了一种基于双层图形透射成像的液面变形和形貌测量方法，包括步骤：
提供图像采集设备，将图像采集设备设置在液面远离地面的一侧，且所述图像采集设备与所述液面具有间隔，并对所述图像采集设备进行标定；
提供具有一定厚度的介质，并将所述介质置于所述液面底部，在垂直于地面所在平面的方向上，所述介质包括相对设置的上层和下层，所述上层和下层均具有信息载体图形；
所述液面静止时，利用图像获取方法通过所述图像采集设备对信息载体图形进行图像采集，得到液面静止时的信息载体图像，所述图像获取方法包括所述图像采集设备采集双层信息载体图像，上层信息载体图形和下层信息载体图形的图像采集到同一图像中，利用几何相位分析技术或者彩色阈值提取法将上层信息载体图形的图像和下层信息载体图形的图像分开到两张图像中，形成上层信息载体参考图像与下层信息载体参考图像；
当所述液面发生形变时，利用所述图像获取方法对信息载体图形进行图像采集，得到液面形变后的信息载体畸变图像，所述信息载体畸变图像包括上层信息载体畸变图像和下层信息载体畸变图像，采集一张信息载体畸变图像时，最终测量的是采集时刻的液面瞬时形貌，当采集随时间变化的序列信息载体畸变图像时，最终测量得到的是液面随时间变化的形貌；
以液面静止时获取的上层信息载体参考图像与下层信息载体参考图像作为参考，对所述上层信息载体畸变图像和所述下层信息载体畸变图像进行处理得到所述上层信息载体图像和所述下层信息载体图像各个时刻的面内U、V方向的位移矢量数据，求解上层信息载体图像各个时刻的全场位移矢量场与下层信息载体图像各个时刻的全场位移矢量场；
按照以下方法，根据上层信息载体图像的全场位移矢量场与下层信息载体图像的全场位移矢量场求解所述液面各个时刻的液体高度数据：
其中，S1为上层信息载体图像的全场面内位移，S2为下层信息载体图像的全场面内位移，d为所述介质的上层与下层之间的厚度，α为图像采集设备所采集图像对应位置的光线与竖直方向的夹角，n1为液体的折射率，n2为所述介质的折射率。
可选的，所述信息载体图形包括条纹、散斑和点阵信息。
可选的，采用数字云纹技术、数字散斑技术或几何相位分析技术对所述对上层信息载体畸变图像和所述下层信息载体畸变图像进行处理得到所述上层信息载体图像和所述下层信息载体图像各个时刻的面内U、V方向的位移矢量数据。
可选的，所述介质包括透明玻璃、或透明塑料，所述介质的折射率与所述液体的折射率相等或不相等。
可选的，所述上层信息载体图像各个时刻的全场位移矢量场，包括所述上层信息载体图像各个时刻的面内U、V方向位移矢量场的矢量和，所述下层信息载体图像各个时刻的全场位移矢量场包括下层信息载体图像各个时刻的面内U、V方向位移矢量场的矢量和。
另一方面，本发明还公开了一种基于双层图形透射成像的液面变形和形貌测量系统，包括图像采集设备、介质、面内位移矢量处理模块以及液面三维形貌处理模块，其中，
所述介质具有一定厚度，设置在液面底部，在垂直于地面所在平面的方向上，所述介质包括相对设置的上层和下层，所述上层和下层均具有信息载体图形；
所述图像采集设备，设置在液面远离地面的一侧，且所述图像采集设备与所述液面具有间隔，并对所述图像采集设备进行标定，与所述面内位移矢量处理模块相耦接，用于所述液面静止时利用图像获取方法对信息载体图形进行图像采集，得到液面静止时的信息载体图像，所述图像获取方法包括：所述图像采集设备采集双层信息载体图像，上层信息载体图形和下层信息载体图形的图像采集到同一图像中，利用几何相位分析技术或者彩色阈值提取法将上层信息载体图形的图像和下层信息载体图形的图像分开到两张图像中，形成上层信息载体参考图像与下层信息载体参考图像；当所述液面发生形变后，对信息载体图形进行图像采集，得到液面形变后的信息载体畸变图像，所述信息载体畸变图像包括上层信息载体畸变图像和下层信息载体畸变图像，并发送至所述面内位移矢量处理模块，采集一张信息载体畸变图像时，最终测量的是采集时刻的液面瞬时形貌，当采集随时间变化的序列信息载体畸变图像时，最终测量得到的是液面随时间变化的形貌；
所述面内位移矢量处理模块，分别与所述图像采集设备和所述液面三维形貌处理模块相耦接，用于以液面静止时获取的上层信息载体参考图像与下层信息载体参考图像作为参考，对所述上层信息载体畸变图像和所述下层信息载体畸变图像进行处理得到所述上层信息载体图像和所述下层信息载体图像各个时刻的面内U、V方向的位移矢量数据，求解上层信息载体图像各个时刻的全场位移矢量场与下层信息载体图像各个时刻的全场位移矢量场，并发送至所述液面三维形貌处理模块；
所述液面三维形貌处理模块，与所述面内位移矢量处理模块相耦接，用于接收所述面内位移矢量处理模块发送的上层信息载体图像的全场位移矢量场与下层信息载体图像的全场位移矢量场数据，按照以下方法求解出液面高度:
其中，S1为上层信息载体图像的全场面内位移，S2为下层信息载体图像的全场面内位移，d为所述介质的上层与下层之间的厚度，α为图像采集设备所采集图像对应位置的光线与竖直方向的夹角，n1为液体的折射率，n2为所述介质的折射率。
可选的，所述信息载体图形包括条纹、散斑和点阵信息。
可选的，所述面内位移矢量处理模块采用数字云纹技术、数字散斑技术或几何相位分析技术对所述对上层信息载体畸变图像和所述下层信息载体畸变图像进行处理得到所述上层信息载体图像和所述下层信息载体图像各个时刻的面内U、V方向的位移矢量数据。
可选的，所述介质包括透明玻璃、或透明塑料，所述介质的折射率与所述液体的折射率相等或不相等。
可选的，所述上层信息载体图像各个时刻的全场位移矢量场，包括所述上层信息载体图像各个时刻的面内U、V方向位移矢量场的矢量和，所述下层信息载体图像各个时刻的全场位移矢量场包括下层信息载体图像各个时刻的面内U、V方向位移矢量场的矢量和。
与现有技术相比，本发明提供的基于双层图形透射成像的液面变形和形貌测量方法和系统，至少实现了如下的有益效果：
1)本发明通过记录液面静止时和液面变形后各时刻的双层信息载体图像，运用相关技术，可以实现对静态和动态液面形貌进行快速高精度地测量；
2)本发明通过分析上下两层信息载体图像畸变与液面高度之间的光学关系，有效提高测量量程的同时提高了测量的灵敏度；
3)本发明测量上下两层信息载体图形反射的光线透射出变形液面后的畸变，结合该上下两层信息载体图形间给定的间隔距离，就可以综合得到液面形貌变化图，测量方法简单，计算效率高，并且便于实验操作。
4)本发明只需要使用单个的图像采集设备，装置简单经济。
当然，实施本发明的任一产品必不特定需要同时达到以上所述的所有技术效果。
通过以下参照附图对本发明的示例性实施例的详细描述，本发明的其它特征及其优点将会变得清楚。
附图说明
被结合在说明书中并构成说明书的一部分的附图示出了本发明的实施例，并且连同其说明一起用于解释本发明的原理。
图1为本发明实施例一所述的液面三维形貌测量方法的流程图；
图2为本发明实施例一所述的液面三维形貌测量方法的原理图；
图3是本发明利用图像获取方法得到的一种双层信息载体图像；
图4是利用彩色阈值提取法得到的上层信息载体图像；
图5是利用彩色阈值提取法得到的下层信息载体图像；
图6是本发明利用图像获取方法得到的又一种双层信息载体图像；
图7和图8是利用几何相位分析技术得到的上层信息载体图像；
图9和图10是利用几何相位分析技术得到的下层信息载体图像；
图11为本发明实施例二所述的液面三维形貌测量系统的结构框图。
具体实施方式
现在将参照附图来详细描述本发明的各种示例性实施例。应注意到：除非另外具体说明，否则在这些实施例中阐述的部件和步骤的相对布置、数字表达式和数值不限制本发明的范围。
以下对至少一个示例性实施例的描述实际上仅仅是说明性的，决不作为对本发明及其应用或使用的任何限制。
对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论，但在适当情况下，所述技术、方法和设备应当被视为说明书的一部分。
在这里示出和讨论的所有例子中，任何具体值应被解释为仅仅是示例性的，而不是作为限制。因此，示例性实施例的其它例子可以具有不同的值。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步讨论。
实施例一：
结合图1，本发明提供了一种基于双层图形透射成像的液面变形和形貌测量方法流程图，具有以下步骤：
步骤101：提供图像采集设备，将图像采集设备设置在液面远离地面的一侧，且图像采集设备与液面具有间隔，并对图像采集设备进行标定；
结合图2，各层信息载体图像的采集可以通过CCD(charge-coupled Device，电荷耦合元件，即图像传感器)相机和/或将图像转变为数字信号功能的图像采集设备实现。图像采集设备(如CCD相机)置于水槽的正上方。
步骤102：提供具有一定厚度的介质，并将介质置于液面底部，在垂直于地面所在平面的方向上，介质包括相对设置的上层和下层，上层和下层均具有信息载体图形；
在一些可选的实施例中，介质包括透明玻璃、或透明塑料，介质的折射率与液体的折射率可以相等，也可以不相等。
结合图2，在实际应用中，将各层信息载体图形置于盛有适量透明液体的玻璃水槽的外面底部。此时，信息载体图形在该方法中就是置于液面底部或下方，就可以应用光的折射定律，即斯涅尔定律。
信息载体图形为：条纹，散斑，点阵等可以根据畸变分析面内位移的图形。
步骤103：液面静止时，利用图像获取方法通过图像采集设备对信息载体图形进行图像采集，得到液面静止时的信息载体图像，图像获取方法包括图像采集设备采集双层信息载体图像，上层信息载体图形和下层信息载体图形的图像采集到同一图像中，利用几何相位分析技术或者彩色阈值提取法将上层信息载体图形的图像和下层信息载体图形的图像分开到两张图像中，形成上层信息载体参考图像与下层信息载体参考图像；
利用图像获取方法时，可以通过控制上下层信息载体图像的颜色不同，通过编程实现颜色阈值提取将上下层信息载体图像分开到两张图像中，形成上层信息载体图像与下层信息载体图像。或者通过控制上下层信息载体图像相位分布动不同，通过几何相位分析技术将上下层信息载体图像分开到两张图像中，形成上层信息载体图像与下层信息载体图像。
具体地可结合图3至图10，图3中采集到的是双层信息载体图像，上层信息载体图形和下层信息载体图形的图像均被采集到了图3中，此时可利用彩色阈值提取法将上层信息载体图形和下层信息载体图形的图像分开得到两张图像，图4和图5，即上层信息载体图像和下层信息载体图像；当然也可以采用几何相位分析技术，图6中是采集到的双层信息载体图像，上层信息载体图形和下层信息载体图形的图像均被采集到了图6中，此时可利用几何相位分析技术将上层信息载体图形和下层信息载体图形的图像分开得到两张图像，图7和图8是上层信息载体图像的两个方向的相位图，图9和图10是下层信息载体图像的两个方向的相位图。
步骤104：当液面发生形变时，利用图像获取方法对信息载体图形进行图像采集，得到液面形变后的信息载体畸变图像，信息载体畸变图像包括上层信息载体畸变图像和下层信息载体畸变图像，采集一张信息载体畸变图像时，最终测量的是采集时刻的液面瞬时形貌，当采集随时间变化的序列信息载体畸变图像时，最终测量得到的是液面随时间变化的形貌；
当然，这里的图像获取方法可参照步骤103中的方法，这里不再赘述。
信息载体畸变图像为液面变形后液面(水)底的信息载体图形将有效光源反射出液面后被图像采集设备采集到的图像。
步骤105：以液面静止时获取的上层信息载体参考图像与下层信息载体参考图像作为参考，对上层信息载体畸变图像和下层信息载体畸变图像进行处理得到上层信息载体图像和下层信息载体图像各个时刻的面内U、V方向的位移矢量数据，求解上层信息载体图像各个时刻的全场位移矢量场与下层信息载体图像各个时刻的全场位移矢量场；
对各层信息载体图像中的畸变运用相关技术进行处理以获取各层信息载体图像上各个时刻的面内U、V方向的位移矢量数据。这里的UV方向与OXY坐标轴方向相同(平行)，只是在场论中习惯于采用‘U、V’叫法。
在本实施例中，液面发生变形后，各层信息载体图形反射的光经过液面的变形区域会发生透射角度的偏移，从而在采集到的各层信息载体图像中产生畸变现象(如：信息载体图像中某一区域会朝着图像的某一侧延伸或收缩等)，故可以根据畸变的方向对发生畸变的信息载体图像运用相关技术获取该发生畸变的各层信息载体图像的面内U、V方向的位移矢量场数据，进一步得到各层信息载体图像全场位移矢量场。
可选的，采用数字云纹技术、数字散斑技术或几何相位分析技术对对上层信息载体畸变图像和下层信息载体畸变图像进行处理得到上层信息载体图像和下层信息载体图像各个时刻的面内U、V方向的位移矢量数据。
数字云纹技术：所有云纹测量技术的基本元件都是栅线。两组栅线叠合在一起产生的云纹条纹会随着栅线几何形状的微小改变而发生宏观移动，云纹法就是通过测量条纹的移动量来分析物体的位移场。栅线是物体变形信息的载体，因此选用合适的栅线可以方便快捷地提取出物体的变形信息。云纹法中的参考栅线，可以是实际的几何栅线，或者是计算机生成的数字栅，栅线的类型通常是根据所要分析问题的类型和分析的方便性来选取。数字云纹技术采用的是数字栅叠加，其云纹条纹的生成不同于传统云纹技术中由两组栅线的光学叠加来实现，而是将数字化的变形栅线图进行二值化后与计算机生成的二值参考栅进行逻辑运算，来模拟栅线叠加形成的几何云纹。
数字散斑技术：数字图像相关法(Digital image correlation，DIC)，又称数字散斑相关法，是将变形前后的两幅数字图像，通过相关计算获取感兴趣区域的变形信息。其基本原理是，对变形前图像中的感兴趣区域进行网格划分，将每个子区域当作刚性运动。再针对每个子区域，通过一定的搜索方法按预先定义的相关函数来进行相关计算，在变形后图像中寻找与该子区域的互相关系数为最大值的区域，即该子区域在变形后的位置，进而获得该子区域的位移。对全部子区域进行计算，即可获得全场的变形信息。
几何相位分析技术：(Geometric phase analysis，GPA)其方法是对变形前后的点阵图像进行傅里叶变换，选择两个非线性矢量，再进行反傅里叶变换，可以得到变形前后点阵图像的相位信息，用变形后的相位信息减去变形前的相位信息得到相位差，根据相位差与应变及位移的关系获得应变场与位移场。
上层信息载体图像各个时刻的全场位移矢量场，包括上层信息载体图像各个时刻的面内U、V方向位移矢量场的矢量和，下层信息载体图像各个时刻的全场位移矢量场包括下层信息载体图像各个时刻的面内U、V方向位移矢量场的矢量和。
步骤106：按照以下方法，根据上层信息载体图像的全场位移矢量场与下层信息载体图像的全场位移矢量场求解液面各个时刻的液体高度数据：
其中，S1为上层信息载体图像的全场面内位移，S2为下层信息载体图像的全场面内位移，d为介质的上层与下层之间的厚度，α为图像采集设备所采集图像对应位置的光线与竖直方向的夹角，n1为液体的折射率，n2为介质的折射率。
可以理解的是，由于介质、液体和空气的折射率均不同，所以图像采集设备对介质上层和下层的信息载体图形进行图像采集时会发生不同程度的折射，由于介质的上层与下层之间具有厚度，所以上层信息载体图像的全场面内位移S1与下层信息载体图像的全场面内位移S2不相等。
进一步地，根据全场位移矢量场求解得到液面的高度分布，需要建立液面的高度于各层信息载体图像面内位移之间的函数关系。具体地说，如图2所示，液面的变形前后的剖面光路图，某点液面到上层信息载体图形之间的距离H满足上述公式，当然计算出高度H就能够了解了液面上各点的液面变形与形貌了。
实施例二：
结合图11，本实施例提供了一种基于双层图形透射成像的液面变形和形貌测量系统，包括图像采集设备201、介质205、面内位移矢量处理模块202以及液面三维形貌处理模块203。
介质205具有一定厚度，设置在液面底部，在垂直于地面所在平面的方向上，介质205包括相对设置的上层和下层，上层和下层均具有信息载体图形；
图像采集设备201，设置在液面远离地面的一侧，且图像采集设备201与液面具有间隔，并对图像采集设备201进行标定，与面内位移矢量处理模块202相耦接，用于液面静止时利用图像获取方法对信息载体图形进行图像采集，得到液面静止时的信息载体图像，图像获取方法包括：图像采集设备采集双层信息载体图像，上层信息载体图形和下层信息载体图形的图像采集到同一图像中，利用几何相位分析技术或者彩色阈值提取法将上层信息载体图形的图像和下层信息载体图形的图像分开到两张图像中，形成上层信息载体参考图像与下层信息载体参考图像；当液面发生形变后，对信息载体图形进行图像采集，得到液面形变后的信息载体畸变图像，信息载体畸变图像包括上层信息载体畸变图像和下层信息载体畸变图像，并发送至面内位移矢量处理模块202，采集一张信息载体畸变图像时，最终测量的是采集时刻的液面瞬时形貌，当采集随时间变化的序列信息载体畸变图像时，最终测量得到的是液面随时间变化的形貌；
面内位移矢量处理模块202，分别与图像采集设备201和液面三维形貌处理模块203相耦接，用于以液面静止时获取的上层信息载体参考图像与下层信息载体参考图像作为参考，对所述上层信息载体畸变图像和所述下层信息载体畸变图像进行处理得到所述上层信息载体图像和所述下层信息载体图像各个时刻的面内U、V方向的位移矢量数据，求解上层信息载体图像各个时刻的全场位移矢量场与下层信息载体图像各个时刻的全场位移矢量场，并发送至液面三维形貌处理模块203；
液面三维形貌处理模块203，与面内位移矢量处理模块202相耦接，用于接收面内位移矢量处理模块202发送的上层信息载体图像的全场位移矢量场与下层信息载体图像的全场位移矢量场数据，按照以下方法求解出液面高度:
其中，S1为上层信息载体图像的全场面内位移，S2为下层信息载体图像的全场面内位移，d为介质205的上层与下层之间的厚度，α为图像采集设备201所采集图像对应位置的光线与竖直方向的夹角，n1为液体204的折射率，n2为介质205的折射率。
可选的，信息载体图形包括条纹、散斑和点阵信息。
可选的，面内位移矢量处理模块202采用数字云纹技术、数字散斑技术或几何相位分析技术对对上层信息载体畸变图像和下层信息载体畸变图像进行处理得到上层信息载体图像和下层信息载体图像各个时刻的面内U、V方向的位移矢量数据。
可选的，介质205包括透明玻璃、或透明塑料，介质205的折射率与液体204的折射率不相等。
可选的，上层信息载体图像各个时刻的全场位移矢量场，包括上层信息载体图像各个时刻的面内U、V方向位移矢量场的矢量和，下层信息载体图像各个时刻的全场位移矢量场包括下层信息载体图像各个时刻的面内U、V方向位移矢量场的矢量和。
在本实施例中，各层信息载体图像是水底各层信息载体图形将有效光源反射出液面被图像采集设备采集到的图像，其中，信息载体图形是包括条纹，散斑，点阵等具备特定灰度梯度的图案。信息载体图像的采集可以通过CCD相机和/或将图像转变为数字信号功能的图像采集设备实现。在实际应用中，将双层信息载体图形置于盛有适量透明液体的玻璃水槽的外面底部。图像采集设备(如CCD相机)置于水槽的正上方。
进一步地，图像采集模块201可以是CCD相机和/或将图像转变为数字信号功能的图像采集设备。
可选的，面内位移矢量处理模块202可以为：由CPU、GPU和/或具有数据处理功能的芯片构成的面内位移矢量处理模块。
可选的，液面三维形貌处理203可以为：由CPU、GPU和/或具有数据处理功能的芯片构成的液面三维形貌处理模块。
在本实施例中，各层信息载体图像是水底各层信息载体图形将有效光源反射出液面被图像采集设备采集到的图像，其中，信息载体图形是包括条纹，散斑，点阵等具备特定灰度梯度的图案。信息载体图像的采集可以通过CCD相机和/或将图像转变为数字信号功能的图像采集设备实现。在实际应用中，将双层信息载体图形置于盛有适量透明液体的玻璃水槽的外面底部。图像采集设备(如CCD相机)置于水槽的正上方。
进一步地，双层信息载体图像采集模块201可以是CCD相机和/或将图像转变为数字信号功能的图像采集设备。
通过上述实施例可知，本发明提供的基于双层图形透射成像的液面变形和形貌测量方法和系统，至少实现了如下的有益效果：
1)本发明通过记录液面静止时和液面变形后各时刻的双层信息载体图像，运用相关技术，可以实现对静态和动态液面形貌进行快速高精度地测量；
2)本发明通过分析上下两层信息载体图像畸变与液面高度之间的光学关系，有效提高测量量程的同时提高了测量的灵敏度；
3)本发明测量上下两层信息载体图形反射的光线透射出变形液面后的畸变，结合该上下两层信息载体图形间给定的间隔距离，就可以综合得到液面形貌变化图，测量方法简单，计算效率高，并且便于实验操作。
4)本发明只需要使用单个的图像采集设备，装置简单经济。
虽然已经通过例子对本发明的一些特定实施例进行了详细说明，但是本领域的技术人员应该理解，以上例子仅是为了进行说明，而不是为了限制本发明的范围。本领域的技术人员应该理解，可在不脱离本发明的范围和精神的情况下，对以上实施例进行修改。本发明的范围由所附权利要求来限定。
Claims (10)
1.一种基于双层图形透射成像的液面变形和形貌测量方法，其特征在于，包括步骤：
提供图像采集设备，将图像采集设备设置在液面远离地面的一侧，且所述图像采集设备与所述液面具有间隔，并对所述图像采集设备进行标定；
提供具有一定厚度的介质，并将所述介质置于所述液面底部，在垂直于地面所在平面的方向上，所述介质包括相对设置的上层和下层，所述上层和下层均具有信息载体图形；
所述液面静止时，利用图像获取方法通过所述图像采集设备对信息载体图形进行图像采集，得到液面静止时的信息载体图像，所述图像获取方法包括所述图像采集设备采集双层信息载体图像，上层信息载体图形和下层信息载体图形的图像采集到同一图像中，利用几何相位分析技术或者彩色阈值提取法将上层信息载体图形的图像和下层信息载体图形的图像分开到两张图像中，形成上层信息载体参考图像与下层信息载体参考图像；
当所述液面发生形变时，利用所述图像获取方法对信息载体图形进行图像采集，得到液面形变后的信息载体畸变图像，所述信息载体畸变图像包括上层信息载体畸变图像和下层信息载体畸变图像，采集一张信息载体畸变图像时，最终测量的是采集时刻的液面瞬时形貌，当采集随时间变化的序列信息载体畸变图像时，最终测量得到的是液面随时间变化的形貌；
以液面静止时获取的上层信息载体参考图像与下层信息载体参考图像作为参考，对所述上层信息载体畸变图像和所述下层信息载体畸变图像进行处理得到所述上层信息载体图像和所述下层信息载体图像各个时刻的面内U、V方向的位移矢量数据，求解上层信息载体图像各个时刻的全场位移矢量场与下层信息载体图像各个时刻的全场位移矢量场；
按照以下方法，根据上层信息载体图像的全场位移矢量场与下层信息载体图像的全场位移矢量场求解所述液面各个时刻的液体高度数据：
其中，S1为上层信息载体图像的全场面内位移，S2为下层信息载体图像的全场面内位移，d为所述介质的上层与下层之间的厚度，α为图像采集设备所采集图像对应位置的光线与竖直方向的夹角，n1为液体的折射率，n2为所述介质的折射率。
2.根据权利要求1所述的基于双层图形透射成像的液面变形和形貌测量方法，其特征在于，所述信息载体图形包括条纹、散斑和点阵信息。
3.根据权利要求1所述的基于双层图形透射成像的液面变形和形貌测量方法，其特征在于，采用数字云纹技术、数字散斑技术或几何相位分析技术对所述对上层信息载体畸变图像和所述下层信息载体畸变图像进行处理得到所述上层信息载体图像和所述下层信息载体图像各个时刻的面内U、V方向的位移矢量数据。
4.根据权利要求1所述的基于双层图形透射成像的液面变形和形貌测量方法，其特征在于，所述介质包括透明玻璃、或透明塑料，所述介质的折射率与所述液体的折射率相等或不相等。
5.根据权利要求1所述的基于双层图形透射成像的液面变形和形貌测量方法，其特征在于，所述上层信息载体图像各个时刻的全场位移矢量场，包括所述上层信息载体图像各个时刻的面内U、V方向位移矢量场的矢量和，所述下层信息载体图像各个时刻的全场位移矢量场包括下层信息载体图像各个时刻的面内U、V方向位移矢量场的矢量和。
6.一种基于双层图形透射成像的液面变形和形貌测量系统，其特征在于，包括图像采集设备、介质、面内位移矢量处理模块以及液面三维形貌处理模块，其中，
所述介质具有一定厚度，设置在液面底部，在垂直于地面所在平面的方向上，所述介质包括相对设置的上层和下层，所述上层和下层均具有信息载体图形；
所述图像采集设备，设置在液面远离地面的一侧，且所述图像采集设备与所述液面具有间隔，并对所述图像采集设备进行标定，与所述面内位移矢量处理模块相耦接，用于所述液面静止时利用图像获取方法对信息载体图形进行图像采集，得到液面静止时的信息载体图像，所述图像获取方法包括：所述图像采集设备采集双层信息载体图像，上层信息载体图形和下层信息载体图形的图像采集到同一图像中，利用几何相位分析技术或者彩色阈值提取法将上层信息载体图形的图像和下层信息载体图形的图像分开到两张图像中，形成上层信息载体参考图像与下层信息载体参考图像；当所述液面发生形变后，对信息载体图形进行图像采集，得到液面形变后的信息载体畸变图像，所述信息载体畸变图像包括上层信息载体畸变图像和下层信息载体畸变图像，并发送至所述面内位移矢量处理模块，采集一张信息载体畸变图像时，最终测量的是采集时刻的液面瞬时形貌，当采集随时间变化的序列信息载体畸变图像时，最终测量得到的是液面随时间变化的形貌；
所述面内位移矢量处理模块，分别与所述图像采集设备和所述液面三维形貌处理模块相耦接，用于以液面静止时获取的上层信息载体参考图像与下层信息载体参考图像作为参考，对所述上层信息载体畸变图像和所述下层信息载体畸变图像进行处理得到所述上层信息载体图像和所述下层信息载体图像各个时刻的面内U、V方向的位移矢量数据，求解上层信息载体图像各个时刻的全场位移矢量场与下层信息载体图像各个时刻的全场位移矢量场，并发送至所述液面三维形貌处理模块；
所述液面三维形貌处理模块，与所述面内位移矢量处理模块相耦接，用于接收所述面内位移矢量处理模块发送的上层信息载体图像的全场位移矢量场与下层信息载体图像的全场位移矢量场数据，按照以下方法求解出液面高度:
其中，S1为上层信息载体图像的全场面内位移，S2为下层信息载体图像的全场面内位移，d为所述介质的上层与下层之间的厚度，α为图像采集设备所采集图像对应位置的光线与竖直方向的夹角，n1为液体的折射率，n2为所述介质的折射率。
7.根据权利要求6所述的基于双层图形透射成像的液面变形和形貌测量系统，其特征在于，所述信息载体图形包括条纹、散斑和点阵信息。
8.根据权利要求6所述的基于双层图形透射成像的液面变形和形貌测量系统，其特征在于，所述面内位移矢量处理模块采用数字云纹技术、数字散斑技术或几何相位分析技术对所述对上层信息载体畸变图像和所述下层信息载体畸变图像进行处理得到所述上层信息载体图像和所述下层信息载体图像各个时刻的面内U、V方向的位移矢量数据。
9.根据权利要求6所述的基于双层图形透射成像的液面变形和形貌测量系统，其特征在于，所述介质包括透明玻璃、或透明塑料，所述介质的折射率与所述液体的折射率相等或不相等。
10.根据权利要求6所述的基于双层图形透射成像的液面变形和形貌测量系统，其特征在于，所述上层信息载体图像各个时刻的全场位移矢量场，包括所述上层信息载体图像各个时刻的面内U、V方向位移矢量场的矢量和，所述下层信息载体图像各个时刻的全场位移矢量场包括下层信息载体图像各个时刻的面内U、V方向位移矢量场的矢量和。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110691868.0A CN113405628B (zh) | 2021-06-22 | 2021-06-22 | 基于双层图形透射成像的液面变形和形貌测量方法和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110691868.0A CN113405628B (zh) | 2021-06-22 | 2021-06-22 | 基于双层图形透射成像的液面变形和形貌测量方法和系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113405628A true CN113405628A (zh) | 2021-09-17 |
| CN113405628B CN113405628B (zh) | 2022-03-29 |
Family
ID=77682384
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110691868.0A Expired - Fee Related CN113405628B (zh) | 2021-06-22 | 2021-06-22 | 基于双层图形透射成像的液面变形和形貌测量方法和系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113405628B (zh) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3075814B2 (ja) * | 1991-11-20 | 2000-08-14 | マツダ株式会社 | 消失性模型を用いた鋳造方法及びその装置 |
| US6236324B1 (en) * | 2000-03-31 | 2001-05-22 | Lawrence V. Nash | Elastic plate liquid level sensor |
| CN101793542A (zh) * | 2010-03-16 | 2010-08-04 | 合肥工业大学 | 高灵敏度光纤光栅液位/流体压力传感器 |
| CN103453850A (zh) * | 2013-08-15 | 2013-12-18 | 北京理工大学 | 基于数字散斑相关技术的透明液面微形貌测量方法和系统 |
| CN203432649U (zh) * | 2013-08-22 | 2014-02-12 | 毛绍方 | 一种液位计 |
| CN104154955A (zh) * | 2014-05-19 | 2014-11-19 | 北京理工大学 | 贮箱液体推进剂液面形貌与剂量动态测量方法和系统 |
| CN104614046A (zh) * | 2015-02-11 | 2015-05-13 | 天津理工大学 | 基于ccd的油水界面位置测量方法 |
| CN109750749A (zh) * | 2019-01-04 | 2019-05-14 | 程志勇 | 一种智能化移动建筑 |
| CN110645917A (zh) * | 2019-09-24 | 2020-01-03 | 东南大学 | 基于阵列式相机的高空间分辨率三维数字图像测量方法 |
-
2021
- 2021-06-22 CN CN202110691868.0A patent/CN113405628B/zh not_active Expired - Fee Related
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3075814B2 (ja) * | 1991-11-20 | 2000-08-14 | マツダ株式会社 | 消失性模型を用いた鋳造方法及びその装置 |
| US6236324B1 (en) * | 2000-03-31 | 2001-05-22 | Lawrence V. Nash | Elastic plate liquid level sensor |
| CN101793542A (zh) * | 2010-03-16 | 2010-08-04 | 合肥工业大学 | 高灵敏度光纤光栅液位/流体压力传感器 |
| CN103453850A (zh) * | 2013-08-15 | 2013-12-18 | 北京理工大学 | 基于数字散斑相关技术的透明液面微形貌测量方法和系统 |
| CN203432649U (zh) * | 2013-08-22 | 2014-02-12 | 毛绍方 | 一种液位计 |
| CN104154955A (zh) * | 2014-05-19 | 2014-11-19 | 北京理工大学 | 贮箱液体推进剂液面形貌与剂量动态测量方法和系统 |
| CN104614046A (zh) * | 2015-02-11 | 2015-05-13 | 天津理工大学 | 基于ccd的油水界面位置测量方法 |
| CN109750749A (zh) * | 2019-01-04 | 2019-05-14 | 程志勇 | 一种智能化移动建筑 |
| CN110645917A (zh) * | 2019-09-24 | 2020-01-03 | 东南大学 | 基于阵列式相机的高空间分辨率三维数字图像测量方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113405628B (zh) | 2022-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cobelli et al. | Global measurement of water waves by Fourier transform profilometry | |
| CN103575227B (zh) | 一种基于数字散斑的视觉引伸计实现方法 | |
| Jurjo et al. | Experimental methodology for the dynamic analysis of slender structures based on digital image processing techniques | |
| CN104457603B (zh) | 一种高温环境下物体变形测量方法 | |
| Qiao et al. | A single-shot phase retrieval method for phase measuring deflectometry based on deep learning | |
| CN104279980A (zh) | 一种基于智能拍照手机的镜面三维面形测量系统 | |
| CN101655357A (zh) | 一种用于二维相位展开的相位梯度相关质量图获取方法 | |
| CN105043301A (zh) | 一种用于三维测量的光栅条纹相位求解方法 | |
| US8306760B1 (en) | Defect detection on optical fiber specimen using 3D surface data | |
| Aparna et al. | Fatigue testing of continuous GFRP composites using digital image correlation (DIC) technique a review | |
| Fu et al. | Non-contact optical dynamic measurements at different ranges: a review | |
| CN115824092B (zh) | 基于彩色复合条纹的相位测量偏折术缺陷检测装置及方法 | |
| TWI407078B (zh) | Micro - lens array surface profile detection system and its detection method | |
| Meng et al. | Defocused calibration for large field-of-view binocular cameras | |
| Liu et al. | A novel orthogonal transmission-virtual grating method and its applications in measuring micro 3-D shape of deformed liquid surface | |
| CN113405628B (zh) | 基于双层图形透射成像的液面变形和形貌测量方法和系统 | |
| CN112697065B (zh) | 基于相机阵列的三维形貌重构方法 | |
| Quan et al. | Determination of surface contour by temporal analysis of shadow moiré fringes | |
| Wang | Development and application of computer-aided fringe analysis | |
| Molimard et al. | Frequency-based image analysis of random patterns: an alternative way to classical stereocorrelation | |
| Zhong et al. | 3D shape measurement of complex transparent liquid surfaces using monoscopic deformed fringe transmissometry | |
| Wang et al. | Multi-lens stereo reconstruction of wave surface in a wave basin | |
| CN105066905B (zh) | 小波变换轮廓术抑噪方法 | |
| CN107543502A (zh) | 即时检测全场厚度的光学装置 | |
| CN111912356A (zh) | 一种基于l-k局域光流的物体三维面形测量方法和装置 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220329 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |