CN105043355A - Side slope micro-deformation monitoring method and side slope micro-deformation monitoring system based on similarity determination criterion - Google Patents

Side slope micro-deformation monitoring method and side slope micro-deformation monitoring system based on similarity determination criterion Download PDF

Info

Publication number
CN105043355A
CN105043355A CN201510242170.5A CN201510242170A CN105043355A CN 105043355 A CN105043355 A CN 105043355A CN 201510242170 A CN201510242170 A CN 201510242170A CN 105043355 A CN105043355 A CN 105043355A
Authority
CN
China
Prior art keywords
moment
image
monitored area
side slope
analyzed
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510242170.5A
Other languages
Chinese (zh)
Other versions
CN105043355B (en
Inventor
赵谦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Dishan Shiju Technology Co ltd
Original Assignee
Xian University of Science and Technology
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 Xian University of Science and Technology filed Critical Xian University of Science and Technology
Priority to CN201510242170.5A priority Critical patent/CN105043355B/en
Publication of CN105043355A publication Critical patent/CN105043355A/en
Application granted granted Critical
Publication of CN105043355B publication Critical patent/CN105043355B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • G01C11/04—Interpretation of pictures
    • G01C11/06—Interpretation of pictures by comparison of two or more pictures of the same area
    • G01C11/12—Interpretation of pictures by comparison of two or more pictures of the same area the pictures being supported in the same relative position as when they were taken
    • G01C11/14—Interpretation of pictures by comparison of two or more pictures of the same area the pictures being supported in the same relative position as when they were taken with optical projection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Multimedia (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

The invention discloses a side slope micro-deformation monitoring method and a side slope micro-deformation monitoring system based on similarity determination criterion, wherein the monitoring method includes the following steps: (1) image collection and synchronous uploading: collecting and uploading side slope surface images in a monitored zone, wherein the monitored zone is divided into K to-be-analyzed zones; (2) imaging receiving and synchronous processing: performing synchronous treatment to the received side slope surface images, including: (201) performing side slope surface image treatment at the initial sampling time and calculating the similarity of the side slope surface image at this time; (202) performing side slope surface image treatment at the next sampling time, and calculating the similarity of the side slope surface image at this time; and (203) determining whether the micro-deformation exists or not. The monitoring system includes K digital projectors, K image photographing devices, an image collection device and an image processing device. The method and the system are reasonable in design, are convenient to achieve, are good in use effect, and can be used for performing simple, real-time and accurate monitoring the side slope micro-deformation.

Description

Based on side slope Light deformation monitoring method and the monitoring system of similarity decision criteria
Technical field
The invention belongs to micro-change monitoring technical field, especially relate to a kind of side slope Light deformation monitoring method based on similarity decision criteria and monitoring system.
Background technology
The side slope surface monitoring of Light deformation and Realtime Alerts most important, that carries out monitoring as Slope of Open Pit Coal Mine surface Light deformation micro-ly becomes the major issue that monitoring and Realtime Alerts concern miner's life security.When reality is monitored side slope surface Light deformation, mainly the rock mass fission produced because of sedimentation is monitored.At present, instrument or method for detecting side slope surface Light deformation have: total powerstation, convergence instrument, optical fiber sensing technology etc., the each tool advantage of above-mentioned detection technique, but also existing defects and deficiency all to some extent, as: when adopting optical fiber sensing technology to carry out the detection of side slope Light deformation, be very difficult by embedded fibers to side slope inside, and the sensitivity detected need carry out modeling analysis according to side slope attribute; Total powerstation and convergence instrument is adopted to carry out side slope Light deformation when detecting, extensive work is done before measurement, comprise and observation station is set, erection prism device etc., and can only gather the three-dimensional data of an observation station in measuring process, need pointwise to gather the distance of side slope surface to instrument, when gathering the irregular slope data of one side, workload is larger at every turn, need take a long time, and be difficult to accomplish real-time, intellectual monitoring.
In recent years, image processing techniques and 3D reconstruction technique achieve fast development.In image processing field, empirical mode decomposition (EMD) method is a kind of signal analysis method with self-adaptation time-frequency resolution capability, and two-dimensional empirical mode decomposition (BEMD) method be the further genralrlization of empirical mode decomposition (EMD) method and its to be a kind of disobeying be disinclined to the adaptive approach of data-driven of basis function, two-dimensional empirical mode decomposition (BEMD) method was also widely used in Image Multiscale analysis in recent years.
Because measuring method has noncontact, precision advantages of higher, therefore it is widely used in three-dimensional measurement field.In three-dimensional measurement field, dynamic, real-time three-dimensional measurement is the difficult point of research.At present, the optical three-dimensional measuring method be widely used has structured light projection method, stereo vision method, laser scanning method, laser interferance method, time-of-flight method etc.Wherein, laser scanning method, due to consuming time longer, is generally difficult to use in dynamic object scanning; Laser interferance method and confocal are to environmental requirement harshness, and equipment complex and expensive, is generally used for the precision measurement of specific industry; And be applicable to the 3D measuring method of dynamic realtime monitoring, mainly contain ToF, holographic imaging method, Structure light method and a few class of stereo vision method.
The defect existed for existing side slope Light deformation monitoring method and deficiency, if 3D method for reconstructing is combined with image processing method, be applied to the monitoring of side slope Light deformation, the monitoring difficulty that the monitoring equipment that just effectively can solve existing side slope Light deformation monitoring method and correspondence exists is large, time-consuming, the problems such as Real-Time Monitoring, monitoring effect be poor can not be realized.In addition, also there is the inaccurate defect of monitoring result in the monitoring equipment of existing side slope Light deformation monitoring method and correspondence, as side slope lastblock rock because of loosen and drop time, because now Light deformation does not appear in slope table face, report to the police if now carry out, probably produce false-alarm.
To sum up, nowadays lack a kind of reasonable in design, realize convenient and the side slope Light deformation monitoring method based on similarity decision criteria that result of use is good and monitoring system, easy, real-time, accurate measurements can be carried out to side slope Light deformation.
Summary of the invention
Technical matters to be solved by this invention is for above-mentioned deficiency of the prior art, a kind of side slope Light deformation monitoring method based on similarity decision criteria is provided, its method step is simple, reasonable in design and realization is convenient, result of use is good, can carry out easy, real-time, accurate measurements to side slope Light deformation.
For solving the problems of the technologies described above, the technical solution used in the present invention is: a kind of side slope Light deformation monitoring method based on similarity decision criteria, is characterized in that the method comprises the following steps:
Step one, image acquisition and synchronized upload: adopt image capture device and according to the sample frequency f preset, the side slope surface image of monitored area is gathered, and the side slope each sampling instant gathered surface image synchronization is sent to image processing equipment;
Described monitored area is region side slope surface needing carry out Light deformation monitoring, and institute monitored area is divided into K region to be analyzed, and the individual described region to be analyzed of K is square and its size is all identical; Wherein, K is positive integer and K >=2; Described side slope surface image is made up of K area image, and K described area image is respectively the image in K described region to be analyzed;
The sampling instant of described image capture device, is denoted as t h, t h=t 1+ (h-1) × Δ t; Wherein, h be positive integer and h=1,2,3 ..., t 1for the initial samples moment of described image capture device, Δ t is time interval of adjacent double sampling before and after described image capture device and Δ t=1/f;
Step 2, image-receptive and synchronously process: after described image processing equipment receives the side slope surface image that each sampling instant of image capture device gathers, process received side slope surface image synchronization, process is as follows:
The side slope surface image procossing that step 201, initial samples moment gather: adopt the side slope surface image of image processing equipment to the image capture device collection of initial samples moment to process, comprise the following steps:
Step 2011, side slope surface image store: store the side slope surface image be currently received, and the sampling instant t to this side slope surface image 1carry out record;
Step 2012, elevation information obtain: call image processing module and process the side slope surface image received in step 2011, obtain the altitude figures of each pixel on described side slope surface image, and on the described side slope surface image obtained, the altitude figures of each pixel is t 1the elevation information of monitored area described in the moment;
The elevation information of described monitored area comprises the elevation information in K described region to be analyzed, the altitude figures of each pixel on the area image that the elevation information in each described region to be analyzed includes this region to be analyzed;
T 1described in moment, the elevation information of monitored area comprises t 1the elevation information in moment K described region to be analyzed; Wherein, t 1moment is numbered the elevation information in the region to be analyzed of k, is denoted as X k(t 1); K is the numbering in described region to be analyzed, k is positive integer and k=1,2 ..., K;
Step 2013, elevation information store: to the t obtained in step 2012 1the elevation information in moment K described region to be analyzed, carries out stores synchronized;
Step 2014, Similarity Measure: according to formula (1), t is calculated 1the image similarity X of monitored area described in the moment 1,2 " K(t 1), wherein X 1,2 " K(t 1) ∈ (0,1]; In formula (1), x k(t 1) be t 1moment is numbered the elevation mean value in the region to be analyzed of k;
The side slope surface image procossing that step 202, next sampling instant gather: adopt the side slope surface image of image processing equipment to the collection of next sampling instant image capture device to process, comprise the following steps:
Step 2021, side slope surface image store: store the side slope surface image be currently received, and the sampling time t to this side slope surface image dcarry out record; Wherein, D is positive integer and D>=2;
Step 2022, elevation information obtain: according to the method described in step 2012, obtain t dthe elevation information of monitored area described in the moment;
T ddescribed in moment, the elevation information of monitored area comprises t dthe elevation information in moment K described region to be analyzed; Wherein, t dmoment is numbered the elevation information in the region to be analyzed of k, is denoted as X k(t d);
Step 2023, elevation information store: to the t obtained in step 2022 dthe elevation information in moment K described region to be analyzed, carries out stores synchronized;
Step 2024, Similarity Measure: according to formula (2), t is calculated dthe image similarity X of monitored area described in the moment 1,2 " K(t d), wherein X 1,2 " K(t d) ∈ (0,1]; In formula (2), x k(t d) be t dmoment is numbered the elevation mean value in the region to be analyzed of k;
Step 203, Light deformation judge: by the t calculated in step 2024 dthe image similarity X of monitored area described in the moment 1,2 " K(t d), respectively with t dbefore moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 " K(t d) and t dbefore moment monitored area described in each moment image similarity between difference when being all not more than ε, illustrate from t 1moment is to t ddescribed in moment there is not Light deformation in monitored area, returns step 202 afterwards, processes the side slope surface image that next sampling instant gathers; Otherwise, t is described dthere is Light deformation in monitored area described in the moment, described image processing equipment controls alarm unit and reports to the police, and to the t of Light deformation occurs dmoment carries out record;
Wherein, ε is the Light deformation alarm threshold value and ε=0.01 ~ 0.1 that preset.
The above-mentioned side slope Light deformation monitoring method based on similarity decision criteria, it is characterized in that: when carrying out elevation information acquisition in step 2012 and step 2022, call image processing module and all according to three frequency color fringe projection method for three-dimensional measurement, obtain the altitude figures of each pixel on described side slope surface image; The side slope surface image gathered in step one is the distortion color fringe image of described monitored area;
Before carrying out image acquisition and synchronized upload in step one, first adopt image processing equipment to generate color fringe image, RGB tri-Color Channels of this color fringe image are generated by sine streak that is low, that neutralize high three kinds of carrier frequency respectively; Afterwards, then be projected on described monitored area by described color fringe image through RGB tri-color channels of digital projector simultaneously; Then, image capture apparatus is adopted to take the distortion color fringe image of described monitored area;
When carrying out elevation information acquisition in step 2012 and step 2022, process is as follows:
Step I, background subtraction and color decoupling zero: subtract each other in the distortion color fringe image of image capture device collection containing the Color Channel of high frequency fringes and the Color Channel containing intermediate frequency striped, obtain the bar graph of high and low frequency compound, decompose with Bidimensional Empirical Mode Decomposition BEMD again, be separated high, middle carrier component; In like manner, by containing in described distortion color fringe image, the Color Channel of low frequency component subtracts each other, in obtaining, the bar graph of low frequency compound, then with Bidimensional Empirical Mode Decomposition BEMD decompose obtain in, low carrier component;
Step I i, phase demodulating: with two-dimentional short time discrete Fourier transform demodulation, the height obtained, neutralize low each carrier component wrapped phase;
Step I ii, phase unwrapping: with change precision go parcel algorithm by low, in and high carrier frequency component complete wrapped phase successively and launch, obtain the expansion phase place of high frequency carrier frequency item, launch the height of monitored area described in phase recovery thus.
The above-mentioned side slope Light deformation monitoring method based on similarity decision criteria, is characterized in that: image processing equipment described in step 203 controls after alarm unit reports to the police, to call elevation information comparison module, to from t 1moment is to t dthe elevation mean value of moment K described analyzed area to be analyzed compares respectively, and according to comparative result, draws t dthe region to be analyzed of moment generation Light deformation; Afterwards, described image processing equipment exports t dthe numbering in the region to be analyzed of moment generation Light deformation;
Wherein, from t 1moment is to t dthe comparative approach of the elevation mean value in moment K described region to be analyzed is all identical; To from t 1moment is to t dwhen the elevation mean value in any one region to be analyzed of moment compares, by t dthe elevation mean value in this region to be analyzed of moment respectively with t 1moment is to t d-1in moment, the elevation mean value in this region to be analyzed of each moment carries out difference comparsion.
The above-mentioned side slope Light deformation monitoring method based on similarity decision criteria, is characterized in that: before the control of image processing equipment described in step 203 alarm unit is reported to the police, and described image processing equipment also needs to call pseudo-distortion judge module, to t dwhether the Light deformation that monitored area described in the moment occurs is that pseudo-distortion judges, process is as follows:
Step I, continuously monitoring: according to the method described in step 2021 to step 2024, to t d+1moment is to t d+Ethe side slope surface image of moment image capture device collection processes, and calculates the image similarity X of monitored area described in each moment 1,2 " K(t f); Wherein, F is positive integer and F=D+1 ~ D+E; Wherein, t d+E-t d+1=Δ T, Δ T are the continuous monitoring time preset; E be positive integer and
Step II, pseudo-distortion judge: by t d+1moment is to t d+Ethe image similarity X of monitored area described in each moment in moment 1,2 " K(t f), respectively with t dthe image similarity X of monitored area described in each moment before moment 1,2 " K(t h) compare: work as X 1,2 " K(t f) and t dthe image similarity X of monitored area described in each moment before moment 1,2 " K(t h) between difference when being all greater than ε, t is described dthe Light deformation that monitored area described in moment occurs is true distortion, and now described image processing equipment control alarm unit is reported to the police; Otherwise, t is described dthe Light deformation that monitored area described in moment occurs is pseudo-distortion.
The above-mentioned side slope Light deformation monitoring method based on similarity decision criteria, is characterized in that: Δ T=2s ~ 3min in step I.
The above-mentioned side slope Light deformation monitoring method based on similarity decision criteria, is characterized in that: the control of image processing equipment described in step 203 alarm unit enters step 3 after reporting to the police;
Continue to monitor after step 3, first generation Light deformation, comprise the following steps:
Step 301, next sampling instant image procossing: according to the method described in step 2021 to step 2024, to t mthe side slope surface image of moment image capture device collection processes, and calculates t mthe image similarity X of monitored area described in the moment 1,2 " K(t m); Wherein, M is positive integer and M>=D+1;
Step 302, Light deformation judge: by the t calculated in step 301 mthe image similarity X of monitored area described in the moment 1,2 " K(t m), respectively with t dmoment is to t min moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 " K(t m) and t dmoment is to t min moment monitored area described in each moment image similarity between difference when being all not more than ε, illustrate from t dmoment is to t mdescribed in moment there is not Light deformation in monitored area, returns step 301 afterwards, processes the side slope surface image that next sampling instant gathers; Otherwise, t is described mthere is Light deformation in monitored area described in the moment, described image processing equipment controls alarm unit and reports to the police, and to the t of Light deformation occurs mmoment and from t 1moment is to t mthe alarm times of alarm unit described in the moment carries out record, enters step 4 afterwards;
Continue to monitor after step 4, next time generation Light deformation, comprise the following steps:
Step 401, next sampling instant image procossing: according to the method described in step 2021 to step 2024, to t nthe side slope surface image of moment image capture device collection processes, and calculates t nthe image similarity X of monitored area described in the moment 1,2 " K(t n); Wherein, N is positive integer and N>=n, and wherein n is positive integer and t nmoment is t nthe generation moment of the last Light deformation in described monitored area before moment;
Step 402, Light deformation judge: by the t calculated in step 401 nthe image similarity X of monitored area described in the moment 1,2 " K(t n), respectively with t nmoment is to t nin moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 " K(t m) and t nmoment is to t nin moment monitored area described in each moment image similarity between difference when being all not more than ε, illustrate from t nmoment is to t ndescribed in moment there is not Light deformation in monitored area, returns step 401 afterwards, processes the side slope surface image that next sampling instant gathers; Otherwise, t is described nthere is Light deformation in monitored area described in the moment, described image processing equipment controls alarm unit and reports to the police, and to the t of Light deformation occurs nmoment and from t 1moment is to t nthe alarm times of alarm unit described in the moment carries out record, enters step 5 afterwards;
Step 5, one or many repeat step 4, continue to monitor described monitored area.
The above-mentioned side slope Light deformation monitoring method based on similarity decision criteria, is characterized in that: before the control of image processing equipment described in step 302 alarm unit is reported to the police, and described image processing equipment also needs to call pseudo-distortion judge module, to t mwhether the Light deformation that monitored area described in the moment occurs is that pseudo-distortion judges, process is as follows:
Step I-11, continuously monitoring: according to the method described in step 2021 to step 2024, to t m+1moment is to t m+Ethe side slope surface image of moment image capture device collection processes, and calculates the image similarity X of monitored area described in each moment 1,2 " K(t f1); Wherein, F1 is positive integer and F1=M+1 ~ M+E; Wherein, t m+E-t m+1=Δ T;
Step II-12, pseudo-distortion judge: by t m+1moment is to t m+Ethe image similarity X of monitored area described in each moment in moment 1,2 " K(t f1), respectively with t dmoment is to t min moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 " K(t f1) and t dmoment is to t min moment monitored area described in each moment image similarity between difference when being all greater than ε, t is described mthe Light deformation that monitored area described in moment occurs is true distortion, and now described image processing equipment control alarm unit is reported to the police; Otherwise, t is described mthe Light deformation that monitored area described in moment occurs is pseudo-distortion;
Before the control of image processing equipment described in step 402 alarm unit is reported to the police, described image processing equipment also needs to call pseudo-distortion judge module, to t nwhether the Light deformation that monitored area described in the moment occurs is that pseudo-distortion judges, process is as follows:
Step I-21, continuously monitoring: according to the method described in step 2021 to step 2024, to t n+1moment is to t n+Ethe side slope surface image of moment image capture device collection processes, and calculates the image similarity X of monitored area described in each moment 1,2 " K(t f2); Wherein, F2 is positive integer and F2=N+1 ~ N+E; Wherein, t m+E-t m+1=Δ T;
Step II-22, pseudo-distortion judge: by t n+1moment is to t n+Ethe image similarity X of monitored area described in each moment in moment 1,2 " K(t f2), respectively with t nmoment is to t nin moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 " K(t f2) and t nmoment is to t nin moment monitored area described in each moment image similarity between difference when being all greater than ε, t is described nthe Light deformation that monitored area described in moment occurs is true distortion, and now described image processing equipment control alarm unit is reported to the police; Otherwise, t is described nthe Light deformation that monitored area described in moment occurs is pseudo-distortion.
The above-mentioned side slope Light deformation monitoring method based on similarity decision criteria, is characterized in that: image processing equipment described in step 302 controls after alarm unit reports to the police, to call elevation information comparison module, to from t dmoment is to t mthe elevation mean value of moment K described analyzed area to be analyzed compares respectively, and according to comparative result, draws t mthe region to be analyzed of moment generation Light deformation; Afterwards, described image processing equipment exports t mthe numbering in the region to be analyzed of moment generation Light deformation;
Wherein, from t dmoment is to t mthe comparative approach of the elevation mean value in moment K described region to be analyzed is all identical; To from t dmoment is to t mwhen the elevation mean value in any one region to be analyzed of moment compares, by t mthe elevation mean value in this region to be analyzed of moment respectively with t dmoment is to t m-1in moment, the elevation mean value in this region to be analyzed of each moment carries out difference comparsion;
Image processing equipment described in step 402 controls after alarm unit reports to the police, to call elevation information comparison module, to from t nmoment is to t nthe elevation mean value of moment K described analyzed area to be analyzed compares respectively, and according to comparative result, draws t nthe region to be analyzed of moment generation Light deformation; Afterwards, described image processing equipment exports t nthe numbering in the region to be analyzed of moment generation Light deformation;
Wherein, from t nmoment is to t nthe comparative approach of the elevation mean value in moment K described region to be analyzed is all identical; To from t nmoment is to t nwhen the elevation mean value in any one region to be analyzed of moment compares, by t nthe elevation mean value in this region to be analyzed of moment respectively with t nmoment is to t n-1in moment, the elevation mean value in this region to be analyzed of each moment carries out difference comparsion.
Simultaneously, it is simple that the present invention also discloses a kind of structure, reasonable in design, input cost is low and use easy and simple to handle, the side slope Light deformation monitoring system based on similarity decision criteria that result of use is good, it is characterized in that: comprise the digital projector that K is laid in K described overlying regions to be analyzed respectively, K respectively to the image capture apparatus that the image in K described region to be analyzed absorbs, the image capture device that the side slope surface image of described monitored area is gathered and the image processing equipment connected with image capture device, the individual described image capture apparatus of K is CCD camera and it all connects with image capture device, K described digital projector all connects with image processing equipment, K described CCD camera is laid in K described overlying regions to be analyzed respectively.
Said system, is characterized in that: described digital projector and CCD camera lay respectively at above the both sides in described region to be analyzed.
The present invention compared with prior art has the following advantages:
1, the side slope Light deformation monitoring method step adopted is simple, reasonable in design and realization is convenient.
2, adopt three frequently color fringe projection method for three-dimensional measurement obtain the altitude figures of each pixel on side slope surface images, obtain that the precision of altitude figures is high and image processing speed is fast.Meanwhile, image-capture field is large, can according to actual needs, and in conjunction with the size of monitored area, the quantity K treating analyzed area adjusts accordingly, and practical operation is very easy.Thus, the present invention, by monitored area being divided into K region to be analyzed, under the prerequisite ensureing monitoring accuracy, can meet the monitoring requirements of different size monitored area.
3, adopt similarity decision criteria to carry out the monitoring of side slope Light deformation, data processing amount is little and data processing speed fast, and monitoring accuracy is high, and carry out precise monitoring to the micro-change process of side slope, measuring accuracy can reach about 0.1mm; Further, monitoring velocity is fast, can complete the side slope Light deformation observation process in each moment in 2 seconds.
4, result of use is good and monitoring result is accurate, can identify as side slope puppet distortion when side slope lastblock rock drops because loosening, avoid false-alarm, thus have the advantages such as monitoring result is accurate, reliability is high, can realize real-time, dynamic monitoring, popularizing application prospect is extensive.
5, the side slope Light deformation monitoring system structure adopted is simple, reasonable in design, equipment installation is laid conveniently and it is easy and simple to handle to use, result of use good, and input cost is lower.
In sum, the present invention is reasonable in design, realization is convenient and result of use is good, can carry out easy, real-time, accurate measurements to side slope Light deformation.
Below by drawings and Examples, technical scheme of the present invention is described in further detail.
Accompanying drawing explanation
Fig. 1 is the method flow block diagram that employing the present invention carries out when side slope Light deformation is monitored.
Fig. 2 is the schematic block circuit diagram of side slope Light deformation monitoring system of the present invention.
Description of reference numerals:
1-image processing equipment; 2-image capture device; 3-digital projector;
4-CCD camera; 5-alarm unit.
Embodiment
A kind of side slope Light deformation monitoring method based on similarity decision criteria as shown in Figure 1, comprises the following steps:
Step one, image acquisition and synchronized upload: adopt image capture device 2 and according to the sample frequency f preset, the side slope surface image of monitored area is gathered, and the side slope each sampling instant gathered surface image synchronization is sent to image processing equipment 1.
Described monitored area is region side slope surface needing carry out Light deformation monitoring, and institute monitored area is divided into K region to be analyzed, and the individual described region to be analyzed of K is square and its size is all identical; Wherein, K is positive integer and K >=2; Described side slope surface image is made up of K area image, and K described area image is respectively the image in K described region to be analyzed.
The sampling instant of described image capture device, is denoted as t h, t h=t 1+ (h-1) × Δ t; Wherein, h be positive integer and h=1,2,3 ..., t 1for the initial samples moment of described image capture device, Δ t is time interval of adjacent double sampling before and after described image capture device and Δ t=1/f.
Step 2, image-receptive and synchronously process: after described image processing equipment 1 receives the side slope surface image that each sampling instant of image capture device 2 gathers, process received side slope surface image synchronization, process is as follows:
The side slope surface image procossing that step 201, initial samples moment gather: the side slope surface image adopting image processing equipment 1 pair of initial samples moment image capture device 2 to gather processes, and comprises the following steps:
Step 2011, side slope surface image store: store the side slope surface image be currently received, and the sampling instant t to this side slope surface image 1carry out record;
Step 2012, elevation information obtain: call image processing module and process the side slope surface image received in step 2011, obtain the altitude figures of each pixel on described side slope surface image, and on the described side slope surface image obtained, the altitude figures of each pixel is t 1the elevation information of monitored area described in the moment;
The elevation information of described monitored area comprises the elevation information in K described region to be analyzed, the altitude figures of each pixel on the area image that the elevation information in each described region to be analyzed includes this region to be analyzed;
T 1described in moment, the elevation information of monitored area comprises t 1the elevation information in moment K described region to be analyzed; Wherein, t 1moment is numbered the elevation information in the region to be analyzed of k, is denoted as X k(t 1); K is the numbering in described region to be analyzed, k is positive integer and k=1,2 ..., K;
Step 2013, elevation information store: to the t obtained in step 2012 1the elevation information in moment K described region to be analyzed, carries out stores synchronized;
Step 2014, Similarity Measure: according to formula (1), t is calculated 1the image similarity X of monitored area described in the moment 1,2 " K(t 1), wherein X 1,2 " K(t 1) ∈ (0,1]; In formula (1), x k(t 1) be t 1moment is numbered the elevation mean value in the region to be analyzed of k;
The side slope surface image procossing that step 202, next sampling instant gather: adopt image processing equipment 1 to process the side slope surface image that next sampling instant image capture device 2 gathers, comprise the following steps:
Step 2021, side slope surface image store: store the side slope surface image be currently received, and the sampling time t to this side slope surface image dcarry out record; Wherein, D is positive integer and D>=2;
Step 2022, elevation information obtain: according to the method described in step 2012, obtain t dthe elevation information of monitored area described in the moment;
T ddescribed in moment, the elevation information of monitored area comprises t dthe elevation information in moment K described region to be analyzed; Wherein, t dmoment is numbered the elevation information in the region to be analyzed of k, is denoted as X k(t d);
Step 2023, elevation information store: to the t obtained in step 2022 dthe elevation information in moment K described region to be analyzed, carries out stores synchronized;
Step 2024, Similarity Measure: according to formula (2), t is calculated dthe image similarity X of monitored area described in the moment 1,2 " K(t d), wherein X 1,2 " K(t d) ∈ (0,1]; In formula (2), x k(t d) be t dmoment is numbered the elevation mean value in the region to be analyzed of k;
Step 203, Light deformation judge: by the t calculated in step 2024 dthe image similarity X of monitored area described in the moment 1,2 " K(t d), respectively with t dbefore moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 " K(t d) and t dbefore moment monitored area described in each moment image similarity between difference when being all not more than ε, illustrate from t 1moment is to t ddescribed in moment there is not Light deformation in monitored area, returns step 202 afterwards, processes the side slope surface image that next sampling instant gathers; Otherwise, t is described dthere is Light deformation in monitored area described in the moment, described image processing equipment 1 controls alarm unit 5 and reports to the police, and to the t of Light deformation occurs dmoment carries out record;
Wherein, ε is the Light deformation alarm threshold value and ε=0.01 ~ 0.1 that preset.
During actual use, can according to specific needs, the value of ε be adjusted accordingly.
In formula (1), min{x 1(t 1), x 2(t 1) ... x k(t 1) represent from x 1(t 1), x 2(t 1) ..., x k(t 1) in get minimum value, max{x 1(t 1), x 2(t 1) ... x k(t 1) represent from x 1(t 1), x 2(t 1) ..., x k(t 1) in get maximal value; In formula (2), min{x 1(t d), x 2(t d) ... x k(t d) represent from x 1(t d), x 2(t d) ..., x k(t d) in get minimum value, max{x 1(t d), x 2(t d) ... x k(t d) represent from x 1(t d), x 2(t d) ..., x k(t d) in get maximal value.
In the present embodiment, x in step 2014 k(t 1) be t 1the altitude figures mean value of all pixels on the area image that moment is numbered the region to be analyzed of k, x in step 2024 k(t d) be t dthe altitude figures mean value of all pixels on the area image that moment is numbered the region to be analyzed of k.
In the present embodiment, when carrying out elevation information acquisition in step 2012 and step 2022, call image processing module and all according to three frequency color fringe projection method for three-dimensional measurement, obtain the altitude figures of each pixel on described side slope surface image; The side slope surface image gathered in step one is the distortion color fringe image of described monitored area.
Before carrying out image acquisition and synchronized upload in step one, first adopt image processing equipment 1 to generate color fringe image, RGB tri-Color Channels of this color fringe image are generated by sine streak that is low, that neutralize high three kinds of carrier frequency respectively; Afterwards, then be projected on described monitored area by described color fringe image through RGB tri-color channels of digital projector 3 simultaneously; Then, image capture apparatus is adopted to take the distortion color fringe image of described monitored area;
When carrying out elevation information acquisition in step 2012 and step 2022, process is as follows:
Step I, background subtraction and color decoupling zero: subtract each other containing the Color Channel of high frequency fringes and the Color Channel containing intermediate frequency striped in the distortion color fringe image that image capture device 2 is gathered, obtain the bar graph of high and low frequency compound, decompose with Bidimensional Empirical Mode Decomposition BEMD again, be separated high, middle carrier component; In like manner, by containing in described distortion color fringe image, the Color Channel of low frequency component subtracts each other, in obtaining, the bar graph of low frequency compound, then with Bidimensional Empirical Mode Decomposition BEMD decompose obtain in, low carrier component;
Step I i, phase demodulating: with two-dimentional short time discrete Fourier transform demodulation, the height obtained, neutralize low each carrier component wrapped phase;
Step I ii, phase unwrapping: with change precision go parcel algorithm by low, in and high carrier frequency component complete wrapped phase successively and launch, obtain the expansion phase place of high frequency carrier frequency item, launch the height of monitored area described in phase recovery thus.
Thus, the present invention adopts three frequency color fringe projection method for three-dimensional measurement to obtain the altitude figures of described monitored area, this method for three-dimensional measurement is based on the three frequency color fringe projection technology of profiling of BEMD, reach the object of the reliable rapid deployment of high precision solution mediation wrapped phase of the lower fringe phase of single image shooting, the texture information of testee can also be recovered simultaneously, only need the expansion phase place accurately obtaining high frequency fringes carrier frequency item, just can recover testee high accuracy three-dimensional profile information.RGB tri-color channels of described digital projector 3 are projected to after on described monitored area simultaneously, and described image capture apparatus (specifically coloured image pickup apparatus) is from the colored deforming stripe figure (i.e. described distortion color fringe image) of another angle shot through described monitored area high modulation.Wherein, background subtraction, color decoupling zero and phase unwrapping are the gordian techniquies of this method for three-dimensional measurement, are the principal elements affecting precision and reliability.Adopt the three-dimensional sensing method that color fringe projection and stereoscopic vision merge, to measure DYNAMIC COMPLEX object.Three frequency color fringe projection method for three-dimensional measurement are a kind of two-way fusion mechanism, namely first binocular stereo vision global registration is instructed by the PHASE DISTRIBUTION of three frequency color fringes, conversely, the meticulous matching and correlation initial phase distortion of recycling binocular fringe gray level figure local multiple features, this distortion derives from reflectivity noise, shape is suddenlyd change and the non-sinusoidal of striped, final stereoscopic vision solution room point position coordinates, reduction testee exact height information.Actual when using, this method for three-dimensional measurement can the elevation information of monitored area described in easy, quick obtaining, and degree of accuracy is high.
During actual use, other method for three-dimensional measurement based on image procossing also can be adopted to obtain the altitude figures of described monitored area.
In the present embodiment, adopt three frequently color fringe projection method for three-dimensional measurement publication number disclosed on 01 18th, 2012 is CN102322822A, the patent No. is for 201110225971.2 and denomination of invention is method for three-dimensional measurement disclosed in the application for a patent for invention file of " a kind of three frequency color fringe projection method for three-dimensional measurement ".
Wherein, similarity decision criteria refers to and carries out the monitoring of side slope Light deformation according to the size of image similarity and change situation.Regard the elevation information of monitored area described in any instant as a fuzzy set, according to fuzzy mathematics theory, the close degree between 2 fuzzy sets can be weighed by similarity (also claiming approach degree).
In the present embodiment, the span of the image similarity of the described monitored area calculated in step 2014 and step 2024 is (0,1], when the value of image similarity is approximately 1, represent that tested side slope surface (i.e. described monitored area) is approximate smooth, all pixels now on described monitored area have maximum similarity; When image similarity value close to 0 time, represent that tested slope table mask has minimum similarity.
In the present embodiment, image processing equipment 1 described in step 203 controls after alarm unit 5 reports to the police, to enter step 3;
Continue to monitor after step 3, first generation Light deformation, comprise the following steps:
Step 301, next sampling instant image procossing: according to the method described in step 2021 to step 2024, to t mthe side slope surface image that moment image capture device 2 gathers processes, and calculates t mthe image similarity X of monitored area described in the moment 1,2 " K(t m); Wherein, M is positive integer and M>=D+1;
Step 302, Light deformation judge: by the t calculated in step 301 mthe image similarity X of monitored area described in the moment 1,2 " K(t m), respectively with t dmoment is to t min moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 " K(t m) and t dmoment is to t min moment monitored area described in each moment image similarity between difference when being all not more than ε, illustrate from t dmoment is to t mdescribed in moment there is not Light deformation in monitored area, returns step 301 afterwards, processes the side slope surface image that next sampling instant gathers; Otherwise, t is described mthere is Light deformation in monitored area described in the moment, described image processing equipment 1 controls alarm unit 5 and reports to the police, and to the t of Light deformation occurs mmoment and from t 1moment is to t mdescribed in moment, the alarm times of alarm unit 5 carries out record, enters step 4 afterwards;
Continue to monitor after step 4, next time generation Light deformation, comprise the following steps:
Step 401, next sampling instant image procossing: according to the method described in step 2021 to step 2024, to t nthe side slope surface image that moment image capture device 2 gathers processes, and calculates t nthe image similarity X of monitored area described in the moment 1,2 " K(t n); Wherein, N is positive integer and N>=n, and wherein n is positive integer and t nmoment is t nthe generation moment of the last Light deformation in described monitored area before moment;
Step 402, Light deformation judge: by the t calculated in step 401 nthe image similarity X of monitored area described in the moment 1,2 " K(t n), respectively with t nmoment is to t nin moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 " K(t m) and t nmoment is to t nin moment monitored area described in each moment image similarity between difference when being all not more than ε, illustrate from t nmoment is to t ndescribed in moment there is not Light deformation in monitored area, returns step 401 afterwards, processes the side slope surface image that next sampling instant gathers; Otherwise, t is described nthere is Light deformation in monitored area described in the moment, described image processing equipment 1 controls alarm unit 5 and reports to the police, and to the t of Light deformation occurs nmoment and from t 1moment is to t ndescribed in moment, the alarm times of alarm unit 5 carries out record, enters step 5 afterwards;
Step 5, one or many repeat step 4, continue to monitor described monitored area.
In the present embodiment, image processing equipment 1 described in step 203 controls before alarm unit 5 reports to the police, and described image processing equipment 1 also needs to call pseudo-distortion judge module, to t dwhether the Light deformation that monitored area described in the moment occurs is that pseudo-distortion judges, process is as follows:
Step I, continuously monitoring: according to the method described in step 2021 to step 2024, to t d+1moment is to t d+Ethe side slope surface image that moment image capture device 2 gathers processes, and calculates the image similarity X of monitored area described in each moment 1,2 " K(t f); Wherein, F is positive integer and F=D+1 ~ D+E; Wherein, t d+E-t d+1=Δ T, Δ T are the continuous monitoring time preset; E be positive integer and
Step II, pseudo-distortion judge: by t d+1moment is to t d+Ethe image similarity X of monitored area described in each moment in moment 1,2 " K(t f), respectively with t dthe image similarity X of monitored area described in each moment before moment 1,2 " K(t h) compare: work as X 1,2 " K(t f) and t dthe image similarity X of monitored area described in each moment before moment 1,2 " K(t h) between difference when being all greater than ε, t is described dthe Light deformation that monitored area described in moment occurs is true distortion, and now described image processing equipment 1 controls alarm unit 5 and reports to the police; Otherwise, t is described dthe Light deformation that monitored area described in moment occurs is pseudo-distortion.
In the present embodiment, Δ T=2s ~ 3min in step I.Formula in, represent and round downwards.
During actual use, can according to specific needs, the value size of Δ T be adjusted accordingly.
In the present embodiment, image processing equipment 1 described in step 302 controls before alarm unit 5 reports to the police, and described image processing equipment 1 also needs to call pseudo-distortion judge module, to t mwhether the Light deformation that monitored area described in the moment occurs is that pseudo-distortion judges, process is as follows:
Step I-11, continuously monitoring: according to the method described in step 2021 to step 2024, to t m+1moment is to t m+Ethe side slope surface image that moment image capture device 2 gathers processes, and calculates the image similarity X of monitored area described in each moment 1,2 " K(t f1); Wherein, F1 is positive integer and F1=M+1 ~ M+E; Wherein, t m+E-t m+1=Δ T;
Step II-12, pseudo-distortion judge: by t m+1moment is to t m+Ethe image similarity X of monitored area described in each moment in moment 1,2 " K(t f1), respectively with t dmoment is to t min moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 " K(t f1) and t dmoment is to t min moment monitored area described in each moment image similarity between difference when being all greater than ε, t is described mthe Light deformation that monitored area described in moment occurs is true distortion, and now described image processing equipment 1 controls alarm unit 5 and reports to the police; Otherwise, t is described mthe Light deformation that monitored area described in moment occurs is pseudo-distortion;
Image processing equipment 1 described in step 402 controls before alarm unit 5 reports to the police, and described image processing equipment 1 also needs to call pseudo-distortion judge module, to t nwhether the Light deformation that monitored area described in the moment occurs is that pseudo-distortion judges, process is as follows:
Step I-21, continuously monitoring: according to the method described in step 2021 to step 2024, to t n+1moment is to t n+Ethe side slope surface image that moment image capture device 2 gathers processes, and calculates the image similarity X of monitored area described in each moment 1,2 " K(t f2); Wherein, F2 is positive integer and F2=N+1 ~ N+E; Wherein, t m+E-t m+1=Δ T;
Step II-22, pseudo-distortion judge: by t n+1moment is to t n+Ethe image similarity X of monitored area described in each moment in moment 1,2 " K(t f2), respectively with t nmoment is to t nin moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 " K(t f2) and t nmoment is to t nin moment monitored area described in each moment image similarity between difference when being all greater than ε, t is described nthe Light deformation that monitored area described in moment occurs is true distortion, and now described image processing equipment 1 controls alarm unit 5 and reports to the police; Otherwise, t is described nthe Light deformation that monitored area described in moment occurs is pseudo-distortion.
In the present embodiment, image processing equipment 1 described in step 203 controls after alarm unit 5 reports to the police, to call elevation information comparison module, to from t 1moment is to t dthe elevation mean value of moment K described analyzed area to be analyzed compares respectively, and according to comparative result, draws t dthe region to be analyzed of moment generation Light deformation; Afterwards, described image processing equipment 1 exports t dthe numbering in the region to be analyzed of moment generation Light deformation;
Wherein, from t 1moment is to t dthe comparative approach of the elevation mean value in moment K described region to be analyzed is all identical; To from t 1moment is to t dwhen the elevation mean value in any one region to be analyzed of moment compares, by t dthe elevation mean value in this region to be analyzed of moment respectively with t 1moment is to t d-1in moment, the elevation mean value in this region to be analyzed of each moment carries out difference comparsion.
For realizing conveniently, preset elevation change alarm threshold value Δ H, to from t 1moment is to t dwhen the elevation mean value in any one region to be analyzed of moment compares, work as t dthe elevation mean value in this region to be analyzed of moment and t 1moment is to t d-1in moment this region to be analyzed of each moment elevation mean value between difference when being all not more than Δ H, illustrate that Light deformation does not occur in this region to be analyzed; Otherwise, illustrate that Light deformation occurs in this region to be analyzed.
Wherein, elevation change alarm threshold value Δ H, according to actual monitoring demand, people is for setting.
Further, image processing equipment 1 described in step 302 controls after alarm unit 5 reports to the police, to call elevation information comparison module, to from t dmoment is to t mthe elevation mean value of moment K described analyzed area to be analyzed compares respectively, and according to comparative result, draws t mthe region to be analyzed of moment generation Light deformation; Afterwards, described image processing equipment 1 exports t mthe numbering in the region to be analyzed of moment generation Light deformation;
Wherein, from t dmoment is to t mthe comparative approach of the elevation mean value in moment K described region to be analyzed is all identical; To from t dmoment is to t mwhen the elevation mean value in any one region to be analyzed of moment compares, by t mthe elevation mean value in this region to be analyzed of moment respectively with t dmoment is to t m-1in moment, the elevation mean value in this region to be analyzed of each moment carries out difference comparsion.
In the present embodiment, to from t dmoment is to t mwhen the elevation mean value in any one region to be analyzed of moment compares, work as t mthe elevation mean value in this region to be analyzed of moment and t dmoment is to t m-1in moment this region to be analyzed of each moment elevation mean value between difference when being all not more than Δ H, illustrate that Light deformation does not occur in this region to be analyzed; Otherwise, illustrate that Light deformation occurs in this region to be analyzed.
Image processing equipment 1 described in step 402 controls after alarm unit 5 reports to the police, to call elevation information comparison module, to from t nmoment is to t nthe elevation mean value of moment K described analyzed area to be analyzed compares respectively, and according to comparative result, draws t nthe region to be analyzed of moment generation Light deformation; Afterwards, described image processing equipment 1 exports t nthe numbering in the region to be analyzed of moment generation Light deformation;
Wherein, from t nmoment is to t nthe comparative approach of the elevation mean value in moment K described region to be analyzed is all identical; To from t nmoment is to t nwhen the elevation mean value in any one region to be analyzed of moment compares, by t nthe elevation mean value in this region to be analyzed of moment respectively with t nmoment is to t n-1in moment, the elevation mean value in this region to be analyzed of each moment carries out difference comparsion.
In the present embodiment, to from t nmoment is to t nwhen the elevation mean value in any one region to be analyzed of moment compares, work as t nthe elevation mean value in this region to be analyzed of moment and t nmoment is to t n-1in moment this region to be analyzed of each moment elevation mean value between difference when being all not more than Δ H, illustrate that Light deformation does not occur in this region to be analyzed; Otherwise, illustrate that Light deformation occurs in this region to be analyzed.
A kind of side slope Light deformation monitoring system based on similarity decision criteria as shown in Figure 2, comprise the digital projector 3 that K is laid in K described overlying regions to be analyzed respectively, K respectively to the image capture apparatus that the image in K described region to be analyzed absorbs, the image capture device 2 that the side slope surface image of described monitored area is gathered and the image processing equipment 1 connected with image capture device 2, the individual described image capture apparatus of K is CCD camera 4 and it all connects with image capture device 2, K described digital projector 3 all connects with image processing equipment 1, K described CCD camera 4 is laid in K described overlying regions to be analyzed respectively.
In the present embodiment, described digital projector 3 and CCD camera 4 lay respectively at above the both sides in described region to be analyzed.
Actual when layings installation is carried out to digital projector 3 and CCD camera 4, above side digital projector 3 being laid in described region to be analyzed, and above opposite side CCD camera 4 being laid in described region to be analyzed.
The above; it is only preferred embodiment of the present invention; not the present invention is imposed any restrictions, every above embodiment is done according to the technology of the present invention essence any simple modification, change and equivalent structure change, all still belong in the protection domain of technical solution of the present invention.

Claims (10)

1., based on a side slope Light deformation monitoring method for similarity decision criteria, it is characterized in that the method comprises the following steps:
Step one, image acquisition and synchronized upload: adopt image capture device (2) and according to the sample frequency f preset, the side slope surface image of monitored area is gathered, and the side slope each sampling instant gathered surface image synchronization is sent to image processing equipment (1);
Described monitored area is region side slope surface needing carry out Light deformation monitoring, and institute monitored area is divided into K region to be analyzed, and the individual described region to be analyzed of K is square and its size is all identical; Wherein, K is positive integer and K >=2; Described side slope surface image is made up of K area image, and K described area image is respectively the image in K described region to be analyzed;
The sampling instant of described image capture device, is denoted as t h, t h=t 1+ (h-1) × Δ t; Wherein, h be positive integer and h=1,2,3 ..., t 1for the initial samples moment of described image capture device, Δ t is time interval of adjacent double sampling before and after described image capture device and Δ t=1/f;
Step 2, image-receptive and synchronously process: after described image processing equipment (1) receives the side slope surface image that image capture device (2) each sampling instant gathers, process received side slope surface image synchronization, process is as follows:
The side slope surface image procossing that step 201, initial samples moment gather: adopt image processing equipment (1) to process the side slope surface image that initial samples moment image capture device (2) gathers, comprise the following steps:
Step 2011, side slope surface image store: store the side slope surface image be currently received, and the sampling instant t to this side slope surface image 1carry out record;
Step 2012, elevation information obtain: call image processing module and process the side slope surface image received in step 2011, obtain the altitude figures of each pixel on described side slope surface image, and on the described side slope surface image obtained, the altitude figures of each pixel is t 1the elevation information of monitored area described in the moment;
The elevation information of described monitored area comprises the elevation information in K described region to be analyzed, the altitude figures of each pixel on the area image that the elevation information in each described region to be analyzed includes this region to be analyzed;
T 1described in moment, the elevation information of monitored area comprises t 1the elevation information in moment K described region to be analyzed; Wherein, t 1moment is numbered the elevation information in the region to be analyzed of k, is denoted as X k(t 1); K is the numbering in described region to be analyzed, k is positive integer and k=1,2 ..., K;
Step 2013, elevation information store: to the t obtained in step 2012 1the elevation information in moment K described region to be analyzed, carries out stores synchronized;
Step 2014, Similarity Measure: according to formula (1), t is calculated 1the image similarity X of monitored area described in the moment 1,2 ... K(t 1), wherein X 1,2 ... K(t 1) ∈ (0,1]; In formula (1), x k(t 1) be t 1moment is numbered the elevation mean value in the region to be analyzed of k;
The side slope surface image procossing that step 202, next sampling instant gather: adopt image processing equipment (1) to process the side slope surface image that next sampling instant image capture device (2) gathers, comprise the following steps:
Step 2021, side slope surface image store: store the side slope surface image be currently received, and the sampling time t to this side slope surface image dcarry out record; Wherein, D is positive integer and D>=2;
Step 2022, elevation information obtain: according to the method described in step 2012, obtain t dthe elevation information of monitored area described in the moment;
T ddescribed in moment, the elevation information of monitored area comprises t dthe elevation information in moment K described region to be analyzed; Wherein, t dmoment is numbered the elevation information in the region to be analyzed of k, is denoted as X k(t d);
Step 2023, elevation information store: to the t obtained in step 2022 dthe elevation information in moment K described region to be analyzed, carries out stores synchronized;
Step 2024, Similarity Measure: according to formula (2), t is calculated dthe image similarity X of monitored area described in the moment 1,2 ... K(t d), wherein X 1,2 ... K(t d) ∈ (0,1]; In formula (2), x k(t d) be t dmoment is numbered the elevation mean value in the region to be analyzed of k;
Step 203, Light deformation judge: by the t calculated in step 2024 dthe image similarity X of monitored area described in the moment 1,2 ... K(t d), respectively with t dbefore moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 ... K(t d) and t dbefore moment monitored area described in each moment image similarity between difference when being all not more than ε, illustrate from t 1moment is to t ddescribed in moment there is not Light deformation in monitored area, returns step 202 afterwards, processes the side slope surface image that next sampling instant gathers; Otherwise, t is described dthere is Light deformation in monitored area described in the moment, described image processing equipment (1) controls alarm unit (5) and reports to the police, and to the t of Light deformation occurs dmoment carries out record;
Wherein, ε is the Light deformation alarm threshold value and ε=0.01 ~ 0.1 that preset.
2. according to the side slope Light deformation monitoring method based on similarity decision criteria according to claim 1, it is characterized in that: when carrying out elevation information acquisition in step 2012 and step 2022, call image processing module and all according to three frequency color fringe projection method for three-dimensional measurement, obtain the altitude figures of each pixel on described side slope surface image; The side slope surface image gathered in step one is the distortion color fringe image of described monitored area;
Before carrying out image acquisition and synchronized upload in step one, first adopt image processing equipment (1) to generate color fringe image, RGB tri-Color Channels of this color fringe image are generated by sine streak that is low, that neutralize high three kinds of carrier frequency respectively; Afterwards, then be projected on described monitored area by described color fringe image through RGB tri-color channels of digital projector (3) simultaneously; Then, image capture apparatus is adopted to take the distortion color fringe image of described monitored area;
When carrying out elevation information acquisition in step 2012 and step 2022, process is as follows:
Step I, background subtraction and color decoupling zero: subtract each other containing the Color Channel of high frequency fringes and the Color Channel containing intermediate frequency striped in the distortion color fringe image that image capture device (2) is gathered, obtain the bar graph of high and low frequency compound, decompose with Bidimensional Empirical Mode Decomposition BEMD again, be separated high, middle carrier component; In like manner, by containing in described distortion color fringe image, the Color Channel of low frequency component subtracts each other, in obtaining, the bar graph of low frequency compound, then with Bidimensional Empirical Mode Decomposition BEMD decompose obtain in, low carrier component;
Step I i, phase demodulating: with two-dimentional short time discrete Fourier transform demodulation, the height obtained, neutralize low each carrier component wrapped phase;
Step I ii, phase unwrapping: with change precision go parcel algorithm by low, in and high carrier frequency component complete wrapped phase successively and launch, obtain the expansion phase place of high frequency carrier frequency item, launch the height of monitored area described in phase recovery thus.
3. according to the side slope Light deformation monitoring method based on similarity decision criteria described in claim 1 or 2, it is characterized in that: after image processing equipment described in step 203 (1) control alarm unit (5) is reported to the police, call elevation information comparison module, to from t 1moment is to t dthe elevation mean value of moment K described analyzed area to be analyzed compares respectively, and according to comparative result, draws t dthe region to be analyzed of moment generation Light deformation; Afterwards, described image processing equipment (1) exports t dthe numbering in the region to be analyzed of moment generation Light deformation;
Wherein, from t 1moment is to t dthe comparative approach of the elevation mean value in moment K described region to be analyzed is all identical; To from t 1moment is to t dwhen the elevation mean value in any one region to be analyzed of moment compares, by t dthe elevation mean value in this region to be analyzed of moment respectively with t 1moment is to t d-1in moment, the elevation mean value in this region to be analyzed of each moment carries out difference comparsion.
4. according to the side slope Light deformation monitoring method based on similarity decision criteria described in claim 1 or 2, it is characterized in that: before image processing equipment described in step 203 (1) control alarm unit (5) is reported to the police, described image processing equipment (1) also needs to call pseudo-distortion judge module, to t dwhether the Light deformation that monitored area described in the moment occurs is that pseudo-distortion judges, process is as follows:
Step I, continuously monitoring: according to the method described in step 2021 to step 2024, to t d+1moment is to t d+Ethe side slope surface image that moment image capture device (2) gathers processes, and calculates the image similarity X of monitored area described in each moment 1,2 ... K(t f); Wherein, F is positive integer and F=D+1 ~ D+E; Wherein, t d+E-t d+1=Δ T, Δ T are the continuous monitoring time preset; E be positive integer and
Step II, pseudo-distortion judge: by t d+1moment is to t d+Ethe image similarity X of monitored area described in each moment in moment 1,2 ... K(t f), respectively with t dthe image similarity X of monitored area described in each moment before moment 1,2 ... K(t h) compare: work as X 1,2 ... K(t f) and t dthe image similarity X of monitored area described in each moment before moment 1,2 ... K(t h) between difference when being all greater than ε, t is described dthe Light deformation that monitored area described in moment occurs is true distortion, and now described image processing equipment (1) control alarm unit (5) is reported to the police; Otherwise, t is described dthe Light deformation that monitored area described in moment occurs is pseudo-distortion.
5., according to the side slope Light deformation monitoring method based on similarity decision criteria according to claim 4, it is characterized in that: Δ T=2s ~ 3min in step I.
6. according to the side slope Light deformation monitoring method based on similarity decision criteria described in claim 1 or 2, it is characterized in that: image processing equipment described in step 203 (1) control alarm unit (5) enters step 3 after reporting to the police;
Continue to monitor after step 3, first generation Light deformation, comprise the following steps:
Step 301, next sampling instant image procossing: according to the method described in step 2021 to step 2024, to t mthe side slope surface image that moment image capture device (2) gathers processes, and calculates t mthe image similarity X of monitored area described in the moment 1,2 ... K(t m); Wherein, M is positive integer and M>=D+1;
Step 302, Light deformation judge: by the t calculated in step 301 mthe image similarity X of monitored area described in the moment 1,2 ... K(t m), respectively with t dmoment is to t min moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 ... K(t m) and t dmoment is to t min moment monitored area described in each moment image similarity between difference when being all not more than ε, illustrate from t dmoment is to t mdescribed in moment there is not Light deformation in monitored area, returns step 301 afterwards, processes the side slope surface image that next sampling instant gathers; Otherwise, t is described mthere is Light deformation in monitored area described in the moment, described image processing equipment (1) controls alarm unit (5) and reports to the police, and to the t of Light deformation occurs mmoment and from t 1moment is to t mthe alarm times of alarm unit described in the moment (5) carries out record, enters step 4 afterwards;
Continue to monitor after step 4, next time generation Light deformation, comprise the following steps:
Step 401, next sampling instant image procossing: according to the method described in step 2021 to step 2024, to t nthe side slope surface image that moment image capture device (2) gathers processes, and calculates t nthe image similarity X of monitored area described in the moment 1,2 ... K(t n); Wherein, N is positive integer and N>=n, and wherein n is positive integer and t nmoment is t nthe generation moment of the last Light deformation in described monitored area before moment;
Step 402, Light deformation judge: by the t calculated in step 401 nthe image similarity X of monitored area described in the moment 1,2 ... K(t n), respectively with t nmoment is to t nin moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 ... K(t m) and t nmoment is to t nin moment monitored area described in each moment image similarity between difference when being all not more than ε, illustrate from t nmoment is to t ndescribed in moment there is not Light deformation in monitored area, returns step 401 afterwards, processes the side slope surface image that next sampling instant gathers; Otherwise, t is described nthere is Light deformation in monitored area described in the moment, described image processing equipment (1) controls alarm unit (5) and reports to the police, and to the t of Light deformation occurs nmoment and from t 1moment is to t nthe alarm times of alarm unit described in the moment (5) carries out record, enters step 5 afterwards;
Step 5, one or many repeat step 4, continue to monitor described monitored area.
7. according to the side slope Light deformation monitoring method based on similarity decision criteria according to claim 6, it is characterized in that: before image processing equipment described in step 302 (1) control alarm unit (5) is reported to the police, described image processing equipment (1) also needs to call pseudo-distortion judge module, to t mwhether the Light deformation that monitored area described in the moment occurs is that pseudo-distortion judges, process is as follows:
Step I-11, continuously monitoring: according to the method described in step 2021 to step 2024, to t m+1moment is to t m+Ethe side slope surface image that moment image capture device (2) gathers processes, and calculates the image similarity X of monitored area described in each moment 1,2 ... K(t f1); Wherein, F1 is positive integer and F1=M+1 ~ M+E; Wherein, t m+E-t m+1=Δ T;
Step II-12, pseudo-distortion judge: by t m+1moment is to t m+Ethe image similarity X of monitored area described in each moment in moment 1,2 ... K(t f1), respectively with t dmoment is to t min moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 ... K(t f1) and t dmoment is to t min moment monitored area described in each moment image similarity between difference when being all greater than ε, t is described mthe Light deformation that monitored area described in moment occurs is true distortion, and now described image processing equipment (1) control alarm unit (5) is reported to the police; Otherwise, t is described mthe Light deformation that monitored area described in moment occurs is pseudo-distortion;
Before image processing equipment described in step 402 (1) control alarm unit (5) is reported to the police, described image processing equipment (1) also needs to call pseudo-distortion judge module, to t nwhether the Light deformation that monitored area described in the moment occurs is that pseudo-distortion judges, process is as follows:
Step I-21, continuously monitoring: according to the method described in step 2021 to step 2024, to t n+1moment is to t n+Ethe side slope surface image that moment image capture device (2) gathers processes, and calculates the image similarity X of monitored area described in each moment 1,2 ... K(t f2); Wherein, F2 is positive integer and F2=N+1 ~ N+E; Wherein, t m+E-t m+1=Δ T;
Step II-22, pseudo-distortion judge: by t n+1moment is to t n+Ethe image similarity X of monitored area described in each moment in moment 1,2 ... K(t f2), respectively with t nmoment is to t nin moment, described in each moment, the image similarity of monitored area compares: work as X 1,2 ... K(t f2) and t nmoment is to t nin moment monitored area described in each moment image similarity between difference when being all greater than ε, t is described nthe Light deformation that monitored area described in moment occurs is true distortion, and now described image processing equipment (1) control alarm unit (5) is reported to the police; Otherwise, t is described nthe Light deformation that monitored area described in moment occurs is pseudo-distortion.
8. according to the side slope Light deformation monitoring method based on similarity decision criteria according to claim 6, it is characterized in that: after image processing equipment described in step 302 (1) control alarm unit (5) is reported to the police, call elevation information comparison module, to from t dmoment is to t mthe elevation mean value of moment K described analyzed area to be analyzed compares respectively, and according to comparative result, draws t mthe region to be analyzed of moment generation Light deformation; Afterwards, described image processing equipment (1) exports t mthe numbering in the region to be analyzed of moment generation Light deformation;
Wherein, from t dmoment is to t mthe comparative approach of the elevation mean value in moment K described region to be analyzed is all identical; To from t dmoment is to t mwhen the elevation mean value in any one region to be analyzed of moment compares, by t mthe elevation mean value in this region to be analyzed of moment respectively with t dmoment is to t m-1in moment, the elevation mean value in this region to be analyzed of each moment carries out difference comparsion;
After image processing equipment described in step 402 (1) controls alarm unit (5) report to the police, call elevation information comparison module, to from t nmoment is to t nthe elevation mean value of moment K described analyzed area to be analyzed compares respectively, and according to comparative result, draws t nthe region to be analyzed of moment generation Light deformation; Afterwards, described image processing equipment (1) exports t nthe numbering in the region to be analyzed of moment generation Light deformation;
Wherein, from t nmoment is to t nthe comparative approach of the elevation mean value in moment K described region to be analyzed is all identical; To from t nmoment is to t nwhen the elevation mean value in any one region to be analyzed of moment compares, by t nthe elevation mean value in this region to be analyzed of moment respectively with t nmoment is to t n-1in moment, the elevation mean value in this region to be analyzed of each moment carries out difference comparsion.
9. one kind utilizes side slope Light deformation monitoring method as claimed in claim 2 to carry out the system of side slope Light deformation, it is characterized in that: comprise the digital projector (3) that K is laid in K described overlying regions to be analyzed respectively, K respectively to the image capture apparatus that the image in K described region to be analyzed absorbs, the image capture device (2) that the side slope surface image of described monitored area is gathered and the image processing equipment (1) connected with image capture device (2), the individual described image capture apparatus of K is CCD camera (4) and it all connects with image capture device (2), K described digital projector (3) all connects with image processing equipment (1), K described CCD camera (4) is laid in K described overlying regions to be analyzed respectively.
10. according to system according to claim 9, it is characterized in that: described digital projector (3) and CCD camera (4) lay respectively at above the both sides in described region to be analyzed.
CN201510242170.5A 2015-05-13 2015-05-13 Side slope Light deformation monitoring method and monitoring system based on similarity decision criteria Active CN105043355B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510242170.5A CN105043355B (en) 2015-05-13 2015-05-13 Side slope Light deformation monitoring method and monitoring system based on similarity decision criteria

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510242170.5A CN105043355B (en) 2015-05-13 2015-05-13 Side slope Light deformation monitoring method and monitoring system based on similarity decision criteria

Publications (2)

Publication Number Publication Date
CN105043355A true CN105043355A (en) 2015-11-11
CN105043355B CN105043355B (en) 2017-05-31

Family

ID=54450093

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510242170.5A Active CN105043355B (en) 2015-05-13 2015-05-13 Side slope Light deformation monitoring method and monitoring system based on similarity decision criteria

Country Status (1)

Country Link
CN (1) CN105043355B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106202908A (en) * 2016-07-05 2016-12-07 西安交通大学 A kind of high slope relaxes the decision method in district
CN107490368A (en) * 2017-08-02 2017-12-19 上海华测导航技术股份有限公司 A kind of side slope setting-out system and method
CN111561917A (en) * 2020-03-30 2020-08-21 同济大学 Road side slope monitoring system
CN114285773A (en) * 2021-11-19 2022-04-05 中国联合网络通信集团有限公司 Flow sampling method, device, equipment and readable storage medium
CN115620489A (en) * 2022-10-24 2023-01-17 深圳市自然资源和不动产评估发展研究中心(深圳市地质环境监测中心) Slope monitoring and early warning method, system, equipment and medium based on live broadcast plug flow

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102052915A (en) * 2010-11-24 2011-05-11 航天科工深圳（集团）有限公司 Geological disaster monitoring system and method
CN102538751A (en) * 2012-01-19 2012-07-04 四川金码科技有限公司 Dynamic and static self-adjusting type settlement/deflection monitoring system and method
CN102788558A (en) * 2012-07-30 2012-11-21 山东师范大学 Three-dimensional deformation measuring system and three-dimensional deformation measuring method combining speckle correlation and speckle interference
CN103018757A (en) * 2012-12-10 2013-04-03 中山大学 Monitoring and display system for position and three-dimensional attitude of floating and semi-submersible type mobile platform
CN202928571U (en) * 2012-11-05 2013-05-08 上海米度测控科技有限公司 Photographic measurement system for long-distance micro-deformation monitoring
US20140253689A1 (en) * 2013-03-08 2014-09-11 Kabushiki Kaisha Topcon Measuring Instrument

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102052915A (en) * 2010-11-24 2011-05-11 航天科工深圳（集团）有限公司 Geological disaster monitoring system and method
CN102538751A (en) * 2012-01-19 2012-07-04 四川金码科技有限公司 Dynamic and static self-adjusting type settlement/deflection monitoring system and method
CN102788558A (en) * 2012-07-30 2012-11-21 山东师范大学 Three-dimensional deformation measuring system and three-dimensional deformation measuring method combining speckle correlation and speckle interference
CN202928571U (en) * 2012-11-05 2013-05-08 上海米度测控科技有限公司 Photographic measurement system for long-distance micro-deformation monitoring
CN103018757A (en) * 2012-12-10 2013-04-03 中山大学 Monitoring and display system for position and three-dimensional attitude of floating and semi-submersible type mobile platform
US20140253689A1 (en) * 2013-03-08 2014-09-11 Kabushiki Kaisha Topcon Measuring Instrument

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106202908A (en) * 2016-07-05 2016-12-07 西安交通大学 A kind of high slope relaxes the decision method in district
CN106202908B (en) * 2016-07-05 2018-10-19 西安交通大学 A kind of determination method in high slope relaxation area
CN107490368A (en) * 2017-08-02 2017-12-19 上海华测导航技术股份有限公司 A kind of side slope setting-out system and method
CN107490368B (en) * 2017-08-02 2020-06-12 上海华测导航技术股份有限公司 Slope lofting system and method
CN111561917A (en) * 2020-03-30 2020-08-21 同济大学 Road side slope monitoring system
CN111561917B (en) * 2020-03-30 2021-10-26 同济大学 Road side slope monitoring system
CN114285773A (en) * 2021-11-19 2022-04-05 中国联合网络通信集团有限公司 Flow sampling method, device, equipment and readable storage medium
CN114285773B (en) * 2021-11-19 2023-11-10 中国联合网络通信集团有限公司 Traffic sampling method, device, equipment and readable storage medium
CN115620489A (en) * 2022-10-24 2023-01-17 深圳市自然资源和不动产评估发展研究中心(深圳市地质环境监测中心) Slope monitoring and early warning method, system, equipment and medium based on live broadcast plug flow

Also Published As

Publication number Publication date
CN105043355B (en) 2017-05-31

Similar Documents

Publication Publication Date Title
US10480939B2 (en) High speed stereoscopic pavement surface scanning system and method
US9970758B2 (en) High speed stereoscopic pavement surface scanning system and method
US10190269B2 (en) High speed stereoscopic pavement surface scanning system and method
CN103471512B (en) A kind of glass plate width detecting system based on machine vision
CN105043355A (en) Side slope micro-deformation monitoring method and side slope micro-deformation monitoring system based on similarity determination criterion
Pages et al. Optimised De Bruijn patterns for one-shot shape acquisition
CN102135236B (en) Automatic non-destructive testing method for internal wall of binocular vision pipeline
KR20130099735A (en) Method and fusion system of time-of-flight camera and stereo camera for reliable wide range depth acquisition
CN103900494B (en) For the homologous points fast matching method of binocular vision 3 D measurement
CN102538705B (en) Secondary-projection-algorithm-based on-line non-contact contour detection system and method of intermediate-thick plate
CN104567726B (en) Vehicle operation troubles detecting system and method
US20130076871A1 (en) High speed photometric stereo pavement scanner
KR20150121179A (en) Real time stereo matching
Staniek Stereo vision method application to road inspection
JP6858415B2 (en) Sea level measurement system, sea level measurement method and sea level measurement program
CN103837087B (en) Pantograph automatic testing method based on active shape model
CN103630088A (en) High-precision tunnel cross section detection method and device based on double laser bands
KR20160121509A (en) Structured light matching of a set of curves from two cameras
CN109341668A (en) Polyphaser measurement method based on refraction projection model and beam ray tracing method
CN103438834A (en) Hierarchy-type rapid three-dimensional measuring device and method based on structured light projection
Beltran et al. A comparison between active and passive 3d vision sensors: Bumblebeexb3 and Microsoft Kinect
KR102004909B1 (en) Tongue rail inspection system for detecting the close contact
TWI632347B (en) Method for integrating three-dimensional image and laser scanning ranging
CN106197320B (en) A kind of scanning of time-sharing multiplex quick three-dimensional and its data processing method
CN113421286B (en) Motion capturing system and method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20221117

Address after: Room 416-1, 4th Floor, Building 3, Fengyue Yunchuang Center, Hanchi 1st Road, Haojing Avenue, Fengdong New Town, Xi'an City, Shaanxi Province 710000

Patentee after: Xi'an Dishan Shiju Technology Co.,Ltd.

Address before: 710054 No. 58, middle section, Yanta Road, Shaanxi, Xi'an

Patentee before: XI'AN University OF SCIENCE AND TECHNOLOGY

TR01 Transfer of patent right