Auto-Tuning of PID Controller: Bachelors of Engineering Mr. Rahul S. Roy Roll No. 12IN2015
Auto-Tuning of PID Controller: Bachelors of Engineering Mr. Rahul S. Roy Roll No. 12IN2015
Bachelors of Engineering
by
Supervisor
2014
Declaration
I declare that this written submission represents my ideas in my own words and where
others ideas or words have been included, I have adequately cited and referenced the original
sources. I also declare that I have adhered to all principles of academic honesty and integrity
and have not misrepresented or fabricated or falsified any idea/data/fact/source in my submission. I understand that any violation of the above will be cause for disciplinary action by the
Institute and can also evoke penal action from the sources which have thus not been properly
cited or from whom proper permission has not been taken when needed.
..............................
(Signature)
..............................
(Mr. Rahul S. Roy and Roll No. 12IN2015)
Date :
Examiners :
1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Supervisors :
1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Chairman
..............................
Date :
Place :
Acknowledgments
With great pleasure, I avail this opportunity to express my profound gratitude and deep
regards to my project co-ordinator Mr.Shirish S. Kulkarni for their spirited guidance, monitoring and constant encouragement throughout the completion of this seminar report. I have
deep sense of admiration for their innate goodness and inexhaustible enthusiasm, it helped me
to work in right direction to attain desired objective.
I am also thankful to Mr.Sharad P. Jadhav, Head of Department of Instrumentation Engineering, RAIT, Nerul for his generous support, devoting their valuable time and helped me
in all possible ways towards successful completion of this work. I thank all those who have
contributed directly or indirectly to this work.
I take this privilege to express my sincere thanks to Dr.Ramesh Vasappanavara, Principal,
RAIT for their support, encouragement and providing the much necessary facilities. I extend
thanks to my friends who have supported in every stage of these report. I cannot end without
thanking my lovely family for their encouragement.
Date
Signature
Abstract
Image registration, or geometric alignment of two-dimensional or three-dimensional
(3-D) image data, is becoming increasingly important in diagnosis, treatment planning, functional studies, computer-guided therapies and in biomedical research. Image registration is an
fundamental and tough task in image processing technique. Especially, CT and MR imaging
for diagnosis and surgical planning indicates that physicians and surgeons gain important information from these modalities, this dissertation work focuses on the implementation of image
registration technique.
In medical image processing field, some techniques are proposed to find a geometrical
transformation that relates the points of an image to their corresponding points of another
image. This report deals with the use of phase correlation for estimation of mutual translation
between images. Log-Polar Transform (LPT) method will be also proposed, thus improving the
accuracy and precision of feature-based registration which will include mutual rotation, scale
and translation. A key feature of LPT methods is the speed offered by the use of FFT routines.
The proposed method will estimates large motions accurately and robustly without the need of
excessive zero padding and over-sampling, thus without sacrificing part of the computational
efficiency for the frequency domain formulation.
Contents
1 Introduction
1.1
1.2
Proposed Method . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.3
Organization Of Report . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2 Literature Review
3.1
3.2
Principle of PCM . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.3
Strategy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3.4
Applications . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
4 Performance Results
12
4.1
Algorithm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
4.2
5 Conclusion
14
Bibliography
15
ii
List of Figures
1.1
iii
List of Tables
3.1
Calculation Table . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
iv
Chapter 1
Introduction
Medical images are increasingly widely used in healthcare and bio-medical research from
wide range of imaging modalities which are now available, such as CT, MRI, PET, SPECT, and
so on. In some clinical scenarios, the information from several different imaging modalities
should be integrated to deduce useful clinical conclusions. Image registration aligns the images
and so establishes correspondence between different features contained on different imaging
modalities, allows monitoring of subtle changes in size or intensity over time or across a population and establishes correspondence between images and physical space in image guided
interventions [2]. Of the multitude of Image registration similarity measures that have been
proposed over the years, mutual information is currently one of the most intensively researched
measures. This attention is a logical consequence of both the favourable characteristics of
the measure and the good registration results reported. Mutual information is an automatic,
intensity-based metric, which does not require the definition of landmarks or features such as
surfaces and can be applied in retrospect. Furthermore, it is one of the few intensity-based
measures well suited to registration of multimodal images.
Chapter 2
Literature Review
Batler et al. registered 2-D projection images using a curve matching method [2]. In these
method samples are generated to sequences of corresponding points to be registered. Corresponding open curves are manually identified, then are matched by searching for the optimal
fit of the local curvatures in the two curves which leads to registration of two-dimensional
projection radiographies. Nicholas Ayache and Andre P. Gueziec have presented a method of
matching of 3D objects extracted from medical images. Crest lines computed on the object
surfaces correspond to meaningful anatomical features, and that they are stable with respect
to rigid transformations. The extraction of the crest lines is done by computing up to third
order derivatives of the image intensity function with appropriate 3D filtering of the volumetric
images, and by the marching lines algorithm. The recovered lines are then approximated by
spline curves, to compute at each point a number of differential invariants. Matching is finally
performed by a new geometric hashing method. The whole chain is now completely automatic,
and provides extremely robust and accurate results, even in the presence of severe occlusions.
Peng Wen [3] , has suggested a medical image registration method using points, contour and
curves is proposed, which has the accuracy of registration based-on points and the robust of
registration based-on lines (including contour and curves). The operator can exactly extract the
features from the images with semi-automatic extraction method.
Audette and Ferrie [4] have reviewed the surface based registration methods in medical
Chapter 3
Image Registration using Phase
Correlation
Phase correlation measures the local movement between two images from their phases
with sub pixel accuracy. These operation results not only higher registration accuracy but also
lower computational complexity.
B(k,l) = A(k,l)exp-j(ka+lb)
(1)
where (k,l) are the fourier domain coordinates, and a and b are the magnitude of the horizontal and vertical shifts that occur between A and B.
To identify a and b using the PCM approach, one computes a normalized cross power
spectrum between A and B to identify the phase correlation matrix
B(k,l)A(k,l)
Q(k,l)= |B(k,l)A(k,l)
| =exp-j(ka+lb)
(2)
Date
Price
Size
Yesterday
big
Today
small
(1)
The shift in the spatial-domain is thus reflected as a phase change in the frequency domain.
Thus by obtaining complex-valued cross-power spectrum expression
Ck,k+1 (f1 , f2 ) = Sk+1 (f1 , f2 )Sk (f1 , f2 )
Ramrao Adik Institute of Technology
(2)
9
(3)
(4)
3.3 Strategy
The illuminating feature is that the problem of finding the exact lateral shift between two
images is recast as finding the rank-one approximation of the normalized phase correlation
matrix, Q. A close inspection reveals that the noise-free model for Q is in fact a rank-one
matrix. Each element in Q can be separated as;
Q(k,l)= exp(-jka)exp(-jlb)
(3)
(4)
and the phase correlation matrix can be rewritten as Q = qa qbH . A straightforward approach
to finding the dominant rank-one subspace of Q is to use the singular value decomposition
(SVD). The linear phase coefficients can then be identified independently in the left and right
dominant singular vectors. From these, estimates of the vertical and horizontal shift can be
derived, even for noninteger translational motion over a large range. To identify the linear
phase in each of the right and left dominant singular vectors of Q.
Ramrao Adik Institute of Technology
10
3.4 Applications
The advantage of computing frequency domain phase correlation is speed. A spatial domain phase correlation requires a significantly larger number of multiplications due to the
added step of the convolution operation. The convolution, as translated in the frequency domain is a simple point by point multiplication leading to a large reduction in the required number of multiplication[12]. The addition of real-time correlation in medical imaging systems to
allow the examination of 3D and 4D images presents significant benefits to all involved, and
thus can help guide medical professionals in the use and adoption of these new technologies.
Image registration using phase correlation finds its applications [7] in medicine combining
data from different modalities e.g. computed tomography (CT) and magnetic resonance imaging (MRI), to obtain more complete information about the patient, monitoring tumor growth,
treatment verification, comparison of the patients data with anatomical atlases[6], in cartography (map updating), remote sensing (multi-spectral classification), integrating information
into geographic information systems (GIS), and in computer vision. It demonstrates superior
performance in identification of cancerous lesions with the use of registration.
11
Chapter 4
Performance Results
The Fourier transform has several properties that can be exploited for image registration. Translation, rotation, reflection and scale all have their counterpart in Fourier domain[3].
Moreover, using the frequency domain, it is possible to achieve excellent robustness against
correlated and frequency dependent noise. An elegant method used to register two images
which are shifted relative to one another is to use phase correlation [9]. Phase correlation relies
on the translation property of the Fourier transform better known as the shift theorem. The
algorithm in figure 4.1 has been implemented using MATLAB software. Performance analysis
has been done at each stage and the results are evaluated and explained below. We tested the
algorithm on different images.
4.1 Algorithm
Given two images f1 and f2 which differ by translation (dx, dy), i.e.;
f2 (x, y) = f1 (x dx, y dy)
(1)
(2)
The principle used of correlation is for matching. In matching, f1 (x, y) is an image containing objects or regions of interested and let f2 (x, y) be that object or region (template). Then, if
12
1
MN
PM 1 PN 1
m=0
m=0
f1 (m, n) f2 (x + m, y + n)
(3)
2D-FFT of both the images is calculated. This step will give two 2-D arrays of FFT
co-efficient I1 and I2 . Then the magnitude of FFT co-efficient is convolved.
IFFT is calculated to find out the matching point (Tpeak) having maximum magnitude.
The co-ordinates of the above obtained peak (Tpeak) will give the translation between
the images.
Thus by removing the translation from the image, will lead to registered image.
13
Chapter 5
Conclusion
In this report, algorithms to spatially register digital images is being incorporated. The
theoretical motivation for the employment of Image registration using phase correlation is provided. The method is derived by expressing a robust matching surface in terms of correlations.
The speed of the method comes from using fast fourier transforms to compute the correlations.
In this approach the transform coefficients are normalized to unit magnitude prior to computing
correlation in the frequency domain. Thus, the correlation is based only on phase information
and is insensitive to changes in image intensity. Although it has the drawback that all transform
components are weighted equally. In principle one should select the spectral pre-filtering so
as to maximize the expected correlation signal-to-noise ratio given the expected second order
moments of the signal and signal noise. The rotation and scale invariant properties of the logpolar transform, along with fast fourier transform and phase correlation allow us to develop
a robust algorithm that works faithfully under geometric distortions like rotation, scale and
translation. Hence this approach is highly effective in registering medical images.
14
Bibliography
[1] Chia-Hsiang Menq Peng Cheng. Real-time continuous image registration enabling ultraprecise 2-d motion tracking. IEEE TRANSACTIONS ON IMAGE PROCESSING,
5:20812090, 2013.
[2] Chen GT Balter JM, Pelizzari CA. Correlation of projection radiographs in radiotherapy
using open curve segments and points. Med Physics, 19:329334, 2010.
[3] Peng Wen. Medical image registration based-on points, contour and curves. International
Conference on Biomedical Engineering and Informatics, pages 132136, 2008.
[4] T.M. Peters M.A. Audette, F.P. Ferrie. An algorithmic overview of surface registration
techniques for medical imaging. Medical image Analysis, 4:201217, 2000.
[5] Prof.Li Minglu Wan Rui. An overview of medical image registration. Proceedings
of the Fifth International Conference on Computational Intelligence and Multimedia
Applications (ICCIMA03), pages 385390, 2003.
[6] B.K. Puri A.D. Edwards K.K. Bhatia, J.V. Hajnal and D. Rueckert. Consistent group
wise non-rigid registration for atlas construction. IEEE International Symposium on
Biomedical Imaging, 1:1518, 2004.
[7] Varshney P.K. Min Xu.
16
17