CN101776759A - Remote sensing image-based area target motion velocity acquiring method and device - Google Patents

Remote sensing image-based area target motion velocity acquiring method and device Download PDF

Info

Publication number
CN101776759A
CN101776759A CN201010106736A CN201010106736A CN101776759A CN 101776759 A CN101776759 A CN 101776759A CN 201010106736 A CN201010106736 A CN 201010106736A CN 201010106736 A CN201010106736 A CN 201010106736A CN 101776759 A CN101776759 A CN 101776759A
Authority
CN
China
Prior art keywords
image
window
velocity
velocity field
target motion
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
CN201010106736A
Other languages
Chinese (zh)
Other versions
CN101776759B (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.)
CENTER FOR EARTH OBSERVATION AND DIGITAL EARTH CHINESE ACADEMY OF SCIENCES
Original Assignee
CENTER FOR EARTH OBSERVATION AND DIGITAL EARTH CHINESE ACADEMY OF SCIENCES
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 CENTER FOR EARTH OBSERVATION AND DIGITAL EARTH CHINESE ACADEMY OF SCIENCES filed Critical CENTER FOR EARTH OBSERVATION AND DIGITAL EARTH CHINESE ACADEMY OF SCIENCES
Priority to CN 201010106736 priority Critical patent/CN101776759B/en
Publication of CN101776759A publication Critical patent/CN101776759A/en
Application granted granted Critical
Publication of CN101776759B publication Critical patent/CN101776759B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

The embodiment of the invention relates to a remote sensing image-based area target motion velocity acquiring method and a remote sensing image-based area target motion velocity acquiring device. The method comprises the following steps: acquiring more than two velocity field changing factors respectively corresponding to matched windows with more than two sizes in a first image and a second image; acquiring the optimal window size from the windows with the more than two sizes in the first image according to the more than two velocity field changing factors; and calculating the area target motion velocity corresponding to the optimal window size. Through the remote sensing image-based area target motion velocity acquiring method and the remote sensing image-based area target motion velocity acquiring device, the optimal window size is acquired from the windows with more than two sizes in the first image through the more than two velocity field changing factors, and then the area target motion velocity corresponding to the optimal window size is calculated, so that the spacial resolution of the velocity field is improved, the gentle change of the velocity field is guaranteed and the accuracy of the area target motion velocity is improved.

Description

Area target motion velocity acquiring method and device based on remote sensing images
Technical field
The embodiment of the invention relates to the remote sensing technology field, especially a kind of area target motion velocity acquiring method and device based on remote sensing images.
Background technology
In recent years, frequently generation activity of the face of land, if too active can the causing in the face of land caused catastrophic consequence, for example: the rubble flow that cause the flood that the lake dam break causes, landslide etc., wherein area target motion is the key character that ground motion changes.Yet, at present area target motion is carried out but famine of data monitored, be difficult to according to Monitoring Data survey region target self characteristics, also increased the difficulty of the disaster alarm that causes by area target motion.
In the prior art, correlation analysis by remote sensing image is monitored area target motion velocity, particularly, image segmentation is become a plurality of graticule mesh, from a plurality of graticule mesh, choose a window, the textural characteristics in the window on two width of cloth images that do not get access to is mutually simultaneously followed the tracks of and the zoning target speed.In actual monitoring situation to area target motion velocity, because non-translational variation is bigger to the less grid influence of size, if non-translational variation takes place the regional aim in having among a small circle, cause the grid matching error, thereby greatly reduce the accuracy of area target motion velocity.
Summary of the invention
The purpose of the embodiment of the invention is to provide a kind of area target motion velocity acquiring method and device based on remote sensing images, improves the accuracy of zoning target speed.
The embodiment of the invention provides a kind of area target motion velocity acquiring method based on remote sensing images, comprising:
The window that is complementary that obtains two above sizes in first image and second image is the changed factor of corresponding plural velocity field respectively, and described first image and described second image are remote sensing images;
Changed factor according to described plural velocity field obtains best window size from the window of described two above sizes of described first image;
Calculate the corresponding area target motion velocity of window size of described the best.
The embodiment of the invention also provides a kind of area target motion velocity deriving means based on remote sensing images, comprising:
First acquisition module is used for obtaining at first image and second image changed factor of the respectively corresponding plural velocity field of the window that is complementary of two above sizes, and described first image and described second image are remote sensing images;
Second acquisition module is used for obtaining best window size according to the changed factor of described plural velocity field from the window of described two above sizes of described first image;
Computing module is used to calculate the corresponding area target motion velocity of window size of described the best.
Area target motion velocity acquiring method and device that the embodiment of the invention provides based on remote sensing images, changed factor by plural velocity field obtains best window size from the window of two above sizes of first image, the corresponding area target motion velocity of the window size of calculating optimum, improved the spatial resolution of velocity field, guarantee the smooth variation of velocity field, improved the accuracy of zoning target speed.
Description of drawings
In order to be illustrated more clearly in the embodiment of the invention or technical scheme of the prior art, to do to introduce simply to the accompanying drawing of required use in embodiment or the description of the Prior Art below, apparently, accompanying drawing in describing below only is some embodiments of the present invention, for those of ordinary skills, under the prerequisite of not paying creative work, can also obtain other accompanying drawing according to these accompanying drawings.
Fig. 1 is the schematic flow sheet that the present invention is based on an embodiment of area target motion velocity acquiring method of remote sensing images;
Fig. 2 is the schematic flow sheet that the present invention is based on another embodiment of area target motion velocity acquiring method of remote sensing images;
Fig. 3 is the synoptic diagram of window and region of search in the step 201 embodiment illustrated in fig. 2;
Fig. 4 is the curved line relation synoptic diagram of window size and FVV in the step 204 embodiment illustrated in fig. 2;
Fig. 5 is the structural representation that the present invention is based on an embodiment of area target motion velocity deriving means of remote sensing images;
Fig. 6 is the structural representation that the present invention is based on another embodiment of area target motion velocity deriving means of remote sensing images.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the invention, the technical scheme in the embodiment of the invention is clearly and completely described, obviously, described embodiment only is the present invention's part embodiment, rather than whole embodiment.Based on the embodiment among the present invention, those of ordinary skills belong to the scope of protection of the invention not making the every other embodiment that is obtained under the creative work prerequisite.
The correlation technique of remote sensing images is to resolve the effective means of area target motion velocity, by two width of cloth remote sensing images being carried out the window coupling, determines two coordinate offset amounts between the window, and then obtains the movement velocity of window area.
Be understandable that according to material balance principle, the movement velocity of regional aim is gradual change in a bigger scope, among a small circle Nei regional aim speed apparently higher than or be lower than on every side that the speed of adjacent domain target is irrational.For the face of land speed of regional aim in adjacent domain that obtains on the evaluation map picture, the changed factor of a velocity field of embodiment of the invention setting (Factor of Velocity Variation, be called for short: FVV), by FVV = Σ 1 N ( Ry 2 Rx 2 + Ry 2 ( V ( x , y ) - V ( x + n , y ) ) 2 + Rx 2 Rx 2 + Ry 2 ( V ( x , y ) - V ( x , y + n ) ) 2 ) N Calculate FVV, wherein, Rx, Ry are respectively remote sensing images pixel separation on directions X and the Y direction in the XY coordinate system,
Figure GSA00000016875400032
Be the weight of directions X velocity contrast,
Figure GSA00000016875400033
Be the weight of Y direction velocity contrast, wherein, weight is represented two pairs of adjacent windows in embodiments of the present invention, and under the identical situation of velocity contrast, the FVV value of the big more correspondence of pixel separation is more little; Window size is n * n pixel (n for more than or equal to 1 positive integer), V (x, y) be the center coordinate (x, the speed of the window of y) locating, V (x+n, y), (x y+n) is and coordinate (x, y) speed of the window of an adjacent n pixel V; N is the sum of the window in the regional aim that participates in the image calculating.In addition, regional aim described in the embodiment of the invention specifically can be to distinguish mutually with ground point target (for example: people, automobile etc.), specifically can refer to face of land active regions, for example: can refer to large-area object, large-area object is specially the large area region that melting of ice and snow takes place, the large area region that the landslide takes place etc.
Fig. 1 is the schematic flow sheet that the present invention is based on an embodiment of area target motion velocity acquiring method of remote sensing images, and as shown in Figure 1, the embodiment of the invention comprises the steps:
Step 101, the window that the is complementary changed factor of corresponding plural velocity field respectively that in first image and second image, obtains two above sizes, wherein, first image and second image are remote sensing images;
Step 102, from the window of two above sizes of first image, obtain best window size according to the changed factor of plural velocity field;
The corresponding area target motion velocity of the window size of step 103, calculating optimum.
In the present embodiment, first image (being also referred to as master image) and second image (being also referred to as from image) are the same area target, mutually remote sensing image data simultaneously not, in the embodiment of the invention the time specifically refer to the shooting time of first image and the shooting time of second image mutually, for example: first image is the remote sensing images of taking on 01 06th, 2007 western part of China Tianshan Mountains regional aim, second image is the remote sensing images of taking on 02 06th, 2008 western part of China Tianshan Mountains regional aim, and then the shooting time of first image and second image is spaced apart 13 months.
The area target motion velocity acquiring method that the embodiment of the invention provides based on remote sensing images, changed factor by plural velocity field obtains best window size from the window of two above sizes of first image, the corresponding area target motion velocity of the window size of calculating optimum, improved the spatial resolution of velocity field, guarantee the smooth variation of velocity field, improved the accuracy of zoning target speed.
Fig. 2 is the schematic flow sheet that the present invention is based on another embodiment of area target motion velocity acquiring method of remote sensing images, Fig. 3 be embodiment illustrated in fig. 2 in the synoptic diagram of window and region of search in the step 201, Fig. 4 is the curved line relation synoptic diagram of window size and FVV in the step 204 embodiment illustrated in fig. 2; In the present embodiment, first image (being also referred to as master image) and second image (being also referred to as from image) are the same area target, mutually remote sensing image data simultaneously not, in the embodiment of the invention the time specifically refer to the shooting time of first image and the shooting time of second image mutually, for example: first image is the remote sensing images of taking on 01 06th, 2007 western part of China Tianshan Mountains regional aim, second image is the remote sensing images of taking on 02 06th, 2008 western part of China Tianshan Mountains regional aim, and then the shooting time interval T of first image and second image is 13 months; As shown in Figure 2, the embodiment of the invention comprises the steps:
Step 201, calculate the window of two above sizes in first image and the related coefficient of corresponding region of search on second image respectively, wherein, corresponding related coefficient set in the region of search of the window of each size in second image;
Wherein, on first image, obtain a plurality of sizes window inequality, find on second image with first image on a plurality of sizes window inequality of being complementary respectively of a plurality of sizes window inequality; In the region of search of second image, calculate the related coefficient of the window corresponding with each size on first image, because default region of search is greater than the size of window, therefore when window slides the calculating related coefficient in the region of search, sliding each time to calculate all to obtain a corresponding related coefficient, the therefore all corresponding related coefficient set in the region of search of window in second image of each size on first image; For example: the window that on first image, obtains six sizes, resolution is respectively 10 * 10,20 * 20,30 * 30,40 * 40,50 * 50,60 * 60, when the size on calculating first image is the related coefficient set on second image of 10 * 10 window, then can on second image, set size and be 15 * 15 region of search; As shown in Figure 3, can pass through similarity index (Similarity Index is called for short: SI)) SI ( u , v ) = Σ y = 1 M Σ x = 1 N ( f ( x , y ) - f ‾ ) ( g ( x + u , y + v ) - g ‾ ( u , v ) ) ( Σ y = 1 M Σ x = 1 N ( f ( x , y ) - f ‾ ) 2 ) 1 / 2 ( Σ y = 1 M Σ x = 1 N ( g ( x + u , y + v ) - g ‾ ( u , v ) ) 2 ) 1 / 2 Calculate the related coefficient described in the present embodiment, f (x, y) be the gray-scale value of the pixel in the window (present embodiment is for convenience of description with this window called after window P) on first image, (x y) is the gray-scale value of window (present embodiment is for convenience of description with this window called after window P ') interior pixel that is complementary with this window P on second image to g.U, v is window P on first image and the coordinate offset amount between the window P ' on second image, f, g (u, v) be respectively the window P on first image and the average gray value of the window P ' interior pixel on second image, N, M are the size of window, by in the region of search, calculating the related coefficient a plurality of and window that first window is complementary, gather thereby further obtain related coefficient.
Step 202, in the set of each related coefficient, obtain the related coefficient of numerical value maximum, with the corresponding window of the related coefficient of the maximum on second image as with first image on the corresponding match window of window;
Step 203, in plural match window the changed factor of the velocity field of zoning target respectively, obtain the changed factor of plural velocity field;
Particularly, FVV = Σ 1 N ( Ry 2 Rx 2 + Ry 2 ( V ( x , y ) - V ( x + n , y ) ) 2 + Rx 2 Rx 2 + Ry 2 ( V ( x , y ) - V ( x , y + n ) ) 2 ) N , Wherein, Rx, Ry are respectively first image pixel separation on directions X and the Y direction in the XY coordinate system, Be the weight of directions X velocity contrast,
Figure GSA00000016875400063
Be the weight of Y direction velocity contrast, n represents the size of described window, for more than or equal to 1 positive integer), V (x, y) be the center coordinate (x, the speed of the window of y) locating, V (x+n, y), (x y+n) is and coordinate (x, y) speed of the window of an adjacent n pixel V; N is the sum that participates in the window of calculating in the image.
Step 204, in the coordinate system that changed factor and window size by velocity field form according to the window size order from small to large of two above sizes of first image successively every the set point number fitting a straight line;
Step 205, when the curvature of determining straight line when setting threshold value, then will set the window size of the pairing window size of threshold value as the best.
As shown in Figure 4, in above-mentioned steps 204 and the step 205, in the coordinate system that the changed factor of velocity field and window size form, when window size is big more, the spatial resolution of corresponding velocity field will reduce, and for example: window size is that 50 * 50 pixels and window size are that the spatial resolution of velocity field of 30 * 30 pixels is low.Thus, be optimal window size by acute variation to the window size of the turning point correspondence of smooth variation with the FVV value, both can guarantee that the degree that velocity field changes is less, can guarantee the spatial resolution of velocity field again; Being specially three points with set point number is that example describes, as shown in Figure 4, on curve, be zero beginning by window size, three points of every vicinity match straight line, when straight line curvature absolute value during less than predetermined threshold value 0.2, thinking that curve begins to become gently, is 0.2 the corresponding window size of the coordinate points window size as the best with this straight line curvature.
The corresponding area target motion velocity of the window size of step 206, calculating optimum;
Particularly, can be according to formula
Figure GSA00000016875400071
The corresponding area target motion velocity of the window size of calculating optimum; Wherein, Vh is the regional aim plane speed on first image, Rx, Ry is first image pixel separation on directions X and the Y direction in the XY coordinate system, T is the time interval of first image and second image taking, Δ x is the coordinate offset amount of window center coordinate on directions X of second image and first image, and Δ y is the coordinate offset amount of window center coordinate on the Y direction of second image and first image.
The area target motion velocity acquiring method that the embodiment of the invention provides based on remote sensing images, changed factor by plural velocity field obtains best window size from the window of two above sizes of first image, the corresponding area target motion velocity of the window size of calculating optimum, improved the spatial resolution of velocity field, guarantee the smooth variation of velocity field, improved the accuracy of zoning target speed.
Fig. 5 is the structural representation that the present invention is based on an embodiment of area target motion velocity deriving means of remote sensing images, and the embodiment of the invention can realize the method flow of above-mentioned embodiment illustrated in figures 1 and 2; As shown in Figure 5, the embodiment of the invention comprises: first acquisition module 51, second acquisition module 52, computing module 53;
Wherein, first acquisition module 51 obtains two above sizes in first image and second image the window that is complementary is the changed factor of corresponding plural velocity field respectively, and wherein, described first image and described second image are remote sensing images; Second acquisition module 52 obtains best window size according to the changed factor of described plural velocity field from the window of described two above sizes of described first image; Computing module 53 calculates the corresponding area target motion velocity of window size of described the best.
The area target motion velocity deriving means that the embodiment of the invention provides based on remote sensing images, second acquisition module 52 obtains best window size by the changed factor of plural velocity field from the window of two above sizes of first image, the corresponding area target motion velocity of window size of computing module 53 calculating optimums, improved the spatial resolution of velocity field, guarantee the smooth variation of velocity field, improved the accuracy of zoning target speed.
Fig. 6 is the structural representation that the present invention is based on another embodiment of area target motion velocity deriving means of remote sensing images, and the embodiment of the invention can realize the method flow of above-mentioned embodiment illustrated in figures 1 and 2; As shown in Figure 6, the embodiment of the invention comprises: first acquisition module 61, second acquisition module 62, computing module 63;
Wherein, first acquisition module 61 obtains two above sizes in first image and second image the window that is complementary is the changed factor of corresponding plural velocity field respectively, and wherein, first image and described second image are remote sensing images; Second acquisition module 62 obtains best window size according to the changed factor of described plural velocity field from the window of described two above sizes of described first image; Computing module 63 calculates the corresponding area target motion velocity of window size of described the best.
Further, first acquisition module 61 can also comprise: first computing unit 611, first acquiring unit 612, second computing unit 613; Wherein, first computing unit 611 calculates the window of two above sizes in described first image and the related coefficient of corresponding region of search target on described second image respectively, wherein, corresponding related coefficient set in the region of search target of the window of each size in described second image; First acquiring unit 612 obtains the related coefficient of numerical value maximum in the set of each related coefficient, with the corresponding window of related coefficient of the described maximum on described second image as with described first image on the corresponding match window of window; Second computing unit 613 is distinguished the changed factor of the velocity field of zoning target in plural match window, obtain the changed factor of plural velocity field.
Further, the changed factor of velocity field is FVV = Σ 1 N ( Ry 2 Rx 2 + Ry 2 ( V ( x , y ) - V ( x + n , y ) ) 2 + Rx 2 Rx 2 + Ry 2 ( V ( x , y ) - V ( x , y + n ) ) 2 ) N , Wherein, FVV is described velocity field changed factor, and Rx, Ry are respectively described master image pixel separation on directions X and the Y direction in the XY coordinate system, Be the weight of described directions X velocity contrast,
Figure GSA00000016875400093
Be the weight of described Y direction velocity contrast, n represents the size of described window, for more than or equal to 1 positive integer), V (x, y) be the center coordinate (x, the speed of the window of y) locating, V (x+n, y), (x y+n) is and coordinate (x, y) speed of the window of an adjacent n pixel V; N is the sum of the window in the regional aim that participates in the image calculating.
Further, second acquisition module 62 can also comprise: fitting a straight line unit 621 and determining unit 622; Wherein, fitting a straight line unit 621 in the coordinate system that forms by described velocity field changed factor and described window size according to the window size order from small to large of described two above sizes of described first image successively every the set point number fitting a straight line; The curvature of determining described straight lines when determining unit 622 is when setting threshold value, then with the window size of the pairing window size of described setting threshold value as the best.
Further, computing module 63 is according to formula
Figure GSA00000016875400094
Calculate the corresponding area target motion velocity of window size of described the best, wherein, Vh is the regional aim plane speed on described first image, Rx, Ry is described first image pixel separation on directions X and the Y direction in the XY coordinate system, T is the time interval of described first image and described second image taking, Δ x is the coordinate offset amount of window center coordinate on described directions X of described second image and described first image, and Δ y is the coordinate offset amount of window center coordinate on described Y direction of described second image and described first image.
The area target motion velocity deriving means that the embodiment of the invention provides based on remote sensing images, second acquisition module 62 obtains best window size by the changed factor of plural velocity field from the window of two above sizes of first image, the corresponding area target motion velocity of window size of computing module 63 calculating optimums, improved the spatial resolution of velocity field, guarantee the smooth variation of velocity field, improved the accuracy of zoning target speed.
The those skilled in the art can be well understood to, and is the convenience described and succinct, and the concrete course of work of the system of foregoing description, equipment, module and unit can not repeat them here with reference to the corresponding process among the preceding method embodiment.
One of ordinary skill in the art will appreciate that: all or part of step that realizes the foregoing description can be finished by the relevant hardware of programmed instruction, aforesaid program can be stored in the computer read/write memory medium, this program is carried out the step that comprises said method embodiment when carrying out; And aforesaid storage medium comprises: various media that can be program code stored such as ROM, RAM, magnetic disc or CD.
It should be noted that at last: above embodiment only in order to technical scheme of the present invention to be described, is not intended to limit; Although with reference to previous embodiment the present invention is had been described in detail, those of ordinary skill in the art is to be understood that: it still can be made amendment to the technical scheme that aforementioned each embodiment put down in writing, and perhaps part technical characterictic wherein is equal to replacement; And these modifications or replacement do not make the essence of appropriate technical solution break away from the spirit and scope of various embodiments of the present invention technical scheme.

Claims (10)

1. the area target motion velocity acquiring method based on remote sensing images is characterized in that, comprising:
The window that is complementary that obtains two above sizes in first image and second image is the changed factor of corresponding plural velocity field respectively, and described first image and described second image are remote sensing images;
Changed factor according to described plural velocity field obtains best window size from the window of described two above sizes of described first image;
Calculate the corresponding area target motion velocity of window size of described the best.
2. method according to claim 1 is characterized in that, the described corresponding respectively plural velocity field changed factor of the window that is complementary that obtains two above sizes in first image and second image comprises:
Calculate the window of two above sizes in described first image and the related coefficient of corresponding region of search on described second image respectively, wherein, corresponding related coefficient set in the region of search of the window of each size in described second image;
In the set of each related coefficient, obtain the related coefficient of numerical value maximum, with the corresponding window of related coefficient of the described maximum on described second image as with described first image on the corresponding match window of window;
In plural match window, distinguish the changed factor of the velocity field of zoning target, obtain the changed factor of plural velocity field.
3. method according to claim 2 is characterized in that the changed factor of described velocity field is FVV = Σ 1 N ( Ry 2 Rx 2 + Ry 2 ( V ( x , y ) - V ( x + n , y ) ) 2 + Rx 2 Rx 2 + Ry 2 ( V ( x , y ) - V ( x , y + n ) ) 2 ) N , Wherein, FVV is described velocity field changed factor, and Rx, Ry are respectively described first image pixel separation on directions X and the Y direction in the XY coordinate system,
Figure FSA00000016875300012
Be the weight of described directions X velocity contrast,
Figure FSA00000016875300013
Be the weight of described Y direction velocity contrast, n represents the size of described window, for more than or equal to 1 positive integer, V (x, y) be the center coordinate (x, the speed of the window of y) locating, V (x+n, y), (x y+n) is and coordinate (x, y) speed of the window of an adjacent n pixel V; N is the sum of the window in the regional aim that participates in the image calculating.
4. method according to claim 1 is characterized in that, describedly obtains best window size according to described plural velocity field changed factor from the window of described two above sizes of described first image and comprises:
In the coordinate system that forms by described velocity field changed factor and described window size according to the window size order from small to large of described two above sizes of described first image successively every the set point number fitting a straight line;
When the curvature of determining described straight line when setting threshold value, then with the window size of the pairing window size of described setting threshold value as the best.
5. according to the arbitrary described method of claim 1~4, it is characterized in that the corresponding area target motion velocity of window size of the described the best of described calculating comprises:
According to formula
Figure FSA00000016875300021
Calculate the corresponding area target motion velocity of window size of described the best, wherein, Vh is the regional aim plane speed on described first image, Rx, Ry is described first image pixel separation on directions X and the Y direction in the XY coordinate system, T is the time interval of described first image and described second image taking, Δ x is the coordinate offset amount of window center coordinate on described directions X of described second image and described first image, and Δ y is the coordinate offset amount of window center coordinate on described Y direction of described second image and described first image.
6. the area target motion velocity deriving means based on remote sensing images is characterized in that, comprising:
First acquisition module is used for obtaining at first image and second image changed factor of the respectively corresponding plural velocity field of the window that is complementary of two above sizes, and described first image and described second image are remote sensing images;
Second acquisition module is used for obtaining best window size according to the changed factor of described plural velocity field from the window of described two above sizes of described first image;
Computing module is used to calculate the corresponding area target motion velocity of window size of described the best.
7. device according to claim 6 is characterized in that, described first acquisition module comprises:
First computing unit, be used for calculating respectively the related coefficient of the window of two above sizes of described first image and corresponding region of search target on described second image, wherein, corresponding related coefficient set in the region of search target of the window of each size in described second image;
First acquiring unit is used for obtaining in the set of each related coefficient the related coefficient of numerical value maximum, with the corresponding window of related coefficient of the described maximum on described second image as with described first image on the corresponding match window of window;
Second computing unit is used for the changed factor in the velocity field of plural match window difference zoning target, obtains plural velocity field changed factor.
8. device according to claim 7 is characterized in that the changed factor of described velocity field is FVV = Σ 1 N ( Ry 2 Rx 2 + Ry 2 ( V ( x , y ) - V ( x + n , y ) ) 2 + Rx 2 Rx 2 + Ry 2 ( V ( x , y ) - V ( x , y + n ) ) 2 ) N , Wherein, FVV is described velocity field changed factor, and Rx, Ry are respectively described first image pixel separation on directions X and the Y direction in the XY coordinate system, Be the weight of described directions X velocity contrast,
Figure FSA00000016875300033
Be the weight of described Y direction velocity contrast, n represents the size of described window, for more than or equal to 1 positive integer, V (x, y) be the center coordinate (x, the speed of the window of y) locating, V (x+n, y), (x y+n) is and coordinate (x, y) speed of the window of an adjacent n pixel V; N is the sum of the window in the regional aim that participates in the image calculating.
9. device according to claim 6 is characterized in that, described second acquisition module comprises:
The fitting a straight line unit is used at the coordinate system that is formed by described velocity field changed factor and described window size according to the window size order from small to large of described two above sizes of described first image successively every the set point number fitting a straight line;
Determining unit, be used for when the curvature of determining described straight line when setting threshold value, then with the window size of the pairing window size of described setting threshold value as the best.
10. according to the arbitrary described device of claim 6~9, it is characterized in that described computing module is used for according to formula
Figure FSA00000016875300041
Calculate the corresponding area target motion velocity of window size of described the best, wherein, Vh is the regional aim plane speed on described first image, Rx, Ry is described first image pixel separation on directions X and the Y direction in the XY coordinate system, T is the time interval of described first image and described second image taking, Δ x is the coordinate offset amount of window center coordinate on described directions X of described second image and described first image, and Δ y is the coordinate offset amount of window center coordinate on described Y direction of described second image and described first image.
CN 201010106736 2010-02-03 2010-02-03 Remote sensing image-based area target motion velocity acquiring method and device Expired - Fee Related CN101776759B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010106736 CN101776759B (en) 2010-02-03 2010-02-03 Remote sensing image-based area target motion velocity acquiring method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010106736 CN101776759B (en) 2010-02-03 2010-02-03 Remote sensing image-based area target motion velocity acquiring method and device

Publications (2)

Publication Number Publication Date
CN101776759A true CN101776759A (en) 2010-07-14
CN101776759B CN101776759B (en) 2012-06-13

Family

ID=42513266

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010106736 Expired - Fee Related CN101776759B (en) 2010-02-03 2010-02-03 Remote sensing image-based area target motion velocity acquiring method and device

Country Status (1)

Country Link
CN (1) CN101776759B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101980301A (en) * 2010-10-28 2011-02-23 北京智安邦科技有限公司 Method and device for acquiring movement speed of target in video image
CN114026436A (en) * 2019-06-25 2022-02-08 索尼集团公司 Image processing apparatus, image processing method, and program

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101980301A (en) * 2010-10-28 2011-02-23 北京智安邦科技有限公司 Method and device for acquiring movement speed of target in video image
CN101980301B (en) * 2010-10-28 2012-08-22 北京智安邦科技有限公司 Method and device for acquiring movement speed of target in video image
CN114026436A (en) * 2019-06-25 2022-02-08 索尼集团公司 Image processing apparatus, image processing method, and program

Also Published As

Publication number Publication date
CN101776759B (en) 2012-06-13

Similar Documents

Publication Publication Date Title
CN102768022B (en) Tunnel surrounding rock deformation detection method adopting digital camera technique
CN103458261B (en) Video scene variation detection method based on stereoscopic vision
CN103810717B (en) A kind of human body behavioral value method and device
CN103645476B (en) A kind of space-time homogeneity filtering method of synthetic-aperture radar differential interferometry graphic sequence
CN104599502A (en) Method for traffic flow statistics based on video monitoring
CN104820996A (en) Target tracking method based on self-adaptive blocks of video
CN105654428A (en) Method and system for image noise reduction
CN103106658A (en) Island or reef coastline rapid obtaining method
CN106372598A (en) Image stabilizing method based on image characteristic detection for eliminating video rotation and jittering
CN101582171A (en) Method and device for creating depth maps
CN103268496A (en) Target identification method of SAR (synthetic aperture radar) images
CN104680553A (en) Bilateral filtration-based variable kernel function target tracing method
CN107464256A (en) A kind of target detection differentiates the correlating method of amendment with possibility
CN105427333A (en) Real-time registration method of video sequence image, system and shooting terminal
CN103985128A (en) Three-dimensional matching method based on color intercorrelation and self-adaptive supporting weight
CN103824305A (en) Improved Meanshift target tracking method
CN102129559B (en) SAR (Synthetic Aperture Radar) image object detection method based on Primal Sketch algorithm
CN108596032B (en) Detection method, device, equipment and medium for fighting behavior in video
CN103927553A (en) Coal and rock recognition method based on multi-scale micro-lamination and contrast ratio joint distribution
CN101776759B (en) Remote sensing image-based area target motion velocity acquiring method and device
CN104484865A (en) Method for removing raindrops in video image
CN103345752B (en) A kind of robot cooperates with mobile phone the method for tracking pedestrians
CN103093472A (en) Optical remote sensing image change detection method based on double dictionary cross sparse representation
CN105425296B (en) Geologic body recognition methods and device
CN105184809A (en) Moving object detection method and moving object detection device

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120613

Termination date: 20140203