JP6966716B2 - Radar image analysis system - Google Patents

Description

The present invention relates to a radar image analysis system that performs image data analysis processing using image data groups of the same observation point acquired by a synthetic aperture radar (SAR) with different acquisition times.

In recent years, there are various information processing systems that analyze images on the ground using satellite images. Some of these analyze the ground using image data acquired by synthetic aperture radar.

Patent Document 1 describes an example of a radar image analysis system that performs image data analysis processing using image data groups of the same observation points acquired by a synthetic aperture radar and having different acquisition times.

In the image processing apparatus (radar image analysis system) described in Patent Document 1, CCD (Coherent Change Detection) is used. The CCD analyzes the change between the two SAR images by using the phase information together with the luminance information of the SAR image.

Further, Patent Document 1 describes the calculation of coherence (value between complex phases) by pairing a large number of synthetic aperture radar images with each other, the generation of a coherent map from the calculated coherence, and the pixels of interest on the coherent map. It is described that the timing at which the coherence changes (the acquisition time of the synthetic aperture radar image) is derived.

In Patent Document 1, when the coherence is calculated for the image data group of the same observation point where the acquisition time of the above-ground area such as the waterside or the forest is different, the coherence is generated in such a place even if there is substantially no change. The problem is that it is calculated with a low value. Further, in the invention of Patent Document 1, whether or not the coherence value calculated from the pair of synthetic aperture radar images fluctuates from the coherence value calculated from the pair of other synthetic aperture radar images is substantially changed. It is used as an index to judge whether or not.

Japanese Unexamined Patent Publication No. 2009-289111

By the way, there is a PS-InSAR (Persistence Scatterer Interferometric SAR) analysis method for SAR image analysis. This analysis uses a large number of SAR images.

In PS-InSAR analysis, which is a type of SAR image analysis, the change in the direct distance (line-of-sight distance) between each measurement target point of the analysis target (object that reflects radar waves) on the ground surface and the synthetic aperture radar. Is shown on the time axis to derive the trend of secular displacement. Aged displacement of the ground surface is generally expressed in terms of annual average velocity (mm / year).

Each measurement target point needs to reflect radio waves from the SAR to the SAR regardless of whether it is a natural object or an artificial object. In addition, in order to carry out good PS-InSAR analysis, it is required that the point to be measured does not change over time in the reflection characteristics (called backscattering characteristics).

In the existing PS-InSAR analysis, when the reflection characteristics of a certain measurement target point change over time, a method is adopted in which the PS-InSAR analysis is set to not be performed for the measurement target point thereafter. Therefore, when PS-InSAR analysis of a SAR image group is performed over a long period of time, in many cases, the number of measurement target points decreases as the months and days follow from the initial image of the analysis result.

If we examine the points to be improved in the existing PS-InSAR analysis, we can find the problem that the secular displacement cannot be obtained for the measurement target points whose reflection characteristics have changed after the date when it is judged to have changed.

PS-InSAR analysis analyzes a large number of image data to obtain the displacement velocity of the ground surface over a wide area with an accuracy of mm / year for each measurement target point. Therefore, for example, the secular displacement of a specific artificial structure can be detected in millimeters for each part by performing time-series analysis of a large number of image data acquired by using a synthetic aperture radar.

However, if the reflection characteristics of the analysis target location change even temporarily, there is a problem that the analysis of that location after the change in the reflection characteristics cannot be derived.

An object of the present invention has been made in view of the above problems, and to provide a radar image analysis system that performs high-precision image analysis by excluding temporary characteristic change factors from image analysis of a synthetic aperture radar. ..

In the radar image analysis system according to the embodiment of the present invention, an image receiving unit that receives an image data group of a synthetic aperture radar at the same observation point having a different acquisition time and one image included in the image data group are selected as a reference image. At the same time, for the secondary image group used together with the reference image, the amount of change in the backward scattering characteristics was calculated for each slave image for each analysis point of the image analysis, and the amount of change in the backward scattering characteristics exceeded a predetermined value. It is characterized by including an image adjustment analysis unit that removes a slave image for each analysis location and analyzes an individual analysis location using the slave image remaining at each analysis location and the reference image.

In the radar image analysis method according to the embodiment of the present invention, a satellite image data group of a synthetic aperture radar at the same observation point having a different acquisition time is received, one image is selected as a reference image from the satellite image data group, and the above-mentioned A slave image group to be used together with the reference image is selected from the satellite image data group, the amount of change in the backward scattering characteristics is calculated for each slave image for each analysis point of the image analysis, and the backward scattering is performed based on the reference image. It is characterized in that a slave image in which the amount of change in characteristics exceeds a predetermined value is removed for each analysis location, and each analysis location is image-analyzed using the slave image remaining at each analysis location and the reference image. ..

The radar image analysis program according to the embodiment of the present invention includes an information processing system, an image receiving unit that receives an image data group of a synthetic aperture radar at the same observation point having a different acquisition time, and one included in the image data group. An image is selected as a reference image, and for the secondary image group used together with the reference image, the amount of change in the backward scattering characteristics is calculated for each slave image for each analysis point of the image analysis, and the amount of change in the backward scattering characteristics is calculated. It is an operation as an image adjustment analysis unit that removes a slave image exceeding a predetermined value for each analysis location and analyzes each analysis location using the slave image remaining at each analysis location and the reference image. It is a feature.

INDUSTRIAL APPLICABILITY According to the present invention, it is possible to provide a radar image analysis system that performs high-precision image analysis by excluding temporary characteristic change factors from image analysis of a synthetic aperture radar.

It is a block diagram which shows the radar image analysis system 1 of 1st Embodiment of this invention. It is a flowchart which shows the operation example of the image analysis system 1 of 1st Embodiment. It is a block diagram which shows the radar image analysis system 2 of the 2nd Embodiment of this invention. It is explanatory drawing which showed the example of the image data selection method for image analysis. It is a block diagram which showed one configuration example of the image analysis system of this invention.

An embodiment of the present invention will be described with reference to the drawings.

[First Embodiment]
FIG. 1 is a configuration diagram showing a radar image analysis system 1 of the first embodiment.

The radar image analysis system 1 performs image data analysis processing using image data groups of the same observation point acquired by the synthetic aperture radar and having different acquisition times.

The radar image analysis system 1 of the first embodiment includes an image reception unit 10 and an image adjustment analysis unit 20.

The image receiving unit 10 receives an image data group of the synthetic aperture radar. The image data group may be accepted in any way, for example, may be downloaded via a communication network, or the image data group may be read from a recording medium on which the image data group is recorded. Further, the radar image analysis system 1 may access and accept image data required for analysis from a network storage or the like each time. The received image data group is recorded in the work memory, storage, and network storage in the radar image analysis system 1 and prepared so that it can be used for analysis.

The image adjustment analysis unit 20 selects a reference image included in the image data group. The reference image may be, for example, an image arbitrarily specified by the operator, or image data whose image acquisition date is close may be automatically selected by accepting input of the date, and the selection method is not limited.

The image adjustment analysis unit 20 receives the selection of the reference image and selects the secondary image group to be used together with the reference image. The reference image and the subordinate image group used for image analysis are collectively called a selected image. This subordinate image is image data of the same observation point acquired at a time different from the acquisition time of the reference image, which is used for image analysis of the reference image together with other subordinate images. As the subordinate image, for example, all of the image data group received by the image receiving unit 10 (excluding the reference image) may be selected, or an arbitrary number of image data groups may be designated by the operator.

Further, the image adjustment analysis unit 20 automatically selects an image data group within the past 10 years from the reference image, for example, if it wants to know the amount of displacement of the specific analysis portion of the reference image in the past 10 years with respect to the image data as the slave image. It may be extracted. Further, the image adjustment analysis unit 20 automatically extracts a new image data group from the reference image, for example, if it wants to know the specific analysis location of the reference image and the current (latest image) displacement amount of the image data to be the slave image. You just have to do it. As described above, the subordinate image group may be appropriately selected according to the desired image analysis result, and the selection method is not particularly limited.

In addition, the image adjustment analysis unit 20 changes the reflection characteristics of each of the selected subordinate images based on the reflection characteristics (backscattering characteristics) of the reference image for each analysis location extracted from the reference image. Calculate whether the specified value is exceeded. The amount of change in the reflection characteristics may be determined by using the intensity and coherence of the image, and if necessary, both of them. After that, the image adjustment analysis unit 20 marks the analysis points of the individual subordinate images whose change amount exceeds a predetermined value, and does not use the image data of the subordinate image at the time of image analysis of the corresponding analysis points, but analyzes the analysis points. Image analysis is proceeded with the image data of the reference image and the image data (group) of other subordinate images (group) for each. The image data (SAR image data) of the corresponding analysis portion of the removed subordinate image is not removed from the storage or the memory, but is only removed from the image analysis of the corresponding analysis portion. For example, even if a slave image is removed when a certain reference image is selected at a certain analysis point, whether or not the slave image is removed when another reference image is selected is determined. Judgment is made based on the amount of change in the reflection characteristics of the analysis points of the other reference image and the corresponding subordinate image.

The image adjustment analysis unit 20 performs image analysis using the selected subordinate image group (the subordinate image whose extracted change amount exceeds a predetermined value is removed) and the reference image for each analysis location. ..

The image analysis may be configured so that the existing method can be arbitrarily selected, and is not particularly limited. For example, the image adjustment analysis unit 20 may perform PS-InSAR analysis.

As described above, by adjusting the subordinate image group used for image analysis together with the reference image for each analysis location, the selected image is reflected from the image analysis of the synthetic aperture radar based on each analysis location of the reference image, even temporarily. Image analysis can be performed by removing the image of each analysis point whose characteristics have changed. As a result, it becomes easier to obtain desired information by removing data fluctuations that can be said to be noise from the image analysis of the synthetic aperture radar. For example, even if an aircraft or animal accidentally enters between the measurement target point and the synthetic aperture radar when acquiring a certain SAR image, it adversely affects the calculation of the accurate line-of-sight distance. , Before starting the PS-InSAR analysis process, the corresponding image can be automatically removed from the analysis of this part.

In a specific example of image analysis, after adjusting the selected image, the image analysis unit 20 quantifies the reflection characteristics such as coherence between the reference image and each of the remaining secondary images for the specified analysis location, and performs PS-InSAR analysis. The tendency (long-term change leveling trend, trend) along the time series is calculated by. At this time, since the slave image set of each analysis point is adjusted, it is possible to calculate a good line-of-sight direction distance at any analysis point, and the tendency of the actual change with time can be clarified.

The image analysis is not limited to PS-InSAR analysis, and desired image analysis such as generation of a coherent map may be performed.

Next, an operation example of the image analysis system 1 of the first embodiment will be shown.

FIG. 2 is a flowchart showing an operation example of the image analysis system 1 of the first embodiment.

First, the image receiving unit 10 receives an image data group of the synthetic aperture radar as an image preparation (S101). Each image data (SAR image data) is image data with different acquisition times (shooting date / time) acquired for the same observation point periodically / irregularly.

Next, the analysis image adjustment unit 20 selects a reference image and a subordinate image group from the image data group as preparation for analysis, and also selects an analysis portion included in the reference image (S202).

The individual analysis points are generally points in the reference image where the brightness is a predetermined amount or more. Image analysis is performed using the image data of the same location in the subordinate image whose brightness is equivalent to the brightness of the reference image in this analysis location. For example, for each pixel of the reference image, a pixel having substantially the same position as the slave image and having the same brightness remains as an analysis point. If the number of analysis points is large, it may be necessary to perform a process of opening an appropriate interval (appropriate pixel) or accept an input for narrowing the analysis range. In addition, individual analysis points and analysis ranges may be accepted from operators and other computer systems.

Next, the analysis image adjusting unit 20 compares the backscattering characteristics of each analysis portion of the selected reference image with the backscattering characteristics of the portion corresponding to the analysis portion of each subordinate image, and compares the backscattering characteristics of each subordinate image and each analysis portion. The amount of change in backscattering characteristics with respect to the reference image is extracted (S103).

The amount of change in each backscattering characteristic may be obtained in any way. For example, the backscattering characteristics of each analysis point of the reference image selected in advance are acquired, and similarly, the amount of change in the backscattering characteristics of the portion corresponding to the analysis part of the subordinate image is acquired for each subordinate image, and then the comparison is performed. The amount of change in each backscattering characteristic may be acquired. In addition, the backscattering characteristics of each pixel are converted into data for all images in advance, and after the reference image is determined, the backscattering characteristics of the analysis points of the reference image and each secondary image are compared with each other. The amount of change in backscattering characteristics may be acquired.

Next, the analysis image adjusting unit 20 determines whether or not to use the analysis portion for image analysis of the analysis portion for each analysis portion of each selected slave image, and removes the unused slave image from the selected image (S104). More specifically, in the analysis image adjustment unit 20, whether the amount of change in the backscattering characteristic exceeds a predetermined value for the same portion of the reference image of the portion corresponding to the analysis portion of the slave image for each selected slave image. Determine. After that, the analysis image adjusting unit 20 removes the secondary image having the backscattering characteristic change amount exceeding the predetermined value from the image analysis of the analysis portion.

That is, the analysis image adjustment unit 20 removes the secondary image from the selected image for image analysis in which the amount of change in the backward scattering characteristic exceeds a predetermined value with reference to the reference image for each analysis location of the image analysis. The remaining secondary image and reference image are reset to the adjusted selection image and sent to the image analysis of each analysis point.

Next, the analysis image adjustment unit 20 performs image analysis using the selected images adjusted for each analysis location (S105). The image analysis result may be output in a mode desired by the operator as appropriate. For example, the image analysis result may be displayed on the screen as an image or a graph, printed, or recorded in a storage.

By operating the radar image analysis system 1 according to the first embodiment as described above, it is possible to perform high-precision image analysis by excluding temporary change factors of reflection characteristics from the image analysis of the synthetic aperture radar. ..

Next, the second embodiment will be described.

[Second Embodiment]
FIG. 3 is a block diagram showing the radar image analysis system 2 of the second embodiment. The radar image analysis system 2 describes the radar image analysis system 1 in some detail. Further, in the second embodiment, the invention will be described by showing some specific processing operations.

The radar image analysis system 2 of the second embodiment includes an image reception unit 10 and an image adjustment analysis unit 20.

The image receiving unit 10 of the present embodiment receives a satellite image data group of the synthetic aperture radar. The image receiving unit 10 receives satellite image data groups of synthetic aperture radars of the same observation point acquired by an observation satellite using the X band band (for example, COSMO-SkyMed) having different acquisition times. From COSMO-SkyMed, satellite image data of the synthetic aperture radar can be acquired at a resolution of 1 m in a designated 10 km square above the ground. In addition, COSMO-SkyMed can also acquire satellite image data of a synthetic aperture radar at a resolution of 3 m on a 40 km square above the ground.

The image adjustment analysis unit 20 includes an analysis image selection unit 21, an analysis location specifying unit 22, a characteristic change amount acquisition unit 23, and an image analysis unit 24.

The analysis image selection unit 21 accepts the designation of the reference image from the satellite image data group. The reference image is not limited to the selection method as described in the first embodiment.

The analysis point specifying unit 22 extracts a pixel group that can be an analysis point from the satellite image data group as an analysis point group. At this time, the analysis point specifying unit 22 accepts an arbitrary range in the image data of the synthetic aperture radar as an analysis range, and extracts a pixel group that can be an analysis point included in the analysis range as an analysis point group. May be good. Further, the analysis location specifying unit 22 may accept the range of the arbitrary structure in the image data as the analysis range. The analysis range may be determined in any way. For example, you may specify 3 or more points on the image data to specify the range, you may specify the range with an appropriate frame, or you can acquire the range and coordinates of the structure from the existing database. , The analysis range may be defined around the range and coordinates. If the analysis range is not specified, the radar image analysis system 2 may operate with the entire satellite image as the analysis range.

The characteristic change amount acquisition unit 23 extracts the change amount of the backscattering characteristic with respect to the reference image of each analysis point included in the analysis range for each slave image group. Extraction of this amount of change may refer to changes in image intensity distribution, coherence, or both, as in the first embodiment.

The image analysis unit 24 removes the secondary image whose amount of change exceeds a predetermined value for each analysis location, and executes image analysis of the analysis range using the remaining secondary image and the reference image.

In this image analysis, for example, the image analysis unit 24 quantifies the coherence with each of the subordinate images with respect to the reference image for each of the designated analysis points, and the tendency (long-term) of each analysis point along the time series. Change leveling trend, trend) should be calculated.

Further, the image analysis unit 24 may extract the displacement difference of the tendency of each analysis point by comparing the tendency of each analysis point along the time series derived by the image analysis. If the displacement difference of the extracted tendency is visualized and depicted on an image map, for example, the variation of the displacement difference between the position coordinates and the analysis point can be extracted.

Further, by performing the same image analysis while moving the reference image along the time series and arranging them side by side, it is possible to obtain an image analysis result in which it is easy to identify the position coordinates having an abnormal tendency. For example, if the image analysis results are superimposed and viewed in order to the future or the past in chronological order, it is easy to obtain analysis data that facilitates mining of new information. For example, it may be possible to discover that the sinking speed of only a part of the level ground is faster than that of the surroundings, and that the tilt angle of a huge rod is different. This mined information may be useful for setting, for example, an important point of interest for close visual confirmation.

For image analysis, PS-InSAR analysis is performed using a reference image and one or more secondary images as selected images, and if pixels of a 1 m resolution SAR image are used as the range of individual analysis points, fluctuations of several millimeters per year can be measured. It is possible. A long-term fluctuation tendency is extracted from this highly accurate fluctuation, and the part where some data acquisition problem temporarily occurs is excluded from the analysis based on the presence or absence of a change in the backscattering characteristics of each analysis location (mutual pixels). do. This makes it possible to derive a more reliable displacement tendency of the analysis range (analysis target).

Further, by limiting the image analysis range to only the range of the arbitrary structure in advance, the image analysis result excluding the unnecessary range can be obtained. According to this mechanism, it is possible to obtain image analysis results that are easy to use for various purposes while reducing computational resources. For example, the image analysis unit 24 outputs an image analysis result that is an analysis result of only an arbitrary structure and is not analyzed at the same time for neighboring structures.

As described above, the radar image analysis system 2 according to the second embodiment can enhance the image analysis of the synthetic aperture radar.

Here, an adjustment method of the selected image of the image adjustment analysis unit 20 will be described with reference to an example.

FIG. 4 is an explanatory diagram showing an example of an image data selection method for image analysis. FIG. 4A shows an example SAR image set and its features. FIG. 4B is an example of a processing flow of the adjustment method of the selected image. In this example, two pixels (two points) in the SAR image indicated by Ax and Bx are set as analysis points.

Images 1 to Image n shown in FIG. 4A are a large number of SAR images of the same location acquired at different times. Ax and Bx in each Image indicate different analysis ranges (pixels) on the image. In addition, all SAR images other than the reference image are used as slave images. After deciding pixel A and pixel B as analysis points, the backscattering characteristics of each secondary image of the analysis point (pixel A or pixel B) are compared with the backscattering characteristics of the reference image to see if they have changed. It is determined for each subordinate image.

In this case, this comparison is made using the intensity distribution and coherence of the analysis site.

If the intensity distribution (or coherence) of the analysis point A changes from the intensity distribution (or coherence) of the reference image (Image m corresponds to that in FIG. 4A), the image adjustment analysis unit 20 is shown in FIG. As shown in (b), when analyzing the analysis point A, image analysis is performed using a data set of other images for analysis excluding Image m.

Further, if the intensity distribution (or coherence) of the analysis point B changes from the intensity distribution (or coherence) of the reference image (Image k corresponds to that in FIG. 4A), the image adjustment analysis unit 20 may perform the image adjustment analysis unit 20. As shown in FIG. 4 (b), when analyzing the analysis point B, image analysis is performed using a data set of images for analysis other than Image k.

When time-series analysis such as PS-InSAR analysis is performed as image analysis using this selected image adjustment method, even if the backscattering characteristics of the analysis location change temporarily, the accuracy of that location is high. You will be able to perform various analyzes. It should be noted that the temporary change in backscattering characteristics occurs, for example, depending on whether or not a car is stopped, the number of cars, and the type of the car on the ground surface, or when an animal, a train, or a flying object passes through (reflection). Since the SAR image has a higher resolution and can be analyzed in millimeters, for example, even if a person is not at the analysis site, the existing image analysis result may be significantly changed. According to the adjustment method of the selected image, it is possible to eliminate such a temporary characteristic change factor and perform high-precision image analysis.

Each part of the radar image analysis system may be realized by using a combination of hardware and software of a computer system as illustrated in FIG. The computer system includes one or more processors and memory tailored to the desired form. Further, in the form of this computer system, in each part, a program for radar image analysis is developed in the above memory, and based on this program, hardware such as one or more processors is operated by an execution instruction group or a code group. , Should be realized. At this time, if necessary, this program may realize each part in cooperation with the functions provided by the software such as the operating system, the microprogram, and the driver.

The program data expanded in the memory appropriately includes an execution instruction group, a code group, a table file, content data, and the like that operate the processor as one or more of the above-mentioned parts.

In addition, this computer system does not necessarily have to be built as a single device, but is built by combining multiple servers / computers / virtual machines, so-called thin clients, distributed computing, and cloud computing. May be good. Further, a part / all parts of the computer system may be replaced with hardware or firmware (for example, one or more LSIs: Large-Scale Integration, FPGA: Field Programmable Gate Array, a combination of electronic elements). Similarly, only a part of each part may be replaced with hardware or firmware.

In addition, this program may be temporarily recorded on a recording medium and distributed. The program recorded on the recording medium is read into a memory via a wired, wireless, or recording medium itself, and operates a processor or the like.

In the present specification, the recording medium includes a storage medium, a memory device, a storage device, and the like of similar terms. Examples of this recording medium include optical disks, magnetic disks, semiconductor memory devices, hard disk devices, tape media, and the like. Further, it is desirable that the recording medium is non-volatile. Further, as the recording medium, a combination of a volatile module (for example, RAM: Random Access Memory) and a non-volatile module (for example, ROM: Read Only Memory) may be used.

The present invention has been described by exemplifying embodiments. However, the specific configuration of the present invention is not limited to the above-described embodiment, and is included in the present invention even if there are changes to the extent that the gist of the present invention is not deviated. For example, changes such as separation and merging of the block configuration of the above-described embodiment and replacement of procedures are free as long as the gist of the present invention and the functions described are satisfied, and the above description does not limit the present invention.

In addition, some or all of the above embodiments may also be described as follows. The following notes do not limit the present invention in any way.
[Appendix 1]
An image reception unit that accepts image data groups of synthetic aperture radar at the same observation point with different acquisition times,
A reference image included in the image data group is selected, and for the slave image group used together with the reference image, a slave image in which the amount of change in the backward scattering characteristics exceeds a predetermined value for each analysis point of the image analysis is selected. An image adjustment analysis unit that performs image analysis of individual analysis points using the remaining secondary image and the reference image after removal.
A radar image analysis system characterized by including.

[Appendix 2]
The radar image analysis according to Appendix 1, wherein the image adjustment analysis unit extracts the amount of change in the backscattering characteristic by using the change in the image intensity distribution, coherence, or both as the change in the backscattering characteristic. system.

[Appendix 3]
As an image analysis, the image analysis unit quantifies the backward scattering characteristics of the designated analysis points with respect to each of the slave images with respect to the reference image, and calculates the tendency of each analysis point along the time series. The radar image analysis system according to Appendix 1 or 2, wherein the radar image analysis system is characterized by the above.

[Appendix 4]
As an image analysis, the image analysis unit quantifies the backward scattering characteristics of the designated analysis points with each of the slave images with respect to the reference image, and calculates the tendency of each analysis point along the time series. The radar image analysis system according to any one of Supplementary note 1 to 3, wherein the displacement difference of the tendency of each analysis point is extracted.

[Appendix 5]
The image adjustment analysis unit
PS-InSAR (Persistence Scatterer Interferometric SAR (Synthetic Aperture Radar)) analysis can be performed using the reference image and multiple subordinate images from 1 as the selected image, and pixels are used as the range of analysis points.
PS- The radar image analysis system according to any one of Supplementary note 1 to 4, which comprises performing InSAR analysis.

[Appendix 6]
The image adjustment analysis unit
An analysis image selection unit that accepts the designation of a reference image from within the image data group,
An analysis location identification unit that extracts a pixel group that can be an analysis location from the image data group as an analysis location group, and an analysis location identification unit.
For the slave image group, a characteristic change amount acquisition unit that extracts the change amount of the backscattering characteristic with respect to the reference image of each analysis point for each slave image, and
For each analysis location, an image analysis unit that removes the secondary image whose amount of change exceeds a predetermined value and analyzes the image using the remaining secondary image and the reference image.
The radar image analysis system according to any one of Supplementary note 1 to 5, wherein the radar image analysis system comprises.

[Appendix 7]
The image receiving unit receives satellite image data groups of synthetic aperture radars at the same observation points with different acquisition times acquired by observation satellites using the X-band band.
The analysis point specifying unit accepts an arbitrary range in the image data of the synthetic aperture radar as an analysis range, and extracts a pixel group that can be an analysis point included in the analysis range as an analysis point group.
The characteristic change amount acquisition unit extracts the amount of change in the backscattering characteristic of each of the analysis points included in the analysis range for each slave image with respect to the reference image for the slave image group.
The image analysis unit removes the secondary image whose amount of change exceeds a predetermined value for each analysis location, and executes image analysis of only an arbitrary range using the remaining secondary image and the reference image. The radar image analysis system according to Appendix 6, wherein the radar image analysis system is characterized by the above.

[Appendix 8]
The analysis location identification unit accepts the range of an arbitrary structure in the image data of the synthetic aperture radar as the analysis range.
The image analysis unit removes the secondary image whose amount of change exceeds a predetermined value for each analysis location, and uses the remaining secondary image and the reference image to image only the range of the arbitrary structure. The radar image analysis system according to Appendix 7, wherein the analysis is performed.

[Appendix 9]
Accepts satellite image data groups of synthetic aperture radars at the same observation point with different acquisition times.
Select a reference image from the satellite image data group and select it.
A subordinate image group to be used together with the reference image is selected from the satellite image data group.
For each analysis point of the image analysis, the subordinate image in which the amount of change in the backward scattering characteristic exceeds a predetermined value is removed based on the reference image, and the remaining subordinate image and the reference image are used for individual analysis. A radar image analysis method characterized by image analysis of a location.

[Appendix 10]
Information processing system,
An image reception unit that accepts image data groups of synthetic aperture radar at the same observation point with different acquisition times,
A reference image included in the image data group is selected, and for the slave image group used together with the reference image, a slave image in which the amount of change in the backward scattering characteristics exceeds a predetermined value for each analysis point of the image analysis is selected. An image adjustment analysis unit that removes and analyzes individual analysis points using the remaining secondary image and the reference image.
A program for radar image analysis characterized by operating as.

1, 2, Radar image analysis system 10 Image reception unit 20 Image adjustment analysis unit 21 Analysis image selection unit 22 Analysis location identification unit 23 Characteristic change amount acquisition unit 24 Image analysis unit

Claims (10)

An image reception unit that accepts image data groups of synthetic aperture radar at the same observation point with different acquisition times,
With selecting an image included in the image data group as a reference image, the subsidiary image set to be used together with the reference image, for each analysis point of each image analysis to calculate a variation of the backscatter characteristics for each sub-image , Image adjustment that removes the secondary image in which the amount of change in the backward scattering characteristics exceeds a predetermined value for each analysis location, and analyzes each analysis location using the secondary image remaining at each analysis location and the reference image. Analysis department and
A radar image analysis system characterized by including. The radar image according to claim 1, wherein the image adjustment analysis unit extracts the amount of change in the backscattering characteristic by using the change in the image intensity distribution, coherence, or both as the change in the backscattering characteristic. Analysis system. As an image analysis, the image adjustment analysis unit quantifies the backward scattering characteristics of each of the slave images with respect to the reference image for the designated analysis points, and the displacement tendency of each analysis point along the time series. The radar image analysis system according to claim 1 or 2, wherein the image is calculated. As an image analysis, the image adjustment analysis unit quantifies the backward scattering characteristics of each of the slave images with respect to the reference image for the designated analysis points, and the displacement tendency of each analysis point along the time series. The radar image analysis system according to any one of claims 1 to 3, wherein the image is calculated and the displacement difference of the displacement tendency of each analysis point is extracted. The image adjustment analysis unit
A reference image and multiple sub-image selection image, employs a pixel as a range of solutions析箇plants,
Based on the presence or absence of changes in the backward scattering characteristics between the reference image included in the selected image and the pixels used as the analysis points of the subordinate image, if there is a change, the pixels at the analysis point of the subordinate image are excluded from PS-. The radar image analysis system according to any one of claims 1 to 4, wherein InSAR (Persistence Scatterer Interferometric SAR (Synthetic Aperture Radar)) analysis is performed. The image adjustment analysis unit
And analyzing the image selection unit that receives designation of the reference image from the image data within the group,
An analysis location identification unit that extracts a pixel group that can be an analysis location from the image data group as an analysis location group, and an analysis location identification unit.
For the slave image group, a characteristic change amount acquisition unit that extracts the change amount of the backscattering characteristic with respect to the reference image of each analysis point for each slave image, and
An image analysis unit that removes the secondary image whose amount of change exceeds a predetermined value for each analysis location and analyzes the image using the secondary image remaining at each analysis location and the reference image. When,
The radar image analysis system according to any one of claims 1 to 5, wherein the radar image analysis system comprises. The image receiving unit receives satellite image data groups of synthetic aperture radars at the same observation points with different acquisition times acquired by observation satellites using the X-band band.
The analysis point specifying unit accepts an arbitrary range in the image data of the synthetic aperture radar as an analysis range, and extracts a pixel group that can be an analysis point included in the analysis range as an analysis point group.
The characteristic change amount acquisition unit extracts the amount of change in the backscattering characteristic of each of the analysis points included in the analysis range for each slave image with respect to the reference image for the slave image group.
The image analysis unit removes the slave image whose amount of change exceeds a predetermined value for each analysis location, and uses the slave image remaining at each analysis location and the reference image. The radar image analysis system according to claim 6, wherein the image analysis of only an arbitrary range is performed. The analysis location identification unit accepts the range of an arbitrary structure in the image data of the synthetic aperture radar as the analysis range.
The image analysis unit removes the slave image whose amount of change exceeds a predetermined value for each analysis location, and uses the slave image remaining at each analysis location and the reference image. The radar image analysis system according to claim 7, wherein the image analysis is performed only in the range of an arbitrary structure. Accepts satellite image data groups of synthetic aperture radars at the same observation point with different acquisition times.
One image is selected as a reference image from the satellite image data group, and the image is selected.
A subordinate image group to be used together with the reference image is selected from the satellite image data group.
The amount of change in the backscattering characteristics is calculated for each subordinate image for each analysis location in the image analysis, and the subordinate image in which the amount of change in the backscattering characteristics exceeds a predetermined value based on the reference image is analyzed for each subordinate image. A radar image analysis method characterized by performing image analysis of individual analysis points using a secondary image remaining at each analysis point and the reference image after removal. Information processing system,
An image reception unit that accepts image data groups of synthetic aperture radar at the same observation point with different acquisition times,
With selecting an image included in the image data group as a reference image, the subsidiary image set to be used together with the reference image, for each analysis point of each image analysis to calculate a variation of the backscatter characteristics for each sub-image , Image adjustment that removes the secondary image in which the amount of change in the backward scattering characteristics exceeds a predetermined value for each analysis location, and analyzes each analysis location using the secondary image remaining at each analysis location and the reference image. Analysis department,
A program for radar image analysis characterized by operating as.

Images