Image Based Measurement of Length and Distance of an Inaccessible Object

Similar documents
Measuring Height of an Object using Accelerometer and Camera in ios and Android Devices

LIGHT: Two-slit Interference

An Image Based Approach to Compute Object Distance

Estimation of Suitable Grade Value for Stopping Sight Distance Computation on Vertical Curves

Synthesis of Planar Mechanisms, Part IX: Path Generation using 6 Bar 2 Sliders Mechanism

1. What is the law of reflection?

Algebra Area of Parallelograms

Part Images Formed by Flat Mirrors. This Chapter. Phys. 281B Geometric Optics. Chapter 2 : Image Formation. Chapter 2: Image Formation

Ch. 25 The Reflection of Light

Construction of 3D Model of a Rectangular Parallelopiped Object from Its 2D Images

Experiment 8 Wave Optics

Fine Classification of Unconstrained Handwritten Persian/Arabic Numerals by Removing Confusion amongst Similar Classes

At the interface between two materials, where light can be reflected or refracted. Within a material, where the light can be scattered or absorbed.

Learning Task: Exploring Reflections and Rotations

Effects Of Shadow On Canny Edge Detection through a camera

An Efficient Single Chord-based Accumulation Technique (SCA) to Detect More Reliable Corners

Human Identification at a Distance Using Body Shape Information

Review for Spring Final Exam Geometry 1. Classify the figure. Name the vertices, edges, and base.

Assignment 1 Due September 6, 2011

Chapter 12 Notes: Optics

Downloaded from UNIT 06 Optics

ON THE VELOCITY OF A WEIGHTED CYLINDER DOWN AN INCLINED PLANE

All forms of EM waves travel at the speed of light in a vacuum = 3.00 x 10 8 m/s This speed is constant in air as well

Edge and local feature detection - 2. Importance of edge detection in computer vision

NEW MONITORING TECHNIQUES ON THE DETERMINATION OF STRUCTURE DEFORMATIONS

Optics II. Reflection and Mirrors

MOVING A VERTEX. How must Kirsten move vertex A to be sure that the area of her new triangle is

Chapter 3 Image Registration. Chapter 3 Image Registration

Light and the Properties of Reflection & Refraction

Appendix E. Plane Geometry

DC2 File Format. 1. Header line 2. Entity line 3. Point line 4. String line

Factor Quadratic Expressions

Name Homeroom. MSMM2014

Physics 101, Lab 1: LINEAR KINEMATICS PREDICTION SHEET

ADVANCED SURVEYING (S-II)

d = (x 2 - x 1 ) 2 + (y 2 - y 1 ) 2 Student Name: Date: Teacher Name: Sunil Dudeja Score:

Answer Key: Three-Dimensional Cross Sections

ACTIVITY TWO CONSTANT VELOCITY IN TWO DIRECTIONS

Horizontal and Vertical Curve Design

6th Grade Vocabulary Mathematics Unit 2

Physical Optics. You can observe a lot just by watching. Yogi Berra ( )

SREE VAHINI INSTITUTE OF SCIENCE AND TECHNOLOGY

pine cone Ratio = 13:8 or 8:5

Developing an intelligent sign inventory using image processing

Diffraction and Interference Lab 7 PRECAUTION

Optics INTRODUCTION DISCUSSION OF PRINCIPLES. Reflection by a Plane Mirror

Experiment 6. Snell s Law. Use Snell s Law to determine the index of refraction of Lucite.

LANE DEPARTURE WARNING SYSTEM FOR VEHICLE SAFETY

Maximum and Minimum Problems

UNIVERSITI MALAYSIA SARAWAK FACULTY OF ENGINEERING CIVIL ENGINEERING DEPARTMENT

A Survey of Light Source Detection Methods

2/26/2016. Chapter 23 Ray Optics. Chapter 23 Preview. Chapter 23 Preview

Geometry Spring Final Review #1, 2014

Interference and Diffraction

4. Describe the correlation shown by the scatter plot. 8. Find the distance between the lines with the equations and.

Checkpoint 1 Define Trig Functions Solve each right triangle by finding all missing sides and angles, round to four decimal places

Simulation of Turbulent Flow over the Ahmed Body

Calculation of the Focal Length of an Offset Satellite Dish Antenna

Chapter 2: Frequency Distributions

PHYSICS. Chapter 34 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT

9. Polarization. 1) General observations [Room 310]

You ll use the six trigonometric functions of an angle to do this. In some cases, you will be able to use properties of the = 46

Reflection of Light. 1)Students will discover how light interacts with certain types of surfaces

MPM 1D Learning Goals and Success Criteria ver1 Sept. 1, Learning Goal I will be able to: Success Criteria I can:

Chapter 6 Review. Extending Skills with Trigonometry. Check Your Understanding

Segmentation of Mushroom and Cap Width Measurement Using Modified K-Means Clustering Algorithm

Subpixel Corner Detection Using Spatial Moment 1)

SUMMARY. CS380: Introduction to Computer Graphics Projection Chapter 10. Min H. Kim KAIST School of Computing 18/04/12. Smooth Interpolation

(a) rotating 45 0 about the origin and then translating in the direction of vector I by 4 units and (b) translating and then rotation.

5B.4 ~ Calculating Sine, Cosine, Tangent, Cosecant, Secant and Cotangent WB: Pgs :1-10 Pgs : 1-7

Monocular Vision Based Autonomous Navigation for Arbitrarily Shaped Urban Roads

Components of Alignment. Horizontal Alignment. Vertical Alignment. Highway Design Project. Vertical Alignment. Vertical Alignment.

Facial Feature Extraction Based On FPD and GLCM Algorithms

The Mathematics of Highway Design

A deformable model driven method for handling clothes

PHYS1004 Problem Sheet - Optics - with Solutions

Chapter 32 Light: Reflection and Refraction. Copyright 2009 Pearson Education, Inc.

Distance in Coordinate Geometry

A NON-TRIGONOMETRIC, PSEUDO AREA PRESERVING, POLYLINE SMOOTHING ALGORITHM

Math 124 Final Examination Autumn Turn off all cell phones, pagers, radios, mp3 players, and other similar devices.

What is it? How does it work? How do we use it?

Prepare a stem-and-leaf graph for the following data. In your final display, you should arrange the leaves for each stem in increasing order.

Preceding vehicle detection and distance estimation. lane change, warning system.

ACTIVITY 2: Reflection of Light

AREA Judo Math Inc.

Geometry Second Semester Final Exam Review

Big Mathematical Ideas and Understandings

diverging. We will be using simplified symbols of ideal lenses:

Slide 2 / 222. Algebra II. Quadratic Functions

Keywords: Thresholding, Morphological operations, Image filtering, Adaptive histogram equalization, Ceramic tile.

DR. UMMA HANY. Current Status: Education: Career Summary:

Geometry SIA #2. Name: Class: Date: Multiple Choice Identify the choice that best completes the statement or answers the question.

Geometry 10 and 11 Notes

Be sure to label all answers and leave answers in exact simplified form.

Geometry Second Semester Review

9. RAY OPTICS AND OPTICAL INSTRUMENTS

Refraction of Light Finding the Index of Refraction and the Critical Angle

Learning Task: Exploring Reflections and Rotations

MET71 COMPUTER AIDED DESIGN

Dispersion (23.5) Neil Alberding (SFU Physics) Physics 121: Optics, Electricity & Magnetism Spring / 17

Transcription:

Image Based Measurement of ength and Distance of an Inaccessible Object Md. Imdadul Islam 1, Mahbubul Alam 2, Sarnali Basak 3, Md. Al-Amin Khandaker 4, S M Raihan Azim 5 Abstract Sometimes it is necessary to measure the length of an object or its distance from a particular location when it is difficult to access the location of the object even not at all accessible. In this paper we proposed an image based measurement technique of such an object and its distance. The technique is very simple, only a digital camera, a wooden bar and three ASERs of intense beam are necessary. Entire measurement is taken ignoring the focal length of the camera but the result provides above 95% accuracy. Such measurement technique; for example is essential for a civil engineer to take approximate measurement before work plan even for a zoologist to measure the size of a wild animal which cannot be caught. Keywords Distance Measurement, Pixel-based ength Measurement, OS, Numerical Integration 1. Introduction If someone is asked to measure the size of an object or its distance from a particular point, the most realistic and easiest thinking is to use a measuring tape. If the length of the object or its distance is very large then civil engineering survey is the optimum solution. Although such survey provides very accurate results but sometimes it is necessary to take an approximate initial measured data for preplan (before mobilization surveying equipment). In this paper work we introduce very low cost equipment, based on image processing technique at the expense of accuracy. Md. Imdadul Islam, Computer Science and Engineering, Mahbubul Alam, Computer Science and Engineering, Sarnali Basak, Computer Science and Engineering, Md. Al-Amin Khandaker, Computer Science and Engineering, S M Raihan Azim, Computer Science and Engineering, 123 Very few research papers are found pertinent to this work. Only four papers are found (not very closed to our analysis) about measurement of length of a fish. The basic principle of these papers are discussed here as the previous literature. The measurement of distance of an object covering this paper is a newly proposed model. In [1] first of all, the edge of the fish is determined which is then splitted into 100 equidistance segments. The center and principle axis of the resultant polygon is determined using moment integration. The result of the moment provides the actual orientation of the fish. The next step is to draw a set of horizontal and vertical lines with respect to the principle axis. Finally the length of the fish is measured along the central axis of the fish taking 8 points on the axis. A cost effective and faster method of measurement of length of fish is proposed in [2]. First of all preprocessing is done on the image to detect the head and tail of fish using canny edge detection, scale-map and sub-segments using a relevance model of [3]. ength of the fish is then measured considering, bending ratio plot, corner factor, focal length (f) and mean object distance (O). Entire procedure is given in block diagram in [3] and the authors of claims 95% accuracy of their result. In [4] the author develops software named AnalyzingDigitalImages to measure the length of an object, using its image. The paper is organized like, section 2 provides the methodology of the proposed model with mathematical formulation, section 3 deals with image based results on several objects, and finally section 4 concludes the entire analysis. 2. Methodology The procedure regarding to the implementation of physical equipment, extraction of the image and the corresponding measurements can be explained with following seven steps: Step 1. Two ASERs of intense beam are placed at 1ft apart on a wooden bar. They are attached on the bar with adhesive material so precisely that the beams emitted are parallel. Step 2. aser beams are projected on an arbitrary object. The object (the major axis of the object) and

the beam should be (approximately) mutually perpendicular. Step 3. Image of the object (including the beam spot) is taken by high resolution digital camera. The camera should be placed on the wooden bar at the midpoint between the ASERs. Step 4. Image is extracted using imtool( file name of the image with extension ) function under MATAB 9.x or above version. Step 5. Using pixel based distance measurement tools, the size of the object is measured according to fig.2 and corresponding condition is: i is the maximum when θ i 0 and this maximum length is the actual length of the object. bar where the separation between them is 1ft like step-1, and the third beam makes an angle α (we recommend to take α below 10 0 ) with the direction of first two parallel beams like fig.3. Step 7. Image of the object (including the three beam spots) is taken by high resolution digital camera. The image is extracted like step-4 under MATAB software. Here both and are known therefore l can be calculated from pixel based measurement tool. Now, the distance of the object, l x. tan ASER-1 ASER-2 ASER-3 aser-1 aser-2 α x θ max (a) θ -10 0 0 0 10 0 (b) θ Fig.3 Distance measurement apparatus l 3. Results In this paper work we first considered a smooth object of uniform and periodic structure. The line of sight (OS) of the camera is placed parallel to horizontal where as the object is perpendicular to the OS of camera. We used high resolution camera of, 5184 3456 pixels and of size 5.7MB. We first select a vertical wall like fig.4 and the snap was taken 10ft ahead of the wall. The width of bricks is measured in number of pixels. The widths are found almost same for all the bricks. The small variation occurs because of angular variation of sides of the object and irregular width of adhesive element between bricks. Fig. 2(a)&(b) Measurement of size of an object To measure the distance between the object and the location of experiment we have to follow the following procedure: Step 6. Three ASER beams are projected on the object two beams are perpendicular to the wooden 124 Fig.4 Non-uniform wall

width of two bricks in pixels International Journal of Advanced Computer Research (ISSN (print): 2249-7277 ISSN (online): 2277-7970) The situation is complexly different if the OS of the camera is slant w.r.t. the horizontal. Fig.5 shows more precise brick wall where the OS of the camera is placed 15 0 with respect to the horizontal plane. The variation of width of the bricks is visualized from the figure. To observe the profile of such variation, the width of bricks is plotted against the distance in number of bricks shown in fig.3. Fitting the data with 3 rd order polynomial, it is found that the deviation in width with distance is almost linear. Fig.6(b). Door with uniform curving Fig.5 (a) Uniform wall et us now apply the above technique on the most uniform structure of Jahangirnagar university campus which is the Shahid Minar. For the case of horizontal OS a snap is taken 50ft away from it and the variation width of bricks from top to bottom is found almost same like fig.7. 88 86 84 82 80 78 76 74 72 Measured data Polynomial fit 70 0 10 20 30 40 50 60 70 80 distance in unit of bricks Fig.5(b) Width of bricks vs the distance of the number of bricks The phenomena of slant OS is also visualized from fig.6 where grill of window and door of straight curving are used. In both cases the snap is taken 10 ft away from the object at an angle of 15 0 w.r.t. horizontal. Fig.6(a) Uniform window grill 125 Fig.7 OS is horizontal Finally three snaps of the Shahid Minar is taken 30ft away from it at 0 0, 15 0 and 30 0 w.r.t. horizontal line. A closed view of the image of 30 0 case is shown in fig.8(a) where the variation in length of bricks are visualized. The profile of pixels/brick vs. distance along the object is plotted for three images shown in fig. 8(b). All the curves are found linear. The variation of width of a brick is found very small within 50ft when the OS of the camera is perpendicular to the object. Since the bricks are squeezed linearly, applying the regression of 3 rd order polynomial we got the constant rate of squeeze, S = 0.965. Therefore the distance measured by the proposed technique should be divided by the factor S.

Fig.8(a) OS is 30 0 with respect to horizontal Fig.10 Height measurement of a person et us introduce some application example pertinent to the paper. The first example is the measurement of the height of a man like fig.10. The measured height of the man, h = (4554.4/762.38) 1 = 5.97ft; whet the actual height is 5ft and 11 inches. Therefore the accuracy is well satisfactory. Our second example is the measurement of height of a column of a building shown in fig.11. Fig.8(b) The variation of width of a brick w.r.t. for different height Applying the ASER spot at 1ft apart on the Shahid Minar we got the distance between them equal to 70 pixels and the height of the tower above the platform is found 4508 pixels shown in fig.9. Now the height of the tower, h = 4508.92/70 = 64.413ft. Using the factor S, the true height will be, 64.413/S=66.749ft. This factor can be ignored for measurement of a small object like fig.10 and 11. Fig.11 Height measurement of a pillar of a building The measured height of the column= (4080/282)/1.5 = 10.79ft; where the exact height = 10.5ft. Our third example is the measurement of length of a curved body based on piece wise linear model like fig.12. Smaller the segments more accurate the result will be like the trapezoidal method of numerical integration. Fig.9 Measurement of height 126 Fig.12 Piece wise linear approximation

video on the front side of the moving car (with two aser spots). Now from the variation of number of pixel between two aser spots (from consecutive video frame) per unit time will give the indication of its velocity. References [1] D.J. White, C. Svellingen, N.J.C. Strachan, Automated measurement of species and length of fish by computer vision, Elsevier Fisheries Research 80 (2006) 203 210 [2] Mohd Rahim, Norhaida Abdullah, Ismail Amin, Mohammad Zakaria, Mustafa Man, Nur Othman, A New Approach in Measuring Fish ength Using FieDI Framework, IAJIT, 2010 [3] Kiranyaz S., Hantao iu, Ferreira M., and Gabbouj M., An Efficient Approach for Boundary Based Corner Detection by Maximizing Bending Ratio and Curvature, in Proceedings of 9th International Symposium on Signal Processing and its Applications, Sharjah, UAE, pp. 1-4, 2007 [4] John Pickle, Measuring ength and Area of Objects in Digital Images Using AnalyzingDigitalImages Software, Concord Academy, March 19, 2008, http://mvh.sr.unh.edu/software/software.htm Fig.13 (a) and (b) Three points model to measure the distance of an object et us now apply three points model to measure the distance of an object (wall of a room is selected to get distinct aser spot on the image) like fig.5. In fig.13(a) the geometry of right angle triangle is found like, = 2028.89/834.19 = 2.43216ft, l = 1f and α=10 0. Now the distance measured by apparatus, x = /tan(α) = 2.43216/0.173 = 13.79ft. The exact length is found 14.5ft. Now the error, ε = (14.5-13.79/14.5) = 4.89%. From fig. 13(b) we got the error ε = 4.63%. 4. Conclusions aser spots are found very blur on the image while taking readings on day light. Therefore very intense aser beam is necessary to take reading on day light. Edge detection technique or contrast enhancement can also be applied (followed by some preprocess) to locate the blur aser spot. The concept of the paper can be applied to measure the length of a straight road from the profile of variation of its width with distance. The entire work can be extended to measure the velocity of a vehicle along a straight road taking 127 Md. Imdadul Islam has completed his B.Sc. and M.Sc Engineering in Electrical and Electronic Engineering from Bangladesh University of Engineering and Technology, Dhaka, Bangladesh in 1993 and 1998 respectively and has completed his Ph.D degree from the Department of Computer Science and Engineering, Jahangirnagar University, Dhaka, Bangladesh in the field of network traffic engineering in 2010. He is now working as a Professor at the Department of Computer Science and Engineering, Jahangirnagar University, Savar, Dhaka, Bangladesh. Previously, he worked as an Assistant Engineer in Sheba Telecom (Pvt.) TD (A joint venture company between Bangladesh and Malaysia, for Mobile cellular and W), from Sept.1994 to July 1996. Dr Islam has a very good field experience in installation of Radio Base Stations and Switching Centers for W. His research field is network traffic, wireless communications, wavelet transform, OFDMA, WCDMA, adaptive filter theory, ANFIS and array antenna systems. He has more than hundred research papers in national and international journals and conference proceedings. He completed his B.Sc. (Hons.) in Computer Science and Engineering from Jahangirnagar University, Savar, Dhaka, Bangladesh in 2006. He

received the Asadul Kabir Gold Medal for securing the highest grade in the B.Sc.(Hons.) final examination 2006 at Jahangirnagar University. He is now working as a ecturer, at the department of Computer Science and Engineering, Jahangirnagar University, Savar, Dhaka, Bangladesh from 13th July 2009. His research interests are wavelet transform, wireless communications, adaptive filter theory, artificial intelligence, image processing and machine learning. Sarnali Basak received her B.Sc. (Hons.) and M.S. degrees in Computer Science and Engineering from Jahangirnagar University in Savar, Dhaka, Bangladesh in 2008 and 2010, respectively. She is working as an assistant professor in the Department of Computer Science and Engineering at University of Asia Pacific, Dhanmondi, Dhaka-1209. She received an NSICT fellowship from the Ministry of Science and Information & Communication Technology, Bangladesh for her MS research in 2011. Ms. Basak has also several research papers in reputed international journals. Her research fields are Image Processing, Pattern Recognition, Artificial Intelligence, Neural Networks and Biometrics. Md. Al-Amin Khandaker completed his Bachelor in Computer Science and Engineering from Jahangirnagar University, Dhaka, Bangladesh in 2012. He is pursuing Masters in Computer Science and Engineering from the same University. He has been working as a Software Engineer at Metatude Asia td, Dhaka (A Dutch offshore software company) since May, 2012. He has expertise in mobile application development. His research interest includes Wireless Communication, Software Engineering, Image Processing and Human Computer Interaction. S M Raihan Azim has completed his B.Sc (Hons.) in Computer Science and Engineering from Jahangirnagar University, Dhaka, Bangladesh in 2012. He is now working as a software Engineer in Vizrt Bangladesh td. His research interests are in Wireless Communication, Artificial Intelligence, Software Engineering and Programming anguages. 128