Marcin Sylwestrzak

Efficient reduction of speckle noise in Optical Coherence Tomography

Speckle pattern, which is inherent in coherence imaging, influences significantly axial and transversal resolution of Optical Coherence Tomography (OCT) instruments. The well known speckle removal techniques are either sensitive to sample motion, require sophisticated and expensive sample tracking systems, or involve sophisticated numerical procedures. As a result, their applicability to in vivo real-time imaging is limited. In this work, we propose to average multiple A-scans collected in a fully controlled way to reduce the speckle contrast. This procedure involves non-coherent averaging of OCT A-scans acquired from adjacent locations on the sample. The technique exploits scanning protocol with fast beam deflection in the direction perpendicular to lateral dimension of the cross-sectional image. Such scanning protocol reduces the time interval between A-scans to be averaged to the repetition time of the acquisition system. Consequently, the averaging algorithm is immune to bulk motion of an investigated sample, does not require any sophisticated data processing to align cross-sectional images, and allows for precise control of lateral shift of the scanning beam on the object. The technique is tested with standard Spectral OCT system with an extra resonant scanner used for rapid beam deflection in the lateral direction. Ultrahigh speed CMOS camera serves as a detector and acquires 200,000 spectra per second. A dedicated A-scan generation algorithm allows for real-time display of images with reduced speckle contrast at 6 frames/second. This technique is applied to in vivo imaging of anterior and posterior segments of the human eye and human skin.

Structural Examination of Easel Paintings with Optical Coherence Tomography

Identification of the order, thickness, composition, and possibly the origin of the paint layers forming the structure of a painting, that is, its stratigraphy, is important in confirming its attribution and history as well as planning conservation treatments. The most common method of examination is analysis of a sample collected from the art object, both visually with a microscope and instrumentally through a variety of sophisticated, modern analytical tools. Because of its invasiveness, however, sampling is less than ideally compatible with conservation ethics; it is severely restricted with respect to the amount of material extirpated from the artwork. Sampling is also rather limited in that it provides only very local information. There is, therefore, a great need for a noninvasive method with sufficient in-depth resolution for resolving the stratigraphy of works of art. Optical coherence tomography (OCT) is a noninvasive, noncontact method of optical sectioning of partially transparent objects, with micrometer-level axial resolution. The method utilizes near-infrared light of low intensity (a few milliwatts) to obtain cross-sectional images of various objects; it has been mostly used in medical diagnostics. Through the serial collection of many such images, volume information may be extracted. The application of OCT to the examination of art objects has been in development since 2003. In this Account, we present a short introduction to the technique, briefly discuss the apparatus we use, and provide a paradigm for reading OCT tomograms. Unlike the majority of papers published previously, this Account focuses on one, very specific, use of OCT. We then consider two examples of successful, practical application of the technique. At the request of a conservation studio, the characteristics of inscriptions on two oil paintings, originating from the 18th and 19th centuries, were analyzed. In the first case, it was possible to resolve some questions concerning the history of the work. From an analysis of the positions of the paint layers involved in three inscriptions in relation to other strata of the painting, the order of events in its history was resolved. It was evident that the original text had been overpainted and that the other inscriptions were added later, thus providing convincing evidence as to the paintings true date of creation. In the second example, a painting was analyzed with the aim of confirming the possibility of forgery of the artists signature, and evidence strongly supporting this supposition is presented. These two specific examples of successful use of the technique on paintings further demonstrate how OCT may be readily adaptable to other similar tasks, such as in the fields of forensic or materials science. In a synergistic approach, in which information is obtained with a variety of noninvasive techniques, OCT is demonstrably effective and offers great potential for further development.

Four-dimensional structural and Doppler optical coherence tomography imaging on graphics processing units

Abstract. The authors present the application of graphics processing unit (GPU) programming for real-time three-dimensional (3-D) Fourier domain optical coherence tomography (FdOCT) imaging with implementation of flow visualization algorithms. One of the limitations of FdOCT is data processing time, which is generally longer than data acquisition time. Utilizing additional algorithms, such as Doppler analysis, further increases computation time. The general purpose computing on GPU (GPGPU) has been used successfully for structural OCT imaging, but real-time 3-D imaging of flows has so far not been presented. We have developed software for structural and Doppler OCT processing capable of visualization of two-dimensional (2-D) data (2000 A-scans, 2048 pixels per spectrum) with an image refresh rate higher than 120 Hz. The 3-D imaging of 100×100 A-scans data is performed at a rate of about 9 volumes per second. We describe the software architecture, organization of threads, and optimization. Screen shots recorded during real-time imaging of a flow phantom and the human eye are presented.

Real-time imaging for Spectral Optical Coherence Tomography with massively parallel data processing

The paper presents the way that colour can serve solving the problem of calibration points indexing in a camera geometrical calibration process. We propose a technique in which indexes of calibration points in a black-and-white chessboard are represented as sets of colour regions in the neighbourhood of calibration points. We provide some general rules for designing a colour calibration chessboard and provide a method of calibration image analysis. We show that this approach leads to obtaining better results than in the case of widely used methods employing information about already indexed points to compute indexes. We also report constraints concerning the technique. Nowadays we are witnessing an increasing need for camera geometrical calibration systems. They are vital for such applications as 3D modelling, 3D reconstruction, assembly control systems, etc. Wherever possible, calibration objects placed in the scene are used in a camera geometrical calibration process. This approach significantly increases accuracy of calibration results and makes the calibration data extraction process easier and universal. There are many geometrical camera calibration techniques for a known calibration scene [1]. A great number of them use as an input calibration points which are localised and indexed in the scene. In this paper we propose the technique of calibration points indexing which uses a colour chessboard. The presented technique was developed by solving problems we encountered during experiments with our earlier methods of camera calibration scene analysis [2]-[3]. In particular, the proposed technique increases the number of indexed points points in case of local lack of calibration points detection. At the beginning of the paper we present a way of designing a chessboard pattern. Then we describe a calibration point indexing method, and finally we show experimental results. A black-and-white chessboard is widely used in order to obtain sub-pixel accuracy of calibration points localisation [1]. Calibration points are defined as corners of chessboard squares. Assuming the availability of rough localisation of these points, the points can be indexed. Noting that differences in distances between neighbouring points in calibration scene images differ slightly, one of the local searching methods can be employed (e.g. [2]). Methods of this type search for a calibration point to be indexed, using a window of a certain size. The position of the window is determined by a vector representing the distance between two previously indexed points in the same row or column. However, experiments show that this approach has its disadvantages, as described below. * E-mail: A.Nowakowski@ire.pw.edu.pl Firstly, there is a danger of omitting some points during indexing in case of local lack of calibration points detection in a neighbourhood (e.g. caused by the presence of non-homogeneous light in the calibration scene). A particularly unfavourable situation is when the local lack of detection effects in the appearance of separated regions of detected calibration points. It is worth saying that such situations are likely to happen for calibration points situated near image borders. Such points are very important for the analysis of optical nonlinearities, and a lack of them can significantly influence the accuracy of distortion modelling. Secondly, such methods may give wrong results in the case of optical distortion with strong nonlinearities when getting information about the neighbouring index is not an easy task. Beside this, the methods are very sensitive to a single false localisation of a calibration point. Such a single false localisation can even result in false indexing of a big set of calibration points. To avoid the above-mentioned problems, we propose using a black-and-white chessboard which contains the coded index of a calibration point in the form of colour squares situated in the nearest neighbourhood of each point. The index of a certain calibration point is determined by colours of four nearest neighbouring squares (Fig.1). An order of squares in such foursome is important. Because the size of a colour square is determined only by the possibility of correct colour detection, the size of a colour square can be smaller than the size of a black or white square. The larger size of a black or white square is determined by the requirements of the exact localisation step which follows the indexing of calibration points [3]. In this step, edge information is extracted from a blackand-white chessboard. This edge information needs larger Artur Nowakowski, Wladyslaw Skarbek Institute of Radioelectronics, Warsaw University of Technology, Nowowiejska 15/19, 00-665 Warszawa, A.Nowakowski@ire.pw.edu.pl Received February 10, 2009; accepted March 27, 2009; published March 31, 2009 http://www.photonics.pl/PLP

Depth-Resolved Multilayer Pigment Identification in Paintings: Combined Use of Laser-Induced Breakdown Spectroscopy (LIBS) and Optical Coherence Tomography (OCT)

A detailed feasibility study on the combined use of laser-induced breakdown spectroscopy with optical coherence tomography (LIBS/OCT), aiming at a realistic depth-resolved elemental analysis of multilayer stratigraphies in paintings, is presented. Merging a high spectral resolution LIBS system with a high spatial resolution spectral OCT instrument significantly enhances the quality and accuracy of stratigraphic analysis. First, OCT mapping is employed prior to LIBS analysis in order to assist the selection of specific areas of interest on the painting surface to be examined in detail. Then, intertwined with LIBS, the OCT instrument is used as a precise profilometer for the online determination of the depth of the ablation crater formed by individual laser pulses during LIBS depth-profile analysis. This approach is novel and enables (i) the precise in-depth scaling of elemental concentration profiles, and (ii) the recognition of layer boundaries by estimating the corresponding differences in material ablation rate. Additionally, the latter is supported, within the transparency of the object, by analysis of the OCT cross-sectional views. The potential of this method is illustrated by presenting results on the detailed analysis of the structure of an historic painting on canvas performed to aid planned restoration of the artwork.

Absolute LIBS stratigraphy with optical coherence tomography

In this contribution preliminary studies on the application of Optical Coherence Tomography (OCT) to absolute depth calibration of Laser Induced Breakdown Spectroscopy (LIBS) data in application to revealing stratigraphy of easel paintings are presented. The procedure of in-situ monitoring of LIBS by means of OCT is described. Numerical method developed for precise extraction of the depth of the LIBS ablation crater is explained. Results obtained with model paintings are discussed.

The application of macro-X-ray fluorescence and optical coherence tomography for examination of parchment manuscripts

Abstract Macro-X-ray fluorescence (XRF) is a newly commercially available research tool very useful in the examination of artwork. Its novelty lies in its ability to create maps of the distribution of chemical elements on scales of a few milimetres. In this contribution, its use together with optical coherence tomography (OCT) in the inspection of an illuminated manuscript is reported for the first time. The former technique is used both for mapping the elemental distribution over large parts of the folios – including illuminated initials – and for quantitative analysis of the composition of the smalt pigment, as well as of changes in the composition of iron–gall ink at different pages. The latter, by providing cross-sectional images of painted details, helps in interpreting the XRF results. All of the results shown relate to the examination of a late sixteenth-century-illuminated parchment manuscript (a gradual) originating from the Convent of the Benedictine Sisters in Lviv in the Ukraine.

Picosecond laser ablation system with process control by optical coherence tomography

In this contribution we describe an apparatus for precise laser ablation of delicate layers, like varnish on pictures. This specific case is very demanding. First of all any changes in colour of remaining varnish layer as well as underneath paint layers are unacceptable. This effect may be induced photochemically or thermically. In the first case strong absorption of the radiation used will eliminate its influence on underlying strata. The thermal effect is limited to so called heat affected zone (HAZ). In addition to colour change, a mechanical damage caused by overheating of the structure adjacent to ablated region should be considered also. All kinds of treads must be carefully eliminated in order to make laser ablation of varnish commonly accepted alternative to chemical and/or mechanical treatments [1]. Since the varnish ablation process is obviously irreversible its effective monitoring is very important to make it safe and trusted. As we showed previously [2-6] optical coherence tomography (OCT) originated from medicine diagnostic method for examination and imaging of cross-sections of weakly absorbing objects can be used for this task. OCT utilises infrared light for non-invasive structure examination and has been under consideration for the examining of objects of art since 2004 [7-10]. In this case the in-depth (axial) resolution is obtained by means of interference of light of high spatial (to ensure sensitivity) and very low temporal coherence (to ensure high axial resolution). In practice, IR sources of bandwidths from 25 to 150 nm are utilised. Resolutions obtained range from 15 down to 2 μm in the media of refracting index equal 1.5. In this contribution we expand application of OCT to space resolved determination of ablation rates, separately for every point of examined area. Such data help in better understanding of the ablation process, fine tuning the laser and finally permit increase of the safety of the ablation process.

Real-time massively parallel processing of spectral optical coherence tomography data on graphics processing units

In this contribution we describe a specialised data processing system for Spectral Optical Coherence Tomography (SOCT) biomedical imaging which utilises massively parallel data processing on a low-cost, Graphics Processing Unit (GPU). One of the most significant limitations of SOCT is the data processing time on the main processor of the computer (CPU), which is generally longer than the data acquisition. Therefore, real-time imaging with acceptable quality is limited to a small number of tomogram lines (A-scans). Recent progress in graphics cards technology gives a promising solution of this problem. The newest graphics processing units allow not only for a very high speed three dimensional (3D) rendering, but also for a general purpose parallel numerical calculations with efficiency higher than provided by the CPU. The presented system utilizes CUDATM graphic card and allows for a very effective real time SOCT imaging. The total imaging speed for 2D data consisting of 1200 A-scans is higher than refresh rate of a 120 Hz monitor. 3D rendering of the volume data build of 10 000 A-scans is performed with frame rate of about 9 frames per second. These frame rates include data transfer from a frame grabber to GPU, data processing and 3D rendering to the screen. The software description includes data flow, parallel processing and organization of threads. For illustration we show real time high resolution SOCT imaging of human skin and eye.

OCT structural examination of 'Madonna dei Fusi' by Leonardo da Vinci

Madonna dei Fusi (‘Madonna of the Yarnwider’) is a spectacular example of Italian Renaissance painting, attributed to Leonardo da Vinci. The aim of this study is to give an account of past restoration procedures. The evidence of a former retouching campaign will be presented with cross-sectional images obtained non-invasively with Optical Coherence Tomography (OCT). Specifically, the locations of overpaintings/retouchings with respect to the original paint layer and secondary varnishes will be given. Additionally, the evidence of a former transfer of the pictorial layer to the new canvas support by detecting the presence of its structure incised into paint layer will be shown.