KR101191458B1 - Structure Deformation Measurement System and method - Google Patents

Structure Deformation Measurement System and method Download PDF

Info

Publication number
KR101191458B1
KR101191458B1 KR20100104837A KR20100104837A KR101191458B1 KR 101191458 B1 KR101191458 B1 KR 101191458B1 KR 20100104837 A KR20100104837 A KR 20100104837A KR 20100104837 A KR20100104837 A KR 20100104837A KR 101191458 B1 KR101191458 B1 KR 101191458B1
Authority
KR
South Korea
Prior art keywords
screen
laser
coordinate system
image
laser irradiator
Prior art date
Application number
KR20100104837A
Other languages
Korean (ko)
Other versions
KR20120043502A (en
Inventor
명현
전해민
Original Assignee
한국과학기술원
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 한국과학기술원 filed Critical 한국과학기술원
Priority to KR20100104837A priority Critical patent/KR101191458B1/en
Publication of KR20120043502A publication Critical patent/KR20120043502A/en
Application granted granted Critical
Publication of KR101191458B1 publication Critical patent/KR101191458B1/en

Links

Images

Abstract

Structure displacement measuring systems and methods are provided. The structure displacement measuring system controls the irradiation angle of the first laser irradiator so that the laser beam irradiated from the first laser irradiator is irradiated inside the second screen, and the laser beam irradiated from the second laser irradiator is inside the first screen. The irradiation angle of the second laser irradiator is controlled to be irradiated, and the relative displacement of the first structure with respect to the second structure using the first image generated by photographing the first screen and the second image generated by photographing the second screen. Estimate This makes it possible to adjust the position of the points inside the screen by controlling the direction of the irradiating lasers for the measurement of the structure displacement, so that the displacement measurement of the structure is always possible in any external situation.

Description

Structure Deformation Measurement System and Method

The present invention relates to a structure displacement measuring system and method, and more particularly to a system and method capable of measuring the displacement of a large structure.

Civil or architectural structures are exposed to external loads such as traffic, earthquakes, gusts, Therefore, in the case of civil engineering or building structure, reasonable and accurate design construction is important, but proper maintenance work is very important for maintaining the optimum usability of the structure and extending the life of the structure. In particular, with the recent increase in large structures such as high-rise buildings and long bridges, the importance of the Structural Health Monitoring System has increased. Such a structural safety diagnosis system can improve the stability of the structure by measuring, analyzing and diagnosing the dynamic behavior of the structures such as bridges and buildings. The structure safety diagnosis system requires many techniques such as structure damage identification method, data acquisition and transmission method, and the like.

Structural safety diagnostic systems typically use inclinometers, acceleration sensors, strain gauges, PZT sensors, etc. to measure structure displacement. Although the measurement of displacement is very important, it has not been studied much because of the huge size of the structure and the difficulty of access. One of the most commonly used methods of structural displacement measurement is the use of contact-type sensors that require a linear variable differential transformer (LVDT) and a stable reference point at the bottom. However, this type of contact sensor is not practical because it requires a reference point to be observed from the outside.

Displacement measurement using vision is also one of the active research areas. Most of the currently developed vision methods are to install a target on the structure to be measured and observe the movement of the target by using a high-performance camera at a remote reference point. This method also requires a reference point and the distance between the target and the camera is far. Because of this, it is difficult to measure due to the weather.

Accordingly, there is a growing need for vision and laser-based mobile robotic devices and methods that do not require a reference point and are robust to external environmental changes such as weather and illumination.

SUMMARY OF THE INVENTION The present invention has been made to solve the above problems, and an object of the present invention is to provide a structure displacement measuring system and method which can adjust the point positions inside the screen by controlling the direction of the laser beams for measuring the structure displacement. In providing.

Another object of the present invention is to provide a structure displacement measuring system and method capable of more accurate structure displacement measurement based on initial offset correction.

In accordance with an aspect of the present invention, a structure displacement measuring system includes: a first laser irradiator installed on a first structure and irradiating a laser beam to a second screen attached to the second structure; A second laser irradiator installed on the second structure and irradiating a laser beam to a first screen attached to the first structure; A first camera photographing the first screen to generate a first image; A second camera photographing the second screen to generate a second image; And controlling an irradiation angle of the first laser irradiator so that the laser beam irradiated from the first laser irradiator is irradiated inside the second screen, and the laser beam irradiated from the second laser irradiator is irradiated inside the first screen. And a control module for controlling the irradiation angle of the second laser irradiator to estimate the relative displacement of the first structure with respect to the second structure using the first image and the second image.

Preferably, the control module corrects the displacement estimated in the estimating step by using a calibration matrix for initial offset correction.

The control module detects a first screen image in the first image, a second screen image in the second image, and a plurality of laser points by the first laser irradiator in the detected second screen image. And a transformation matrix for detecting a laser point by the second laser irradiator from the detected first screen image, and converting a coordinate system of the first laser irradiator into a coordinate system of a first screen and a coordinate system of the first screen. Converts the coordinate system of the second screen into a coordinate system of the second screen, converts the coordinate system of the second laser irradiator into a coordinate system of the second screen, and converts the coordinate system of the second screen into the coordinate system of the second screen. A second equation of motion is generated by using a transformation matrix that is converted into a coordinate system of the first screen, and the second equation of motion It is preferable to estimate the displacement using the second equation of motion.

The control module may estimate the displacement by applying the first and second kinematic equations to an iterative calculation method such as Newton-Labson method or Extended Kalman Filtering (EKF) method.

The control module calculates a first laser irradiation angle at which the positional difference between the laser point formed on the second screen and the second target point is minimized by the first laser irradiator, and calculates the first laser irradiation angle according to the first laser irradiation angle. The laser irradiation angle of the laser irradiator is controlled, and the second laser irradiation angle is calculated by calculating a second laser irradiation angle that minimizes the position difference between the laser point formed on the first screen and the first target point by the second laser irradiator. It is preferable to control the laser irradiation angle of the said 2nd laser irradiator accordingly.

The first laser irradiator may irradiate two laser beams onto the second screen, and the second laser irradiator may irradiate one laser beam onto the first screen.

The control module controls the center of the two laser points by the first laser irradiator to be coupled to the second target point of the second screen, and the laser point by the second laser irradiator is controlled by the first screen of the first screen. It is preferable to control so as to be in one target point.

The first target point may be the center of the first screen, and the second target point may be the center of the second screen.

Meanwhile, the structure displacement measuring method according to the present invention includes a first irradiation step of irradiating a laser beam to a second screen attached to the second structure from a first laser irradiator installed on the first structure; A second irradiation step of irradiating a laser beam to a first screen attached to the first structure from a second laser irradiator installed on the second structure; A first control step of controlling an irradiation angle of the first laser irradiator such that the laser beam irradiated from the first laser irradiator is irradiated into the second screen; A second control step of controlling an irradiation angle of the second laser irradiator such that the laser beam irradiated from the second laser irradiator is irradiated into the first screen; A first generation step of photographing the first screen to generate a first image; A second generation step of generating a second image by photographing the second screen; And estimating a relative displacement of the first structure relative to the second structure using the first image and the second image.

The structure displacement measuring method may further include correcting the displacement estimated in the estimating step by using a calibration matrix for initial offset correction.

The estimating step may include detecting a first screen image from the first image;
Detecting a second screen image in the second image; Detecting a plurality of laser points by the first laser irradiator from the detected second screen image; Detecting a laser point by the second laser irradiator on the detected first screen image; Generating a first equation of motion using a transformation matrix for converting the coordinate system of the first laser irradiator into a coordinate system of the first screen and a transformation matrix for converting the coordinate system of the first screen to the coordinate system of the second screen; Generating a second equation of motion using a transformation matrix for converting the coordinate system of the second laser irradiator into a coordinate system of the second screen and a transformation matrix for converting the coordinate system of the second screen to the coordinate system of the first screen; And estimating the displacement by using the first and second motion equations.

The estimating of the displacement may include estimating the displacement by applying the first and second equations of motion to an iterative calculation method such as Newton-Labson method or Extended Kalman Filtering (EKF) method.

The first control step may include: calculating a first laser irradiation angle that minimizes a position difference between a laser point formed on the second screen and a second target point by the first laser irradiator; And controlling the laser irradiation angle of the first laser irradiator according to the calculated first laser irradiation angle, wherein the second control step comprises: a laser formed on the first screen by the second laser irradiator Calculating a second laser irradiation angle that minimizes the positional difference between the point and the first target point; And controlling the laser irradiation angle of the second laser irradiator according to the calculated second laser irradiation angle.

The first irradiating step irradiates two laser beams from the first laser irradiator to the second screen, and the second irradiating step irradiates one laser beam from the second laser irradiator to the first screen. can do.

The first control step may include controlling the center of two laser points by the first laser irradiator to the second target point of the second screen, and the second control step may be performed by the second laser irradiator. The laser point may be controlled to be coupled to the first target point of the first screen.

The first target point may be the center of the first screen, and the second target point may be the center of the second screen.

On the other hand, the structure displacement control system according to the present invention, is installed on the first structure, the first laser irradiator for irradiating a laser beam to a second screen attached to the second structure; A second laser irradiator installed on the second structure and irradiating a laser beam to a first screen attached to the first structure; A first camera photographing the first screen to generate a first image; A second camera photographing the second screen to generate a second image; And a transformation matrix for detecting a plurality of laser points by the first laser irradiator in the second image, and converting a coordinate system of the first laser irradiator into a coordinate system of a first screen, and converting a coordinate system of the first screen into the second system. A first equation of motion is generated using a transformation matrix that is converted into a coordinate system of a screen, a laser point is detected by the second laser irradiator in the first image, and the coordinate system of the second laser irradiator is converted into a coordinate system of a second screen. A second kinematic equation is generated using a transform matrix for converting to and a transform matrix for converting the coordinate system of the second screen to the coordinate system of the first screen, and using the first and second motion equations, And a control module for estimating the relative displacement of the first structure relative to the second structure.

Preferably, the control module corrects the displacement estimated in the estimating step by using a calibration matrix for initial offset correction.

On the other hand, the structure displacement control method according to the invention, the first irradiation step of irradiating a laser beam to the second screen attached to the second structure from the first laser irradiator installed on the first structure; A second irradiation step of irradiating a laser beam to a first screen attached to the first structure from a second laser irradiator installed on the second structure; A first generation step of photographing the first screen to generate a first image; A second generation step of generating a second image by photographing the second screen; Detecting a plurality of laser points by the first laser irradiator in the second image; Generating a first equation of motion using a transformation matrix for converting the coordinate system of the first laser irradiator into a coordinate system of the first screen and a transformation matrix for converting the coordinate system of the first screen to the coordinate system of the second screen; Detecting a laser point by the second laser irradiator in the first image; Generating a second equation of motion using a transformation matrix for converting the coordinate system of the second laser irradiator into a coordinate system of the second screen and a transformation matrix for converting the coordinate system of the second screen to the coordinate system of the first screen; And estimating a relative displacement of the first structure relative to the second structure using the first and second motion equations.

The structure displacement control method may further include correcting the displacement estimated in the estimating step by using a correction matrix for initial offset correction.

On the other hand, the computer-readable recording medium according to the present invention, the first irradiation step of irradiating a laser beam from the first laser irradiator installed on the first structure to the second screen attached to the second structure; A second irradiation step of irradiating a laser beam to a first screen attached to the first structure from a second laser irradiator installed on the second structure; A first control step of controlling an irradiation angle of the first laser irradiator such that the laser beam irradiated from the first laser irradiator is irradiated into the second screen; A second control step of controlling an irradiation angle of the second laser irradiator such that the laser beam irradiated from the second laser irradiator is irradiated into the first screen; A first generation step of photographing the first screen to generate a first image; A second generation step of generating a second image by photographing the second screen; And estimating a relative displacement of the first structure relative to the second structure using the first image, the irradiation angle of the first laser irradiator , the second image and the irradiation angle of the second laser irradiator; Contains programs that can run.

On the other hand, the computer-readable recording medium according to the present invention, the first irradiation step of irradiating a laser beam from the first laser irradiator installed on the first structure to the second screen attached to the second structure; A second irradiation step of irradiating a laser beam to a first screen attached to the first structure from a second laser irradiator installed on the second structure; A first generation step of photographing the first screen to generate a first image; A second generation step of generating a second image by photographing the second screen; Detecting a plurality of laser points by the first laser irradiator in the second image; Generating a first equation of motion using a transformation matrix for converting the coordinate system of the first laser irradiator into a coordinate system of the first screen and a transformation matrix for converting the coordinate system of the first screen to the coordinate system of the second screen; Detecting a laser point by the second laser irradiator in the first image; Generating a second equation of motion using a transformation matrix for converting the coordinate system of the second laser irradiator into a coordinate system of the second screen and a transformation matrix for converting the coordinate system of the second screen to the coordinate system of the first screen; And estimating a relative displacement of the first structure relative to the second structure using the first and second motion equations. Contains programs that can run.

As described above, according to the present invention, it is possible to adjust the position of the point to the inside of the screen by controlling the direction of the laser irradiation for the structure displacement measurement, it is always possible to measure the displacement of the structure in any external situation.

In addition, more accurate structural relative displacement measurements can be made based on initial offset calibration.

In addition, according to the present invention, the distance between the screen and the camera can be shortened, so that unlike the conventional vision method, the external environment is hardly affected.

1 is a view showing the concept of a laser stance fixed structure displacement measuring system,
2 is a detailed block diagram of a laser stance fixed structure displacement measuring system;
3 shows a simulation result showing the effect of the initial offset correction,
4 is a view illustrating a concept of a laser attitude controlled structure displacement measuring system;
5 shows the laser irradiator-A in detail;
6 is a detailed block diagram of a laser attitude controlled structure displacement measuring system;
7 to 11 are views showing the results of simulating the measurement performance of the structure displacement of the controlled structure displacement measurement system shown in FIG. 4, and
12 and 13 are diagrams showing simulation results capable of comparing the controlled structure displacement measuring system and the fixed structure displacement measuring system.

Hereinafter, with reference to the drawings will be described the present invention in more detail.

1. Laser attitude "fixed" structure displacement measuring system

1 is a view showing the concept of a laser stance fixed structure displacement measurement system. Laser attitude "Fixed" structure displacement measuring system is a system in which the direction of the laser irradiated on the screen for measuring the structure displacement is fixed. Hereinafter, for convenience of notation and description, it will be abbreviated as "fixed structure displacement measuring system".

As shown in FIG. 1, a fixed structure displacement measurement system includes a measurement module-A 110, a measurement module-B 120, and a calculation module (not shown). Measurement module-A 110 is installed in Structure-A (not shown) and measurement module-B 120 is installed in Structure-B (not shown).

Measurement module-A 110 includes a screen-A 111, a laser irradiator-A 113, and a camera-A 115, and measurement module-B 120 includes a screen-B 121, a laser irradiator. -B 123 and camera-B 125.

The laser irradiator-A 113 irradiates a pair of parallel laser beams to the screen-B 121. Accordingly, a pair of laser points are formed on the screen-B 121 attached to the structure-B. The camera-B 125 transfers the image-B generated by photographing the screen-B 121 to the calculation module.

The laser irradiator-B 123 irradiates one laser beam to the screen-A 111. Accordingly, one laser point is formed on the screen-A 111 attached to the structure-A. The camera-A 115 transfers the image-A generated by photographing the screen-A 111 to the calculation module.

The calculation module estimates the structure displacement by using the image-A received from the camera-A 115 and the image-B received from the camera-B 125. Here, "displacement structure" is a relative displacement, physically coordinate system of the screen -A (111) attached to the structure -A relative to the coordinate system of the screen Σ B -B (121) attached to the structure -B Σ A Means.

The structure displacement p can be expressed as [x, y, z, θ, φ, ψ]. here,

1) x is the linear displacement in the x direction from screen-B 121 to screen-A 111,

2) y is the linear displacement in the y direction from screen-B 121 to screen-A 111,

3) z is the linear displacement in the z direction from Screen-B 121 to Screen-A 111,

4) θ is the rotational displacement of the x coordinate axis of the screen-A (111) relative to the x coordinate axis of the screen-B 121,

5) φ is the rotational displacement of the y coordinate axis of the Screen-A (111) with respect to the y coordinate axis of the Screen-B 121, and

6) ψ is a rotational displacement of the z coordinate axis of Screen-A 111 with respect to the z coordinate axis of Screen-B 121.

Hereinafter, a calculation module for estimating the structure displacement p = [x, y, z, θ, φ, ψ] using the image-A and the image-B will be described in detail with reference to FIG. 2. 2 is a detailed block diagram of a computing module.

As illustrated in FIG. 2, the calculation module 200 includes an image corrector 210, a screen detector 220, a laser point detector 230, and a displacement estimator 240.

The image corrector 210 corrects and removes the distortion added to the image-A received from the camera-A 115. Similarly, for the image-B received from the camera-B 125, the image corrector 210 corrects and removes the distortion.

The screen detector 220 detects the Screen-A (111) image from the corrected image-A and detects the Screen-B (121) image from the corrected image-B by using a Corner Detection Algorithm. do.

The laser point detector 230 detects the position A O of the laser point projected from the detected screen-A 111 image. In addition, the laser point detector 230 detects positions B O and B Y of two laser points projected from the detected screen-B 121 image.

The displacement estimator 240 uses the positions A O, B O, and B Y of the laser points detected by the laser point detector 230 to displace the structure p = [x, y, z, θ, φ, ψ ] Is estimated.

2. Structure displacement estimation process

Hereinafter, the process of estimating the structure displacement p = [x, y, z, θ, φ, ψ] by the displacement estimator 240 will be described in detail.

Referring to FIG. 1, the geometric relationship between the laser point positions m = [ A O, B O, B Y] T and the structure displacement p = [x, y, z, θ, φ, ψ] T is shown. Kinematic Equation can be derived.

To this end, the coordinate system of the screen -B (121) Σ B requires the coordinate system of the screen from -A (111) Σ conversion matrix for converting the coordinates to A A T B, which is expressed as Equation 1 below.

Figure 112010069417013-pat00001

Here, T (x, y, z) is a translation matrix about the x, y, and z axes, R (x, θ) is a rotation matrix about the x-axis, and R (y, φ) is the rotation matrix about the y axis, and R (z, ψ) is the rotation matrix about the z axis.

The position A O of the laser point formed on the screen-A 111 may be calculated according to Equation 2 below using the transformation matrix A T B.

Figure 112010069417013-pat00002

Here, [0, L] means the x coordinate and the y coordinate of the laser irradiator-B 123 on the screen-B 121. And Z BA means the distance from the screen-B 121 to the screen-A (111).

In the same manner as in Equation 2, B O and B Y can also be calculated. On the other hand, since the screen-A 111 and the screen-B 121 are xy planes, the z coordinate is zero. Accordingly, the z coordinates of the laser points on screen-A 111 or screen-B 121 are zero.

If the calculated positions A 0, B O and B Y of the laser points are respectively divided with respect to the x coordinate and the y coordinate, a motion equation such as Equation 3 below may be derived.

Figure 112010069417013-pat00003

Here, s θ is sinθ, c θ is cosθ, A O x is the coordinate, A O of A O x y is the y-coordinate, B O x is the x-coordinate of B O, B O y of A O is y of B O coordinate, y a x is the x coordinate of y a, y a y represents the y coordinate of the a y.

By applying the Newton-Raphson method shown in Equation 4 below using the equation of motion M shown in Equation 3, displacement p = [x, y, z, θ, φ, ψ] It is possible to estimate.

Figure 112010069417013-pat00004

Where Jp (= ∂M / ∂p) is Jacobian (Jacoby function: Jacobian) of the equation M, J + p is pseudo-inverse of Jacobian,

Figure 112010069417013-pat00005
Denotes laser point positions estimated by M, and m (k) denotes actual observed laser point positions.

In the above, the Newton-Lapson technique is used to estimate the displacement p, but it is only one implementation, so it is possible to use another iterative calculation technique. For example, it is also possible to estimate the displacement p = [x, y, z, θ, φ, ψ] by applying the motion equation M shown in Equation 3 to the Extended Kalman Filtering (EKF) technique.

3. Initial Offset  correction( Calibration of initial offset )

Hereinafter, in the fixed structure displacement measuring system described above and the controlled structure displacement measuring system described below, a method of correcting offset due to initial displacement will be described in detail.

For the initial offset calibration, we first calculate the initial displacement p 0 = [x 0 , y 0 , z 0 , θ 0 , φ 0 , ψ 0 ] T and use it to generate a calibration matrix. Matrix A T cal is shown in Equation 5 below.

Figure 112010069417013-pat00006

Where T (x, y, z) is a translation matrix for the X, Y, and Z axes, R (x, θ) is a rotation matrix for the X axis, and R (y, [theta] is the rotation matrix about the Y axis, and R (z, ψ) is the rotation matrix about the z axis.

Thereafter, as shown in Equation 6 below, the inverse of the calibration matrix A T cal is multiplied by the transformation matrix A T B to obtain a corrected transformation matrix T new .

Figure 112010069417013-pat00007

Using the corrected transformation matrix T new, we can obtain the displacement p new whose initial offset is corrected, and the corrected displacement p new = [x ', y', z ', θ', φ ', ψ'] As shown in Equation 7, it can be expressed using the elements of T new .

Figure 112010069417013-pat00008

3 shows a simulation result showing the effect of the initial offset correction. In FIG. 3, the solid line is the displacement estimation result when the initial offset is corrected, and the dotted line is the displacement estimation result when the initial offset is not corrected.

As shown in FIG. 3, it can be seen that the solid line starts from 0, due to the correction of the initial offset.

4. Laser attitude "controlled" structure displacement measuring system

In the fixed structure displacement measuring system described above, laser points are required to be present in the screen as a prerequisite for estimating structure displacement. However, due to environmental impacts, it may be difficult for laser points to be irradiated into the screen, which needs to be compensated for.

Hereinafter, a structure displacement measuring system capable of adjusting point positions by controlling the irradiation direction of lasers by 2-DOF driving is proposed.

4 is a view illustrating a concept of a laser attitude controlled structure displacement measuring system. Laser attitude "Controlled" structure displacement measuring system is a system that can control the direction of the laser irradiated to the screen to measure the structure displacement, and ultimately adjust the position of the laser point on the screen. For convenience, abbreviated as "controlled structure displacement measuring system".

As shown in FIG. 4, the controlled structure displacement measurement system includes a measurement module-A 410, a measurement module-B 420, and a control module (not shown). Measurement module-A 410 is installed in Structure-A (not shown) and measurement module-B 420 is installed in Structure-B (not shown).

Measurement module-A 410 has a screen-A 411, a laser irradiator-A 413 and a camera-A 415, and measurement module-B 420 has a screen-B 421, a laser irradiator A B-423 and a camera-B 425.

The laser irradiator-A 413 irradiates a pair of parallel laser beams. The laser irradiator-A 413 is rotatable about the x and z axes. The left side of FIG. 5 shows a view of the laser irradiator-A 413 as viewed from above, and the right side of FIG. 5 shows a view of the laser irradiator-A 413 as viewed from the right.

The laser irradiator-A 413 is drive controlled by the control module. Specifically, the control module is arranged such that the center of the pair of laser points by the laser irradiator-A 413 is centered on the target point-B in the screen-B 421 attached to the structure-B. Control the laser irradiation angle.

The camera-B 425 transfers the image-B generated by photographing the screen-B 421 to the calculation module.

The laser irradiator-B 423 irradiates one laser beam. The laser irradiator-B 423 is also rotatable about the x and z axes.

The laser irradiator-B 423 is also drive controlled by the control module. Specifically, the control module is configured such that the laser irradiation angle of the laser irradiator-B 423 is such that the laser point by the laser irradiator-B 423 is coupled to the target point-A in the screen-B 421 attached to the structure-B. To control.

The camera-A 415 transfers the image-A generated by photographing the screen-A 411 to the calculation module.

The control module controls the laser irradiation angles of the laser irradiators 413 and 423 by using 'image-A received from the camera-A 415' and 'image-B received from the camera-B 425'. And estimate the displacement of the structure.

Hereinafter, a control module that controls the laser irradiation angle using the image-A and the image-B and estimates the structure displacement p = [x, y, z, θ, φ, ψ] in detail with reference to FIG. 6. Explain. 6 is a detailed block diagram of the control module.

As shown in FIG. 6, the control module 600 includes an image corrector 610, a screen detector 620, a laser point detector 630, a displacement estimator 640, a point position difference calculator 650, A laser irradiation angle calculator 660, a laser irradiation angle controller 670, a transformation matrix generator 680, and a calibration matrix generator 690 are provided.

The image corrector 610, the screen detector 620, and the laser point detector 630 may be implemented in the same manner as the image corrector 210, the screen detector 220, and the laser point detector 230 illustrated in FIG. 2. Therefore, detailed description thereof will be omitted.

The point position difference calculator 650 calculates the position difference between the center of the pair of laser points by the laser irradiator-A 413 and the target point-B. In addition, the point position difference calculator 650 calculates a position difference between the laser point and the target point-A by the laser irradiator-B 423.

Here, target point-A is the center of screen-A, and target point-B is the center of screen-B.

The laser irradiation angle calculation unit 660 is a laser irradiation of the laser irradiator-A (413) to minimize the positional difference between the 'center of a pair of laser points by the laser irradiator-A (413)' and 'target point-B' Calculate the angle. In addition, the laser irradiation angle calculator 660 calculates the laser irradiation angle of the laser irradiator-B 423 which minimizes the positional difference between the 'laser point by the laser irradiator-B 423' and the 'target point-A'. do.

The laser irradiation angle controller 670 controls the laser irradiator-A 413 and the laser irradiator-B 423, respectively, so that the laser irradiation is performed according to the laser irradiation angle calculated by the laser irradiation angle calculator 660.

A transformation matrix generating unit 680, a laser irradiator -A (413) coordinate system Σ A laser irradiator -B (423) generates a "to the coordinate system of the screen -A Σ A transformation matrix A T A to convert ', and for for the coordinate system Σ B generates "the transformation matrix B T B to convert to the coordinate system of the screen Σ B -B '.

Hereinafter, a process of generating the transformation matrix A T A and the transformation matrix B T B by the transformation matrix generator 680 will be described in detail.

Equation 8 below including the laser irradiation angles of the laser irradiator-A 413 can be assumed.

Figure 112010069417013-pat00009

Where B t is the center of screen-B 421 to target point-B. B t x is the x-coordinate of B t, B t y is the y coordinate of B t. Θ A and ψ A mean laser irradiation angles such that the center of the laser points by the laser irradiator-A 413 is located at the center of the screen-B 421 which is the target point-B.

In addition, [0, 0] means the x coordinate and the y coordinate of the laser irradiator-A 413 in the screen-A 411, and Z A'B is the laser irradiator-A 413 and the screen-B 421. ) Is the distance between.

Since the z coordinate is zero in the coordinate system Σ B of the screen-B 421, the z coordinate of the target point B t is zero. Meanwhile, in Equation 8, the transformation matrix B T A may be obtained by calculating an inverse of A T B shown in Equation 1.

Meanwhile, similarly to Equation 8, the following Equation 9 including the laser irradiation angles of the laser irradiator-B 423 may be assumed.

Figure 112010069417013-pat00010

Where A t is the center of screen-A 411 to target point-A. A t x t A is the x-coordinate, y is the y-coordinate of A t A t. Θ B and ψ B mean laser irradiation angles such that the laser point irradiated from the laser irradiator-B 423 is located at the center of the screen-A 411 which is the target point-A.

In addition, [0, -L] means the x coordinate and the y coordinate of the laser irradiator-B 423 in the screen-B 421, and Z B'A denotes the laser irradiator-B 423 and the screen-A ( 411) distance.

Since the z coordinate is zero in the coordinate system Σ A of the screen-A 411, the z coordinate of the target point A t is zero.

When the above Equations 8 and 9 are combined into one, the relationship between the target points and the laser irradiation angles may be expressed by Equation 10.

Figure 112010069417013-pat00011

here,

Figure 112010069417013-pat00012
to be. And, θ is a rotation angle in a coordinate system Σ AB to B of the screen -A (411) screen coordinate system Σ -B (421) at the time of A in the x-axis, θ B is to the x axis screen -B It is the rotation angle from the coordinate system (Sigma) B of (421) to the coordinate system (Sigma) B ' of the laser irradiator-B (423).

Meanwhile,

Figure 112010069417013-pat00013
Similarly,
Figure 112010069417013-pat00014
ego,
Figure 112010069417013-pat00015
Is,
Figure 112010069417013-pat00016
to be.

By applying the Newton-Rabson technique or EKF using Equation M r shown in Equation 10, it is possible to estimate the laser irradiation angles that can locate the laser point or the center of the laser points at the target point.

Coordinates of laser irradiation -A (413) Σ A 'screen -A (411) Σ coordinate system transformation matrix to convert A A T A of the' possible is calculated using the θ ψ A and A, equation (11) below, and same.

Figure 112010069417013-pat00017

Further, in the same way as above, the coordinate system of the applicator -B (423) using the θ and ψ B B. Σ B to calculate the "a transformation matrix B T B to convert to the coordinate system of the laser B Σ screen -B (421) ' Of course you can.

Meanwhile, the calibration matrix generator 690 generates a calibration matrix A T cal for initial offset calibration. The calibration matrix generator 690 generates a calibration matrix according to the above-described "3. initial offset calibration".

The displacement estimator 640 may include 1) positions of three laser points detected by the laser point detector 630, and 2) transformation matrices A T A , B T generated by the transformation matrix generator 680. B ′ 3) The structure displacement is estimated using the calibration matrix A T cal generated by the calibration matrix generator 690.

Hereinafter, a method of estimating structure displacement by the displacement estimating unit 640 will be described in detail.

First, using the positions of the three laser points detected by the laser point detector 630 and the transformation matrices A T A and B T B generated by the transformation matrix generator 680, the displacement p = The equation of motion including [x, y, z, θ, φ, ψ] is derived.

Specifically, the position A O B of the laser point irradiated from the laser irradiator-B 423 to the screen-A 411 is expressed by Equation 12 below.

Figure 112010069417013-pat00018

Further, the positions B O A and B Y A of the laser points irradiated from the laser irradiator-A 413 to the screen-B 421 are expressed by Equation 13 below.

Figure 112010069417013-pat00019

In Equations 12 and 13, the z coordinates of A O B , B O A and B Y A are zero.

Integrating equations (12) and (13) into one, a motion equation including the displacement p = [x, y, z, θ, φ, ψ] for the controlled structure displacement measurement system is derived as shown in equation (14) below. .

Figure 112010069417013-pat00020

Subsequently, the displacement estimator 640 repeatedly applies the Newton-Rapson method shown in Equation 4 using the equation M vs. in Equation 14, whereby displacement p = [x, y, z, θ, φ , ψ] is estimated. In addition, the displacement estimating unit 640 may estimate the displacement p = [x, y, z, θ, φ, ψ] by applying the motion equation M vs shown in Equation 14 to the EKF technique.

The displacement estimator 640 obtains the corrected transform matrix T new by substituting the correction matrix A T cal generated by the correction matrix generator 690 into Equation 6 described above, and corrects the corrected transform matrix T new. By using Equation 7, the displacement p new = [x ', y', z ', θ', φ ', ψ'] obtained by correcting the initial offset is obtained.

5. Simulation Results

7 to 11 are diagrams showing the results of simulating the measurement performance of the structure displacement of the controlled structure displacement measurement system shown in FIG.

7 to 11, the solid line represents the actual structure displacement, the long dotted line represents the estimated structure displacement when the initial offset calibration is performed, and the short dotted line represents the estimated structure displacement when the initial offset calibration is not performed. Indicates. Meanwhile, the distance between Structure-A and Structure-B was assumed to be 100m.

According to the simulation result using the Newton-Rapson method shown in FIG. 7, it can be seen that the estimated structure displacement converges to the actual structure displacement in a relatively short time when the initial offset correction is performed.

In addition, FIG. 8 shows a simulation result of applying the Newton-Rabson technique when noise is added artificially to the detected laser point positions. In this case, the results are relatively satisfactory.

On the other hand, Figure 9 shows the results of estimating the structure displacement by applying the EKF technique in the same conditions as in Figure 8, it can be seen that the excellent performance even when the EKF technique is applied.

In addition, FIG. 10 shows a simulation result of applying the Newton-Rabson technique when the x coordinate of the structure displacement is artificially changed. In this case, it can be seen that the convergence to the actual structure displacement is changed well.

On the other hand, Figure 11 shows the results of estimating the structure displacement by applying the EKF technique under the same conditions as in Figure 10, it can be seen that the excellent performance even when the EKF technique is applied.

The upper part of FIG. 12 shows the position coordinates of the radar points in the screen when the Newton-Lapson technique is applied in the controlled structure displacement measuring system, and the lower part of FIG. The position coordinates of the point are shown. 12 it can be seen that the displacement of the laser point is more stable in the controlled structure displacement measurement system.

In addition, FIG. 13 illustrates a case in which the EKF technique is applied under the same conditions as in FIG. 12, and the same result as in FIG. 12 can be obtained.

Meanwhile, the structure-A and the structure-B in which the measurement module-A and the measurement module-B mentioned in the above embodiments are installed may be separate structures as well as a single structure.

While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it is to be understood that the invention is not limited to the disclosed exemplary embodiments, but, on the contrary, It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the present invention.

110, 410: Measurement module-A 111, 411: Screen-A
113, 413: Laser irradiator-A 115, 415: Camera-A
120, 420: Measurement module-B 121, 421: Screen-B
123, 423: laser irradiator-B 125, 425: camera-B
200: operation module 210: image correction unit
220: screen detection unit 230: laser point detection unit
240: displacement estimation unit 600: control module
610: image correction unit 620: screen detection unit
630: laser point detection unit 640: displacement estimation unit
650: point position difference calculation unit 660: laser irradiation angle calculation unit
670: laser irradiation angle control unit 680: transformation matrix generation unit
690: calibration matrix generator

Claims (22)

A first laser irradiator installed on the first structure and irradiating a laser beam to a second screen attached to the second structure;
A second laser irradiator installed on the second structure and irradiating a laser beam to a first screen attached to the first structure;
A first camera photographing the first screen to generate a first image;
A second camera photographing the second screen to generate a second image; And
The irradiation angle of the first laser irradiator is controlled so that the laser beam irradiated from the first laser irradiator is irradiated inside the second screen, and the laser beam irradiated from the second laser irradiator is irradiated inside the first screen. And a control module for controlling an irradiation angle of the second laser irradiator and estimating a relative displacement of the first structure with respect to the second structure using the first image and the second image. Structure displacement measuring system.
The method of claim 1,
The control module,
And correcting the estimated displacement using a calibration matrix for initial offset calibration.
The method of claim 1,
The control module,
Detecting a first screen image in the first image, detecting a second screen image in the second image, detecting a plurality of laser points by the first laser irradiator in the detected second screen image,
Detecting a laser point by the second laser irradiator on the detected first screen image,
Generating a first equation of motion using a transformation matrix for converting the coordinate system of the first laser irradiator into a coordinate system of the first screen and a transformation matrix for converting the coordinate system of the first screen to the coordinate system of the second screen,
Generating a second equation of motion using a transformation matrix for converting the coordinate system of the second laser irradiator into a coordinate system of the second screen and a transformation matrix for converting the coordinate system of the second screen to the coordinate system of the first screen,
And using the first and second motion equations to estimate the displacement.
The method of claim 3, wherein
The control module,
And applying the first and second kinematic equations to an iterative calculation technique of one of Newton-Rapson and Extended Kalman Filtering (EKF) techniques to estimate the displacement.
The method of claim 1,
The control module,
The first laser irradiation angle is calculated by the first laser irradiator to minimize the position difference between the laser point formed on the second screen and the second target point, the laser irradiation of the first laser irradiator according to the first laser irradiation angle Control the angle,
The second laser irradiation angle is calculated by the second laser irradiator to minimize the position difference between the laser point formed on the first screen and the first target point, and the laser irradiation of the second laser irradiator according to the second laser irradiation angle Structure displacement measuring system, characterized in that for controlling the angle.
The method of claim 1,
The first laser irradiator,
Irradiate two laser beams onto the second screen,
The second laser irradiator,
And irradiating one laser beam onto the first screen.
The method according to claim 6,
The control module,
Control the center of the two laser points by the first laser irradiator to be joined to the second target point of the second screen,
And control the laser point by the second laser irradiator to the first target point of the first screen.
8. The method of claim 7,
The first target point is the center of the first screen,
And the second target point is the center of the second screen.
A first irradiation step of irradiating a laser beam from a first laser irradiator installed on the first structure to a second screen attached to the second structure;
A second irradiation step of irradiating a laser beam to a first screen attached to the first structure from a second laser irradiator installed on the second structure;
A first control step of controlling an irradiation angle of the first laser irradiator such that the laser beam irradiated from the first laser irradiator is irradiated into the second screen;
A second control step of controlling an irradiation angle of the second laser irradiator such that the laser beam irradiated from the second laser irradiator is irradiated into the first screen;
A first generation step of photographing the first screen to generate a first image;
A second generation step of generating a second image by photographing the second screen; And
Estimating a relative displacement of the first structure with respect to the second structure using the first image and the second image.
The method of claim 9,
And correcting the displacement estimated in the estimating step by using a calibration matrix for initial offset correction.
The method of claim 9,
The estimating step,
Detecting a first screen image in the first image;
Detecting a second screen image in the second image;
Detecting a plurality of laser points by the first laser irradiator from the detected second screen image;
Detecting a laser point by the second laser irradiator on the detected first screen image;
Generating a first equation of motion using a transformation matrix for converting the coordinate system of the first laser irradiator into a coordinate system of the first screen and a transformation matrix for converting the coordinate system of the first screen to the coordinate system of the second screen;
Generating a second equation of motion using a transformation matrix for converting the coordinate system of the second laser irradiator into a coordinate system of the second screen and a transformation matrix for converting the coordinate system of the second screen to the coordinate system of the first screen; And
Estimating the displacement by using the first equation and the second equation of motion.
12. The method of claim 11,
Estimating the displacement,
And applying the first and second equations of motion to an iterative calculation method of any one of Newton-Labson and Extended Kalman Filtering (EKF) techniques to estimate the displacement.
The method of claim 9,
The first control step,
Calculating a first laser irradiation angle at which a position difference between a laser point formed on the second screen and a second target point is minimized by the first laser irradiator; And
Controlling the laser irradiation angle of the first laser irradiator according to the calculated first laser irradiation angle;
The second control step,
Calculating a second laser irradiation angle by which the position difference between the laser point formed on the first screen and the first target point is minimized by the second laser irradiator; And
And controlling the laser irradiation angle of the second laser irradiator in accordance with the calculated second laser irradiation angle.
The method of claim 9,
The first irradiation step,
Irradiating two laser beams from the first laser irradiator to the second screen,
The second irradiation step,
And irradiating one laser beam from the second laser irradiator onto the first screen.
The method of claim 14,
The first control step,
Control the center of the two laser points by the first laser irradiator to be joined to the second target point of the second screen,
The second control step,
And controlling the laser point by the second laser irradiator to the first target point of the first screen.
16. The method of claim 15,
The first target point is the center of the first screen,
And the second target point is the center of the second screen.
A first laser irradiator installed on the first structure and irradiating a laser beam to a second screen attached to the second structure;
A second laser irradiator installed on the second structure and irradiating a laser beam to a first screen attached to the first structure;
A first camera photographing the first screen to generate a first image;
A second camera photographing the second screen to generate a second image; And
A transformation matrix for detecting a plurality of laser points by the first laser irradiator in the second image, and converting a coordinate system of the first laser irradiator into a coordinate system of a first screen, and converting a coordinate system of the first screen to the second screen A first equation of motion is generated using a transformation matrix that transforms the coordinate system into a coordinate system, a laser point is detected by the second laser irradiator in the first image, and the coordinate system of the second laser irradiator is used as the coordinate system of the second screen. A second motion equation is generated using a transform matrix for transforming and a transform matrix for transforming a coordinate system of the second screen into a coordinate system of the first screen. The second motion equation is generated using the first motion equation and the second motion equation. And a control module for estimating the relative displacement of the first structure relative to the second structure. Above measurement system.
18. The method of claim 17,
The control module,
And correcting the estimated displacement using a calibration matrix for initial offset calibration.
A first irradiation step of irradiating a laser beam from a first laser irradiator installed on the first structure to a second screen attached to the second structure;
A second irradiation step of irradiating a laser beam to a first screen attached to the first structure from a second laser irradiator installed on the second structure;
A first generation step of photographing the first screen to generate a first image;
A second generation step of generating a second image by photographing the second screen;
Detecting a plurality of laser points by the first laser irradiator in the second image;
Generating a first equation of motion using a transformation matrix for converting the coordinate system of the first laser irradiator into a coordinate system of the first screen and a transformation matrix for converting the coordinate system of the first screen to the coordinate system of the second screen;
Detecting a laser point by the second laser irradiator in the first image;
Generating a second equation of motion using a transformation matrix for converting the coordinate system of the second laser irradiator into a coordinate system of the second screen and a transformation matrix for converting the coordinate system of the second screen to the coordinate system of the first screen; And
Estimating a relative displacement of the first structure relative to the second structure using the first and second motion equations.
20. The method of claim 19,
And correcting the displacement estimated in the estimating step by using a calibration matrix for initial offset correction.
A first irradiation step of irradiating a laser beam from a first laser irradiator installed on the first structure to a second screen attached to the second structure;
A second irradiation step of irradiating a laser beam to a first screen attached to the first structure from a second laser irradiator installed on the second structure;
A first control step of controlling an irradiation angle of the first laser irradiator such that the laser beam irradiated from the first laser irradiator is irradiated into the second screen;
A second control step of controlling an irradiation angle of the second laser irradiator such that the laser beam irradiated from the second laser irradiator is irradiated into the first screen;
A first generation step of photographing the first screen to generate a first image;
A second generation step of generating a second image by photographing the second screen; And
Estimating a relative displacement of the first structure relative to the second structure by using the first image and the second image; and a computer-readable recording medium containing a program capable of performing the same.
A first irradiation step of irradiating a laser beam from a first laser irradiator installed on the first structure to a second screen attached to the second structure;
A second irradiation step of irradiating a laser beam to a first screen attached to the first structure from a second laser irradiator installed on the second structure;
A first generation step of photographing the first screen to generate a first image;
A second generation step of generating a second image by photographing the second screen;
Detecting a plurality of laser points by the first laser irradiator in the second image;
Generating a first equation of motion using a transformation matrix for converting the coordinate system of the first laser irradiator into a coordinate system of the first screen and a transformation matrix for converting the coordinate system of the first screen to the coordinate system of the second screen;
Detecting a laser point by the second laser irradiator in the first image;
Generating a second equation of motion using a transformation matrix for converting the coordinate system of the second laser irradiator into a coordinate system of the second screen and a transformation matrix for converting the coordinate system of the second screen to the coordinate system of the first screen; And
Estimating a relative displacement of the first structure with respect to the second structure using the first and second motion equations.
KR20100104837A 2010-10-26 2010-10-26 Structure Deformation Measurement System and method KR101191458B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR20100104837A KR101191458B1 (en) 2010-10-26 2010-10-26 Structure Deformation Measurement System and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR20100104837A KR101191458B1 (en) 2010-10-26 2010-10-26 Structure Deformation Measurement System and method

Publications (2)

Publication Number Publication Date
KR20120043502A KR20120043502A (en) 2012-05-04
KR101191458B1 true KR101191458B1 (en) 2012-10-18

Family

ID=46263669

Family Applications (1)

Application Number Title Priority Date Filing Date
KR20100104837A KR101191458B1 (en) 2010-10-26 2010-10-26 Structure Deformation Measurement System and method

Country Status (1)

Country Link
KR (1) KR101191458B1 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101247966B1 (en) 2012-10-23 2013-04-02 주식회사 동우기술단 Method of safety diagnostic with telescope-camera for structure
KR101646981B1 (en) 2015-02-27 2016-08-09 한국건설기술연구원 Data processing system for evaluating safety of structure, and method for the same
KR101914942B1 (en) * 2018-10-08 2018-11-06 정홍석 Apparatus and method for measuring alignment state of the roll
KR101917619B1 (en) 2018-02-23 2018-11-13 (주)카이센 System for measuring bridge deflection
KR101930878B1 (en) * 2018-06-22 2018-12-19 정홍석 Apparatus and method for measuring alignment state of the roll

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8937725B2 (en) * 2012-06-14 2015-01-20 Nikon Corporation Measurement assembly including a metrology system and a pointer that directs the metrology system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003232615A (en) 2002-02-08 2003-08-22 Jfe Engineering Kk Method and device for measuring position of measurement point, method of measuring curved shape, and screen means
JP2009092535A (en) 2007-10-10 2009-04-30 Ono Sokki Co Ltd Optical displacement gauge
KR20100103928A (en) 2009-03-16 2010-09-29 한국과학기술원 Structural displacement measurement system and method using thereof

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003232615A (en) 2002-02-08 2003-08-22 Jfe Engineering Kk Method and device for measuring position of measurement point, method of measuring curved shape, and screen means
JP2009092535A (en) 2007-10-10 2009-04-30 Ono Sokki Co Ltd Optical displacement gauge
KR20100103928A (en) 2009-03-16 2010-09-29 한국과학기술원 Structural displacement measurement system and method using thereof

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101247966B1 (en) 2012-10-23 2013-04-02 주식회사 동우기술단 Method of safety diagnostic with telescope-camera for structure
KR101646981B1 (en) 2015-02-27 2016-08-09 한국건설기술연구원 Data processing system for evaluating safety of structure, and method for the same
KR101917619B1 (en) 2018-02-23 2018-11-13 (주)카이센 System for measuring bridge deflection
KR101930878B1 (en) * 2018-06-22 2018-12-19 정홍석 Apparatus and method for measuring alignment state of the roll
KR101914942B1 (en) * 2018-10-08 2018-11-06 정홍석 Apparatus and method for measuring alignment state of the roll
WO2020075955A1 (en) * 2018-10-08 2020-04-16 주식회사 피씨엠솔루션 Apparatus for detecting relative positioning information between rolls, and method for measuring roll alignment state by using same
US11192158B2 (en) 2018-10-08 2021-12-07 Pcm Solution Co., Ltd. Apparatus for detecting relative positioning information between rolls, and method for measurement roll alignment state by using same

Also Published As

Publication number Publication date
KR20120043502A (en) 2012-05-04

Similar Documents

Publication Publication Date Title
KR101217186B1 (en) Displacement measurement system of structure and method of the same
KR101191458B1 (en) Structure Deformation Measurement System and method
US9989973B2 (en) Measurement control system for multi-shaft supported air floatation platform
Du et al. Online robot calibration based on hybrid sensors using Kalman Filters
Myung et al. Paired structured light for structural health monitoring robot system
CN111207774A (en) Method and system for laser-IMU external reference calibration
WO2020073929A1 (en) Operation trolley, and robotic arm optical target positioning device and positioning method thereof
JP2017209781A5 (en)
JP2006322937A (en) Determination method of 3d coordinates of object surface
US20180194010A1 (en) Mechanism-parameter-calibration method for robotic arm system
JP2022046691A (en) State checking device for built-in object, operation checking device and state checking method for built-in object
JP4302830B2 (en) Robot calibration method and apparatus
Park et al. Deformation monitoring of a building structure using a motion capture system
EP3567340A1 (en) Visual inspection arrangement
KR101560578B1 (en) Apparatus and method for controling direction error of gimbal apparatus using image processing
KR101573641B1 (en) 6-DOF Displacement Measurement System and Method
KR101029751B1 (en) Structural displacement measurement system and method using thereof
KR101340555B1 (en) Apparatus and method for generating base view image
Lu et al. Six-axis position measurement system for levitated motion stages
Hartmann et al. High accurate pointwise (geo-) referencing of a k-tls based multi-sensor-system
KR101283932B1 (en) Method for measuring direction error of gimbal platform and apparatus thereof
Isa et al. Kinematic error analysis of stage tracking using stereo vision
KR20210157747A (en) Displacement Measurement System of structure using image processing unit and laser
KR101436097B1 (en) Non-Contacting Method for Measuring 6-DOF Motion Based on Laser Sensor
Zöller et al. Preview of Driving Surface Unevenness in Wheeled Mobile Driving Simulators

Legal Events

Date Code Title Description
A201 Request for examination
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20161004

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20170928

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20181002

Year of fee payment: 7