Image Processing by Digital Filter Using Matlab

Download as pdf or txt
Download as pdf or txt
You are on page 1of 7
 
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: editor@ijaiem.org Volume 3, Issue 9, September 2014 ISSN 2319 – 4847 Volume 3, Issue 9, September 2014
Page 126
 
ABSTRACT
 
This paper highlights different image processing stages, components of a typical digital image processing system and strategies  concerning digital filter design for the applications of image smoothing. The aim of this paper is to study the implementation of  two dimensional FIR digital filter using weighted least square method for the application of image demising. The Oxford  dictionary defines the word image as the optical appearance of something produced in a mirror or through a lens. Image may be  other types of radiant energy and devices . However, optical images are most common and are most important. Amount of light (radiant) energy received at a point of a scene by an observer or by an image sensor varies with direction and distance of that  point. The term scene stands for a landscape or a view as seen by a spectator. The scene, indoor, we see around us is, in general,  three dimensional, i.e. each point in the scene may have different depth from the observer. However, there are scenes which can  be considered as two dimensional, i.e. points comprising the scene have equal depth.
KEYWORDS-COMPONENT:
 Matlab software, SIM 300 GSM/GPRS Modem Some examples of two-dimensional scenes are:
 
flat terrain seen from satellite situated at a high altitude.
 
the shadow of objects created by X-ray,
 
flat polished surface of object or thin slices seen through microscopic, and
 
documents printed on paper. In fact, radiant energy is recorded at corresponding points on a plane to form an image . Hence, the brightness and colour recorded in an image may be represented as a function of several variables. The simplest kind of intensity (optical) image we can think of is a black-and-white (B/W) image. These image are most common and are represented by a function of two variables
 y xg
,
 where
 y xg
,
 is the grayness or brightness or intensity of the image at the spatial coordinate
 y x
,
. A multi-spectral image is a vector-valued function having number of components equal to that of spectral bands and is represented by
 y xg y xg y xg
n
,,,,,,
21
 
 at each
 y x
,
. The colour image is a special case of multi-spectral image: and if we consider intensity measured in three wavelengths corresponding to red, green and blue, then
 y xg
,
 =
 y xg y xg y xg
 BG R
,,,,,,
 
. A temporal argument is added to represent time-varying image
 y xg
,,
 (also called image sequences).
1.
 
INTRODUCTION
The process of receiving and analyzing visual information by the human species is referred to as sight, perception or understanding. Similarly, the process of receiving and analyzing visual information by digital computer is called digital image processing and scene analysis. The term ‘image’ rather than’ picture is used here, because the computer stores and  processes numerical image of a scene. Although ‘pattern recognition’ and ‘image processing’ have a lot in common, yet they were developed as separate disciplines . Two broad classes of techniques, viz., processing and analysis have evolved in the field of digital image processing and analysis. Processing of an image includes improvement in its appearance and efficient representation. So the field consists of not only feature extraction, analysis and recognition of image, but also coding, filtering, enhancement and restoration.
IMAGE PROCESSING BY DIGITAL FILTER USING MATLAB
Netranand Pathak
1
, Ketandeep Jamwal
2
 , Mukund Madhav
3
 , Aditya Madhav
4
 
1
M-Tech, Electrical Engineering (Power system) GALGOTIA UNIVERSITY-NOIDA
1-3
 B-Tech, Electronics & Communication Engineering, BHSBIET-Lehragaga
 
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: editor@ijaiem.org Volume 3, Issue 9, September 2014 ISSN 2319 – 4847 Volume 3, Issue 9, September 2014
Page 127
 
Enhancement Restoration Segmentation Feature extraction Primitive and relation extraction Statistical/str uctural analysis
Description and interpretation
Image sensing
Scene
Digitization Compression Mass storage
Transmission Improved image Fig.1
 Different stages of image processing and analysis scheme (Information flow along the outgoing link from each block)
Digital image
2.
 
METHODOLOGY
The entire process of image processing and analysis starting from the receiving of visual information to the giving out of description of the scene, may be divided into three major stages. Which are also considered as major sub-areas
(a)
 
Discretization and representation:
 This is concerned with converting visual information into a discrete form and making it compatible for computer processing, approximating visual information to save storage space as well as time requirement in subsequent processing.
(b)
 
Processing:
It is concerned with improving image quality by filtering etc, compressing data to save storage and channel capacity during transmission.
(c) Analysis:
 It is mainly related to extracting the image features, quantifying shapes, registration and recognition. In the initial stage, the input is a scene (visual information), and the output is corresponding digital image. In the secondary stage, both the input and the output are images where the output is an improved version of the input. And, in the final stage, the input is still an image but the output is a description of the contents of that image.
 
3.
 
RESULTS
Given the impulse response of a digital filter, its frequency response is used to study its behavior to different frequency components. In case of symmetric 2-D FIR digital filter, the magnitude response is same along both the frequency axes, i.e.,
1
 
 and 
2
 
. Therefore, maximum magnitude of ripples along the two-frequency axes is same in case of 2-D symmetric FIR digital filter. In the design of FIR digital filters using WLS techniques, the update of weights is important  because it affects the convergence of the procedure. There are no analytical methods to ensure that the iterative update schemes in weights are convergence as pointed out by many workers. The convergence problem of course, is still inevitable in the proposed method. The constant parameter
 determines the update of weights whereas
 and
 
 
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: editor@ijaiem.org Volume 3, Issue 9, September 2014 ISSN 2319 – 4847 Volume 3, Issue 9, September 2014
Page 128
 
determine the update of G. Eventually, all the parameters affect the convergent properties of the algorithm such as the convergent speed and the achievable peak ripple magnitude (prm). There are no analytical schemes for obtaining these  parameters, and the relationship between these parameters are still unknown. However, from the experience of numerous design examples in the literature, we can make a reasonable choice of these parameters. In this low-pass filter with size 5
5 L = 40,
 p
 
 p = 0.4
 and
s
 
 = 0.6
 is used for the illustration of the choice.
1.
 
(a) The constant
 
 : It is an acceleration factor in the course of optimization when
 
 > 1. As seen from Fig. 4 when
 
 = 2 or 2.4 the procedure convergence faster significantly than the other cases with a small increase of minimum  prm. Therefore, the choice of
 
 = 2 or 2.4 is better under the consideration of convergent speed and convergent error. However, in this work, we choose
 
 = 2.
2.
 
(b) The constant
 
:
 
 can also be viewed as an acceleration factor for the algorithm. The larger
 
 value may result in faster convergence. However, it may cause divergence. For,
 
 > 0.5 it has been found that the proposed algorithm fails to converge for several cases. For
 
 < 0.4 the algorithm seems converges eventually, whereas the convergent speed is slow or the final prm is relatively large. However, experiments show that
 
 = 0.45 is a better choice for the tradeoff between convergent speed and convergent error, as illustrated in Fig. 4. Based on these observations, in this design example, the algorithm diverges when
 
 > 0.62 except
 
 = 0.68. The number of iterations corresponding to the
 
 values, which yield divergence, are meaningless.
3.
 
c) The constant
 
:It has been observed that in few cases, the prm of the (k+1)th iteration is large than that if the k 
th
 iteration for some K. This may cause slow convergence. To attack the problem, we try to inhibit the increase of prm by decreasing the update amount of G, when the situation occurs. In this work,
 
 = 1 is chosen initially since it gives the maximum convergent speed. Once pr m
1
> pr m
 has been detected,
 
 is changed to a large value. In other words, two values
 
 = 1 and
 
=
 
*
 are used in turn indicate that
 
7050*
 
 is a good consideration of convergent speed and convergent error. Therefore, in this work, we choose
 
*
= 50.
 
4.F
REQUENCY
R
ESPONSE
E
VALUATION OF THE
I
MPLEMENTED
F
ILTER
The symmetric 2-D FIR filter is designed using two dimensional weighted least squares method for different filter specifications, i.e., pass-band frequency, stop-band frequency, tolerance and order of the filter. The frequency response comparison of these filters with different specifications is illustrated in Figs.3-7. The frequency responses so obtained are normalized frequency responses for the purposes of comparison. Based on literature and our experience of working with the present problem as well as depending upon the requirements of the present application for which the filter is to be designed, the tolerance limit selected in the present case is 0.01. If the value for the tolerance limit selected is smaller than the selected choice, the filter response improves but at the same time, computational time becomes much larger and vice-versa. So the selection of the tolerance limit, is in fact, a trade off between the computational time and filter response. From the results so obtained, first of all, it has been observed that with the modified algorithm, the frequency response of the implemented filter is always circular
Fig.2
 Performance analysis of the implemented two dimensional digital FIR filter using weighted least squares method (a) Frequency response of a circularly symmetric low pass filter (Order, N = 5, Tolerance = 0.01, Pass-band radius = 0.4
 
 
International Journal of Application or Innovation in Engineering& Management (IJAIEM)
Web Site: www.ijaiem.org Email: editor@ijaiem.org Volume 3, Issue 9, September 2014 ISSN 2319 – 4847 Volume 3, Issue 9, September 2014
Page 129
 
and Stop-band radius = 0.6
); (b) Number of iterations required for convergence versus
 plot when
 = 2;
*
 = 50 (c)  Number of iterations required for convergence versus
 plot when
 = 2;
 = 0.45; symmetric as illustrated in the Figs 3-7. The transition band obtained is quite sharp in almost all the Figs.3-7 with the implemented technique as compared to that of other design methods reported in the literature. Also, as shown in the Fig.3, the frequency response of low pass symmetric FIR filter is obtained for the specifications: (Order of the filter, N = 5; pass-band cut-off frequency = 0.11
; stop-band cut-off frequency = 0.99
; and Tolerance = 0.01). It has been observed from the Fig.4 that the frequency response is circular symmetric but the magnitude of ripples in the pass-band as well as stop-band is extremely high and it is much beyond the tolerance limits. Hence the filter with these specifications is not suitable for the present problem. It has been further observed that with the increase of pass-band cut-off frequency from 0.11
 to 0.16
, the magnitude of ripples decrease sharply and consequently the filter response improves greatly as verified from Figs.4-7. In this manner the final performance parameters selected for the design of filter are: (Order of the filter, N = 5; pass-band cut-off frequency = 0.16
; stop-band cut-off frequency = 0.99
; and Tolerance = 0.01). On the other hand, the ripple magnitude can also be reduced by increasing the order of the filter but at the same time computational complexity and hence convergence time increases, hence the use of higher order filter beyond an order 5 is also not an optimal solution. Instead we made the use of varying the pass-band cut-off frequency while keeping the stop-band cutoff frequency on the maximum side, i.e., 0.99
 for a fifth order filter in the present problem. Further, the transition band is kept purposely larger keeping stop-band cut-off frequency at 0.99
, so that the design time is reduced and also to get a wider range for the selection of pass-band frequency so the image information is not lost substantially and at the same time all types of images can be handled with a single filter having same performance parameters. It has also been observed that the sharp ripples occur in the transition band in all the Figs.3-7. However, at the same time with use of the selected parameters, the ripples in the Fig 3 are found to be much smoother and are equally distributed in the entire frequency domain.
Fig.3
 Frequency response of a circularly symmetric low-pass filter using weighted least square method for Order of the filter, N = 5, PB cut-off frequency = 0.16
, SB cut-off frequency = 0.99
, and Tolerance = 0.01

Reward Your Curiosity

Everything you want to read.
Anytime. Anywhere. Any device.
No Commitment. Cancel anytime.
576648e32a3d8b82ca71961b7a986505