Fractional-order TV-L 2 model for image denoising
|
|
- Noah Merritt
- 6 years ago
- Views:
Transcription
1 Cent. Eur. J. Phys. 11(10) DOI: /s Central European Journal of Physics Fractional-order TV-L 2 model for image denoising Research Article Dali Chen 1, Shenshen Sun 1, Congrong Zhang 1, YangQuan Chen 2, Dingyu Xue 1 1 Information Science and Engineering, Northeastern University, Wenhua Road 3-11, Heping Districe Shenyang, Liaoning, CHINA 2 MESA Lab, University of California, Merced, 5200 North Lake Road Merced, CA 95343, USA Received 31 January 2013; accepted 26 April 2013 Abstract: This paper proposes a new fractional order total variation (TV) denoising method, which provides a much more elegant and effective way of treating problems of the algorithm implementation, ill-posed inverse, regularization parameter selection and blocky effect. Two fractional order TV-L 2 models are constructed for image denoising. The majorization-minimization (MM) algorithm is used to decompose these two complex fractional TV optimization problems into a set of linear optimization problems which can be solved by the conjugate gradient algorithm. The final adaptive numerical procedure is given. Finally, we report experimental results which show that the proposed methodology avoids the blocky effect and achieves state-of-the-art performance. In addition, two medical image processing experiments are presented to demonstrate the validity of the proposed methodology. PACS (2008): y, Xx, Tt Keywords: image denoising fractional calculus total variation majorization-minimization algorithm Versita sp. z o.o. 1. Introduction Since the work of Rudin-Osher-Fatemi [1], the TV-L 2 model has been used in many applications such as image restoration [2], image decomposition [3], image decompression [4] and image segmentation [5]. The typical TV-L 2 model, for example ROF model, has been proved to be able to achieve a good trade-off between edge preservation and noise removal [1]. However, it tends to produce the so-called blocky (staircase) effects on the images because it favors a piecewise constant solution in bounded chendali@ise.neu.edu.cn variation (BV) space [6]. In order to deal with blocky effects, the modification of TV-L 2 model, which generalizes the differential order in regularization term, has aroused the more and more attentions of numerous scholars. There are two types of generalization of the differentiation in regularization term. The first one deals with higherorder differentiation. In 2000, You and Kaveh introduced a fourth-order partial differential equation (PDE)-based denoising model in which the regularized solution is obtained by solving the minimization of potential function of second-order derivative of the image [7]. Although this model reduces the blocky effects, it tends to cause the sign of uplifting effect and formation of artifacts around edges. For this problem, some improved fourth-order PDE models were proposed by replacing the Laplacian operator of 1414
2 Dali Chen, Shenshen Sun, Congrong Zhang, YangQuan Chen, Dingyu Xue diffusivity function by the different operators such as the gradient [8]. In this paper, our interest focus on the second generalization which deals with fractional order differentiation. Fractional calculus is a rapidly growing mathematical discipline, which provides an important tool for nonlocal field theories [9]. Recently, fractional calculus has been greatly studied in computer vision [10 14]. The main reason for this development is the expectation that the use of this theory will lead to a much more elegant and effective way of treating problems of blocky effect and detail information protection. Specially, the fractional order total variation (TV) models play an important role for image denoising, inpainting and motion estimation [15 19]. However, the algorithms to compute these TV-based models are developed on the continuous domain. Some of them use the methods based on duality and others solve the associated Euler-Lagrange equation, which is a non-linear partial differential equation (PDE). Note that these methods have to be discretized for the application on digital image, and the discretization method thus needs to be considered. The majorization-minimization (MM) algorithm is an iterative optimization method which is formulated in the discrete domain. The main idea of MM is to solve a sequence of simpler optimization problems instead of a complex one [20]. Recently, the majorization-minimization algorithm has been used to solve the TV deconvolution and denoising optimization problems [21, 22]. In our paper, a new image denoising method is proposed, in which the MM algorithm is used to solve the fractional order TV-L 2 (TV α -L 2 ) denoising model. The proposed denoising method is able to avoid the blocky effect and achieve state-of-the-art performance. This paper is organized as follows. Section 1 introduces prior work and our motivation, and some preliminary knowledges used in this paper are presented in Section 2. In Section 3, we firstly introduce two fractional order TV-L 2 denoising models, which are based on the different TV regularizers. Then, the majorizors of two fractional order TV regularizers are given in one uniform formula. Finally, the numerical procedure is shown for the readers. Experimental evaluation is presented in Section 4 and the paper is concluded in Section Preliminaries 2.1. Majorization-Minimization algorithm The Majorization-Minimization algorithm is an iterative optimization method which decomposes an arbitrarily complex optimization problem into a set of simpler ones. We consider the optimization problem described by arg min E(x). (1) x RN Instead of solving Eq. (1) directly, the Majorization- Minimization algorithm solves a sequence of simpler optimization problems x k+1 = arg min x R N H k(x), (2) where the functions H k (x), k = 0, 1, 2,... are a sequence of majorizors of the cost function E(x), which should satisfy the following two conditions: H k (x) E(x), H k (x k ) = E(x k ). In addition, the Majorization-Minimization algorithm has a useful property [22]. Consider E(x) = E 1 (x) + E 2 (x), and Q k (x) is the majorizor of E 1 (x). Then the majorizor of E(x) can be obtained by the function Q k (x) + E 2 (x). x 2.2. Matrix approximate method An efficient matrix approximate method for fractional order derivative was proposed by [23], which leads to significant simplification of the numerical solution. For easy reading, we reproduce the relevant results. Let us assume that the given signal f is sampled from its continuous version at an uniformly spaced grid size of h. Thereby, the discrete form f i = f(i h), for i = 0, 1,..., m. In image processing applications, the grid size h is one, so we set h = 1 for easy description. From Grünwald-Letnikov fractional derivative definition [24], the discrete formula of the fractional order derivative of the digital signal can be defined by the following formula D α x m k=0 w (α) k f i k, (3) ( ) where w (α) α k = ( 1) k represent the coefficients of the k polynomial (1 z) α. The coefficients can also be obtained recursively from w (α) 0 = 1, w (α) k = (1 α + 1 k ) w (α) k 1, k = 1, 2,.... (4) Using the matrix approximate method, we can rewrite Eq. (4) by the following form D α x f ψ f, (5) 1415
3 Fractional-order TV-L 2 model for image denoising where f = [f 0, f 1,..., f m ] T and ψ is a matrix defined by the following formula w (α) w (α) 1 w (α) ψ = (6).. w m (α) w (α) m 1... w (α) 0 When α = 1, the matrix ψ is a sparse banded matrix which consists of only two diagonals, the main diagonal and one lower diagonal. So, only two points are used to measure the gradient information. When α is not integer, the matrix ψ is a lower triangular matrix and all the points before k are used to calculus the fractional order derivative of the kth point. This long-term memory is an important characteristic of fractional differentiation and also is a prominent difference between fractional differentiation and integer order differentiation. Because of the long-term memory, fractional differentiation becames an important tool for modelling characteristic phenomena in various application. 3. Fractional-order TV-L 2 model 3.1. Model description Let f(x, y) = u(x, y) + v(x, y) denote the observed noisy image, where (x, y) T denotes the location with a rectangular image domain Ω R 2 and v is white Gaussian noise. TV-L 2 denoising model estimates the desired clean image u(x, y) by solving the following finite-dimensional optimization problem û = arg min u { E(u) = f u λtv(u)}, (7) where f u 2 2 is the data fidelity term, TV(u) is the regularization term, u ν is ν-norm of u, and λ is regularization parameter which controls the degree of smoothing. In this paper, we consider a fractional order TV p -L 2 (TV α p- L 2 )model, defined as û = arg min u { E(u) = f u λtvα p(u) }, (8) where p = 1 or 2. For better description, the discrete definition of the regularization term TV α p(u) is used, which is obtained by the following formula TV α p(u) = n Ω D α u(n) p, (9) where D α is the fractional-order derivative operator. When p = 1, D α u(n) 1 = D α hu(n) + D α v u(n) ; when p = 2, D α u(n) 2 = (D α hu(n)) 2 + (D α v u(n)) 2. D α h and D α v are linear operators corresponding to horizontal and vertical fractional order derivative Majorizor of TV α p(u) The key to apply Majorization-Minimization algorithm to the TV α p-l 2 denoising model is to find the majorizor of TV α p(u). Using the fact that t t2 2 t k + t k 2, t, t k R, (10) it follows that the function Q (k,1) (u) defined by Q (k,1) (u) = 1 2 TVα 1 (u k ) + 1 ( (Dα h u(n))2 2 D α n Ω h u k(n) + (Dα v u(n)) 2 Dv α u k (n) ) (11) satisfies Q (k,1) (u) TV α 1 (u) and Q (k,1) (u k ) = TV α 1 (u k ). Therefore, Q (k,1) (u) is the majorizor of TV α 1 (u). For easy description, we let D α h and D α v denote matrices such that D α hu and D α v u respectively yield the fractional order horizontal and vertical derivative. Based on this, we can rewrite Eq. (11) by the following compact formula Q (k,1) (u) = 1 2 TVα 1 (u k ) ut (D α ) T V 1 (k,1)d α u = C α 1 (u k ) ut (D α ) T V 1 (k,1)d α u, (12) where D α = [(D α h) T, (D α v ) T ] T, V (k,1) = diag(λ h k, Λv k ), Λh k = diag( D α hu k ) and Λ v k = diag( Dα v u k ). In the same way, we define the function ( Q (k,2) (u) =TV α 2 (u k ) + 1 (D α hu(n)) 2 (D α hu k (n)) 2 2 n Ω (D α hu k (n)) 2 + (D α v u k (n)) 2 ) + (Dα v u(n)) 2 (D α v u k (n)) 2. (13) (D α hu k (n)) 2 + (D α v u k (n)) 2 Since the inequality t t k + (t t k ) 2 t k for any t 0, the function Q (k,2) (u) satisfies Q (k,2) (u) TV α 2 (u) and Q (k,2) (u k ) = TV α 2 (u k ). Therefore, Q (k,2) (u) is the majorizor of TV α 2 (u) and has a compact form described by Q (k,2) (u) = TV α 2 (u k ) ut (D α ) T V 1 (k,2)d α u = C α 2 (u k ) ut (D α ) T V 1 (k,2)d α u, (14) 1416
4 Dali Chen, Shenshen Sun, Congrong Zhang, YangQuan Chen, Dingyu Xue where V (k,2) = diag(v k, v k ) and v k = diag( (D α hu k ) 2 + (D α v u k ) 2 ). Integrating Eq. (12) and Eq. (14) into one uniform formula, we can obtain the majorizor of TV α p(u) by the following function Q (k,p) (u) = C α p(u k ) ut (D α ) T V 1 (k,p)d α u. (15) 3.3. Numerical scheme This section gives details on the employed numerical procedure for the proposed TV α p(u)-l 2 denoising approach. Based on the property of Majorization-Minimization algorithm, the majorizor for the cost function E(u) is constructed by adding the data term f u 2 2 to Eq. (15), which is given by the function H (k,p) (u) = f u ut (D α ) T V 1 (k,p)d α u+c α p(u k ). (16) Since C α p(u k ) is a constant independent of u, it is irrelevant for the optimization problem. We thus estimate u by solving a sequence of optimization problems u k+1 = arg min u { f u λ 2 ut (D α ) T V 1 (k,p)d α u }. (17) Since Eq. (16) is a quadratic function, an explicit solution to Eq. (17) leads to a linear system ( I + λ(d α ) T V 1 (k,p)d α) u k+1 = f, (18) which can be solved by the conjugate gradient algorithm. However, some values of D α u p will go to zero as the iterations progress, and therefore the matrix V(k,p) 1 will go to infinity. This makes the solution of Eq. (18) more difficult. We tackle this difficulty by adopting the spatial smoothing strategy which has been largely used to overcome the ill-posed problem. It is common to smooth the image prior to differentiation, e.g. by convolving the updated image with the Gaussian kernel K σ (x, y) of standard deviation σ, u (x, y) = (K σ u)(x, y). The regularization parameter λ in Eq. (18) is related to the noise statistics and controls the degree of filtering of the obtained solution. When λ is too large, details and texture can be over smoothed; when λ is too small, it yields overly oscillatory estimates owing to either noise or discontinuities. Therefore, the selection of the regularization parameter is a critical issue to which much attention has been devoted. Many popular approaches have been proposed in the regularization and Bayesian frameworks. In this paper, the method proposed in [22] is used to adjust λ, λ p k = MNσ 2 /2TV α p(u k ). Here M and N are respectively the row and column of given image, and σ is the standard deviation of noise. To summarize, our entire noise removal algorithm in a form of a pseudo-code is done in following 1. Let the input image be f and set k = 0 and u k = f. Initialize K σ, σ, M, N and iteration number K iter. 2. Set u k = K σ u k and λ p k = MNσ 2 /2TV α p(u k ). 3. Compute A k = ( I + λ p k (Dα ) T V 1 (k,p) Dα) according to Eq. (12) to (15). 4. Compute u k+1 by using the conjugate gradient algorithm to solve the linear system A k u k+1 = f. 5. If k = K iter, stop; else, set k = k +1 and go to step Experiments and analysis 4.1. Restraint of block effect In this section, some experiments are given to assess the capability of reducing block effect of our proposed fractional order TV-L 2 model. Consider a one-dimension signal, y(t) = 4 sin(2πt) + 8 sin(3πt), which is contaminated by the additive white Gaussian noise (AWGN) with standard deviation (SD) of 0.5. The corrupted signal is shown in Fig. 1(a). We use the traditional TV-L 2 model, in which the differential order of regular term is one, to process the contaminated signal and the result at the 200th iteration is shown in Fig. 1(b). For comparison, the signal processed by the proposed fractional order TV-L 2 model with α = 2.8 is shown in Fig. 1(c). Note that the TV α 1 and TVα 2 have the same form when the given signal is one dimension signal. In order to ensure the fair comparison, the same iteration number K iter = 100 and regular parameter λ = 10 are used in this experiment. The difference between the traditional TV-L 2 model and proposed fractional order TV-L 2 model is obvious: while the traditional TV-L 2 model approximates the observed signal with a step signal, the proposed fractional order model with a piecewise planar signal which looks more natural and does not produce false edges. In conclusion, the proposed fractional order TV-L 2 model can reduce blocky effect effectively comparing with the traditional TV-L 2 model Analysis of denoising performance The aim of this section is to analyze the denoising performance of the proposed fractional order TV p -L 2 model. For 1417
5 Fractional-order TV-L 2 model for image denoising Figure 1. Noisy signal (a) with 0.5 SD Gaussian noise denoised with the traditional TV-L 2 model (b) and fractional order TV-L 2 model (c). Table 1. PSNR quantitative comparison among five denoising models. determine the quality of a processed image. calculated by the following formula It can be Image SD IP-M F-O-PDE IF-O-PDE ROF TV α 2 -L Barbara Lena Peppers this purpose, three famous test images, named Barbara, Lena and Peppers, are shown in Fig. 2, which have been corrupted by the additive white Gaussian noise (AWGN) with standard deviation (SD) of 10, 20 and 30 respectively. Five diffusion models are used to process the given noisy images, which are the improved Perona and Malik (IP-M) model [25], fourth order (F-O) PDE model [7], improved fourth order (IF-O) PDE model[8], ROF model [1] and proposed fractional order TV 2 -L 2 (TV α 2 - L 2 ) model. In this experiment, we set σ =5and t =0.25 for IP-M, F-O-PDE and IF-O-PDE models, and we set λ =0.07, τ =0.1 and K = 20 for ROF model. These parameters are able to ensure the best denoising performance of the corresponding denoising model. For TV α 2 -L2 model, we set α =1.8. The denoising results are shown in Table 1. In the table, the first column shows three given images and the second column shows their noisy images corrupted by the additive white Gaussian noise (AWGN) with standard deviation (SD) of 10, 20 and 30 respectively. In order to quantify the denoised image, we consider the peak-signal-to-noise ratio (PSNR), which has been largely used in literature and commonly applied to VL 2 PSNR = 10 log MN 10 N M, (19) i=1 j=1 (u(i, j) û(i, j))2 where u is the original image, û is the denoised image and V L is the maximal gray level of the image. The PSNRs of the denoised images processed by the different denoising models are listed under the corresponding denoising model. The bigger PSNR is, the better denoising performance is. For easy observation, we use bold font for the biggest PSNR values in the table. From the table, it is obvious that the PSNR values of our proposed model are bigger than that of the other four models, so we can conclude that our TV α 2 -L2 model outperforms the other models. In order to further verify the denoising performance of our proposed model, Fig. 3 shows the denoised images processed by four denoising models. Here, we use the different iteration number K iter in the different denoising model. For TV α 2 -L2 model, we set K iter = 10 because the method converges at 10 iterations. Similarly, we set K iter = 100 for the ROF model to ensure the converged result. In order to obtain the best result, we set K iter = 13 for the IP-M model and K iter = 10 for the IF-O PDE model. The first figure is the famous image "Goldhill" corrupted by the additive white Gaussian noise with standard deviation (SD) of 10. In order to show the clear comparison, the partial enlarged view of the noisy image is shown in the second figure. The denoised result of our proposed model is shown in the right-top figure. The second row shows the denoised images processed by IP-M model, IF-O PDE model and ROF model. In the left-bottom figure, although the detail information is preserved, there are a lot of noises unremoved. From the middle-bottom figure, it is obvious that the denoised image produces many false edges, and the result in the right-bottom figure looks blocky. Only the result of our model shown in the right-top figure looks 1418
6 Dali Chen, Shenshen Sun, Congrong Zhang, YangQuan Chen, Dingyu Xue Figure 2. Barbara (a), Lena (b) and Peppers (c). Figure 3. Goldhill with 10 SD Gaussian noise (a) and its partial view (b) processed by TVα2 -L2 model (c), IP-M model (d), IF-O PDE model (e) and ROF model (f). natural and does not produce false edges. So we conclude that our proposed fractional order model is able to achieve a good trade-off between edge preservation and noise removal. Table 2 shows the PSNR values of denoised images processed by the TVα1 -L2 and TVα2 -L2 model. The first row is the given image and the second row is the standard deviation (SD) of additive white Gaussian noise of noisy image. The PSNR values obtained by the TVα2 -L2 model are listed in the third row and the values obtained by the TVα1 -L2 model are listed in the last row. Here we set α = 1.6. It is obvious that the PSNR values of TVα2 -L2 model are bigger than that of TVα1 -L2 model. So we conclude that the denoising performance of TVα2 -L2 model is better. From the other α values, we can get the same conclusion. 1419
7 Fractional-order TV-L 2 model for image denoising Table 2. PSNR quantitative comparison between TV α 1 -L2 model and TV α 2 -L2 model Barbara Lena Peppers Goldhill Model TV L TV L Figure 4. Relation between the iteration number and PSNR value. In this experiment, the proposed TV α 2 -L2 model is used to process the image "Lena" corrupted by the additive white Gaussian noise with standard deviation (SD) of 10, and the PSNR value is recorded at each iteration. Fig. 4 shows the relation between the iteration number and PSNR value. The horizontal axis of this figure is iteration number and the vertical axis is PSNR value. The different colorful curve denotes the result obtained from the different α, for example, the red curve is the result obtained by the TV α 2 -L2 model with α =1.4. From the figure, we can see that the PSNR values are stable after 15 iterations, which indicates our proposed TV α 2 -L2 denoising algorithm is stable. In addition, it can be seen that the steady-state PSNR value of our proposed fractional order TV-L 2 model is bigger than that of the traditional TV-L 2 model, which shows that the fractional order TV-L 2 model has better denoising performance. In addition, in order to decide the value of the fractional order α, we study the relation between the PSNR and α. The Fig. 5 shows the relation between the PSNR and α on Lena and Peppers image corrupted by the additive white Gaussian noise with standard deviation (SD) of 10. The left image is the result of Lena image and the right one is the result of Peppers image. From the figure, we can obtain the following conclusions. Firstly, the PSNR reaches a maximum between α =1and α =2. Secondly, the PSNR at α =1is lower than PSNR at α>1, which is owing to the blocky effect. Finally, the PSNR decreases rapidly as α tends to zero and when α =0, Figure 5. Relation between the fractional order α and PSNR on Lena (left) and Peppers (right) image corrupted by the additive white Gaussian noise with standard deviation (SD) of 10. the PSNR reaches a minimum. According to these facts, we can select the fractional order α between α =1.4 and α = Lung nodule segmentation experiment Image denoising technology not only can improve the observability of the image, but also can increase the accuracy of some postprocessing technologies. In this experiment, the lung nodule CT image is used as the test image, which is downloaded from the LIDC medical imaging database [26]. The given CT image is obtained from Siemens low-dose spiral CT, whose peak voltage ranges from 120 kvp to 140 kvp and tube current ranges from 1420
8 Dali Chen, Shenshen Sun, Congrong Zhang, YangQuan Chen, Dingyu Xue Figure 6. Lung nodule CT image with Gaussian noise and the partial view of the segmented results of the undenoised image and the images denoised by the TV 1 2 -L2 and TV L2 denoising methods. Figure 7. PET image with Gaussian noise and the partial view of the segmented results of the undenoised image and the images denoised by the TV 1 2 -L2 and TV L2 denoising methods. 40 mas to 388 mas. In order to assess our denoising method, the lung nodule CT image is contaminated by the additive white Gaussian noise and the typical regiongrowing segmentation method [27] is used to segment the lung nodule from the noisy image and denoised image. The results are shown in Fig. 6. The first figure shows the given noisy lung nodule CT image. For easy observation, the partial enlarged view of the segmented result of the noisy image is shown in the second figure, from which we can see that the lung nodule cannot be segmented. The third and forth figures show the segmented results of the images denoised by the TV 1 2 -L2 and TV L2 denoising methods, respectively. Although these two segmented results both show the right location of lung nodule, the detected edge of TV L2 method is more accurate than that of TV 1 2 -L2 method. This demonstrates that the proposed fractional order TV-L 2 denoising method is helpful for improving the accuracy of the postprocessing technologies in the lung nodule CT image segmentation Cardiac muscular segmentation experiment In the medical image processing, the segment of the cardiac muscular tissue is an important step. The purpose of this experiment is to assess the effect of our method on the the cardiac muscular segmentation. In this experiment, the test image is the cardiac muscular PET image downloaded from the DICOM7 medical image database ( The additive white Gaussian noise is added in the test image. We use the typical iteration threshold segmentation algorithm [27] to segment the cardiac muscular tissue from the noisy image and the denoised images respectively. For easy observation, the pseudo-color technology is used on the segmented images, and the results are shown in Fig. 7. The first figure shows the noisy image and the red part of which is the cardiac muscular tissue. The partial enlarged view of the segmented result of the noisy image is shown in the second figure, in which some noises remain. The third and forth figures show the segmented results of the images denoised by the TV 1 2 -L2 and TV L2 denoising methods, respectively. Although the segmented results both are satisfactory, the result of TV L2 denoising method looks much more smooth than that of TV 1 2 -L2 method. Therefore, we can conclude that the proposed fractional order TV-L 2 denoising method is helpful for improving the accuracy of the postprocessing technologies in the cardiac muscular PET image segmentation. 1421
9 Fractional-order TV-L 2 model for image denoising 5. Conclusion A new adaptive fractional order TV-based denoising method was proposed in this paper. The main contributions are as following: (1) two fractional order TV-L 2 (TV α p-l 2 ) models are constructed for image denoising; (2) the majorization-minimization algorithm was used to solve the TV α p-l 2 model, which provides an effective tool to solve fractional TV optimization problem; (3) the majorizors of two fractional order TV regularizers are obtained in one uniform formula; (4) the experiments demonstrate that the proposed methodology is able to avoid the blocky effect and achieve state-of-the-art performance. Future works involve extending the proposed methodology to the other total variational models such as optical flow model. Acknowledgments This work was supported by National Natural Science Foundation of China (No ), Scientific Research Fund of Liaoning Provincial Education Department (L ) and Fundamental Research Funds for the Central Universities (N ). The authors would like to thank Mário A. T. Figueiredo, Ivan Selesnick and Igor Podlubny for their outstanding prior work. References [1] L. Ruding, S. Osher, E. Fatemi, Physica D 60, 259 (1992) [2] J. F. Aujol, J. Math. Imaging Vis. 34, 307 (2009) [3] C. Vogel, M. Oman, IEEE T. Image Process. 7, 813 (1998) [4] F. Alter, S. Durand, J. Froment, J. Math. Imaging Vis. 23, 199 (2005) [5] F. Li, C. Shen, C. Li, J. Math. Imaging Vis. 37, 98 (2010) [6] J. Zhang, Z. Wei, L. Xiao, J. Math. Imaging Vis. 43, 39 (2012) [7] Y. L. You, M. Kaveh, IEEE T. Image Process. 9, 1723 (2000) [8] M. Hajiaboli, IPSJ Transactions on Computer Vision and Application 2, 94 (2010) [9] R. Herrmann, Fractional Calculus: An Introduction for Physicists (World Scientific, New Jersey, 2011) [10] S. C. Liu, S. Chang, IEEE T. Image Process. 6, 1176 (1997) [11] S. Didas, B. Burgeth, A. Imiya, J. Weickert, Scale Space and PDE Methods in Computer Vision 3459, 13 (2005) [12] B. Ninness, IEEE T. Inf. Theory 44, 32 (1998) [13] I. Petras, D. Sierociuk, I. Podlubny, IEEE T. Signal Proces. 60, 5561 (2012) [14] Y. F. Pu, J. L. Zhou, X. Yuan, IEEE T. Image Process. 19, 491 (2010) [15] B. Jian, X. C. Feng, IEEE T. Image Process. 16, 2492 (2007) [16] P. Guidotti, J. V. Lambers, J. Math. Imaging Vis. 33, 25 (2009) [17] E. Cuesta, M. Kirane, S. A. Malik, Signal Process. 92, 553 (2012) [18] M. Janev, S. Pilipović, T. Atanacković, R. Obradović, N. Ralević, Math. Comput. Model. 54, 729 (2011) [19] D. Chen, H. Sheng, Y. Q. Chen, D. Y. Xue, Phil. Trans. R. Soc. A., DOI: /rsta [20] D. Hunter, K. Lange, The American Statistician 58, 30 (2004) [21] M. A. T. Figueiredo, J. M. Bioucas Dias, R. D. Nowak, IEEE T. Image Process. 16, 2980 (2007) [22] J. P. Oliveira, J. M. Bioucas Dias, M. A. T. Figueiredo, Signal Process. 89, 1683 (2009) [23] I. Podlubny, Fractional Calculus and Applied Analysis 3, 359 (2000) [24] I. Podlubny, Fractional Differential Equations (Academic Press, San Diego, 1999) [25] P. Perona, J. Malik, IEEE T. Pattern Anal. 12, 629 (1990) [26] S. G. Armato, et al., Med. Phys. 38, 915 (2011) [27] R. C. Gonzalez, R. E. Woods, Digital Image Processing, 2nd edition (Addison-Wesley, Massachusetts, 1992) 1422
ENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT
ENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT PRASHANT ATHAVALE Abstract. Digital images are can be realized as L 2 (R 2 objects. Noise is introduced in a digital image due to various reasons.
More informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
More informationA Nonlocal p-laplacian Equation With Variable Exponent For Image Restoration
A Nonlocal p-laplacian Equation With Variable Exponent For Image Restoration EST Essaouira-Cadi Ayyad University SADIK Khadija Work with : Lamia ZIAD Supervised by : Fahd KARAMI Driss MESKINE Premier congrès
More informationMathematical Problems in Image Processing
Gilles Aubert Pierre Kornprobst Mathematical Problems in Image Processing Partial Differential Equations and the Calculus of Variations Second Edition Springer Foreword Preface to the Second Edition Preface
More informationConvex Hodge Decomposition of Image Flows
Convex Hodge Decomposition of Image Flows Jing Yuan 1, Gabriele Steidl 2, Christoph Schnörr 1 1 Image and Pattern Analysis Group, Heidelberg Collaboratory for Image Processing, University of Heidelberg,
More informationLINEARIZED BREGMAN ITERATIONS FOR FRAME-BASED IMAGE DEBLURRING
LINEARIZED BREGMAN ITERATIONS FOR FRAME-BASED IMAGE DEBLURRING JIAN-FENG CAI, STANLEY OSHER, AND ZUOWEI SHEN Abstract. Real images usually have sparse approximations under some tight frame systems derived
More informationIMAGE RESTORATION: TOTAL VARIATION, WAVELET FRAMES, AND BEYOND
IMAGE RESTORATION: TOTAL VARIATION, WAVELET FRAMES, AND BEYOND JIAN-FENG CAI, BIN DONG, STANLEY OSHER, AND ZUOWEI SHEN Abstract. The variational techniques (e.g., the total variation based method []) are
More informationPoisson Image Denoising Using Best Linear Prediction: A Post-processing Framework
Poisson Image Denoising Using Best Linear Prediction: A Post-processing Framework Milad Niknejad, Mário A.T. Figueiredo Instituto de Telecomunicações, and Instituto Superior Técnico, Universidade de Lisboa,
More informationErkut Erdem. Hacettepe University February 24 th, Linear Diffusion 1. 2 Appendix - The Calculus of Variations 5.
LINEAR DIFFUSION Erkut Erdem Hacettepe University February 24 th, 2012 CONTENTS 1 Linear Diffusion 1 2 Appendix - The Calculus of Variations 5 References 6 1 LINEAR DIFFUSION The linear diffusion (heat)
More informationNovel integro-differential equations in image processing and its applications
Novel integro-differential equations in image processing and its applications Prashant Athavale a and Eitan Tadmor b a Institute of Pure and Applied Mathematics, University of California, Los Angeles,
More informationEnergy-Based Image Simplification with Nonlocal Data and Smoothness Terms
Energy-Based Image Simplification with Nonlocal Data and Smoothness Terms Stephan Didas 1, Pavel Mrázek 2, and Joachim Weickert 1 1 Mathematical Image Analysis Group Faculty of Mathematics and Computer
More informationTWO-PHASE APPROACH FOR DEBLURRING IMAGES CORRUPTED BY IMPULSE PLUS GAUSSIAN NOISE. Jian-Feng Cai. Raymond H. Chan. Mila Nikolova
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume 00, Number 0, Xxxx XXXX pp. 000 000 TWO-PHASE APPROACH FOR DEBLURRING IMAGES CORRUPTED BY IMPULSE PLUS GAUSSIAN NOISE Jian-Feng
More informationRIGOROUS MATHEMATICAL INVESTIGATION OF A NONLINEAR ANISOTROPIC DIFFUSION-BASED IMAGE RESTORATION MODEL
Electronic Journal of Differential Equations, Vol. 2014 (2014), No. 129, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu RIGOROUS MATHEMATICAL
More informationNON-LINEAR DIFFUSION FILTERING
NON-LINEAR DIFFUSION FILTERING Chalmers University of Technology Page 1 Summary Introduction Linear vs Nonlinear Diffusion Non-Linear Diffusion Theory Applications Implementation References Page 2 Introduction
More informationEE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6)
EE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6) Gordon Wetzstein gordon.wetzstein@stanford.edu This document serves as a supplement to the material discussed in
More informationVariational Methods in Image Denoising
Variational Methods in Image Denoising Jamylle Carter Postdoctoral Fellow Mathematical Sciences Research Institute (MSRI) MSRI Workshop for Women in Mathematics: Introduction to Image Analysis 22 January
More informationConvex Hodge Decomposition and Regularization of Image Flows
Convex Hodge Decomposition and Regularization of Image Flows Jing Yuan, Christoph Schnörr, Gabriele Steidl April 14, 2008 Abstract The total variation (TV) measure is a key concept in the field of variational
More informationPOISSON noise, also known as photon noise, is a basic
IEEE SIGNAL PROCESSING LETTERS, VOL. N, NO. N, JUNE 2016 1 A fast and effective method for a Poisson denoising model with total variation Wei Wang and Chuanjiang He arxiv:1609.05035v1 [math.oc] 16 Sep
More informationTensor-based Image Diffusions Derived from Generalizations of the Total Variation and Beltrami Functionals
Generalizations of the Total Variation and Beltrami Functionals 29 September 2010 International Conference on Image Processing 2010, Hong Kong Anastasios Roussos and Petros Maragos Computer Vision, Speech
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationA Variational Approach to Reconstructing Images Corrupted by Poisson Noise
J Math Imaging Vis c 27 Springer Science + Business Media, LLC. Manufactured in The Netherlands. DOI: 1.7/s1851-7-652-y A Variational Approach to Reconstructing Images Corrupted by Poisson Noise TRIET
More informationRecent developments on sparse representation
Recent developments on sparse representation Zeng Tieyong Department of Mathematics, Hong Kong Baptist University Email: zeng@hkbu.edu.hk Hong Kong Baptist University Dec. 8, 2008 First Previous Next Last
More informationPartial Differential Equations and Image Processing. Eshed Ohn-Bar
Partial Differential Equations and Image Processing Eshed Ohn-Bar OBJECTIVES In this presentation you will 1) Learn what partial differential equations are and where do they arise 2) Learn how to discretize
More informationTRACKING SOLUTIONS OF TIME VARYING LINEAR INVERSE PROBLEMS
TRACKING SOLUTIONS OF TIME VARYING LINEAR INVERSE PROBLEMS Martin Kleinsteuber and Simon Hawe Department of Electrical Engineering and Information Technology, Technische Universität München, München, Arcistraße
More informationBERNSTEIN FILTER: A NEW SOLVER FOR MEAN CURVATURE REGULARIZED MODELS. Yuanhao Gong. National University of Singapore
BERNSTEIN FILTER: A NEW SOLVER FOR MEAN CURVATURE REGULARIZED MODELS Yuanhao Gong National University of Singapore ABSTRACT The mean curvature has been shown a proper regularization in various ill-posed
More informationSPARSE SIGNAL RESTORATION. 1. Introduction
SPARSE SIGNAL RESTORATION IVAN W. SELESNICK 1. Introduction These notes describe an approach for the restoration of degraded signals using sparsity. This approach, which has become quite popular, is useful
More informationPDE-based image restoration, I: Anti-staircasing and anti-diffusion
PDE-based image restoration, I: Anti-staircasing and anti-diffusion Kisee Joo and Seongjai Kim May 16, 2003 Abstract This article is concerned with simulation issues arising in the PDE-based image restoration
More informationAdaptive Corrected Procedure for TVL1 Image Deblurring under Impulsive Noise
Adaptive Corrected Procedure for TVL1 Image Deblurring under Impulsive Noise Minru Bai(x T) College of Mathematics and Econometrics Hunan University Joint work with Xiongjun Zhang, Qianqian Shao June 30,
More informationTotal Variation Image Edge Detection
Total Variation Image Edge Detection PETER NDAJAH Graduate School of Science and Technology, Niigata University, 8050, Ikarashi 2-no-cho, Nishi-ku, Niigata, 950-28, JAPAN ndajah@telecom0.eng.niigata-u.ac.jp
More informationTotal Variation Theory and Its Applications
Total Variation Theory and Its Applications 2nd UCC Annual Research Conference, Kingston, Jamaica Peter Ndajah University of the Commonwealth Caribbean, Kingston, Jamaica September 27, 2018 Peter Ndajah
More informationIntroduction to Nonlinear Image Processing
Introduction to Nonlinear Image Processing 1 IPAM Summer School on Computer Vision July 22, 2013 Iasonas Kokkinos Center for Visual Computing Ecole Centrale Paris / INRIA Saclay Mean and median 2 Observations
More informationNONLINEAR DIFFUSION PDES
NONLINEAR DIFFUSION PDES Erkut Erdem Hacettepe University March 5 th, 0 CONTENTS Perona-Malik Type Nonlinear Diffusion Edge Enhancing Diffusion 5 References 7 PERONA-MALIK TYPE NONLINEAR DIFFUSION The
More informationError Analysis for H 1 Based Wavelet Interpolations
Error Analysis for H 1 Based Wavelet Interpolations Tony F. Chan Hao-Min Zhou Tie Zhou Abstract We rigorously study the error bound for the H 1 wavelet interpolation problem, which aims to recover missing
More informationA Dual Formulation of the TV-Stokes Algorithm for Image Denoising
A Dual Formulation of the TV-Stokes Algorithm for Image Denoising Christoffer A. Elo, Alexander Malyshev, and Talal Rahman Department of Mathematics, University of Bergen, Johannes Bruns gate 12, 5007
More informationEfficient Beltrami Filtering of Color Images via Vector Extrapolation
Efficient Beltrami Filtering of Color Images via Vector Extrapolation Lorina Dascal, Guy Rosman, and Ron Kimmel Computer Science Department, Technion, Institute of Technology, Haifa 32000, Israel Abstract.
More informationFast Multilevel Algorithm for a Minimization Problem in Impulse Noise Removal
Fast Multilevel Algorithm for a Minimization Problem in Impulse Noise Removal Raymond H. Chan and Ke Chen September 25, 7 Abstract An effective 2-phase method for removing impulse noise was recently proposed.
More informationAnalysis of Numerical Methods for Level Set Based Image Segmentation
Analysis of Numerical Methods for Level Set Based Image Segmentation Björn Scheuermann and Bodo Rosenhahn Institut für Informationsverarbeitung Leibnitz Universität Hannover {scheuerm,rosenhahn}@tnt.uni-hannover.de
More informationNonlinear diffusion filtering on extended neighborhood
Applied Numerical Mathematics 5 005) 1 11 www.elsevier.com/locate/apnum Nonlinear diffusion filtering on extended neighborhood Danny Barash Genome Diversity Center, Institute of Evolution, University of
More informationInverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 2007 Technische Universiteit Eindh ove n University of Technology
Inverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 27 Introduction Fredholm first kind integral equation of convolution type in one space dimension: g(x) = 1 k(x x )f(x
More informationDesign of Image Adaptive Wavelets for Denoising Applications
Design of Image Adaptive Wavelets for Denoising Applications Sanjeev Pragada and Jayanthi Sivaswamy Center for Visual Information Technology International Institute of Information Technology - Hyderabad,
More informationarxiv: v1 [math.na] 3 Jan 2019
arxiv manuscript No. (will be inserted by the editor) A Finite Element Nonoverlapping Domain Decomposition Method with Lagrange Multipliers for the Dual Total Variation Minimizations Chang-Ock Lee Jongho
More informationOn the convergence properties of the modified Polak Ribiére Polyak method with the standard Armijo line search
ANZIAM J. 55 (E) pp.e79 E89, 2014 E79 On the convergence properties of the modified Polak Ribiére Polyak method with the standard Armijo line search Lijun Li 1 Weijun Zhou 2 (Received 21 May 2013; revised
More informationLinear Diffusion and Image Processing. Outline
Outline Linear Diffusion and Image Processing Fourier Transform Convolution Image Restoration: Linear Filtering Diffusion Processes for Noise Filtering linear scale space theory Gauss-Laplace pyramid for
More informationMinimizing Isotropic Total Variation without Subiterations
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Minimizing Isotropic Total Variation without Subiterations Kamilov, U. S. TR206-09 August 206 Abstract Total variation (TV) is one of the most
More informationA Localized Linearized ROF Model for Surface Denoising
1 2 3 4 A Localized Linearized ROF Model for Surface Denoising Shingyu Leung August 7, 2008 5 Abstract 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 1 Introduction CT/MRI scan becomes a very
More informationDETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION
Journal of Fractional Calculus and Applications, Vol. 6(1) Jan. 2015, pp. 83-90. ISSN: 2090-5858. http://fcag-egypt.com/journals/jfca/ DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL
More informationA FAMILY OF NONLINEAR DIFFUSIONS CONNECTING PERONA-MALIK TO STANDARD DIFFUSION. Patrick Guidotti
DISCRETE AND CONTINUOUS DYNAMICAL SYSTEMS SERIES S Volume 5, Number 3, June 202 doi:0.3934/dcdss.202.5.xx pp. X XX A FAMILY OF NONLINEAR DIFFUSIONS CONNECTING PERONA-MALIK TO STANDARD DIFFUSION Patrick
More informationUsing Entropy and 2-D Correlation Coefficient as Measuring Indices for Impulsive Noise Reduction Techniques
Using Entropy and 2-D Correlation Coefficient as Measuring Indices for Impulsive Noise Reduction Techniques Zayed M. Ramadan Department of Electronics and Communications Engineering, Faculty of Engineering,
More informationA Denoising-Decomposition Model Combining TV Minimisation and Fractional Derivatives
East Asian Journal on Applied Mathematics Vol. 8, No. 3, pp. 447-462 doi: 10.4208/eajam.130917.150218 August 2018 A Denoising-Decomposition Model Combining TV Minimisation and Fractional Derivatives Donghuan
More informationRestoration of images with rotated shapes
Numer Algor (008) 48:49 66 DOI 0.007/s075-008-98-y ORIGINAL PAPER Restoration of images with rotated shapes S. Setzer G. Steidl T. Teuber Received: 6 September 007 / Accepted: 3 January 008 / Published
More informationVariational Methods in Signal and Image Processing
Variational Methods in Signal and Image Processing XU WANG Texas A&M University Dept. of Electrical & Computer Eng. College Station, Texas United States xu.wang@tamu.edu ERCHIN SERPEDIN Texas A&M University
More informationImage processing and Computer Vision
1 / 1 Image processing and Computer Vision Continuous Optimization and applications to image processing Martin de La Gorce martin.de-la-gorce@enpc.fr February 2015 Optimization 2 / 1 We have a function
More informationNoisy Word Recognition Using Denoising and Moment Matrix Discriminants
Noisy Word Recognition Using Denoising and Moment Matrix Discriminants Mila Nikolova Département TSI ENST, rue Barrault, 753 Paris Cedex 13, France, nikolova@tsi.enst.fr Alfred Hero Dept. of EECS, Univ.
More informationA Total Fractional-Order Variation Model for Image Restoration with Non-homogeneous Boundary Conditions and its Numerical Solution
A Total Fractional-Order Variation Model for Image Restoration with Non-homogeneous Boundary Conditions and its Numerical Solution Jianping Zhang and Ke Chen arxiv:509.04237v [cs.cv] 6 Sep 205 Abstract
More informationParameter Identification in Partial Differential Equations
Parameter Identification in Partial Differential Equations Differentiation of data Not strictly a parameter identification problem, but good motivation. Appears often as a subproblem. Given noisy observation
More informationANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE INTO CARTOON PLUS TEXTURE. C.M. Elliott and S.A.
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 6, Number 4, December 27 pp. 917 936 ANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE
More informationUPRE Method for Total Variation Parameter Selection
UPRE Method for Total Variation Parameter Selection Youzuo Lin School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ 85287 USA. Brendt Wohlberg 1, T-5, Los Alamos National
More informationAdaptive Primal Dual Optimization for Image Processing and Learning
Adaptive Primal Dual Optimization for Image Processing and Learning Tom Goldstein Rice University tag7@rice.edu Ernie Esser University of British Columbia eesser@eos.ubc.ca Richard Baraniuk Rice University
More informationDenoising of NIRS Measured Biomedical Signals
Denoising of NIRS Measured Biomedical Signals Y V Rami reddy 1, Dr.D.VishnuVardhan 2 1 M. Tech, Dept of ECE, JNTUA College of Engineering, Pulivendula, A.P, India 2 Assistant Professor, Dept of ECE, JNTUA
More informationHigh Order Numerical Methods for the Riesz Derivatives and the Space Riesz Fractional Differential Equation
International Symposium on Fractional PDEs: Theory, Numerics and Applications June 3-5, 013, Salve Regina University High Order Numerical Methods for the Riesz Derivatives and the Space Riesz Fractional
More informationResearch Article Denoising Algorithm Based on Generalized Fractional Integral Operator with Two Parameters
Hindawi Publishing Corporation Discrete Dynamics in Nature and Society Volume 212, Article ID 529849, 14 pages doi:11155/212/529849 Research Article Denoising Algorithm Based on Generalized Fractional
More informationNon-negative Quadratic Programming Total Variation Regularization for Poisson Vector-Valued Image Restoration
University of New Mexico UNM Digital Repository Electrical & Computer Engineering Technical Reports Engineering Publications 5-10-2011 Non-negative Quadratic Programming Total Variation Regularization
More informationResearch Article Nonlocal Mean Image Denoising Using Anisotropic Structure Tensor
Advances in Optical Technologies Volume 2013, Article ID 794728, 6 pages http://dx.doi.org/10.1155/2013/794728 Research Article Nonlocal Mean Image Denoising Using Anisotropic Structure Tensor Xi Wu, 1,
More informationImage Deblurring in the Presence of Impulsive Noise
Image Deblurring in the Presence of Impulsive Noise Leah Bar Nir Sochen Nahum Kiryati School of Electrical Engineering Dept. of Applied Mathematics Tel Aviv University Tel Aviv, 69978, Israel April 18,
More informationMULTI-SCALE IMAGE DENOISING BASED ON GOODNESS OF FIT (GOF) TESTS
MULTI-SCALE IMAGE DENOISING BASED ON GOODNESS OF FIT (GOF) TESTS Naveed ur Rehman 1, Khuram Naveed 1, Shoaib Ehsan 2, Klaus McDonald-Maier 2 1 Department of Electrical Engineering, COMSATS Institute of
More informationSOS Boosting of Image Denoising Algorithms
SOS Boosting of Image Denoising Algorithms Yaniv Romano and Michael Elad The Technion Israel Institute of technology Haifa 32000, Israel The research leading to these results has received funding from
More informationA Tensor Variational Formulation of Gradient Energy Total Variation
A Tensor Variational Formulation of Gradient Energy Total Variation Freddie Åström, George Baravdish and Michael Felsberg Linköping University Post Print N.B.: When citing this work, cite the original
More informationA BAYESIAN APPROACH TO NONLINEAR DIFFUSION BASED ON A LAPLACIAN PRIOR FOR IDEAL IMAGE GRADIENT
A BAYESIAN APPROACH TO NONLINEAR DIFFUSION BASED ON A LAPLACIAN PRIOR FOR IDEAL IMAGE GRADIENT Aleksandra Pižurica, Iris Vanhamel, Hichem Sahli, Wilfried Philips and Antonis Katartzis Ghent University,
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationScale Space Analysis by Stabilized Inverse Diffusion Equations
Scale Space Analysis by Stabilized Inverse Diffusion Equations Ilya Pollak, Alan S. Willsky and Hamid Krim Laboratory for Information and Decision Systems, Massachusetts Institute of Technology, 77 Massachusetts
More informationPatch similarity under non Gaussian noise
The 18th IEEE International Conference on Image Processing Brussels, Belgium, September 11 14, 011 Patch similarity under non Gaussian noise Charles Deledalle 1, Florence Tupin 1, Loïc Denis 1 Institut
More informationMotion Estimation (I) Ce Liu Microsoft Research New England
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationA discretized Newton flow for time varying linear inverse problems
A discretized Newton flow for time varying linear inverse problems Martin Kleinsteuber and Simon Hawe Department of Electrical Engineering and Information Technology, Technische Universität München Arcisstrasse
More informationA Partial Differential Equation Approach to Image Zoom
A Partial Differential Equation Approach to Image Zoom Abdelmounim Belahmidi and Frédéric Guichard January 2004 Abstract We propose a new model for zooming digital image. This model, driven by a partial
More informationImage enhancement by non-local reverse heat equation
Image enhancement by non-local reverse heat equation Antoni Buades 1, Bartomeu Coll 2, and Jean-Michel Morel 3 1 Dpt Matematiques i Informatica Univ. Illes Balears Ctra Valldemossa km 7.5, 07122 Palma
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationSPARSITY-ASSISTED SIGNAL SMOOTHING (REVISITED) Ivan Selesnick
SPARSITY-ASSISTED SIGNAL SMOOTHING (REVISITED) Ivan Selesnick Electrical and Computer Engineering Tandon School of Engineering, New York University Brooklyn, New York, USA ABSTRACT This paper proposes
More informationOpen Access Measuring Method for Diameter of Bearings Based on the Edge Detection Using Zernike Moment
Send Orders for Reprints to reprints@benthamscience.ae 114 The Open Automation and Control Systems Journal, 2015, 7, 114-119 Open Access Measuring Method for Diameter of Bearings Based on the Edge Detection
More informationUnravel Faults on Seismic Migration Images Using Structure-Oriented, Fault-Preserving and Nonlinear Anisotropic Diffusion Filtering
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2019 SGP-TR-214 Unravel Faults on Seismic Migration Images Using Structure-Oriented,
More informationA Linearly Convergent First-order Algorithm for Total Variation Minimization in Image Processing
A Linearly Convergent First-order Algorithm for Total Variation Minimization in Image Processing Cong D. Dang Kaiyu Dai Guanghui Lan October 9, 0 Abstract We introduce a new formulation for total variation
More informationTruncation Strategy of Tensor Compressive Sensing for Noisy Video Sequences
Journal of Information Hiding and Multimedia Signal Processing c 2016 ISSN 207-4212 Ubiquitous International Volume 7, Number 5, September 2016 Truncation Strategy of Tensor Compressive Sensing for Noisy
More informationInvestigating the Influence of Box-Constraints on the Solution of a Total Variation Model via an Efficient Primal-Dual Method
Article Investigating the Influence of Box-Constraints on the Solution of a Total Variation Model via an Efficient Primal-Dual Method Andreas Langer Department of Mathematics, University of Stuttgart,
More informationNonlinear regularizations of TV based PDEs for image processing
Nonlinear regularizations of TV based PDEs for image processing Andrea Bertozzi, John Greer, Stanley Osher, and Kevin Vixie ABSTRACT. We consider both second order and fourth order TV-based PDEs for image
More informationCache Oblivious Stencil Computations
Cache Oblivious Stencil Computations S. HUNOLD J. L. TRÄFF F. VERSACI Lectures on High Performance Computing 13 April 2015 F. Versaci (TU Wien) Cache Oblivious Stencil Computations 13 April 2015 1 / 19
More informationA fast nonstationary preconditioning strategy for ill-posed problems, with application to image deblurring
A fast nonstationary preconditioning strategy for ill-posed problems, with application to image deblurring Marco Donatelli Dept. of Science and High Tecnology U. Insubria (Italy) Joint work with M. Hanke
More informationITK Filters. Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms
ITK Filters Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms ITCS 6010:Biomedical Imaging and Visualization 1 ITK Filters:
More informationNonlinear Diffusion. 1 Introduction: Motivation for non-standard diffusion
Nonlinear Diffusion These notes summarize the way I present this material, for my benefit. But everything in here is said in more detail, and better, in Weickert s paper. 1 Introduction: Motivation for
More informationA Riemannian Framework for Denoising Diffusion Tensor Images
A Riemannian Framework for Denoising Diffusion Tensor Images Manasi Datar No Institute Given Abstract. Diffusion Tensor Imaging (DTI) is a relatively new imaging modality that has been extensively used
More informationIntroduction to Sparsity in Signal Processing
1 Introduction to Sparsity in Signal Processing Ivan Selesnick Polytechnic Institute of New York University Brooklyn, New York selesi@poly.edu 212 2 Under-determined linear equations Consider a system
More informationImage Noise: Detection, Measurement and Removal Techniques. Zhifei Zhang
Image Noise: Detection, Measurement and Removal Techniques Zhifei Zhang Outline Noise measurement Filter-based Block-based Wavelet-based Noise removal Spatial domain Transform domain Non-local methods
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 402 Iterative TV-regularization of Grey-scale Images Martin Fuchs and Joachim Weickert Saarbrücken
More informationNotes on Regularization and Robust Estimation Psych 267/CS 348D/EE 365 Prof. David J. Heeger September 15, 1998
Notes on Regularization and Robust Estimation Psych 67/CS 348D/EE 365 Prof. David J. Heeger September 5, 998 Regularization. Regularization is a class of techniques that have been widely used to solve
More informationImage Deblurring in the Presence of Impulsive Noise
Image Deblurring in the Presence of Impulsive Noise Leah Bar Nir Sochen Nahum Kiryati School of Electrical Engineering Dept. of Applied Mathematics Tel Aviv University Tel Aviv, 69978, Israel Abstract
More informationarxiv: v1 [cs.cv] 18 Sep 2017
Noname manuscript No. will be inserted by the editor) Variational Methods for Normal Integration Yvain Quéau Jean-Denis Durou Jean-François Aujol arxiv:709.05965v [cs.cv] 8 Sep 07 the date of receipt and
More informationMichal Kuneš
A (Zernike) moment-based nonlocal-means algorithm for image denoising Michal Kuneš xkunes@utia.cas.cz ZOI UTIA, ASCR, Friday seminar 13.03.015 Introduction - Uses nonlocal (NL) means filter - Introduce
More informationBayesian-based Wavelet Shrinkage for SAR Image Despeckling Using Cycle Spinning *
Jun. 006 Journal of Electronic Science and Technology of China Vol.4 No. ayesian-based Wavelet Shrinkage for SAR Image Despeckling Using Cycle Spinning * ZHANG De-xiang,, GAO Qing-wei,, CHEN Jun-ning.
More informationImage Cartoon-Texture Decomposition and Feature Selection using the Total Variation Regularized L 1 Functional
Image Cartoon-Texture Decomposition and Feature Selection using the Total Variation Regularized L 1 Functional Wotao Yin 1, Donald Goldfarb 1, and Stanley Osher 2 1 Department of Industrial Engineering
More informationA memory gradient algorithm for l 2 -l 0 regularization with applications to image restoration
A memory gradient algorithm for l 2 -l 0 regularization with applications to image restoration E. Chouzenoux, A. Jezierska, J.-C. Pesquet and H. Talbot Université Paris-Est Lab. d Informatique Gaspard
More information