The reconstruction and quantitative characterization of organogenesis during the development of living species is an area of great interest in the field of bio-medicine: What are the processes underlying animal embryogenesis where a single cell (a zygote) transforms into a multi-cellular organism comprising a rich diversity of cells in a detailed spatiotemporal organization?. Recent advances in microscopy techniques combined with the application of biological markers have now enabled the capture of in-vivo 3D image sequences of animal models, at the scale of individual cells, as they undergo embryogenesis. Robust and efficient computational techniques are required to meet the challenge of quantitatively processing and analyzing the huge amount of data provided by such in-vivo observations. In this context, this PhD Thesis contributes to the design of mathematical methods that aim automatically to digitally reconstruct the embryogenesis of animal models (mainly zebrafish) from 3D+time data acquired with state-of-the-art microscopy technologies. To accomplish this long-standing goal, the multiscale integration of cell dynamics must be combined with the molecular dynamics that occur within that spatio-temporal context. Towards this goal, on the one hand, cell dynamics can be explored from lineage tree information combined with information concerning cell shape and structure; and on the other, molecular dynamics can be gleaned from the analysis of spatio-temporal patterns of gene expression data. This PhD thesis comprises four parts that cover different aspects of the multiscale digital reconstruction of embryogenesis; it relies upon the use of various distinct methodologies, ranging from image processing to quantitative data analysis for biological interpretations. The first two parts explore the reconstruction of cell dynamics during zebrafish embryogenesis. In the first part, entitled 'In-toto reconstruction of early zebrafish embryogenesis from multiharmonic imaging', a dedicated image processing pipeline (mitosis tracking, membranes segmentation) for digitizing the first three hours of zebrafish development has been implemented. This strategy is employed to produce the lineage trees in a group of six embryos including complete division timings, coordinates, and cell shapes with minute temporal accuracy and micrometer spatial resolution. The analysis of this data shows the existence of a metasynchronous mitotic wave during the first cleavages in zebrafish development. In the second part 'Towards the spatio-temporal atlas of gene expression in zebrafish embryo development', a prototype of the image registration and data quantification workflows needed to generate a gene expression atlas of zebrafish development from fluorescence in-situ hybridization images is shown. The viability of the proposed approach is illustrated with a proof-ofconcept example that aligns and quantifies five gene expression patterns on a zebrafish template at the beginning of gastrulation. In the next two parts, new generic tools for reconstructing cell dynamics are presented. The third part, entitled 'Spatio-temporal mathematical morphology processing for in-vivo embryogenesis imaging', explores the extension of morphological image processing to 4D datasets. The definition of the basic spatio-temporal structuring elements (hypersegment, hypersphere, hypertube, hyperhorn) allows to present several applications (filtering, tracking, segmentation) dedicated to the analysis of 4D datasets of embryogenesis acquired with different microscopy techniques. Working directly in 4D spaces provides more spatio-temporal coherence than classical methods. Finally, in the fourth part -entitled 'Measuring the cell lineage tree'- a quantitative method that characterizes cell lineage tree data is presented. Because of the lack of mathematical formalisms in this field (usually only qualitative descriptions are provided), the Morphogenetic Entropy measure is defined. It allows one to quantify the informational content of a cell inside a lineage tree, concerning pluripotency and behavior during cell division (proliferation, diversification, stem cell mode). Morphogenetic Entropy measurements can be broadly applied to research on embryogenesis and stem cells. Overall, this PhD Thesis presents mathematical methods that when applied to in-vivo images from animal models (e.g. zebrafish, sea urchin) allow the quantitative description of embryogenesis processes at a cellular scale. The proposed methods analyze structural information, cell lineage trees and gene expression data. From a broader perspective, the use of these methods should contribute to fill the gap between genetics and epigenetics by providing a deeper understanding of embryogenesis.

Los filtros adaptativos constituyen una herramienta importante en el ámbito del procesa­do digital de señales estadísticas. Frente a los filtros convencionales, los filtros adaptativos permite operar en entornos donde los estadísticos de las señales que participan en la ope­ración de filtrado no son conocidos a priori. Por consiguiente, son capaces de adaptarse a las condiciones cambiantes del entorno, lo que ha permitido no solo mejorar las presta­ciones ofrecidas por los filtros de coeficientes fijos sino que también el desarrollo de apli­

caciones nuevas.

Para poder adaptarse a escenarios desconocidos, los filtros adaptativos están dotados con un algoritmo que les permite modificar el valor de sus coeficientes en función de la infor­mación presente en el entorno. La manera en que el algoritmo estima el valor de los coefi­cientes determinará las características del filtro tales como su velocidad de convergencia, complejidad computacional, robustez numérica, etc.

El objetivo de este proyecto ha sido implementar en hardware un filtro adaptativo en for­ma de array sistólico destinado a la identificación de sistemas. Disponer el algoritmo en forma sistólica consiste en descomponer las operaciones que constituyen el algoritmo de adaptación en pequeñas unidades de cálculo que operaran de manera paralela para actua­lizar el vector de pesos del filtro. Para poder cumplir dicho objetivo se han estudiado las diferentes alternativas que en la actualidad son capaces de solucionar el problema adapta­tivo y en las que, además, es posible descomponer sus operaciones en forma de array sistólico. Entre todos los algoritmos estudiados se escogió el algoritmo QRD-RLS inverso, ya que supone un buen compromiso entre velocidad de convergencia, complejidad compu­tacional y sencíllez de la estructura.

Se han realizado dos implementaciones, ambas descritas utilizando el lenguaje VHDL. La primera consistió en una estructura no sintetizable en la que se empleó aritmética en co­ma flotante. Dicha estructura sirvió como paso previo a la transformación del array sistóli­co en una estructura sintetizable, ya que facilitó aspectos como la cuantificación y la depu­ración. La segunda implementación, realizada utilizando aritmética en punto fijo, supone la materialización del algoritmo en una estructura sintetizable. En ambas implementacio­nes se pueden modificar los parámetros del filtro a partir de un conjunto de constantes y sin necesidad de tener que recodificar la estructura.

A partir de las implementaciones realizadas, se han cubierto dos objetivos. Por un lado se ha creado una plataforma de trabajo que facilitará ajustar el filtro de tal forma que sea posible adecuarlo al escenario real elegido. Por otro, lado se ha comprobado que la imple­mentación en punto fijo del algoritmo es viable y que el algoritmo preserva sus caracterís­ticas funcionales y su estabilidad numérica.

Palabras clave

Filtros adaptativos, identificación de sistemas, algoritmos adaptativos, array sistólico, QRD-RLS inverso, QRD-RLS, RLS, LMS, VHDL.