A Novel Video Logging Method based on the Self-Focus Lens Array

Similar documents
Far field intensity distributions of an OMEGA laser beam were measured with

Analysis on Drill String Vibration Signal of Stick Slip and Bit Bouncing

Design of a digital holographic interferometer for the. ZaP Flow Z-Pinch

Corrosion Assessment of Offshore Oil Pipeline Based on Ultrasonic. Technique

OPTICAL SYSTEMS OBJECTIVES

Optical Performance of Nikon F-Mount Lenses. Landon Carter May 11, Measurement and Instrumentation

Demosaicing Algorithm for Color Filter Arrays Based on SVMs

IMAGE FORMATION. Light source properties. Sensor characteristics Surface. Surface reflectance properties. Optics

Astronomical Cameras

A Geometric Correction Method of Plane Image Based on OpenCV

Notes from Lens Lecture with Graham Reed

Synopsis of paper. Optomechanical design of multiscale gigapixel digital camera. Hui S. Son, Adam Johnson, et val.

LENSES. INEL 6088 Computer Vision

R.B.V.R.R. WOMEN S COLLEGE (AUTONOMOUS) Narayanaguda, Hyderabad.

Bias errors in PIV: the pixel locking effect revisited.

Research on Optical Fiber Flow Test Method With Non-Intrusion

Chapters 1 & 2. Definitions and applications Conceptual basis of photogrammetric processing

Fiberoptic and Waveguide Sensors

Adaptive Fingerprint Binarization by Frequency Domain Analysis

Introduction to DSP ECE-S352 Fall Quarter 2000 Matlab Project 1

EMC ANALYSIS OF ANTENNAS MOUNTED ON ELECTRICALLY LARGE PLATFORMS WITH PARALLEL FDTD METHOD

is a method of transmitting information from one place to another by sending light through an optical fiber. The light forms an electromagnetic

REAL TIME SURFACE DEFORMATIONS MONITORING DURING LASER PROCESSING

P Forsmark site investigation. RAMAC and BIPS logging in borehole HFM11 and HFM12

Chapter 18 Optical Elements

6.1 - Introduction to Periodic Functions

7. Michelson Interferometer

USE OF FT IN IMAGE PROCESSING IMAGE PROCESSING (RRY025)

Application Research on Hydraulic Coke Cutting Monitoring System Based on Optical Fiber Sensing Technology

Real-Time Scanning Goniometric Radiometer for Rapid Characterization of Laser Diodes and VCSELs

CCD Automatic Gain Algorithm Design of Noncontact Measurement System Based on High-speed Circuit Breaker

Research on 3-D measurement system based on handheld microscope

Visible Light Communication-based Indoor Positioning with Mobile Devices

An infrared-curtain-based sensor for non-contact recognition of raw-milk gas-liquid two-phase flows Li Boyi 1, Cheng Xiaodong 2

An Indian Journal FULL PAPER. Trade Science Inc. Parameters design of optical system in transmitive star simulator ABSTRACT KEYWORDS

Projection. Announcements. Müller-Lyer Illusion. Image formation. Readings Nalwa 2.1

Numerical simulation of a gradient-index fibre probe and its properties of light propagation

Geometry of Aerial Photographs

Suppression of Pulse Interference in Partial Discharge Measurement Based on Phase Correlation and Waveform Characteristics

Section 2 ADVANCED TECHNOLOGY DEVELOPMENTS

ENHANCEMENT OF SYNTHETIC APERTURE FOCUSING TECHNIQUE (SAFT) BY ADVANCED SIGNAL PROCESSING

CS 443: Imaging and Multimedia Cameras and Lenses

Bearing fault detection of wind turbine using vibration and SPM

A Compact Dual-Polarized Antenna for Base Station Application

MINIATURE X-RAY SOURCES AND THE EFFECTS OF SPOT SIZE ON SYSTEM PERFORMANCE

Computer Generated Holograms for Optical Testing

Study on a Single-Axis Fabry-Perot Fiber-Optic Accelerometer and its Signal Demodulation Method

GEOMETRICAL OPTICS Practical 1. Part I. BASIC ELEMENTS AND METHODS FOR CHARACTERIZATION OF OPTICAL SYSTEMS

ACOUSTIC BEHAVIOR OF MULTIPHASE FLOW CONDITIONS IN A VERTICAL WELL

CSE 527: Introduction to Computer Vision

International Conference on Information Sciences, Machinery, Materials and Energy (ICISMME 2015)

Image Measurement of Roller Chain Board Based on CCD Qingmin Liu 1,a, Zhikui Liu 1,b, Qionghong Lei 2,c and Kui Zhang 1,d

Used in Image Acquisition Area CCD Driving Circuit Design

3D radar imaging based on frequency-scanned antenna

Lecture 3 Complex Exponential Signals

Image Formation and Capture

A Recognition of License Plate Images from Fast Moving Vehicles Using Blur Kernel Estimation

定射角定方位射孔技术. Set Firing Angle Oriented Perforating Technology AMSTERDAM XI AN WUHUA JUNENG BLASTING EQUIPMENT CO., LTD

A 3D Profile Parallel Detecting System Based on Differential Confocal Microscopy. Y.H. Wang, X.F. Yu and Y.T. Fei

ARRAY PROCESSING FOR INTERSECTING CIRCLE RETRIEVAL

Ch 24. Geometric Optics

Background. Computer Vision & Digital Image Processing. Improved Bartlane transmitted image. Example Bartlane transmitted image

Projection. Readings. Szeliski 2.1. Wednesday, October 23, 13

UNIT-II : SIGNAL DEGRADATION IN OPTICAL FIBERS

A train bearing fault detection and diagnosis using acoustic emission

FRAUNHOFER AND FRESNEL DIFFRACTION IN ONE DIMENSION

Modeling and Synthesis of Aperture Effects in Cameras

Corrections for downhole NMR logging

Ihor TROTS, Andrzej NOWICKI, Marcin LEWANDOWSKI

Artifacts Reduced Interpolation Method for Single-Sensor Imaging System

ECC419 IMAGE PROCESSING

Chapter 25. Optical Instruments

Automatic inspection system for measurement of lens field curvature by means of computer vision

Fiber-optic Michelson Interferometer Sensor Fabricated by Femtosecond Lasers

Measurement of channel depth by using a general microscope based on depth of focus

Very short introduction to light microscopy and digital imaging

A Novel Algorithm for Hand Vein Recognition Based on Wavelet Decomposition and Mean Absolute Deviation

Experiment 1: Fraunhofer Diffraction of Light by a Single Slit

APPLICATION NOTE

Image Formation and Capture. Acknowledgment: some figures by B. Curless, E. Hecht, W.J. Smith, B.K.P. Horn, and A. Theuwissen

Be aware that there is no universal notation for the various quantities.

Supplementary Materials

Design of the Wide-view Collimator Based on ZEMAX

Coherent Laser Measurement and Control Beam Diagnostics

ELECTROMAGNETIC FIELD APPLICATION TO UNDERGROUND POWER CABLE DETECTION

Overview. Pinhole camera model Projective geometry Vanishing points and lines Projection matrix Cameras with Lenses Color Digital image

ON THE CREATION OF PANORAMIC IMAGES FROM IMAGE SEQUENCES

E X P E R I M E N T 12

Development of a new multi-wavelength confocal surface profilometer for in-situ automatic optical inspection (AOI)

A COMPACT DUAL INVERTED C-SHAPED SLOTS ANTENNA FOR WLAN APPLICATIONS

Faraday Rotators and Isolators

The Beam Characteristics of High Power Diode Laser Stack

Gravitational Lensing Experiment

Design of Omni-Directional Tilt Sensor Based on Machine Vision

Shaft Vibration Monitoring System for Rotating Machinery

Research Article Analysis and Design of Leaky-Wave Antenna with Low SLL Based on Half-Mode SIW Structure

1.Discuss the frequency domain techniques of image enhancement in detail.

Direct observation of beamed Raman scattering

MUSKY: Multispectral UV Sky camera. Valentina Caricato, Andrea Egidi, Marco Pisani and Massimo Zucco, INRIM

Active Aperture Control and Sensor Modulation for Flexible Imaging

Transcription:

Advances in Petroleum Exploration and Development Vol. 5, No. 2, 2013, pp. 75-82 DOI:10.3968/j.aped.1925543820130502.1393 ISSN 1925-542X [Print] ISSN 1925-5438 [Online] www.cscanada.net www.cscanada.org A Novel Video Logging Method based on the Self-Focus Lens Array JIA Huiqin [a],* ; LI Zhouli [a] [a] Key Laboratory of Education Ministry for Photoelectric Logging and Detecting of Oil and Gas, Xi'an Shiyou University, Xi'an, China. * Corresponding author. Supported by National Natural Science Foundation of China (No.51244009), the Education Department of Shaanxi Province of China (No. 2012JC21) and the Education Department Laboratory fund of Shaanxi Province of China (No. 12JS078). Received 21 April 2013; accepted 26 June 2013 Abstract At present, down-hole video logging method is used to observe the bottom and lateral well wall image, which places video camera on the bottom of logging instrument. This method can acquire the bottom image clearly. It s difficult to obtain lateral image because optic axis of lenses is placed along with the well axis, and it s impossible to place an existing camera along the radial direction because of restriction of borehole diameter or pipe diameter and object distance of camera, etc. A method and instrument for acquiring lateral image is presented on this paper, multiple self-focus lenses are placed along radial direction, used special relay lens transmitting the multiple imaging to an image sensor, and formed one image, then transmitted the image to ground. In order to use optic spectral properties, light cone circling the self-focus lenses is used to transfer the image, the minimum overlap radius for measurement boreholethe is analyzed,the seal problem for the optical system is designed.the lateral well wall image is acquired through researching the method of lateral multiple lens, then the phase correction method is used to fuse the image from different angles of lateral well wall. The video well logging instrument is developed using the above method, which can real-timely display the crack opening, filling substance, porosity and rock component of down-hole casing and borehole.so this method can provide an interpretation tools for the pipe internal or downhole phenomena. Key words: Lateral video logging; Image acquisition; Self-focus lens array; Circular light source; Phase correction method; Image fusion Jia, H. Q., & Li, Z. L. (2013). A Novel Video Logging Method based on the Self-Focus Lens Array. Advances in Petroleum Exploration and Development, 5(2), 75-82. Available from: URL: http://www.cscanada. net/index.php/aped/article/view/j.aped.1925543820130502.1393 DOI: http://dx.doi.org/10.3968/j.aped.1925543820130502.1393 INTRODUCTION The downhole video logging method commonly used in the logging instrument, which was placed on the bottom, used to observe the bottom and lateral images, to collect the image at the bottom of borehole [1-3]. When used for lateral imaging, due to the optical axis of the lens along the borehole axis direction, it s hard to obtain the lateral well wall image, because of the hole diameter or pipe diameter and the existing camera distance limit, radial direction place an existing camera is impossible. The main reason is that these instruments designed based on the normal camera technology,the diameters of normal lens is very large, place the camera in the bottom of well, which try to acquire the image on the bottom, but can not obtain the lateral well wall image. Since this proposal placed the main light source and imaging lens optical axis parallel to the axis of the borehole wall and sleeve through the plane at the bottom of the instrument (or convex) window irradiation bottom of the well, and received the reflected light. Obviously, such a scheme for identifying wells whereabouts matter, auxiliary salvage contributed. Since, near the reflection light in the viewing angle range, and distant reflection less attenuation of light and, in principle, cannot be placed to meet the direction of the well axis, and the well logging interpretation require for lateral imaging. The normal optical system just like the human eye, looking straight ahead, to see the front of the image, 75 Copyright Canadian Research & Development Center of Sciences and Cultures

A Novel Video Logging Method based on the Self-Focus Lens Array for the lateral can only get fuzzy and unclear images. In order to get a clear lateral borehole image, it s needed to disposit the lens facing to the lateral sidewall [4-5]. Our group study with Self-Focus Lens Array (SFLA) as objective lens, which placed in the lateral multiple objective lens, light cone is used to transfer image device, light source uses circular one, the light strength of which can be controlled by light source controller. Lateral multilens imaging system can overcome the deficiencies of the vertical camera difficult to get a clear image of the lateral borehole wall. Using the proposal, and compared to the previous proposal, image distortion can be eliminated. Also using the image mosaic method can form a complete image of a sectional lateral without the need for the mechanism rotates, the other can be avoided due to the lens movement appears ghosting [6-7]. 1. DESIGN PRINCIPLE FOR VIDEO LOGGING INSTRUMENT 1.1 Basic Principle of Video Logging Instrument The block diagram of Video logging instrument is shown in the Figure 1, which consist of two parts: downhole and ground part. The downhole consist of power module, Light Source Control Module, Light source, Optical system, DSP control module, CCD and modulation module. (1) Light source control module is used to adjust light strength. Light source is used to supply the light for Optical system. Light source is designed as circular one, and the LED is distributed around the light source board, the light strength is controlled by the digital potentiometer based on the I 2 C bus. (2) CCD is used to convert the light signal to electronic one, which transfer the borehole image to DSP control module, which is used to convert the differential LVDS signal to the TTL/CMOS one, and then transfer to the modulation module by the cameralink interface [8-9]. (3) Modulation module realizes the following functions: convert the cameralink interface to the normal well information coding format. Demodulation module realizes the functions of converting the normal well information coding format to the USB data format. (4) Power module is used to supply the power to light source control module and CCD. Which can supply the ±5V, ±3.3V to CCD,light source control module and modulation module. (5) DSP is used to realize the image compression because if sending the original image to ground, it will occupy the cable transmission bandwidth, simultaneously it can be used to control the light source strength of light source. (6) Optical system is the core of this system. The image is clear or not decides it. The composition block diagram is shown in the Figure 2. (7) Software is used to display the underground image, Interpret the crack opening, filling substance, porosity and rock component of down-hole casing. Which is realized by the graphical display software LabVIEW 8.6. Figure 1 Basic Diagram for Video Logging Instrument 1.2 Design of Optical System SFLA is the most important component of Optical system, and the most important properties of it is smallness, the diameter of it is very small. Light cones is used to transfer light. In the laterally placed some SFLA, in the well axis placed only one self-focus lens. Copyright Canadian Research & Development Center of Sciences and Cultures 76

JIA Huiqin; LI Zhouli (2013). Advances in Petroleum Exploration and Development, 5(2), 75-82 Optical system is divided into two cavities, the above cavity set two light cone fixed plates, which placed up and down, and the center of which place one light cone parallel to the well axis, it can be defined as Horizontal Light Cone (HLC), and 90 degree bending light cone around the LCFP, it can be defined as Vertical Light Cone (VLC) and the small side of which coupled to SFLA, the side of SFLA face to the tested well wall through the hole on the shell of optical system (SPS). The small end of HLC is coupled to self-focus lens, the other side of self-focus lens face to the bottom of the borehole.vlc face the well wall, can obtain the image of well wall. The two Light cone fixed plate (LCFP) is sheathed on the outer shell of the optical system. And the Circular Glass Shell (CGS) is used to isolate the optical system from the outside well environment. (b) The Orthographic Views of VLC Figure 3 Orthographic Views of Light Cone Figure 2 Front View of Video Logging Instrument The most important part for single lens imaging system is coupling and image transfer. Light cones are used to transfer light. It s very difficult to mechanical processing. The orthographic views of 90 degree bending Light cones is shown in Figure 3. (a) The Orthographic Views of HLC 1.3 Minimum Overlap Radius for Measurement Borehole SFLA is used to obtain the whole information of borehole, so it will cause the overlap between the two self-focus lens. With the increase of numerical aperture, can see the scale will be bigger, but the gradient refractive index lens edge will appear serious fuzzy phenomenon, caused very big effect on imaging, and so select the appropriate gradient refractive index lens is particularly important. The method for computation of minimum overlap radius is as follows: The eight self-focus lens is equally distributed on the SPS, so the angle among each eight self-focus lens are 45.The radius of SPS is d,θ is the field angle of selffocus lens, R2 is the distance from the edge of self-focus lens to sidewall, r1 is the radius of well. θ 1 is the angel between axial direction of self-focus lens and, is 22.5. r 1 value is obtained according to the cosine theorem.there will compute under three cases, the result is shown in the Table 1. Seen from the table.1, when the parameter of d is determined, the numerical aperture N A= sin(θ ) is larger, the overlapping radius is small, so obtained image information is much. With the increase of the numerical aperture N A, so the range of seen information is large, but the edge of selffocus lens is fuzzy, it influences the imaging, so it s important to select the right self-focus lens [10]. 77 Copyright Canadian Research & Development Center of Sciences and Cultures

A Novel Video Logging Method based on the Self-Focus Lens Array Image region of single lens on the CCD area is shown in the figure 6. Seen from figure 6, there are nine Imaging region, In which the center of the circular region for the well axis is the direction of the imaging area, other eight ones is the lateral bearing imaging area, all the images come from the self-focus lens through the light cones. Figure 4 Schematic Diagram of Computation the Minimum Overlapping Radius Table 1 r 1 Value Under Different Cases 10 12 14 d (mm) sin(θ) r 1 (mm) 0.46 72.2613 0.552 41.5204 0.6 34.9051 0.46 86.7136 0.552 49.8245 0.6 41.8861 0.46 101.1659 0.552 58.1286 0.6 48.8672 1.4 Seal Design for Optical System It needs to consider the seal problem because the logging instrument is placed in the downhole, there will have muds or other materials that can pollute the lens, the circular ring glass window is designed to avoid the lens to be polluted. The mechanical structure is shown in the Figure 5.The field angle of one self-focus lens is 2 θ. Figure 6 nine Imaging region of CCD 2. IMAGE PROCESSING AND INFORMATION INTERPRETATION 2.1 The Borehole Image Processing Principle Seen from Figure 5, the borehole image from the different angles has some overlap, in order to cut the overlap part, and form the registration image for well wall. Then it s needed to process the image from the SFLA. In order to realize image registration, it s needed to realize image translation, rotation and scaling. The process is shown in the Figure 7. Downhole Image from the CCD is firstly needed to registration, Panoramic and fusion; finally form the GIS map for the downhole. 2θ Figure 5 Distruibued Relation of VLC Figure 7 Flow Diagram of Image Fusion Copyright Canadian Research & Development Center of Sciences and Cultures 78

JIA Huiqin; LI Zhouli (2013). Advances in Petroleum Exploration and Development, 5(2), 75-82 2.2 Process of Image Translation Assumption the two images: f 1 (x,y), f 2 (x,y), and there are only translation relation between the two image, and the translation amount is (x 0, y 0 ). f 1 (x,y) = f 2 (x-x 0,y-y 0 ) (1) According to the translation property of Fourier transform, there are such relations: F 2 (u,v) = e -j(ux 0 +vy 0 ) F 1 (u,v) (2) Where F 1 (u,v) and F 2 (u,v) are the Fourier transform of f 1 (x,y) and f 2 (x,y) respectively. The cross spectrum of them is as follows: * F1( uvf, ) 2 ( uv, ) j( ux0+ vy0) = e * 3) F1( uvf, ) 2 ( uv, ) where F * 2(u,v) is the complex conjugate of F 2 (u,v), the fourier transform of e -j(ux 0 +vy 0 ) is a two-dimensional pulse function, the function name is δ(x-x 0,y-y 0 ). And phase correlation is to get the inverse fourier transform of equation (3), then search the peak value,to decide the registration parameter(x 0,y 0 ). The step of image translation is as follows: (1) Acquire images contained with the different angels through the optical system; (2) execute the FFT algorithm to the acquired images; A = fft(f 1 x,y)b = fft(f 2 x,y)) (4) (3) compute the cross power spectrum of images; (4) compute the inverse cross power spectrum aiming at the step (3); C = A.*conj(B)/norm(A.*conj(B)) (5) (5) Search the maximum amplitude, then which is the translation amount. Figure 8 is result of estimating the parameters through the phase correlation method. Figure 8(a) is the reference image, Figure 8(b) is registration Image, Figure 8(c) is Fourier Transform of Cross Power Spectrum, Figure 8(d) is splicing image. (b) Registration Image (c) Fourier Transform of Cross Power Spectrum (d) Splicing Image Figure 8 stimation of Registration Parameter Using Phase Correlation Method (a) Reference Image 2.3 Process of Image Rotation The steps of image rotation and scaling is as follows: Assumption the two images: f 1 (x,y) and f 2 (x,y), is obtained from f 1 (x,y) to scale the image with λ parameters and rotate one with θ 0 counterclockwise, and meets: 79 Copyright Canadian Research & Development Center of Sciences and Cultures

A Novel Video Logging Method based on the Self-Focus Lens Array 1 1 f2( xy, ) = f1[ λ ( xcosθ0 + ysin θ0), λ ( sinθ0 + ycos θ0)] (6) After logarithmic polar coordinate transformation, then f2(ln ρθ, ) = f1(ln ρ ln λθ, θ0) (7) Image Descartes coordinate rotation and scaling has been converted to the logarithmic polar coordinate on the translation relation, in the log polar coordinate system using phase correction method. Figure 9 realizes the image registration, Figure 9(b) is the registration image after Figure 9(a) rotates 17.36, and zooms 1.37 times, at the same time which cut in the center at the same size, realize the image registration using the phase correction method. 2.4 Mosaic Image After Registration (a) Raw Image (a) Reference Image (b) Zoom and Rotation Image (b) Registration Image (c) Logarithmic Polar Coordinates of Figure (a) (c) After Registration Image (d) Logarithmic Polar Coordinates of Figure (b) Figure 9 Mosaic Effect Diagram After Registration (d) Mosaic Image Figure 10 Mosaic Image After Registration Copyright Canadian Research & Development Center of Sciences and Cultures 80

JIA Huiqin; LI Zhouli (2013). Advances in Petroleum Exploration and Development, 5(2), 75-82 Figure 10 is the mosaic image for the two downhole image. Figure 10(a) is the reference image, Figure 10(b) is registration Image, Figure 10(c) is the image after registration, and Figure 10(d) is mosaic image. 2.5 Fusion Image Using Phase Corection Method Image fuses using the pixel maximum method, assumption that the size of two images is m n, the image after fusion is f, then aiming at the pixel Gray scale value,choose the maximum one to execute the image fusion, the equation is: f(m,n) = max{a(m,n),b(m,n)} (8) m,n is the row and column of image.when fusion is executed, comparison of the size of pixel value corresponding place in the raw image, and taken the maximum of pixel value as pixel value in the position (m,n) of image F. Taken the radial 8 self-focus lens images for image registration. The difference between downhole 1 and downhole 2 is diameter of borehole. Figure 11 is Mosaic Experimental Result of downhole 1. Figure 11(a) is the reference image, Figure 11(b) is registration image, Figure 11(c) is a Mosaic image, Figure 11(d) is the fusion image. (d) Fusion Image Figure 11 Mosaic Experimental Result of Downhole 1 Figure 12 is Mosaic Experimental Result of downhole 2. Figure 12(a) is the reference image, Figure 12(b) is registration Image, and Figure 12(c) is a Mosaic image, Figure 12(d) is the fusion image. (a) Reference Image (a) Reference Image (b) Registration Image (b) Registration Image (c) Mosaic Image (c) Mosaic Image 81 Copyright Canadian Research & Development Center of Sciences and Cultures

A Novel Video Logging Method based on the Self-Focus Lens Array (3) Because the image is obtained from different angles and the acquired image has overlap. In order to form the whole well wall information map, phase correction method is used to realize image fusion, and get the good result. (d) Fusion Image Figure 12 Mosaic Experimental Result of Downhole 2 Seen from Table 2, it shows the Estimated Parameters for downhole 1 and downhole 2. Using this parameters can make the image from different angles to one GIS map, which can display the information directly, which provide a novel method for well logging such as crack opening, filling substance, porosity and rock component of down-hole casing, or casing casing damage inspection. Table 2 Two Experimental Result Image Estimated Parameters [ x0, y 0, θ, 0 λ ] downhole 1 [ 113,37,17.27,1.336 ] downhole 2 [ 80,18,17.31,1.346 ] CONCLUSIONS This paper researches to use the SFLA to acquire the signals in the process of well logging instrument design, using this method can obtain the image of well wall. The main work of this paper is below: (1) Because the size of SFLA is about 1.5 mm.it is very suitable for acquiring the well wall information, because the diameter of borehole is about the 77 mm. (2) Using circular light source provides a new light principle for downhole image acquiring. Through test, circular light source has good effect on the acquiring the well wall image information, and improve the definition of well wall image. REFERENCES [1] Li, Z. L. et al. (2011). Design and Implementation of the Annular Light Source and Its Controller for Downhole Lateral Multi-Lens Imaging System. Journal of Xi an Shiyou University, 4, 4. [2] Li, Z. L., Jia, H. Q., & Wang, P. et al. (2010). A Multiple Lens Compound Imaging Video Logging Device. China Patents, 60877011.11.2010. [3] Ren, A L. (2010). Research and Implementation on the Control Method of the Underground Lateral and Annular Light Source Intensity. Electronic Test, 12, 4. [4] Wang, W., & Lu, P. Z. (2006). An Overview of Digital Image Mosaic. Journal of Chinese Computer Systems, 27, 1348. [5] Li, Z. X., Mao, Y. B., & Wang, Z. Q. (2005). A Method of Image Mosaicing Using Log Polar Coordinate Mapping. Journal of Image and Graphics, 10, 60. [6] Hu, S. J., Ge, X. W., & Su, L. (2007). Panoramic Image Mosaic Based on Domain Correlation. Journal of Hefei University of Technology, 30, 60. [7] Chen, Q., Dfrise, M., & Deconinck, F. (1994). Symmetric Phase-only Matched Filtering of Fourier-Mellin Transform for Image Registration and Recognition. IEEE Transactions on Pattern Analysis and Machine Intelligence, 16, 1154. [8] Fang, J. W., Zhai, C., & Jin, Y. (2008). Simulation and Optimization of Phase Correlation Algorithm in Image Mosaics. Computer Applications and Software, 16, 47. [9] Wei, C. T., Wu, P., Zhang, Z. X., & Zhang, J. Q. (2011). An Image Matching Algorithm Based on Improved Log-Polar Image Transform. Bulletin of Surveying and Mapping, 15, 54. [10] Zheng, Z. B., & Ye, Z. F. (2006). Image Registration Algorithm Based on Phase-Correlation. Journal of Data Acquisition & Processing, 21, 444. Copyright Canadian Research & Development Center of Sciences and Cultures 82