|
|
Publications sur Estimation bayesienne
Résultat de la recherche dans la liste des publications :
Article |
1 - Invariant Bayesian estimation on manifolds.
I. H. Jermyn. Annals of Statistics, 33(2): pages 583--605, avril 2005.
Mots-clés : Estimation bayesienne, MAP, MMSE, Invariant, Metrique, Jeffrey's.
@ARTICLE{jermyn_annstat05,
|
| author |
= |
{Jermyn, I. H.}, |
| title |
= |
{Invariant Bayesian estimation on manifolds}, |
| year |
= |
{2005}, |
| month |
= |
{avril}, |
| journal |
= |
{Annals of Statistics}, |
| volume |
= |
{33}, |
| number |
= |
{2}, |
| pages |
= |
{583--605}, |
| url |
= |
{http://dx.doi.org/10.1214/009053604000001273}, |
| pdf |
= |
{ftp://ftp-sop.inria.fr/ariana/Articles/jermyn_annstat05.pdf}, |
| keyword |
= |
{Estimation bayesienne, MAP, MMSE, Invariant, Metrique, Jeffrey's} |
| } |
Abstract :
A frequent and well-founded criticism of the maximum em a posteriori (MAP) and minimum mean squared error (MMSE) estimates of a continuous parameter param taking values in a differentiable manifold paramspace is that they are not invariant to arbitrary `reparametrizations' of paramspace. This paper clarifies the issues surrounding this problem, by pointing out the difference between coordinate invariance, which is a em sine qua non for a mathematically well-defined problem, and diffeomorphism invariance, which is a substantial issue, and then provides a solution. We first show that the presence of a metric structure on paramspace can be used to define coordinate-invariant MAP and MMSE estimates, and we argue that this is the natural way to proceed. We then discuss the choice of a metric structure on paramspace. By imposing an invariance criterion natural within a Bayesian framework, we show that this choice is essentially unique. It does not necessarily correspond to a choice of coordinates. In cases of complete prior ignorance, when Jeffreys' prior is used, the invariant MAP estimate reduces to the maximum likelihood estimate. The invariant MAP estimate coincides with the minimum message length (MML) estimate, but no discretization or approximation is used in its derivation. |
|
 haut de la page
Article de conférence |
1 - Parametric blind deconvolution for confocal laser scanning microscopy.
P. Pankajakshan et B. Zhang et L. Blanc-Féraud et Z. Kam et J.C. Olivo-Marin et J. Zerubia. Dans Proc. 29th International Conference of IEEE EMBS (EMBC-07), pages 6531-6534, août 2007.
Mots-clés : Microscopie confocale, Blind Deconvolution, Poisson noise, Variation totale, Algorithme EM, Estimation bayesienne.
Copyright : ©2007 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or redistribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from the IEEE.
@INPROCEEDINGS{Pankajakshan07a,
|
| author |
= |
{Pankajakshan, P. and Zhang, B. and Blanc-Féraud, L. and Kam, Z. and Olivo-Marin, J.C. and Zerubia, J.}, |
| title |
= |
{Parametric blind deconvolution for confocal laser scanning microscopy}, |
| year |
= |
{2007}, |
| month |
= |
{août}, |
| booktitle |
= |
{Proc. 29th International Conference of IEEE EMBS (EMBC-07)}, |
| pages |
= |
{6531-6534}, |
| pdf |
= |
{http://ieeexplore.ieee.org/iel5/4352184/4352185/04353856.pdf?tp=&isnumber=&arnumber=4353856}, |
| keyword |
= |
{Microscopie confocale, Blind Deconvolution, Poisson noise, Variation totale, Algorithme EM, Estimation bayesienne} |
| } |
Abstract :
In this paper, we propose a method for the
iterative restoration of fluorescence Confocal Laser Scanning
Microscopic (CLSM) images and parametric estimation of the
acquisition system’s Point Spread Function (PSF). The CLSM is
an optical fluorescence microscope that scans a specimen in 3D
and uses a pinhole to reject most of the out-of-focus light. However,
the quality of the images suffers from two basic physical
limitations. The diffraction-limited nature of the optical system,
and the reduced amount of light detected by the photomultiplier
cause blur and photon counting noise respectively. These images
can hence benefit from post-processing restoration methods
based on deconvolution. An efficient method for parametric
blind image deconvolution involves the simultaneous estimation
of the specimen 3D distribution of fluorescent sources and
the microscope PSF. By using a model for the microscope
image acquisition physical process, we reduce the number of
free parameters describing the PSF and introduce constraints.
The parameters of the PSF may vary during the course of
experimentation, and so they have to be estimated directly from
the observed data. A priori model of the specimen is further
applied to stabilize the alternate minimization algorithm and to
converge to the solutions. |
|
haut de la page
4 Rapports de recherche et Rapports techniques |
1 - Space non-invariant point-spread function and its estimation in fluorescence microscopy.
P. Pankajakshan et L. Blanc-Féraud et Z. Kam et J. Zerubia. Research Report 7157, INRIA, décembre 2009.
Mots-clés : Confocal Laser Scanning Microscopy, point spread function, Estimation bayesienne, Estimation MAP, Deconvolution, fluorescence microscopy.
@TECHREPORT{ppankajakshan09c,
|
| author |
= |
{Pankajakshan, P. and Blanc-Féraud, L. and Kam, Z. and Zerubia, J.}, |
| title |
= |
{Space non-invariant point-spread function and its estimation in fluorescence microscopy}, |
| year |
= |
{2009}, |
| month |
= |
{décembre}, |
| institution |
= |
{INRIA}, |
| type |
= |
{Research Report}, |
| number |
= |
{7157}, |
| url |
= |
{http://hal.archives-ouvertes.fr/inria-00438719/en/}, |
| keyword |
= |
{Confocal Laser Scanning Microscopy, point spread function, Estimation bayesienne, Estimation MAP, Deconvolution, fluorescence microscopy} |
| } |
Résumé :
Dans ce rapport de recherche, nous rappelons brièvement comment la nature limitée de diffraction de l'objectif d'un microscope optique, et le bruit
intrinsèque peuvent affecter la résolution d'une image observée. Un algorithme de déconvolution aveugle a été proposé en vue de restaurer les fréquences manquants au delà de la limite de diffraction. Cependant, sous d'autres conditions, l'approximation du systéme imageur l'imagerie sans aberration n'est plus valide et donc les aberrations de la phase du front d'onde émergeant d'un médium ne sont plus ignorées. Dans la deuxième partie de
ce rapport de recherche, nous montrons que la distribution d'intensité originelle et la localisation d'un objet peuvent être retrouvées uniquement en obtenant de la phase du front d'onde
réfracté, à partir d'images d'intensité observées. Nous démontrons cela par obtention de la fonction de ou a partir d'une microsphère imagée. Le bruit et l'infl uence de la taille de la
microsphère peuvent être diminués et parfois complètement supprimes des images observées en utilisant un estimateur maximum a posteriori. Néanmoins, a cause de l'incohérence du système d'acquisition, une récupération de phase a partir d'intensités observées n'est possible que si la restauration de la phase est contrainte. Nous avons utilisé l'optique géométrique
pour modéliser la phase du front d'onde réfracté, et nous avons teste l'algorithme sur des images simulées. |
Abstract :
In this research report, we recall briefly how the diffraction-limited nature of an optical microscope's objective, and the intrinsic noise can affect the observed images' resolution. A blind deconvolution algorithm can restore the lost frequencies beyond the diffraction limit. However, under other imaging conditions, the approximation of aberration-free imaging, is not applicable, and the phase aberrations of the emerging wavefront from a specimen immersion medium cannot be ignored any more. We show that an object's location and its original intensity distribution can be recovered by retrieving the refracted wavefront's phase from the observed intensity images. We demonstrate this by retrieving the point-spread function from an imaged microsphere. The noise and the influence of the microsphere size can be mitigated and sometimes completely removed from the observed images by using a maximum a posteriori estimate. However, due to the incoherent nature of the acquisition system, phase retrieval from the observed intensities will be possible only if the phase is constrained. We have used geometrical optics to model the phase of the refracted wavefront, and tested the algorithm on some simulated images. |
|
2 - A Binary Tree-Structured MRF Model for Multispectral Satellite Image Segmentation.
G. Scarpa et G. Poggi et J. Zerubia. Rapport de Recherche 5062, INRIA, France, décembre 2003.
Mots-clés : Estimation bayesienne, Classification, Champs de Markov, Modeles hierarchiques.
@TECHREPORT{Scarpa03,
|
| author |
= |
{Scarpa, G. and Poggi, G. and Zerubia, J.}, |
| title |
= |
{A Binary Tree-Structured MRF Model for Multispectral Satellite Image Segmentation}, |
| year |
= |
{2003}, |
| month |
= |
{décembre}, |
| institution |
= |
{INRIA}, |
| type |
= |
{Research Report}, |
| number |
= |
{5062}, |
| address |
= |
{France}, |
| url |
= |
{https://hal.inria.fr/inria-00071522}, |
| pdf |
= |
{https://hal.inria.fr/file/index/docid/71522/filename/RR-5062.pdf}, |
| ps |
= |
{https://hal.inria.fr/docs/00/07/15/22/PS/RR-5062.ps}, |
| keyword |
= |
{Estimation bayesienne, Classification, Champs de Markov, Modeles hierarchiques} |
| } |
Résumé :
Dans ce rapport, nous proposons un modèle markovien a priori structuré à arbre binaire (le TS-MRF) pour la segmentation d'images satellitaires multispectrales. Ce modèle permet de représenter un champ bidimensionnel par une séquence de champs de Markov binaires, chacun correspondant à un noeud de l'arbre. Pour avoir une bonne classification, on peut adapter le modèle TS-MRF à la structure intrinsèque des données, en définissant un MRF, à plusieurs paramètres, très flexible. Bien que l'on définisse le modèle global sur tout l'arbre, l'optimisation et l'estimation peuvent être poursuivis en considérant un noeud à la fois, à partir de la racine jusqu'aux feuilles, avec une réduction significative de la complexité. En effet, on a montré expérimentalement que l'algorithme global est beaucoup plus rapide qu'un algorithme conventionnel fondé sur le modèle markovien d'Ising, en particulier quand le nombre des bandes spectrales est très grand. Grâce à la procédure d'optimisation séquentielle, ce modèle permet aussi de déterminer le nombre des classes présentes dans l'image satellitaire, dans le cadre d'une classification non supervisée, à travers une condition d'arrêt définie localement pour chaque noeud. Nous avons effectué des expériences sur une image SPOT de la baie de Lannion, pour laquelle nous disposons d'une vérité terrain, et nous avons trouvé que le modèle proposé fournit de meilleurs résultats que certains autres modèles de Markov et que d'autres méthodes variationnelles. |
Abstract :
In this work we detail a tree-structured MRF (TS-MRF) prior model useful for segmentation of multispectral satellite images. This model allows a hierarchical representation of a 2-D field by the use of a sequence of binary MRFs, each corresponding to a node in the tree. In order to get good performances, one can fit the intrinsic structure of the data to the TS-MRF model, thereby defining a multi-parameter, flexible, MRF. Although a global MRF model is defined on the whole tree, optimization as well estimation can be carried out by working on a single node at a time, from the root down to the leaves, with a significant reduction in complexity. Indeed the overall algorithm is proved experimentally to be much faster than a comparable algorithm based on a conventional Ising MRF model, especially when the number of bands becomes very large. Thanks to the sequential optimization procedure, this model also addresses the cluster validation problem of unsupervised segmentation, through the use of a stopping condition local to each node. Experiments on a SPOT image of the Lannion Bay, a ground-truth of which is available, prove the superior performance of the algorithm w.r.t. some other MRF based algorithms for supervised segmentation, as well as w.r.t. some variational methods. |
|
3 - On Bayesian Estimation in Manifolds.
I. H. Jermyn. Rapport de Recherche 4607, Inria, France, novembre 2002.
Mots-clés : Évenement rare, Estimation bayesienne, Invariant.
@TECHREPORT{4607,
|
| author |
= |
{Jermyn, I. H.}, |
| title |
= |
{On Bayesian Estimation in Manifolds}, |
| year |
= |
{2002}, |
| month |
= |
{novembre}, |
| institution |
= |
{Inria}, |
| type |
= |
{Research Report}, |
| number |
= |
{4607}, |
| address |
= |
{France}, |
| url |
= |
{https://hal.inria.fr/inria-00071978}, |
| pdf |
= |
{https://hal.inria.fr/file/index/docid/71978/filename/RR-4607.pdf}, |
| ps |
= |
{https://hal.inria.fr/docs/00/07/19/78/PS/RR-4607.ps}, |
| keyword |
= |
{Évenement rare, Estimation bayesienne, Invariant} |
| } |
Résumé :
Il est fréquemment dit que les estimées au sens du maximum a posteriori (MAP) et du minimum de l'erreur quadratique moyenne (MMSE) d'un paramètre continu ne sont pas invariantes relativement aux «reparamètrisations» de l'espace des paramètres . Ce rapport clarifie les questions autour de ce problème, en soulignant la différence entre l'invariance aux changements de coordonnées, qui est une condition sine qua non pour un problème mathématiq- uement bien défini, et l'invariance aux difféomorphismes, qui est une question significative, et fournit une solution. On montre d'abord que la présence d'une structure métrique sur peut être utilisée pour définir les estimées aux sens du MAP et du MMSE qui sont invariantes aux changements de coordonnées, et on explique pourquoi cela est la fa on naturelle et nécessaire pour le faire. Le problème de l'estimation et les quantités géométriques qui y sont associées sont tous définis d'une fa on clairement invariante aux changements de coordonnées. On montre que la même estimée au sens du MAP est obtenue en utilisant soit la `maximisation d'une densité' soit une fonction de perte delta, définie de fa on invariante. Puis, on discute le choix d'une métrique pour . En imposant un critère d'invariance qui est naturel dans le cadre bayesien, on montre que ce choix est unique. Il ne correspond pas nécessairement à un choix de coordonnées. L'estimée au sens du MAP qui en résulte coincide avec l'estimée fondée sur la longueur minimum de message (MML), mais la demonstration n'utilise pas de discrétisation ou d'approximation. |
Abstract :
It is frequently stated that the maximum a posteriori (MAP) and minimum mean squared error (MMSE) estimates of a continuous parameter are not invariant to arbitrary «reparametrizations» of the parameter space . This report clarifies the issues surrounding this problem, by pointing out the difference between coordinate invariance, which is a sine qua non for a mathematically well-defined problem, and diffeomorphism invariance, which is a substantial issue, and provides a solution. We first show that the presence of a metric structure on can be used to define coordinate-invari- ant MAP and MMSE estimates, and we argue that this is the natural and necessary way to proceed. The estimation problem and related geometrical quantities are all defined in a manifestly coordinate-invariant way. We show that the same MAP estimate results from `density maximization' or from using an invariantly-defined delta function loss. We then discuss the choice of a metric structure on . By imposing an invariance criterion natural within a Bayesian framework, we show that this choice is essentially unique. It does not necessarily correspond to a choice of coordinates. The resulting MAP estimate coincides with the minimum message length (MML) estimate, but no discretization or approximation is used in its derivation. |
|
4 - Satellite image deconvolution using complex wavelet packets.
A. Jalobeanu et L. Blanc-Féraud et J. Zerubia. Rapport de Recherche 3955, Inria, juin 2000.
Mots-clés : Deconvolution, Estimation bayesienne, Paquet d'ondelettes.
@TECHREPORT{jalo00,
|
| author |
= |
{Jalobeanu, A. and Blanc-Féraud, L. and Zerubia, J.}, |
| title |
= |
{Satellite image deconvolution using complex wavelet packets}, |
| year |
= |
{2000}, |
| month |
= |
{juin}, |
| institution |
= |
{Inria}, |
| type |
= |
{Research Report}, |
| number |
= |
{3955}, |
| url |
= |
{https://hal.inria.fr/inria-00072694}, |
| pdf |
= |
{https://hal.inria.fr/file/index/docid/72694/filename/RR-3955.pdf}, |
| ps |
= |
{https://hal.inria.fr/docs/00/07/26/94/PS/RR-3955.ps}, |
| keyword |
= |
{Deconvolution, Estimation bayesienne, Paquet d'ondelettes} |
| } |
Résumé :
La déconvolution des images satellitaires floues et bruitées est un problème inverse mal posé. L'inversion directe entraîne une amplification inacceptable du bruit. Généralement, soit le problème est régularisé lors de l'inversion, soit le bruit est filtré après déconvolution et décomposition dans le domaine de la transformée en ondelettes. Nous avons developpé dans ce rapport la deuxième solution, en seuillant les coefficients d'une nouvelle transformée en paquets d'ondelettes complexes, les fonctions de seuillage étant estimées de manière automatique. L'utilisation de paquets d'ondelettes complexes rend cette méthode invariante par translation, et tient compte des directions, tout en restant d'une complexité O(N). Les résultats obtenus présentent à la fois des textures nettes et un très bon rapport signal/bruit dans les zones homogènes. Par rapport aux algorithmes concurrents, la méthode que nous proposons est plus rapide, invariante par rotation, et tient compte de la directionnalité des détails et des textures de l'image pour mieux les restaurer. Les images déconvoluées de cette manière peuvent être utilisées telles quelles (la restauration peut être intégrée directement dans la chaîne d'acquisition). Mais elles peuvent également constituer le point de départ d'une méthode de régularisation adaptative, permettant d'obtenir des contours plus francs. |
Abstract :
The deconvolution of blurred and noisy satellite images is an ill-posed inverse problem. The direct inversion leads to unacceptable noise amplificatio- n. Usually, either the problem is regularized during the inversion process, or the noise is filtered after deconvolution and decomposition in the wavelet transform domain. Herein, we have developed the second solution, by thresholding the coefficients of a new complex wavelet packet transform; the thresholding functions are automatically estimated. The use of complex wavelet packets enables translation invariance, and takes into account the directions, while remaining of complexity O(N). The obtained results exhibit both correctly restored textures and a high SNR in homogeneous areas. Compared to concurrent algorithms, the proposed method is faster, rotation invariant and takes into account the directions of the details and textures of the image to restore them better. The images deconvolved this way can be used as they are (the restoration step proposed here can be directly inserted in the acquisition chain). But they also can provide a starting point of an adaptive regularization method, enabling one to obtain sharper edges. |
|
haut de la page
Ces pages sont générées par 
|
Rapport d'activité 2011
bibtex
