Image analysis and processing methods in verifying the correctness of performing low-invasive esthetic medical procedures
© Koprowski et al.; licensee BioMed Central Ltd. 2013
Received: 5 April 2013
Accepted: 5 June 2013
Published: 9 June 2013
Efficacy and safety of various treatments using fractional laser or radiofrequency depend, to a large extent, on precise movement of equipment head across the patient’s skin. In addition, they both depend on uniform distribution of emitted pulses throughout the treated skin area. The pulses should be closely adjacent but they should not overlap. Pulse overlapping results in amplification of irradiation dose and carries the danger of unwanted effects.
Images obtained in infrared mode (Flir SC5200 thermovision camera equipped with photon detector) were entered into Matlab environment. Thermal changes in the skin were forced by CO2RE laser. Proposed image analysis and processing methods enable automatic recognition of CO2RE laser sites of action, making possible to assess the correctness of performed cosmetic procedures.
80 images were acquired and analyzed. Regions of interest (ROI) for the entire treatment field were determined automatically. In accordance with the proposed algorithm, laser-irradiated L i areas (ROI) were determined for the treatment area. On this basis, error values were calculated and expressed as percentage of area not covered by any irradiation dose (δ o ) and as percentage area which received double dose (δ z ). The respective values for the analyzed images were δ o =17.87±10.5% and δ z =1.97±1.5%, respectively.
The presented method of verifying the correctness of performing low-invasive esthetic medical (cosmetic) procedures has proved itself numerous times in practice. Advantages of the method include: automatic determination of coverage error values δ o and δz, non-invasive, sterile and remote-controlled thermovisual mode of measurements, and possibility of assessing dynamics of patient’s skin temperature changes.
KeywordsImage processing Esthetic medical procedures Laser Fully automatic
Esthetic medicine market is among the most dynamically developing sectors of industry across the world. Especially popular prove treatments that are low-invasive. According to the American Society of Plastic Surgeons, in the USA alone 13.8 million of low-invasive procedures were performed in 2011, at the estimated value of ca. USD 12.2 billion. Considering trends in population demographics in both developed and developing countries, the tendency of low-invasive esthetic medical procedures to grow in popularity will likely become more dynamic. Owing to the new developments in the medical equipment market it has now become possible to obtain satisfactory results of esthetic medicine procedures, together with a relatively low treatment invasiveness which also means shorter recovery time. Such minimum invasiveness, modest intensity of adverse side effects, coupled with satisfactory effects of treatment is possible, however, only after optimizing procedure parameters.
Action of fractional lasers causes focal ablation of epidermis, whereas radiofrequency-based procedures result in local overheating of both epidermis and corium. At tissue level these agents cause remodeling of collagen fibers and stimulate epidermal regeneration [1–7].
Energy of laser radiation or radiofrequency energy is delivered to patient’s skin in the form of pulses which affect a definite tissue area. Efficacy and safety of treatment using a fractional laser or radiofrequency depend, to a great extent, on precise movement of the therapeutic equipment head across the patient’s skin . In addition, they both depend on uniform distribution of emitted pulses throughout the treated skin area. The pulses should be closely adjacent but they should not overlap [8, 9]. Pulse overlapping results in dose amplification and carries the danger of undesired effects  and .
Known methods for monitoring the efficiency and safety of laser esthetic procedures include optoacoustic [11, 12] and optodynamic methods . They are based on using high-sensitivity cameras (detectors) that allow imaging fluorescence phenomena in real time, as well as measuring fluorescence intensity while performing concomitant spectral analysis. Use of these methods has its limitations due to, for example, speed of image acquisition or impossibility of treating the patient in a totally uninvasive manner with increased positioning distances.
In the presented study we aimed to verify the correctness of performing laser-mediated esthetical medical procedures. This was achieved based on automatic calculation of the degree of coverage of the treated area by CO2RE laser-sent pulses. The study was performed using the proposed method of analyzing images obtained in infrared mode.
We analyzed thermovisual image sequences collected from 15 patients using a Flir SC5200 thermal imaging camera equipped with photon detector. In total, 80 images were analyzed (from patients’ right and left cheek, chin, forehead and nose). Thermal changes in human skin were induced using CO2RE laser. All patients were adequately prepared prior to the cosmetic procedure and thermovisual measurements. The error of thermovisual measurement method was minimized by taking into account 1) false estimation of the object’s emissivity, 2) radiation originating from the surroundings and reflected from the object, 3) atmospheric attenuation and scattering and own atmospheric emission, 4) changes in emission from camera optical components, 5) errors intrinsic to methodology of adopted measurement course, 6) air current convection, 7) emotional state of patient, 8) patient’s dress, 9) thermal conductivity of limited and diffuse heat sources, 10) skin vascularization, 11) meals eaten by patient within preceding 24 hrs, 12) crossed radiation, 13) patient’s movements prior to and during examination, 14) undisclosed diseases, and 15) faults in the algorithm. The cited measurement errors can be easily diminished or totally eliminated by assuring constant room temperature (no air movement due to drafts or air conditioning, solar irradiation, radiators, etc.) and securing time necessary to acclimatize a patient in such room (typically half an hour). This is an essential requirement for the majority of thermovisual measurememts.
After minimizing errors due to these causes we carried out measurements according to the methodology presented below.
Infrared images generated by a thermal scanning camera (Flir SC5200) equipped with photon detector were entered into Matlab environment. The camera has indium antimonide (InSb) detector, with 3–5 μm spectral range, 320×256 pixel resolution and 30×30 μm pixel pitch. Thermal changes in patient’s skin were induced by CO2RE laser (CO2 type – 10600 nm wave-length, pulsed laser beam emission mode, 4.5 J/cm2 energy density, 1–150 mJ impulse energy, 16.7 kHz impulse frequency, 20–3000 μs impulse duration, 10 mm 2 maximum scanned area, 120 μm or 150 μm dot size) [6, 14–16].
Application of individual irradiation doses and, thus, induction of thermal changes, was performed manually, by sequentially applying the equipment head to patient’s face. The irradiation procedure was performed by an expert physician who attempted to cover as much as feasible of the whole analyzed area . Synchronization of laser triggering and image acquisition had been programmed. Synchronization error due to operating system delay, data transmission timing and number of stills per second did not exceed 0.1 second. The thermal scanning equipment was positioned during the procedure at ca. 30 cm from patient’s face.
Input image L(m,n) (where m denotes rows and n denotes columns) at 320×256 pixels, following an increase in resolution to M×N=480×640 pixels, was filtered using a median filter with h mask dimensions M h ×N h =3×3. Increase of image resolution was achieved using the nearest neighbor method, avoiding thus new pixel values [17–19]. The generated image L MED was subjected to processing aimed at detecting regions of interest (ROI). Automatic ROI detection denotes a phase in the image analysis and processing, during which assignment of skin area receiving irradiation dose takes place. This method should work correctly in analyzing temperature changes in areas subjected to therapeutic procedure. Due to close relationship between forcing of temperature change and time of skin temperature reaction onset, this process was followed in more detail.
Binarization (L B ) may result also in other smaller areas. Removal of smaller, erroneously indicated areas was achieved by labeling procedure (marking of clusters). The largest cluster (area) was then chosen. In all of the acquired thermovisual images this area was indicated correctly. This area (ROI) and L R image created on its basis were the subject of subsequent analysis.
Interpretation of errors δ o and δ z defined in this way is straightforward. It defines correctness of the performed procedure (forcing). Increased value of δ o means that the operator did not drive the laser head uniformly, leaving untreated areas. This is not harmful to the patient but requires additional corrective treatment. On the other hand, increased values of δ z error indicate harm, as patient receives in these areas a double dose of irradiation.
Practical use of this algorithm is demonstrated in the next section.
Determination of δ o and δ z errors is affected also by other elements specific for the algorithm itself or/and for procedure methodology. These elements include: 1) error due to camera placement with respect to patient’s skin, 2) error due to non-perpendicular placement of laser head with respect to patient’s skin during the procedure of applying irradiation dose, 3) local disturbances in skin thermoregulation, 4) presence of perspiration, 5) interference such as, e. g., patient’s hair falling accidentally onto forehead during the procedure.
In practice, the last two elements predominate: one is caused by skin reaction to temperature and the other by improperly secured hair.
The obtained results as well as δ o and δ z error values are affected by personal habits of the operator (technician). It has been noticed that these error values strongly depend on individual habits of the technician and, only to a lesser degree, on the shape of facial area subjected to treatment. Differences in δ o and δ z error values for two technicians may vary in a broad range. Due to this reason an automated system of laser triggering has been proposed. The system is based on tracking in visible light the skin areas subjected to treatment (a CCD camera is placed in the laser head). Treated skin areas are memorized using visible light. Following manual relocation (in any direction) of the laser head by a fixed distance, the laser is automatically triggered. Such a system allows minimizing values of δ o and δ z errors. In addition, these errors stay independent of the individual habits of an operator. The system is patent-protected  and shall be described in detail in future papers.
Comparison of results with other methods
Contemporary generation of lasers do not offer yet a qualitative analysis of procedures performed using them. It is assumed that an expert laser operator performs the cosmetic procedure correctly, without causing overlapping of irradiation doses. Various practical methods allowing laser beam control and visible image analysis have been reported to date, especially in patent claim literature . Also known have been descriptions of visualization methods using visible light and accomplished by various types of cameras placed, e. g., in laser head . None of these solutions offers, however, analysis of the correctness of procedures performed with laser equipment. Neither temperature fields nor their degree of homogeneity have been assessed. In only a few reports temperature fields and their distribution within the skin were analyzed e. g. [29–31]. In the model proposed by Frahm et al. , model simulations of superficial temperature correlated with the measured skin surface temperature (ρ>0.90, p<0.001). Reported were studies comparing three Infrared Thermal Detection Systems. In this case correlations between ITDS and oral temperatures were similar for OptoTherm (ρ=0.43) and FLIR (ρ=0.42), but significantly lower for Wahl (ρ=0.14, p<0.001). Among numerous references pertaining to application of thermovision in medicine only a few e. g. [30–38] have dealt with quantitative (not qualitative) measurements. As an example, Bagavathiappan et al.  reported a temperature difference of 0.7–1°C as statistically significant. Based on this one can conclude that thermovisual analysis of human skin does require taking into account numerous factors which interfere with measurement. In the case of the algorithm presented herein only a minute skin fragment is analyzed. An expert laser operator has full control over this fragment, and is capable of minimizing the effect of additional factors to a negligible level.
The presented method of verifying the correctness of performing laser-mediated esthetic medical procedures has repeatedly proven itself in practice. Its advantages include: 1) automatic determination of δ o and δ z error values, 2) non-invasive sterile and remote-controlled thermovisual measurements, 3) possibility of learning how to assess procedure correctness through training, 4) assessment of dynamics of patient’s skin temperature changes, and 5) assessment of correct choice of irradiation dose, treatment length and individual equipment setting.
The described method has been currently used in esthetic medical procedures performed at the Silesian Medical College in Katowice, Poland.
No outside funding was received for this study.
- Adrian RM: Pulsed carbon dioxide and long pulse 10–ms erbium-YAG laser resurfacing: a comparative clinical and histological study. J Cutan Laser Ther 1999, 1: 197–202. 10.1080/14628839950516670View ArticleGoogle Scholar
- Fitzpatrick RE, Rostan EF, Marchell N: Collagen tightening by carbon dioxide versus erbium:YAG laser. Lasers Surg Med 2000, 27: 395–403. 10.1002/1096-9101(2000)27:5<395::AID-LSM1000>3.0.CO;2-4View ArticleGoogle Scholar
- Horton S, Alster TS: Preoperative and postoperative considerations for carbon dioxide laser resurfacing. Cutis 1999, 64: 399–406.Google Scholar
- Ross EV, Sajben FP, Hsia J, Barnette D, Miller CH, McKinlay JR: Non ablative skin remodeling: selective dermal heating with mid-infrared laser and contact cooling combination. Lasers Surg Med 2000, 26: 186–195. 10.1002/(SICI)1096-9101(2000)26:2<186::AID-LSM9>3.0.CO;2-IView ArticleGoogle Scholar
- Koch RJ: Radiofrequency non-ablative tissue tightening. Facial Plast Surg Clin North Am 2004, 12(3):339–346. 10.1016/j.fsc.2004.02.007View ArticleGoogle Scholar
- Abraham MT, Ross EV: Current concepts in nonablative radiofrequency rejuvenation of the lower face and neck. Facial Plast Surg 2005, 21(1):65–73. 10.1055/s-2005-871765View ArticleGoogle Scholar
- Carroll L, Humphreys TR: Laser tissue interactions. Clin Dermatol 2006, 24: 2–7. 10.1016/j.clindermatol.2005.10.019View ArticleGoogle Scholar
- Manuskiatti W, Triwongwaranat D, Varothai S, Eimpunth S, Wanitphakdeedecha R: Efficacy and safety of carbon-dioxide ablative fractional resurfacing device for treatment of atrophic acne scars in. J Am Acad Dermatol 2009, 63(2):274–283.View ArticleGoogle Scholar
- Steiner R: New laser technology and future applications. Med Laser Appl 2006, 21: 131–140. 10.1016/j.mla.2006.03.007View ArticleGoogle Scholar
- Alexiades-Armenakas MR, Dover JS, Arndt KA: The spectrum of laser skin resurfacing: Nonablative, fractional and ablative resurfacing. J Am Acad Dermatol 2008, 58(5):719–737. 10.1016/j.jaad.2008.01.003View ArticleGoogle Scholar
- Buehler A, Kacprowicz M, Taruttis A, Ntziachristos V: Real-time handheld multispectral optoacoustic imaging. Opt Lett 2013, 38(9):1404–1406. 10.1364/OL.38.001404View ArticleGoogle Scholar
- Ntziachristos V: Clinical translation of optical and optoacoustic imaging. Philos Trans A Math Phys Eng Sci 2011, 369(1955):4666–4678. 10.1098/rsta.2011.0270View ArticleGoogle Scholar
- Cencič B, Grad L, Možina J, Jezeršek M: Optodynamic monitoring of laser tattoo removal. J Biomed Opt 2012, 17(4):047003. 10.1117/1.JBO.17.4.047003View ArticleGoogle Scholar
- Foster KR, Zhang H, Osepchuk JM: Thermal response of tissues to millimeter waves: implications for setting exposure guidelines. Health Phys 2010, 99(6):806–810. 10.1097/HP.0b013e3181db29e6View ArticleGoogle Scholar
- Foster KR: Thermographic detection of breast cancer. IEEE Eng Med Biol Mag 1998, 17(6):10. 10.1109/51.734241View ArticleGoogle Scholar
- Foster KR, Morrissey JJ: Thermal aspects of exposure to radiofrequency energy: Report of a workshop. Int J Hyperthermia 2011, 27(4):307–319. 10.3109/02656736.2010.545965View ArticleGoogle Scholar
- Korzyńska A, Hoppe A, Strojny W, Wertheim D: Investigation of a combined texture and contour method for segmentation of light microscopy cell images. Proceedings of The Second IASTED International Conference on Biomedical Engineering 2004, 234–239.Google Scholar
- Gonzalez RC, Woods RE: Digital Image Processing Using Matlab. In Prentice Hall. Printed in United States of America; 2008.Google Scholar
- Koprowski R, Wojaczynska-Stanek K, Wrobel Z: Automatic segmentation of characteristic areas of the human head on thermographic images. Machine Graphics and Vision 2007, 16(3–4):251–274.Google Scholar
- Otsu N: A threshold selection method from gray-level histograms. IEEE Trans Sys Man Cyber 1979, 9(1):62–66.MathSciNetView ArticleGoogle Scholar
- Porwik P, Wróbel K, Doroz R: The Polish Coins Denomination Counting by Using Oriented Circular Hough Transform. Advances in Intelligent and Soft Computing 2009, 57: 569–576. 10.1007/978-3-540-93905-4_66View ArticleGoogle Scholar
- Porwik P, Para T: Some handwritten signature parameters in biometric recognition process. Proceedings of the International Conference on Information Technology Interfaces, (ITI2007), Dubrovnik 2007, 185–190.View ArticleGoogle Scholar
- Wróbel K, Doroz R: The new method of signature recognition based on least squares contour alignment. International Conference On Biometrics And Kansei Engineering 2009, 80–83.Google Scholar
- Sonka M, Michael Fitzpatrick J: Medical Image Processing and Analysis. In Handbook of Medical Imaging. Belligham: SPIE; 2000.Google Scholar
- Koprowski R, Wróbel Z: Image Processing in Optical Coherence Tomography Using Matlab. Katowice, Poland: University of Silesia; 2011.Google Scholar
- Koprowski R, Wróbel Z, Wilczyński S: A system to help performing low-invasive aesthetic medical procedures. Urząd Patentowy Rzeczypospolitej Polskiej 2012. Patent number P-398896 (submission date 20.04.2012)Google Scholar
- Wynne JJ, Gomory SH, Felsensteln JM: Laser dermablator and dermablation. United States Patent 2000. Patent Number 6,165,170, Date of Patent Dec. 26, 2000Google Scholar
- Altshuler GB, O’Shea L, Lazanicka OM: Dermatological treatment with visualization. United States Patent 2007. Patent Number 7,220,254 B2, Date of Patent May. 22, 2007Google Scholar
- Frahm KS, Andersen OK, Arendt-Nielsen L, Mørch CD: Spatial temperature distribution in human hairy and glabrous skin after infrared CO2 laser radiation. Biomed Eng Online 2010, 9: 69. 10.1186/1475-925X-9-69View ArticleGoogle Scholar
- Vogel A, Dlugos C, Nuffer R, Birngruber R: Optical properties of human sclera, and their consequences for transscleral laser applications. Lasers Surg Med 1991, 11(4):331–340. 10.1002/lsm.1900110404View ArticleGoogle Scholar
- Yeo C, Son T, Park J, Lee YH, Kwon K, Nelson JS, Jung B: Development of compression-controlled low-level laser probe system: towards clinical application. Lasers Med Sci 2010, 25(5):699–704. 10.1007/s10103-010-0779-8View ArticleGoogle Scholar
- Nguyen AV, Cohen NJ, Lipman H, Brown CM, Molinari NA, Jackson WL, Kirking H, Szymanowski P, Wilson TW, Salhi BA, Roberts RR, Stryker DW, Fishbein DB: Comparison of 3 Infrared Thermal Detection Systems and Self-Report for Mass Fever Screening. Emerg Infect Dis 2010, 16(11):1710–1717. 10.3201/eid1611.100703View ArticleGoogle Scholar
- Bagavathiappan S, Saravanan T, Philip J, Jayakumar T, Raj B, Karunanithi R, Panicker TMR, Korath MP, Jagadeesan K: Infrared thermal imaging for detection of peripheral vascular disorders. J Med Phys 2009, 34(1):43–47. 10.4103/0971-6203.48720View ArticleGoogle Scholar
- Herry CL, Frize M: Quantitative assessment of pain-related thermal dysfunction through clinical digital infrared thermal imaging. Biomed Eng Online 2004, 3: 19. 10.1186/1475-925X-3-19View ArticleGoogle Scholar
- Bichinho GL, Gariba MA, Sanches IJ, Gamba HR, Cruz FPF, Nohama PN: A Computer Tool for the Fusion and Visualization of Thermal and Magnetic Resonance Images. J Digit Imaging 2009, 22(5):527–534. 10.1007/s10278-007-9046-3View ArticleGoogle Scholar
- Jones BF: A reappraisal of the use of infrared thermal image analysis in medicine. IEEE Trans Med Imag 1998, 17(6):1019–1027. 10.1109/42.746635View ArticleGoogle Scholar
- Jones BF, Plassmann P: Digital infrared thermal imaging of human skin. IEEE Eng Med Biol Mag 2002, 21(6):41–48. 10.1109/MEMB.2002.1175137View ArticleGoogle Scholar
- Brioschi ML, Macedo JF, Macedo RAC: Skin thermometry: new concepts. J Vasc Br 2003, 2(2):151–160.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.