Translator Disclaimer
Open Access
30 August 2022 Registration of polarimetric images for in vivo skin diagnostics
Lennart Jütte, Gaurav Sharma, Harshkumar Patel, Bernhard Roth
Author Affiliations +

Significance: Mueller matrix (MM) polarimetry is a promising tool for the detection of skin cancer. Polarimetric in vivo measurements often suffer from misalignment of the polarimetric images due to motion, which can lead to false results.

Aim: We aim to provide an easy-to-implement polarimetric image data registration method to ensure proper image alignment.

Approach: A feature-based image registration is implemented for an MM polarimeter for phantom and in vivo human skin measurements.

Results: We show that the keypoint-based registration of polarimetric images is necessary for in vivo skin polarimetry to ensure reliable results. Further, we deliver an efficient semiautomated method for the registration of polarimetric images.

Conclusions: Image registration for in vivo polarimetry of human skin is required for improved diagnostics and can be efficiently enhanced with a keypoint-based approach.



Imaging of the human skin with polarimetric techniques has gained importance in recent years, as the incidence of melanoma, which is the deadliest type of skin cancer, is rising throughout all skin types.1 Beyond the classical skin cancer screening done by a dermatologist with a dermoscope,2 new modalities are being implemented for the early detection of melanoma: polarimetric imaging devices such as Stokes3 or Mueller polarimeters4,5 and machine learning-based computer-aided diagnosis systems6 are among the recently reported approaches. Furthermore, the early detection of melanoma is increasingly profiting from automated registration and analysis.7 In addition, adaptive optics such as adjustable liquid lenses allow for the implementation of autofocus in imaging systems without moving parts.8 However, because the up to 36 images necessary for Mueller matrix (MM) determination9 are typically obtained in acquisition times of 20 to 30 s, depending on the type of the setup,10 unintended body movements of the patient can become a significant limitation. Patient motion can lead to misalignment of the polarimetric images.11 Further, it can induce motion blur, in particular because laser power levels have to be reduced due to laser safety and therefore relatively high exposure times are needed to acquire the MMs. In addition, the motion can add up in cases in which a summation process of the polarimetric data is needed to enhance the signal-to-noise ratio (SNR).11 For example, the registration of scanned polarimetric images for ophthalmology was investigated by Nourrit et al.11 The above-mentioned problems in the MM acquisition can be partly overcome by hardware upgrades, e.g., camera and lens upgrades. However, such upgrades typically rely on complex and costly components.12 Therefore, software solutions represent an alternative approach and promise improved results.

In this work, we investigate the need and potential benefits of registering polarimetric images for application in vivo for dermoscopy. Registration of such data is particularly complex because the skin rarely shows sufficient key points with strong contrast (except for mole evaluation13). In addition, the contrast of key points usually changes considerably between the polarization states of the lighting and the analyzer. This issue and the possibility of key points moving out of the field of view make it difficult to detect the same set of key points within the up to 36 images typically acquired in MM polarimetry (MMP). To evaluate the potential benefits of this method, we determine the MM and its polarimetric parameters with and without the polarimetric image registration for the case of a printed melanoma phantom, healthy skin with a drop of honey, and a benign nevus, respectively. The instrumentation and registration schemes are described in Sec. 2. The results of registration and polarimetric analysis are presented in Sec. 3. Section 4 concludes the study.




Instrumentation and Experimental Procedure

With MMP, the information about the polarization changing properties of a sample can be obtained.14 The MM combines the complete polarization properties of a sample within a 4×4 matrix.15 The sample is illuminated at different polarization states. The MM can be calculated from the intensity of the reflected or transmitted light.16 For in vivo skin imaging, only the reflection modality is relevant. Mathematically, the MM is a transformation matrix for Stokes vectors that describe the polarization state of light. The polarization state of the light after interaction with a sample S0 is calculated from the Stokes vector of the incident light Si and the MM entries Mij. of the sample as17

Eq. (1)


These polarization states are described using Stokes vectors18 as follows:

Eq. (2)


The indices of the intensity values I stand for the polarization states as shown in Table 1.

Table 1

Explanation of indices and the correlated polarization states.

PolarizationHorizontalVerticalLinear +45 degLinear −45 degRight circularLeft circular

The Stokes vector changes after the interaction of incident light with the sample, which is recorded to measure the MM. For the used setup, either 16 (needed states: H, V, P, and R) or 36 (needed states: H, V, P, M, R, and L) different images are acquired.8 Although the acquisition of the MM based on 16 images is faster and therefore results in less motion error, the usage of 36 images improves the SNR. More detailed descriptions of the formalism can be found in the literature.17 To reduce calibration measurement errors, in this work, we use six polarization states for the acquisition of the MM from the required 36 intensity measurements to increase the measurement accuracy and the SNR in comparison with the calculation of the MM from 16 images only.19 The calibration of the polarimeter is tested by measuring samples of known MMs such as commercially available polarizers, retarders, and diffusors as well as air. In general, the calibration of our polarimetric device involves two main steps. First, the relative orientation of the optical elements (i.e., linear polarizers and liquid crystal retarders) is set. Second, the physical behavior of the optical elements (i.e., the phase-shift of the liquid crystal retarders) needs to be calibrated by employing quarter- and half-wave-plates. The detailed calibration steps can be found in the literature.20 The results of the calibration procedures performed in this work are very consistent with the expected matrices. An increase in measurement time to improve the SNR needs to be carefully considered, especially in vivo measurements in which movement disorder can change the result.

The experimentally obtained MM entries do not show a direct relation to physical properties of the sample. Therefore, for a more detailed interpretation of the measurements, a polar decomposition is usually carried out. In this work, we use the common polar decomposition of Lu and Chipman21 in the pixel-by-pixel way with the rows of the pixel array indexed as i and the columns indexed as j:

Eq. (3)

Mexp,ij  =  MΔ,ij·MR,ij·MD,ij.

In this decomposition Mexp is the experimentally obtained MM, and MΔ, MR, and MD represent the pure depolarizer, retarder, and diattenuator properties, respectively.21

The key parameters resulting from the polar decomposition are Δ, the depolarization power; R, the total retardance; D, the diattenuation; and P, the total polarizance. These physical properties of the sample are calculated from the following equations:21

Eq. (4)


Eq. (5)


Eq. (6)


Eq. (7)


We use our in-house MM polarimeter22 for the measurements reported in this work. A sketch of the optical elements of the system is shown in Fig. 1.

Fig. 1

Sketch of the MM polarimeter used. The IR distance sensor for the automatic focus system is placed below the path of the beam. A continuous wave (cw) laser is coupled into an optical fiber with an attached vibrational motor for speckle reduction. Together with a fixed linear polarizer orientated at +45  deg with respect to the reference coordinate system, the liquid crystal retarders (LCR) 1 and 2 form the polarization state generator. The polarization state analyzer (PSA) contains two LCRs and a fixed linear polarizer orientated at 45  deg with respect to the reference coordinate system.


The distance sensor (DT35-B15851, Sick AG, Waldkirch, Germany) measures the distance from a point close to (but outside) the polarimeter’s field of view to avoid interference with its illumination, as visualized in Fig. 2.

Fig. 2

CAD design of autofocus system for a MM polarimeter.


With the liquid lens, it is possible to adjust the focus in real-time based on the distance information provided by the distance sensor and on a function derived from the calibration procedure, which relies on the collection of value pairs of working distance and the liquid lens current. In detail, the liquid lens current is adjusted for different working distances until the imaged plane object is in maximum focus.23 The automatic focus works well for most skin parts within a theoretical working distance range of 30 cm to 8 m.8 However, the working distance is limited by the illumination intensity (and therefore noise) and the desired image resolution. With the automatic focus, it is possible to fully open the imaging aperture and accept the reduced depth of field, as larger openings reduce the acquisition time due to shorter possible exposure times.

To minimize patient motion, we designed a mobile MM acquisition system based on an ergonomic arm (Ergotron LX ARM 45-241-026, Ergotron, St. Paul, Minnesota) that allows patients to lie on a patient bed. Figure 3 shows a CAD design of the complete measurement system.

Fig. 3

CAD design of the in vivo skin MMP imaging setup. The patient movement is minimized through the possibility of the patient lying down.


The mobile MM acquisition system allows the operator to image essentially all parts of the human skin. For better visualization, the appropriate laser safety measures are not shown in this representation.


Image Processing

Polarimetric datasets are obtained with the MMP only under unintentional movement (e.g., breathing and shaking). In the following, we describe the different image registration techniques used in this work to correct for the misalignment caused by that motion. First, registration assessment is conducted on the basis of subjective criteria. In addition, an objective criterion is applied with the comparison of the numerical values of each polarimetric parameter, respectively.


Image Registration of Polarimetric Images

Image registration is the process of aligning two or more images of the same scene.24 It is often used as a precondition for other image processing applications.25 Determining an effective approach to image registration depends on the application case. Careful selection of a point transformation model is required to provide reference points between the images. In addition, a method for comparing information to identify the parameters necessary to correctly align images is needed. In general, there are two well-known methods for automatic image registration: feature-based and intensity-based registration algorithms.26 By contrast, manual image registration relies on control point mapping registration algorithms.27

Local features and their descriptors constitute the basis for many computer vision algorithms.28 Their applications include image registration,29 object detection30 and classification,31 tracking,32 motion estimation,33 and content-based image retrieval.34 Local features refer to a pattern or distinct structure that exists in an image, such as a point, edge, or small patch of the image. They are usually associated with an image patch that differs from its immediate environment in texture, color or intensity. Although the distinction from its environment is important, the representation of the content of the characteristic is usually not relevant. Examples of local features are blobs, corners, and edge pixels.24 Due to the lack of matching points in skin images, feature matching is not a reliable method for batch image processing alone.

The Image Processing Toolbox™ by MATLAB (MATLAB, 2021. version 9.11.0 (R2021b), Natick, Massachusetts: The MathWorks Inc.) provides tools for point mapping to determine the required transformation parameters to align images. In point mapping, the user selects points in a pair of images that identify the same reference. Then a geometric mapping is deduced from the positions of these control points.24 In selecting control points, a high level of accuracy is required. To align the target image, we need to select at least two pairs of matching points between the target image and the reference image. More pairs of matching points improve the registration result. This increases the processing time for batch image processing.

We propose a semiautomated key point-based registration method for polarimetric images to take into account the inevitable movement during acquisition. We combine the registration of both feature-based (automatic) and control point mapping (manual). Figure 4 shows a visual representation of this process.

Fig. 4

Flow chart of semiautomated keypoint-based registration with an example of a heart symbol printed on paper.


The process starts with an algorithm based on KAZE35 features. In the case of insufficient pairs of matching points, it will switch to the manual mapping algorithm of the control points for this particular image. The Computer Vision Toolbox™ provides several methods for detecting corner and blob features and includes various descriptors. In this method, feature detection selects regions of an image with unique content, such as blobs. The feature detection finds possible points for subsequent processing that do not necessarily refer to physical structures. Finding features that remain locally invariant so that they are identifiable even in the presence of rotation or scaling is essential to feature detection.

Feature extraction involves the calculation of a descriptor, which is typically done on regions centered around detected features. By the means of image processing, descriptors transform a local pixel neighborhood into a compact vector representation. This new representation allows for comparison between neighborhoods, regardless of changes in scale or direction. We chose the best feature detector and descriptor from the perspective of our application criteria and the nature of our data by testing the available feature detector and descriptor combinations. Table 2 shows a comparison of the number of matching pairs between two images for different blob detectors and descriptors.

Table 2

Comparison of the number of matching pairs between two images for different blob detector and descriptor combinations. We exemplarily show the values for two images of the polarimetric raw data of a nevus (compare with Fig. 15) and a honey drop on healthy skin (compare with Fig. 13).

Blob detectorsDescriptorNevusHoney drop on healthy skin
Features detectedMatching pairs between reference and target imagesFeatures detectedMatching pairs between reference and target images
Reference imageTarget imageReference imageTarget image
Note: The detector and descriptor combinations employed in this article are highlighted in bold.

If there are sufficient valid matches, the false matches are removed by selecting only the strongest among them. This work is based on the M-estimator sample consensus (MSAC),36 a variant of the random sample consensus (RANSAC)37 algorithm. This method finds a geometric transform, separates the correct matches from the spurious matches, and then applies the computed geometric transform to the image being compared.36 For this method, a minimum of two matching pairs is required. If the requirement is not met, then it allows the operator to manually select the two matching pairs between the references and the particular target image. The manual selection of control points is shown in Fig. 5.

Fig. 5

Screenshot of the manual control point selection process. The user selects the control points by clicking on a feature in the top-left image and then the corresponding feature in the top-right image. The lower row images are used for overview. The blue box region is shown enlarged in the upper row of images.


After alignment of all images, the MM is calculated. However, the edges of the spatially resolved MM entries are usually not valid and show distracting extreme values. To achieve results without distracting edges, we crop all images to the reference overlap zone of all aligned images or the area of interest.



Figure 6 shows the samples studied in a first step that were captured by a smartphone (upper row) and by a noncontact dermoscope23 (lower row).

Fig. 6

Clinical images of (a) healthy skin with hair and (b) a benign nevus. Dermoscopy images of (c) healthy skin with hair and (d) a benign nevus as a standard skin imaging reference to MM imaging.


We implement our approach on healthy skin with hair and a mole with surrounding skin, respectively. In addition, we apply a tiny drop of honey to the healthy skin as a phantom for small skin patches with varying optical activity.38 Aligning these samples is particularly important due to their strong contrast resulting from the hair and the nevus. The clinical and dermoscopic images serve as reference for the polarimetric skin imaging.


Registration of Polarimetric Data for In Vivo Skin Imaging

First, we show that in polarimetric data the visibility of features in the skin is usually polarization sensitive, as shown in Fig. 7.

Fig. 7

The visibility of features in human skin can depend on the polarization states. Two features with strong polarization dependency are marked exemplarily.


From the comparison of the two images from the polarimetric raw data in Fig. 7, it is obvious that the visibility of features in human skin can depend on the combination of the state of the polarization state generator and the PSA.

In the following, we demonstrate that movement during data acquisition can lead to false MMs and derived polarimetric parameters. To do this, we use an icon as shown in Fig. 8 printed on paper as a lesion phantom.

Fig. 8

Phantom designed based on ABCD criteria.


The phantom design is based on the ABCD criteria for melanoma diagnosis. It is asymmetric and has irregular borders and a diameter >5  mm. The phantom is moved throughout the acquisition. The translational and rotational movement within the focal plane, realized by manually moving the target on a rail with one additional rotational degree of freedom throughout the acquisition, leads to an error of alignment as shown in Fig. 9.

Fig. 9

The detected matching features are at different positions in two exemplary overlayed raw images, as indicated by the yellow lines connecting corresponding matching points. This error of orientation is due to motion during the acquisition.


In the following the spatially resolved polarimetric parameters retardance, polarizance, depolarization power, and diattenuation are calculated. In the remainder of this work, we refer to our feature-based registration method as alignment (a) and the standard intensity-based image registration method as alignment (b). As displayed in Fig. 10, the results differ largely depending on whether the polarimetric data were registered or not and which registration method was applied.

Fig. 10

Polarimetric parameters of the phantom. (a) Not moved, (b) moved and unaligned, (c) feature-based alignment, and (d) intensity-based alignment.


The previous comparison shows that registration has two immediate positive effects. First, image alignment leads to an improved image quality, as the features of the phantom are more clearly visible. Second, alignment corrects spatially resolved polarimetric parameters. The trend of retardance, polarization, and diattenuation of the phantom is reversed if the images are not aligned. The aligned and cropped data parameters indicate that the alignment has been achieved, but that the values on the phantom edges are still different compared with the static reference sample, as seen in the polarizance panels in Fig. 10, third and fourth rows. As shown in Fig. 10, it is also obvious that the polarimetric parameters from the raw data that have been registered with the feature-based method shows a greater agreement with the polarimetric data from the data at rest (ground truth) compared with the application of the standard intensity-based method. Further, the feature-based method is less prone to artifacts on the edges of the phantom.

As shown in Table 3, it is obvious that the feature-based registration leads to results that are closest to the ground truth for R, P, and Δ. For D, both alignment methods give a similar numerical value that does not match with the ground truth.

Table 3

Comparison of the numerical values of the determined polarimetric parameters from the center of the phantom.

Alignment (a)0.90770.09880.18200.05990.43140.07630.31030.0726
Alignment (b)2.18260.40730.25670.19560.55740.61530.30710.2743
Note: Significant values are highlighted in bold.

Next, a polarimetric in vivo measurement of a honey drop on healthy skin of a human arm with hair, as shown in Fig. 11, is carried out.

Fig. 11

Exemplary raw data of the honey drop on healthy skin as a phantom for spatially defined small changes in polarization on human skin. The contours of the honey drop are marked.


The diameter of the honey drop is 2  mm. Figure 12 shows the degree of unintentional motion during data collection by comparing the position of detected features in two images of the polarimetric data.

Fig. 12

Exemplary resulting misalignment of the honey drop on healthy skin due to motion during the acquisition. Two exemplary images of the raw data are overlayed. The yellow arrows connect corresponding matching points. The region of the honey drop is marked in both images.


Figure 12 shows the degree of unintentional motion and the need for registration of the polarimetric data. The spatially resolved polarimetric parameters for a drop of honey on healthy skin with hair are displayed in Fig. 13 for the raw data as well as after feature-based and intensity-based alignment.

Fig. 13

Polarimetric parameters for a drop of honey on healthy skin with hair. (a) unaligned, (b) feature-based alignment, and (c) intensity-based alignment. The region of the honey drop is marked. The contours are best identified in the case of the feature-based alignment due to the greatest contrast.


Feature-based registration of polarimetric images has an immediate effect on spatially resolved polarimetric parameters. After registration, the hair can be distinguished from the skin, and the range of values has changed. In addition, the tiny drop of honey is most distinguishable from the surrounding skin in the case of feature-based registration, especially in the Δ-patch, showing that such polarimetric changes in small skin patches can only be restored with a suitable image registration method.

Subsequently, the influence of in-focal plane motion on the polarimetric parameters for in vivo nevus assessment was investigated. An in vivo measurement with our MMP is performed on a benign nevus of a volunteer. The degree of misalignment is obvious from Fig. 14, as the matching feature points are detected at slightly different positions.

Fig. 14

Exemplary resulting misalignment of matching feature points of the nevus due to motion during the acquisition. Two exemplary images of the raw data are overlayed. The yellow arrows connect corresponding matching points.


It is clear from Fig. 14 that the features of the nevus do not overlap in the polarimetric images due to motion and that the polarimetric data must therefore be registered. In Fig. 15, the polarimetric parameters of a nevus of unaligned and aligned data are shown.

Fig. 15

Polarimetric parameters for a mole. (a) unaligned, (b) feature-based alignment, and (c) intensity-based alignment.


The feature-based registration of polarimetric data allows for the in vivo assessment of nevi. By contrast, unaligned or intensity-based aligned polarimetric data lead to artifacts in the polarimetric parameters that interfere with dermatological interpretation. The inconsistencies of polarimetric parameters at the edge of the nevus, as obvious from the retardance patch in the second row in Fig. 15, suggest that the image registration can be further improved. It is necessary to examine, in future work, whether the inconsistencies at the edges of the nevus are caused by the occurrence of biological and structural changes or by a slight residual misalignment.


Autofocus for Polarimetric In Vivo Skin Imaging

Finally, we investigated the effect of movement in the direction of the optical axis of the polarimetric system and the associated changes in the polarimetric parameters. For this, the focal plane was changed during the acquisition of the polarimetric parameters displayed in Fig. 16.

Fig. 16

Polarimetric parameters from data of the phantom with motion in the direction of the optical axis without autofocus. The induced artifacts are obvious in comparison with the phantom at rest in Fig. 9, first row.


In comparison with Fig. 10, first row, it is clear that the change of the focus plane during acquisition leads to a bad focus on some of the images of the polarimetric data, which results in artifacts on the edges of the phantom. Consequently, an automatic focus system is required to take into account the skin movement in the direction of the optical axis during acquisition.



We implemented a semiautomatic feature-based registration method and studied the impact of different registration methods on polarimetric analysis. In addition, we found mapping functions that are suitable to reliably detecting the corresponding skin features as key points in polarimetric data. The body movement can induce motion blur because of relatively high exposure times. Furthermore, the movement of the skin in the 36 images acquired leads to false polarimetric parameters, preventing reliable diagnosis. In rare cases, some key points leave the field of view due to movement. To avoid false MM results, the characteristics of the 36 images are aligned using image processing techniques. In addition, the MMP is mounted on a movable arm that enables the system to be placed easily in the skin area studied while the patient movement is minimized with the patient bed setting. The results indicate that our methodology is a good solution to the problem. Further research into the feature detection technique is needed to enable a fully automated and reliable registration of in vivo skin polarimetry data. The results emphasize the importance of adequate image registration techniques for in vivo skin polarimetry. Thus, equipped with image registration, MM polarimetry could be a valuable asset for dermatology. The performance of image registration is most important when the skin area of interest is particularly small, as characteristics of polarimetric activity can be lost without data registration. Our proposed method is suitable for the alignment of polarimetric images overcoming the problems of intensity and threshold matching. In the future, we will use the approach for in vivo measurements of inflammatory skin diseases and melanoma skin cancer.


The authors declare no conflicts of interest.


This research was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy within the Cluster of Excellence PhoenixD (EXC 2122, Project ID 390833453). This work has been supported by iToBoS (Intelligent Total Body Scanner for Early Detection of Melanoma), project funded by the European Union’s Horizon 2020 research and innovation program, under grant agreement no. 965221

Code, Data, and Materials Availability

Data are available upon reasonable request.



F. Bray et al., “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: A Cancer J. Clin., 68 (6), 394 –424 (2018). Google Scholar


N. G. Marghoob, K. Liopyris and N. Jaimes, “Dermoscopy: a review of the structures that facilitate melanoma detection,” J. Am. Osteopathic Assoc., 119 (6), 380 –390 (2019). Google Scholar


D. C. Louie et al., “Degree of optical polarization as a tool for detecting melanoma: proof of principle,” J. Biomed. Opt., 23 (12), 125004 (2018). JBOPFO 1083-3668 Google Scholar


D. Fricke et al., “Non-contact fast Mueller matrix measurement system for investigation of inflammatory skin diseases,” Progr. Biomed. Opt. Imaging, 20 (13), 6 (2019). Google Scholar


N. T. Luu et al., “Characterization of Mueller matrix elements for classifying human skin cancer utilizing random forest algorithm,” J. Biomed. Opt., 26 (7), 075001 (2021). JBOPFO 1083-3668 Google Scholar


L. R. Soenksen et al., “Using deep learning for dermatologist-level detection of suspicious pigmented skin lesions from wide-field images,” Sci. Transl. Med., 13 (581), eabb3652 (2021). STMCBQ 1946-6234 Google Scholar


K. Korotkov and R. Garcia, “Computerized analysis of pigmented skin lesions: a review,” Artif. Intell. Med., 56 (2), 69 –90 (2012). AIMEEW 0933-3657 Google Scholar


D. Fricke, Entwicklung einer dermatoskopischen Kamera mit Müller-Matrix-Bildgebung, Institutionelles Repositorium der Leibniz Universität Hannover, Hannover (2020). Google Scholar


J. S. Baba et al., “Development and calibration of an automated Mueller matrix polarization imaging system,” J. Biomed. Opt., 7 (3), 341 –349 (2002). JBOPFO 1083-3668 Google Scholar


T. Huang et al., “Fast Mueller matrix microscope based on dual DoFP polarimeters,” Opt. Lett., 46 (7), 1676 –1679 (2021). OPLEDP 0146-9592 Google Scholar


V. Nourrit et al., “Nonlinear registration for scanned retinal images: application to ocular polarimetry,” Appl. Opt., 47 (29), 5341 –5347 (2008). APOPAI 0003-6935 Google Scholar


J. C. Ramella-Roman et al., “Development and characterization of a snapshot Mueller matrix polarimeter for the determination of cervical cancer risk in the low resource setting,” Progr. Biomed. Opt. Imaging, 19 (19), 33 (2018). Google Scholar


C. N. E. Anagnostopoulos, D. D. Vergados and P. Mintzias, “Image registration of follow-up examinations in digital dermoscopy,” in 13th IEEE Int. Conf. BioInf. and BioEng., 1 –4 (2013). Google Scholar


M. Peyvasteh et al., “Evolution of raw meat polarization-based properties by means of Mueller matrix imaging,” J. Biophotonics, 14 (5), e202000376 (2021). Google Scholar


L. Jütte et al., “Mueller matrix-based approach for the ex vivo detection of riboflavin-treated transparent biotissue,” Appl. Sci., 11 (23), 11515 (2021). Google Scholar


D. Fricke et al., “Mueller matrix measurement of electrospun fiber scaffolds for tissue engineering,” Polymers, 11 (12), 2062 (2019). Google Scholar


W. S. Bickel and W. M. Bailey, “Stokes vectors, Mueller matrices, and polarized scattered light,” Am. J. Phys., 53 (5), 468 –478 (1985). AJPIAS 0002-9505 Google Scholar


D. Ivanov et al., “Colon cancer detection by using Poincaré sphere and 2D polarimetric mapping of ex vivo colon samples,” J. Biophotonics, 13 (8), e202000082 (2020). Google Scholar


B. D. Cameron, Y. Li and A. Nezhuvingal, “Determination of optical scattering properties in turbid media using Mueller matrix imaging,” J. Biomed. Opt., 11 (5), 054031 (2006). JBOPFO 1083-3668 Google Scholar


E. Compain, S. Poirier and B. Drevillon, “General and self-consistent method for the calibration of polarization modulators, polarimeters, and Mueller-matrix ellipsometers,” Appl. Opt., 38 (16), 3490 –3502 (1999). APOPAI 0003-6935 Google Scholar


S.-Y. Lu and R. A. Chipman, “Interpretation of Mueller matrices based on polar decomposition,” J. Opt. Soc. Am. A, 13 (5), 1106 (1996). JOAOD6 0740-3232 Google Scholar


D. Fricke et al., “Mueller matrix analysis of collagen and gelatin containing samples towards more objective skin tissue diagnostics,” Polymers, 12 (6), 1400 (2020). Google Scholar


D. Fricke et al., “Non-contact dermatoscope with ultra-bright light source and liquid lens-based autofocus function,” Appl. Sci., 9 (11), 2177 (2019). Google Scholar


B. Zitová and J. Flusser, “Image registration methods: a survey,” Image Vis. Comput., 21 (11), 977 –1000 (2003). Google Scholar


M. Irani and S. Peleg, “Improving resolution by image registration,” CVGIP: Graph. Models Image Process., 53 (3), 231 –239 (1991). Google Scholar


F. P. M. Oliveira and J. M. R. S. Tavares, “Medical image registration: a review,” Comput. Methods Biomech. Biomed. Eng., 17 (2), 73 –93 (2014). Google Scholar


S. Kumar et al., “Robust image registration technique for SAR images,” in First Int. Conf. Ind. and Inf. Syst., 519 –524 (2006). Google Scholar


J. Li and N. M. Allinson, “A comprehensive review of current local features for computer vision,” Neurocomputing, 71 (10–12), 1771 –1787 (2008). NRCGEO 0925-2312 Google Scholar


Y. Lu et al., “A novel image registration approach via combining local features and geometric invariants,” PloS One, 13 (1), e0190383 (2018). POLNCL 1932-6203 Google Scholar


S. Nigam et al., “An effective local feature descriptor for object detection in real scenes,” in IEEE Conf. Inf. & Commun. Technol., 244 –248 (2013). Google Scholar


X. Zhou et al., “Image Classification Using Super-Vector Coding of Local Image Descriptors,” 141 –154 Springer Berlin Heidelberg, Berlin, Heidelberg (2010). Google Scholar


F. Pernici and A. Del Bimbo, “Object tracking by oversampling local features,” IEEE Trans. Pattern Anal. Mach. Intell., 36 (12), 2538 –2551 (2014). ITPIDJ 0162-8828 Google Scholar


H. Huang et al., “Correlation and local feature based cloud motion estimation,” in ACM Conf., 1 –9 (2012). Google Scholar


S. Murala, R. P. Maheshwari and R. Balasubramanian, “Local tetra patterns: a new feature descriptor for content-based image retrieval,” IEEE Trans. Image Process.: A Publ. IEEE Signal Process. Soc., 21 (5), 2874 –2886 (2012). Google Scholar


P. F. Alcantarilla, A. Bartoli and A. J. Davison, KAZE Features, 214 –227 Heidelberg (2012). Google Scholar


S. Wu, W. Zeng and H. Chen, “A sub-pixel image registration algorithm based on SURF and M-estimator sample consensus,” Pattern Recognit. Lett., 140 261 –266 (2020). PRLEDG 0167-8655 Google Scholar


M. A. Fischler and R. C. Bolles, “Random sample consensus,” Commun. ACM, 24 (6), 381 –395 (1981). CACMA2 0001-0782 Google Scholar


R. Espinosa-Luna et al., “Polarimetric applications to identify bee honey,” 82870F (2012). Google Scholar


Lennart Jütte received his BS and MS degrees in nanotechnology from Leibniz University Hannover in 2019 and 2021, respectively, where he is currently pursuing his PhD in applied optics at the Hannover Centre for Optical Technologies. His current research interests include dermoscopy, biophotonics, and polarimetry.

Bernhard Roth received his PhD in atomic and particle physics from the University of Bielefeld, Bielefeld, Germany, in 2001 and his State Doctorate (Habilitation) degree in experimental quantum optics from the University of Duesseldorf, Duesseldorf, Germany, in 2007. From 2002 to 2007, he was a Research Group Leader at the University of Duesseldorf. Since 2012, he has been the scientific and managing director at the Hanover Centre for Optical Technologies, and since 2014, he has been a professor of physics at Leibniz University Hanover, Hanover, Germany,. His scientific activities include applied and fundamental research in laser development and spectroscopy, polymer optical sensing, micro- and nanooptics fabrication, and optical technology for illumination, information technology, and the life sciences.

Biographies of the other authors are not available.

© The Authors. Published by SPIE under a Creative Commons Attribution 4.0 International License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Lennart Jütte, Gaurav Sharma, Harshkumar Patel, and Bernhard Roth "Registration of polarimetric images for in vivo skin diagnostics," Journal of Biomedical Optics 27(9), 096001 (30 August 2022).
Received: 2 May 2022; Accepted: 12 August 2022; Published: 30 August 2022


Image registration


In vivo imaging



Image processing

Back to Top