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

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 through 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 in first image and second image, obtains two above sizes is the changed factor of corresponding plural velocity field respectively, and the changed factor of said velocity field does $\mathrm{FVV}=\frac{\underset{1}{\overset{N}{\mathrm{\Σ}}}\sqrt{(\frac{{\mathrm{Ry}}^2}{{\mathrm{Rx}}^2+{\mathrm{Ry}}^2}{(V(x,y)-V(x+n,y))}^2+\frac{{\mathrm{Rx}}^2}{{\mathrm{Rx}}^2+{\mathrm{Ry}}^2}{(V(x,y)-V(x,y+n))}^2)}}{N},$ Wherein, FVV is said velocity field changed factor, and Rx, Ry are respectively said first image pixel separation on directions X and the Y direction in the XY coordinate system,

Be the weight of said directions X velocity contrast,

Be the weight of said Y direction velocity contrast, n representes the size of said 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 of participating in the image calculating, and said first image and said second image are remote sensing images;

Changed factor according to said plural velocity field obtains best window size from the window of said two above sizes of said first image;

Calculate the corresponding area target motion velocity of window size of said 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 the changed factor of said velocity field does $\mathrm{FVV}=\frac{\underset{1}{\overset{N}{\mathrm{\Σ}}}\sqrt{(\frac{{\mathrm{Ry}}^2}{{\mathrm{Rx}}^2+{\mathrm{Ry}}^2}{(V(x,y)-V(x+n,y))}^2+\frac{{\mathrm{Rx}}^2}{{\mathrm{Rx}}^2+{\mathrm{Ry}}^2}{(V(x,y)-V(x,y+n))}^2)}}{N},$ Wherein, FVV is said velocity field changed factor, and Rx, Ry are respectively said first image pixel separation on directions X and the Y direction in the XY coordinate system, Be the weight of said directions X velocity contrast,

Be the weight of said Y direction velocity contrast, n representes the size of said 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 of participating in the image calculating, and said first image and said second image are remote sensing images;

Second acquisition module is used for obtaining best window size according to the changed factor of said plural velocity field from the window of said two above sizes of said first image;

Computing module is used to calculate the corresponding area target motion velocity of window size of said the best.

Area target motion velocity acquiring method and device that the embodiment of the invention provides based on remote sensing images; Changed factor through 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; Improve the spatial resolution of velocity field, guaranteed 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; Obviously, the 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 property, 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

To combine the accompanying drawing in the embodiment of the invention below, the technical scheme in the embodiment of the invention is carried out clear, intactly description, 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 are not making the every other embodiment that is obtained under the creative work prerequisite, all belong to the scope of the present invention's protection.

The correlation technique of remote sensing images is to resolve the effective means of area target motion velocity, through two width of cloth remote sensing images being carried out the window coupling, confirms two coordinate offset amounts between the window, and then obtains the movement velocity of window area.

It is understandable that according to material balance principle, the movement velocity of regional aim is gradual change in a bigger scope, the regional aim speed among a small circle 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), through $\mathrm{FVV}=\frac{\underset{1}{\overset{N}{\mathrm{\Σ}}}\sqrt{(\frac{{\mathrm{Ry}}^2}{{\mathrm{Rx}}^2+{\mathrm{Ry}}^2}{(V(x,y)-V(x+n,y))}^2+\frac{{\mathrm{Rx}}^2}{{\mathrm{Rx}}^2+{\mathrm{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,

Be the weight of directions X velocity contrast,

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 of participating 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 through 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; Improve the spatial resolution of velocity field, guaranteed 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 with first image on the related coefficient of the corresponding window of each size; Because preset 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)) $\mathrm{SI}(u,v)=\frac{\underset{y=1}{\overset{M}{\mathrm{\Σ}}}\underset{x=1}{\overset{N}{\mathrm{\Σ}}}(f(x,y)-\stackrel{\‾}{f})(g(x+u,y+v)-\stackrel{\‾}{g}(u,v))}{{\left(\underset{y=1}{\overset{M}{\mathrm{\Σ}}}\underset{x=1}{\overset{N}{\mathrm{\Σ}}}{(f(x,y)-\stackrel{\‾}{f})}^2\right)}^{1/2}{\left(\underset{y=1}{\overset{M}{\mathrm{\Σ}}}\underset{x=1}{\overset{N}{\mathrm{\Σ}}}{(g(x+u,y+v)-\stackrel{\‾}{g}(u,v))}^2\right)}^{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 for describe convenient with this window called after window P) on first image; (x y) is the gray-scale value of the window that is complementary with this window P on second image (present embodiment is convenient with this window called after window P ' for describing) interior pixel to g.U, v are window P and the coordinate offset amounts between the window P ' on second image on first image, f; G (u; V) be respectively window P and the average gray value of the window P ' interior pixel on second image on first image, N, M are the size of window; Through 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 maximum related coefficient of numerical value, 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, $\mathrm{FVV}=\frac{\underset{1}{\overset{N}{\mathrm{\Σ}}}\sqrt{(\frac{{\mathrm{Ry}}^2}{{\mathrm{Rx}}^2+{\mathrm{Ry}}^2}{(V(x,y)-V(x+n,y))}^2+\frac{{\mathrm{Rx}}^2}{{\mathrm{Rx}}^2+{\mathrm{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,

Be the weight of Y direction velocity contrast, n representes the size of said 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 participating in the window of calculating in the image.

Step 204, every successively according to the window size order from small to large of two above sizes of first image in the coordinate system that changed factor and window size by velocity field form at a distance from the set point number fitting a straight line;

Step 205, when the curvature of confirming 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 was 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 corresponding window size of the turning point of smooth variation with the FVV value, both can guarantee that the degree of velocity field variation 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 is 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 the corresponding area target motion velocity of window size of formula calculating optimum; Wherein, Vh is the regional aim plane speed on first image; Rx, Ry are first image pixel separation on directions X and the Y direction in the XY coordinate system, and 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 through 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; Improve the spatial resolution of velocity field, guaranteed 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, said first image and said second image are remote sensing images; Second acquisition module 52 obtains best window size according to the changed factor of said plural velocity field from the window of said two above sizes of said first image; Computing module 53 calculates the corresponding area target motion velocity of window size of said 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 through 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; Improve the spatial resolution of velocity field, guaranteed 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 said second image are remote sensing images; Second acquisition module 62 obtains best window size according to the changed factor of said plural velocity field from the window of said two above sizes of said first image; Computing module 63 calculates the corresponding area target motion velocity of window size of said 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 said first image and the related coefficient of corresponding region of search target on said second image respectively; Wherein, corresponding related coefficient set in the region of search target of the window of each size in said second image; First acquiring unit 612 obtains the maximum related coefficient of numerical value in the set of each related coefficient, with the corresponding window of related coefficient of the said maximum on said second image as with said 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 does $\mathrm{FVV}=\frac{\underset{1}{\overset{N}{\mathrm{\Σ}}}\sqrt{(\frac{{\mathrm{Ry}}^2}{{\mathrm{Rx}}^2+{\mathrm{Ry}}^2}{(V(x,y)-V(x+n,y))}^2+\frac{{\mathrm{Rx}}^2}{{\mathrm{Rx}}^2+{\mathrm{Ry}}^2}{(V(x,y)-V(x,y+n))}^2)}}{N},$ Wherein, FVV is said velocity field changed factor, and Rx, Ry are respectively said master image pixel separation on directions X and the Y direction in the XY coordinate system,

Be the weight of said directions X velocity contrast, Be the weight of said Y direction velocity contrast, n representes the size of said 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 of participating in the image calculating.

Further, second acquisition module 62 can also comprise: fitting a straight line unit 621 and definite unit 622; Wherein, fitting a straight line unit 621 is every successively at a distance from the set point number fitting a straight line according to the window size order from small to large of said two above sizes of said first image in the coordinate system that is formed by said velocity field changed factor and said window size; When confirming curvature that said straight lines are confirmed in unit 622 when setting threshold value, then with the window size of the pairing window size of said setting threshold value as the best.

Further; Computing module 63 calculates the corresponding area target motion velocity of window size of said the best according to formula

; Wherein, Vh is the regional aim plane speed on said first image; Rx; Ry is said first image pixel separation on directions X and the Y direction in the XY coordinate system; T is the time interval of said first image and said second image taking, and Δ x is the coordinate offset amount of window center coordinate on said directions X of said second image and said first image, and Δ y is the coordinate offset amount of window center coordinate on said Y direction of said second image and said 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 through 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; Improve the spatial resolution of velocity field, guaranteed the smooth variation of velocity field, improved the accuracy of zoning target speed.

The those skilled in the art can be well understood to, for the convenience described with succinct, the concrete course of work of the system of foregoing description, equipment, module and unit can repeat no more at this 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 accomplished through the relevant hardware of programmed instruction; Aforesaid program can be stored in the computer read/write memory medium; This program 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.

What should explain at last is: above embodiment is only in order to explaining technical scheme of the present invention, but not to its restriction; Although with reference to previous embodiment the present invention has been carried out detailed explanation, 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 are revised or replacement, do not make the spirit and the scope of the essence disengaging various embodiments of the present invention technical scheme of relevant art scheme.

Claims (8)

1. the area target motion velocity acquiring method based on remote sensing images is characterized in that, comprising:

The window that is complementary that in first image and second image, obtains two above sizes is the changed factor of corresponding plural velocity field respectively, and the changed factor of said velocity field does $\mathrm{FVV}=\frac{\underset{1}{\overset{N}{\mathrm{\Σ}}}\sqrt{(\frac{{\mathrm{Ry}}^2}{{\mathrm{Rx}}^2+{\mathrm{Ry}}^2}{(V(x,y)-V(x+n,y))}^2+\frac{{\mathrm{Rx}}^2}{{\mathrm{Rx}}^2+{\mathrm{Ry}}^2}{(V(x,y)-V(x,y+n))}^2)}}{N},$ Wherein, FVV is said velocity field changed factor, and Rx, Ry are respectively said first image pixel separation on directions X and the Y direction in the XY coordinate system,

Be the weight of said directions X velocity contrast,

Be the weight of said Y direction velocity contrast, n representes the size of said 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 of participating in the image calculating, and said first image and said second image are remote sensing images;

Changed factor according to said plural velocity field obtains best window size from the window of said two above sizes of said first image;

Calculate the corresponding area target motion velocity of window size of said the best.

2. method according to claim 1 is characterized in that, the said corresponding respectively plural velocity field changed factor of the window that is complementary that in first image and second image, obtains two above sizes comprises:

Calculate the window of two above sizes in said first image and the related coefficient of corresponding region of search on said second image respectively, wherein, corresponding related coefficient set in the region of search of the window of each size in said second image;

In the set of each related coefficient, obtain the maximum related coefficient of numerical value, with the corresponding window of related coefficient of the said maximum on said second image as with said 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 1 is characterized in that, saidly from the window of said two above sizes of said first image, obtains best window size according to said plural velocity field changed factor and comprises:

Every successively according to the window size order from small to large of said two above sizes of said first image in the coordinate system that forms by said velocity field changed factor and said window size at a distance from the set point number fitting a straight line;

When the curvature of confirming said fitting a straight line when setting threshold value, then with the window size of the corresponding window size of the coordinate points of said setting threshold value as the best.

4. according to each described method in the claim 1～3, it is characterized in that the corresponding area target motion velocity of window size of the said the best of said calculating comprises:

Calculate the corresponding area target motion velocity of window size of said the best according to formula

; Wherein, Vh is the regional aim plane speed on said first image; Rx; Ry is said first image pixel separation on directions X and the Y direction in the XY coordinate system; T is the time interval of said first image and said second image taking, and Δ x is the coordinate offset amount of window center coordinate on said directions X of said second image and said first image, and Δ y is the coordinate offset amount of window center coordinate on said Y direction of said second image and said first image.

5. 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 the changed factor of said velocity field does $\mathrm{FVV}=\frac{\underset{1}{\overset{N}{\mathrm{\Σ}}}\sqrt{(\frac{{\mathrm{Ry}}^2}{{\mathrm{Rx}}^2+{\mathrm{Ry}}^2}{(V(x,y)-V(x+n,y))}^2+\frac{{\mathrm{Rx}}^2}{{\mathrm{Rx}}^2+{\mathrm{Ry}}^2}{(V(x,y)-V(x,y+n))}^2)}}{N},$ Wherein, FVV is said velocity field changed factor, and Rx, Ry are respectively said first image pixel separation on directions X and the Y direction in the XY coordinate system,

Be the weight of said directions X velocity contrast, Be the weight of said Y direction velocity contrast, n representes the size of said 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 of participating in the image calculating, and said first image and said second image are remote sensing images;

Second acquisition module is used for obtaining best window size according to the changed factor of said plural velocity field from the window of said two above sizes of said first image;

Computing module is used to calculate the corresponding area target motion velocity of window size of said the best.

6. device according to claim 5 is characterized in that, said first acquisition module comprises:

First computing unit; Be used for calculating respectively the related coefficient of the window of two above sizes of said first image and corresponding region of search target on said second image; Wherein, corresponding related coefficient set in the region of search target of the window of each size in said second image;

First acquiring unit is used for obtaining the maximum related coefficient of numerical value in the set of each related coefficient, with the corresponding window of related coefficient of the said maximum on said second image as with said first image on the corresponding match window of window;

Second computing unit is used for obtaining plural velocity field changed factor at the plural match window changed factor of the velocity field of zoning target respectively.

7. device according to claim 5 is characterized in that, said second acquisition module comprises:

The fitting a straight line unit is used at the coordinate system that is formed by said velocity field changed factor and said window size every successively at a distance from the set point number fitting a straight line according to the window size order from small to large of said two above sizes of said first image;

Confirm the unit, be used for when the curvature of confirming said fitting a straight line when setting threshold value, then with the window size of the corresponding window size of the coordinate points of said setting threshold value as the best.

8. according to each described device in the claim 5～7; It is characterized in that; Said computing module is used for calculating according to formula

the corresponding area target motion velocity of window size of said the best; Wherein, Vh is the regional aim plane speed on said first image; Rx, Ry are said first image pixel separation on directions X and the Y direction in the XY coordinate system, and T is the time interval of said first image and said second image taking; Δ x is the coordinate offset amount of window center coordinate on said directions X of said second image and said first image, and Δ y is the coordinate offset amount of window center coordinate on said Y direction of said second image and said first image.

Images