| 1. |
- Borga, Magnus, 1965-
(författare)
-
Learning Multidimensional Signal Processing
- 1998
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- The subject of this dissertation is to show how learning can be used for multidimensional signal processing, in particular computer vision. Learning is a wide concept, but it can generally be defined as a system’s change of behaviour in order to improve its performance in some sense.Learning systems can be divided into three classes: supervised learning, reinforcement learning and unsupervised learning. Supervised learning requires a set of training data with correct answers and can be seen as a kind of function approximation. A reinforcement learning system does not require a set of answers. It learns by maximizing a scalar feedback signal indicating the system’s performance. Unsupervised learning can be seen as a way of finding a good representation of the input signals according to a given criterion.In learning and signal processing, the choice of signal representation is a central issue. For high-dimensional signals, dimensionality reduction is often necessary. It is then important not to discard useful information. For this reason, learning methods based on maximizing mutual information are particularly interesting.A properly chosen data representation allows local linear models to be used in learning systems. Such models have the advantage of having a small number of parameters and can for this reason be estimated by using relatively few samples. An interesting method that can be used to estimate local linear models is canonical correlation analysis (CCA). CCA is strongly related to mutual information. The relation between CCA and three other linear methods is discussed. These methods are principal component analysis (PCA), partial least squares (PLS) and multivariate linear regression (MLR). An iterative method for CCA, PCA, PLS and MLR, in particular low-rank versions of these methods, is presented.A novel method for learning filters for multidimensional signal processing using CCA is presented. By showing the system signals in pairs, the filters can be adapted to detect certain features and to be invariant to others. A new method for local orientation estimation has been developed using this principle. This method is significantly less sensitive to noise than previously used methods.Finally, a novel stereo algorithm is presented. This algorithm uses CCA and phase analysis to detect the disparity in stereo images. The algorithm adapts filters in each local neighbourhood of the image in a way which maximizes the correlation between the filtered images. The adapted filters are then analysed to find the disparity. This is done by a simple phase analysis of the scalar product of the filters. The algorithm can even handle cases where the images have different scales. The algorithm can also handle depth discontinuities and give multiple depth estimates for semi-transparent images.
|
|
| 2. |
- Brun, Anders, 1976-
(författare)
-
Manifolds in Image Science and Visualization
- 2007
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- A Riemannian manifold is a mathematical concept that generalizes curved surfaces to higher dimensions, giving a precise meaning to concepts like angle, length, area, volume and curvature. A glimpse of the consequences of a non-flat geometry is given on the sphere, where the shortest path between two points – a geodesic – is along a great circle. Different from Euclidean space, the angle sum of geodesic triangles on the sphere is always larger than 180 degrees.Signals and data found in applied research are sometimes naturally described by such curved spaces. This dissertation presents basic research and tools for the analysis, processing and visualization of such manifold-valued data, with a particular emphasis on future applications in medical imaging and visualization.Two-dimensional manifolds, i.e. surfaces, enter naturally into the geometric modelling of anatomical entities, such as the human brain cortex and the colon. In advanced algorithms for processing of images obtained from computed tomography (CT) and ultrasound imaging (US), images themselves and derived local structure tensor fields may be interpreted as two- or three-dimensional manifolds. In diffusion tensor magnetic resonance imaging (DT-MRI), the natural description of diffusion in the human body is a second-order tensor field, which can be related to the metric of a manifold. A final example is the analysis of shape variations of anatomical entities, e.g. the lateral ventricles in the brain, within a population by describing the set of all possible shapes as a manifold.Work presented in this dissertation include: Probabilistic interpretation of intrinsic and extrinsic means in manifolds. A Bayesian approach to filtering of vector data, removing noise from sampled manifolds and signals. Principles for the storage of tensor field data and learning a natural metric for empirical data.The main contribution is a novel class of algorithms called LogMaps, for the numerical estimation of logp (x) from empirical data sampled from a low-dimensional manifold or geometric model embedded in Euclidean space. The logp (x) function has been used extensively in the literature for processing data in manifolds, including applications in medical imaging such as shape analysis. However, previous approaches have been limited to manifolds where closed form expressions of logp (x) have been known. The introduction of the LogMap framework allows for a generalization of the previous methods. The application of LogMaps to texture mapping, tensor field visualization, medial locus estimation and exploratory data analysis is also presented.
|
|
| 3. |
- Cros, Olivier
(författare)
-
Structural properties of the mastoid using image analysis and visualization
- 2017
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- The mastoid, located in the temporal bone, houses an air cell system whose cells have a variation in size that can go far below current conventional clinical CT scanner resolution. Therefore, the mastoid air cell system is only partially represented in a CT scan. Where the conventional clinical CT scanner lacks level of minute details, micro-CT scanning provides an overwhelming amount of ne details. The temporal bone being one of the most complex in the human body, visualization of micro-CT scanning of this boneawakens the curiosity of the experimenter, especially with the correct visualization settings.This thesis first presents a statistical analysis determining the surface area to volume ratio of the mastoid air cell system of human temporal bone, from micro-CT scanning using methods previously applied for conventional clinical CT scans. The study compared current results with previous studies, with successive downsampling the data down to a resolution found in conventional clinical CT scanning. The results from the statistical analysis showed that all the small mastoid air cells, that cannot be detected in conventional clinical CT scans, do heavily contribute to the estimation of the surface area, and in consequence to the estimation of the surface area to volume ratio by a factor of about 2.6. Such a result further strengthens the idea of the mastoid to play an active role in pressure regulation and gas exchange.Discovery of micro-channels through specific use of a non-traditional transfer function was then reported, where a qualitative and a quantitative pre-analysis were performed and reported. To gain more knowledge about these micro-channels, a local structure tensor analysis was applied where structures are described in terms of planar, tubular, or isotropic structures. The results from this structural tensor analysis suggest these microchannels to potentially be part of a more complex framework, which hypothetically would provide a separate blood supply for the mucosa lining the mastoid air cell system.The knowledge gained from analysing the micro-channels as locally providing blood to the mucosa, led to the consideration of how inflammation of the mucosa could impact the pneumatization of the mastoid air cell system. Though very primitive, a 3D shape analysis of the mastoid air cell system was carried out. The mastoid air cell system was first represented in a compact form through a medial axis, from which medial balls could be used. The medial balls, representative of how large the mastoid air cells can be locally, were used in two complementary clustering methods, one based on the size diameter of the medial balls and one based on their location within the mastoid air cell system. From both quantitative and qualitative statistics, it was possible to map the clusters based on pre-defined regions already described in the literature, which opened the door for new hypotheses concerning the effect of mucosal inflammation on the mastoid pneumatization.Last but not least, discovery of other structures, previously unreported in the literature, were also visually observed and briefly discussed in this thesis. Further analysis of these unknown structures is needed.
|
|
| 4. |
- Eklund, Anders, 1981-
(författare)
-
Computational Medical Image Analysis : With a Focus on Real-Time fMRI and Non-Parametric Statistics
- 2012
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Functional magnetic resonance imaging (fMRI) is a prime example of multi-disciplinary research. Without the beautiful physics of MRI, there wouldnot be any images to look at in the first place. To obtain images of goodquality, it is necessary to fully understand the concepts of the frequencydomain. The analysis of fMRI data requires understanding of signal pro-cessing, statistics and knowledge about the anatomy and function of thehuman brain. The resulting brain activity maps are used by physicians,neurologists, psychologists and behaviourists, in order to plan surgery andto increase their understanding of how the brain works.This thesis presents methods for real-time fMRI and non-parametric fMRIanalysis. Real-time fMRI places high demands on the signal processing,as all the calculations have to be made in real-time in complex situations.Real-time fMRI can, for example, be used for interactive brain mapping.Another possibility is to change the stimulus that is given to the subject, inreal-time, such that the brain and the computer can work together to solvea given task, yielding a brain computer interface (BCI). Non-parametricfMRI analysis, for example, concerns the problem of calculating signifi-cance thresholds and p-values for test statistics without a parametric nulldistribution.Two BCIs are presented in this thesis. In the first BCI, the subject wasable to balance a virtual inverted pendulum by thinking of activating theleft or right hand or resting. In the second BCI, the subject in the MRscanner was able to communicate with a person outside the MR scanner,through a virtual keyboard.A graphics processing unit (GPU) implementation of a random permuta-tion test for single subject fMRI analysis is also presented. The randompermutation test is used to calculate significance thresholds and p-values forfMRI analysis by canonical correlation analysis (CCA), and to investigatethe correctness of standard parametric approaches. The random permuta-tion test was verified by using 10 000 noise datasets and 1484 resting statefMRI datasets. The random permutation test is also used for a non-localCCA approach to fMRI analysis.
|
|
| 5. |
- Knutsson, Hans
(författare)
-
Filtering and reconstruction in image processing
- 1982
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Image processing is a broad field posing a wide range of problems. The Work presented in this dissertation is mainly concerned with filter design subject to different criteria and constraints.The first part describes the development of a new radiographic reconstruction method designated Ectomography. The method is novel in that it allows reconstruction of an arbitrarily thick layer of an object using limited viewing angle.The subject of the second part is estimation and filtering of local image information. Quadrature filters are designed enabling accurate orientation and frequency estimates. The extracted information is shown to provide a good basis for efficient image enhancement and coding procedures.
|
|
| 6. |
- Landelius, Tomas
(författare)
-
Reinforcement Learning and Distributed Local Model Synthesis
- 1997
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Reinforcement learning is a general and powerful way to formulate complex learning problems and acquire good system behaviour. The goal of a reinforcement learning system is to maximize a long term sum of instantaneous rewards provided by a teacher. In its extremum form, reinforcement learning only requires that the teacher can provide a measure of success. This formulation does not require a training set with correct responses, and allows the system to become better than its teacher.In reinforcement learning much of the burden is moved from the teacher to the training algorithm. The exact and general algorithms that exist for these problems are based on dynamic programming (DP), and have a computational complexity that grows exponentially with the dimensionality of the state space. These algorithms can only be applied to real world problems if an efficient encoding of the state space can be found.To cope with these problems, heuristic algorithms and function approximation need to be incorporated. In this thesis it is argued that local models have the potential to help solving problems in high-dimensional spaces and that global models have not. This is motivated with the biasvariance dilemma, which is resolved with the assumption that the system is constrained to live on a low-dimensional manifold in the space of inputs and outputs. This observation leads to the introduction of bias in terms of continuity and locality.A linear approximation of the system dynamics and a quadratic function describing the long term reward are suggested to constitute a suitable local model. For problems involving one such model, i.e. linear quadratic regulation problems, novel convergence proofs for heuristic DP algorithms are presented. This is one of few available convergence proofs for reinforcement learning in continuous state spaces.Reinforcement learning is closely related to optimal control, where local models are commonly used. Relations to present methods are investigated, e.g. adaptive control, gain scheduling, fuzzy control, and jump linear systems. Ideas from these areas are compiled in a synergistic way to produce a new algorithm for heuristic dynamic programming where function parameters and locality, expressed as model applicability, are learned on-line. Both top-down and bottom-up versions are presented.The emerging local models and their applicability need to be memorized by the learning system. The binary tree is put forward as a suitable data structure for on-line storage and retrieval of these functions.
|
|
| 7. |
- Lundström, Claes, 1973-
(författare)
-
Efficient Medical Volume Visualization : An Approach Based on Domain Knowledge
- 2007
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Direct Volume Rendering (DVR) is a visualization technique that has proved to be a very powerful tool in many scientific visualization applications. Diagnostic medical imaging is one domain where DVR could provide clear benefits in terms of unprecedented possibilities for analysis of complex cases and highly efficient work flow for certain routine examinations. The full potential of DVR in the clinical environment has not been reached, however, primarily due to limitations in conventional DVR methods and tools.This thesis presents methods addressing four major challenges for DVR in clinical use. The foundation of all methods is to incorporate the domain knowledge of the medical professional in the technical solutions. The first challenge is the very large data sets routinely produced in medical imaging today. To this end a multiresolution DVR pipeline is proposed, which dynamically prioritizes data according to the actual impact in the rendered image to be reviewed. Using this prioritization the system can reduce the data requirements throughout the pipeline and provide high performance and visual quality in any environment.Another problem addressed is how to achieve simple yet powerful interactive tissue classification in DVR. The methods presented define additional attributes that effectively captures readily available medical knowledge. The task of tissue detection is also important to solve in order to improve efficiency and consistency of diagnostic image review. Histogram-based techniques that exploit spatial relations in the data to achieve accurate and robust tissue detection are presented in this thesis.The final challenge is uncertainty visualization, which is very pertinent in clinical work for patient safety reasons. An animation method has been developed that automatically conveys feasible alternative renderings. The basis of this method is a probabilistic interpretation of the visualization parameters.Several clinically relevant evaluations of the developed techniques have been performed demonstrating their usefulness. Although there is a clear focus on DVR and medical imaging, most of the methods provide similar benefits also for other visualization techniques and application domains.
|
|
| 8. |
- Rydell, Joakim, 1979-
(författare)
-
Advanced MRI Data Processing
- 2007
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Magnetic resonance imaging (MRI) is a very versatile imaging modality which can be used to acquire several different types of images. Some examples include anatomical images, images showing local brain activation and images depicting different types of pathologies. Brain activation is detected by means of functional magnetic resonance imaging (fMRI). This is useful e.g. in planning of neurosurgical procedures and in neurological research. To find the activated regions, a sequence of images of the brain is collected while a patient or subject alters between resting and performing a task. The variations in image intensity over time are then compared to a model of the variations expected to be found in active parts of the brain. Locations with high correlation between the intensity variations and the model are considered to be activated by the task.Since the images are very noisy, spatial filtering is needed before the activation can be detected. If adaptive filtering is used, i.e. if the filter at each location is adapted to the local neighborhood, very good detection performance can be obtained. This thesis presents two methods for adaptive spatial filtering of fMRI data. One of these is a modification of a previously proposed method, which at each position maximizes the similarity between the filter response and the model. A novel feature of the presented method is rotational invariance, i.e. equal sensitivity to activated regions in different orientations. The other method is based on bilateral filtering. At each position, this method averages pixels which are located in the same type of brain tissue and have similar intensity variation over time.A method for robust correlation estimation is also presented. This method automatically detects local bursts of noise in a signal and disregards the corresponding signal segments when the correlation is estimated. Hence, the correlation estimate is not affected by the noise bursts. This method is useful not only in analysis of fMRI data, but also in other applications where correlation is used to determine the similarity between signals.Finally, a method for correcting artifacts in complex MR images is presented. Complex images are used e.g. in the Dixon technique for separate imaging of water and fat. The phase of these images is often affected by artifacts and therefore need correction before the actual water and fat images can be calculated. The presented method for phase correction is based on an image integration technique known as the inverse gradient. The method is shown to provide good results even when applied to images with severe artifacts.
|
|
| 9. |
- Sigfridsson, Andreas, 1978-
(författare)
-
Multidimensional MRI of Myocardial Dynamics : Acquisition, Reconstruction and Visualization
- 2009
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Methods for measuring deformation and motion of the human heart in-vivo are crucial in the assessment of cardiac function. Applications ranging from basic physiological research, through early detection of disease to follow-up studies, all rely on the quality of the measurements of heart dynamics. This thesis presents new improved magnetic resonance imaging methods for acquisition, image reconstruction and visualization of cardiac motion and deformation.As the heart moves and changes shape during the acquisition, synchronization to the heart dynamics is necessary. Here, a method to resolve not only the cardiac cycle but also the respiratory cycle is presented. Combined with volumetric imaging, this produces a five-dimensional data set with two cyclic temporal dimensions. This type of data reveals unique physiological information, such as interventricular coupling in the heart in different phases of the respiratory cycle.The acquisition can also be sensitized to motion, measuring not only the magnitude of the magnetization but also a signal proportional to local velocity or displacement. This allows for quantification of the motion which is especially suitable for functional study of the cardiac deformation. In this work, an evaluation of the influence of several factors on the signal-to-noise ratio is presented for in-vivo displacement encoded imaging. Additionally, an extension of the method to acquire multiple displacement encoded slices in a single breath hold is also presented.Magnetic resonance imaging is usually associated with long scan times, and many methods exist to shorten the acquisition time while maintaining acceptable image quality. One class of such methods involves acquiring only a sparse subset of k-space. A special reconstruction is then necessary in order to obtain an artifact-free image. One family of these reconstruction techniques tailored for dynamic imaging is the k-t BLAST approach, which incorporates data-driven prior knowledge to suppress aliasing artifacts that otherwise occur with the sparse sampling. In this work, an extension of the original k-t BLAST method to two temporal dimensions is presented and applied to data acquired with full coverage of the cardio-respiratory cycles. Using this technique, termed k-t2 BLAST, simultaneous reduction of scan time and improved spatial resolution is demonstrated. Further, the loss of temporal fidelity when using the k-t BLAST approach is investigated, and an improved reconstruction is proposed for the application of cardiac function analysis.Visualization is a crucial part of the imaging chain. Scalar data, such as regular anatomical images, are straightforward to display. Myocardial strain and strain-rate, however, are tensor quantities which do not lend themselves to direct visualization. The problem of visualizing the tensor field is approached in this work by combining a local visualization that displays all degrees of freedom for a single tensor with an overview visualization using a scalar field representation of the complete tensor field. The scalar field is obtained by iterated adaptive filtering of a noise field, creating a continuous geometrical representation of the myocardial strain-rate tensor field.The results of the work presented in this thesis provide opportunities for improved imaging of myocardial function, in all areas of the imaging chain; acquisition, reconstruction and visualization.
|
|
| 10. |
- Svensson, Björn, 1979-
(författare)
-
A Multidimensional Filtering Framework with Applications to Local Structure Analysis and Image Enhancement
- 2008
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Filtering is a fundamental operation in image science in general and in medical image science in particular. The most central applications are image enhancement, registration, segmentation and feature extraction. Even though these applications involve non-linear processing a majority of the methodologies available rely on initial estimates using linear filters. Linear filtering is a well established cornerstone of signal processing, which is reflected by the overwhelming amount of literature on finite impulse response filters and their design.Standard techniques for multidimensional filtering are computationally intense. This leads to either a long computation time or a performance loss caused by approximations made in order to increase the computational efficiency. This dissertation presents a framework for realization of efficient multidimensional filters. A weighted least squares design criterion ensures preservation of the performance and the two techniques called filter networks and sub-filter sequences significantly reduce the computational demand.A filter network is a realization of a set of filters, which are decomposed into a structure of sparse sub-filters each with a low number of coefficients. Sparsity is here a key property to reduce the number of floating point operations required for filtering. Also, the network structure is important for efficiency, since it determines how the sub-filters contribute to several output nodes, allowing reduction or elimination of redundant computations.Filter networks, which is the main contribution of this dissertation, has many potential applications. The primary target of the research presented here has been local structure analysis and image enhancement. A filter network realization for local structure analysis in 3D shows a computational gain, in terms of multiplications required, which can exceed a factor 70 compared to standard convolution. For comparison, this filter network requires approximately the same amount of multiplications per signal sample as a single 2D filter. These results are purely algorithmic and are not in conflict with the use of hardware acceleration techniques such as parallel processing or graphics processing units (GPU). To get a flavor of the computation time required, a prototype implementation which makes use of filter networks carries out image enhancement in 3D, involving the computation of 16 filter responses, at an approximate speed of 1MVoxel/s on a standard PC.
|
|
