Summary
气穴微泡使用连接到变焦镜头的高速摄像机进行成像。阐述了实验设置，并采用图像分析计算气穴面积。使用 ImageJ 完成图像分析。
Abstract
提出了一种用于成像气穴气泡并计算气穴面积的实验和图像分析技术。这里提出的高速成像实验技术和图像分析方案，也可以应用于其他研究领域的成像显微气泡;因此，它具有广泛的应用。我们将其应用于牙科超声波刻度器周围的图像气穴。图像气穴来描述它，并了解如何将其用于各种应用，这一点非常重要。围绕牙科超声波刻度器发生的气穴可以用作一种新的牙菌斑去除方法，比目前的牙周治疗技术更有效，造成的损害更少。我们提出了一种使用高速摄像机和变焦镜头成像牙科超声波刻度器尖端周围发生的气穴气泡云的方法。我们还使用机器学习图像分析计算气穴面积。开源软件用于图像分析。提供的图像分析易于复制，不需要编程经验，并且可以很容易地修改以适应用户的应用。
Introduction
成像气泡的运动对于各种应用都很重要，因为它控制着系统的流体动力学。有许多应用，这非常有用：在流化床反应器1，2，,2或用于清洁气穴气泡3,3，4。成像气泡的目的是更好地了解气泡动力学或气泡云的方向和运动。这可以通过观察图像结构以及使用图像分析来获取定量信息（如气泡的大小）来实现。
气穴气泡是当压力低于饱和压力值5时在流体中发生的气体或蒸汽实体。当声场以超声波频率应用于流体时，它们可能发生。它们反复生长和坍塌，在坍塌时可以释放能量，以高速微喷气机和冲击波6，7,的形式。这些可以通过剪切力将表面上的颗粒清除，并导致表面清洁8。气穴气泡正在调查不同行业的表面清洁，如半导体，食品和伤口清洁,9，10，11，12。1011,129它们还可用于清洁牙齿和生物材料中的牙菌斑，如牙科植入物12、13。,13气穴发生在目前使用的牙科仪器，如超声波刻度仪和内音文件，并显示出潜在的附加清洁过程与这些仪器14。
气穴气泡的振荡发生在几微秒内，因此需要高速摄像机以每秒数千帧的成像来捕捉其运动。我们演示了一种围绕牙科超声波刻度器进行微泡气穴成像的方法。目的是了解气穴如何围绕不同的超声波刻度器变化，因此它可以被优化为清洁牙菌斑的新方法。
以前用于研究气穴的方法包括气精，它使用发光来检测气穴发生的地方15，16。15,然而，这是一种间接技术，它无法实时可视化气穴气泡。因此，它无法准确确定它在仪器上发生的确切地点，并且除非与其他成像技术相结合，否则无法获得有关气泡动力学的信息。高速成像不仅可以成像气穴气泡的产生和坍塌，还可以成像气穴发生的类型：气穴云、微流和微喷头6,6、7、18。,18这些提供了有关气穴如何清洁表面的更多信息。
我们提出了一种使用高速摄像机成像气穴微泡并计算气穴发生平均面积的方法。此方法使用围绕不同牙科超声波刻度器尖端发生的气穴示例进行演示，尽管实验和图像分析步骤可用于其他应用，例如成像其他宏和微泡。
Protocol
1. 仪器设置
- 选择要成像的仪器或对象。在这个实验中，对超声波刻度仪进行了成像。气穴气泡发生在水中超声波刻度器的尖端周围。
- 选择一个微型定位阶段，通过 XYZ 平移和旋转对仪器进行成像。放在实验室插孔上。将仪器手柄连接到微定位级
- 选择光学透明的水容器进行成像。这些实验中使用的容器是用玻璃显微镜幻灯片创建的。
- 选择具有旋转平台的 XY 阶段。放在实验室插孔上。将水容器放在舞台上，并加注过滤水（反渗透或蒸馏）。
2. 高速摄像机设置
- 选择具有所需帧速率和分辨率的高速摄像机以及带光纤导光的高强度光源。
- 将微定位滑动板连接到高速摄像机机身，并将其连接到三脚架支架上。
- 选择具有所需分辨率和焦距的镜头并将其连接到相机。对于此实验，以 8.4 μm/像素的分辨率使用变焦镜头。
- 向成像罐加注水，并定位仪器的尖端以所需方向在水箱中成像。
- 在连接相机并在软件中加载实时视图后，使用低放大倍率对焦到超声波刻度器的尖端，必要时重新定位光源。将仪器和光源放在摄像机前面并对焦。调整到所需的帧速率和亮度。
注：在高帧速率、短快门速度和/或高放大倍数下进行成像需要更高的光照强度。可以在反射模式或传输模式下提供照明。在此协议中，使用高强度冷照明装置在传输模式（亮场）中提供照明。 - 为高速摄像机设置最佳帧速率和快门速度。在这个实验中，帧速率为6400 fps，快门速度为262纳秒。快速移动气泡（如气穴气泡）需要短快门速度，以确保它们聚焦。
- 调整变焦镜头的放大倍率和光源的强度，使背景为白色，而不会过度曝光。
3. 校准
- 记录尖端的位置（在 x-y 级旋转，仪器的旋转角度，以进行可重复性）。
- 为确保每次重复的视场一致，请选择参考点并记下坐标。在这种情况下，参考点是超声波刻度器的尖端。然后，它可以在视图领域内同一位置的未来实验中重新定位。
- 如果像素大小未知，则在设定的放大倍率下使用具有 10 μm 标记的颗粒图像，并使用图像分析软件（如斐济）计算分辨率。
4. 高速视频录制
- 在不气穴的情况下对仪器进行成像。在计算气穴气泡区域时，将从图像分析中的气穴图像中减去。以 TIFF 等格式保存视频，这样不会丢失图像质量。
- 成像使用气穴操作的仪器。确保有足够的帧进行精确分析，例如，每次 5 次重复，每次 500 帧。
5. 图像处理
- 从 ImageJ 网站（https://imagej.net/Fiji）下载斐济19。 提供了 ImageJ 宏代码，该代码会自动执行下面描述的图像分析步骤，还可以更改以适应应用程序。宏的单个步骤在步骤 5.3-5.5 中描述。
- 裁剪图像以移除因照明不均匀而导致的任何较暗区域（如果需要）。确保所有图像都裁剪到图像中的相同大小和相同点。
- 通过使用自动阈值之一自动设置阈值，将图像转换为二进制图像。在此示例中，使用最小自动阈值。
- 运行填充孔命令，从气泡中删除任何被错误地分段的黑色像素。
- 计算堆栈的直方图，以显示与缩放器对应的像素数和每个帧中的气穴。
- 在这种情况下，与气泡对应的像素为白色，并且值为 255。保存这些测量值。
- 对于仪器在没有气泡的情况下运行的视频，重复步骤 5.3-5.6。
- 仅根据直方图的结果计算超声波刻度器尖端的均值。
- 从缩放器周围的气泡视频计算的每个区域中减去仪器的均值区域。气泡区域留有测量。
- 使用斐济的图像计算器从缩放器的二进制图像中减去具有气泡的缩放器的二进制图像来可视化。
- 计算气泡区域的平均值和标准偏差。
- 将值从像素数转换为面积（本例中为 μm2），乘以像素大小平方。通过以与成像相同的放大倍率对高速摄像机进行成像，并使用 ImageJ 设置比例来计算每个像素的大小。
- 绘制数据。如果比较不同的条件，也可以进行统计分析，以显示气泡面积的任何显著差异。
6. ImageJ 宏
- 在 ImageJ/Fiji 菜单中，转到 插件 > 新 > 宏。确保在语言菜单下选中 IJ1 宏，然后复制并粘贴以下代码。单击运行以执行宏 （补充文件）。
Representative Results
图1中显示了测试的超声波刻度器尖端之一的图像分析步骤。FSI 1000 尖端和 10P 尖端在水箱内成像，冷却水关闭（图 2）。气穴发生在尖端 FSI 1000 的弯曲附近，最大功率，靠近尖端 10P 中的自由端（图 3 和图 4）。FSI 1000 尖端的气穴平均面积为 0.1 ± 0.07 mm 2，10P 尖端的均值为 0.50 ±0.25 mm 2（图 5）。2
图1：高速成像设置和图像分析 步骤 （a） 研究中使用的高速成像设置的示意图。（b） 研究中使用的图像分析步骤的示意图，仅显示刻度器尖端左侧的原始图像，并带气穴，然后相互二进制和减去，以计算气穴云的面积。 请单击此处查看此图的较大版本。
图2：比较不同提示高速 图像静止显示气穴发生在测试的两个超声波刻度器尖端周围 （a） FSI 1000 （b） 10P. 请点击这里查看这个数字的较大版本。
图 3：提示 10P 高速图像：高速图像静止提示 10P，以每秒 6400 帧拍摄的视频。 在尖端的自由端周围可以看到气穴。 请单击此处查看此图的较大版本。
图 4：提示 FSI1000 高速图像：提示 FSI 1000 的高速图像静止图像，以每秒 6400 帧的速度拍摄的视频。 在尖端中间可以看到气穴。 请单击此处查看此图的较大版本。
图5：气穴区图像分析结果。 FSI 1000 和 10P 超声波刻度仪尖端周围气穴平均面积，使用所述图像分析技术计算。误差条表示标准偏差。 请单击此处查看此图的较大版本。
补充文件。 请点击这里下载此文件。
Discussion
本文所述的技术能够对具有高空间和时间分辨率的快速移动微泡进行成像。它有可能有利于广泛的科学学科，如化学工程，牙科和医学。工程应用包括用于清洁表面的成像气穴气泡，或用于流化床反应器中的成像气泡。生物医学应用包括围绕医疗和牙科仪器的成像气穴，以及使用气穴气泡从硬和软组织中成像生物膜脱毛。在这项研究中，我们演示了该技术，通过两个不同的牙科超声波刻度器尖端成像气穴。本研究中测试的两个尖端之间，气穴量不同，在尖端 10P 的自由端周围观察到更多的气穴云。这以前与振动振幅20有关。高速视频显示，FSI 1000 尖端的振动较小，这可能是为什么围绕此尖端的气穴较少的原因。
图像分析方法的一个局限性是，删除刻度器面积的图像减法技术并不完全准确，因为刻度器是振荡的，因此减法可能会将刻度器的一些区域错误地分割为气泡。但是，这一点通过从大量帧（n=2000）中平均面积来解释。对于要减去的对象是静止的应用程序中，这不是问题。对于要减去的移动对象的算数具有更高方差的算物，我们建议在减去结果之前同步两个视频中的移动。在目前的研究中，我们没有同步振荡，但由于振动较低，我们可以假设在两次测量中振荡彼此非常对应。
图像阈值是准确的，因为亮场照明提供了具有良好对比度的统一背景。确保背景是统一的，并且不包含任何其他可能错误地分段的对象，这一点至关重要。可以使用其他自动阈值来修改阈值方法，以适应应用程序。手动阈值（用户设置阈值）也是可能的，但不建议这样做，因为它会降低结果的可重复性，因为不同的用户将选择不同的阈值。
图像分析已用于许多其他气泡成像研究。他们还使用类似的背光方法，在气泡和背景之间获得最佳对比度，并阈值对气泡进行分割 21、22、23、24。,22,23,24本研究显示的方法也可以推广到许多不同的气泡成像应用，这些应用不仅限于高速成像。高速成像已用于在水中产生的气穴气泡，也围绕仪器，如内音文件和超声波刻度器12，25，26，27，28。12,25,26,27,28例如，Rivas等人和Macedo等人使用连接到显微镜上的高速摄像机，由冷光源提供照明，用于用气穴进行图像清洁，并在内向文件17、29周围进行图像气穴。明亮的场照明在背景和气泡之间提供了更多的对比度，因此可以使用简单的分割技术，如阈值，正如Rivas等人所证明的，用于成像和量化气穴侵蚀和清洁随着时间的推移29。由于灰度4，30的变异性较高，暗场照明使阈,值更加困难。图像分析已用于其他研究，以收集有关气泡1，2,的信息。Vyas等人使用机器学习方法在超声波刻度器20周围分割气穴气泡。本文中描述的方法更快，因为它使用简单的阈值，因此计算密集度较低，并且可以分析在刻度器上方和下方发生的气泡。但是，当前论文中使用的阈值方法只有在背景均匀时才准确。如果在成像过程中无法获得统一的背景，可以使用其他图像处理技术，例如使用背景减法使用滚球半径来校正不均匀的照明，使用中值或高斯滤镜进行滤波以消除噪声，或者也使用基于机器学习的技术20，31。,31
最后，提出了一种高速成像和分析方案，用于成像和计算微观运动物体的面积。通过成像超声波刻度器周围的气穴气泡，我们演示了这种方法。它可用于其他牙科仪器（如牙内音文件）周围的成像气穴，并可轻松适应其他非牙科气泡成像应用。
Disclosures
作者没有什么可透露的。
Acknowledgments
作者感谢工程和物理科学研究理事会EP/P015743/1的资助。
Materials
| Name | Company | Catalog Number | Comments |
| 0.25x attachment | Navitar | 1-50011 | |
| 12x with 12mm fine focus Long distance microscope zoom lens |
Navitar | 1-50486 | |
| 2x adaptor with f mount | Navitar | 1-62922 | |
| Cavitron Plus Ultrasonic Scaler | Dentsply Sirona | 8184003 | |
| Cavitron Ultrasonic Insert FSI 1000FSI 1000 | Dentsply Sirona | UCAFTHD | |
| Fibre light guide. 8mm fibre bundle 1500mm length. Focussing lens assembly for Hayashi light, 1/4"-20 tripod thread for mounting. |
Hayashi | LGC1- 8L1500 |
|
| Geared head | Manfrotto | MN405 | 7.5kg load capacity |
| HDF7010 High-Power LED Endoscope light source. 150W LED provides cold output equivalent to 250W Xenon. |
Hayashi | LA-HDF710 | |
| Heavy weight Tripod | Manfrotto | MN475B | Geared centre column, 12kg load capacity |
| High Speed Camera | Photron | 103526 | FASTCAM Mini AX200 900K M3 (16GB memory) |
| High-Precision Rotation Stage | Thorlabs | PR01/M | |
| Laboratory jacks | Camlab | 1194083 | |
| Micropositioning sliding plate | Manfrotto | SKU 454 | |
| Micropositioning stage 3D | Thorlabs | PT3/M | |
| Micropositioning stage rotation | Thorlabs | OCT-XYR1/M | OCT-XYR1/M - XY Stage with Solid Top Plate |
| NEWTRON P5 XS Ultrasonic Scaler | Acteon | F62118 | |
| Ultrasonic Insert 10P | Acteon | F00253 |
References
- Asegehegn, T. W., Schreiber, M., Krautz, H. J. Investigation of bubble behavior in fluidized beds with and without immersed horizontal tubes using a digital image analysis technique. Journal of Power Technologies. 210 (3), 248-260 (2011).
- Busciglio, A., Vella, G., Micale, G., Rizzuti, L. Analysis of the bubbling behaviour of 2D gas solid fluidized beds: Part I. Digital image analysis technique. Chemical Engineering Journal. 140 (1), 398-413 (2008).
- Versluis, M.
- Matsumoto, H., Yoshimine, Y., Akamine, A. Visualization of irrigant flow and cavitation induced by Er: YAG laser within a root canal model. Journal of Endodontics. 37 (6), 839-843 (2011).
- Young, F. R. Cavitation. , World Scientific. (1999).
- Brennen, C. E. Cavitation and Bubble Dynamics. , Cambridge University Press. (2013).
- Leighton, T. The acoustic bubble. , Academic Press. (2012).
- Verhaagen, B., Rivas, D. F. Measuring cavitation and its cleaning effect. Ultrasonics Sonochemistry. 29, 619-628 (2016).
- Oulahal-Lagsir, N., Martial-Gros, A., Boistier, E., Blum, L., Bonneau, M. The development of an ultrasonic apparatus for the non-invasive and repeatable removal of fouling in food processing equipment. Letters in Applied Microbiology. 30 (1), 47-52 (2000).
- Gale, G. W., Busnaina, A. A. Roles of cavitation and acoustic streaming in megasonic cleaning. Particulate Science and Technology. 17 (3), 229-238 (1999).
- Erriu, M., et al. Microbial biofilm modulation by ultrasound: Current concepts and controversies. Ultrasonics Sonochemistry. 21, 15-22 (2014).
- Van der Sluis, L., Versluis, M., Wu, M., Wesselink, P. Passive ultrasonic irrigation of the root canal: a review of the literature. International Endodontic Journal. 40 (6), 415-426 (2007).
- Vyas, N., Sammons, R. L., Addison, O., Dehghani, H., Walmsley, A. D. A quantitative method to measure biofilm removal efficiency from complex biomaterial surfaces using SEM and image analysis. Scientific Reports. 6, 32694 (2016).
- Walmsley, A. D., Lea, S. C., Felver, B., King, D. C., Price, G. J. Mapping cavitation activity around dental ultrasonic tips. Clinical Oral Investigations. 17 (4), 1227-1234 (2013).
- Price, G. J., Tiong, T. J., King, D. C. Sonochemical characterisation of ultrasonic dental descalers. Ultrasonics Sonochemistry. 21, 2052-2060 (2014).
- Felver, B., King, D. C., Lea, S. C., Price, G. J., Damien Walmsley, A. Cavitation occurrence around ultrasonic dental scalers. Ultrasonics Sonochemistry. 16, 692-697 (2009).
- Macedo, R. G., et al. Sonochemical and high-speed optical characterization of cavitation generated by an ultrasonically oscillating dental file in root canal models. Ultrasonics Sonochemistry. 21, 324-335 (2014).
- Reuter, F., Lauterborn, S., Mettin, R., Lauterborn, W. Membrane cleaning with ultrasonically driven bubbles. Ultrasonics Sonochemistry. 37, 542-560 (2017).
- Schindelin, J., et al. Fiji: an open-source platform for biological-image analysis. Nature Methods. 9 (7), 676-682 (2012).
- Vyas, N., et al. High-speed Imaging of Cavitation around Dental Ultrasonic Scaler Tips. PLoS One. 11 (3), 0149804 (2016).
- Ahmed, F. S., Sensenich, B. A., Gheni, S. A., Znerdstrovic, D., Al Dahhan, M. H. Bubble dynamics in 2D bubble column: comparison between high-speed camera imaging analysis and 4-point optical probe. Chemical Engineering Communications. 202 (1), 85-95 (2015).
- Honkanen, M. Reconstruction of three-dimensional bubble surface from high-speed orthogonal imaging of dilute bubbly flow. Proceedings of Computational Methods in Multiphase Flow V, New Forest, UK. , 469-480 (2009).
- do Amaral, C. E., et al. Image processing techniques for high-speed videometry in horizontal two-phase slug flows. Flow Measurement and Instrumentation. 33, 257-264 (2013).
- Lau, Y., Deen, N., Kuipers, J. Development of an image measurement technique for size distribution in dense bubbly flows. Chemical Engineering Science. 94, 20-29 (2013).
- Matsumoto, Y., Yoshizawa, S. Behaviour of a bubble cluster in an ultrasound field. International Journal for Numerical Methods in Fluids. 47 (6-7), 591-601 (2005).
- Peeters, H. H., Iskandar, B., Suardita, K., Suharto, D. Visualization of removal of trapped air from the apical region of the straight root canal models generating 2-phase intermittent counter flow during ultrasonically activated irrigation. Journal of Endodontics. 40 (6), 857-861 (2014).
- Halford, A., et al. Synergistic effect of microbubble emulsion and sonic or ultrasonic agitation on endodontic biofilm in vitro. Journal of Endodontics. 38 (11), 1530-1534 (2012).
- Kauer, M., Belova-Magri, V., Cairós, C., Linka, G., Mettin, R. High-speed imaging of ultrasound driven cavitation bubbles in blind and through holes. Ultrasonics Sonochemistry. 48, 39-50 (2018).
- Rivas, D. F., et al. Localized removal of layers of metal, polymer, or biomaterial by ultrasound cavitation bubbles. Biomicrofluidics. 6 (3), 034114 (2012).
- Pishchalnikov, Y. A., et al. Cavitation Bubble Cluster Activity in the Breakage of Kidney Stones by Lithotripter Shockwaves. Journal of Endourology. 17 (7), 435-446 (2003).
- Sternberg, S. R.
