WO2004056262A1 - Apparatus and method for the reproduction of measuring conditions based on the body surface or surface veins and the contacting pressure - Google Patents

Apparatus and method for the reproduction of measuring conditions based on the body surface or surface veins and the contacting pressure Download PDF

Info

Publication number
WO2004056262A1
WO2004056262A1 PCT/CN2003/000813 CN0300813W WO2004056262A1 WO 2004056262 A1 WO2004056262 A1 WO 2004056262A1 CN 0300813 W CN0300813 W CN 0300813W WO 2004056262 A1 WO2004056262 A1 WO 2004056262A1
Authority
WO
WIPO (PCT)
Prior art keywords
measurement
image
axis
contact pressure
human body
Prior art date
Application number
PCT/CN2003/000813
Other languages
French (fr)
Chinese (zh)
Inventor
Kexin Xu
Xiaodong Hu
Jingying Jiang
Original Assignee
Tianjin Sunshine Optics Technologies Co., Ltd.
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 Tianjin Sunshine Optics Technologies Co., Ltd. filed Critical Tianjin Sunshine Optics Technologies Co., Ltd.
Priority to AU2003272847A priority Critical patent/AU2003272847A1/en
Publication of WO2004056262A1 publication Critical patent/WO2004056262A1/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/27Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/14Coupling media or elements to improve sensor contact with skin or tissue
    • A61B2562/146Coupling media or elements to improve sensor contact with skin or tissue for optical coupling

Definitions

  • the present invention relates to a measurement device and method, and in particular to a device and method for reproducing measurement conditions based on human body surface or surface texture features and contact pressure. Background technique
  • the currently available methods are: (1) measuring multiple groups of infrared absorption spectra for a specific part of a specific individual; (2) measuring the component content of each group of spectra at a corresponding time using an invasive method; (3) using the above steps to obtain Establish a mathematical model of the infrared absorption spectrum and component content data; (4) Use this mathematical model to analyze the component content corresponding to the infrared absorption spectrum obtained from subsequent actual measurements.
  • certain restrictions must be placed on the measurement locations and conditions, that is, the actual measurement conditions of the subsequent component content are consistent with the measurement conditions established by the mathematical model.
  • the current implementation methods can be roughly divided into two categories: One is the use of specific limit devices, such as: hand molds, finger cots, limit stops, etc., so that the measurement location of the measured object can basically be It is limited to a certain area. The second is to make specific marks at or near the measurement site, use image acquisition and processing methods to calculate the relative position, and move the measurement device to the corresponding position.
  • specific limit devices such as: hand molds, finger cots, limit stops, etc.
  • the present invention establishes a measurement condition reproduction device based on human body surface layer or surface texture features and contact pressure, and proposes a corresponding implementation method.
  • the surface texture features are the surface texture of the skin such as fingerprints and palm prints.
  • the surface layer refers to the texture features reflected by blood vessels, bones and pores near the surface of the skin.
  • the measurement position will be the saved image as a template image, and the contact pressure information will be attached.
  • the human body measurement part After the human body measurement part is repositioned, first collect the surface layer or surface texture image of the measurement part, and then match and identify the current image and the template image by image processing to obtain the displacement of the measurement part in the xy plane and Offset of the z-axis rotation angle.
  • the image is matched with the template image using the following method: (1)
  • the wavelet transform and exhaustion method are used to determine the rotation angle.
  • the wavelet transform reduces the data amount of the texture image at the measurement site and accelerates the process of rotation angle correction;
  • a specific spectrum measurement device as shown in Figure 1: It consists of fiber optic probes in the x, y, and z directions. Linear displacement and z-axis angle servo positioning device, stage, lighting device, camera, frame grabber, spectrum detection system and computer processing system; contact pressure sensor and temperature sensor are set at the measurement site and the probe, and the pressure, Temperature signal processing and acquisition device.
  • the work of the device includes the following steps:
  • step (2) Perform spectral measurement under the conditions of step (1), and create a specific mathematical model of spectral measurement for the measurement location based on the spectral data;
  • the device makes the optical fiber probe reach the corresponding position, and then adjusts the position of the measuring element in the z-axis direction according to the contact pressure between the optical fiber probe and the measurement part, and measures the temperature of the contact part, as shown in FIG. 3;
  • step (2) Use optical fiber probe and spectral detection system to perform spectral detection, and measure the measurement contact at that time Perform temperature analysis on the temperature of the part, the measured spectral data and the mathematical model of spectral measurement established in step (2);
  • the surface texture or surface texture of the human body can be the texture of the skin surface, such as fingerprints and palm prints, or the texture features of blood vessels, bones and pores near the surface of the skin, which are fixed under certain conditions;
  • the probe is repeatedly positioned in the vertical z-axis
  • the wavelet transform is used to reduce the data amount of the texture image at the measurement site and accelerate the process of rotation angle correction
  • the servo system is used to realize the automatic positioning of the optical fiber probe.
  • Figure 1 Block diagram of measurement condition reproduction system
  • Figure 4 Structure of an image positioning system based on palmprint recognition
  • Figure 5 A schematic diagram of the structure of an image positioning system based on fingerprint recognition
  • Figure 6a a front view of the structure of a ring-shaped edge lighting device
  • Figure 6b Side view of the structure of the ring-shaped edge lighting device
  • Figure 7 Schematic diagram of position shift between two images
  • FIG. 8 Image matching processing flow
  • Figure 10 The palmprint image is currently collected
  • Figure 11 Phase correlation results of the currently acquired image after rotation correction processing with the template image
  • Figure 12 Superposition results of the current acquisition image after rotation and translation correction with the template image
  • Figure 13a Fiber optic probe in contact with the measurement site Schematic diagram of pressure detection
  • Figure 13b Top view of the contact pressure detection between the optical fiber probe and the measurement site
  • Figure 14 Schematic diagram of the contact temperature between the optical fiber probe and the measurement site
  • Figure 16 Comparison of the CV values of the spectral measurement results before and after the system is reproduced using the measurement conditions; Figure 16- (1) Test object 1;
  • a lighting device is installed below the stage with a hole, the device is provided with a hole attached to the stage, and the axis coincides;
  • the servo positioning device is composed of four parts: an X-axis servo positioning device, a y-axis servo positioning device, z-axis servo positioning device, z-axis rotation servo positioning device; z-axis rotation servo positioning device is installed on z-axis servo positioning device; optical fiber probe is installed on Z- axis rotation servo positioning device, and the center axis coincides with its rotation axis; servo
  • the displacement mechanism in the positioning device is driven by a stepper motor and controlled by a servo controller, which is controlled by a computer.
  • the texture features of the human body's surface layer or surface are relatively obvious, and the location is convenient for non-invasive infrared spectrum detection. According to the characteristics of the measurement parts on the fingers and palms, the structure of the measurement condition reproduction system is also different.
  • the optical fiber probe 8 and the CCD The camera 1 is located at the same location 2 as the inside of the palm, that is, the image matching information is the texture of the palm area, and the measurement position of the infrared spectrum is the palm of the hand.
  • the stage 3 is driven by the X-axis servo positioning device 4 and the y-axis servo positioning device 5.
  • the X-axis servo positioning device is fixed on the base plate 7.
  • the y-axis servo positioning device is orthogonally fixed to the X-axis servo positioning device on the X-axis servo positioning device.
  • the stage is fixed on the y-axis servo positioning device.
  • the servo positioning device for driving the optical fiber probe 8 is composed of two parts step by step: a z-axis servo positioning device 6 and a z-axis rotary servo positioning device 9.
  • the z-axis servo positioning device is fixed on the base plate by a vertical connection.
  • the z-axis servo positioning device is installed on the z-axis servo positioning device.
  • the optical fiber probe is mounted on the Z- axis rotary servo positioning device, and the center axis coincides with its rotation axis.
  • the displacement mechanism in the servo positioning device is driven by a stepping motor and controlled by a servo controller, which is controlled by a computer. .
  • the optical fiber probe 8 and the CCD camera 1 are located at the measured part 2 such as both sides of a finger, that is, the information of image matching is the fingerprint, and the measurement position of the infrared spectrum is The back of the finger, such as the nail area, is fixed on the stage 3.
  • the servo positioning device driving the optical fiber probe is composed of four parts step by step: X-axis servo positioning device 4, y-axis servo positioning device 5, Z- axis servo positioning device 6, z-axis rotary servo positioning device 9.
  • the X-axis servo positioning device is connected to the bottom plate 7.
  • the y-axis servo positioning device is orthogonally fixed to the X-axis servo positioning device on the X-axis servo positioning device.
  • the z-axis servo positioning device is mounted on the y-axis servo positioning device through a vertical connection, and the Z- axis rotates.
  • the servo positioning device is mounted on the z-axis servo positioning device.
  • the optical fiber probe 8 is mounted on the z-axis rotation servo positioning device 9, and the center axis thereof coincides with its rotation axis.
  • the displacement mechanism in the servo positioning device is driven by a stepping motor and controlled by a servo controller, which is controlled by a computer.
  • the palm and fingers are placed on a carrier with the front side facing down, and a through hole 10 is provided on the carrier to facilitate the collection of palm prints or fingerprint images.
  • a ring-shaped edge lighting device 11 is installed under the carrier. The device also has a through hole, the same size as the hole of the carrier, and the axes coincide.
  • Illumination of the measurement site is achieved by means of circular edge illumination.
  • a light emitting diode (LED) light source fixing structure is fixed below the palm and finger bearing surface, as shown in Figs. 6a and 6b.
  • the material of the fixing structure is plexiglass, and the circumference of the circular hole 10 is 24 grooves 12 are evenly arranged, and each groove can be embedded with a light emitting diode (LED).
  • the bottom surface of the groove 12 is an inclined surface, which is higher near the center to ensure that the LED is illuminated obliquely upward.
  • the plexiglass on the side of the center hole is roughened, so that the diffuse reflection of the LED illuminates the measurement site.
  • the CCD camera 1 is installed directly below the center axis of the circular hole, as shown in Figs. 4 and 5.
  • the CCD camera is fixedly installed on the bottom plate directly below the fingerprint.
  • the optical fiber probe In the palmprint positioning system, it is required that the optical fiber probe should be directly under the palm in the test state, which conflicts with the front image of the palmprint.
  • the optical fiber probe is mounted on the z-axis rotary servo positioning device, and the CCD camera is mounted on the z- axis servo positioning device.
  • the Z axis servo positioning device can be translated to different positions under the drive of the X axis servo positioning device and the y axis servo positioning device, that is, the CCD camera is moved to the stage during the image acquisition stage. Just below the circular hole, move the fiber-optic probe directly under the circular hole of the stage in the spectrum measurement stage.
  • the CCD camera outputs a standard video signal and inputs it to a frame grabber in a computer system.
  • Figure 7 is a schematic diagram of the offsets of the two images.
  • the origins of the two coordinate axes in the figure are the centers of the two images. It can be seen that there is a displacement offset and a rotational offset of the xy plane in the two images.
  • the image capture card collects palm print and fingerprint images, it needs to enter the image processing process to determine the position offset of the two images, and then drive the servo positioning system with the position offset to make the optical fiber probe reach the corresponding position.
  • the image processing flow is divided into three stages: (1) the determination of the rotation angle using wavelet transform and exhaustion method; (2) the use of phase correlation to detect the displacement change of pixel accuracy; (3) the use of quadratic surface fitting for sub- Pixel-precision displacement change detection.
  • FIG. 8 is an image processing flowchart, which is described in detail below:
  • the resolution of the set image is
  • wavelet transform is performed on the currently acquired image and the template image before determining the rotation angle of the exhaustion method, reducing the image resolution to 32x32, which greatly reduces the time of this image processing stage.
  • the wavelet multi-scale of the image represents a pyramid-like structure, which is very beneficial for fast image matching from coarse to fine.
  • the overall characteristics of the image are basically maintained in the scale space, that is, the low-frequency portion, and the amount of data in the scale space portion after extraction is significantly reduced compared to the original image.
  • Phase correlation is an image matching algorithm with little geometric distortion. Law. Because the geometric distortion has a large impact on the high-frequency components of the image and a small impact on the low-frequency components, a phase correlation algorithm based on the Fourier spectrum with a low-pass filter characteristic is used in this device, which can greatly reduce the geometric distortion on the matching performance Impact.
  • ⁇ and are the results of Fourier transform of two images (currently acquired image and template image), respectively.
  • the phase spectrum contains the position shift information of the two images, and it is a power spectrum whose spectral amplitude is 1 in the full frequency domain.
  • the phase correlation function is a delta pulse function located at the position offset (x0, yO) of the two figures, which is also called the correlation peak.
  • the value is 1, otherwise it is 0. Therefore, the phase correlation calculation results of the two images are used in this device to determine the displacement amount of the images.
  • phase matching is used for image matching to determine the displacement offset with the following characteristics: (1) Large displacement detection range. When there is only a displacement change between the two images, the maximum displacement offset of the xy plane can be detected to reach half the image width. (2) Sharp related peaks. When the two images are completely correlated, the calculated ⁇ pulse function has a very sharp correlation peak, which can achieve an accurate matching of the positions of the center points of the two images. (3) Little dependence on image gray. When there is a gray level difference between the two images, the position of the calculated ⁇ pulse function does not change, but there is a difference in amplitude, so it has a strong ability to resist image occlusion. (4) The amount of rotation offset between the two images has a greater impact on the results of image matching.
  • the rotation angle offset is greater than 5 °, even if the two images are identical, the correlation peak attenuation is zero, that is, the displacement offset of the two images in the xy plane cannot be determined based on the correlation peak. Therefore, in the previous step, the rotation angle is determined by judging the size of the correlation peak, and the image is rotated by this angle to perform phase correlation, and finally the position offset of the two images is obtained.
  • the position offset of the currently acquired image and the template image can be determined, and the minimum unit of the offset is 1 pixel.
  • the following sub-pixel matching algorithm based on quadric surface fitting is used.
  • a surface fitting method is adopted. This method has the advantages of fast and easy calculation, and the accuracy can reach sub-pixels.
  • the idea of the surface fitting method is: take the best matching point at the pixel level as the center, perform surface fitting according to the similarity measure, and then calculate the exact position of the extreme point by corresponding mathematical methods.
  • This device uses the correlation coefficient of phase correlation as the similarity measurement feature, selects the quadric surface as the fitting function, and uses the multivariate least squares regression method in the calculation to determine the exact location of the extreme point.
  • the quadric fitting function uses the formula:
  • PC (x, y) ax + by + cxy + dx + ey + f
  • PC (x, y) is a phase correlation value corresponding to the position (x, y).
  • the above function can be written as:
  • the device adopts the multivariate least squares regression method in the fitting calculation, which makes the calculation simple and accurate.
  • the vector X is used as the regression coefficient, and it is assumed that the value of the random variable B depends on the independent variable in the matrix A.
  • the regression coefficient is obtained as the coefficient of the fitting function. After the coefficients of the fitting function are obtained, the precise position of the image offset with sub-pixel accuracy can be obtained using the following formula.
  • FIG. 9 as a template palmprint image 10 is currently acquired palmprint, the resolution of the image is two 256 X 256. It can be seen that there is a significant angular and displacement shift in the two images.
  • Increase angle The amount is set to 0.1 degree, and the rotation angle is detected by applying wavelet transform and exhaustion method, and the rotation angle of the palm image of the current set relative to the template palm image is 13.9 degrees counterclockwise. After the currently collected palm print image is rotated counterclockwise by 13.9 degrees, phase correlation matching processing is performed with the template palm print image.
  • Figure 11 is a flat distribution map of the correlation function obtained by processing, and the horizontal offset of the two images is obtained according to the position of the maximum correlation peak.
  • the amount is -20 pixels, and the vertical offset is 64 pixels.
  • FIG. 12 it can be seen from the superimposed image that the above-mentioned detection of the rotation angle and displacement offset using phase correlation is correct and effective.
  • the horizontal offset of the two images is -20.08 pixels, vertical
  • the offset in the direction is 64.20 pixels.
  • the rotation detection accuracy of the measurement condition reproduction device can reach 0.1 degrees, and the translation detection accuracy can reach 1/50 of a pixel.
  • the relative position of the measurement site in the xy plane that is, the offset and rotation angle of the xy axis
  • the xy plane servo positioning device to drive the optical fiber probe so that the probe axis coincides with the center point of the measurement site.
  • the optical fiber probe is rotated by a certain angle through the z-axis rotation positioning device, which is the rotation angle obtained by the image processing process.
  • the optical fiber probe is driven by the z-direction approaching servo device to bring the optical fiber probe into contact with the measurement site.
  • the pressure sensor After the optical fiber probe contacts the measurement part, the pressure sensor will output the contact pressure signal and input it into the signal acquisition card.
  • the computer compares the collected contact pressure signal with a preset value, and then adjusts the relative position of the optical fiber probe in the z direction. Position, so that the contact pressure between the optical fiber probe and the measurement site is basically maintained near a preset pressure value.
  • the measurement of contact pressure is shown in Figure 13.
  • a circular-shaped gas path bracket 13 is fixed on the top of the cylindrical optical fiber probe, and a silicone tube 14 is installed in the annular area between the bracket and the optical fiber probe 8.
  • the inside of the silicone tube is filled with air and connected with a
  • the silicon pressure sensor is connected to form a closed gas path.
  • the computer can The card collects the contact temperature of the two, and uses the data as parameters for mathematical model establishment and infrared spectrum analysis.
  • the temperature sensor 15 is a thin-film thermistor, and is attached to the measurement surface of the optical fiber probe, as shown in FIG. 14. There are two bundles of optical fibers installed in the optical fiber probe, one bundle is distributed in the center area of the probe, and the other bundle is distributed in a ring area of the probe. Four temperature sensors are evenly pasted in the middle area between the two optical fibers, and ten temperature sensors are evenly pasted along the circumference in the peripheral area of the optical fiber bundle 16.
  • the signal acquisition card collects the temperature of the 14 detection points into the computer system and performs average processing to ensure the accuracy of the temperature measurement.
  • the infrared optical path detection system starts to work.
  • the relative position of the measurement site is required to be kept basically fixed.
  • the CCD camera collects an image of the contact area of the optical fiber probe and the measurement site-the position of the round hole on the stage, and uses it as a motion state to determine
  • the computer collects an image of the same area after every five wavelengths of spectral measurement, and performs a subtraction operation with the template image discriminated by the motion state. Whether the measurement site has moved during the process.
  • Fig. 15 is a processing flow for judging a motion state.
  • the spectrum measurement can be entered.
  • Measurement conditions The measurement of the spectral data when the reproduction device is not working: First, adjust the center axis of the optical fiber probe and the center axis of the round hole of the stage to coincide, and then manually adjust the optical fiber probe to be vertical The orientation position makes the optical fiber probe have a certain contact pressure with the palm. At the same time, the z-axis rotary servo positioning mechanism is reset to the zero position.
  • the palm of the measured object is lifted after each spectrum measurement is completed, and the servo is reset when it is placed on the stage.
  • the positioning device does not work, that is, the position of the optical fiber probe is fixed; the palm is repeatedly placed on the carrier 10 times, each time 10 sets of spectral data were obtained.
  • Measurement conditions Measurement of the spectral data during the operation of the reproduction device: Determine and save the measurement conditions before the measurement of the first group of spectral data, such as: template palm print image, contact pressure, and the palm of the measured object remains after each spectral measurement is completed.
  • the evaluation index of the statistical graph is the ratio of the maximum fluctuation of the light intensity corresponding to each wavelength point in the 10 sets of spectral data to the light intensity, which is abbreviated as CV.
  • the dashed line in the figure indicates the statistical result of the spectral data when the measurement condition reproduction device is not working, and the solid line in the figure indicates the statistical result of the spectral data when the measurement condition reproduction device is in operation.
  • the change in the position of the measurement site has a greater impact on the spectral measurement, and corresponding measures must be adopted to limit the change in position;
  • the present invention is based on The human body surface layer or surface texture feature and contact pressure measurement condition reproduction device can greatly reduce the influence of the change of the position of the measurement location on the spectrum measurement, and improve the measurement stability and repeatability of the spectrum.

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Biomedical Technology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Biophysics (AREA)
  • Chemical & Material Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

The present invention relates to an apparatus and a method for reproduction of measuring conditions based on the body surface or surface veins and the contacting pressure. The apparatus comprises an optical fiber measuring head, servo positioning devices for straight displacements in the x, y, z directions and angle around the z-axis, an object stage, a lighting device, a camera, an image collecting card, a spectrum detecting system and a computer processing system; the apparatus is simultaneously mounted with a contacting pressure sensor, a temperature sensor and a processing and collecting device for pressure and temperature signals. The steps of the method are as follows: (1) collecting, matching and positioning the image of body surface or surface veins; (2) performing the positioning for measuring position of the xy plane displacement and the rotating angle around z-axis and the like according to the result of the image matching; (3) performing the positioning in z direction according to the contacting pressure of the position of the measuring element and the body surface or surface veins; (4) during the operation of the measuring element, collecting the image of the body surface or surface veins characters, monitoring and measuring the movement state of the measured position.

Description

基于人体表层或表面纹理特征和接触压力的  Based on human body surface or surface texture features and contact pressure
测量条件复现装置与方法  Measuring condition reproduction device and method
技术领域 Technical field
本发明涉及一种测量装置和方法, 特别涉及一种基于人体表层或表 面紋理特征和接触压力的测量条件复现装置与方法。 背景技术  The present invention relates to a measurement device and method, and in particular to a device and method for reproducing measurement conditions based on human body surface or surface texture features and contact pressure. Background technique
在人体内成分 (如: 血糖浓度等)的无创的红外光谱检测过程中, 由 于人体组织结构方面如皮肤的组织与厚度、 血管的分布、 脂肪的成份与 厚度等的差异和外界环境的波动, 要求测量条件如测量的部位、 接触压 力等保持恒定。  During the non-invasive infrared spectrum detection of human body components (such as blood glucose concentration, etc.), due to differences in human tissue structure such as skin tissue and thickness, blood vessel distribution, fat composition and thickness, and fluctuations in the external environment, It is required to keep the measurement conditions such as the measurement site and the contact pressure constant.
在基于红外光谱的检测系统中, 即使人体内某待测成分的含量基本 恒定, 不同人体或同一人体不同部位的光谱测量数据也有可能呈现出较 大的差异。 因此目前可行的方法是: (1)针对特定个体的特定部位测量出 多组红外吸收光谱; (2)用有创的方法测量出每组光谱对应时刻的成分含 量; (3)利用以上步骤得到的红外吸收光谱和成分含量的数据建立数学模 型; (4)利用该数学模型分析出后续实际测量得到的红外吸收光谱所对应 的成分含量。 为了保证该数据模型的可重用性及最终测量结果的可靠性, 必须对测量部位和条件进行一定的限制, 即后续成分含量实际测量的条 件与数学模型建立的测量条件保持一致。  In a detection system based on infrared spectroscopy, even if the content of a component to be measured in the human body is basically constant, the spectral measurement data of different human bodies or different parts of the same human body may show a large difference. Therefore, the currently available methods are: (1) measuring multiple groups of infrared absorption spectra for a specific part of a specific individual; (2) measuring the component content of each group of spectra at a corresponding time using an invasive method; (3) using the above steps to obtain Establish a mathematical model of the infrared absorption spectrum and component content data; (4) Use this mathematical model to analyze the component content corresponding to the infrared absorption spectrum obtained from subsequent actual measurements. In order to ensure the reusability of the data model and the reliability of the final measurement results, certain restrictions must be placed on the measurement locations and conditions, that is, the actual measurement conditions of the subsequent component content are consistent with the measurement conditions established by the mathematical model.
为了保持测量条件恒定, 目前的实现方法可大致分为两大类: 一是 采用特定的限位装置, 如: 手模、 指套、 限位挡块等, 使得被测对象的 测量部位基本能够限定在某个区域; 二是在测量部位或附近区域作特定 的标记, 利用图像采集和处理的方法, 计算相对位置, 并移动测量装置 到相应位置。 现有技术存在以下问题:  In order to keep the measurement conditions constant, the current implementation methods can be roughly divided into two categories: One is the use of specific limit devices, such as: hand molds, finger cots, limit stops, etc., so that the measurement location of the measured object can basically be It is limited to a certain area. The second is to make specific marks at or near the measurement site, use image acquisition and processing methods to calculate the relative position, and move the measurement device to the corresponding position. The following problems exist in the prior art:
' (1)不能自动、 有效地识别被测对象是否为同一个体, 必须人为干预; (2)测量位置的复现精度较低, 最高只能到像素级; (4)对测量部位的旋'转现象缺乏有效的解决方案; '(1) It is not possible to automatically and effectively identify whether the measured object is the same individual, and human intervention is required; (2) The reproduction accuracy of the measurement position is low, and it can only reach the pixel level; (4) There is no effective solution to the rotation phenomenon of the measurement site;
(5)对测量过程中发生的测量部位移动所造成的测量误差未能进行考虑 和解决。 . 发明内容  (5) The measurement error caused by the movement of the measurement site during the measurement process cannot be considered and resolved. Summary of the Invention
为了克服上述现有技术中的不足, 本发明建立基于人体表层或表面 紋理特征和接触压力的测量条件复现装置, 并提出相应的实现方法。  In order to overcome the above-mentioned shortcomings in the prior art, the present invention establishes a measurement condition reproduction device based on human body surface layer or surface texture features and contact pressure, and proposes a corresponding implementation method.
本发明的方法步骤如下:  The method steps of the present invention are as follows:
(1)人体表层或表面紋理特征的图像拾取、 匹配和定位;  (1) Image pickup, matching and positioning of human body surface or surface texture features;
(2)根据表层或表面紋理图像匹配的结果进行 xy平面位移和 z轴旋转 角度等测量位置的定位;  (2) Positioning of measurement positions such as xy plane displacement and z-axis rotation angle according to the results of surface or surface texture image matching;
(3)依据测量元件与人体表层或表面纹理所在部位的接触压力进行 z方 向的定位;  (3) Positioning in the z direction according to the contact pressure between the measuring element and the surface layer or surface of the human body;
(4)在测量元件工作的过程中对人体表层或表面纹理特征继续采集图 像, 监测测量部位的运动状态。  (4) During the operation of the measuring element, continue to collect images of human body surface or surface texture features, and monitor the movement status of the measurement part.
其中: among them:
对人体测量部位的表层或表面纹理进行图像采集, 表面纹理特征 是皮肤表面紋路如: 指纹和掌纹等, 表层是指皮肤表层附近血管、 骨骼 及毛孔等体现出的纹理特征, 依据人为预设的测量位置将保存图像作为 模板图像, 并附加接触压力信息。  Image collection of the surface layer or surface texture of the human body measurement surface. The surface texture features are the surface texture of the skin such as fingerprints and palm prints. The surface layer refers to the texture features reflected by blood vessels, bones and pores near the surface of the skin. The measurement position will be the saved image as a template image, and the contact pressure information will be attached.
人体测量部位重新放置后, 首先采集测量部位的表层或表面纹理图 像, 然后通过图像处理的方法对当前图像与模板图像进行匹配识别, 得 到当前测量部位与模板图像建立时测量部位在 xy平面位移和 z轴旋转角 度的偏移量。  After the human body measurement part is repositioned, first collect the surface layer or surface texture image of the measurement part, and then match and identify the current image and the template image by image processing to obtain the displacement of the measurement part in the xy plane and Offset of the z-axis rotation angle.
依据上述得到的 xy平面位移和 z轴旋转角度的偏移量驱动伺服装置 使测量元件到达相应位置, 然后依据测量元件与测量部位的接触压力调 整测量元件在 z轴方向的位置;  Drive the servo device according to the XY plane displacement and the z-axis rotation angle deviation obtained above to make the measuring element reach the corresponding position, and then adjust the position of the measuring element in the z-axis direction according to the contact pressure between the measuring element and the measurement part;
启动测量元件进入工作状态, 继续采集测量部位的表层或表面紋理 图像, 通过图像相减运算得到测量部位的运动状态。  Start the measurement element into the working state, continue to collect the surface layer or surface texture image of the measurement part, and obtain the movement state of the measurement part by image subtraction operation.
图像与模板图像进行匹配是采用下述的方法: (1)应用小波变换及穷尽法进行旋转角度的确定, 小波变换减少测量部 位纹理图像的数据量, 加速旋转角度校正的进程; The image is matched with the template image using the following method: (1) The wavelet transform and exhaustion method are used to determine the rotation angle. The wavelet transform reduces the data amount of the texture image at the measurement site and accelerates the process of rotation angle correction;
(2)利用相位相关进行像元精度的位移变化检测, 减小图像失真和光 不均匀对图像匹配的影响;  (2) Use phase correlation to detect displacement changes of pixel accuracy to reduce the impact of image distortion and light unevenness on image matching;
(3)利用曲面拟合技术进行亚像元精度的位移变化检测, 使得图像匹配 的精度达到亚像元级。  (3) Use surface fitting technology to detect displacement changes of sub-pixel accuracy, so that the accuracy of image matching reaches sub-pixel level.
根据本发明的一种基于人体表层或表面纹理特征和接触压力的测量 条件复现方法我们设计了具体的光谱测量装置: 如图 1 所示: 它由光纤 测头, x、 y、 z 方向的直线位移和 z 轴角度伺服定位装置, 载台, 照明 装置, 摄像机, 图像采集卡, 光谱检测系统和计算机处理系统组成; 在 测量部位与测头设置有接触压力传感器和温度传感器, 和对压力、 温度 信号的处理及采集装置。  According to a method for reproducing the measurement conditions based on human body surface or surface texture features and contact pressure, we have designed a specific spectrum measurement device: as shown in Figure 1: It consists of fiber optic probes in the x, y, and z directions. Linear displacement and z-axis angle servo positioning device, stage, lighting device, camera, frame grabber, spectrum detection system and computer processing system; contact pressure sensor and temperature sensor are set at the measurement site and the probe, and the pressure, Temperature signal processing and acquisition device.
该装置的工作包括下列步骤:  The work of the device includes the following steps:
(1)建立模板模型, 其步骤如图 2所示: 驱动摄像机到初始位置, 采集 人体表层或表面的纹理信息, 将图像保存到模板图像库, 驱动光纤测头 到图像中心, 驱动光纤测头与测量部位接触, 并保持设定的接触压力, 测量接触部位的温度;  (1) Establish a template model, the steps are shown in Figure 2: Drive the camera to the initial position, collect the texture information of the human body surface or surface, save the image to the template image library, drive the optical fiber probe to the image center, and drive the optical fiber probe Make contact with the measurement site, maintain the set contact pressure, and measure the temperature of the contact site;
(2)在步骤 (1)的条件下进行光谱测量, 并依据光谱数据创建该测量部 位特定的光谱测量数学模型;  (2) Perform spectral measurement under the conditions of step (1), and create a specific mathematical model of spectral measurement for the measurement location based on the spectral data;
(3)人体测量部位重新放置后, 首先驱动摄像机到初始位置, 采集人体 测量部位的表层或表面纹理图像, 然后判断模板图像是否存在, 如果不 存在显示错误信息并结束, 否则选择相应的模板图像, 然后通过图像处 理的方法对当前图像与模板图像进行匹配识别, 得到当前测量 ¾位与模 板图像建立时测量部位在 xy平面位移和 z轴旋转角度的偏移量及最大相 关峰的幅值, 如果最大相关峰的幅值超过限定值, 提示模板图像选择错 误并结束, 否则继续判断 xy平面的偏移量是否超出驱动机构的位移量, 如果超限, 提示超限信息后结束, 否则驱动伺服装置使光纤测头到达相 应位置, 然后依据光纤测头与测量部位的接触压力调整测量元件在 z 轴 方向的位置, 并测量接触部位的温度, 如图 3所示;  (3) After the human body measurement part is repositioned, first drive the camera to the initial position, collect the surface or surface texture image of the human body measurement part, and then judge whether the template image exists. If there is no error message displayed and end, otherwise select the corresponding template image Then, the current image and the template image are matched and identified through image processing to obtain the current measurement position and the offset of the measurement position in the xy plane and the z-axis rotation angle and the amplitude of the maximum correlation peak when the template image is established If the amplitude of the maximum correlation peak exceeds the limit, it prompts the template image to be selected incorrectly and ends. Otherwise, it continues to determine whether the offset of the xy plane exceeds the displacement of the drive mechanism. If it exceeds the limit, it ends after prompting the limit information, otherwise it drives the servo. The device makes the optical fiber probe reach the corresponding position, and then adjusts the position of the measuring element in the z-axis direction according to the contact pressure between the optical fiber probe and the measurement part, and measures the temperature of the contact part, as shown in FIG. 3;
(4)利用光纤测头和光谱检测系统进行光谱检测, 测量出当时测量接触 部位的温度, 与测量得到的光谱数据与步骤 (2)建立的光谱测量数学模型 进行数据分析; (4) Use optical fiber probe and spectral detection system to perform spectral detection, and measure the measurement contact at that time Perform temperature analysis on the temperature of the part, the measured spectral data and the mathematical model of spectral measurement established in step (2);
(5)在光谱检测的过程中继续采集测量部位的表层或表面.纹理图像, 通过 图像相减运算得到测量部位的运动状态, 如果在测量过程中发现测量部 位发生移动, 终止光谱检测, 装置重新进入测量条件复现流程。  (5) During the spectrum detection process, continue to collect the surface layer or surface of the measurement site. The texture image is obtained by the subtraction operation of the measurement site. If the movement of the measurement site is found during the measurement process, the spectrum detection is terminated and the device is restarted. Enter the measurement condition reproduction process.
本发明具有以下优点:  The invention has the following advantages:
(1)依据人体表层或表面的纹理特征实现测头在 xy平面位移和旋转角度 的重复定位;  (1) Repetitive positioning of the probe in the xy plane displacement and rotation angle according to the texture features of the human body surface or surface;
(2)人体表层或表面紋理特征可以是皮肤表面紋路如指纹和掌纹, 也可 以是皮肤表层附近血管、 骨骼及毛孔等体现出的纹理特征, 在一定条件 下是固定不变的;  (2) The surface texture or surface texture of the human body can be the texture of the skin surface, such as fingerprints and palm prints, or the texture features of blood vessels, bones and pores near the surface of the skin, which are fixed under certain conditions;
(3)依据光纤测头与测量部位的接触压力实现测头在垂直方向 z轴的重 复定位;  (3) According to the contact pressure between the optical fiber probe and the measurement site, the probe is repeatedly positioned in the vertical z-axis;
(4)采集测量部位的温度信息, 并将 引入到红外光谱分析过程中; (4) Collect the temperature information of the measurement site and introduce it into the infrared spectrum analysis process;
(5)在红外光谱的采集过程中监测人体测量部位的运动状态, 保证光谱 测量的正确性; (5) Monitor the movement status of the body measurement part during the acquisition of the infrared spectrum to ensure the accuracy of the spectrum measurement;
(6)利用小波变换减少测量部位紋理图像的数据量, 加速旋转角度校正 的进程;  (6) The wavelet transform is used to reduce the data amount of the texture image at the measurement site and accelerate the process of rotation angle correction;
(7)利用相位相关技术实现图像匹配, 减小图像失真和光照不均匀对图 像匹配的影响;  (7) Use phase correlation technology to achieve image matching, reduce the impact of image distortion and uneven illumination on image matching;
(8)利用曲面拟合技术使得图像匹配的精度达到亚像元级;  (8) Use surface fitting technology to make the image matching accuracy reach sub-pixel level
(9)采用环状边缘照明和正面直接摄像的方式简化纹理图像的采集; (9) Simplify the collection of texture images by means of ring-shaped edge lighting and direct front camera;
(10)利用伺服系统实现光纤测头的自动定位。 附图说明 (10) The servo system is used to realize the automatic positioning of the optical fiber probe. BRIEF DESCRIPTION OF THE DRAWINGS
图 1 : 测量条件复现系统框图;  Figure 1: Block diagram of measurement condition reproduction system;
图 2: 初始测量条件确定流程;  Figure 2: Initial measurement condition determination process;
图 3: 测量条件复现处理流程;  Figure 3: Measurement condition reproduction processing flow;
图 4: 基于掌纹识别的图像定位系统结构简图;  Figure 4: Structure of an image positioning system based on palmprint recognition;
图 5: 基于指纹识别的图像定位系统结构简图; 图 6a: 环状边缘照明装置结构主视图; Figure 5: A schematic diagram of the structure of an image positioning system based on fingerprint recognition; Figure 6a : a front view of the structure of a ring-shaped edge lighting device;
图 6b: 环状边缘照明装置结构侧视图;  Figure 6b: Side view of the structure of the ring-shaped edge lighting device;
图 7: 两幅图像位置偏移示意图;  Figure 7: Schematic diagram of position shift between two images;
图 8: 图像匹配处理流程;  Figure 8: Image matching processing flow;
图 9: 模板掌纹图像;  Figure 9: Template palm print image;
图 10: 当前采集掌纹图像;  Figure 10: The palmprint image is currently collected;
图 11 : 将当前采集图像旋转校正处理后与模板图像的相位相关结果; 图 12: 将当前采集图像经旋转与平移校正后与模板图像的叠合结果; 图 13a: 光纤测头与测量部位接触压力的检测示意图;  Figure 11: Phase correlation results of the currently acquired image after rotation correction processing with the template image; Figure 12: Superposition results of the current acquisition image after rotation and translation correction with the template image; Figure 13a: Fiber optic probe in contact with the measurement site Schematic diagram of pressure detection;
图 13b: 光纤测头与测量部位接触压力的检测俯视图;  Figure 13b: Top view of the contact pressure detection between the optical fiber probe and the measurement site;
图 14: 光纤测头与测量部位接触温度的检测示意图;  Figure 14: Schematic diagram of the contact temperature between the optical fiber probe and the measurement site;
图 15: 运动状态判别处理流程;  Figure 15: Motion state discrimination processing flow;
图 16: 采用测量条件复现系统前后光谱测量结果 CV值比较; 图 16-(1)测试对象 1 ;  Figure 16: Comparison of the CV values of the spectral measurement results before and after the system is reproduced using the measurement conditions; Figure 16- (1) Test object 1;
图 16-(2)测试对象 2;  Figure 16- (2) Test object 2;
图 16-(3)测试对象 3。 具体实施方式  Figure 16- (3) Test object 3. detailed description
下面结合附图和具体实施方式对本发明作进一步详细的论述。  The present invention is discussed in further detail below with reference to the drawings and specific embodiments.
1.测量条件复现系统的结构  1. Structure of measurement condition reproduction system
带有孔的载台的下方安装照明装置, 该装置设置有与载台相附的孔, 且轴线重合; 伺服定位装置由四部分组合而成: X轴伺服定位装置、 y轴 伺服定位装置、 z轴伺服定位装置、 z轴旋转伺服定位装置; z轴旋转伺 服定位装置安装在 z轴伺服定位装置上; 光纤测头安装在 Z轴旋转伺服 定位装置上, 并且中心轴与其旋转轴重合; 伺服定位装置中的位移机构 由步进电机驱动完成, 并由伺服控制器控制, 伺服控制器由计算机控制。 A lighting device is installed below the stage with a hole, the device is provided with a hole attached to the stage, and the axis coincides; the servo positioning device is composed of four parts: an X-axis servo positioning device, a y-axis servo positioning device, z-axis servo positioning device, z-axis rotation servo positioning device; z-axis rotation servo positioning device is installed on z-axis servo positioning device; optical fiber probe is installed on Z- axis rotation servo positioning device, and the center axis coincides with its rotation axis; servo The displacement mechanism in the positioning device is driven by a stepper motor and controlled by a servo controller, which is controlled by a computer.
在人体表层或表面的紋理特征中, 指纹和掌纹的纹理特征比较明显, 而且所处的部位便于进行无创的红外光谱检测。 依据手指和手掌上测量 部位的特点, 测量条件复现系统的结构也有所区别。  Among the texture features of the human body's surface layer or surface, the texture features of fingerprints and palm prints are relatively obvious, and the location is convenient for non-invasive infrared spectrum detection. According to the characteristics of the measurement parts on the fingers and palms, the structure of the measurement condition reproduction system is also different.
如图 4所示: 在基于掌纹识别的图像定位系统中, 光纤测头 8和 CCD 摄像机 1 同处被测部位 2如手掌内侧, 即图像匹配的信息为掌心区域的 紋理, 红外光谱的测量位置为手的掌心。 载台 3被 X轴伺服定位装置 4、 y轴伺服定位装置 5驱动, X轴伺服定位装置固定在底板上 7, y轴伺服 定位装置与 X轴正交固定在 X轴伺服定位装置上, 载台固定在 y轴伺服 定位装置上。 驱动光纤测头 8的伺服定位装置由两部分逐级组合而成: 即 z轴伺服定位装置 6、 z轴旋转伺服定位装置 9。 z轴伺服定位装置通 过垂直连接件固定在底板上 7, z轴旋转伺服定位装置安装在 z轴伺服定 位装置上。 光纤测头安装在 Z轴旋转伺服定位装置上, 并且中心轴与其 旋转轴重合。 伺服定位装置中的位移机构由步进电机驱动完成, 并由伺 服控制器控制, 伺服控制器由计算机控制。 . As shown in Figure 4: In an image positioning system based on palmprint recognition, the optical fiber probe 8 and the CCD The camera 1 is located at the same location 2 as the inside of the palm, that is, the image matching information is the texture of the palm area, and the measurement position of the infrared spectrum is the palm of the hand. The stage 3 is driven by the X-axis servo positioning device 4 and the y-axis servo positioning device 5. The X-axis servo positioning device is fixed on the base plate 7. The y-axis servo positioning device is orthogonally fixed to the X-axis servo positioning device on the X-axis servo positioning device. The stage is fixed on the y-axis servo positioning device. The servo positioning device for driving the optical fiber probe 8 is composed of two parts step by step: a z-axis servo positioning device 6 and a z-axis rotary servo positioning device 9. The z-axis servo positioning device is fixed on the base plate by a vertical connection. The z-axis servo positioning device is installed on the z-axis servo positioning device. The optical fiber probe is mounted on the Z- axis rotary servo positioning device, and the center axis coincides with its rotation axis. The displacement mechanism in the servo positioning device is driven by a stepping motor and controlled by a servo controller, which is controlled by a computer. .
如图 5所示: 在基于指紋识别的图像定位系统中,光纤测头 8和 CCD 摄像机 1分处被测部位 2如手指的两侧, 即图像匹配的信息为指紋, 红 外光谱的测量位置为手指的背面如指甲部位, 载台 3 固定不动。 驱动光 纤测头的伺服定位装置由四部分逐级组合而成: X轴伺服定位装置 4、 y 轴伺服定位装置 5、 Z轴伺服定位装置 6、 z轴旋转伺服定位装置 9。 X轴 伺服定位装置与底板 7相连, y轴伺服定位装置与 X轴正交固定在 X轴伺 服定位装置上, z轴伺服定位装置通过垂直连接件安装在 y轴伺服定位装 置上, Z轴旋转伺服定位装置安装在 z轴伺服定位装置上。 光纤测头 8安 装在 z轴旋转伺服定位装置 9上, 并且中心轴与其旋转轴重合。 伺服定 位装置中的位移机构由步进电机驱动完成, 并由伺服控制器控制, 伺服 控制器由计算机控制。 As shown in Figure 5: In an image positioning system based on fingerprint recognition, the optical fiber probe 8 and the CCD camera 1 are located at the measured part 2 such as both sides of a finger, that is, the information of image matching is the fingerprint, and the measurement position of the infrared spectrum is The back of the finger, such as the nail area, is fixed on the stage 3. The servo positioning device driving the optical fiber probe is composed of four parts step by step: X-axis servo positioning device 4, y-axis servo positioning device 5, Z- axis servo positioning device 6, z-axis rotary servo positioning device 9. The X-axis servo positioning device is connected to the bottom plate 7. The y-axis servo positioning device is orthogonally fixed to the X-axis servo positioning device on the X-axis servo positioning device. The z-axis servo positioning device is mounted on the y-axis servo positioning device through a vertical connection, and the Z- axis rotates. The servo positioning device is mounted on the z-axis servo positioning device. The optical fiber probe 8 is mounted on the z-axis rotation servo positioning device 9, and the center axis thereof coincides with its rotation axis. The displacement mechanism in the servo positioning device is driven by a stepping motor and controlled by a servo controller, which is controlled by a computer.
在以上两种应用的图像定位系统中, 手掌和手指放置在一个载台上, 正面朝下, 载台上有一通孔 10, 以便于掌纹或指紋图像的采集。 在载台 的下方安装了一个环状边缘照明装置 11, 该装置同样开有一通孔, 大小 与载台的孔相同, 轴线重合。  In the image positioning system of the above two applications, the palm and fingers are placed on a carrier with the front side facing down, and a through hole 10 is provided on the carrier to facilitate the collection of palm prints or fingerprint images. A ring-shaped edge lighting device 11 is installed under the carrier. The device also has a through hole, the same size as the hole of the carrier, and the axes coincide.
2. 照明的实现 '  2. Lighting realization ''
测量部位的照明采用环状边缘照明的方式实现。 为了便于掌紋和指紋 的图像釆集, 在手掌和手指的下方有一圆孔 10, 如图 4和图 5所示。 在 手掌和手指承载面的下方固定一个发光二极管 (LED)光源的固定结构,如 图 6a和图 6b所示。 该固定结构的材料为有机玻璃, 在圆孔 10的圆周上 均匀排布 24个凹槽 12, 每一凹槽可嵌入一个发光二极管 (LED), 凹槽 12 的底面为一斜面, 靠近中心处较高, 保证 LED向斜上方照射。 中心孔侧 面的有机玻璃打磨粗糙, 使得 LED漫反射照射测量部位。 Illumination of the measurement site is achieved by means of circular edge illumination. In order to facilitate the image collection of palm prints and fingerprints, there are circular holes 10 under the palms and fingers, as shown in FIGS. 4 and 5. A light emitting diode (LED) light source fixing structure is fixed below the palm and finger bearing surface, as shown in Figs. 6a and 6b. The material of the fixing structure is plexiglass, and the circumference of the circular hole 10 is 24 grooves 12 are evenly arranged, and each groove can be embedded with a light emitting diode (LED). The bottom surface of the groove 12 is an inclined surface, which is higher near the center to ensure that the LED is illuminated obliquely upward. The plexiglass on the side of the center hole is roughened, so that the diffuse reflection of the LED illuminates the measurement site.
3. 图像采集的实现  3. Implementation of image acquisition
为了避免因斜成像引起的图像畸变和随之产生的定位误差, CCD摄 像机 1安装在圆孔中心轴的正下方, 如图 4和图 5所示。  In order to avoid image distortion caused by oblique imaging and the subsequent positioning errors, the CCD camera 1 is installed directly below the center axis of the circular hole, as shown in Figs. 4 and 5.
在指紋定位系统中, 由于光纤测头在指甲的上方, 与 CCD摄像机的 位置不相冲突, 因此 CCD摄像机固定安装在指紋正下方的底板上。  In the fingerprint positioning system, because the optical fiber probe is above the nail and does not conflict with the position of the CCD camera, the CCD camera is fixedly installed on the bottom plate directly below the fingerprint.
在掌纹定位系统中, 要求光纤测头在测试状态下应位于手掌的正下 方, 这与掌纹的正面成像相冲突。 ·为了解决这一矛盾, 在将光纤测头安 装在 z轴旋转伺服定位装置的同时, 将 CCD摄像机安装在 Z轴伺服定位 装置上。 由于伺服定位装置是逐级组合而成, Z轴伺服定位装置可在 X轴 伺服定位装置和 y轴伺服定位装置的驱动下平移到不同的位置, 即在图 像采集阶段将 CCD摄像机移动到载台圆孔的正下方, 在光谱测量的阶段 将光纤测头移动到载台圆孔的正下方。 In the palmprint positioning system, it is required that the optical fiber probe should be directly under the palm in the test state, which conflicts with the front image of the palmprint. In order to solve this contradiction, the optical fiber probe is mounted on the z-axis rotary servo positioning device, and the CCD camera is mounted on the z- axis servo positioning device. As the servo positioning device is combined step by step, the Z axis servo positioning device can be translated to different positions under the drive of the X axis servo positioning device and the y axis servo positioning device, that is, the CCD camera is moved to the stage during the image acquisition stage. Just below the circular hole, move the fiber-optic probe directly under the circular hole of the stage in the spectrum measurement stage.
CCD摄像机输出标准视频信号, 输入到计算机系统中的图像采集卡。 The CCD camera outputs a standard video signal and inputs it to a frame grabber in a computer system.
4. 图像处理流程与算法 4. Image processing flow and algorithm
在进行光纤测头重复定位之前, 需要确定当前采集图像与模板图像的 偏移量。 图 7为两幅图像偏移量的示意图, 图中两个坐标轴的原点分别 为两幅图像的中心, 可看出两幅图像存在 xy平面的位移偏移量和旋转偏 移量。 图像采集卡采集掌纹和指纹图像后, 需要进入图像处理流程, 确 定两幅图像的位置偏移量, 然后以该位置偏移量驱动伺服定位系统, 使 得光纤测头到相应的位置。 图像处理流程分为三个阶段: (1)应用小波变 换及穷尽法进行旋转角度的确定; (2)利用相位相关进行像元精度的位移 变化检测; (3)利用二次曲面拟合进行亚像元精度的位移变化检测。  Before repeating the positioning of the optical fiber probe, the offset between the currently acquired image and the template image needs to be determined. Figure 7 is a schematic diagram of the offsets of the two images. The origins of the two coordinate axes in the figure are the centers of the two images. It can be seen that there is a displacement offset and a rotational offset of the xy plane in the two images. After the image capture card collects palm print and fingerprint images, it needs to enter the image processing process to determine the position offset of the two images, and then drive the servo positioning system with the position offset to make the optical fiber probe reach the corresponding position. The image processing flow is divided into three stages: (1) the determination of the rotation angle using wavelet transform and exhaustion method; (2) the use of phase correlation to detect the displacement change of pixel accuracy; (3) the use of quadratic surface fitting for sub- Pixel-precision displacement change detection.
图 8为图像处理流程图, 下面详细描述:  FIG. 8 is an image processing flowchart, which is described in detail below:
4.1应用小波变换及穷尽法进行旋转角度的确定  4.1 Determination of rotation angle using wavelet transform and exhaustion method
在多次测量条件下, 人体测量部位难以保证位置的完全一致, 其中包 括旋转角度和平面位移。 由于进行基于相位相关的图像匹配对旋转的角 度比较敏感, 因此在确定 xy平面位移之前需要确定当前采集图像与模板 图像的旋转角度。 Under multiple measurement conditions, it is difficult to ensure that the position of the body measurement site is completely consistent, including rotation angle and plane displacement. Because image matching based on phase correlation is more sensitive to the angle of rotation, it is necessary to determine the current acquired image and template before determining the xy plane displacement. The rotation angle of the image.
为了使得图像定位系统具有较高的精度, 釆集图像的分辨率为 In order to make the image positioning system have higher accuracy, the resolution of the set image is
768x582。 在手掌和手指的载台上有简单的限位装置, 手掌和手指的位移 偏移量在一个较小的范围内, 使得截取中央分辨率为 256x256的图像进 行图像匹配就可满足定位的要求。 然而, 基于相位相关的图像匹配需要 进行傅立叶变换, 运算量很大, 而且穷尽法需要进行多次相位相关运算。 图像分辨率即使为 256x256, 数据量仍很大, 直接对该分辨率的图像进 行穷尽法的相位相关需要很长的时间。 为了减少相位相关图像匹配运算 的数据量, 在进行穷尽法的旋转角度确定之前对当前采集图像和模板图 像进行小波变换, 将图像的分辨率降低到 32x32, 大大地降低了该图像 处理阶段的时间。 图像的小波多尺度表示一种金字塔式的结构, 这对于 进行由粗到细的快速图像匹配是非常有利的。 由小波变换的理论, 图像 经小波变换之后, 在尺度空间基本保持图像整体特征, 即低频部分, 且 经抽采后尺度空间部分相对原图数据量明显减少。 768x582. There are simple limit devices on the stage of the palms and fingers. The displacement of the palms and fingers is within a small range, which makes it possible to meet the positioning requirements by capturing images with a central resolution of 256x256 and performing image matching. However, image matching based on phase correlation requires a Fourier transform, which is computationally intensive, and the exhaustive method requires multiple phase correlation operations. Even if the image resolution is 256x256, the amount of data is still very large, and it takes a long time to directly perform the exhaustive phase correlation on the image with this resolution. In order to reduce the amount of data in the phase-dependent image matching operation, wavelet transform is performed on the currently acquired image and the template image before determining the rotation angle of the exhaustion method, reducing the image resolution to 32x32, which greatly reduces the time of this image processing stage. . The wavelet multi-scale of the image represents a pyramid-like structure, which is very beneficial for fast image matching from coarse to fine. According to the theory of wavelet transform, after the image is transformed by wavelet, the overall characteristics of the image are basically maintained in the scale space, that is, the low-frequency portion, and the amount of data in the scale space portion after extraction is significantly reduced compared to the original image.
采用穷尽法确定两幅图像的旋转角度偏移量的具体步骤如下: (1)假 设当前釆集图像相对于模板图像旋转变化的角度 A范围为 [Amin,AmaX], 根据旋转角度测量精度的要求, 确定角度穷尽的步长 ΔΑ; (2)将当前采集 图像旋转 Amin, 求其傅立叶变换, 并与模板图像进行相位相关, 求得对 应的相关峰信噪比,并保存到一个列表中; (3)将当前采集图像旋转 Amin+ΔΑ, 求其傅立叶变换, 并与模板图像进行相位相关, 求得对应的 相关峰信噪比, 保存在上面所述的列表中; (4)以步长为 ΔΑ逐渐增加当 前采集图像的旋转角度, 直至旋转角度达到 Amax, 在每一步求其傅立 叶变换, 并与模板图像进行相位相关, 求得对应的相关峰信噪比, 同样 保存在上面所述的列表中; (5)对列表进行遍历, 比较每一旋转角度对应 的相关峰信噪比, 相关峰信噪比最大值对应的旋转角度为两幅图像的角 度偏移量。 The specific steps of using the exhaustive method to determine the rotation angle offset of the two images are as follows: (1) Assume that the range A of the rotation change of the current set image relative to the template image is [Amin, Am aX ]. It is required to determine the step angle ΔΑ where the angle is exhausted; (2) Rotate the current acquired image by Amin, find its Fourier transform, and perform phase correlation with the template image, find the corresponding peak-to-noise ratio of the relevant peak, and save it to a list; (3) Rotate the current acquired image by Amin + ΔΑ, find its Fourier transform, and perform phase correlation with the template image, find the corresponding peak-to-noise ratio of the relevant peak, and save it in the list described above; (4) use the step size Gradually increase the rotation angle of the currently acquired image for ΔΑ until the rotation angle reaches Amax, find its Fourier transform at each step, and perform phase correlation with the template image to obtain the corresponding correlation peak signal-to-noise ratio, which is also stored in the above-mentioned In the list; (5) traverse the list, compare the correlation peak signal-to-noise ratio corresponding to each rotation angle, and the rotation angle corresponding to the maximum correlation peak signal-to-noise ratio Offset angle of the two images.
4.2利用相位相关进行像元精度的位移变化检测  4.2 Detection of pixel-accurate displacement changes using phase correlation
利用 4.1所述的小波变换和穷尽法确定当前釆集图像和模板图像的旋 转角度之后, 需要再次利用相位相关处理进行图像匹配, 确定两幅图像 的 xy平面的位移偏移量。 相位相关是图像一种受几何失真小的匹配算 法。 因为几何失真对图像的高频分量影响大, 而对低频分量影响小, 在 本装置中采用具有低通滤波器特性、 以傅立叶频谱为基础的相位相关算 法, 可以大大减小几何失真对匹配性能的影响。 After using the wavelet transform and exhaustion method described in 4.1 to determine the rotation angle of the current set image and template image, it is necessary to use phase correlation processing again for image matching to determine the displacement offset of the xy plane of the two images. Phase correlation is an image matching algorithm with little geometric distortion. Law. Because the geometric distortion has a large impact on the high-frequency components of the image and a small impact on the low-frequency components, a phase correlation algorithm based on the Fourier spectrum with a low-pass filter characteristic is used in this device, which can greatly reduce the geometric distortion on the matching performance Impact.
下式为相位相关运算的基本公式。 '
Figure imgf000011_0001
The following formula is the basic formula of the phase correlation operation. '
Figure imgf000011_0001
F ξ,η) Όξ,η)  F ξ, η) Όξ, η)
其中, ^和 分别为两幅图像 (当前采集图像和模板图像)的傅立叶 变换的结果。 由 (1)式及傅立叶变换的理论可知, 该相位谱包含了两幅图 像的位置平移信息, 而且它是一个频谱幅度在全频域内为 1 的功率谱。 对 (1)式进行逆傅立叶变换可知, 相位相关函数是一个位于两图位置偏移 (x0, yO)处的 δ脉冲函数, 也称之为相关峰。 当两幅图像完全相似时, 其 值为 1, 反之为 0。 因此, 在本装置中利用两幅图像的相位相关运算结果 确定图像的位移偏移量。 Among them, ^ and are the results of Fourier transform of two images (currently acquired image and template image), respectively. According to the formula (1) and the theory of Fourier transform, the phase spectrum contains the position shift information of the two images, and it is a power spectrum whose spectral amplitude is 1 in the full frequency domain. According to the inverse Fourier transform of (1), it can be seen that the phase correlation function is a delta pulse function located at the position offset (x0, yO) of the two figures, which is also called the correlation peak. When the two images are completely similar, the value is 1, otherwise it is 0. Therefore, the phase correlation calculation results of the two images are used in this device to determine the displacement amount of the images.
在本装置中, 利用相位相关进行图像匹配, 从而确定位移偏移量时 具有以下特点: (1)较大的位移检测范围。当两幅图像仅存在位移变化时, 能够检出 xy平面最大的位移偏移量达到图像宽度的一半。 (2)尖锐的相 关峰。 当两幅图像完全相关时, 计算得到的 δ脉冲函数有非常尖锐的相关 峰, 可实现两幅图像中心点位置的精确匹配。 (3)对图像灰度依赖小。 当 两幅图像之间存在灰度差时, 计算得到的 δ脉冲函数的位置不发生变化, 只是幅度上有所区别, 因此具有很强的抗图像遮挡的能力。 (4)两幅图像 之间的旋转偏移量对图像匹配的结果影响较大。 当两幅图像存在一个旋 转角度偏移量时, 其傅立叶变换的频谱也随之旋转。 当旋转角度偏移量 大于 5°时, 即使两幅图像完全相同, 其相关峰值衰减为零, 即不能依据 相关峰来确定两幅图像在 xy平面的位移偏移量。 因此在上一步骤通过判 断相关峰的大小确定旋转角度, 并对图像旋转该角度后进行相位相关, 最终得到两幅图像的位置偏移量。  In this device, phase matching is used for image matching to determine the displacement offset with the following characteristics: (1) Large displacement detection range. When there is only a displacement change between the two images, the maximum displacement offset of the xy plane can be detected to reach half the image width. (2) Sharp related peaks. When the two images are completely correlated, the calculated δ pulse function has a very sharp correlation peak, which can achieve an accurate matching of the positions of the center points of the two images. (3) Little dependence on image gray. When there is a gray level difference between the two images, the position of the calculated δ pulse function does not change, but there is a difference in amplitude, so it has a strong ability to resist image occlusion. (4) The amount of rotation offset between the two images has a greater impact on the results of image matching. When there is a rotation angle offset between the two images, the spectrum of the Fourier transform is also rotated. When the rotation angle offset is greater than 5 °, even if the two images are identical, the correlation peak attenuation is zero, that is, the displacement offset of the two images in the xy plane cannot be determined based on the correlation peak. Therefore, in the previous step, the rotation angle is determined by judging the size of the correlation peak, and the image is rotated by this angle to perform phase correlation, and finally the position offset of the two images is obtained.
通过上述图像处理过程, 可以确定当前采集图像和模板图像的位置偏 移量, 该偏移量的最小单位为 1个像素。 为了进一步精确确定位置偏移 量, 提高系统的定位精度, 采用下述所述的基于二次曲面拟合的亚像元 匹配算法。 Through the above image processing process, the position offset of the currently acquired image and the template image can be determined, and the minimum unit of the offset is 1 pixel. To further accurately determine the position offset To improve the positioning accuracy of the system, the following sub-pixel matching algorithm based on quadric surface fitting is used.
4.3利用二次曲面拟合进行亚像元精度的位移变化检测  4.3 Sub-pixel precision displacement change detection using quadric surface fitting
在本系统中, 采用曲面拟合方法, 该方法具有快速、容易计算的优点, 精度能达到亚像元。 曲面拟合法的思想是: 以像元级上的最佳匹配点为 中心, 按相似性度量进行曲面拟合, 然后通过相应的数学方法计算得到 极值点的精确位置。 本装置采用相位相关的相关系数作为相似性度量特 征, 选择二次曲面作为拟合函数, 在计算中采用多变量最小二乘回归法 确定极值点的精确位置。  In this system, a surface fitting method is adopted. This method has the advantages of fast and easy calculation, and the accuracy can reach sub-pixels. The idea of the surface fitting method is: take the best matching point at the pixel level as the center, perform surface fitting according to the similarity measure, and then calculate the exact position of the extreme point by corresponding mathematical methods. This device uses the correlation coefficient of phase correlation as the similarity measurement feature, selects the quadric surface as the fitting function, and uses the multivariate least squares regression method in the calculation to determine the exact location of the extreme point.
二次曲面拟合函数采用公式为:  The quadric fitting function uses the formula:
PC(x,y) = ax + by + cxy + dx + ey + f PC (x, y) = ax + by + cxy + dx + ey + f
其中, PC(x,y)为对应于位置 (x,y)的相位相关值。 上述函数可以写成 如下形式: Among them, PC (x, y) is a phase correlation value corresponding to the position (x, y). The above function can be written as:
AX = B 式中,
Figure imgf000012_0001
AX = B where,
Figure imgf000012_0001
本装置在拟合计算中采用多变量最小二乘回归法, 使得计算简单、 准 确。 在计算过程中, 将向量 X作为回归系数, 并假设随机变量 Β的取值 依赖于矩阵 Α中的自变量, 回归系数的求取即为拟合函数的系数。 在求 得拟合函数的系数之后, 可以利用下式求得亚像元精度的图像偏移的精 确位置。  The device adopts the multivariate least squares regression method in the fitting calculation, which makes the calculation simple and accurate. In the calculation process, the vector X is used as the regression coefficient, and it is assumed that the value of the random variable B depends on the independent variable in the matrix A. The regression coefficient is obtained as the coefficient of the fitting function. After the coefficients of the fitting function are obtained, the precise position of the image offset with sub-pixel accuracy can be obtained using the following formula.
2db - ce 2ae― dc  2db-ce 2ae― dc
X = - y  X =-y
c2 - 4ab c2 - ab 图 9为模板掌紋图像, 图 10为当前采集掌紋图像, 两幅图像的分辨 率为 256 X 256。 可看出两幅图像存在明显的角度和位移偏移。 将角度增 量设置为 0.1度, 应用小波变换及穷尽法进行旋转角度检测, 得到当前 釆集掌纹图像相对于模板掌紋图像的旋转角度为逆时针 13.9度。 将当前 采集掌紋图像逆时针旋转 13.9度后, 与模板掌紋图像进行相位相关匹配 处理, 图 11为处理得到相关函数的平面分布图, 依据最大相关峰的位置 得到两幅图像在水平方向的偏移量为 -20像素, 垂直方向的偏移量为 64 像素。 将当前采集掌纹图像图逆时针旋转 13.9度, 并在水平和垂直方向 分别平移 -20像素和 64像素, 然后与模板掌纹图像进行叠加。 如图 12所 示, 从叠加图像可看出上述利用相位相关的旋转角度和位移偏移量的检 测是正确、 有效的。 在利用相位相关进行像元精度的位移变化检测的基 础上, 再利用二次曲面拟合进行亚像元精度的位移变化检测, 得到两幅 图像在水平方向的偏移量为 -20.08像素, 垂直方向的偏移量为 64.20像 素。 通过进一步评价, 测量条件复现装置的旋转检测精度可达到 0.1度, 平移检测精度可达到 1/50像元。 c 2 - 4ab c 2 - ab FIG. 9 as a template palmprint image 10 is currently acquired palmprint, the resolution of the image is two 256 X 256. It can be seen that there is a significant angular and displacement shift in the two images. Increase angle The amount is set to 0.1 degree, and the rotation angle is detected by applying wavelet transform and exhaustion method, and the rotation angle of the palm image of the current set relative to the template palm image is 13.9 degrees counterclockwise. After the currently collected palm print image is rotated counterclockwise by 13.9 degrees, phase correlation matching processing is performed with the template palm print image. Figure 11 is a flat distribution map of the correlation function obtained by processing, and the horizontal offset of the two images is obtained according to the position of the maximum correlation peak. The amount is -20 pixels, and the vertical offset is 64 pixels. Rotate the currently collected palm print image map 13.9 degrees counterclockwise, and translate -20 pixels and 64 pixels in the horizontal and vertical directions, respectively, and then superimpose the template palm print image. As shown in FIG. 12, it can be seen from the superimposed image that the above-mentioned detection of the rotation angle and displacement offset using phase correlation is correct and effective. On the basis of using phase correlation to detect displacement changes in pixel accuracy, and then using quadric surface fitting to detect displacement changes in sub-pixel accuracy, the horizontal offset of the two images is -20.08 pixels, vertical The offset in the direction is 64.20 pixels. Through further evaluation, the rotation detection accuracy of the measurement condition reproduction device can reach 0.1 degrees, and the translation detection accuracy can reach 1/50 of a pixel.
5.接触参数的测量与光纤测头 z向位置的调整  5. Measurement of contact parameters and adjustment of the z-direction position of the optical fiber probe
图像处理流程完成之后将得到测量部位在 xy平面的相对位置即 xy轴 的偏移量和旋转角度, 然后依靠 xy平面伺服定位装置驱动光纤测头, 使 测头轴线与测量部位的中心点重合, 并通过 z轴旋转定位装置将光纤测 头旋转一定角度, 该角度为图像处理流程得到的旋转角度, 最后通过 z 向逼近伺服装置驱动光纤测头, 使光纤测头与测量部位相接触。  After the image processing process is completed, the relative position of the measurement site in the xy plane, that is, the offset and rotation angle of the xy axis, and then rely on the xy plane servo positioning device to drive the optical fiber probe so that the probe axis coincides with the center point of the measurement site. The optical fiber probe is rotated by a certain angle through the z-axis rotation positioning device, which is the rotation angle obtained by the image processing process. Finally, the optical fiber probe is driven by the z-direction approaching servo device to bring the optical fiber probe into contact with the measurement site.
光纤测头与测量部位接触后, 压力传感器将输出接触压力信号, 并输 入到信号采集卡中, 计算机将采集到的接触压力信号与预设值相比较, 然后调整光纤测头在 z向的相对位置, 使得光纤测头与测量部位的接触 压力基本保持在预设压力值附近。 接触压力的测量如图 13所示。 在圆柱 形光纤测头的顶端固定一圆环状的气路托架 13, 在托架与光纤测头 8之 间的环形区域内安装一根硅胶管 14, 硅胶管内部充满空气, 并与一个硅 压力传感器相连, 构成一个封闭的气路。 依据测头与测量部位的接触压 力实现测头在垂直方向 z轴的重复定位。 当人体的测量部位与光纤测头 相接触时, 由于测量部位压挤充气的硅胶管, 引起封闭气路中的压力发 生变化, 即可得到人体测量部位与光纤测头的接触压力。  After the optical fiber probe contacts the measurement part, the pressure sensor will output the contact pressure signal and input it into the signal acquisition card. The computer compares the collected contact pressure signal with a preset value, and then adjusts the relative position of the optical fiber probe in the z direction. Position, so that the contact pressure between the optical fiber probe and the measurement site is basically maintained near a preset pressure value. The measurement of contact pressure is shown in Figure 13. A circular-shaped gas path bracket 13 is fixed on the top of the cylindrical optical fiber probe, and a silicone tube 14 is installed in the annular area between the bracket and the optical fiber probe 8. The inside of the silicone tube is filled with air and connected with a The silicon pressure sensor is connected to form a closed gas path. Repeated positioning of the probe in the vertical z-axis is achieved based on the contact pressure between the probe and the measurement location. When the measurement part of the human body is in contact with the optical fiber probe, the pressure in the closed gas path changes due to the pressure of the inflatable silicone tube at the measurement part, and the contact pressure of the human measurement part and the optical fiber probe can be obtained.
光纤测头与测量部位相接触后, 计算机能够通过温度传感器和信号采 集卡将两者的接触温度采集, 并将其数据作为数学模型建立和红外光谱 分析的参数。 温度传感器 15为薄膜式热敏电阻, 粘贴在光纤测头的测量 面, 如图 14所示。 光纤测头内安装有两束光纤, 一束分布在测头的中心 区域, 另一束分布在测头的一个圆环区域。 在两束光纤之间的中间区域 均匀贴有 4个温度传感器, 在光纤束 16的外围区域沿圆周均匀贴有 10 个温度传感器。 信号采集卡将 14个检测点的温度一起采集到计算机系统 中, 并进行平均处理, 保证温度测量的准确性。 After the optical fiber probe is in contact with the measurement site, the computer can The card collects the contact temperature of the two, and uses the data as parameters for mathematical model establishment and infrared spectrum analysis. The temperature sensor 15 is a thin-film thermistor, and is attached to the measurement surface of the optical fiber probe, as shown in FIG. 14. There are two bundles of optical fibers installed in the optical fiber probe, one bundle is distributed in the center area of the probe, and the other bundle is distributed in a ring area of the probe. Four temperature sensors are evenly pasted in the middle area between the two optical fibers, and ten temperature sensors are evenly pasted along the circumference in the peripheral area of the optical fiber bundle 16. The signal acquisition card collects the temperature of the 14 detection points into the computer system and performs average processing to ensure the accuracy of the temperature measurement.
6.运动状态判别  6. Exercise state discrimination
当光纤测头运动到相应位置., 并与测量部位相接触后, 红外光路检测 系统幵始工作。 在红外光谱检测的过程中, 要求测量部位的相对位置基 本保持固定。 为了达到这一目的, 在驱动光纤测头到计算后的位置并启 动光谱测量的同时, CCD摄像机采集光纤测头与测量部位接触区域 -载台 圆孔位置的图像, 并将其作为运动状态判别的模板图像; 在光谱测量的 过程中, 计算机每进行 5个波长的光谱测量后, 采集同一区域的图像, 并与运动状态判别的模板图像进行相减运算, 依据运算的结果判断在光 谱测量的过程中测量部位是否发生了移动。 如果在光谱测量过程中两幅 图像相减的结果基本为零, 表明测量部位的相对位置未发生变化, 继续 光谱的测量循环, 直至光谱测量结束并显示测量结束提示信息, 否则表 明检测到测量部位的相对位置发生了变化, 中止光谱测量, 提示测量数 据无效, 并要求用户重新进入测量条件复现流程, 保证了光谱测量的正 确性。 图 15为运动状态判别处理流程。  When the optical fiber probe moves to the corresponding position and comes into contact with the measurement site, the infrared optical path detection system starts to work. In the process of infrared spectrum detection, the relative position of the measurement site is required to be kept basically fixed. In order to achieve this, while driving the optical fiber probe to the calculated position and starting the spectral measurement, the CCD camera collects an image of the contact area of the optical fiber probe and the measurement site-the position of the round hole on the stage, and uses it as a motion state to determine In the process of spectral measurement, the computer collects an image of the same area after every five wavelengths of spectral measurement, and performs a subtraction operation with the template image discriminated by the motion state. Whether the measurement site has moved during the process. If the result of the subtraction of the two images during the spectrum measurement process is basically zero, it indicates that the relative position of the measurement site has not changed. Continue the spectrum measurement cycle until the end of the spectrum measurement and display the end of measurement prompt message, otherwise it indicates that the measurement site is detected The relative position has changed, and the spectrum measurement was suspended, indicating that the measurement data is invalid, and requiring the user to re-enter the measurement condition reproduction process to ensure the accuracy of the spectrum measurement. Fig. 15 is a processing flow for judging a motion state.
7.红外光谱测量分析结果  7. Infrared spectrum measurement and analysis results
利用上述测量条件复现方法确定光纤测头的相对位置后, 可进入光谱 测量。 为了比较采用测量条件复现系统前后光谱测量的结果, 针对掌心 部位进行了一系列测试。 测试方案有两种: (1)测量条件复现装置不工作 时光谱数据的测量: 首先将光纤测头的中心轴与载台圆孔的中心轴调整 为重合, 然后手动调整光纤测头在垂直方向的位置, 使得光纤测头与掌 心有一定的接触压力, 同时 z轴旋转伺服定位机构复位到零位, 被测对 象的手掌在每次光谱测量完成后抬起, 重新放置到载台时伺服定位装置 不工作, 即光纤测头的位置固定不变; 手掌重复放置到载台 10次, 每次 获得 10组光谱数据。 (2)测量条件复现装置工作时光谱数据的测量: 在第 1组光谱数据测量之前确定并保存测量条件如: 模板掌紋图像、 接触压 力, 被测对象的手掌在每次光谱测量完成后仍然抬起, 但是在每次重新 放置后, 在进入光谱数据测量之前, 利用测量条件复现装置复现第一组 光谱数据测量前保存的测量条件, 最后进行光谱数据的测量; 同样将手 掌重复放置到载台 10次, 共获得 10组光谱数据。 After the relative position of the optical fiber probe is determined by using the above-mentioned measurement condition reproduction method, the spectrum measurement can be entered. In order to compare the results of the spectral measurements before and after using the measurement condition reproduction system, a series of tests were performed on the palm area. There are two kinds of test schemes: (1) Measurement conditions: The measurement of the spectral data when the reproduction device is not working: First, adjust the center axis of the optical fiber probe and the center axis of the round hole of the stage to coincide, and then manually adjust the optical fiber probe to be vertical The orientation position makes the optical fiber probe have a certain contact pressure with the palm. At the same time, the z-axis rotary servo positioning mechanism is reset to the zero position. The palm of the measured object is lifted after each spectrum measurement is completed, and the servo is reset when it is placed on the stage. The positioning device does not work, that is, the position of the optical fiber probe is fixed; the palm is repeatedly placed on the carrier 10 times, each time 10 sets of spectral data were obtained. (2) Measurement conditions Measurement of the spectral data during the operation of the reproduction device: Determine and save the measurement conditions before the measurement of the first group of spectral data, such as: template palm print image, contact pressure, and the palm of the measured object remains after each spectral measurement is completed. Lift up, but after each relocation, before entering the spectral data measurement, use the measurement condition reproduction device to reproduce the measurement conditions saved before the first set of spectral data measurement, and finally perform the measurement of the spectral data; also place the palm repeatedly Go to the carrier 10 times and obtain 10 sets of spectral data.
利用以上测试方案, 对 3名不同的人体进行了光谱数据的测量, 统计 结果如图 16-(1)、 图 16-(2)、 图 16- (3)所示。 统计图的评价指标为 10组 光谱数据中每个波长点对应光强的最大波动与光强的比值, 縮写为 CV。 图中虚线表示的是测量条件复现装置不工作时光谱数据的统计结果, 图 中实线表示的是测量条件复现装置工作时光谱数据的统计结果。  Using the above test scheme, the spectral data of three different humans were measured. The statistical results are shown in Figure 16- (1), Figure 16- (2), and Figure 16- (3). The evaluation index of the statistical graph is the ratio of the maximum fluctuation of the light intensity corresponding to each wavelength point in the 10 sets of spectral data to the light intensity, which is abbreviated as CV. The dashed line in the figure indicates the statistical result of the spectral data when the measurement condition reproduction device is not working, and the solid line in the figure indicates the statistical result of the spectral data when the measurement condition reproduction device is in operation.
通过比较不同人体的 CV值, 可得出以下结论: (1)测量部位的位置变 化对光谱测量有较大的影响,必须采用相应的措施来限定位置的变化; (2) 本发明提出的基于人体表层或表面纹理特征和接触压力的测量条件复现 装置能大大减小测量部位位置改变对光谱测量的影响, 提高光谱的测量 稳定性和重复性。  By comparing the CV values of different human bodies, the following conclusions can be drawn: (1) The change in the position of the measurement site has a greater impact on the spectral measurement, and corresponding measures must be adopted to limit the change in position; (2) The present invention is based on The human body surface layer or surface texture feature and contact pressure measurement condition reproduction device can greatly reduce the influence of the change of the position of the measurement location on the spectrum measurement, and improve the measurement stability and repeatability of the spectrum.

Claims

权 利 要 求 Rights request
1.一种基于人体表层或表面纹理特征和接触压力的测量条件复现方 法, 其特征是主要包括以下步骤: . 1. A method for reproducing measurement conditions based on human body surface or surface texture features and contact pressure, which is characterized in that it mainly includes the following steps:
(1)人体表层或表面纹理特征的图像拾取、 匹配和定位;  (1) Image pickup, matching and positioning of human body surface or surface texture features;
(2)根据表层或表面纹理图像匹配的结果进行 xy平面位移和 z轴旋转角 度等测量位置的定位;  (2) positioning of measurement positions such as xy plane displacement and z-axis rotation angle according to the results of surface layer or surface texture image matching;
(3)依据测量元件与人体表层或表面纹理所在部位的接触压力进行 z方 向的定位;  (3) Positioning in the z direction according to the contact pressure between the measuring element and the surface layer or surface of the human body;
(4)在测量元件工作的过程中对人体表层或表面紋理特征采集图像, 监 测测量部位的运动状态。  (4) During the operation of the measuring element, collect images of human body surface or surface texture features, and monitor the movement state of the measurement part.
2.如权利要求 1 所述的一种基于人体表层或表面紋理特征和接触压力 的测量条件复现的方法, 其特征是所述的步骤为: '  2. A method of reproducing measurement conditions based on human body surface or surface texture features and contact pressure according to claim 1, wherein the steps are:
(1)对人体测量部位的表层或表面纹理进行图像采集, 表面纹理特征是 皮肤表面紋路如: 指纹和掌紋等, 表层是指皮肤表层附近血管、 骨骼及 毛孔等体现出的纹理特征, 依据人为预设的测量位置将保存图像作为模 板图像, 并附加接触压力信息;  (1) Image collection of the surface layer or surface texture of the body measurement site. The surface texture is characterized by the surface texture of the skin such as fingerprints and palm prints. The surface layer refers to the texture features reflected by blood vessels, bones and pores near the surface of the skin. The preset measurement position will save the image as a template image and attach the contact pressure information;
(2)人体测量部位重新放置后, 首先采集测量部位的表层或表面纹理图 像, 然后通过图像处理的方法对当前图像与模板图像进行匹配识别, 得 到当前测量部位与模板图像建立时测量部位在 xy平面位移和 z轴旋转角 度的偏移量;  (2) After the human body measurement part is repositioned, first collect the surface layer or surface texture image of the measurement part, and then match and identify the current image and the template image through image processing to obtain the measurement part at xy when the current measurement part and the template image are established. Offset of plane displacement and z-axis rotation angle;
(3)依据步骤 (2)得到的 xy平面位移和 Z轴旋转角度的偏移量驱动伺服 装置使测量元件到达相应位置, 然后依据测量元件与测量部位的接触压 力调整测量元件在 z轴方向的位置; (3) Drive the servo device according to the displacement of the xy plane displacement and the Z axis rotation angle obtained in step (2) to make the measuring element reach the corresponding position, and then adjust the measuring element in the z axis direction according to the contact pressure between the measuring element and the measurement part position;
(4)启动测量元件进入工作状态, 继续采集测量部位的表层或表面纹理 图像, 通过图像相减运算得到测量部位的运动状态。 - (4) Start the measuring element into the working state, continue to collect the surface layer or surface texture image of the measurement part, and obtain the movement state of the measurement part through the image subtraction operation. -
3.如权利要求 1、 2所述的人体表层或表面紋理特征和接触压力的测量 条件复现的方法, 其特征是所述的当前图像与模板图像进行匹配的方法 是: (1)应用小波变换及穷尽法进行旋转角度的确定, 小波变换减少测量部 位纹理图像的数据量, 加速旋转角度校正的进程; 3. The method for reproducing the measurement conditions of human body surface layer or surface texture features and contact pressure according to claim 1, 2, wherein the method of matching the current image with the template image is: (1) The wavelet transform and exhaustion method are used to determine the rotation angle. The wavelet transform reduces the data amount of the texture image at the measurement site and accelerates the process of rotation angle correction;
(2)利用相位相关进行像元精度的位移变化检测, 减小图像失真和光照 不均匀对图像匹配的影响;  (2) Use phase correlation to detect pixel-accurate displacement changes to reduce the effects of image distortion and uneven illumination on image matching;
(3)利用曲面拟合技术进行亚像元精度的位移变化检测, 使得图像匹配 的精度达到亚像元级。  (3) Use surface fitting technology to detect displacement changes of sub-pixel accuracy, so that the accuracy of image matching reaches sub-pixel level.
4.一种基于人体表层或表面紋理特征和接触压力的测量条件复现方法 的光谱测量装置, 它是由光纤测头, x、 y、 z方向的直线位移和 z轴角 度伺服定位装置、 载台、 照明装置、 摄像机、 图像采集卡、 光谱检测系 统和计算机处理系统组成; 其特征是测量部位与测头设置有接触压力传 感器、 温度传感器和对压力、 温度信号的处理及采集装置。  4. A spectroscopic measuring device based on a method for reproducing the conditions of a human body's surface layer or surface texture and a measurement condition of a contact pressure, which is a fiber-optic probe, a linear displacement in the x, y, and z directions and a z-axis angle servo positioning device. It consists of a table, a lighting device, a camera, an image acquisition card, a spectrum detection system, and a computer processing system. It is characterized in that a contact pressure sensor, a temperature sensor, and a processing and acquisition device for pressure and temperature signals are provided at the measurement site and the probe.
5.如权利要求 4 所述的一种基于人体表层或表面紋理特征和接触压力 的测量条件复现方法的光谱测量装置, 其特征是所述的装置中, 带有孔 的载台的下方安装照明装置, 该装置设置有与载台相附的孔, 且轴线重 合; 伺服定位装置由四部分组合而成: X轴伺服定位装置、 y轴伺服定位 装置、 Z轴伺服定位装置、 z轴旋转伺服定位装置; z轴旋转伺服定位装 置安装在 z轴伺服定位装萱上; 光纤测头安装在 z轴旋转伺服定位装置 上, 并且中心轴与其旋转轴重合; 伺服定位装置中的位移机构由步进电 机驱动完成, 并由伺服控制器控制, 伺服控制器由计算机控制。 The spectrum measurement device based on the method for reproducing the measurement conditions of the human body surface layer or surface texture and the contact pressure according to claim 4, wherein the device is installed below a carrier with holes. Lighting device, the device is provided with a hole attached to the stage, and the axis coincides; the servo positioning device is composed of four parts: X-axis servo positioning device, y-axis servo positioning device, Z- axis servo positioning device, z-axis rotation Servo positioning device; The z-axis rotary servo positioning device is installed on the z-axis servo positioning device; the optical fiber probe is installed on the z-axis rotary servo positioning device, and the central axis coincides with its rotation axis; The displacement mechanism in the servo positioning device consists of steps The input motor is driven and controlled by a servo controller, which is controlled by a computer.
.  .
6.如权利要求 4 所述的一种基于人体表层或表面紋理特征和接触压力 的测量条件复现方法的光谱测量装置, 其特征是所述的照明装置是照明 采用环状边缘照明的方式实现: 在载台的下方设置有发光二极管光源的 固定结构, 该固定结构的材料为有机玻璃, 在圆孔的圆周上均匀排布 24 个凹槽, 每一凹槽可嵌入一个发光二极管 (LED), 凹槽的底面为一斜面, 靠近中心处较高, 保证 LED向斜上方照射; 中心孔侧面的有机玻璃打磨 粗糙, 使得 LED漫反射照射测量部位。 6. A spectrum measurement device based on a method for reproducing the measurement conditions of human body surface or surface texture features and contact pressure according to claim 4, characterized in that said illumination device is realized by means of ring-shaped edge illumination. : A fixed structure of light emitting diode light source is arranged below the stage. The material of the fixed structure is plexiglass. 24 grooves are evenly arranged on the circumference of the circular hole, and each groove can be embedded with a light emitting diode (LED). The bottom surface of the groove is an inclined surface, which is higher near the center to ensure that the LED is illuminated obliquely upwards; the plexiglass on the side of the center hole is roughened so that the LED diffuse reflection irradiates the measurement site.
7.如权利要求 4 所述的一种基于人体表层或表面纹理特征和接触压力 的测量条件复现方法的光谱测量装置, 其特征是所述的装置中, CCD摄 像机安装在载台的下方, 采用正面直接摄像的方式; 当测头与摄像机在 同一侧时, 将光纤测头安装在 Z轴旋转伺服定位装置的同时, 将 CCD摄 像机安装在 z轴伺服定位装置上。 7. A spectrum measurement device based on a method for reproducing the measurement conditions of human body surface or surface texture features and contact pressure according to claim 4, wherein in the device, a CCD camera is installed below the stage, Adopt direct front camera; when the probe is on the same side as the camera, the optical fiber probe is installed in the Z- axis rotation servo positioning device, and the CCD is taken. The camera is mounted on the z-axis servo positioning device.
8.如权利要求 4 所述的一种基于人体表层或表面纹理特征和接触压力 的测量条件复现方法的光谱测量装置, 其特征是所述的接触压力传感器 是在光纤测头的顶端固定一环状的气路托架, 在托架与光纤测头之间的 环形区域内安装具有弹性的硅胶管, 管内部充满空气, 并与压力传感器 相连, 构成一个封闭的气路, 依据测头与测量部位的接触压力实现测头 在垂直方向 z轴的重复定位。  The spectrum measurement device based on the method for reproducing the measurement conditions of the human body surface layer or surface texture and the contact pressure according to claim 4, wherein the contact pressure sensor is fixed at the top of the optical fiber probe. An annular gas path bracket is installed with a flexible silicone tube in the annular area between the bracket and the optical fiber probe. The tube is filled with air and connected to a pressure sensor to form a closed gas path. The contact pressure of the measurement site enables repeated positioning of the probe in the vertical z-axis.
9.如权利要求 4 所述的一种基于人体表层或表面纹理特征和接触压力 的测量条件复现方法的光谱测量装置, 其特征是所述的温度传感器为薄 膜式热敏电阻, 粘贴在光纤测头的测量面。  The spectrum measurement device based on the method for reproducing the measurement conditions of human body surface layer or surface texture characteristics and contact pressure according to claim 4, wherein the temperature sensor is a thin film thermistor and is attached to an optical fiber Measuring surface of the probe.
10.如权利要求 4所述的一种基于人体表层或表面紋理特征和接触压力 的测量条件复现方法的光谱测量装置, 其特征在于它包括下列步骤:  10. A spectrum measurement device based on a method for reproducing a measurement condition of a human body surface layer or surface texture feature and contact pressure according to claim 4, characterized in that it comprises the following steps:
(1)建立模板模型, 其步骤为: 驱动摄像机到初始位置, 采集人体表层 或表面的紋理信息, 将图像保存到模板图像库, 驱动光纤测头到图像中 心, 驱动光纤测头与测量部位接触, 并保持设定的接触压力, 测量接触 部位的温度;  (1) Create a template model, the steps are as follows: drive the camera to the initial position, collect the texture information of the human body surface or surface, save the image to the template image library, drive the optical fiber probe to the image center, and drive the optical fiber probe to contact the measurement site , And maintain the set contact pressure, and measure the temperature of the contact area;
(2)在步骤 (1)的条件下进行光谱测量, 并依据光谱数据创建该测量部 位特定的光谱测量数学模型;  (2) Perform spectral measurement under the conditions of step (1), and create a specific mathematical model of spectral measurement for the measurement location based on the spectral data;
(3)人体测量部位重新放置后, 首先采集人体测量部位的表层或表面纹 理图像, 然后通过图像处理的方法对当前图像与模板图像进行匹配识别, 得到当前测量部位与模板图像建立时测量部位在 xy平面位移和 z轴旋转 角度的偏移量, 然后驱动伺服装置使光纤测头到达相应位置, 然后依据 光纤测头与测量部位的接触压力调整测量元件在 z轴方向的位置;  (3) After the human body measurement part is repositioned, first collect the surface layer or surface texture image of the human body measurement part, and then match and identify the current image and the template image through the image processing method, and obtain that the measurement part is located when the current measurement part and the template image are established. xy plane displacement and z axis rotation angle offset, then drive the servo device to make the optical fiber probe reach the corresponding position, and then adjust the position of the measuring element in the z axis direction according to the contact pressure between the optical fiber probe and the measurement part;
(4)利用光纤测头和光谱检测系统进行光谱检测, 测量出当时测量接触 部位的温度, 与测量得到的光谱数据代入到步骤 (2)建立的光谱测量数学 模型进行数据分析;  (4) Use optical fiber probes and spectral detection systems to perform spectral detection, measure the temperature of the measurement contact site at that time, and substitute the measured spectral data into step (2) The mathematical model of spectral measurement established for data analysis;
(5)在光谱检测的过程中继续采集测量部位的表层或表面纹理图像, 通 过图像相减运算得到测量部位的运动状态, 如果在测量过程中发现测量 部位发生移动, 终止光谱检测, 装置重新进入测量条件复现流程。  (5) During the spectrum detection process, continue to collect the surface or surface texture image of the measurement part, and obtain the movement state of the measurement part through the image subtraction operation. If the measurement part is found to move during the measurement process, the spectral detection is terminated, and the device re-enters Measurement condition reproduction process.
PCT/CN2003/000813 2002-11-22 2003-09-24 Apparatus and method for the reproduction of measuring conditions based on the body surface or surface veins and the contacting pressure WO2004056262A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
AU2003272847A AU2003272847A1 (en) 2002-11-22 2003-09-24 Apparatus and method for the reproduction of measuring conditions based on the body surface or surface veins and the contacting pressure

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN02148848.7 2002-11-22
CN 02148848 CN1203809C (en) 2002-11-22 2002-11-22 Measurement condition reproducing device and method based on body's surface texture characteristic and contact pressure

Publications (1)

Publication Number Publication Date
WO2004056262A1 true WO2004056262A1 (en) 2004-07-08

Family

ID=4751550

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2003/000813 WO2004056262A1 (en) 2002-11-22 2003-09-24 Apparatus and method for the reproduction of measuring conditions based on the body surface or surface veins and the contacting pressure

Country Status (3)

Country Link
CN (1) CN1203809C (en)
AU (1) AU2003272847A1 (en)
WO (1) WO2004056262A1 (en)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE602007010192D1 (en) * 2006-01-26 2010-12-16 Toshiba Kk X-ray CT device, as well as methods for phantom assignment
BR112012028410A2 (en) * 2010-05-08 2016-11-16 Univ California image registration system and image registration method
US8479581B2 (en) 2011-05-03 2013-07-09 General Electric Company Device and method for measuring pressure on wind turbine components
CN103892805B (en) * 2014-03-26 2016-09-07 太仓微芯电子科技有限公司 The compress control method of bio signal detection and system
CN104809357B (en) * 2015-05-21 2017-10-03 北京理工大学 The analysis method and device of a kind of multiple dimensioned touching act of three-dimensional rough surface
CN106360941B (en) * 2016-11-17 2019-08-13 广州杰赛科技股份有限公司 A kind of nail localization method and nail art device
CN108229255A (en) * 2016-12-21 2018-06-29 江苏迪康姆生物技术有限公司 A kind of palmmprint acquires all-in-one machine
CN107242855A (en) * 2017-06-05 2017-10-13 天津大学 A kind of biological tissue's dynamic modulation spectral measurement device and method
CN109984725B (en) * 2017-12-29 2022-07-12 天津先阳科技发展有限公司 Contact pressure interference suppression method and device in diffuse reflection measurement and measurement method
CN109269706A (en) * 2018-10-26 2019-01-25 东南大学 A kind of view-based access control model sensor legged type robot foot end multi-dimensional force sensing device
CN113925486A (en) * 2020-07-14 2022-01-14 先阳科技有限公司 Tissue component measurement method, tissue component measurement device, electronic apparatus, tissue component measurement system, and storage medium
CN114468994B (en) * 2021-02-11 2023-02-28 先阳科技有限公司 Tissue component measuring method and device and wearable equipment
CN115420208B (en) * 2022-11-04 2023-03-24 之江实验室 Texture sensor based on optical fiber knot sensitive structure and elastic shifting piece

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09108204A (en) * 1995-10-24 1997-04-28 Kdk Corp Measurement position reproduction method, measurement position reproduction device, and body fluid component concentration measuring device using it
CN1165556A (en) * 1995-08-07 1997-11-19 株式会社京都第一科学 Method and apparatus for measuring light projection concentration
CN1201379A (en) * 1995-10-31 1998-12-09 京都第一科学株式会社 Measuring condition setting jig, measuring condition setting method and biological information measuring instrument
CN1211735A (en) * 1997-08-22 1999-03-24 株式会社京都第一科学 Method of and apparatus for organism measurement
CN1224600A (en) * 1997-08-06 1999-08-04 仓敷纺织株式会社 Pressure regulating device and method of regulating pressure applied to living body during biodata measurement
CN1044528C (en) * 1996-09-05 1999-08-04 马开荣 Thin-film thermosensitive resistance infrared detector

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1165556A (en) * 1995-08-07 1997-11-19 株式会社京都第一科学 Method and apparatus for measuring light projection concentration
JPH09108204A (en) * 1995-10-24 1997-04-28 Kdk Corp Measurement position reproduction method, measurement position reproduction device, and body fluid component concentration measuring device using it
CN1201379A (en) * 1995-10-31 1998-12-09 京都第一科学株式会社 Measuring condition setting jig, measuring condition setting method and biological information measuring instrument
CN1044528C (en) * 1996-09-05 1999-08-04 马开荣 Thin-film thermosensitive resistance infrared detector
CN1224600A (en) * 1997-08-06 1999-08-04 仓敷纺织株式会社 Pressure regulating device and method of regulating pressure applied to living body during biodata measurement
CN1211735A (en) * 1997-08-22 1999-03-24 株式会社京都第一科学 Method of and apparatus for organism measurement

Also Published As

Publication number Publication date
CN1416780A (en) 2003-05-14
CN1203809C (en) 2005-06-01
AU2003272847A1 (en) 2004-07-14

Similar Documents

Publication Publication Date Title
WO2004056262A1 (en) Apparatus and method for the reproduction of measuring conditions based on the body surface or surface veins and the contacting pressure
JP3579686B2 (en) Measuring position reproducing method, measuring position reproducing device, and optical measuring device using the same
CN104132613B (en) Noncontact optical volume measurement method for complex-surface and irregular objects
CN101703394B (en) Radial pulse wave detection device
US7015473B2 (en) Method and apparatus for internal feature reconstruction
Scholz Reliability and validity of the WATSMART™ three-dimensional optoelectric motion analysis system
US20180058840A1 (en) Conformance Test Artifact for Coordinate Measuring Machine
US20120253165A1 (en) Intraocular pressure detecting device and detecting method thereof
JP2009521052A (en) Biological information detection using sweep-type imaging device
EP0766063A2 (en) Lens parameter measurement using optical sectioning
CN111076855A (en) Glass stress detection device
EP3317651A1 (en) Interferometric roll-off measurement using a static fringe pattern
CN217429994U (en) Nicotinic acid skin reaction image processing system
TW201017092A (en) Three-dimensional model reconstruction method and system thereof
JP2000283824A (en) Container information measuring apparatus
TWI687200B (en) Optical pulse image measuring device and method for analyzing change of pulse waveform
CN112858106A (en) System and method for measuring liquid drop contact angle
TWI352190B (en)
Duckworth et al. A clinically affordable non-contact wound measurement device
CN115227286A (en) Device and method for puncture self-adaptive guiding of ultrasonic probe and blood sampling box
KR102049040B1 (en) Apparatus for diagnosing biometric roughness
JPH05253210A (en) Wrinkle measuring system
JPH04138133A (en) Coil position detecting method and mr device
CN109758121A (en) A kind of non-intuitive detection device of manpower vascular surface and method
CN2681128Y (en) Screw detecting machine

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
122 Ep: pct application non-entry in european phase
NENP Non-entry into the national phase

Ref country code: JP

WWW Wipo information: withdrawn in national office

Country of ref document: JP