Optical transfer diagnosis differentiating benign and malignant pigmented lesions in a simulated primary care practice.

BACKGROUND: The detection of melanoma poses a substantial challenge, particularly for primary care providers (PCPs) who may have limited training in discriminating between suspicious and benign melanocytic lesions. The noninvasive optical transfer diagnosis (OTD) method was designed to be used by PCPs in their decision-making process. OBJECTIVES: To assess the potential of the OTD method by developing, training and validating an OTD indication algorithm for automated discrimination between benign melanocytic lesions and malignant lesions, based on a set of 712 lesions. METHODS: The authors performed in vivoOTD capture and subsequent analysis of 712 pigmented lesions. Of the lesions, 415 were clinically and dermoscopically benign and 297 were dermoscopically suspicious or equivocal. After image capture, all suspicious or equivocal lesions were biopsied and examined histopathologically. RESULTS: Of the 297 suspicious or equivocal lesions, histopathological findings revealed 80 to be malignant (64 melanomas, 13 basal cell carcinomas and 3 squamous cell carcinomas). OTD misdiagnosed one of the 80 malignant lesions as benign (sensitivity, 99%). OTD specificity was 93% for the dermoscopically benign lesions, 73% for all lesions included in the study and 36% for the clinically suspicious but histopathologically benign lesions. CONCLUSIONS: High sensitivity and specificity, as provided by OTD in this preliminary study, would help PCPs reduce the number of referrals for dermatology consultation, excision or biopsy. Further studies are planned for screening patients in a primary care setting, with comparisons of OTD results with biopsy or dermoscopy results.

Simultaneous polarimeter retrievals of microphysical aerosol and ocean color parameters from the "MAPP" algorithm with comparison to high-spectral-resolution lidar aerosol and ocean products.

We present an optimal-estimation-based retrieval framework, the microphysical aerosol properties from polarimetry (MAPP) algorithm, designed for simultaneous retrieval of aerosol microphysical properties and ocean color bio-optical parameters using multi-angular total and polarized radiances. Polarimetric measurements from the airborne NASA Research Scanning Polarimeter (RSP) were inverted by MAPP to produce atmosphere and ocean products. The RSP MAPP results are compared with co-incident lidar measurements made by the NASA High-Spectral-Resolution Lidar HSRL-1 and HSRL-2 instruments. Comparisons are made of the aerosol optical depth (AOD) at 355 and 532 nm, lidar column-averaged measurements of the aerosol lidar ratio and Ångstrøm exponent, and lidar ocean measurements of the particulate hemispherical backscatter coefficient and the diffuse attenuation coefficient. The measurements were collected during the 2012 Two-Column Aerosol Project (TCAP) campaign and the 2014 Ship-Aircraft Bio-Optical Research (SABOR) campaign. For the SABOR campaign, 73% RSP MAPP retrievals fall within ±0.04 AOD at 532 nm as measured by HSRL-1, with an R value of 0.933 and root-mean-square deviation of 0.0372. For the TCAP campaign, 53% of RSP MAPP retrievals are within 0.04 AOD as measured by HSRL-2, with an R value of 0.927 and root-mean-square deviation of 0.0673. Comparisons with HSRL-2 AOD at 355 nm during TCAP result in an R value of 0.959 and a root-mean-square deviation of 0.0694. The RSP retrievals using the MAPP optimal estimation framework represent a key milestone on the path to a combined lidar+polarimeter retrieval using both HSRL and RSP measurements.

Long-term comparisons of UV index values derived from a NILU-UV instrument, NWS, and OMI in the New York area.

A comparison is presented of UV index (UVI) values obtained under different cloud conditions from a Norwegian Institute for Air Research UV (NILU-UV) instrument, the ozone monitoring instrument (OMI) onboard the Aura satellite, and the National Weather Service (NWS) model for the time period of 2010-2014. The NILU-UV irradiance meter is a ground-based, multi-channel, moderate bandwidth filter instrument. UVI values derived from measurements by a NILU-UV instrument deployed in the New York area (40.74°N, -74.03°E) to monitor the erythemal UV radiation from 2010 to present is compared to UVI values derived from OMI measurements and predicted by the NWS model. OMI overestimated the UVI values by 13.06% for all cloud conditions compared with the UVI values derived from measurements by the NILU-UV instrument. The heavier the cloud cover, the higher the overestimation. The mean relative difference between the UVI derived from the NWS model and from NILU-UV measurements was 11.54%. The UVI prediction by NWS was also overestimated under cloudy conditions. Under overcast conditions the NWS predictions of UVI values differ significantly from those derived from NILU-UV measurements, yielding a correlation of only 0.8025 and a mean relative difference of 28.25%.

Comparisons of three NILU-UV instruments deployed at the same site in the New York area.

The Norwegian Institute for Air Research ultraviolet (NILU-UV) irradiance meter is a ground-based, multichannel, moderate bandwidth filter instrument that measures irradiances at ultraviolet (UV) and visible wavelengths with five channels in the UV (302, 312, 320, 340, and 380 nm) and one channel in the visible (400-700 nm) part of the spectrum. Minute-by-minute irradiances recorded in these channels are used to infer the total ozone column (TOC) amount, and a radiation modification factor (RMF) designed to have a value close to 100 under cloud-free conditions. The performance of three NILU-UV instruments deployed side-by-side in the New York area (40.74°N, -74.03°E) is assessed, and derived TOC values are compared with those derived from the ozone monitoring instrument (OMI) deployed on NASA's AURA satellite. Based on about three years of data, it was found that the three instruments yielded similar TOC values that were in close agreement with those derived from the OMI. The relative difference in TOC values derived from the three NILU-UV instruments was generally less than 2.5%. Cloud cover affects the accuracy of the inferred TOC, but reliable values can be obtained in the presence of clouds, although the accuracy deteriorates under heavy overcast conditions with RMF values smaller than 65 (low cloud transmittance).

Comparison of discrete ordinate and Monte Carlo simulations of polarized radiative transfer in two coupled slabs with different refractive indices.

A comparison is presented of two different methods for polarized radiative transfer in coupled media consisting of two adjacent slabs with different refractive indices, each slab being a stratified medium with no change in optical properties except in the direction of stratification. One of the methods is based on solving the integro-differential radiative transfer equation for the two coupled slabs using the discrete ordinate approximation. The other method is based on probabilistic and statistical concepts and simulates the propagation of polarized light using the Monte Carlo approach. The emphasis is on non-Rayleigh scattering for particles in the Mie regime. Comparisons with benchmark results available for a slab with constant refractive index show that both methods reproduce these benchmark results when the refractive index is set to be the same in the two slabs. Computed results for test cases with coupling (different refractive indices in the two slabs) show that the two methods produce essentially identical results for identical input in terms of absorption and scattering coefficients and scattering phase matrices.

Modeling optical properties of human skin using Mie theory for particles with different size distributions and refractive indices.

We used size distributions of volume equivalent spherical particles with complex refractive index to model the inherent optical properties (IOPs) in four different layers of human skin at ten different wavelengths in the visible and near-infrared spectral bands. For each layer, we first computed the size-averaged absorption coefficient, scattering coefficient, and asymmetry factor for the collection of particles in a host medium using Mie theory and compared these IOPs in each layer with those obtained from a bio-optical model (BOM). This procedure was repeated, using an optimization scheme, until satisfactory agreement was obtained between the IOPs obtained from the particle size distribution and those given by the BOM. The size distribution as well as the complex refractive index of the particles, obtained from this modeling exercise, can be used to compute the phase matrix, which is an essential input to model polarized light transport in human skin tissue.

Retrieval of the physiological state of human skin from UV-Vis reflectance spectra - a feasibility study.

We test the feasibility of using an accurate radiative transfer model for the coupled air-tissue system in conjunction with a classic inversion scheme based on Bayesian optimal estimation theory for retrieval of parameters describing the physiological state of human skin. To that end, we analyse ultraviolet and visible reflectance spectra from human skin measured before, immediately after, and on each day for two weeks after photodynamic treatment with the hexyl ester of ALA and exposure to red light (632 nm). For the first time, we show that it is possible to perform a simultaneous retrieval of the melanosome concentration in both the basal and the upper layers of the epidermis.

Optical detection and monitoring of pigmented skin lesions.

A method is presented for discriminating between malignant and benign pigmented skin lesions based on multispectral and multi-angle images. It is discussed how to retrieve maps of physiology properties and morphometric parameters from recorded images using a bio-optical model, radiative transfer calculations, and nonlinear inversion, and how to employ automated zooming to extract lesion and surrounding masks. Training and validation of a classification scheme for separation between benign and malignant tissue yielded sensitivity/specificity ranging from 97%/97% for application to a small dataset comprised of lesions not used for training and validation to 99%/93% for application to a larger dataset.

Ozone depletion and its consequences for the fluence of carcinogenic sunlight.

A slight reduction of the ozone level over the northern hemisphere in the period 1969-1986 has been reported [D. Lindley, Nature (Lond.), 323: 293, 1988]. Ozone measurements performed in Oslo are in agreement with this. However, the ozone level for 1987 and 1988 was above normal, and no negative or positive trend is apparent for the last 10 years. The consequences of an ozone reduction for the fluence rate of carcinogenically effective sunlight was evaluated on the basis of recent action spectra for mutagenesis in cells, carcinogenesis in mice, and erythema induction in humans. Depending on the choice of action spectrum we find amplification factors (defined as percentage increase in yearly fluence of carcinogenically efficient sunlight per percentage reduction of the ozone level) between 1.1 and 1.3 at latitudes between 0 and 20 degrees and between 0.9 and 1.1 for Northern Europe. These estimates are significantly lower than 2.0, which is the value found when the calculations are based on the DNA absorption spectrum (R. B. Setlow, Proc. Natl. Acad. Sci. USA, 71:3363-3366, 1974).

Biological UV-doses and the effect of an ozone layer depletion.

Effective UV-doses were calculated based on the integrated product of the biological action spectrum (the one proposed by IEC, which extends to 400 nm, was adopted) and the spectral irradiance. The calculations include absorption and scattering of UV-radiation in the atmosphere, both for normal ozone conditions as well as for a depleted ozone layer. For Scandinavian latitudes the effective annual UV-dose increases by approximately 4% per degrees of latitude towards the Equator. An ozone depletion of one percent increases the annual UV-dose by approximately 1% at 60 degrees N (increases slightly at lower latitudes). A large depletion of 50% over Scandinavia (60 degrees N) would give these countries an effective UV-dose similar to that obtained, with normal ozone conditions, at a latitude of 40 degrees N (California or the Mediterranean countries). The Antarctic ozone hole increases the annual UV-dose by 20 to 25% which is a similar increase as that attained by moving 5 to 6 degrees of latitude nearer the Equator. The annual UV-dose at higher latitudes is mainly determined by the summer values of ozone. Both the ozone values and the effective UV-doses vary from one year to another (within +/- 4%). No positive or negative trend is observed for Scandinavia from 1978 to 1988.

Time-resolved polarimetry over water waves: relating glints and surface statistics.

The phenomenon of sunglint, well known in satellite remote sensing, lacks a fundamental characterization under controlled laboratory conditions. Exploiting an apparatus specifically assembled for the purpose, we examine the signal collected by a photopolarimeter, pointed at a wavy water surface with measurable statistics and illuminated by a laser source. We also analyze the wave slope distributions, retrieved with an imaging system, and correlate them with the time series of glints. More particularly, we investigate the link between the occurrence of glints and that of the slopes from which they originate. In this context, the results obtained by applying the Hilbert-Huang transform technique to the slope time series are compared with those obtained through a traditional Fourier transform. This novel study first identifies the individual atomic glints as Fresnel reflection originating from a single wave facet. It then discusses the periodic character of a sequence of glints generated by a gravity wave state, as opposed to the erratic behavior of glints typical of capillary wave states. In mixed gravity-capillary conditions, it is shown that the glint properties are governed mainly by the capillary regime.

Use of shadows to retrieve water vapor in hazy atmospheres.

Techniques aimed at retrieving water vapor from satellite data of reflected near-infrared solar radiation have progressed significantly in recent years. These techniques rely on observation of water vapor attenuation of near-infrared solar radiation reflected by the Earth's surface. Ratios of measured radiances at wavelengths inside and outside water vapor absorbing channels are used for retrieval purposes. These ratios partially remove the dependence of surface reflectance on wavelength and are used to retrieve the total column water vapor amount. Hazy atmospheric conditions, however, introduce errors into this widely used technique. A new method based on radiance differences between clear and nearby shadowed surfaces, combined with ratios between water vapor absorbing and window regions, is presented that improves water vapor retrievals under hazy atmospheric conditions. Radiative transfer simulations are used to demonstrate the advantage offered by this technique.

Radiative transfer in nonuniformly refracting layered media: atmosphere-ocean system.

We have applied the discrete-ordinate method to solve the radiative-transfer problem pertaining to a system consisting of two strata with different indices of refraction. The refraction and reflection at the interface are taken into account. The relevant changes (as compared with the standard problem with a constant index of refraction throughout the medium) in formulation and solution of the radiative-transfer equation, including the proper application of interface and boundary conditions, are described. Appropriate quadrature points (streams) and weights are chosen for the interface-continuity relations. Examples of radiative transfer in the coupled atmosphere-ocean system are provided. To take into account the region of total reflection in the ocean, additional angular quadrature points are required, compared with those used in the atmosphere and in the refractive region of the ocean that communicates directly with the atmosphere. To verify the model we have tested for energy conservation. We also discuss the effect of the number of streams assigned to the refractive region and the total reflecting region on the convergence. Our results show that the change in the index of refraction between the two strata significantly affects the radiation field. The radiative-transfer model we present is designed for application to the atmosphere-ocean system, but it can be applied to other systems that need to consider the change in the index of refraction between two strata.

Derivation of total ozone abundance and cloud effects from spectral irradiance measurements.

We describe a method to infer total ozone abundance and effective cloud transmission from global (diffuse plus direct) spectral irradiance measurements taken at the Earth's surface. The derivation of total ozone abundance relies on the comparison of measured irradiance ratios at two wavelengths in the UV part of the spectrum with a synthetic chart of this ratio computed for a variety of ozone abundances. One of these wavelengths should be appreciably absorbed by ozone (e.g., 305 nm) compared with the other one (e.g., 340 nm). This synthetic ratio (and therefore also the inferred total ozone abundance) is insensitive to the value of the surface albedo used in the model computations. Comparison with independent in situ and remote (from ground and space) determinations of total ozone abundance shows that measurements of global irradiances provide a reliable means of inferring the total column ozone amount for clear as well as cloudy sky conditions. Computer simulations are used to demonstrate that the ozone abunance inferred from global irradiance measurements is quite insensitive to cloud effects, whereas the use of the scattered irradiance only or the zenith sky intensity (measured routinely in the Dobson network on overcast days) requires substantial corrections for cloud effects. Effective cloud transmission is estimated from the data by comparing the measured irradiance at a wavelength where ozone absorption is minimal (e.g., 350 nm) to the clear-sky value. Irradiances generated by a plane-parallel radiation model as a function of cloud optical thickness are used to estimate an equivalent stratified cloud optical depth. These estimates of cloud transmission and optical depth are sensitive to ground reflection, implying that the accurate determination of cloud attenuation requires precise knowledge of the surface albedo.

Numerically stable algorithm for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered media.

We summarize an advanced, thoroughly documented, and quite general purpose discrete ordinate algorithm for time-independent transfer calculations in vertically inhomogeneous, nonisothermal, plane-parallel media. Atmospheric applications ranging from the UV to the radar region of the electromagnetic spectrum are possible. The physical processes included are thermal emission, scattering, absorption, and bidirectional reflection and emission at the lower boundary. The medium may be forced at the top boundary by parallel or diffuse radiation and by internal and boundary thermal sources as well. We provide a brief account of the theoretical basis as well as a discussion of the numerical implementation of the theory. The recent advances made by ourselves and our collaborators-advances in both formulation and numerical solution-are all incorporated in the algorithm. Prominent among these advances are the complete conquest of two illconditioning problems which afflicted all previous discrete ordinate implementations: (1) the computation of eigenvalues and eigenvectors and (2) the inversion of the matrix determining the constants of integration. Copies of the FORTRAN program on microcomputer diskettes are available for interested users.

Image reconstructions in diffraction tomography from limited transmitted field data sets.

Diffraction tomography reconstructions of objects from limited transmitted field data sets are discussed together with theoretical analyses and results of numerical experiments. It is shown that limited data sets, representing only a small part of the complete data sets, can be used for reconstructions in diffraction tomography with satisfactory accuracy. We also find that, in diffraction tomography based on the hybrid filtered backpropagation and the first-Rytov approximation, the use of limited data sets can provide a larger range of validity than the use of complete data sets, the reason being that limited data sets pose less-severe phase-unwrapping problems.

Ultraviolet radiation at high latitudes and the risk of skin cancer.

Solar irradiance in the spectral region 280 to 800 nm was measured with a double monochromator at 2 locations in Norway, Tromsø (69.7 degrees N) and Longyearbyen (78.2 degrees N). During the observational (midnight sun) period in Longyearbyen, the maximum UVB irradiance recorded was less than 0.3 W/m2, and no radiation was detected for wavelengths below 300 nm. Such low levels are believed to be a consequence of the low solar elevation angle and the high ozone content of the Arctic ozone layer, which absorbs the incident UV light. With levels between 280 and 350 DU over the period of study, Tromsø and Longyearbyen recorded only one-ninth of the calculated UVB radiation at the equator. There is therefore a considerably higher risk of radiation damage to the skin in equatorial regions (controlling for skin type), a finding that agrees with the statistical evidence for a 7-8 times higher rate of skin cancer in the white population of equatorial countries.

Application of the extended boundary condition method to particles with sharp edges: a comparison of two surface integration approaches.

We discuss, test, and compare two surface integration approaches that have been proposed for applying the extended boundary condition method (EBCM) to particles with sharp edges. One is based on approximating surface parameterization by a smooth function. By investigating the accuracy of this approach we find a quantitative condition for the radius of curvature of the approximate particle surface at the edge. The second approach is based on a special quadrature scheme for performing surface integration in the EBCM. For the simple test case of a cubic particle we find that the numerical advantages of the second method outweigh those of the first method, resulting in an overall reduction of computation time by a factor of 2. We conclude that the second method is preferable to the first when one is dealing with regularly shaped particles, for which the special quadrature scheme is reasonably simple to implement, and with particles with a relatively small number of sharp edges.

Application of the extended boundary condition method to homogeneous particles with point-group symmetries.

The numerical evaluation of surface integrals is the most time-consuming part of the extended boundary condition method (EBCM) for calculating the T matrix. An efficient implementation of the method is presented for homogeneous particles with discrete geometric symmetries and is applied to regular polyhedral prisms of finite length. For such prisms, an efficient quadrature scheme for computing the surface integrals is developed. Exploitation of these symmetries in conjunction with the new quadrature scheme leads to a reduction in CPU time by 3 orders of magnitude from that of a general EBCM implementation with no geometry-specific adaptations. The improved quadrature scheme and the exploitation of symmetries account for, respectively, 1 and 2 orders of magnitude in the total reduction of the CPU time. Test results for scattering by rectangular parallelepipeds and hexagonal plates are shown to agree well with corresponding results obtained by use of the discrete-dipole approximation. A model application for various polyhedral prisms is presented.

Validity of the diffusion approximation in bio-optical imaging.

Accurate numerical simulations based on rigorous radiative transfer theory are used to assess the validity of the diffusion approximation that is frequently used in bio-optical imaging. These simulations show that the error is large for a non-index-matched boundary between air and tissue. This weakness of the diffusion approximation underscores the need to understand how diffusion theory can be used to extract accurate values of tissue optical properties. A validity criterion for the diffusion approximation is established on the basis of the single-scattering albedo a and the asymmetry factor g for a slab with index-matched boundaries.