CN106910190A - For the squatter building scene Recognition system of land resources monitoring - Google Patents

For the squatter building scene Recognition system of land resources monitoring Download PDF

Info

Publication number
CN106910190A
CN106910190A CN201710102202.0A CN201710102202A CN106910190A CN 106910190 A CN106910190 A CN 106910190A CN 201710102202 A CN201710102202 A CN 201710102202A CN 106910190 A CN106910190 A CN 106910190A
Authority
CN
China
Prior art keywords
region
contrast images
benchmark image
piece
treat
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.)
Pending
Application number
CN201710102202.0A
Other languages
Chinese (zh)
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.)
XIAN TIANHE DEFENCE TECHNOLOGY Co Ltd
Original Assignee
XIAN TIANHE DEFENCE TECHNOLOGY Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by XIAN TIANHE DEFENCE TECHNOLOGY Co Ltd filed Critical XIAN TIANHE DEFENCE TECHNOLOGY Co Ltd
Priority to CN201710102202.0A priority Critical patent/CN106910190A/en
Publication of CN106910190A publication Critical patent/CN106910190A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • G06T7/001Industrial image inspection using an image reference approach
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20056Discrete and fast Fourier transform, [DFT, FFT]

Landscapes

  • Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)

Abstract

The disclosure is directed to a kind of squatter building scene Recognition system for land resources monitoring.The system includes:The first two field picture of predeterminated position is obtained as benchmark image and using the two field picture after the first two field picture as treating contrast images;To benchmark image and treat that contrast images carry out piecemeal treatment the benchmark image and to treat that contrast images are divided into corresponding multiple block regions respectively;To benchmark image and in treating contrast images, corresponding a pair of block region carries out Fourier transformation and is calculated crosspower spectrum, then inverse Fourier transform is carried out to crosspower spectrum obtains a spatial transform matrix, and calculates the maximum of spatial transform matrix;Take benchmark image and treat corresponding piece of region of maximum of the minimum spatial transform matrix of top n in the maximum of the corresponding spatial transform matrix in corresponding each pair block region in contrast images as the region of variation treated in contrast images.The disclosure can improve the precision and accuracy of squatter building scene Recognition, it is to avoid system is reported by mistake and failed to report.

Description

For the squatter building scene Recognition system of land resources monitoring
Technical field
This disclosure relates to land and resource information monitoring technology field, more particularly to a kind of squatter building for land resources monitoring Scene Recognition system.
Background technology
With the fast development of China's economy, soil imbalance between supply and demand becomes increasingly conspicuous, and illegal building occupies cultivated land phenomenon, city Illegal land used or illegal land used build stealing for phenomenon and mineral resources and adopt illegal mining phenomenon and occur repeatedly.At present in land resources Monitoring aspect mainly monitors land use change survey situation using Satellite Remote Sensing technological means, during by before and after different years Between soil Remote Sensing Imagery Change verify each place illegal land situation.However, satellite monitoring is more to be carried out from State-level Macroscopic supervision, is related to zonal monitoring resource, then more to be supervised using video monitoring system.
At present, the video monitoring system of land resources there are problems that in practical engineering application, such as in complicated field Under scape (such as field arable land, forest), due to by dynamic background (such as leaf is rocked, the rustle of leaves in the wind, illuminance abrupt variation) for disobeying The influence of scene Recognition is built, recognition accuracy is low, add that scene is complicated and changeable, cause the video monitoring system of land resources easily to be sent out It is raw to report by mistake and fail to report, larger burden is caused to subsequent treatment, cause some unnecessary erroneous judgements.
Therefore, it is necessary to providing a kind of new technical scheme improves one or more problem present in such scheme.
It should be noted that information is only used for strengthening the reason of background of this disclosure disclosed in above-mentioned background section Solution, therefore can include not constituting the information to prior art known to persons of ordinary skill in the art.
The content of the invention
The purpose of the disclosure is to provide a kind of squatter building scene Recognition system for land resources monitoring, and then is at least existed Overcome to a certain extent due to one or more problem caused by the limitation of correlation technique and defect.
Other characteristics and advantage of the disclosure will be apparent from by following detailed description, or partially by the disclosure Practice and acquistion.
According to the first aspect of the embodiment of the present disclosure, there is provided a kind of squatter building scene Recognition system for land resources monitoring System, the system includes:
Image collection module, for obtaining the first two field picture of land resources region to be monitored predeterminated position as reference map Picture, and the two field picture after first two field picture is obtained as treating contrast images;
Image block module, for the benchmark image and treating that contrast images carry out piecemeal treatment with by the base respectively Quasi- image and treat that contrast images are divided into corresponding multiple block regions;
Image processing module, for corresponding a pair of block region to be carried out in Fu to the benchmark image and in treating contrast images Leaf transformation is simultaneously calculated crosspower spectrum, then one spatial transform matrix of inverse Fourier transform acquisition is carried out to the crosspower spectrum simultaneously Calculate the maximum of the spatial transform matrix;And
Area determination module, for obtaining the benchmark image and corresponding each pair block region is corresponding in treating contrast images Corresponding piece of region conduct of maximum of the minimum spatial transform matrix of top n in the maximum of the spatial transform matrix The region of variation treated in contrast images;Wherein N is more than or equal to 2 and less than the described piece of integer of the total number in region.
In a kind of exemplary embodiment of the disclosure, described image processing module is used for:
The benchmark image is calculated according to below equation and the i-th mutual work(in corresponding a pair of block region in contrast images is treated Rate is composed:
Wherein, F1i(u, v) is f1iThe Fourier transformation of (x, y), f1i(x, y) represents i-th piece of area in the benchmark image Domain;
F2i(u, v) is f2iThe Fourier transformation of (x, y), f2i(x, y) treats i-th piece of region in contrast images described in representing,It is F2iThe conjugation of (u, v);
x0,y0Respectively it is described treat contrast images in i-th piece of region relative to i-th piece of region in the benchmark image In the side-play amount in x and y directions.
In a kind of exemplary embodiment of the disclosure, the value of the N is 5.
In a kind of exemplary embodiment of the disclosure, described image processing module is used for:
To the benchmark image and treat that contrast images carry out uniform non-overlapping piecemeal treatment and form multiple block regions respectively.
In a kind of exemplary embodiment of the disclosure, the benchmark image is 1920* with the resolution ratio for treating contrast images 1080, the size in described piece of region is 120*120.
In a kind of exemplary embodiment of the disclosure, each described piece of region has station location marker information;The region is true Cover half block, is used for:
Obtain the benchmark image and treat the corresponding spatial transform matrix in corresponding each pair block region in contrast images Maximum in corresponding piece of region of maximum of the minimum spatial transform matrix of top n treat contrast images as described Region of variation include:
Obtain the benchmark image and treat the corresponding spatial transform matrix in corresponding each pair block region in contrast images Maximum in the minimum spatial transform matrix of top n corresponding piece of region of maximum and record the N number of piece of region Station location marker information;
The N number of piece of area maps are treated contrast images by station location marker information according to the N number of piece of region to described As the region of variation.
In a kind of exemplary embodiment of the disclosure, the system also includes false-alarm judge module, is used for:
The gray level co-occurrence matrixes in each the block region in the region of variation are calculated, and calculates the gray level co-occurrence matrixes Characteristic information;
Judge whether the characteristic information meets a default criterion, the corresponding institute of this feature information is removed if meeting State a block region in region of variation;
The region of variation behind a block region as object variations region will be removed.
In a kind of exemplary embodiment of the disclosure, the characteristic information of the gray level co-occurrence matrixes includes contrast, energy At least one of with entropy.
In a kind of exemplary embodiment of the disclosure, the false-alarm judge module is used for:
Whether whether the contrast is judged more than first threshold, and judges the energy value less than Second Threshold, while Judge the entropy whether more than the 3rd threshold value;
When the contrast is more than the first threshold, the energy value is big less than the Second Threshold and the entropy When the 3rd threshold value meets simultaneously, the characteristic information meets the default criterion.
In a kind of exemplary embodiment of the disclosure, the first threshold is 0.6, and the Second Threshold is 0.12, described 3rd threshold value is 2.1.
The technical scheme provided by this disclosed embodiment can include the following benefits:
In a kind of embodiment of the disclosure, by the above-mentioned squatter building scene Recognition system monitored for land resources, by base Quasi- image and treat that contrast images piecemeal forms corresponding piece of region, frequency is transformed into from time domain by regarding block region as 2D signal Domain simultaneously obtains two crosspower spectrums in image corresponding blocks region so as to further infer that out whether scene is sent out using Fourier transformation Changing.Thus it is possible, on the one hand, improve the precision and accuracy rate to the complex scene identification for changing;On the other hand, can be with Realization accurately recognizes early warning to squatter building scene, further avoid land resources monitoring system and wrong report occurs and fails to report, So as to avoid causing larger burden and caused some unnecessary erroneous judgements to subsequent treatment.
It should be appreciated that the general description of the above and detailed description hereinafter are only exemplary and explanatory, not The disclosure can be limited.
Brief description of the drawings
Accompanying drawing herein is merged in specification and constitutes the part of this specification, shows the implementation for meeting the disclosure Example, and it is used to explain the principle of the disclosure together with specification.It should be evident that drawings in the following description are only the disclosure Some embodiments, for those of ordinary skill in the art, on the premise of not paying creative work, can also basis These accompanying drawings obtain other accompanying drawings.
Fig. 1 schematically shows the squatter building scene Recognition system monitored for land resources in disclosure exemplary embodiment and shows It is intended to;
Fig. 2 schematically shows another squatter building scene Recognition system for land resources monitoring in disclosure exemplary embodiment System schematic diagram;
Fig. 3 schematically shows the another squatter building scene Recognition system for land resources monitoring in disclosure exemplary embodiment System schematic diagram;
Fig. 4 schematically shows another squatter building scene Recognition system for being used for land resources monitoring in disclosure exemplary embodiment System schematic diagram;
Fig. 5 A schematically show the benchmark image schematic diagram in disclosure exemplary embodiment Scene one;
Fig. 5 B schematically show in disclosure exemplary embodiment Scene one and treat contrast images schematic diagram;
Fig. 5 C schematically show the change identification test result schematic diagram of disclosure exemplary embodiment Scene one;
Fig. 6 A schematically show the benchmark image schematic diagram in disclosure exemplary embodiment Scene two;
Fig. 6 B schematically show in disclosure exemplary embodiment Scene two and treat contrast images schematic diagram;
Fig. 6 C schematically show the change identification test result schematic diagram of disclosure exemplary embodiment Scene two.
Fig. 7 A schematically show the benchmark image schematic diagram in disclosure exemplary embodiment Scene three;
Fig. 7 B schematically show in disclosure exemplary embodiment Scene three and treat contrast images schematic diagram;
Fig. 7 C schematically show the change identification test result schematic diagram of disclosure exemplary embodiment Scene three.
Specific embodiment
Example embodiment is described more fully with referring now to accompanying drawing.However, example embodiment can be with various shapes Formula is implemented, and is not understood as limited to example set forth herein;Conversely, thesing embodiments are provided so that the disclosure will more Fully and completely, and by the design of example embodiment those skilled in the art is comprehensively conveyed to.Described feature, knot Structure or characteristic can be combined in one or more implementation methods in any suitable manner.
Additionally, accompanying drawing is only the schematic illustrations of the disclosure, it is not necessarily drawn to scale.Identical accompanying drawing mark in figure Note represents same or similar part, thus will omit repetition thereof.Some block diagrams shown in accompanying drawing are work( Energy entity, not necessarily must be corresponding with physically or logically independent entity.These work(can be realized using software form Energy entity, or these functional entitys are realized in one or more hardware modules or integrated circuit, or at heterogeneous networks and/or place These functional entitys are realized in reason device device and/or microcontroller device.
A kind of squatter building scene Recognition system for land resources monitoring is provide firstly in this example embodiment.With reference to Shown in Fig. 1, the system 100 can include image collection module 101, image block module 102, the and of image processing module 103 Area determination module 104.Wherein:
Described image acquisition module 101, the first two field picture for obtaining land resources region to be monitored predeterminated position is made On the basis of image, and obtain the two field picture after first two field picture as treating contrast images;
Described image piecemeal module 102, for the benchmark image and treat contrast images carry out respectively piecemeal treatment with The benchmark image and will treat that contrast images are divided into corresponding multiple block regions;
Described image processing module 103, for the benchmark image and treating corresponding a pair of block region in contrast images Carry out Fourier transformation and be calculated crosspower spectrum, then inverse Fourier transform is carried out to the crosspower spectrum to obtain spatial domain change Change matrix and calculate the maximum of the spatial transform matrix;
The area determination module 104, for obtaining the benchmark image and treating corresponding each pair block area in contrast images Corresponding piece of the maximum of the minimum spatial transform matrix of top n in the maximum of the corresponding spatial transform matrix in domain Region is used as the region of variation treated in contrast images;Wherein N is more than or equal to 2 and less than the described piece of total number in region Integer.
By the above-mentioned squatter building scene Recognition system monitored for land resources, on the one hand, improve to changing The precision and accuracy rate of complex scene identification;On the other hand, it is possible to achieve early warning is accurately recognized to squatter building scene, enters one Step avoids land resources monitoring system be there is wrong report and fails to report, so as to avoid causing larger burden to subsequent treatment and lead Some the unnecessary erroneous judgements for causing.
Below, the unit of the said system in this example embodiment will be carried out in more detail referring to figs. 1 to Fig. 4 Explanation.
Described image acquisition module 101 is used to obtain the first two field picture conduct of land resources region to be monitored predeterminated position Benchmark image, and the two field picture after first two field picture is obtained as treating contrast images.
Exemplary, the predeterminated position is same presetting bit under complex scene, and the presetting bit can be as the case may be It is selected.To certain specific presetting bit, the first two field picture that will be gathered in the presetting bit is used as benchmark image, the one of subsequent acquisition Two field picture is used as pending image.Can be adopted by the image processing equipment such as camera or photoelectric monitor equipment during specific collection Collection image.
Described image piecemeal module 102 is used for the benchmark image and treats that contrast images carry out piecemeal treatment to incite somebody to action respectively The benchmark image and treat that contrast images are divided into corresponding multiple block regions.
In this example embodiment, it is described to the benchmark image and treat contrast images carry out respectively piecemeal treatment include: To the benchmark image and treat that contrast images carry out uniform non-overlapping piecemeal treatment and form multiple block regions respectively.
It is exemplary, the benchmark image and treat that the resolution ratio of contrast images is 1920*1080, described piece of region it is big It is small to be 120*120.Because pending image resolution ratio is 1920*1080, size is larger, and treatment effect is carried out to entire image It is not good, small change can be flooded, while region of variation can not be positioned.Therefore, to pending two images Piecemeal is carried out, the size of block is taken here for 120*120, so, for the image of 1920*1080 sizes, the number of the block being divided into It is 16*9=144.
Described image processing module 103 is used for the benchmark image and corresponding a pair of block region is entered in treating contrast images Row Fourier transformation is simultaneously calculated crosspower spectrum, then inverse Fourier transform one spatial transform of acquisition is carried out to the crosspower spectrum Matrix simultaneously calculates the maximum of the spatial transform matrix.
In this example embodiment, described to the benchmark image and in treating contrast images, corresponding a pair of block region is carried out Fourier transformation is simultaneously calculated crosspower spectrum and can include:The benchmark image is calculated according to below equation and contrast images are treated In i-th crosspower spectrum in corresponding a pair of block region:
Wherein, F1i(u, v) is f1iThe Fourier transformation of (x, y), f1i(x, y) represents i-th piece of area in the benchmark image Domain;F2i(u, v) is f2iThe Fourier transformation of (x, y), f2i(x, y) treats i-th piece of region in contrast images described in representing,It is F2iThe conjugation of (u, v);x0,y0Respectively it is described treat contrast images in i-th piece of region relative to the reference map Side-play amount of the i-th piece of region in x and y directions as in.
It is exemplary, to i-th piece of region in benchmark image and it is corresponding treat contrast images in i-th piece of region enter respectively Row Fourier transformation, then calculates crosspower spectrum, inverse Fourier transform is carried out to crosspower spectrum and obtains a spatial transform matrix simultaneously The maximum of the spatial transform matrix is calculated, the calculating of matrix maximum refers to existing mature technology, repeats no more.
Ideally, the above-mentioned inverse Fourier transform that carried out to crosspower spectrum can obtain an impulse function, the impact letter Several is not 0 on corresponding translation position, is 0 in other positions.And in practice, due to being influenceed by various noises, gained The spatial transform matrix occurrence for arriving is not 0 phenomenon substantially, but still has a peak-peak, the corresponding position of the peak value Translational movement as between two images.Under normal circumstances, the translational movement that can be calculated according to intersection crosspower spectrum will scheme As carrying out registration, subsequently to carry out the application of image mosaic etc.However, whether being indifferent to image in this example embodiment There is translation, i.e., do not carry out registration, it is that it is paid close attention to it is important that whether correspondence two images occur whether the change of target, i.e. scene become Change, for example whether the building for having squatter building occurs.In this example embodiment, according to signal processing theory, if two signals do not have Change, the amplitude of its crosspower spectrum is essentially 1;, whereas if two signals to be compared change, then its amplitude Greatly reduce.Here, regard image as 2D signal, calculate the intersection crosspower spectrum of two images corresponding blocks, the reflection of its amplitude Situation of change between correspondence image.It is theoretical according to this, it is possible to use the amplitude of the intersection crosspower spectrum of two images becomes Change to judge to treat whether contrast images change, so that the accuracy rate and precision that identification judges are greatly improved.
The area determination module 104 is used to obtain the benchmark image and treats corresponding each pair block region in contrast images Corresponding piece of area of maximum of the minimum spatial transform matrix of top n in the maximum of the corresponding spatial transform matrix Domain is used as the region of variation treated in contrast images;Wherein N is more than or equal to 2 and whole less than described piece of total number in region Number.
With reference to shown in Fig. 2, each described piece of region has station location marker information, such as unique call number or coding Deng.The area determination module 104 obtains the benchmark image and treats the corresponding institute in corresponding each pair block region in contrast images Corresponding piece of region of maximum of the spatial transform matrix of top n minimum in the maximum of spatial transform matrix is stated as institute Stating the region of variation treated in contrast images can include step 201~step 202.Wherein:
Step 201:Obtain the benchmark image and treat the corresponding spatial domain in corresponding each pair block region in contrast images Corresponding piece of region of maximum of the minimum spatial transform matrix of top n and to record this N number of in the maximum of transformation matrix The station location marker information in block region.
It is exemplary, it is above-mentioned to carry out the spatial transform matrix equal with block area size of generation after inverse Fourier transform.So Afterwards, the maximum of spatial transform matrix is counted, maximum and its corresponding piece of index in region, i.e. position identification information is preserved.
Step 202:The N number of piece of area maps are treated right by the station location marker information according to the N number of piece of region described in Than in image as the region of variation.
Exemplary, the value of the N is 5, and other numerical value are may also be certainly.To the corresponding spatial transform in all pieces of regions The maximum of matrix carries out ascending sort, minimum corresponding piece of region of 5 maximums before selection, according to station location marker information In mapping that to treat contrast images, obtain corresponding blocks region in the position in treating contrast images is as contrast images are treated The region that region of variation, i.e. scene change.The region of variation changes maximum region, Ye Jihou in having reacted two images Select region of variation block.
With reference to shown in Fig. 2, rocked to process dynamic background such as leaf, the influence of the change such as the rustle of leaves in the wind, illuminance abrupt variation, Also removed with the characteristic value difference of the gray level co-occurrence matrixes in false-alarm block region using object block region in the present exemplary embodiment Most of false-alarm is improving the accuracy of identification.On the basis of above-described embodiment, in the present exemplary embodiment, the system 100 can with false-alarm judge module 105, for performing following steps S301~S303, with reference to shown in Fig. 3, wherein:
Step S301:The gray level co-occurrence matrixes in each the block region in the region of variation are calculated, and calculates the gray scale The characteristic information of co-occurrence matrix.
Exemplary, the characteristic information of the gray level co-occurrence matrixes includes at least one of contrast, energy and entropy. The gray level co-occurrence matrixes can reflect that image, on direction, adjacent spaces, the integrated information of amplitude of variation, is the office for analyzing image The basis of portion's pattern and their queueing disciplines, is a kind of method being quantitatively described to image texture.Gray level co-occurrence matrixes side Method reflects texture with conditional probability, is the performance of the Gray Correlation of adjacent pixel, and it is defined as in image at a distance of two for D The joint probability distribution that individual gray-scale pixels occur simultaneously, is a kind of second-order statistic, and pixel relationship can be according to different textures Attribute is selected.
Use PδGray level co-occurrence matrixes are represented, it is a matrix of L*L (L is the gray level of image), Pδ(i, j), i, j= 1,2 ..., L represents the times or frequency that two pixels with spatial relation δ and gray scale respectively i and j occur.Here δ=(Dx,Dy), pixel relationship can be selected according to different texture features, i.e. Dx,DyCan be with unrestricted choice.For thinner Texture can take pel spacing for 1, δ=(± 1,0) is horizontal sweep;δ=(0, ± 1) it is vertical scanning;δ=(1, -1), δ =(- 1,1) are 45 degree of scannings;Scanned for 180 degree δ=(1,1), δ=(1,1).Once locus determines, you can generation gray scale Co-occurrence matrix.
If have a larger value according to the unit that the gray level co-occurrence matrixes can be seen that near diagonal, explanatory diagram as Pixel has similar pixel value, is had than larger value if the deviation from cornerwise element, and pixels illustrated gray scale has locally Large change.In order to obtain more textural characteristics, we calculate contrast, energy and entropy here, and specific calculation is as follows:
First, contrast:
Texture rill is deeper, and its contrast is bigger, and visual effect is more clear;Conversely, contrast is small, then rill is shallow, effect It is fuzzy.Gray scale difference is the big pixel of contrast to more, and the value of contrast is bigger.Away from cornerwise unit in gray level co-occurrence matrixes Element value is bigger, and contrast con is bigger.So the bigger images of con are more clear.Contrast computing formula is as follows:
Wherein, P (i, j) represents the i-th row, the value of jth row in gray level co-occurrence matrixes.
2nd, energy
Energy Asm is the quadratic sum of gray level co-occurrence matrixes element value, so also referred to as energy, reflects gradation of image point Cloth uniformity coefficient and texture fineness degree.Asm values show greatly the texture pattern of a kind of more uniform and regular change.It is calculated as follows:
Wherein, P (i, j) represents the i-th row, the value of jth row in gray level co-occurrence matrixes.
3rd, entropy
The physical meaning of entropy Ent is exactly the rule degree of object, and more orderly entropy is smaller, and more unordered entropy is bigger.Entropy is same herein Represent image information content, when in gray level co-occurrence matrixes all elements have maximum randomness, space co-occurrence matrix in all values When almost equal, in gray level co-occurrence matrixes during element dispersed distribution, entropy is larger.It illustrates the non-uniform degree of texture in image Or complexity.It is calculated as follows:
Wherein, P (i, j) represents the i-th row, the value of jth row in gray level co-occurrence matrixes.
Step S302:Judge whether the characteristic information meets a default criterion, this feature letter is removed if meeting Cease a block region in the corresponding region of variation.
With reference to shown in Fig. 4, in the present exemplary embodiment, the false-alarm judge module 105 judges in above-mentioned steps S302 Whether the characteristic information meets a default criterion can include:
Step S401:Whether whether the contrast is judged more than first threshold, and judges the energy value less than second Threshold value, while judging the entropy whether more than the 3rd threshold value;
In the present exemplary embodiment, understand that the textural characteristics in object block region and false-alarm block region are not by statistical analysis Together, and gray level co-occurrence matrixes can describe the textural characteristics of image.Therefore, the gray level co-occurrence matrixes in each block region are calculated, and Its characteristic information is counted, statistical analysis is carried out to a large amount of scene images, draw different characteristic information (such as contrast, energy and entropy) Threshold value, formulate criterion be the default criterion, so can further remove possible false-alarm block region, obtain Treat final object variations region in contrast images.
Step S402:When the contrast is more than the first threshold, the energy value less than the Second Threshold and The entropy meets more than the 3rd threshold value simultaneously when, the characteristic information meets the default criterion.
Exemplary, the first threshold can be 0.6, and the Second Threshold can be 0.12, and the 3rd threshold value can Think 2.1.So when the contrast being calculated is more than 0.6, the energy value is less than 0.12 and the entropy is more than 2.1 when meeting simultaneously, then it is assumed that a corresponding block region is false-alarm block region.
Step S303:The region of variation behind a block region as object variations region will be removed.Namely when Three above condition meets simultaneously when, corresponding false-alarm block region is removed.This embodiment scheme is easy and effective, it is easy to accomplish, solution Easily there is model modification difficulty, is easily influenceed by dynamic background interference using background modeling and difference in conventional method of having determined, improves The precision of scene Recognition and the degree of accuracy.
In order to verify the said system in this example embodiment in actual different scenes squatter building behavioral value it is accurate Property, inventor is tested the construction activity in two groups of different scenes in practice.
Fig. 5 A show the benchmark image in scene one, and Fig. 5 B show in scene one and treat contrast images.Fig. 5 C are field The change identification test result schematic diagram of scape one.Be can be seen that compared with benchmark image from Fig. 5 A and Fig. 5 B, treat that contrast images occur Crowd activity and the change of brick wall, additionally, also leaf is rocked and slight illumination variation.For practical application, crowd The region of variation of activity and brick wall needs to detect, and rocking for leaf, it is background region of variation, it is necessary to filter.Adopt The test result obtained with the said system in this example embodiment is as shown in Figure 5 C.From the results, it was seen that this example reality Rectangle frame region in applying the said system in mode and can detecting the region of variation of brick wall and groups of people in scene, i.e. Fig. 5 C, Rocking etc. for leaf is considered that false-alarm has been removed.
Fig. 6 A show the benchmark image in scene two, and Fig. 6 B show in scene two and treat contrast images.Fig. 6 C are field The change identification test result schematic diagram of scape two.Be can be seen that compared with benchmark image from Fig. 6 A and Fig. 6 B, treat that contrast images occur Crowd activity.Additionally, the change of calico and the rustle of leaves in the wind on also big ceiling.For practical application, crowd activity region It is to need the region of variation that detects, and the change of calico and the rustle of leaves in the wind that region is then background variation zone on big ceiling Domain is, it is necessary to filter.The testing result obtained using the said system in this example embodiment as shown in Figure 6 C, can be with from result Find out, the said system in this example embodiment can detect the rectangle frame area in the region of variation of people in scene, i.e. Fig. 6 C Domain, the change of calico and the rustle of leaves in the wind that region is considered that false-alarm is removed on big ceiling.
Additionally, in order to verify the robustness that the said system in this example embodiment changes to global illumination, to reality Two groups of images in scene three are tested.Fig. 7 A show the benchmark image in scene three, and Fig. 7 B are shown in scene three Treat contrast images.Fig. 7 C are the change identification test result schematic diagram of scene three.Be can be seen that and base from the contrast of Fig. 7 A and Fig. 7 B Quasi- image is compared, and treats that contrast images occur the change of brick wall and illumination, in addition with rocking for leaf.Using this example embodiment party The testing result that said system in formula is obtained as seen in figure 7 c, from the results, it was seen that above-mentioned in this example embodiment System can detect the multiple rectangle frame regions in the region of variation of brick wall in scene, i.e. Fig. 7 C.Although scene illumination there occurs Larger change, but the region of variation of brick wall still can be detected, there is not false-alarm.It can thus be seen that this example is implemented Said system in mode compares robust to global illumination change.
The embodiment of the present invention also has following beneficial technique effect:The embodiment of the present invention goes out from signal processing theory first Hair, the change occurred using the frequency domain correlated judgment image of signal, compared with traditional background modeling method, the embodiment of the present invention The change of region generation is concerned with, it is insensitive to the slight change in scene, solve conventional method easily dry by dynamic background Disturb the problem of influence.Secondly, the embodiment of the present invention is based on Fourier transformation theory, easy and effective, it is easy to which hardware is realized, is actual Using there is provided possible.In addition, the embodiment of the present invention using the gray level co-occurrence matrixes of target area calculate different characteristic values with Target area and false-alarm region are distinguished, for rear end provides preliminary early warning, can be greatly reduced based on the most of false-alarm of this removal Front-end image sends the bandwidth and cost pressure caused to rear end.Finally, the embodiment of the present invention can also be realized to scene change The coarse positioning in region, be rear end carry out further identification search space is greatly reduced, be that final quick judgement lays the foundation.
Although it should be noted that be referred in above-detailed the system for action executing some modules or Person's unit, but this division is not enforceable.In fact, according to embodiment of the present disclosure, above-described two or The feature and function of more multimode or unit can embody in a module or unit.Conversely, above-described one The feature and function of individual module or unit can be further divided into being embodied by multiple modules or unit.As module Or the part that shows of unit can be or may not be physical location, you can with positioned at a place, or can also divide Cloth is on multiple NEs.Some or all of module therein can be according to the actual needs selected to realize wooden disclosure side The purpose of case.Those of ordinary skill in the art are without creative efforts, you can to understand and implement.
Through the above description of the embodiments, those skilled in the art is it can be readily appreciated that example described herein is implemented Mode can be realized by software, it is also possible to be realized by way of software is with reference to necessary hardware.Therefore, according to the disclosure The technical scheme of implementation method can be embodied in the form of software product, and the software product can store non-volatile at one Property storage medium (can be CD-ROM, USB flash disk, mobile hard disk etc.) in or network on, including some instructions are causing a calculating Equipment (can be personal computer, server or network equipment etc.) performs the systemic-function according to disclosure implementation method.
Those skilled in the art will readily occur to its of the disclosure after considering specification and putting into practice invention disclosed herein Its embodiment.The application is intended to any modification, purposes or the adaptations of the disclosure, these modifications, purposes or Person's adaptations follow the general principle of the disclosure and including the undocumented common knowledge in the art of the disclosure Or conventional techniques.Description and embodiments are considered only as exemplary, and the true scope of the disclosure and spirit are by appended Claim is pointed out.

Claims (10)

1. it is a kind of for land resources monitoring squatter building scene Recognition system, it is characterised in that the system includes:
Image collection module, for obtaining the first two field picture of land resources region to be monitored predeterminated position as benchmark image, And the two field picture after first two field picture is obtained as treating contrast images;
Image block module, for the benchmark image and treating that contrast images carry out piecemeal treatment with by the reference map respectively Picture and treat that contrast images are divided into corresponding multiple block regions;
Image processing module, for corresponding a pair of block region to carry out Fourier's change to the benchmark image and in treating contrast images Crosspower spectrum is changed and be calculated, then inverse Fourier transform is carried out to the crosspower spectrum and obtained a spatial transform matrix and calculate The maximum of the spatial transform matrix;And
Area determination module, for obtaining the benchmark image and treating that corresponding each pair block region is corresponding described in contrast images Corresponding piece of region of maximum of the minimum spatial transform matrix of top n is used as described in the maximum of spatial transform matrix Treat the region of variation in contrast images;Wherein N is more than or equal to 2 and less than the described piece of integer of the total number in region.
2. squatter building scene Recognition system according to claim 1, it is characterised in that described image processing module, is used for:
The benchmark image is calculated according to below equation and i-th cross-power in corresponding a pair of block region in contrast images is treated Spectrum:
e j 2 π ( ux 0 + vy 0 ) = F 1 i ( u , v ) · F 2 i * ( u , v ) | F 1 i ( u , v ) · F 2 i * ( u , v ) | ;
Wherein, F1i(u, v) is f1iThe Fourier transformation of (x, y), f1i(x, y) represents i-th piece of region in the benchmark image;
F2i(u, v) is f2iThe Fourier transformation of (x, y), f2i(x, y) treats i-th piece of region in contrast images described in representing,It is F2iThe conjugation of (u, v);
x0,y0Respectively it is described treat in contrast images i-th piece of region relative to i-th piece of region in the benchmark image in x and The side-play amount in y directions.
3. squatter building scene Recognition system according to claim 1, it is characterised in that the value of the N is 5.
4. the squatter building scene Recognition system according to any one of claims 1 to 3, it is characterised in that described image piecemeal module, For:
To the benchmark image and treat that contrast images carry out uniform non-overlapping piecemeal treatment and form multiple block regions respectively.
5. squatter building scene Recognition system according to claim 4, it is characterised in that the benchmark image and treat contrast images Resolution ratio is 1920*1080, and the size in described piece of region is 120*120.
6. squatter building scene Recognition system according to claim 4, it is characterised in that each described piece of region has station location marker Information;The area determination module, is used for:
Obtain the benchmark image and treat that the corresponding spatial transform matrix in corresponding each pair block region is most in contrast images Corresponding piece of region of maximum of the spatial transform matrix of top n minimum and the N number of piece of position in region is recorded in big value Identification information;
Station location marker information according to the N number of piece of region using the N number of piece of area maps to it is described treat in contrast images as The region of variation.
7. squatter building scene Recognition system according to claim 6, it is characterised in that the system also judges mould including false-alarm Block, is used for:
The gray level co-occurrence matrixes in each the block region in the region of variation are calculated, and calculates the feature of the gray level co-occurrence matrixes Information;
Judge whether the characteristic information meets a default criterion, the corresponding change of this feature information is removed if meeting Change a block region in region;
The region of variation behind a block region as object variations region will be removed.
8. squatter building scene Recognition system according to claim 7, it is characterised in that the characteristic information of the gray level co-occurrence matrixes Including at least one of contrast, energy and entropy.
9. squatter building scene Recognition system according to claim 8, it is characterised in that the false-alarm judge module, is used for:
Whether whether the contrast is judged more than first threshold, and judges the energy value less than Second Threshold, while judging Whether the entropy is more than the 3rd threshold value;
When the contrast is more than the first threshold, the energy value is more than institute less than the Second Threshold and the entropy The 3rd threshold value is stated while when meeting, the characteristic information meets the default criterion.
10. squatter building scene Recognition system according to claim 9, it is characterised in that the first threshold is 0.6, described the Two threshold values are 0.12, and the 3rd threshold value is 2.1.
CN201710102202.0A 2017-02-24 2017-02-24 For the squatter building scene Recognition system of land resources monitoring Pending CN106910190A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710102202.0A CN106910190A (en) 2017-02-24 2017-02-24 For the squatter building scene Recognition system of land resources monitoring

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710102202.0A CN106910190A (en) 2017-02-24 2017-02-24 For the squatter building scene Recognition system of land resources monitoring

Publications (1)

Publication Number Publication Date
CN106910190A true CN106910190A (en) 2017-06-30

Family

ID=59208766

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710102202.0A Pending CN106910190A (en) 2017-02-24 2017-02-24 For the squatter building scene Recognition system of land resources monitoring

Country Status (1)

Country Link
CN (1) CN106910190A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112233169A (en) * 2020-10-27 2021-01-15 江苏科博空间信息科技有限公司 Based on unmanned aerial vehicle area of offending land calculation system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103034860A (en) * 2012-12-14 2013-04-10 南京思创信息技术有限公司 Scale-invariant feature transform (SIFT) based illegal building detection method
CN103491351A (en) * 2013-09-29 2014-01-01 东南大学 Intelligent video monitoring method for illegal buildings
CN103546726A (en) * 2013-10-28 2014-01-29 东南大学 Method for automatically discovering illegal land use

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103034860A (en) * 2012-12-14 2013-04-10 南京思创信息技术有限公司 Scale-invariant feature transform (SIFT) based illegal building detection method
CN103491351A (en) * 2013-09-29 2014-01-01 东南大学 Intelligent video monitoring method for illegal buildings
CN103546726A (en) * 2013-10-28 2014-01-29 东南大学 Method for automatically discovering illegal land use

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
景丽等: "基于相位相关的图像区域复制篡改检测与定位", 《计算机应用》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112233169A (en) * 2020-10-27 2021-01-15 江苏科博空间信息科技有限公司 Based on unmanned aerial vehicle area of offending land calculation system

Similar Documents

Publication Publication Date Title
CN106897694A (en) For the squatter building scene recognition method of land resources monitoring
Romeo et al. Crop row detection in maize fields inspired on the human visual perception
Shao et al. Application of a fast linear feature detector to road extraction from remotely sensed imagery
CN104766099B (en) Distribution type fiber-optic circumference vibration signal processing and recognition methods based on image
CN109029381A (en) A kind of detection method of tunnel slot, system and terminal device
EP2297701A2 (en) Method for video analysis
CN105335955A (en) Object detection method and object detection apparatus
CN104636764B (en) A kind of image latent writing analysis method and its device
CN106156758B (en) A kind of tidal saltmarsh method in SAR seashore image
CN101551852A (en) Training system, training method and detection method
CN103747271A (en) Video tamper detection method and device based on mixed perceptual hashing
CN104463104B (en) A kind of stationary vehicle target rapid detection method and device
CN112883802A (en) Method for identifying destructive event of pipeline optical fiber vibration safety early warning system
CN106295657A (en) A kind of method extracting human height's feature during video data structure
CN110348434A (en) Camera source discrimination method, system, storage medium and calculating equipment
CN108960012A (en) Feature point detecting method, device and electronic equipment
CN102629325B (en) Image characteristic extraction method, device thereof, image copy detection method and system thereof
CN113111706B (en) SAR target feature unwrapping and identifying method for azimuth continuous deletion
CN106910190A (en) For the squatter building scene Recognition system of land resources monitoring
CN106375756B (en) It is a kind of to remove the detection method distorted for the single object of monitor video
CN102024244A (en) Method and device for embedding and detecting watermarks based on image characteristic region
CN112529837A (en) Remote sensing image change detection algorithm based on coupling discrimination feature self-learning network
CN117058548A (en) Remote sensing image target detection method and device, storage medium and electronic equipment
CN116612097A (en) Method and system for predicting internal section morphology of wood based on surface defect image
CN111126508A (en) Hopc-based improved heterogeneous image matching method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20170630