Usually, image compression techniques that use only one transform exhibit some poor properties. For instance, the Discrete Cosine Transform (DCT) cannot efficiently represent high frequency components, resulting in blurred images. The Multiple Bases Representation (MBR) compression technique, which uses two or more transforms, is found to be superior to the single transform techniques in terms of representation efficiency. However, some bits in the MBR representation are needed to track the basis information. The MBR image quality is deteriorated by discontinuities at block boundaries, as is the standard DCT transform. In this thesis, test images are distorted by MBR compression using a Recursive Residual Projection algorithm. This algorithm is a sub-optimal method to find the best basis vector subset for representing images based on multiple orthogonal bases. The MBR distorted images are reconstructed by the iterative method of Projection onto Convex Sets (POCS). Many constraints that form convex sets are reviewed and examined. Due to the high distortion at the block boundaries, some constraints are introduced particularly to reduce artifacts at the boundaries. Some constraints add energy to the reconstructed images while others remove energy. Thus, the initial vectors play a key role in the performance of the POCS method for better MBR reconstruction. This thesis also determines the most appropriate initial vector for each constraint. Finally, the composite projections associated with the sign, minimum decreasing and norm-of-slope constraints are used to improve the reconstruction of the MBR distorted images and the effect of ordering of the projections is investigated. / Master of Science

No description available.

L’approximation par harmoniques sphériques (SPN) simplifiées de l’équation de transfert radiatif a été proposée comme un modèle fiable de propagation de la lumière dans les tissus biologiques. Cependant, peu de solutions analytiques ont été trouvées pour ce modèle. De telles solutions analytiques sont d’une grande valeur pour valider les solutions numériques des équations SPN, auxquelles il faut recourir dans le cas de tissus avec des géométries courbes complexes. Dans la première partie de cette thèse, des solutions analytiques pour deux géométries courbes sont présentées pour la première fois, à savoir pour la sphère et pour le cylindre. Pour les deux solutions, les conditions aux frontières générales tenant compte du saut d’indice de réfraction à l’interface du tissus et de son milieu environnant, telles qu’applicables à l’optique biomédicale, sont utilisées. Ces solutions sont validées à l’aide de simulations Monte Carlo basées sur un maillage de discrétisation du milieu. Ainsi, ces solutions permettent de valider rapidement un code numérique, par exemple utilisant les différences finies ou les éléments finis, sans nécessiter de longues simulations Monte Carlo. Dans la deuxième partie de cette thèse, la reconstruction itérative pour l’imagerie par tomographie optique diffuse par fluorescence est proposée sur la base d’une fonction objective et de son terme de régularisation de type Lq-Lp. Pour résoudre le problème inverse d’imagerie, la discrétisation du modèle de propagation de la lumière est effectuée en utilisant la méthode des différences finies. La reconstruction est effectuée sur un modèle de souris numérique en utilisant un maillage multi-échelle. Le problème inverse est résolu itérativement en utilisant une méthode d’optimisation. Pour cela, le gradient de la fonction de coût par rapport à la carte de concentration de l’agent fluorescent est nécessaire. Ce gradient est calculé à l’aide d’une méthode adjointe. Des mesures quantitatives utilisées en l’imagerie médicale sont utilisées pour évaluer la performance de l’approche de reconstruction dans différentes conditions. L’approche Lq-Lp montre des performances quantifiées élevées par rapport aux algorithmes traditionnels basés sur des fonction coût de type somme de carrés de différences. / Abstract : The simplified spherical harmonics (SPN) approximation to the radiative transfer equation has been proposed as a reliable model of light propagation in biological tissues. However, few analytical solutions have been found for this model. Such analytical solutions are of great value to validate numerical solutions of the SPN equations, which must be resorted to when dealing with media with complex curved geometries. In the first part of this thesis, analytical solutions for two curved geometries are presented for the first time, namely for the sphere and for the cylinder. For both solutions, the general refractiveindex mismatch boundary conditions, as applicable in biomedical optics, are resorted to. These solutions are validated using mesh-based Monte Carlo simulations. So validated, these solutions allow in turn to rapidly validate numerical code, based for example on finite differences or on finite elements, without requiring lengthy Monte Carlo simulations. provide reliable tool for validating numerical simulations. In the second part, iterative reconstruction for fluorescence diffuse optical tomography imaging is proposed based on an Lq-Lp framework for formulating an objective function and its regularization term. To solve the imaging inverse problem, the discretization of the light propagation model is performed using the finite difference method. The framework is used along with a multigrid mesh on a digital mouse model. The inverse problem is solved iteratively using an optimization method. For this, the gradient of the cost function with respect to the fluorescent agent’s concentration map is necessary. This is calculated using an adjoint method. Quantitative metrics resorted to in medical imaging are used to evaluate the performance of the framework under different conditions. The results obtained support this new approach based on an Lq-Lp formulation of cost functions in order to solve the inverse fluorescence problem with high quantified performance.

This thesis reports the work done towards design and fabrication of a versatile and low cost, frequency domain DOT (Diffuse Optical Tomography) Imager. A design which uses only a single fiber for the source and a single fiber bundle for the detector is reported. From near the source, to diametrically opposite to the source, the detected intensity of scattered light varies by three to four orders in magnitude, depending on the tissue/phantom absorption and scattering properties. The photo multiplier tube’s (PMT’s) gain is controlled to operate it in the linear range, thus increasing the dynamic range of detection. Increasing the dynamic range by multi channel data acquisition is also presented. Arresting the oscillations of a stepper using a negative torque braking method is also adopted in this application for increasing the speed of data acquisition. The finite element method (FEM) for obtaining photon density solution to the transport equation and the model based iterative image reconstruction (MPBIIR) algorithm are developed for verifying the experimental prototype. Simulation studies presented towards the end of this thesis work provide insight into the nature of measurements. The optical absorption reconstructed images from the simulation, verified the validity of implementation of the reconstruction method for further reconstructions from data gathered from the developed imager. A single iteration of MOBIIR to segment the region of interest (ROI) using an homogeneous measurement estimate is presented. Using the single iteration MOBIIR to obtain a relatively more accurate starting value for the optical absorption coefficient, and the reconstruction results for data obtained from tissue mimicking solid epoxy-resin phantom with a single in-homogeneity inclusion is also presented to demonstrate the imager prototype.