American Journal of Signal Processing

p-ISSN: 2165-9354    e-ISSN: 2165-9362

2017;  7(2): 39-43



The Application of Numerical Approximation Methods upon Digital Images

Ali Jassim Mohamed Ali

Department of Physics, College of Science, Mustansiriyah University, Baghdad, Iraq

Correspondence to: Ali Jassim Mohamed Ali, Department of Physics, College of Science, Mustansiriyah University, Baghdad, Iraq.


Copyright © 2017 Scientific & Academic Publishing. All Rights Reserved.

This work is licensed under the Creative Commons Attribution International License (CC BY).


Numerous problems arise in diverse areas of science and engineering, as well as from the physical, computer, biological, economic, and even social sciences. The need for choosing such an application is more clearly and concisely demonstrate how shall the numerical technique be applied in such real-life situations. Numerical techniques, explore the required theory to get numerical solutions to the mathematical problem, especially when computer deals with such huge and complicated calculations. By assumption that the root is exist in image's gray level interval, an attempt to explorer the application of some numerical algorithms is tested here by using some numerical analysis algorithms. An estimation to the error bounds for each approximation process had been calculated. Bisection, Newton Raphson, Secant and False Position methods are some of these methods which have been used here upon some digital images. Among the various used approximation methods and according to subjective and quantitative evaluation results, one can be noted that the Bisection method is the best approximation technique. Besides, image's histogram plays the main role in sequencing such methods of approximations.

Keywords: Bisection method, Newton Raphson, Secant method, False Position method, Absolute Error, Relative Error

Cite this paper: Ali Jassim Mohamed Ali, The Application of Numerical Approximation Methods upon Digital Images, American Journal of Signal Processing, Vol. 7 No. 2, 2017, pp. 39-43. doi: 10.5923/j.ajsp.20170702.01.

1. Introduction

One of the primary objectives of numerical analysis is to decide which algorithms can be implemented to solve mathematical problems on computer [1]. These algorithms accurately approximate the solution of problems which the latter can't be solved exactly and after that a typical techniques for estimating error bounds for the approximation process can be handled [2]. Many researches focused their attention toward using such methods to solve their problems. C. Thinzar and N. Aye made use of the digital image correlation that is based on the sub-pixel accuracy with the aid of Newton Raphson method to detect storm direction for managing and to prevent the possible land's attacked region [3]. D. Biolè, etc. used Secant method in addition to other techniques to measure apparent contact angles of non-axisymmetric drops [4]. Md. G. Moazzam, etc. used a robust method for solving transcendental equations. The presented method by them gave better results rather than the well-known False Position method and the Bisection method for solving such equations [5]. W. Chu, etc. demonstrated a new registration algorithm based on Newton Raphson iteration process to align images with rigid body transformation [6]. So, an attempt to detect error propagation which can be occurred by applying such numerical methods to digital image had been utilized here. These methods are:

1.1. Bisection Method

It is one of the most basic problems of numerical approximation. It is also known as the root-finding problem or Binary-search method. It is a process of dividing a set continually in half to search for the solution to a certain problem [2]. It states that if f(x) is a continuous function defined on the interval a<x<b, with f(a) and f(b) of opposite sign, then there is at least one root in the interval between a and b [5]. By setting a1=a and b1=b and take the midpoint of [a,b] such that [7, 8]
So the next steps must be taken in the considerations:
• If f(p1)=0; a root was found at p1.
• If f(a1).f(b1)<0, then there is a root between a1 and p1.
• If f(a1).f(b1)>0, then there is a root between p1 and b1. And so on.

1.2. Newton Raphson Method

It is one of the well-known numerical approaches for solving a root-finding problem [2]. This technique supposes that po be an approximation to the root (i.e. iteration guess) and belongs to the interval [a,b] in such that f '(po)≠0 and |p-po| is so small and hence [7, 9]

1.3. Secant Method

This method uses a secant line joining two points that cut curve's function and can be presented in the following expression [2, 8]:

1.4. False Position Method

In this technique, one uses results that are known to be false to converge to the true root. This method chooses an initial approximations po and p1 such that f(po).f(p1)<0. The new approximation value will be then obtained in the way as described in the secant method mentioned above. After that and in order to decide which secant line to be used, the product of f(p2) and f(p1) should be taken and one of the following considerations shall be verified [7]:
• Choose p3 as a line joining (p1,f(p1)) and (p2,f(p2)) if f(p1).f(p2)<0.
• Choose p3 as a line joining (po,f(po)) and (p2,f(p2)) if f(p1).f(p2)>0.
In order to estimate the approximation errors, two methods had been used here; these are the absolute and relative errors as shown in the following expressions [2]:
Where p* is an approximation value for p.

2. Result and Discussion

For simplicity, an assumption had been made here that the root in the interval (i.e. image's data) is unique and the stopping criteria shall be as follows:
where ε is the tolerance. Figure (1) shows the used images with their fitting histograms' curves. The fitting process had been executed by using Table Curve 2D V 5.01 program.
Figure 1. (a) Vegetable image, (b) fitting histogram for (a), (c) Eye Image, (d) fitting histogram for (c), (e) Balloon image, (f) fitting histogram for (e)
In general and as a result of application, numerical approximation methods distorted the images in alternatives ratios depending on image's histogram. An image with sharp chromatic gradient is better than other used images which have a wide range of color gradients, this can be seen through noticing Figs. (2-4). In addition to that, one can noticed that there is no clear sequence for algorithms' effects upon single image. In most cases, one can find that the best numerical methods' performance is the Bisection method with lowest errors presented in Table (1).
Table 1. Approximation Errors for all used images
Figure 2. Histograms for vegetable image before and after numerical methods' applications
Figure 3. Histograms for eye image before and after numerical methods' applications
Figure 4. Histograms for balloon image before and after numerical methods' applications
A distorted case occurred in eye image in the application of the Newton Raphson method with a shift toward lower values of the image's gray level while it isn't the case for other used images; the distortion case resulted near the ends of images' gray levels for all used algorithms in general, this can be recognized through noticing approximation errors explained in the above table. Image's histogram restricts the sequence of algorithms in their applications; the symmetric behavior for both Secant and False Position methods resulted in the case of vegetable and balloon images in contraction to that occurred to the case in eye image which the latter witnessed a strong distortion state in the case of Secant method technique.

3. Conclusions

By assumption that the root is exist in image's gray level interval, an attempt to explorer the application of some numerical algorithms is tested here by using some numerical analysis algorithms. As a consequence result, one can noticed that the used numerical algorithms distorted the image after all with different distortion ratios depend upon image's histogram. In addition to that, the sequence for techniques' effects is not clear upon such used images but in general one can be concluded that the Bisection method can be recognized as the best one with least absolute and relative errors.


[1]  Z. Xu. (2010) Numerical Analysis, ACMS 40390. [Online]. Available:
[2]  R. L. Burden and J. D. Faires, Numerical Analysis, Ninth Edition, Brooks/ Cole, Cengage Learning, 2011.
[3]  C. Thinzar, and N. Aye, "Detection the storm movement by sub pixel registration approach of Newton Raphson method" International Journal of e-Education, e-Business, e-Management and e-Learning, Vol. 4, No. 1, 2014.
[4]  D. Biolè, M. Wang, and V. Bertola, "Assessment of direct image processing methods to measure the apparent contact angle of liquid drops" Experimantal Thermal and Fluid Science, Vol. 76, 2016. Available:
[5]  Md. G. Moazzam, A. Chakraborty, and Md. A. Bhuiyan, "A roubust method for solving transcendental equations" IJCSI International Journal of Computer Science Issues, Vol. 9, Issue 6, No. 2, 2012.
[6]  W. Chu, L. Ma, J. Song, and T. Vorburger, "An iterative image registration algorithm by optimizing similarity measurement" Journal of Research of the National Institute of Standards and Technology, Vol. 115, No. 1, 2010.
[7]  J. Kiusalaas, Numerical Methods in Engineering with Matlab, Second Edition, Cambridge University Press, 2010.
[8]  S. C. Chapra, Applied Numerical Methods with Matlab ® for Engineers and Scientists, Third Edition, McGraw-Hill Companies, Inc., 2012.
[9]  S. T. Karris, Analysis Using Matlab ® and Excel ®, Third Edition, Orchard Publications, 2007.