EP1008105A1 - Method and arrangement for finding the centroid of a body - Google Patents

Method and arrangement for finding the centroid of a body

Info

Publication number
EP1008105A1
EP1008105A1 EP98901151A EP98901151A EP1008105A1 EP 1008105 A1 EP1008105 A1 EP 1008105A1 EP 98901151 A EP98901151 A EP 98901151A EP 98901151 A EP98901151 A EP 98901151A EP 1008105 A1 EP1008105 A1 EP 1008105A1
Authority
EP
European Patent Office
Prior art keywords
image
segment
marker
centre
length
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
EP98901151A
Other languages
German (de)
English (en)
French (fr)
Inventor
Thorleif Josefsson
Göran FRIDH
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qualisys AB
Original Assignee
Qualisys AB
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Qualisys AB filed Critical Qualisys AB
Publication of EP1008105A1 publication Critical patent/EP1008105A1/en
Withdrawn legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/66Analysis of geometric attributes of image moments or centre of gravity
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/64Analysis of geometric attributes of convexity or concavity
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30196Human being; Person

Definitions

  • the present invention relates to a method and arrangement for determining at least one dimension data for a body, specially data relating to centroid coordinate of an image of said body reproduced on a sensor that produces a signal containing pixel data for said reproduced image.
  • markers are attached to the object to be analysed.
  • the object provided with the markers was first filmed and then manually analysed and digitalised to determine the correct position of the markers. This was a time-consuming procedure.
  • CCD sensor which is a light sensor
  • a CCD sensor consisting of lines of charged coupled sensors arranged as a matrix, i.e. arranged in an X and Y coordinate system, for one or several colours, converts the light (from the optical element) projected on it, by electronically scanning in Y direction each line of X sensors and producing a television (video) signal. Then, the signals may be analysed in different ways to detect the position of the markers attached to the object.
  • U.S. Patent 5,077,294 describes a motion analyse system using a camera converting the image of a body having illuminated markers placed thereon into a plurality of image frames, each comprising a plurality of image lines. Each image line comprises a plurality of pixels. Then an image processor is used to receive each image line and detect the occurrence of a marker in each image line. This is done by comparing each pixel amplitude to a threshold value and indicating when a pixel amplitude value exceeds the threshold value. Then an offset value representing the difference between the pixel amplitude value and the threshold value is added to a running total of amplitude values for the detected marker.
  • the offset value is to calculate the moment value for each pixel in each marker based on the position of the pixel in the scan line.
  • the calculated moment line value is then added to a running total of amplitude values for the detected marker.
  • the amplitude and the moment sums for each series of pixels for each marker in each image line are passed to a signal processor, which groups the amplitude and the moment sums for each marker together, then the groups' values are used to calculate a centroid value for each marker.
  • Another parameter is determined by a method characterized in steps of: interpolating said signal and calculating a length of a segment representing a section of said image, and using said calculated length to determine the area of the image by using summing all length.
  • the invention also relates to an arrangement for determining at least one dimension data for a body, comprising: means for collecting image data for said body, converting means for producing a digital signal containing pixel data for said reproduced image, comparator unit for generating difference between a pixel signal level and a threshold value, computation unit for calculating a centre point for length of each segment representing a section of said image, calculating a centroid coordinate and/or an area of said image and/or a radius of said image.
  • Fig. 1 is a schematic diagram of a simple motion analyse system according to the invention.
  • Fig. 2 is a block diagram over the primary functions of the camera according to fig. 1.
  • Fig. 3 shows an example of a marker used in a system according to present invention, having a substantially spherical form.
  • Fig. 4 is a diagram showing the analog and digital signal to be processed according to the present invention.
  • Fig. 5 is a schematic section through a digitalised marker according to the present invention.
  • Fig. 6 is a schematic dataflow diagram of a preferred camera used in an arrangement according to the present invention.
  • Fig. 7 is a schematic view of a marker image for roundness test.
  • an analogous system uses a conventional video signal from a camera as an input.
  • the X and Y coordinates for the marker which is separated from the surroundings in intensity of light are calculated.
  • the aim is to measure a movement of the marker as exact as possible, i.e. the inaccuracy, which is a result of the video signal consisting of a set of finite number of dots to be minimised.
  • the video signal consists of a number of lines, which are scanned in chronological order.
  • a marker generates an image, which extends over one or several lines.
  • the marker image on a line is called a segment.
  • the time is measured partly from the beginning of the line to the beginning of the segment QQ and partly from the beginning of the line to the end of the segment (X e ).
  • the mean value of these two periods is a measure for the position of a segment in the space, in horizontal direction (if the lines are horizontal) while the serial number of the line (S) is a measure for position of the segment in the vertical direction.
  • the length 1 of the segments is then X e -X s .
  • the ]T sign indicates that the summation is carried out over all segments being a member of the marker image.
  • the time points X s and X can be measured (when scanning the pixel value in horizontal direction) with an electronic counting device connected to an oscillator.
  • the counter starts in the beginning av the line and it is read when the start and end of a segment are reached.
  • the oscillator frequency due to technical and economical reasons is limited. In the digital case the problem may be that the image elements cannot be as small as required.
  • a comparator is provided fed with video signals, which starts an integrator, which generates a linear potential slope, which starts from a potential V a at time X s to reach V b .
  • the slope is than sampled and measured when the counter changes between two values.
  • the time when the slope passes a predetermined threshold value is calculated and used to define the time X ss .
  • the difference between X,. and X ss is may be a constant and it is determined by the integrator and the delays in the comparators.
  • the time X,.,. is easily calculated from the measured points on the potential slope at the time for the change of the counter provided that at least two points on the slope are measured. For example, if two voltage values are measured, V, at t, and V 2 at t j and V 0 is between V, and V 2 , X ss is interpolated by formula 3:
  • the time X e is measured in same way.
  • a linear slope is used, because the calculations become simpler, however, other curves may be used.
  • a simple schematic system according to the present invention is shown in fig. 1.
  • the system comprises at least one camera 10 directed to an object, in this case a human body 12, to be analysed and at least one marker 11 attached to the object 12.
  • the camera 10 may be connected to a computer device 13, to further process the received signals from the camera 10.
  • FIG. 3 An example of a substantially spherical marker 25 is shown in fig. 3, having a reflecting surface.
  • the camera 10 may be provided with means to emit IR-light, which is then reflected by the marker 25.
  • the special design of the marker allows the cameras from different angels apprehend a circular form for a correct positioning.
  • the camera 10 is equipped with a CCD unit operating substantially in same way as described above.
  • the block diagram of fig. 2 schematically illustrates some primary parts of the camera 10 intended to execute the method according to the invention.
  • the camera 10 includes optical elements 14 such as lenses and other focusing means (not shown) to project the image 16 of a marking device 11 onto a CCD unit 15.
  • optical elements 14 such as lenses and other focusing means (not shown) to project the image 16 of a marking device 11 onto a CCD unit 15.
  • the surface of the CCD unit is then scanned and the image signals including pixel information are converted to a suitable video signal by means of a converting unit 17, which may be integrated in the CCD unit.
  • the video signal representing image lines is then serially or in parallel sent to a processing unit 18.
  • the processing unit digitalises the received video signal, for example using an A/D-converter 19. This signal may also be digitalised in unit 17.
  • the processing unit may be connected to a memory unit, not shown, containing a set of instructions for controlling the processing unit.
  • the image elements may be arranged in lines, for example by means of a low-pass filter to provide some continuos signal, which can be processed as the analogous signal.
  • each image element is measured individually and from the measured values, a value is interpolated, determining when a threshold T is passed, analogously to the Basic Theory part. In this case the integrator and the comparator are eliminated.
  • Graphs representing the digital and video signals are presented in the diagram of fig. 4 using references D and A, respectively.
  • the axes of the diagram represent voltage level and pixel number in a digital or video signal.
  • the digitalized signal is then passed over to a comparison unit 20, which interpolates individual sample values about the predetermined threshold value T, also called video level, which may be obtained from a memory unit 21.
  • T also called video level
  • the object is to determine when the amplitude of the signal passes the value T.
  • Each passage presents a start and stop coordinate of each segment with a high resolution, which can be about 30 x number of pixels on a row.
  • a computation unit 22 following calculation is executed:
  • V, and V 2 are the signal levels of preceding and succeeding pixels, respectively, received from the comparatione unit 21.
  • formula 4 may be regarded as a special case of formula 3.
  • the pixel number may be obtained from a counter (not shown). Depending on the components used, levels V, and V 2 may be measured having resolution of 10 bits, the pixel number (MSB) 9 bits and (T-N,)/(N -N ⁇ ) 7 bits. Then the centre point x' of the marker is computed in a computation unit 22 by means of previous values stored in a memory unit 23, using formula (5):
  • l k is the length of the segment k (i.e., X ek -X sk ), according to fig. 5,
  • S is the serial number of the image element, and x k is the centre of the segment k.
  • the formulas (1) and (2) are substituted by formulas (5) and (6), respectively.
  • formula (1) and (2) alone do not contribute to obtaining an exact value as desired.
  • the n power of l k is calculated.
  • Fig. 5 schematically illustrates the digitalized two dimensional image 26 of the marker according to fig. 3.
  • the reason for computation of the square of the l k is to obtain a more accurate value, specially in circular images, as the lengths of the segments, i.e. ⁇ k , l k and l k+1 , shifts very fast the longer the segment is from the centre region of the circular image.
  • n does not need to be an integer.
  • the accuracy improvement in vertical direction (y') becomes spectacular for a video signal, as the number of line segments in a marker is few.
  • the l k 's may then be stored in the memory unit 23, for further calculations.
  • A i ⁇ . ⁇ , which yields formula (7):
  • results are then transformed to an interface unit to further transmission to the computer unit 13, in which the computed values, x' and r, can be used to show the position of the marker on a screen for simulation, positioning and other applications.
  • the dataflow in a preferred embodiment of a camera 10 is illustrated in fig. 6, in which parts functioning as above are designated with same designation numeral.
  • the camera 10 is provided withe s light source 27, such as an IR-flash, surrounding the lens opening on the camera casing.
  • the CCD sensor 15, is arranged in communication with the lens 14.
  • the signal from CCD is provided to a video processing unit 28, which after processing the video signal from the CCD sensor 15 forwards the processed signal to an analog/digital converter 19.
  • the CCD sensor 15, video processing unit 28 and A/D-converter are synchronised by a common timing signal outputted from timing control device 29, e.g. including an oscillator and a counter.
  • the result from the A/D converter is provided to a memory unit 30, preferably of FIFO memory (First In, First Out).
  • a register 31 is arranged for supplying a preceding data.
  • a comparator 20 generating a write signal, preferably two pixel values per translations are stored in the memory unit 30.
  • the number of the pixel and line number of image line is provided by the counter of the timing control unit 29.
  • the data stored in the memory unit 30 is then used by a processing unit 32, for instance for centroid calculation as described above.
  • an instruction set memory unit 33 is arranged for storing the operation instructions
  • a data memory unit 34 for storing processed data
  • a communication unit 35 is arranged for transferring the processed data, for example to a computer.
  • the system must carry out a "roundness" test, so that only images that pass the test, may be analysed.
  • the roundness test according to the invention may be employed in any other application requiring similar test or procedure.
  • a first step the symmetry of the digital image of the marker is tested to see if a rectangle, which exactly enclose all segments is a square, i.e.
  • the processing unit 18 may be integrated in the camera 10, or in a peripheral unit. The number of units and their function in the processing unit 18 may also vary.
  • the system is not limited to a camera using CCD sensors.
  • Video signals from a camera can be used in a scanning apparatus to track and determine the position of a marker for centroid calculations.
  • the form of the marker may vary due to the application area.
  • a trapezoid formed marker may be used.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Image Analysis (AREA)
EP98901151A 1997-01-13 1998-01-13 Method and arrangement for finding the centroid of a body Withdrawn EP1008105A1 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
SE9700065 1997-01-13
SE9700065A SE512426C2 (sv) 1997-01-13 1997-01-13 Förfarande och anordning för att bestämma centrumpunktskoordinaten för ett objekt
PCT/SE1998/000032 WO1998030976A1 (en) 1997-01-13 1998-01-13 Method and arrangement for finding the centroid of a body

Publications (1)

Publication Number Publication Date
EP1008105A1 true EP1008105A1 (en) 2000-06-14

Family

ID=20405397

Family Applications (1)

Application Number Title Priority Date Filing Date
EP98901151A Withdrawn EP1008105A1 (en) 1997-01-13 1998-01-13 Method and arrangement for finding the centroid of a body

Country Status (5)

Country Link
EP (1) EP1008105A1 (ja)
JP (1) JP2001508209A (ja)
AU (1) AU5685098A (ja)
SE (1) SE512426C2 (ja)
WO (1) WO1998030976A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000253293A (ja) * 1999-02-26 2000-09-14 Minolta Co Ltd デジタルカメラ
WO2010073174A1 (en) * 2008-12-25 2010-07-01 Koninklijke Philips Electronics N.V. System and method for image capturing

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5072294A (en) * 1989-06-07 1991-12-10 Loredan Biomedical, Inc. Method and apparatus for analyzing a body having a marker located thereon
JP2854359B2 (ja) * 1990-01-24 1999-02-03 富士通株式会社 画像処理システム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO9830976A1 *

Also Published As

Publication number Publication date
SE9700065L (sv) 1998-07-14
WO1998030976A1 (en) 1998-07-16
JP2001508209A (ja) 2001-06-19
SE512426C2 (sv) 2000-03-13
AU5685098A (en) 1998-08-03
SE9700065D0 (sv) 1997-01-13

Similar Documents

Publication Publication Date Title
EP0951696B1 (en) Method and arrangement for determining the position of an object
US6549288B1 (en) Structured-light, triangulation-based three-dimensional digitizer
CN104837436B (zh) 三维扫描中的运动补偿
EP1596158B1 (en) Three-dimensional shape input device
US6256099B1 (en) Methods and system for measuring three dimensional spatial coordinates and for external camera calibration necessary for that measurement
US7502018B2 (en) Projector, electronic whiteboard system using projector and method of acquiring indicated point
US20040246229A1 (en) Information display system, information processing apparatus, pointing apparatus, and pointer cursor display method in information display system
EP1192414A1 (en) Method and system for measuring the relief of an object
US4628469A (en) Method and apparatus for locating center of reference pulse in a measurement system
US20020118274A1 (en) Three-dimensional measuring method and system
WO2018168757A1 (ja) 画像処理装置、システム、画像処理方法、物品の製造方法、プログラム
JP3493403B2 (ja) 3次元計測装置
CN105717502B (zh) 一种基于线阵ccd的高速激光测距装置
WO1998030976A1 (en) Method and arrangement for finding the centroid of a body
Gruen et al. Real-time photogrammetry at the digital photogrammetric station (DIPS) of ETH Zurich
JP2681745B2 (ja) レーザ光を利用したスペックルパターンによる被計測物の上下および横移動量の測定方法。
JPH10512953A (ja) 光学的なレンジ及び速度検出システム
JPH076777B2 (ja) パターンの輪郭検出方法及びこの方法を用いた測長装置
Shortis et al. Calibration of CCD cameras for field and frame capture modes
JPS60194302A (ja) 被写体計測装置
Greaves State of the Art in Automated Motion Tracking and Analysis Systems
JP2001033245A (ja) 平面上の点の位置測定方法
JPH10215467A (ja) 位置測定方法及び位置測定システム
Stone et al. A CCD system for reading circles in wide-angle astrometry
Baca Spatial reconstruction of marker trajectories from high-speed video image sequences

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 19990809

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE CH DE DK ES FI FR GB GR IE IT LI LU MC NL PT SE

GRAH Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOS IGRA

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20030614