Anatomical 3D Modeling Using IR Sensors and Radiometric Processing Based on Structure from Motion: Towards a Tool for the Diabetic Foot Diagnosis.

Medical infrared thermography has proven to be a complementary procedure to physiological disorders, such as the diabetic foot. However, the technique remains essentially based on 2D images that display partial anatomy. In this context, a 3D thermal model provides improved visualization and faster inspection. This paper presents a 3D reconstruction method associated with temperature information. The proposed solution is based on a Structure from Motion and Multi-view Stereo approach, exploiting a set of multimodal merged images. The infrared images were obtained by automatically processing the radiometric data to remove thermal interferences, segment the RoI, enhance false-color contrast, and for multimodal co-registration under a controlled environment and a ∆T < 2.6% between the RoI and thermal interferences. The geometric verification accuracy was 77% ± 2%. Moreover, a normalized error was adjusted per sample based on a linear model to compensate for the curvature emissivity (error ≈ 10% near to 90°). The 3D models were displayed with temperature information and interaction controls to observe any point of view. The temperature sidebar values were assigned with information retrieved only from the RoI. The results have proven the feasibility of the 3D multimodal construction to be used as a promising tool in the diagnosis of diabetic foot.

Simulation of an optimized technique based on DS-CDMA for simultaneous transmission of multichannel biosignals.

Telemedicine is becoming increasingly, with applications in many areas of healthcare, such as home telecare of the elderly, diagnosis at a distance and robotic surgery. The simultaneous transmission of several leads of biomedical signals should be considered in telemedicine, given the many benefits it brings. Code division multiple access (CDMA) is a multiple access technique that enables users to transmit independent information simultaneously within the same bandwidth. The direct sequence CDMA (DS-CDMA) is a variant of the CDMA technique in which a pseudorandom sequence having a higher bandwidth than the information signal is used to modulate the information signal directly. Biomedical signals are confidential; thus, their transmission must be secured. In this paper we propose a protocol similar to DS-CDMA for the simultaneous transmission of all of the leads of some multichannel biomedical signals. We assigned orthogonal codes to different leads of a signal. The convolution of each lead with the code gives a signal spread over a broad frequency band. All of the spread signals are then mixed to produce a single composite signal. This composite signal is frequency modulated, amplified and transmitted. At the reception, inverse functions to the previous are developed to perform demodulation, demultiplexing and extraction of the physiological signals transmitted. We used the discrete Walsh functions as codes. The results obtained are satisfactory, even in situations where the noise disturbances are significant.

Qualitative estimation of pelvic organ interactions and their consequences on prostate motion: study on a deceased person.

In an attempt to have better targeting of the prostate during radiotherapy it is necessary to understand the mechanical interactions between bladder, rectum, and prostate and estimate their consequences on prostate motion. For this, the volumes of bladder, rectum, and lungs were modified concomitantly on a deceased person. A CT acquisition was performed for each of these different pelvic configurations (36 acquisitions). An increase in the volume of the bladder or lungs induces a compression of tissues of the pelvic area from its supero-anterior (S-A) to infero-posterior (I-P) side. Conversely, an increase of rectum volume induces a compression from the I-P to the S-A side of the pelvic region. These compressive actions can be added or subtracted from each other, depending on their amplitudes and directions. Prostate motion occurs when a movement of the rectum is observed (this movement depends, itself, on lungs and bladder volume). The maximum movement of prostate is 9 mm considering maximal bladder or rectal action, and 11 mm considering maximum lung action. In some other cases, opposition of compressive effects can lead to stasis of the prostate. Based on the volumes of bladder, rectum, and lungs, it is possible to qualitatively estimate the movement of organs of the pelvic area. The best way to reduce prostate movement is to recommend the patient to have an empty rectum, with either full bladder and/or full lungs.

3-D reconstruction of microcalcification clusters using stereo imaging: algorithm and mammographic unit calibration.

The three-dimensional (3-D) shape of microcalcification clusters is an important indicator in early breast cancer detection. In fact, there is a relationship between the cluster topology and the type of lesion (malignant or benign). This paper presents a 3-D reconstruction method for such clusters using two 2-D views acquired during standard mammographic examinations. For this purpose, the mammographic unit was modeled using a camera with virtual optics. This model was used to calibrate the acquisition unit and then to reconstruct the clusters in the 3-D space after microcalcification segmentation and matching. The proposed model is hardware independent since it is suitable for digital mammographic units with different geometries and with various physical acquisition principles. Three-dimensional reconstruction results are presented here to prove the validity of the method. Tests were first performed using a phantom with a well-known geometry. The latter contained X-ray opaque glass balls representing microcalcifications. The positions of these balls were reconstructed with a 16.25-microm mean accuracy. This very high inherent algorithm accuracy is more than enough for a precise 3-D cluster representation. Further validation tests were carried out using a second phantom including a spherical cluster. This phantom was built with materials simulating the behavior of both mammary tissue and microcalcifications toward Xrays. The reconstructed shape was effectively spherical. Finally, reconstructions were carried out for real clusters and their results are also presented.

Correspondences between microcalcification projections on two mammographic views acquired with digital systems.

In this paper, we have proposed an algorithm for automatic matching of MC projections viewed on two mammograms of the same breast. The implemented algorithm consists in three steps. From five morphological features of the MC, a similarity function was built between each MC of the first image of the pair, and each MC of the second image. These values quantified the resemblance between each pair of MC and permitted, for a given MC of the first image, to sort the MC of the second image which could be matched to it. From the geometry of the system providing the pairs of images being analysed, we derived some geometrical constraints that must be satisfied by corresponding MC. In order to take into account the fact that due to breast deformation, the corresponding MC is often off the classical epipolar line constructed on the basis of stereovision assumption, we instead consider an epipolar strip on both side of the approximated epipolar line. The MC of the second image which was out of that strip was eliminated. Then a coefficient was applied to the remaining MC that took into account their distances to the approximated epipolar line. Finally, the selection procedure was used to pick out the right pair. In order to test our algorithm, we compared the result it yielded with those coming from two operators who matched a number of MC with confidence. A concordance of 78.43% was obtained between confident manual matching and automatic matching. Since this algorithm was designed to be part of a tool for 3D reconstruction of microcalcification clusters, we reconstructed some clusters with manual matching and automatic matching and compared the shapes of the clusters obtained in these two ways. The resemblance in some cases was very good and average in a number of others. This suggests that our algorithm should be improved. Therefore, apart from ongoing effort to introduce other constraints, we believe that taking into account the third view provided by the imaging system could be of great help. We shall soon explore this possibility. Overall, we believe that despite the failure of our tool in some cases it can already at this stage be used with some confidence.

Assessment of a Monte-Carlo simulation of SPECT recordings from a new-generation heart-centric semiconductor camera: from point sources to human images.

Geant4 application for tomographic emission (GATE), a Monte-Carlo simulation platform, has previously been used for optimizing tomoscintigraphic images recorded with scintillation Anger cameras but not with the new-generation heart-centric cadmium-zinc-telluride (CZT) cameras. Using the GATE platform, this study aimed at simulating the SPECT recordings from one of these new CZT cameras and to assess this simulation by direct comparison between simulated and actual recorded data, ranging from point sources to human images. Geometry and movement of detectors, as well as their respective energy responses, were modeled for the CZT 'D.SPECT' camera in the GATE platform. Both simulated and actual recorded data were obtained from: (1) point and linear sources of (99m)Tc for compared assessments of detection sensitivity and spatial resolution, (2) a cardiac insert filled with a (99m)Tc solution for compared assessments of contrast-to-noise ratio and sharpness of myocardial borders and (3) in a patient with myocardial infarction using segmented cardiac magnetic resonance imaging images. Most of the data from the simulated images exhibited high concordance with the results of actual images with relative differences of only: (1) 0.5% for detection sensitivity, (2) 6.7% for spatial resolution, (3) 2.6% for contrast-to-noise ratio and 5.0% for sharpness index on the cardiac insert placed in a diffusing environment. There was also good concordance between actual and simulated gated-SPECT patient images for the delineation of the myocardial infarction area, although the quality of the simulated images was clearly superior with increases around 50% for both contrast-to-noise ratio and sharpness index. SPECT recordings from a new heart-centric CZT camera can be simulated with the GATE software with high concordance relative to the actual physical properties of this camera. These simulations may be conducted up to the stage of human SPECT-images even if further refinement is needed in this setting.

Validation of a digital mammographic unit model for an objective and highly automated clinical image quality assessment.

In mammography, image quality assessment has to be directly related to breast cancer indicator (e.g. microcalcifications) detectability. Recently, we proposed an X-ray source/digital detector (XRS/DD) model leading to such an assessment. This model simulates very realistic contrast-detail phantom (CDMAM) images leading to gold disc (representing microcalcifications) detectability thresholds that are very close to those of real images taken under the simulated acquisition conditions. The detection step was performed with a mathematical observer. The aim of this contribution is to include human observers into the disc detection process in real and virtual images to validate the simulation framework based on the XRS/DD model. Mathematical criteria (contrast-detail curves, image quality factor, etc.) are used to assess and to compare, from the statistical point of view, the cancer indicator detectability in real and virtual images. The quantitative results given in this paper show that the images simulated by the XRS/DD model are useful for image quality assessment in the case of all studied exposure conditions using either human or automated scoring. Also, this paper confirms that with the XRS/DD model the image quality assessment can be automated and the whole time of the procedure can be drastically reduced. Compared to standard quality assessment methods, the number of images to be acquired is divided by a factor of eight.

Compared performance of high-sensitivity cameras dedicated to myocardial perfusion SPECT: a comprehensive analysis of phantom and human images.

UNLABELLED: Differences in the performance of cadmium-zinc-telluride (CZT) cameras or collimation systems that have recently been commercialized for myocardial SPECT remain unclear. In the present study, the performance of 3 of these systems was compared by a comprehensive analysis of phantom and human SPECT images. METHODS: We evaluated the Discovery NM 530c and DSPECT CZT cameras, as well as the Symbia Anger camera equipped with an astigmatic (IQ x SPECT) or parallel-hole (conventional SPECT) collimator. Physical performance was compared on reconstructed SPECT images from a phantom and from comparable groups of healthy subjects. RESULTS: Classifications were as follows, in order of performance. For count sensitivity on cardiac phantom images (counts x s(-1) x MBq(-1)), DSPECT had a sensitivity of 850; Discovery NM 530c, 460; IQ x SPECT, 390; and conventional SPECT, 130. This classification was similar to that of myocardial counts normalized to injected activities from human images (respective mean values, in counts x s(-1) x MBq(-1): 11.4 ± 2.6, 5.6 ± 1.4, 2.7 ± 0.7, and 0.6 ± 0.1). For central spatial resolution: Discovery NM 530c was 6.7 mm; DSPECT, 8.6 mm; IQ x SPECT, 15.0 mm; and conventional SPECT, 15.3 mm, also in accordance with the analysis of the sharpness of myocardial contours on human images (in cm(-1): 1.02 ± 0.17, 0.92 ± 0.11, 0.64 ± 0.12, and 0.65 ± 0.06, respectively). For contrast-to-noise ratio on the phantom: Discovery NM 530c had a ratio of 4.6; DSPECT, 4.1; IQ x SPECT, 3.9; and conventional SPECT, 3.5, similar to ratios documented on human images (5.2 ± 1.0, 4.5 ± 0.5, 3.9 ± 0.6, and 3.4 ± 0.3, respectively). CONCLUSION: The performance of CZT cameras is dramatically higher than that of Anger cameras, even for human SPECT images. However, CZT cameras differ in that spatial resolution and contrast-to-noise ratio are better with the Discovery NM 530c, whereas count sensitivity is markedly higher with the DSPECT.

Voxel-based quantitative analysis of brain images from ¹8F-FDG PET with a block-matching algorithm for spatial normalization.

OBJECTIVE: Statistical Parametric Mapping (SPM) is widely used for the quantitative analysis of brain images from ¹8F fluorodeoxyglucose positron emission tomography (FDG PET). SPM requires an initial step of spatial normalization to align all images to a standard anatomic model (the template), but this may lead to image distortion and artifacts, especially in cases of marked brain abnormalities. This study aimed at assessing a block-matching (BM) normalization algorithm, where most transformations are not directly computed on the overall brain volume but through small blocks, a principle that is likely to minimize artifacts. METHODS: Large and/or small hypometabolic areas were artificially simulated in initially normal FDG PET images to compare the results provided by statistical tests computed after either SPM or BM normalization. RESULTS: Results were enhanced by BM, compared with SPM, with regard to (i) errors in the estimation of large defects volumes (about 2-fold lower) because of a lower image distortion, and (ii) rates of false-positive foci when numerous or extended abnormalities were simulated. These observations were strengthened by analyses of FDG PET examinations from epileptic patients. CONCLUSIONS: Results obtained with the BM normalization of brain FDG PET appear more precise and robust than with SPM normalization, especially in cases of numerous or extended abnormalities.

Computation of realistic virtual phantom images for an objective lesion detectability assessment in digital mammography.

Image quality assessment is required for an optimal use of mammographic units. On the one hand, there are objective image quality assessment methods based on the measurement of technical parameters such as modulation transfer function (MTF), noise power spectrum (NPS) or detection quantum efficiency (DQE) describing performances of digital detectors. These parameters are, however, without direct relationship with lesion detectability in clinical practice. On the other hand, there are image quality assessment methods involving time consuming procedures, but presenting a direct relationship with lesion detectability. This contribution describes an X-ray source/digital detector model leading to the simulation of virtual contrast-detail phantom (CDMAM) images. The virtual image computation method requires the acquisition of only few real images and allows for an objective image quality assessment presenting a direct relationship with lesion detectability. The transfer function of the proposed model takes as input physical parameters (MTF* and noise) measured under clinical conditions on mammographic units. As presented in this contribution, MTF* is a modified MTF taking into account the effects due to X-ray scatter in the breast and magnification. Results obtained with the structural similarity index prove that the simulated images are quite realistic in terms of contrast and noise. Tests using contrast detail curves highlight the fact that the simulated and real images lead to very similar data quality in terms of lesion detectability. Finally, various statistical tests show that quality factors computed for both the simulated images and the real images are very close for the two data sets.

Digestive activity evaluation by multichannel abdominal sounds analysis.

This paper introduces a complete methodology for abdominal sounds analysis, from signal acquisition to statistical data analysis. The goal is to evaluate if and how phonoenterograms can be used to detect different functioning modes of the normal gastrointestinal tract, both in terms of localization and of time evolution during the digestion. After the description of the acquisition protocol and the employed instrumentation, several signal processing steps are presented: wavelet denoising and segmentation, artifact suppression, and source localization. Next, several physiological features are extracted from the processed signals issued from a database of 14 healthy volunteers, recorded during 3 h after a standardized meal. Data analysis is performed using a multifactorial statistical method. Based on the introduced approach, we show that the abdominal regions of healthy volunteers present statistically significant phonoenterographic characteristics, which evolve differently during the normal digestion. The most significant feature allowing us to distinguish regions and time differences is the number of recorded sounds, but important information is also carried by sound amplitudes, frequencies, and durations. Depending on the considered feature, the sounds produced by different abdominal regions (especially stomach, ileocaecal, and lower abdomen regions) present a specific distribution over space and time. This information, statistically validated, is usable in further studies as a comparison term with other normal or pathological conditions.

Mosaicing of bladder endoscopic image sequences: distortion calibration and registration algorithm.

Cancers located on the internal wall of bladders can be detected in image sequences acquired with endoscopes. The clinical diagnosis and follow-up can be facilitated by building a unique panoramic image of the bladder with the images acquired from different viewpoints. This process, called image mosaicing, consists of two steps. In the first step, consecutive images are pairwise registered to find the local transformation matrices linking geometrically consecutive images. In the second step, all images are placed in a common and global coordinate system. In this contribution, a mutual information-based similarity measure and a stochastic gradient optimization method were implemented in the registration process. However, the images have to be preprocessed in order to register the data in a robust way. Thus, a simple correction method of the distortions affecting endoscopic images is presented. After the placement of all images in the global coordinate system, the parameters of the local transformation matrices are all adjusted to improve the visual aspect of the panoramic images. Phantoms are used to evaluate the global mosaicing accuracy and the limits of the registration algorithm. The mean distances between ground truth positions in the mosaiced image range typically in 1-3 pixels. Results given for in vivo patient data illustrate the ability of the algorithm to give coherent panoramic images in the case of bladders.

ECG data compression using Jacobi polynomials.

Data compression is a frequent signal processing operation applied to ECG. We present here a method of ECG data compression utilizing Jacobi polynomials. ECG signals are first divided into blocks that match with cardiac cycles before being decomposed in Jacobi polynomials bases. Gauss quadratures mechanism for numerical integration is used to compute Jacobi transforms coefficients. Coefficients of small values are discarded in the reconstruction stage. For experimental purposes, we chose height families of Jacobi polynomials. Various segmentation approaches were considered. We elaborated an efficient strategy to cancel boundary effects. We obtained interesting results compared with ECG compression by wavelet decomposition methods. Some propositions are suggested to improve the results.