Submitted by Renata Lopes (renatasil82@gmail.com) on 2017-02-23T19:26:50Z No. of bitstreams: 1 francianeconceicaopeters.pdf: 4730497 bytes, checksum: 201c60342a8bf9edc9b308fa50fafa54 (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2017-02-24T12:08:43Z (GMT) No. of bitstreams: 1 francianeconceicaopeters.pdf: 4730497 bytes, checksum: 201c60342a8bf9edc9b308fa50fafa54 (MD5) / Made available in DSpace on 2017-02-24T12:08:43Z (GMT). No. of bitstreams: 1 francianeconceicaopeters.pdf: 4730497 bytes, checksum: 201c60342a8bf9edc9b308fa50fafa54 (MD5) Previous issue date: 2010-01-27 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A construção de imagens associadas à distribuição de condutividades no interior de um meio condutor a partir de injeção de corrente elétrica e medidas de potencial no contorno externo do corpo é uma técnica conhecida como tomografia por impedância elétrica. É um problema inverso que tem sido estudado visando aplicações biomédicas, monitoramento de processos industriais e investigação geofísica. Em alguns casos, é possível levar em consideração informações conhecidas sobre o domínio do corpo no processo de construção da imagem, recaindo no problema da detecção de inclusões que é o problema efetivamente tratado neste trabalho. Este problema pode ser resolvido por meio da minimização de uma função da diferença entre potenciais medidos no contorno e calculados para uma dada distribuição de condutividades. O presente trabalho desenvolve uma estratégia para a solução deste problema baseada na parametrização da geometria do contorno das inclusões cujas formas e dimensões se pretende determinar. O problema de minimização é resolvido por meio do Método de Levenberg-Marquardt e o problema direto via Método dos Elementos de Contorno. Para avaliar o desempenho da estratégia proposta são apresentados resultados numéricos envolvendo contornos definidos por splines, problemas com a presença de ruído nas medidas, avaliação de protocolos de injeção de corrente e medição de potencial elétrico e ainda uma aplicação voltada ao monitoramento cardíaco. / The images reconstruction of the conductivity distribution inside a conductive body based on electrical current injection and potential measurements on the outer boundary of this body is a technique known as electrical impedance tomography. This is an inverse problem that has been studied in biomedical applications, industrial process monitoring and geophysics investigation. In some cases, it is possible to take into account in the reconstruction process, informations about the body, leading to the problem of identifying inclusions, that is the problem actually treated in this work. This inverse problem can be solved by the minimization of a function, defined as the difference between the measured potentials and the computed ones for a given conductivity distribution. The present work describes a strategy to solve this problem based on the parametrization of the inclusions boundary, whose shape and size is intended to be determined. The minimization problem is solved via Levenberg-Marquardt Method and the forward one is solved via Boundary Elements Method. In order to evaluate the performance of the proposed strategy, numerical experiments with inclusions of boundaries defined by splines, problems with noisy data, current injection and potential measurement protocols and an application of the strategy to the cardiac function monitoring are presented.

CAPES / Em ensaios não destrutivos por ultrassom, o sinal obtido a partir de um sistema de aquisição de dados real podem estar contaminados por ruído e os ecos podem ter atrasos de tempo subamostrais. Em alguns casos, esses aspectos podem comprometer a informação obtida de um sinal por um sistema de aquisição. Para lidar com essas situações, podem ser utilizadas técnicas de estimativa de atraso temporal (Time Delay Estimation ou TDE) e também técnicas de reconstrução de sinais, para realizar aproximações e obter mais informações sobre o conjunto de dados. As técnicas de TDE podem ser utilizadas com diversas finalidades na defectoscopia, como por exemplo, para a localização precisa de defeitos em peças, no monitoramento da taxa de corrosão em peças, na medição da espessura de um determinado material e etc. Já os métodos de reconstrução de dados possuem uma vasta gama de aplicação, como nos NDT, no imageamento médico, em telecomunicações e etc. Em geral, a maioria das técnicas de estimativa de atraso temporal requerem um modelo de sinal com precisão elevada, caso contrário, a localização dessa estimativa pode ter sua qualidade reduzida. Neste trabalho, é proposto um esquema alternado que estima de forma conjunta, uma referência de eco e atrasos de tempo para vários ecos a partir de medições ruidosas. Além disso, reinterpretando as técnicas utilizadas a partir de uma perspectiva probabilística, estendem-se suas funcionalidades através de uma aplicação conjunta de um estimador de máxima verossimilhança (Maximum Likelihood Estimation ou MLE) e um estimador máximo a posteriori (MAP). Finalmente, através de simulações, resultados são apresentados para demonstrar a superioridade do método proposto em relação aos métodos convencionais. / Abstract (parágrafo único): In non-destructive testing (NDT) with ultrasound, the signal obtained from a real data acquisition system may be contaminated by noise and the echoes may have sub-sample time delays. In some cases, these aspects may compromise the information obtained from a signal by an acquisition system. To deal with these situations, Time Delay Estimation (TDE) techniques and signal reconstruction techniques can be used to perform approximations and also to obtain more information about the data set. TDE techniques can be used for a number of purposes in the defectoscopy, for example, for accurate location of defects in parts, monitoring the corrosion rate in pieces, measuring the thickness of a given material, and so on. Data reconstruction methods have a wide range of applications, such as NDT, medical imaging, telecommunications and so on. In general, most time delay estimation techniques require a high precision signal model, otherwise the location of this estimate may have reduced quality. In this work, an alternative scheme is proposed that jointly estimates an echo model and time delays for several echoes from noisy measurements. In addition, by reinterpreting the utilized techniques from a probabilistic perspective, its functionalities are extended through a joint application of a maximum likelihood estimator (MLE) and a maximum a posteriori (MAP) estimator. Finally, through simulations, results are presented to demonstrate the superiority of the proposed method over conventional methods.

CNPq; FUNTEF-PR / Sensores wire-mesh (WMS) são capazes de gerar imagens da seção transversal de escoamentos multifásicos e tem sido amplamente utilizados para investigar fenômenos de escoamentos em plantas piloto. Tais dispositivos são capazes de medir parâmetros de escoamento tais como distribuição da fração de fase (por exemplo fração de gás ou líquido) e visualizar escoamentos multifásicos com alta resolução temporal e espacial. Sendo portanto, uma ferramenta importante para investigações de escoamentos mais detalhadas. No entanto, seu princípio de medição é baseado em eletrodos intrusivos posicionados dentro do tubo onde o escoamento flui. A resolução da imagem gerada pelo sensor é dada pelo número de cruzamentos entre os fios transmissores e receptores. Em muitos processos, no entanto, efeitos de intrusividade de tal sensor pode ser uma limitação no seu uso. Por isso, um número reduzido de fios poderia permitir uma expansão do campo de aplicações do sensor wire-mesh. Por essa razão, o presente trabalho sugere um método de reconstrução de imagem para aumentar a resolução dos dados de um sensor wire-mesh com um número de eletrodos menor que o ótimo. Desta forma, os efeitos de intrusividade no processo investigado poderiam ser reduzidos. O método de reconstrução é baseado em uma abordagem estatística de regularização e é conhecido como Maximum a Posteriori (MAP). Dados de escoamento de um WMS 16x16 são usados para determinar um modelo gaussiano multivariável do escoamento, o qual são empregados como regularização na reconstrução. Uma matriz de sensitividade é estimada pelo método de elementos finitos (FEM) para incorporar o algoritmo MAP. Dados experimentais são usados para validar o método proposto, sendo comparado com interpolação do tipo spline. Resultados experimentais mostram que a reconstrução por MAP possui um desempenho melhor do que interpolação do tipo spline, alcançando desvios de fração de vazio dentro de uma faixa de ± 10% na grande maioria dos pontos de operação. A validação foi executada em um loop de escoamento horizontal água/gás em regime intermitente (golfada). / Wire-mesh sensors (WMS) are able to generate cross-sectional images of multiphase flow and have been widely used to investigate flow phenomena in pilot plant studies. Such devices are able to measure flow parameters such as phase fraction (e.g. gas/liquid fraction) distribution and visualize multiphase flows with high temporal and spatial resolution. Hence, being important tool for detailed flow investigation. However, its sensing principle is based on intrusive electrodes placed inside the pipe where a multiphase flow streams. The image resolution generated by the sensor is given by the number of crossing points formed by the transmitter and receptor wires. In many processes, however, the intrusive effect of such sensor might be a limitation on its use. Therefore, a reduced number of wires could possibly increase the application field of wire-mesh sensors. For this reason, the present work presents an image reconstruction method to increase resolution of WMS data with less than optimal number of electrode wires. In this way, a reduction of intrusive effects on the process under investigation may be achieved. The reconstruction method is based on statistical view of regularization and is known as Maximum a Posterior (MAP). 16x16 WMS flow data are used to determine a Multivariate Gaussian flow model, which in turn is used as regularization in the reconstruction. A sensitive matrix is estimated by finite element method (FEM) to incorporate MAP algorithm. Experimental data are used to validate the proposed method, which is compared with spline interpolation. Experimental results show that the MAP reconstruction performs better than interpolation and achieves deviation in gas void fraction estimation in the range of ±10% in the vast majority of operating points. The tests were performed in a horizontal water-gas flow loop operating at intermittent (slug) flow regime.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Neste trabalho, três técnicas para resolver numericamente problemas inversos de transporte de partículas neutras a uma velocidade para aplicações em engenharia nuclear são desenvolvidas. É fato conhecido que problemas diretos estacionários e monoenergéticos de transporte são caracterizados por estimar o fluxo de partículas como uma função-distribuição das variáveis independentes de espaço e de direção de movimento, quando os parâmetros materiais (seções de choque macroscópicas), a geometria, e o fluxo incidente nos contornos do domínio (condições de contorno), bem como a distribuição de fonte interior são conhecidos. Por outro lado, problemas inversos, neste trabalho, buscam estimativas para o fluxo incidente no contorno, ou a fonte interior, ou frações vazio em barras homogêneas. O modelo matemático usado tanto para os problemas diretos como para os problemas inversos é a equação de transporte independente do tempo, a uma velocidade, em geometria unidimensional e com o espalhamento linearmente anisotrópico na formulação de ordenadas discretas (SN). Nos problemas inversos de valor de contorno, dado o fluxo emergente em um extremo da barra, medido por um detector de nêutrons, por exemplo, buscamos uma estimativa precisa para o fluxo incidente no extremo oposto. Por outro lado, nos problemas inversos SN de fonte interior, buscamos uma estimativa precisa para a fonte armazenada no interior do domínio para fins de blindagem, sendo dado o fluxo emergente no contorno da barra. Além disso, nos problemas inversos SN de fração de vazio, dado o fluxo emergente em uma fronteira da barra devido ao fluxo incidente prescrito no extremo oposto, procuramos por uma estimativa precisa da fração de vazio no interior da barra, no contexto de ensaios não-destrutivos para aplicações na indústria. O código computacional desenvolvido neste trabalho apresenta o método espectronodal de malha grossa spectral Greens function (SGF) para os problemas diretos SN em geometria unidimensional para gerar soluções numéricas precisas para os três problemas inversos SN descritos acima. Para os problemas inversos SN de valor de contorno e de fonte interior, usamos a propriedade da proporcionalidade da fuga de partículas; ademais, para os problemas inversos SN de fração de vazio, oferecemos a técnica a qual nos referimos como o método físico da bissecção. Apresentamos resultados numéricos para ilustrar a precisão das três técnicas, conforme descrito nesta tese. / In this work, three techniques for numerically solving one-speed neutral particle inverse transport problems for nuclear engineering applications are developed. It is well known that direct steady-state monoenergetic transport problems are characterized by estimating the flux of particles as a distribution function of space and direction-of-motion independent variables, when the material parameters (cross sections), the geometry, and the incoming flux at the boundaries of the domain (boundary conditions), as well as the interior source distribution are known. Conversely, inverse problems, in this work, seek for estimates to the incident boundary flux, or interior source, or void fractions in homogeneous slabs. The mathematical model used for direct and inverse problems is the time-independent one-speed slab-geometry transport equation with linearly anisotropic scattering in the discrete ordinates (SN) formulation. In the boundary-value inverse problems, given the existing flux at one boundary of the slab, as measured by a neutron detector, for example, we seek for accurate estimate for the incident flux at the opposite boundary. On the other hand, in the interior source inverse SN problems, we seek for accurate estimate for the interior source stored within the slab for shielding purpose, given the exiting flux at the boundary of the slab. Furthermore, as with the void fraction inverse SN problems, given the exiting flux at one boundary of the slab due to prescribed incident flux at the opposite boundary, we seek for accurate estimate of the void fraction within the slab in the context of non-destructive testing applications in industry. The computer code developed in this work presents the coarse-mesh spectral Greens function (SGF) nodal method for direct SN problems in slab geometry to generate accurate numerical solutions to the three inverse SN problems described above. For the boundary-value and interior source inverse SN problems, we use the proportionality property of the leakage of particles; moreover, for the void fraction inverse SN problems, we offer the technique that we refer to as the physical bisection method. We present numerical results to illustrate the accuracy of the three techniques, as described in this dissertation.

Nesta tese mostraremos uma aplicação do método SGF, cf., spectral Greens function, para gerar fluxos angulares nas interfaces dos nodos na formulação de ordenadas discretas (SN) para equação de transporte de partículas neutras em geometria cartesiana unidimensional e a uma velocidade, usando quadraturas de alta ordem. Utilizando o método de malha grossa SGF, primeiramente determinamos as constantes arbitrárias da solução geral analítica das equações SN em cada nodo espacial. Então usamos a fórmula de quadratura angular para estimar uma expressão para o fluxo escalar de nêutrons e substituímos no termo de fonte de espalhamento isotrópico na equação de transporte. Resolvemos analiticamente a equação unidimensional de transporte de nêutrons com a fonte de espalhamento aproximada desta maneira e geramos os valores para o fluxo angular no interior de cada nodo espacial. Como o método SGF gera soluções numéricas de malha grossa completamente livres de erro de truncamento espacial, esperamos que o esquema de reconstrução analítica proposto tenha alta precisão para os fluxos angulares, considerando as condições de continuidade nas interfaces dos nodos do domínio espacial. Esta técnica caracteriza um problema inverso, pois a partir da solução de malha grossa do método SGF, podemos reconstruir o fluxo angular de nêutrons em qualquer ponto do domínio. / We offer in this work an application of the SGF method, cf., spectral Greens function, to generate the angular fluxes at the region interfaces of multilayer slabs, using high-order angular quadrature sets in the one-speed discrete ordinates (SN) formulation of the neutral particle transport equation. Using the SGF coarse-mesh numerical solution, we first determine the arbitrary constants of the analytical general solution of the SN equations within each spatial node. Then, we use the angular quadrature formula to estimate the expression of the scalar flux distribution, that we substitute into the isotropic scattering source term of the transport equation. We solve analytically the slab-geometry transport equation, with the scattering source so approximated, in order to generate the angular flux profile within each spatial node. As the SGF method generates coarse-mesh numerical solution, which is completely free from spatial truncation errors, we expect that the offered approximate analytical reconstruction scheme be accurate enough for the localized angular flux distribution, considering the node interface continuity conditions within the domain. This technique is thought of as an inverse problem since from the SGF coarse-mesh nodal solution we are able to reconstruct the angular flux profile at any point inside the domain.

Apresenta-se um estudo em Teoria de Domínios das equações integrais da forma geral f (x) = h(x)+g Z b(x) a(x) g(x, y, f (y))dy com h, a e b definidas para x ∈ [a0,b0], a0 ≤a(x)≤b(x)≤b0 e g definida para x, y ∈ [a0,b0], cujo lado direito define uma contração sobre o espaço métrico de funções reais contínuas limitadas. O ponto de partida desse trabalho é a reescrita da Análise Intervalar para Teoria de Domínios do problema de valor incial em equações diferenciais ordinárias que possuem solução como ponto fixo do operador de Picard. Com o conjunto dos números reais interpretados pelo Domínio Intervalar, as funções reais são estendidas para operarem no domínio de funçoes intervalares de variável real. Em particular, faz-se a extensão canônica do campo vetorial em relação à segunda variável. Nesse contexto, pela primeira vez tem-se o estudo das equações integrais de Fredholm e Volterra sobre o domínio de funções intervalares de variável real definida pelo operador integral intervalar com a participação da extensão canônica de g em relação à terceira variável. Adicionando ao domínio de funções intervalares sua função medição, efetua-se a análise da convergência do operador intervalar de Fredholm e Volterra em Teoria de Domínios com o cálculo da sua derivada informática em relação à medição no seu ponto fixo. Com a representação das funções intervalares em função passo constante a partir da partição do intervalo [a0,b0], reescrevese o algoritmo da Análise Intervalar em Teoria de Domínios com a introdução do cálculo da aproximação da extensão canônica de g e com o comprimento do intervalo da partição tendendo para zero. Estende-se essa abordagem mais completa do estudo das equações integrais na resolução de problemas de valores iniciais e valor de contorno em equações diferenciais ordinárias e parciais. Uma vez que para uma pequena variação do campo vetorial v ou do valor inicial y0 da equação diferencial f ′(x) = v(x, f (x)) com a condição inicial f (x0) = y0, pode-se ter uma solução tão próxima da solução f da equação quanto possível, formaliza-se pela primeira vez em Teoria de Domínios um algoritmo na resolução do problema inverso em que, conhecendo a função f , determina-se uma equação diferencial ordinária com o cálculo de um campo vetorial v tal que o operador de Picard associado mapeia f tão próxima quanto possível a ela mesma. / We present a study in Domain Theory of integral equations of the form f (x) = h(x)+g Z b(x) a(x) g(x, y, f (y))dy for a0 ≤ a(x) ≤ b(x) ≤ b0 with h, a, b defined for x ∈ [a0,b0] and g defined for x, y ∈ [a0,b0], in which the right-hand side defines a contraction on the metric space of continuous realvalued functions on [a0,b0]. The starting point of this work is to revisit Interval Analysis in Domain Theory for the initial-value problem in ordinary differential equations where a solution is expressed as a fixed point of the Picard operator. With the set of real numbers interpreted as the interval domain, real-valued functions are extended to work in the space of interval-valued functions of the real variable domain. In particular, the vector field is extended in the second argument. Under these conditions, for the first time Fredholm and Volterra integral equations have solutions expressed as fixed points of a contraction mapping in terms of the splitting on interval-valued functions of the real variable domain. The measurement for interval-valued functions of the real variable domain is considered where we can asssess the convergence properties of the interval integral operator by means of the informatic derivative. The proposed techniques are applied to more general methods in ordinary differencial equations (ODEs) and partial differential equations (PDEs). For the first time, an algorithm is proposed to provide solutions to the inverse problem for Odinary Differential Equation where, given a function f , it is found a vector field v that defines a Picard operator which maps the solution f as close as possible to itself, such that the ODE f ′(x) = v(x, f (x)) admits f as either an exact or, as closely as desired, an approximate solution.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Neste trabalho, três técnicas para resolver numericamente problemas inversos de transporte de partículas neutras a uma velocidade para aplicações em engenharia nuclear são desenvolvidas. É fato conhecido que problemas diretos estacionários e monoenergéticos de transporte são caracterizados por estimar o fluxo de partículas como uma função-distribuição das variáveis independentes de espaço e de direção de movimento, quando os parâmetros materiais (seções de choque macroscópicas), a geometria, e o fluxo incidente nos contornos do domínio (condições de contorno), bem como a distribuição de fonte interior são conhecidos. Por outro lado, problemas inversos, neste trabalho, buscam estimativas para o fluxo incidente no contorno, ou a fonte interior, ou frações vazio em barras homogêneas. O modelo matemático usado tanto para os problemas diretos como para os problemas inversos é a equação de transporte independente do tempo, a uma velocidade, em geometria unidimensional e com o espalhamento linearmente anisotrópico na formulação de ordenadas discretas (SN). Nos problemas inversos de valor de contorno, dado o fluxo emergente em um extremo da barra, medido por um detector de nêutrons, por exemplo, buscamos uma estimativa precisa para o fluxo incidente no extremo oposto. Por outro lado, nos problemas inversos SN de fonte interior, buscamos uma estimativa precisa para a fonte armazenada no interior do domínio para fins de blindagem, sendo dado o fluxo emergente no contorno da barra. Além disso, nos problemas inversos SN de fração de vazio, dado o fluxo emergente em uma fronteira da barra devido ao fluxo incidente prescrito no extremo oposto, procuramos por uma estimativa precisa da fração de vazio no interior da barra, no contexto de ensaios não-destrutivos para aplicações na indústria. O código computacional desenvolvido neste trabalho apresenta o método espectronodal de malha grossa spectral Greens function (SGF) para os problemas diretos SN em geometria unidimensional para gerar soluções numéricas precisas para os três problemas inversos SN descritos acima. Para os problemas inversos SN de valor de contorno e de fonte interior, usamos a propriedade da proporcionalidade da fuga de partículas; ademais, para os problemas inversos SN de fração de vazio, oferecemos a técnica a qual nos referimos como o método físico da bissecção. Apresentamos resultados numéricos para ilustrar a precisão das três técnicas, conforme descrito nesta tese. / In this work, three techniques for numerically solving one-speed neutral particle inverse transport problems for nuclear engineering applications are developed. It is well known that direct steady-state monoenergetic transport problems are characterized by estimating the flux of particles as a distribution function of space and direction-of-motion independent variables, when the material parameters (cross sections), the geometry, and the incoming flux at the boundaries of the domain (boundary conditions), as well as the interior source distribution are known. Conversely, inverse problems, in this work, seek for estimates to the incident boundary flux, or interior source, or void fractions in homogeneous slabs. The mathematical model used for direct and inverse problems is the time-independent one-speed slab-geometry transport equation with linearly anisotropic scattering in the discrete ordinates (SN) formulation. In the boundary-value inverse problems, given the existing flux at one boundary of the slab, as measured by a neutron detector, for example, we seek for accurate estimate for the incident flux at the opposite boundary. On the other hand, in the interior source inverse SN problems, we seek for accurate estimate for the interior source stored within the slab for shielding purpose, given the exiting flux at the boundary of the slab. Furthermore, as with the void fraction inverse SN problems, given the exiting flux at one boundary of the slab due to prescribed incident flux at the opposite boundary, we seek for accurate estimate of the void fraction within the slab in the context of non-destructive testing applications in industry. The computer code developed in this work presents the coarse-mesh spectral Greens function (SGF) nodal method for direct SN problems in slab geometry to generate accurate numerical solutions to the three inverse SN problems described above. For the boundary-value and interior source inverse SN problems, we use the proportionality property of the leakage of particles; moreover, for the void fraction inverse SN problems, we offer the technique that we refer to as the physical bisection method. We present numerical results to illustrate the accuracy of the three techniques, as described in this dissertation.

CNPq / Este trabalho tem por objetivo apresentar uma nova proposta aos algoritmos de Demosaicing existentes, utilizando uma abordagem mais flexível quanto ao uso do Color filter array (CFA) em imagens coloridas de único sensor. O algoritmo proposto tem como base a estrutura de problemas inversos, cujo funcionamento utiliza um modelo de operação matriz-vetor que é adaptável ao CFA empregado. A partir deste conceito, o algoritmo trata o problema de Demosaicing como o de minimização de função custo, com um termo referente à dependência da estimativa com os dados obtidos e com o modelo de captura, o outro termo é relacionado aos conhecimentos observados em imagens que podem ser explorados para uma estimativa mais precisa, tal elemento é chamado de Prior. A proposta estabelecida tem como base algoritmos de regularização com foco na alta correlação presente entre os canais de cor (R, G e B), e na suavidade local de regiões uniformes, essa base formaliza o Prior empregado no trabalho. A minimização da proposta é atingida iterativamente através do IRLS-CG, que é a combinação de dois algoritmos de minimização eficientes, que apresenta rápidas respostas, e a capacidade de trabalhar com a norma L1 em conjunto com a norma L2. Com o intuito de atestar a qualidade do algoritmo proposto, foi elaborado um experimento em que o mesmo foi testado com diferentes CFAs e em situação com ruído gaussiano de 35dB e sem ruído algum em imagens da base de dados da Kodak, e os resultados comparados com algoritmos do estado-da-arte, no qual o desempenho da proposta apresentou resultados excelentes, inclusive em CFAs que destoam do padrão Bayer, que é o mais comumente usado na atualidade. / This research presents a new proposal to Demosaicing algorithms, using a more flexible approach to deal with the Color filter array (CFA) in single sensor color imaging. The proposed algorithm is structured in the inverse problems model, whose functions employ a CFA adaptive matrix-vector operational model. From this concept, the Demosaicing problem is treated as a cost function minimization with two terms, one referring to the dependence between the estimation and the data provided by the acquisition model, and other term related to features observed in images, which can be explored to form a more precise estimation, this last term is known as Prior. The established proposal is applied in regularization algorithms with focus on the high correlation among color channels (R, G, and B), and in the local smoothness of uniform regions. Both characteristics organize the Prior employed in this work. The minimization proposed is iteratively achieved through IRLS-CG, which is the combination of two efficient minimization algorithms, that presents quick responses, and the capacity to deal with L1 and L2 norm at the same time. The quality of the proposed algorithm is verified in an experiment in which varous CFA were used and a situation with 35dB gaussian noise and another one with no noise applied to the Kodak dataset, and the results were compared with state-of-the-art algorithms, in which the performance of the proposed Prior showed excellent results, including when the CFA is different from Bayer’s, which is the most commonly used pattern.

CNPq / Este trabalho tem por objetivo apresentar uma nova proposta aos algoritmos de Demosaicing existentes, utilizando uma abordagem mais flexível quanto ao uso do Color filter array (CFA) em imagens coloridas de único sensor. O algoritmo proposto tem como base a estrutura de problemas inversos, cujo funcionamento utiliza um modelo de operação matriz-vetor que é adaptável ao CFA empregado. A partir deste conceito, o algoritmo trata o problema de Demosaicing como o de minimização de função custo, com um termo referente à dependência da estimativa com os dados obtidos e com o modelo de captura, o outro termo é relacionado aos conhecimentos observados em imagens que podem ser explorados para uma estimativa mais precisa, tal elemento é chamado de Prior. A proposta estabelecida tem como base algoritmos de regularização com foco na alta correlação presente entre os canais de cor (R, G e B), e na suavidade local de regiões uniformes, essa base formaliza o Prior empregado no trabalho. A minimização da proposta é atingida iterativamente através do IRLS-CG, que é a combinação de dois algoritmos de minimização eficientes, que apresenta rápidas respostas, e a capacidade de trabalhar com a norma L1 em conjunto com a norma L2. Com o intuito de atestar a qualidade do algoritmo proposto, foi elaborado um experimento em que o mesmo foi testado com diferentes CFAs e em situação com ruído gaussiano de 35dB e sem ruído algum em imagens da base de dados da Kodak, e os resultados comparados com algoritmos do estado-da-arte, no qual o desempenho da proposta apresentou resultados excelentes, inclusive em CFAs que destoam do padrão Bayer, que é o mais comumente usado na atualidade. / This research presents a new proposal to Demosaicing algorithms, using a more flexible approach to deal with the Color filter array (CFA) in single sensor color imaging. The proposed algorithm is structured in the inverse problems model, whose functions employ a CFA adaptive matrix-vector operational model. From this concept, the Demosaicing problem is treated as a cost function minimization with two terms, one referring to the dependence between the estimation and the data provided by the acquisition model, and other term related to features observed in images, which can be explored to form a more precise estimation, this last term is known as Prior. The established proposal is applied in regularization algorithms with focus on the high correlation among color channels (R, G, and B), and in the local smoothness of uniform regions. Both characteristics organize the Prior employed in this work. The minimization proposed is iteratively achieved through IRLS-CG, which is the combination of two efficient minimization algorithms, that presents quick responses, and the capacity to deal with L1 and L2 norm at the same time. The quality of the proposed algorithm is verified in an experiment in which varous CFA were used and a situation with 35dB gaussian noise and another one with no noise applied to the Kodak dataset, and the results were compared with state-of-the-art algorithms, in which the performance of the proposed Prior showed excellent results, including when the CFA is different from Bayer’s, which is the most commonly used pattern.

Nesta tese mostraremos uma aplicação do método SGF, cf., spectral Greens function, para gerar fluxos angulares nas interfaces dos nodos na formulação de ordenadas discretas (SN) para equação de transporte de partículas neutras em geometria cartesiana unidimensional e a uma velocidade, usando quadraturas de alta ordem. Utilizando o método de malha grossa SGF, primeiramente determinamos as constantes arbitrárias da solução geral analítica das equações SN em cada nodo espacial. Então usamos a fórmula de quadratura angular para estimar uma expressão para o fluxo escalar de nêutrons e substituímos no termo de fonte de espalhamento isotrópico na equação de transporte. Resolvemos analiticamente a equação unidimensional de transporte de nêutrons com a fonte de espalhamento aproximada desta maneira e geramos os valores para o fluxo angular no interior de cada nodo espacial. Como o método SGF gera soluções numéricas de malha grossa completamente livres de erro de truncamento espacial, esperamos que o esquema de reconstrução analítica proposto tenha alta precisão para os fluxos angulares, considerando as condições de continuidade nas interfaces dos nodos do domínio espacial. Esta técnica caracteriza um problema inverso, pois a partir da solução de malha grossa do método SGF, podemos reconstruir o fluxo angular de nêutrons em qualquer ponto do domínio. / We offer in this work an application of the SGF method, cf., spectral Greens function, to generate the angular fluxes at the region interfaces of multilayer slabs, using high-order angular quadrature sets in the one-speed discrete ordinates (SN) formulation of the neutral particle transport equation. Using the SGF coarse-mesh numerical solution, we first determine the arbitrary constants of the analytical general solution of the SN equations within each spatial node. Then, we use the angular quadrature formula to estimate the expression of the scalar flux distribution, that we substitute into the isotropic scattering source term of the transport equation. We solve analytically the slab-geometry transport equation, with the scattering source so approximated, in order to generate the angular flux profile within each spatial node. As the SGF method generates coarse-mesh numerical solution, which is completely free from spatial truncation errors, we expect that the offered approximate analytical reconstruction scheme be accurate enough for the localized angular flux distribution, considering the node interface continuity conditions within the domain. This technique is thought of as an inverse problem since from the SGF coarse-mesh nodal solution we are able to reconstruct the angular flux profile at any point inside the domain.