WO2019106850A1 - Ｓａｒ画像解析システム、画像処理装置、画像処理方法および画像処理プログラム - Google Patents

Ｓａｒ画像解析システム、画像処理装置、画像処理方法および画像処理プログラム Download PDF

Info

Publication number
WO2019106850A1
WO2019106850A1 PCT/JP2017/043377 JP2017043377W WO2019106850A1 WO 2019106850 A1 WO2019106850 A1 WO 2019106850A1 JP 2017043377 W JP2017043377 W JP 2017043377W WO 2019106850 A1 WO2019106850 A1 WO 2019106850A1
Authority
WO
WIPO (PCT)
Prior art keywords
phase
correlation
image processing
reflection point
pixel
Prior art date
Application number
PCT/JP2017/043377
Other languages
English (en)
French (fr)
Inventor
大地 田中
宝珠山　治
Original Assignee
日本電気株式会社
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 日本電気株式会社 filed Critical 日本電気株式会社
Priority to US16/767,805 priority Critical patent/US11709254B2/en
Priority to JP2019556532A priority patent/JP7020494B2/ja
Priority to PCT/JP2017/043377 priority patent/WO2019106850A1/ja
Publication of WO2019106850A1 publication Critical patent/WO2019106850A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021—SAR image post-processing techniques
    • G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis

Definitions

  • the present invention relates to a SAR image analysis system, an image processing apparatus, an image processing method, and an image processing program.
  • Patent Document 1 describes an example of a SAR (Synthetic Aperture Rader) image analysis system.
  • SAR Synthetic Aperture Rader
  • a model fitting is performed on the phase observed at each stable reflection point candidate, and a model in which displacement occurs linearly with respect to time is assumed as a fitting model.
  • An object of the present invention is to provide a technique for solving the above-mentioned problems.
  • an image processing apparatus is Phase array storage means for storing a phase array of each pixel across a plurality of images; Clustering means for clustering the respective pixels into a plurality of clusters based on the phase array; A phase correlation calculating unit that calculates a correlation between a phase change in each of the plurality of clusters and a phase array in each of the pixels; Equipped with
  • an image processing method is A phased array storage step of storing the phased array of each pixel across a plurality of images in a phased array storage means; Clustering the respective pixels into a plurality of clusters based on the phase array; Calculating a correlation between a phase change in each of the plurality of clusters and a phase array in each of the pixels; including.
  • an image processing program is A phased array storage step of storing the phased array of each pixel across a plurality of images in a phased array storage means; Clustering the respective pixels into a plurality of clusters based on the phase array; Calculating a correlation between a phase change in each of the plurality of clusters and a phase array in each of the pixels; On a computer.
  • the SAR image analysis system is An image processing apparatus for identifying stable reflection points based on a plurality of SAR images; A SAR image analyzer for analyzing a SAR image using data of the stable reflection point; SAR image analysis system comprising The image processing apparatus at least Phase array storage means for storing a phase array of each pixel across the plurality of SAR images; Clustering means for clustering the respective pixels into a plurality of clusters based on the phase array; A phase correlation calculating unit that calculates a correlation between a phase change in each of the plurality of clusters and a phase array in each of the pixels; Equipped with
  • the SAR image analysis method is An image processing step for identifying a stable reflection point based on a plurality of SAR images; A SAR image analysis step of analyzing a SAR image using data of the stable reflection point; SAR image analysis method including
  • the image processing step includes at least A phase array storage step of storing the phase array of each pixel across the plurality of SAR images in a phase array storage unit; Clustering the respective pixels into a plurality of clusters based on the phase array; Calculating a correlation between a phase change in each of the plurality of clusters and a phase array in each of the pixels; including.
  • the present invention it is possible to generate an evaluation index of likelihood of a stable reflection point without being influenced by a large phase change for an object whose displacement is nonlinear, an object whose elevation is high, and an object whose displacement is large.
  • the image processing apparatus 100 is an apparatus that generates an evaluation index of likelihood of a stable reflection point from the phase array of each pixel across a plurality of images.
  • the image processing apparatus 100 includes a phase array storage unit 101, a clustering unit 102, and a phase correlation calculation unit 103.
  • the phased array storage unit 101 stores the phased array of each pixel across a plurality of images.
  • the clustering unit 102 clusters each pixel into a plurality of clusters based on the phase array.
  • the phase correlation calculation unit 103 calculates the correlation between the phase change in each of the plurality of clusters and the phase arrangement of each pixel.
  • the correlation between the phase change in each cluster of a plurality of clusters and the phase arrangement at each pixel is calculated, and therefore, for objects whose displacement is nonlinear, objects whose elevation is high, and objects whose displacement is large
  • An evaluation index of stable reflection point likeness can be generated without being influenced by a large phase change.
  • FIG. 17 is a diagram showing a functional configuration of the SAR image analysis system 1700 of the base technology described in Patent Document 1. As shown in FIG.
  • the SAR image analysis system 1700 includes a SAR image storage unit 1701, a phase identification unit 1702, a stable reflection point candidate identification unit 1703, a model fitting unit 1704, and a model error identification unit 1705; A stable reflection point identification unit 1706 and a displacement / altitude identification unit 1707 are provided.
  • the SAR image analysis system 1700 having such a configuration operates as follows. First, the phase specifying unit 1702 inputs, from the SAR image storage unit 1701, a plurality of SAR images obtained by capturing the same place at different times, and specifies changes in the phase of each pixel.
  • the stable reflection point candidate identification unit 1703 receives a plurality of SAR images obtained by photographing the same place at different times, and identifies stable reflection point candidates that are pixels having a high possibility of stable reflection.
  • stable reflection point candidates for example, in Patent Document 1, stable reflection is achieved in that a plurality of SAR images captured at different times in the same place are input and the change in reflection intensity is not large across the plurality of SAR images. It is identified as a point candidate.
  • the model fitting unit 1704 receives the stable reflection point candidate specified by the stable reflection point candidate specifying unit as an input, and performs model fitting on the phase observed in each stable reflection point candidate.
  • the phase on the SAR image is an amount that changes according to the altitude, displacement, the state of the atmosphere, and the like. That is, the fitted model represents how the altitude / displacement of the object to be imaged, the state of the atmosphere, and the like change depending on the time when the SAR image is taken.
  • Patent Document 1 assumes a model in which displacement occurs linearly with respect to time.
  • the model error identification unit 1705 identifies how much the fitting performed by the model fitting unit deviates.
  • the stable reflection point identification unit 1706 receives the error identified by the model error identification unit as an input, and identifies a point with a small error as the stable reflection point.
  • the displacement / altitude identification unit 1707 identifies displacement and elevation on these stable reflection points as input. As described above, stable reflection points less affected by noise, and precise displacement and elevation at each stable reflection point are output.
  • Such an SAR image analysis system 1700 has a problem that displacement measurement can not be successfully performed on an object whose displacement is nonlinear, an object whose elevation is high, and an object whose displacement is large.
  • the first cause is that it is determined whether or not a non-linear displacement object is not a stable reflection point because it is evaluated whether it is stable based on the magnitude of the error from the linear displacement model.
  • the problem is that Another cause is that the phase fitting itself is computationally difficult.
  • the phase is an amount like a position on the circumference, and for example, even if it causes a phase change of one turn or a phase change of two turns, it returns to the origin when observed after the change, It is indistinguishable from the case where there is no phase change at all. Therefore, when the phase change is large, such as when there is a large displacement or when the elevation is high, the correct phase change is not known, and the model can not be fitted well.
  • Each phase of a pixel which is at the same point as each pixel of one SAR image to be is included as each element of the array.
  • the degree of stable reflection point candidates to be evaluated is divided into point groups with similar phase changes, and it is evaluated whether the average phase change and the phase change of the point to be evaluated are correlated or not. It can be determined whether it is not affected by noise. Since the above process does not include the process of fitting the model, stable reflection points can be specified even for an object for which it is difficult to apply a phase model, such as a high-rise building.
  • FIG. 2 is a block diagram showing a functional configuration of the image processing apparatus 200 according to the present embodiment.
  • the image processing apparatus 200 includes, as the minimum configuration, a phase array storage unit 201, a clustering unit 202, and a phase correlation calculation unit 203.
  • the image processing apparatus 200 may also include a SAR image storage unit 204, a stable reflection point candidate identification unit 205, and a stable reflection point identification unit 206.
  • the phase array storage unit 201 stores a phase array which is an array in which the phases of the pixels in the entire image are recorded.
  • the clustering unit 202 divides pixels into clusters having similar phases using a clustering method. As a clustering method, for example, although the method described in Non-Patent Document 1 may be used, it is not limited thereto.
  • the phase correlation calculation unit 203 calculates the correlation between the representative phase change in each cluster and the phase change in each of the pixels based on the clustering result and the phase arrangement of each of the pixels in the entire image.
  • the stable reflection point candidate identification unit 205 identifies a pixel that can be a stable reflection point.
  • the clustering unit 202 clusters only stable reflection point candidates.
  • the clustering unit 202 clusters pixels of the entire image.
  • the stable reflection point identification unit 206 stabilizes pixels having a correlation exceeding a predetermined threshold based on the correlation between each pixel calculated in the phase correlation calculation unit 203 and each cluster.
  • a stable reflection point is specified on condition that it is a reflection point.
  • the phased array storage unit 201 stores the phased array in each of the pixels of the entire image. In storing the phase array, the phase array storage unit 201 generates, for each pixel, an array having the phase of each pixel of the plurality of SAR images as an element. This phase may be determined, for example, as a phase difference between a reference SAR image and another SAR image. In this case, among the plurality of input SAR images, a reference SAR image is determined in advance. Then, the phase array storage unit 201 obtains the phase difference between the reference SAR image and each of the other SAR images, and stores them as a phase array in each pixel.
  • the phased array storage unit 201 may store phases of a plurality of SAR images in time-series order or the like as the phased array without defining the reference SAR image. Each phase is stored, for example, in the form of a complex number normalized so that the absolute value is "1".
  • the phased array storage unit 201 may remove a change in phase that may change depending on the elevation of the object using a DEM (Digital Elevation Model) or the like. It may be removed using, for example, a known method.
  • DEM Digital Elevation Model
  • the phase correlation calculation unit 203 calculates the correlation between the representative phase change of each cluster and the phase arrangement of each pixel of the entire image.
  • the calculated correlation can be used as an evaluation index for determining that the stronger the reflection point seems to be.
  • FIG. 3A is a diagram showing an outline of a SAR image analysis system 300 including the image processing apparatus 200 according to the present embodiment.
  • the SAR satellite 350 captures a SAR image.
  • the SAR satellite 350 captures an SAR image of a surface or a structure located on the ground, such as a high-rise building or a high-rise apartment.
  • the SAR image captured by the SAR satellite 350 is stored, for example, in the SAR image storage unit 204.
  • the SAR satellite 350 may transmit the captured SAR image to the SAR image analysis device 310.
  • the SAR image analysis device 310 acquires a plurality of SAR images captured by the SAR satellite 350.
  • a plurality of SAR images acquired by the SAR image analysis apparatus 310 are, for example, captured at the same spot, the same region, the same range, etc. of the ground surface under different imaging conditions (date, weather, SAR satellite used for imaging, etc.) It is a SAR image.
  • the SAR image analysis device 310 may acquire the SAR image stored in the SAR image storage unit 204.
  • the SAR image analysis device 310 may acquire an SAR image directly from the SAR satellite 350.
  • the SAR image analysis device 310 specifies a stable reflection point from the acquired SAR image.
  • the SAR image analyzer 310 further identifies the phase at each of the identified stable reflection points.
  • FIG. 3B is a diagram showing a functional configuration of the SAR image analysis system 300 including the image processing apparatus 200 according to the present embodiment.
  • the SAR image analysis system 300 includes a SAR image storage unit 204, an SAR image analysis device 310, a stable reflection point storage unit 330, and a displacement / altitude storage unit 340.
  • the SAR image analysis device 310 includes the image processing device 200 of this embodiment, the stable reflection point identification unit 206, and the displacement / altitude calculation unit 311.
  • the SAR image analysis device 310 operates as follows.
  • the image processing apparatus 200 stores SAR image of a plurality of SAR images captured at different times in the same place in order to analyze the difference in phase observed at the same point (pixel) across a plurality of SAR images. Acquired from the unit 204.
  • the image processing apparatus 200 according to the present embodiment generates a correlation between a representative phase change of a clustered cluster and a phase array of each pixel as an evaluation index of the likelihood of a stable reflection point.
  • the stable reflection point identification unit 206 identifies a stable reflection point based on the acquired evaluation index.
  • the identified stable reflection point is stored, for example, in the stable reflection point storage unit 330.
  • the displacement / altitude calculation unit 311 inputs the stable reflection point identified by the stable reflection point identification unit 206, and calculates the displacement and elevation at each stable reflection point.
  • the calculated displacement and altitude are stored, for example, in the displacement / altitude storage unit 340.
  • FIG. 4 is a block diagram showing a functional configuration of the phase correlation calculation unit 203 according to the present embodiment.
  • the phase statistics computing unit 431 identifies, for each cluster, a statistic representing a representative phase change of each cluster. For example, the average of the phase arrangement in each of the pixels constituting each cluster is used as the statistic indicating the representative phase change of each cluster.
  • the weighted average of the phase alignment may be specified as a representative phase change of each cluster using a weight function that changes according to the distance.
  • w is a monotonically non-increasing positive function (d) is, for example with respect to the positive input d, the position vector of the respective pixels r i, the phase Assuming that a complex column vector storing each element of the array is s i , a vector s C (r) representing a typical phase change of the cluster C at the location r is specified as the following equation (1) .
  • a statistic representing a typical phase change of each cluster a matrix having, as elements, phase differences for all pairs of elements of the phase array, is calculated for the phase array in each pixel constituting the cluster, The average value of these matrices may be calculated and used.
  • the matrix which records the phase difference for all pairs of elements of the phased array if tentatively determined for the pixel i, s i with a column vector s i of matches an element of the phase sequence of the pixel i s i * It can be determined as However, * is a symbol showing complex conjugate transposition.
  • * is a symbol showing complex conjugate transposition.
  • the weighted average of the matrix s i s i * may be a statistic representing a representative phase change of each cluster. That is, using the same letter definition as used in defining s C (r), the statistic representing a typical phase change of a cluster is in the form of a matrix weighted average S c (r) , Is specified as in the following equation (2).
  • other statistics such as a mode value or median value may be taken instead of the average value of vectors or matrices.
  • the statistics of complex numbers are directly calculated on a vector or matrix of complex numbers in the above, it is also possible to once take out the argument, specify the statistic for the argument, and return it to the complex number. Further, in the average calculation, in general, the absolute value of each element of the vector or matrix to be obtained does not become “1”, but it may be renormalized so that the absolute value becomes “1”.
  • the phase correlation calculation unit 432 calculates the correlation between the representative phase change of the cluster and the phase change of each pixel of the entire image. It can be evaluated that the pixel with the stronger correlation has less noise and is more likely to be a stable reflection point. In calculating the correlation, for example, the following calculation is performed on each pixel of the entire image. Temporarily, let i be the index of the pixel specifying the correlation, s i be the column vector containing the elements of the phase array in the pixel i, and r i be the position vector of the pixel i. As a statistic representing a typical phase change of the cluster C, it is assumed that the vector s C (r) defined in the equation (1) is used.
  • the phase correlation operation unit 432 calculates the correlation corr (i, C) of the representative phase change of the cluster C and the phase change at the pixel i by the following equation (3) for all pixels i of the whole image. Identify for each.
  • the matrix s C (r) defined in the equation (2) is used as a statistic representing a typical phase change of the cluster C
  • the phase correlation operation unit 432 calculates the typical phase change of the cluster C.
  • the correlation corr (i, C) of the phase change in the pixel i is specified for each of all the pixels i by the following equation (4).
  • v_one is a column vector such that all elements are 1.
  • a value obtained by appropriately applying conversion by a monotonous non-decreasing function to a positive input such as a square root or a square may be used to a value specified as Expression (4).
  • the correlation between the representative phase change of each of the clusters described above and the phase change of each of the pixels in the entire image is obtained by the number of clusters for each of the pixels in the entire image.
  • the phase correlation operation unit 432 determines the correlation between each pixel and a plurality of clusters. Are integrated into one, and the value is output as a correlation in each pixel.
  • the grouping process is realized, for example, by specifying a cluster including each of the pixels of the entire image and using only the correlation with the representative phase change of the cluster.
  • the largest correlation may be extracted and used from among the plurality of correlations identified with each cluster.
  • a cluster closest to each of the pixels of the entire image may be specified, and only correlation with a representative phase change of the cluster may be used.
  • FIG. 5A is a diagram showing the configuration of the SAR image storage unit 204 according to the present embodiment.
  • the SAR image storage unit 204 accumulates SAR images captured and transmitted by a satellite or the like.
  • the configuration of the SAR image storage unit 204 is not limited to that shown in FIG. 5A.
  • the SAR image storage unit 204 stores a reflection point ID (image ID) 512 included in the SAR image of each SAR image ID 511.
  • a luminance 513, a position 514, a phase 515 and the like are stored in association with each reflection point ID (image ID) 512.
  • the phased array storage unit 201 stores a phased array 522 composed of a set of a time stamp and a phase of each SAR image corresponding to each reflection point ID (pixel ID) 521.
  • FIG. 5C is a diagram showing the configuration of the stable reflection point candidate table 530 according to the present embodiment.
  • the stable reflection point candidate table 530 is used by the stable reflection point candidate identification unit 205 to identify stable reflection point candidates among reflection points (pixels) in a plurality of SAR images.
  • the configuration of the stable reflection point candidate table 530 is not limited to FIG. 5C.
  • the stable reflection point candidate table 530 is a stable reflection point candidate condition 532 for specifying a stable reflection point candidate in association with each reflection point ID (pixel ID) 531, and the stability of the determination result by the stable reflection point candidate condition 532
  • the reflection point candidate flag 533 is stored.
  • the stable reflection point candidate flag 533 is assumed to be a stable reflection point candidate in the case of “1” and not a stable reflection point candidate in the case of “0”.
  • the stable reflection point candidate condition 532 includes the stability of luminance, the stability of position, the stability of phase, and the like, but the determination is made based on at least one condition not limited thereto.
  • the cluster table 540 stores, in correspondence with the cluster ID 541, the reflection point ID (pixel ID) 542 classified into the cluster, the phase array 543 thereof, and the stable reflection point candidate flag 544.
  • the reflection point (pixels) of all “1” s of the stable reflection point candidate flags 544 are clustered.
  • FIG. 5E is a diagram showing the configuration of the phase statistic table 550 according to the present embodiment.
  • the phase statistic table 550 is used by the phase statistic calculator 431 to calculate a representative phase change of each cluster.
  • the configuration of the phase statistic table 550 is not limited to FIG. 5E.
  • FIG. 5F is a diagram showing the configuration of the phase correlation table 560 according to the present embodiment.
  • the phase correlation table 560 is used by the phase correlation operation unit 432 to calculate the correlation between the phase array of each reflection point (pixel) and the phase statistic (phase change) of each cluster.
  • the configuration of the phase correlation table 560 is not limited to that shown in FIG. 5F.
  • the phase correlation table 560 calculates and stores the phase correlation values of all combinations between the phase array group 561 of each reflection point (pixel) and the phase statistic (phase change) group 562 of each cluster.
  • FIG. 5G is a view showing the configuration of the stable reflection point table 570 according to the present embodiment.
  • the stable reflection point table 570 is used by the stable reflection point identification unit 206 to identify a stable reflection point based on the phase correlation value from the phase correlation calculation unit 203.
  • the configuration of the stable reflection point table 570 is not limited to FIG. 5G.
  • the stable reflection point table 570 stores a stable reflection point candidate flag 572, a phase correlation calculation result 573, and a stable reflection point flag 574 in association with each reflection point ID (pixel ID) 571.
  • the stable reflection point flag 574 is a stable reflection point in the case of “1” and not a stable reflection point in the case of “0”.
  • FIG. 6 is a block diagram showing the hardware configuration of the image processing apparatus 200 according to the present embodiment. Although FIG. 6 shows a configuration in which the stable reflection point candidate identification unit 205 is provided and the stable reflection point identification unit 206 is not provided, the present invention is not limited to this.
  • a CPU (Central Processing Unit) 610 is a processor for arithmetic control, and implements the functional configuration units in FIG. 2 and FIG. 4 by executing a program. Note that there may be a plurality of CPUs 610 corresponding to the respective functions.
  • a ROM (Read Only Memory) 620 stores fixed data and programs such as initial data and programs.
  • the network interface 630 controls communication with other devices such as the SAR image storage unit 204 via the network.
  • a random access memory (RAM) 640 is a random access memory used by the CPU 610 as a work area for temporary storage. In the RAM 640, an area for storing data necessary for realizing the present embodiment is secured.
  • the SAR image data 641 is data of a SAR image acquired from the SAR image storage unit 204.
  • the phase array data 642 is data of the phase array of each reflection point (pixel) acquired from the phase array storage unit 201.
  • the stable reflection point candidate data 643 is data of stable reflection point candidates of each reflection point (pixel) acquired from the stable reflection point candidate table 530.
  • the cluster data 644 is data of a cluster into which each reflection point (pixel) acquired from the cluster table 540 is classified.
  • the storage 650 stores a database, various parameters, or the following data or programs necessary for realizing the present embodiment.
  • the phased array storage unit 201 stores the phased array of each reflection point (pixel) as shown in FIG. 5B.
  • the stable reflection point candidate table 530 stores, as shown in FIG. 5C, whether or not each reflection point (pixel) is a stable reflection point candidate.
  • the cluster table 540 stores reflection points (pixels) classified into each cluster as shown in FIG. 5D.
  • the phase statistics table 550 stores phase statistics (phase change) of each cluster as shown in FIG. 5E.
  • the phase correlation table 560 stores, as shown in FIG. 5F, the correlation between the phase array of each reflection point (pixel) and the phase statistic (phase change) of each cluster.
  • the storage 650 stores the following programs.
  • the image processing device control program 651 is a program that controls the overall processing of the image processing device 200 of the present embodiment.
  • the phase alignment storage module 652 is a module that stores the phase alignment of each reflection point (pixel) across a plurality of SAR images.
  • the stable reflection point candidate identification module 653 is a module for identifying stable reflection point candidates based on the stability of data of each reflection point (pixel).
  • the clustering module 654 is a module that clusters each reflection point (pixel) into a plurality of clusters based on the phase array of each reflection point (pixel).
  • the phase statistics computing module 655 is a module for computing a representative phase feature (phase change) of each cluster.
  • the phase correlation operation module 656 is a module that calculates the correlation between the phase array of each reflection point (pixel) and the phase statistic (phase change) of each cluster.
  • the input / output interface 660 provides an interface for controlling data input / output with an input / output device.
  • the input / output interface 660 is connected to a display unit 661 that displays the state of the image processing apparatus 200 and an operation unit 662 that inputs an operation to the image processing apparatus 200.
  • a storage medium 663 storing a SAR image may be connected.
  • FIG. 7 is a flowchart showing the processing procedure of the image processing apparatus 200 according to the present embodiment. This flowchart is executed by the CPU 610 in FIG. 6 using the RAM 640, and implements the functional components in FIG. 2 and FIG.
  • step S701 the image processing apparatus 200 generates and accumulates a phase array of each reflection point (pixel) from a plurality of SAR images.
  • step S703 the image processing apparatus 200 clusters each reflection point (pixel) of the SAR image into a plurality of clusters based on the phase array.
  • step S705 the image processing apparatus 200 calculates a statistic representing a representative phase change of each cluster from the phase array of reflection points (pixels) in each cluster.
  • step S 707 the image processing apparatus 200 calculates the correlation between the statistic representing the representative phase change of each cluster and the phase array of each reflection point (pixel) of the entire image.
  • the image processing apparatus 200 outputs the calculated correlation as it is or collectively at step S709. It can be judged that the reflection point (pixel) having a strong correlation is less affected by noise and seems to be a stable reflection point suitable for displacement analysis, which is a stable evaluation index.
  • the correlation between the phase change in each cluster of the plurality of clusters and the phase arrangement of each pixel is used to evaluate the stable reflection point, an object whose displacement is nonlinear or an object whose elevation is high For a target with a large displacement, the likelihood of a stable reflection point can be evaluated without being affected by a large phase change.
  • the pixel is once divided into clusters having similar phase changes, and statistics on representative phase changes in each of the clusters are calculated. Thereby, the phase change which removed the influence by noise is calculated
  • the phase change including noise in each pixel of the whole image and the noise in each cluster are included. It is understood how uncorrelated phase changes are correlated. Therefore, a pixel with strong correlation can be identified as a stable reflection point with little noise.
  • the image processing apparatus according to the present embodiment differs from the second embodiment in that it has a noise removing unit that removes noise from the calculated phase statistics.
  • the other configurations and operations are similar to those of the second embodiment, and therefore the same configurations and operations are denoted by the same reference numerals and the detailed description thereof is omitted.
  • FIG. 8 is a block diagram showing a functional configuration of the phase correlation calculation unit 803 of the image processing apparatus 800 according to the present embodiment.
  • the configuration of the image processing apparatus 800 is obtained by replacing the phase correlation calculation unit 203 in FIG. 2 with a phase correlation calculation unit 803.
  • the phase correlation calculation unit 803 in the present embodiment includes a phase statistics calculation unit 431, a noise removal unit 833, and a phase correlation calculation unit 432. That is, the phase correlation calculation unit 203 according to the second embodiment has only the point that the noise removal unit 833 is provided between the phase statistics calculation unit 431 and the phase correlation calculation unit 432 in the present embodiment. It is different from
  • the phase statistic calculator 431 calculates a statistic representing a typical phase change of a cluster.
  • the noise removing unit 833 removes the influence of noise from the statistic calculated by the phase statistic computing unit 431.
  • the phase correlation operation unit 432 calculates the correlation with each reflection point (pixel) using the statistic whose noise influence has been removed by the noise removal unit 833.
  • the noise removing unit 833 removes the influence of noise from the statistic related to the representative phase change of the cluster calculated by the phase statistic calculating unit 431.
  • the statistics concerning the representative phase change of each cluster is calculated in the phase statistics calculation unit 431, the influence of noise to some extent is eliminated.
  • the robustness against noise can be further enhanced by using prior knowledge on the phase alignment that can be specified in each of the pixels in each cluster and removing noise that does not match prior knowledge as noise. it can.
  • the matrix S C (r) which is a statistic calculated by equation (2), corresponds to one eigenvector corresponding to a typical phase change of a cluster and noise by eigenvalue decomposition It can be seen that it has the property of being divided into the remaining eigenvectors. Therefore, for example, the eigenvalue decomposition of the matrix S C (r) which is the feature quantity calculated by equation (2), and the largest eigenvector is the eigenvector corresponding to the representative phase change of the cluster, Let eigenvectors be attributed to noise. Thus, by separating the matrix S C (r) into the matrix corresponding to the largest eigenvector and the other matrices, it is possible to obtain statistics from which the influence of noise is removed.
  • the eigenvalue decomposition of the matrix S C (r) which is the feature quantity calculated by equation (2), and the number of eigenvectors having large eigenvalues are assumed to be the eigenvectors corresponding to the representative phase change of the cluster , And other eigenvectors are caused by noise.
  • the matrix S C (r) by separating the matrix S C (r) into a matrix corresponding to a number of eigenvectors having large eigenvalues and the other matrices, it is possible to obtain a statistic in which the influence of noise is removed.
  • the influence of noise may be removed by assuming a model that is gradual with respect to time with respect to displacement.
  • the matrix S C (r) which is a statistic calculated by Equation (2)
  • phase comparison is performed between the elements of the phase array corresponding to the image close to the imaging time as the image acquisition time.
  • the absolute value of the component being performed is increased.
  • the absolute value of the component whose phase comparison is performed between the elements of the phase array corresponding to the image distant from the photographing time becomes small.
  • noise with respect to components other than the diagonal components of the matrix S C (r), there is a property that the absolute value is increased regardless of whether the time is distant or not.
  • FIG. 9 is a diagram showing the configuration of the noise removal table 910 according to the present embodiment.
  • the noise removal table 910 is used by the noise removal unit 833 to remove noise of the phase statistics (phase change) calculated by the phase statistics calculation unit 431.
  • the configuration of the noise removal table 910 is not limited to that shown in FIG.
  • FIG. 10 is a flowchart showing the processing procedure of the image processing apparatus 800 according to the present embodiment. This flowchart is executed by the CPU 610 in FIG. 6 to which the noise removal table 910 is added using the RAM 640, and implements the functional configuration units in FIG. 2 and FIG. In FIG. 10, the same steps as in FIG. 7 will be assigned the same step numbers, and redundant description will be omitted.
  • Steps S207 to S705 operate in the same manner as in FIG. Thereafter, in step S1006, the image processing apparatus 800 removes the influence of noise from the statistic in the cluster calculated in step S705.
  • step S 707 the image processing apparatus 800 matches the representative phase change of each cluster and the phase change of the evaluation target in the same manner as in FIG. 7 using the statistic from which the influence of noise is removed. Calculate the correlation. Therefore, as the correlation is stronger, the influence of noise is smaller and it can be judged as a stable reflection point suitable for displacement analysis.
  • the present embodiment it is possible to specify stable reflection points necessary for performing displacement analysis on an object causing nonlinear displacement, an object having a high altitude, and an object having a large displacement.
  • the present embodiment includes a portion for removing the influence of noise as a configuration, and the stable reflection point can be specified more robustly.
  • the image processing apparatus according to the present embodiment differs from the second embodiment in that the correlation operation unit and the phase statistics operation unit of the phase correlation calculation unit are arranged in reverse.
  • the other configurations and operations are similar to those of the second embodiment, and therefore the same configurations and operations are denoted by the same reference numerals and the detailed description thereof is omitted.
  • FIG. 11 is a block diagram showing the functional configuration of the phase correlation calculation unit 1103 of the image processing apparatus 1100 according to this embodiment.
  • the configuration of the image processing apparatus 1100 is obtained by replacing the phase correlation calculation unit 203 in FIG. 2 with a phase correlation calculation unit 1103.
  • the phase correlation calculation unit 1103 in the present embodiment includes a phase correlation calculation unit 1131 and a phase statistics calculation unit 1132.
  • the phase correlation operation unit 1131 calculates the phase array stored in the phase array storage unit 201 in each pixel of the entire image and the phase array at each reflection point (pixel) forming each cluster clustered by the clustering unit 202. Calculate the correlation.
  • the phase statistics amount calculation unit 1132 calculates, for example, within a cluster, the correlation with the phase change calculated at each reflection point (pixel) of the entire image at each reflection point (pixel) constituting each cluster. Calculate statistics such as average value of correlation.
  • the phase correlation calculation unit 1103 calculates statistics of the correlation of the phase change between each of the evaluation targets and each of the pixels constituting each cluster.
  • the phase correlation operation unit 1131 calculates the correlation of the phase change between each of the pixels in the entire image and each of the pixels constituting the cluster.
  • the correlation between the phase changes of the pixels is calculated, for example, as follows. Temporarily, the index of the pixel for which the correlation is to be specified is i, the cluster is C, and the index of the pixel in cluster C is j ⁇ C.
  • the elements of the phased-array in the pixel i in the column vector s i if containing the elements in the phase sequence observed in the pixel j in column vector s j, the correlation corr phase change between the pixel i and the pixel j (i, j) is calculated by the following equation (5).
  • a monotonous non-decreasing function may be applied to a positive number. For example, it may be squared. Alternatively, a function that is “0 (zero)” may be applied to inputs below a certain threshold. By applying such a function, it becomes possible to perform evaluation by removing the influence of the extremely weakly correlated pixel.
  • the phase statistic calculator 1132 receives the correlation between the phase array of each pixel of the entire image and the phase array of each pixel constituting each cluster, which is calculated by the phase correlation calculator 1131. Then, a statistic representing the correlation between the phase arrangement at each pixel of the entire image and the phase change as a whole at the pixels constituting each cluster is calculated. Assuming that i of the pixel for which correlation is to be calculated and the cluster of which the correlation is to be calculated is C, the phase statistic calculator 1132 calculates a statistic corr (i, C) representing the correlation with the entire cluster, For corr (i, j) obtained in (5), statistics are calculated for all pixels j ⁇ C included in the cluster C.
  • the average value for jrC of corr (i, j) is calculated as in the following equation (6).
  • other statistics such as maximum value, mode value and median value may be used.
  • weighted evaluation may be performed using a monotonous non-increasing load function according to the distance.
  • w (d) is a monotonically non-increasing function with respect to d
  • r i and r j are specified as the positions of pixels i and j, respectively, as shown in equation (7) below
  • You may The statistic representing the correlation between the phase change in each of the pixels of the entire image and the phase change in the pixels constituting each cluster calculated as described above is calculated by the phase correlation calculation unit 203 in the second embodiment.
  • the correlation between the phase change in each of the identified pixels and the representative phase change of each cluster can be treated almost the same.
  • FIG. 12A is a diagram showing the configuration of the phase correlation table 1210 according to the present embodiment.
  • the phase correlation table 1210 is used by the phase correlation operation unit 1131 to calculate the correlation between the phase array of each pixel of the entire image and the phase array of each pixel constituting the cluster.
  • the configuration of the phase correlation table 1210 is not limited to that shown in FIG. 12A.
  • the phase correlation table 1210 calculates phase correlation values of all combinations between the phase array group 561 of each reflection point (pixel) and the phase array group 1262 of each reflection point (pixel) classified into each cluster. It is memorized.
  • FIG. 12B is a diagram showing the configuration of the phase statistic table 1220 according to the present embodiment.
  • the phase statistic table 1220 is obtained by correlating the phase array of each pixel calculated by the phase statistic calculator 1132 in the phase correlation calculator 1131 with the phase array of each pixel constituting each cluster, and It is used to calculate a statistic that represents the correlation with the phase change of each cluster.
  • the configuration of the phase statistic table 1220 is not limited to FIG. 12B.
  • the phase statistics table 1220 stores correlation statistics 1222, 1223, ... of each cluster in correspondence with the reflection point ID (pixel ID) 1221.
  • the correlation statistic 1222, 1223, ... of each cluster is a correlation group between each reflection point ID (pixel ID) 1221 and each reflection point (pixel) in each cluster and a correlation statistic calculated from the correlation group It stores correlation averages and correlation weighted averages.
  • FIG. 13 is a flowchart showing the processing procedure of the image processing apparatus 1100 according to this embodiment. This flowchart is executed by the CPU 610 in FIG. 6 having the phase correlation table 1210 and the phase statistics table 1220 using the RAM 640, and implements the functional configuration units in FIG. 2 and FIG. In FIG. 13, the same steps as in FIG. 7 will be assigned the same step numbers, and redundant description will be omitted. Steps S701 and S703 operate in the same manner as in FIG.
  • step S1305 the image processing apparatus 1100 calculates the correlation between the phase array at each pixel of the entire image and the phase array at each pixel constituting each cluster.
  • step S1307 the image processing apparatus 1100 calculates correlation statistics for each cluster with respect to the correlation between the phase array in each pixel and the phase array in each pixel constituting each cluster, which is calculated in step S1305.
  • step S1309 the image processing apparatus 1100 outputs the calculated correlation statistic as it is or collectively. The stronger the statistic about the correlation in the cluster finally calculated, the smaller the influence of noise can be determined as a stable reflection point suitable for displacement analysis.
  • each cluster is set for each of the pixels in the entire image.
  • the evaluation can be performed by removing the influence of pixels having extremely different phase changes among the pixels to be configured. Such an evaluation method is useful, for example, when an algorithm of clustering does not divide a cluster sufficiently and is identified in a state where a plurality of clusters are mixed.
  • the image processing apparatus according to the present embodiment is different from the second to fourth embodiments in that the calculation result of the phase correlation calculation unit is fed back to the stable reflection point candidate specifying unit.
  • the other configurations and operations are similar to those of the second embodiment, and therefore the same configurations and operations are denoted by the same reference numerals and the detailed description thereof is omitted.
  • FIG. 14 is a block diagram showing a functional configuration of the image processing apparatus 1400 according to the present embodiment.
  • the same functional components as in FIG. 2 will be assigned the same reference numerals and overlapping descriptions will be omitted.
  • the phase correlation calculation unit 203 and the stable reflection point candidate identification unit 205 in FIG. 2 are replaced with a phase correlation calculation unit 1403 and a stable reflection point candidate identification unit 1405, and those are connected by a feedback line.
  • the stable reflection point candidate identification unit 1405 identifies a pixel that can be a stable reflection point.
  • the stable reflection point candidate identification unit 1405 in the present embodiment differs from the stable reflection point candidate identification unit 205 in the second embodiment in that stable reflection point candidates are again identified based on the correlation calculated by the phase correlation calculation unit 1403. , Has the ability to re-identify.
  • the clustering unit 202 again performs clustering with the stable reflection candidate specified by the stable reflection point candidate specifying unit 1405 again as an input.
  • the phase correlation calculation unit 1403 calculates the correlation with the phase change representing each cluster clustered again by the clustering unit 202 with respect to the phase array in each pixel of the entire image. As described above, the correlation calculated by the phase correlation calculation unit 1403 is fed back to the stable reflection point candidate identification unit 1405, and identification of stable reflection point candidates is performed again.
  • the phase correlation calculating unit 1403 is the phase correlation calculating unit 203 in the second embodiment described above, or the phase correlation calculating unit 803 in the third embodiment, or the phase correlation calculating unit 1103 in the fourth embodiment. Either may be used.
  • FIG. 15 is a diagram showing the configuration of the stable reflection point candidate table 1530 according to the present embodiment.
  • the stable reflection point candidate table 1530 is used by the stable reflection point candidate identification unit 1405 to repeat identification of stable reflection point candidates among the respective reflection points (pixels) in a plurality of SAR images based on correlation. .
  • the same components as in FIG. 5C are denoted by the same reference numerals, and redundant description will be omitted.
  • the configuration of the stable reflection point candidate table 1530 is not limited to FIG.
  • the stable reflection point candidate table 1530 is changed based on the phase correlation calculation result 1532 fed back from the phase correlation calculation unit 1403 and the phase correlation calculation result 1532 in association with each reflection point ID (pixel ID) 531.
  • the reflection point candidate flag 1533 is stored.
  • the stable reflection point candidate flag 1533 is a stable reflection point candidate in the case of “1” and is not a stable reflection point candidate in the case of “0”.
  • the reflection point (pixel) where the reflection point ID (pixel ID) 531 is “A0001-0002” is initially a stable reflection point candidate but is excluded from the stable reflection point candidates based on the phase correlation calculation result 1532 There is.
  • FIG. 16 is a flowchart showing the processing procedure of the image processing apparatus 1400 according to this embodiment. This flowchart is executed by the CPU 610 in FIG. 6 having the stable reflection point candidate table 1530 using the RAM 640, and implements the functional configuration unit in FIG. In FIG. 14, the same steps as in FIG. 7 will be assigned the same step numbers, and redundant description will be omitted. Step S701 operates in the same manner as in FIG.
  • step S1603 the image processing apparatus 1400 specifies a pixel that can be a stable reflection point as a stable reflection point candidate.
  • step S1605 the image processing apparatus 1400 clusters stable reflection point candidates based on the phase array.
  • step S1607 the image processing apparatus 1400 calculates the correlation between the representative phase change in each cluster and the phase array in each pixel of the entire image.
  • step S1607 one of (S705, S707) in FIG. 7, (S705, S1006, S707) in FIG. 10 or (S1305, S1307) in FIG. 13 is executed.
  • step S1609 the image processing apparatus 1400 determines whether the identified cluster or the correlation of the phase with the identified cluster satisfies the end condition.
  • the termination condition for example, it may be used whether or not a predetermined number of repetitions have been performed.
  • all the stable reflection point candidates satisfy the condition as the stable reflection point when the correlation in the pixel which becomes the stable reflection point candidate is stronger than a certain level, that is, stable reflection point candidate identification It may be determined that the end condition is satisfied, assuming that the unit 1405 correctly identifies the stable reflection point.
  • the image processing apparatus 1400 feeds back the current phase correlation calculation result to the stable reflection point candidate identification unit 1405 in step S1611, and updates the identification of stable reflection point candidates again from step S1603. Or repeat and repeat the process. If the end condition is satisfied, the image processing apparatus 1400 outputs the phase correlation calculation result of the narrowed stable reflection point candidate in step S1613 and ends the process.
  • the likelihood of a stable reflection point can be more accurately evaluated by performing identification of stable reflection point candidates again in accordance with the strength of correlation specified for each of the pixels in the entire image.
  • the correlation can be calculated.
  • clusters of stable reflection point candidates are extracted by extracting clusters for stable reflection point candidates judged to have less noise, and by making identification results of the stable reflection point candidates more accurate. It is possible to reduce the influence of noise when correlating the phase change representative of the phase change with the phase change in each of the pixels of the entire image.
  • the present invention can be applied to a program for analyzing displacement of an object such as a high-rise building or a bridge which can have nonlinear displacement or large displacement using a SAR image. Moreover, if a sound wave is used, it is applicable also to synthetic aperture sonar. By extracting pixels with stable phase in a cluster according to the present invention, it is possible to find objects that have been buried in noise only by reflection intensity.
  • the present invention is also applicable to phase images captured by a Time Of Flight camera or the like. By applying the present invention, it becomes possible to distinguish whether the fluctuation of the observation value is caused by noise even for the object whose position moves a little, and to extract the solid shape of the object robustly. Will be able to In addition, the present invention can be applied to phase images using coherent light interference.
  • the present invention may be applied to a system configured of a plurality of devices or to a single device. Furthermore, the present invention is also applicable to the case where an image processing program for realizing the functions of the embodiments is supplied to a system or apparatus directly or remotely. Therefore, in order to realize the functions of the present invention on a computer, a program installed on the computer, a medium storing the program, and a WWW (World Wide Web) server for downloading the program are also included in the scope of the present invention. .
  • a non-transitory computer readable medium storing a program that causes a computer to execute at least the processing steps included in the above-described embodiment is included in the scope of the present invention.
  • Phase array storage means for storing a phase array of each pixel across a plurality of images; Clustering means for clustering the respective pixels into a plurality of clusters based on the phase array; A phase correlation calculating unit that calculates a correlation between a phase change in each of the plurality of clusters and a phase array in each of the pixels;
  • An image processing apparatus comprising: (Supplementary Note 2) The phase correlation calculation means First phase statistic calculation means for calculating representative phase change statistics within each cluster; First phase correlation operation means for calculating the correlation between the representative phase change and the phase array in each of the pixels;
  • the image processing apparatus according to claim 1 further comprising: (Supplementary Note 3) The image processing apparatus according to claim 2, wherein said first phase statistic computing means uses an average of phase arrangements of a plurality of pixels in each cluster as the representative statistic of phase change.
  • phase correlation calculation means further comprises noise removal means for removing the influence of noise on said representative phase change in a cluster.
  • the noise removal means eigenvalue decomposes the representative phase change statistic calculated by the first phase statistic calculation means, and sets the largest eigenvector as the representative phase change from which the influence of noise is removed.
  • the noise removal means eigen value decomposes the representative phase change statistics calculated by the first phase statistics calculation means, and the several eigenvectors from the largest eigenvector are removed from the influence of noise.
  • the phase correlation calculation means Second phase correlation operation means for calculating the correlation between the phase arrangement of the pixels in each cluster and the phase arrangement of each of the pixels; Second phase statistics computing means for computing representative correlation statistics within each cluster;
  • the image processing apparatus according to appendix 10 wherein the second phase statistic computing means performs weighting corresponding to the position of each pixel. (Supplementary Note 12) 11.
  • the image processing apparatus according to any one of appendices 1 to 11, wherein the phase correlation calculation means includes a process of putting together the calculated correlations.
  • the image processing apparatus according to any one of appendices 1 to 12, further comprising stable reflection point candidate specifying means for specifying a candidate pixel that becomes a stable reflection point where reflection is stable from each pixel across the plurality of images.
  • the stable reflection point candidate specifying unit updates the specification of the candidate reflection point candidate stable reflection point again based on the correlation calculated by the phase correlation calculation unit.
  • Image processing method including: (Appendix 18) A phased array storage step of storing the phased array of each pixel across a plurality of images in a phased array storage means; Clustering the respective pixels into a plurality of clusters based on the phase array; Calculating a correlation between a phase change in each of the plurality of clusters and a phase array in each of the pixels;
  • An image processing program that causes a computer to execute.
  • An image processing apparatus for identifying stable reflection points based on a plurality of SAR images; A SAR image analyzer for analyzing a SAR image using data of the stable reflection point; SAR image analysis system comprising The image processing apparatus at least Phase array storage means for storing a phase array of each pixel across the plurality of SAR images; Clustering means for clustering the respective pixels into a plurality of clusters based on the phase array; A phase correlation calculating unit that calculates a correlation between a phase change in each of the plurality of clusters and a phase array in each of the pixels; SAR image analysis system with (Supplementary Note 20) An image processing step for identifying a stable reflection point based on a plurality of SAR images; A SAR image analysis step of analyzing a SAR image using data of the stable reflection point; SAR image analysis method including The image processing step includes at least A phase array storage step of storing the phase array of each pixel across the plurality of SAR images in a phase array storage unit; Clustering

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Image Analysis (AREA)

Abstract

本発明は、変位が非線形である対象や、標高が高い対象、変位が大きな対象について、大きな位相変化に影響されずに安定反射点らしさの評価指標を生成する画像処理装置に関する。本画像処理装置は、複数の画像に渡る各画素の位相配列を蓄積する位相配列蓄積部と、位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリング部と、複数のクラスタの各クラスタ内での位相変化と各画素における位相配列との相関を算出する位相相関算出部と、を備える。

Description

ＳＡＲ画像解析システム、画像処理装置、画像処理方法および画像処理プログラム
　本発明は、ＳＡＲ画像解析システム、画像処理装置、画像処理方法および画像処理プログラムに関する。
　上記技術分野において、特許文献１には、ＳＡＲ（Synthetic Aperture Rader)画像解析システムの一例が記載されている。特許文献１のＳＡＲ画像解析システムでは、各安定反射点候補において観測された位相に対して、モデルの当て嵌めを行い、当て嵌めるモデルとして変位が時間に対して線形に起こるというモデルを仮定する。次に、モデル当て嵌めが、どの程度ずれているかを特定し、モデルに当てはまらない位相の変化が全て標高や変位、大気の状態に起因しないノイズによるものとし、誤差の少なかった点を安定反射点として特定する。最後に、安定反射点を入力として、それらの点上で変位と標高を特定する。以上により、ノイズの影響が少ない安定反射点と、各安定反射点における精密な変位・標高が出力される。
特表２００３－５００６５８号公報（国際特許公開ＷＯ００／７２０４５）
　しかしながら、上記文献に記載の技術では、変位が非線形である対象や、標高が高い対象、変位が大きな対象について、変位計測が上手くできないという問題がある。
　本発明の目的は、上述の課題を解決する技術を提供することにある。
　上記目的を達成するため、本発明に係る画像処理装置は、
　複数の画像に渡る各画素の位相配列を蓄積する位相配列蓄積手段と、
　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリング手段と、
　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出手段と、
　を備える。
　上記目的を達成するため、本発明に係る画像処理方法は、
　複数の画像に渡る各画素の位相配列を位相配列蓄積手段に蓄積する位相配列蓄積ステップと、
　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリングステップと、
　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出ステップと、
　を含む。
　上記目的を達成するため、本発明に係る画像処理プログラムは、
　複数の画像に渡る各画素の位相配列を位相配列蓄積手段に蓄積する位相配列蓄積ステップと、
　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリングステップと、
　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出ステップと、
　をコンピュータに実行させる。
　上記目的を達成するため、本発明に係るＳＡＲ画像解析システムは、
　複数のＳＡＲ画像に基づいて安定反射点を特定するための画像処理装置と、
　前記安定反射点のデータを用いてＳＡＲ画像を解析するＳＡＲ画像解析装置と、
　を備えるＳＡＲ画像解析システムであって、
　前記画像処理装置は、少なくとも、
　　前記複数のＳＡＲ画像に渡る各画素の位相配列を蓄積する位相配列蓄積手段と、
　　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリング手段と、
　　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出手段と、
　を備える。
　上記目的を達成するため、本発明に係るＳＡＲ画像解析方法は、
　複数のＳＡＲ画像に基づいて安定反射点を特定するための画像処理ステップと、
　前記安定反射点のデータを用いてＳＡＲ画像を解析するＳＡＲ画像解析ステップと、
　を含むＳＡＲ画像解析方法であって、
　前記画像処理ステップは、少なくとも、
　　前記複数のＳＡＲ画像に渡る各画素の位相配列を位相配列蓄積手段に蓄積する位相配列蓄積ステップと、
　　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリングステップと、
　　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出ステップと、
　を含む。
　本発明によれば、変位が非線形である対象や、標高が高い対象、変位が大きな対象について、大きな位相変化に影響されずに安定反射点らしさの評価指標を生成することができる。
本発明の第１実施形態に係る画像処理装置の機能構成を示すブロック図である。 本発明の第２実施形態に係る画像処理装置の機能構成を示すブロック図である。 本発明の第２実施形態に係る画像処理装置を含むＳＡＲ画像解析システムの概要を示す図である。 本発明の第２実施形態に係る画像処理装置を含むＳＡＲ画像解析システムの機能構成を示す図である。 本発明の第２実施形態に係る位相相関算出部の機能構成を示すブロック図である。 本発明の第２実施形態に係るＳＡＲ画像記憶部の構成を示す図である。 本発明の第２実施形態に係る位相配列蓄積部の構成を示す図である。 本発明の第２実施形態に係る安定反射点候補テーブルの構成を示す図である。 本発明の第２実施形態に係るクラスタテーブルの構成を示す図である。 本発明の第２実施形態に係る位相統計量テーブルの構成を示す図である。 本発明の第２実施形態に係る位相相関テーブルの構成を示す図である。 本発明の第２実施形態に係る安定反射点テーブルの構成を示す図である。 本発明の第２実施形態に係る画像処理装置のハードウェア構成を示すブロック図である。 本発明の第２実施形態に係る画像処理装置の処理手順を示すフローチャートである。 本発明の第３実施形態に係る位相相関算出部の機能構成を示すブロック図である。 本発明の第３実施形態に係るノイズ除去テーブルの構成を示す図である。 本発明の第３実施形態に係る画像処理装置の処理手順を示すフローチャートである。 本発明の第４実施形態に係る位相相関算出部の機能構成を示すブロック図である。 本発明の第４実施形態に係る位相相関テーブルの構成を示す図である。 本発明の第４実施形態に係る位相統計量テーブルの構成を示す図である。 本発明の第４実施形態に係る画像処理装置の処理手順を示すフローチャートである。 本発明の第５実施形態に係る画像処理装置の機能構成を示すブロック図である。 本発明の第５実施形態に係る安定反射点候補テーブルの構成を示す図である。 本発明の第５実施形態に係る画像処理装置の処理手順を示すフローチャートである。 前提技術のＳＡＲ画像解析システムの機能構成を示す図である。
　以下に、図面を参照して、本発明の実施の形態について例示的に詳しく説明する。ただし、以下の実施の形態に記載されている構成要素は単なる例示であり、本発明の技術範囲をそれらのみに限定する趣旨のものではない。
　［第１実施形態］
　本発明の第１実施形態としての画像処理装置１００について、図１を用いて説明する。画像処理装置１００は、複数の画像に渡る各画素の位相配列から安定反射点らしさの評価指標を生成する装置である。
　図１に示すように、画像処理装置１００は、位相配列蓄積部１０１と、クラスタリング部１０２と、位相相関算出部１０３と、を含む。位相配列蓄積部１０１は、複数の画像に渡る各画素の位相配列を蓄積する。クラスタリング部１０２は、位相配列に基づいて、各画素を複数のクラスタにクラスタリングする。位相相関算出部１０３は、複数のクラスタの各クラスタ内での位相変化と各画素における位相配列との相関を算出する。
　本実施形態によれば、複数のクラスタの各クラスタ内での位相変化と各画素における位相配列との相関を算出するので、変位が非線形である対象や、標高が高い対象、変位が大きな対象について、大きな位相変化に影響されずに安定反射点らしさの評価指標を生成することができる。
　［第２実施形態］
　次に、本発明の第２実施形態に係る画像処理装置について説明するが、本実施形態の画像処理装置の特徴を明瞭とするため、まず、前提技術として、ＳＡＲ画像解析システムの一例を説明する。
　《前提技術の説明》
　図１７は、特許文献１に記載されている、前提技術のＳＡＲ画像解析システム１７００の機能構成を示す図である。
　図１７に示すように、ＳＡＲ画像解析システム１７００は、ＳＡＲ画像記憶部１７０１と、位相特定部１７０２と、安定反射点候補特定部１７０３と、モデル当て嵌め部１７０４と、モデル誤差特定部１７０５と、安定反射点特定部１７０６と、変位・標高特定部１７０７と、から構成される。
　このような構成を有するＳＡＲ画像解析システム１７００はつぎのように動作する。まず、位相特定部１７０２は、同じ場所を異なる時期に撮影した複数のＳＡＲ画像をＳＡＲ画像記憶部１７０１から入力し、各々の画素における位相の変化を特定する。安定反射点候補特定部１７０３は、同じ場所を異なる時期に撮影した複数のＳＡＲ画像を入力とし、反射が安定している可能性の高い画素である、安定反射点候補を特定する。安定反射点候補としては、例えば、特許文献１では、同じ場所を異なる時期に撮影した複数のＳＡＲ画像を入力とし、複数のＳＡＲ画像にわたって反射の強度の変化が大きくないような点を、安定反射点候補として特定している。次に、モデル当て嵌め部１７０４は、安定反射点候補特定部で特定された安定反射点候補を入力とし、各安定反射点候補において観測された位相に対して、モデルの当て嵌めを行う。ここで、ＳＡＲ画像上の位相は、標高や変位、大気の状態等に応じて変化をする量である。すなわち、当て嵌めるモデルは、撮影対象の標高・変位、大気の状態等が、ＳＡＲ画像を撮影する時間によって、どのように変化するかを表すものである。特許文献１では、変位が時間に対して線形に起こるというモデルを仮定している。次に、モデル誤差特定部１７０５では、モデル当て嵌め部で行った当て嵌めが、どの程度ずれているかを特定する。特許文献１では、モデルに当てはまらない位相の変化が全て標高や変位、大気の状態に起因しないノイズによるものだと仮定し、モデルからの誤差の大きさはノイズの大きさを表すとしている。次に、安定反射点特定部１７０６は、モデル誤差特定部で特定された誤差を入力とし、誤差の少なかった点を安定反射点として特定する。最後に、変位・標高特定部１７０７は、前記の安定反射点を入力として、それらの点上で変位と標高を特定する。以上により、ノイズの影響が少ない安定反射点と、各安定反射点における精密な変位・標高が出力される。
　（前提技術の課題）
　しかしながら、このようなＳＡＲ画像解析システム１７００では、変位が非線形である対象や、標高が高い対象、変位が大きな対象について、変位計測が上手くできないという問題がある。その第１の原因は、線形な変位のモデルからの誤差の大きさを元にして安定かどうかを評価しているために、非線形な変位をしている物が安定反射点ではないと識別されてしまうことにある。もう１つの原因は、位相の当て嵌め自体が計算的に困難だということにある。位相は円周上の位置のような量であり、例えば一周分の位相変化を起こしても、あるいは２周分の位相変化を起こしても、変化後に観測すると原点に戻ってきているために、位相変化が全くない場合との区別がつかない。このため、大きな変位がある場合や標高が高い場合等、位相の変化が大きくなる場合に、正しい位相変化が分からなくなり、モデルが上手く当て嵌められなくなる。
　《本実施形態の説明》
　上記課題を解決するために、本実施形態に係る画像処理装置は、同じ場所を異なる時期に撮影した複数のＳＡＲ画像を入力として画素の各々における位相配列を作成する。その位相配列を入力として、画素を位相が同期しているクラスタに分割する。そして、クラスタリング結果とＳＡＲ画像上の画素の各々における位相配列を入力として、クラスタ各々の代表的な位相変化と、ＳＡＲ画像上の画素の各々における位相変化との相関を算出する。ここで、位相配列とは入力となるＳＡＲ画像のうち基準となる１枚に対して、その各々の画素上で作る配列であり、同じ場所を異なる時期に撮影した複数のＳＡＲ画像に渡り、基準となる１枚のＳＡＲ画像の各々の画素と同一地点であるような画素の各位相を配列の各要素として持つ。位相変化が類似した点群に分割し、その平均的な位相変化と評価対象となる点の位相変化が相関しているのかどうかを評価することにより、評価対象となる安定反射点候補がどの程度ノイズの影響を受けていないのかを判断することができる。以上の過程にはモデルの当て嵌めを行う処理は含まれていないため、高層ビルなどの位相のモデルを当てはめにくい対象についても安定反射点を特定することができる。
　《画像処理装置の機能構成》
　図２は、本実施形態に係る画像処理装置２００の機能構成を示すブロック図である。
　図２を参照すると、本実施形態の画像処理装置２００は、最小構成として、位相配列蓄積部２０１と、クラスタリング部２０２と、位相相関算出部２０３とを備える。また、画像処理装置２００は、ＳＡＲ画像記憶部２０４と、安定反射点候補特定部２０５、安定反射点特定部２０６を含んでもよい。
　位相配列蓄積部２０１は、ＳＡＲ画像記憶部２０４からの複数のＳＡＲ画像に基づいて、画像全体の画素の各々における位相を記録した配列である位相配列を蓄積する。クラスタリング部２０２は、クラスタリング手法をもちいて、画素を位相が類似しているクラスタに分割する。クラスタリング手法としては、例えば、非特許文献１に記載の方法を用いてよいが、これに限定されない。位相相関算出部２０３は、クラスタリング結果と画像全体の画素各々における位相配列に基づいて、各クラスタ内の代表的な位相変化と、画素の各々における位相変化との、相関を算出する。
　安定反射点候補特定部２０５は、安定反射点になりうる画素を特定する。安定反射点候補特定部２０５を備える場合、クラスタリング部２０２は、安定反射点候補のみをクラスタリングする。安定反射点候補特定部２０５を備えてない場合、クラスタリング部２０２は、画像全体の画素をクラスタリングする。安定反射点特定部２０６を備える場合、安定反射点特定部２０６は、位相相関算出部２０３において算出された各画素と各クラスタとの相関を元に、所定の閾値を上回る相関を持つ画素を安定反射点とするなどを条件として、安定反射点を特定する。
　続いて、本実施形態における画像処理装置の各機能構成について説明する。
　（位相配列蓄積部）
　位相配列蓄積部２０１は、画像全体の画素の各々における位相配列を蓄積する。位相配列の蓄積に当たって、位相配列蓄積部２０１は、画素ごとに、複数のＳＡＲ画像の各画素における位相を要素とする配列を生成する。この位相は、例えば、基準となるＳＡＲ画像と他のＳＡＲ画像との間での位相差として求めてもよい。この場合は、入力となる複数のＳＡＲ画像のうち、基準となるＳＡＲ画像があらかじめ定められる。そして、位相配列蓄積部２０１は、基準となるＳＡＲ画像と他のＳＡＲ画像の各々との位相差を求めて、それらを各画素における位相配列として蓄積する。また、別の例として、位相配列蓄積部２０１は、基準となるＳＡＲ画像を定めずに、複数のＳＡＲ画像の位相を時系列の順などに格納したものを位相配列としてもよい。各々の位相は、例えば絶対値が“１”となるように正規化された複素数の形式で格納される。また、位相配列蓄積部２０１は、ＤＥＭ(Digital Elevation Model)等を用いて、対象の標高に依存して変化しえる位相の変化を除去してもよいし、大気の揺らぎに依存した位相変化を公知の手法を用いるなどして除去してもよい。
　（クラスタリング部）
　クラスタリング部２０２は、画像全体の画素各々を、類似した位相を持つクラスタに分割する。例えば、非特許文献１に記載の方法を用いれば、画像全体の画素の各々を類似した位相を持つクラスタに分割することができる。
　（位相相関算出部）
　位相相関算出部２０３は、各クラスタの代表的な位相変化と、画像全体の各画素における位相配列との相関を算出する。算出された相関は、強いほど安定反射点らしいと判断するための評価指標として用いることができる。
　《ＳＡＲ画像解析システム》
　図３Ａは、本実施形態に係る画像処理装置２００を含むＳＡＲ画像解析システム３００の概要を示す図である。
　ＳＡＲ画像解析システム３００は、画像処理装置２００を含むＳＡＲ画像解析装置３１０とＳＡＲ衛星３５０とを含む。ＳＡＲ画像解析システム３００は、ＳＡＲ画像を処理することにより、地上の構造物、例えば、高層ビルや高層マンションなどの傾きや地盤沈下などによるビルの沈み込みなどを判定する。
　ＳＡＲ衛星３５０は、ＳＡＲ画像を撮影する。ＳＡＲ衛星３５０は、地表や、地上にある建造物、例えば、高層ビルや高層マンションなどのＳＡＲ画像を撮影する。ＳＡＲ衛星３５０が撮影したＳＡＲ画像は、例えば、ＳＡＲ画像記憶部２０４に保存される。また、ＳＡＲ衛星３５０は、撮影したＳＡＲ画像をＳＡＲ画像解析装置３１０に送信してもよい。
　ＳＡＲ画像解析装置３１０は、ＳＡＲ衛星３５０が撮影した複数のＳＡＲ画像を取得する。ここで、ＳＡＲ画像解析装置３１０が取得する複数のＳＡＲ画像は、例えば、地表の同一地点や同一領域、同一範囲などを異なる撮影条件（日時、天候、撮影に用いたＳＡＲ衛星など）で撮影したＳＡＲ画像である。ＳＡＲ画像解析装置３１０は、ＳＡＲ画像記憶部２０４に保存されたＳＡＲ画像を取得してもよい。また、ＳＡＲ画像解析装置３１０は、ＳＡＲ衛星３５０から直接ＳＡＲ画像を取得してもよい。そして、ＳＡＲ画像解析装置３１０は、取得したＳＡＲ画像から安定反射点を特定する。ＳＡＲ画像解析装置３１０は、さらに、特定した安定反射点の各々における位相を特定する。
　図３Ｂは、本実施形態に係る画像処理装置２００を含むＳＡＲ画像解析システム３００の機能構成を示す図である。
　図３Ｂに示したように、ＳＡＲ画像解析システム３００は、ＳＡＲ画像記憶部２０４、ＳＡＲ画像解析装置３１０、安定反射点記憶部３３０および変位・標高記憶部３４０を有する。ここで、ＳＡＲ画像解析装置３１０は、本実施形態の画像処理装置２００、安定反射点特定部２０６および変位・標高算出部３１１を有する。
　そして、ＳＡＲ画像解析装置３１０は、次のように動作する。本実施形態の画像処理装置２００は、複数のＳＡＲ画像間にわたる同一地点（画素）で観測された位相の違いを解析するために、同じ場所を異なる時期に撮影した複数のＳＡＲ画像をＳＡＲ画像記憶部２０４から取得する。本実施形態の画像処理装置２００は、クラスタリングされたクラスタの代表的な位相変化と各画素の位相配列との相関を、安定反射点らしさの評価指標として生成する。安定反射点特定部２０６は、取得した評価指標に基づいて安定反射点を特定する。特定された安定反射点は、例えば、安定反射点記憶部３３０に記憶される。
　次に、変位・標高算出部３１１は、安定反射点特定部２０６で特定された安定反射点を入力とし、各安定反射点における変位および標高を算出する。算出した変位および標高は、例えば、変位・標高記憶部３４０に記憶される。
　（位相相関算出部）
　図４は、本実施形態に係る位相相関算出部２０３の機能構成を示すブロック図である。
　図４を参照すると、位相相関算出部２０３は、位相統計量演算部４３１と、位相相関演算部４３２とを有する。
　（位相統計量演算部）
　位相統計量演算部４３１は、各クラスタに対して、各クラスタの代表的な位相変化を表す統計量を特定する。各クラスタの代表的な位相変化を表す統計量としては、例えば、各クラスタを構成する画素各々における位相配列の平均を用いる。または、距離に応じて変化する荷重関数を用いて位相配列の荷重平均を、各クラスタの代表的な位相変化として特定してもよい。すなわち、仮に相関を特定する対象の画素をi、クラスタをCとし、w(d)を例えば正の入力dに対して単調非増加な正の関数とし、画素各々の位置ベクトルをｒi、位相配列の各要素を格納した複素の列ベクトルをｓiとすると、場所rにおけるクラスタCの代表的な位相変化を表すベクトルｓC(r)は、下記の数式（１）のように特定される。
Figure JPOXMLDOC01-appb-M000001
　また、各クラスタの代表的な位相変化を表す統計量として、クラスタを構成する画素各々における位相配列に対して、その位相配列の要素の全てのペアに対する位相差を要素として持つ行列を算出し、これら行列の平均値を算出して用いてもよい。例えば、位相配列の要素の全ペアに対する位相差を記録した行列は、仮に画素iに対して求めるならば、画素iにおける位相配列の要素を収めた列ベクトルｓiを用いてｓi ｓi *として求めることができる。ただし、*は複素共役転置を表す記号である。各クラスタの代表的な位相変化を表す統計量としては、仮にクラスタCに対して特定するのであれば、クラスタCを構成する各画素iに対して求められた行列ｓi ｓi *の平均値を用いればよい。
　また、数式（１）において定義したのと同様に、行列ｓi ｓi *の加重平均を、各クラスタの代表的な位相変化を表す統計量であるとしてもよい。すなわち、ｓC(r)を定義する際に用いたのと同様の文字の定義を用いて、クラスタの代表的な位相変化を表す統計量は、行列の加重平均ＳC(r)の形式では、下記の数式（２）のように特定される。
Figure JPOXMLDOC01-appb-M000002
　また、クラスタの代表的な位相変化の特定に際して、ベクトルや行列の平均値ではなく、最頻値や中央値等、他の統計量を取っても構わない。また、上では複素数のベクトルや行列に対して複素数の統計量を直接計算したが、一旦、偏角を取り出し、偏角に対する統計量を特定してから、複素数に戻してもよい。また、平均計算では一般的に、求まるベクトルや行列の各要素の絶対値は“１”にならないが、絶対値が“１”になるように正規化し直しても構わない。
　（位相相関演算部）
　位相相関演算部４３２は、クラスタの代表的な位相変化と、画像全体の画素各々における位相変化との相関を演算する。この相関が強い画素ほど、ノイズが少なく、より安定反射点らしいと評価できる。相関を演算するにあたっては、例えば、画像全体の画素各々に対して下記のような計算を行う。仮に、相関を特定する画素のインデックスをiとし、画素iにおける位相配列の要素を収めた列ベクトルをｓiとし、画素iの位置ベクトルをｒiとする。クラスタCの代表的な位相変化を表す統計量として、数式（１）で定義したベクトルｓC(r)を用いるとする。このとき、位相相関演算部４３２は、クラスタCの代表的な位相変化と画素iにおける位相変化の相関corr(i,C)を、下記の数式（３）によって、画像全体の全ての画素iの各々に対して特定する。
Figure JPOXMLDOC01-appb-M000003
　または、クラスタＣの代表的な位相変化を表す統計量として、数式（２）で定義した行列ｓC(r)を用いるとすれば、位相相関演算部４３２はクラスタCの代表的な位相変化と画素iにおける位相変化の相関corr(i,C)を、下記の数式（４）によって、全ての画素iの各々に対して特定する。ただし、v_oneは全ての要素が1であるような列ベクトルである。
Figure JPOXMLDOC01-appb-M000004
　また、相関としては、数式（４）のようにして特定された値に、平方根や二乗等、正の入力に対して単調非減少な関数による変換を適宜適用したものを用いても構わない。
　以上のようにして算出された、各クラスタの代表的な位相変化と、画像全体の画素の各々における位相変化との相関は、画素の各々におけるノイズの量が少ないほど強くなるため、ノイズの少ない画素である安定反射点を特定するのに役に立つ。すなわち、相関が高い画素ほど安定反射点らしいと判断することができる。ただし、上述した各クラスタの代表的な位相変化と、画像全体の画素の各々における位相変化との相関は、画像全体の画素の各々に対してクラスタの数だけ求まる。
　そこで、これらの相関を画像全体の画素の各々が安定反射点かどうかを特定するのに用いるために、位相相関演算部４３２は、画素各々に対して、複数のクラスタとの間で特定した相関を１つにまとめ、その値を画素各々における相関として出力する。このまとめる処理については、例えば、画像全体の画素の各々に対してそれを含むクラスタを特定しておき、そのクラスタの代表的な位相変化との相関のみを用いるようにすることで実現される。または、画像全体内の画素の各々に対して、各クラスタとの間で特定された複数の相関の中から、最大の相関を取り出して用いてもよい。あるいは、画像全体の画素の各々に対して最近接するクラスタを特定しておき、そのクラスタの代表的な位相変化との相関のみを用いるようにしてもよい。
　《各データの記憶構成》
　以下、図５Ａ～図５Ｇを参照して、本実施形態で使用される各データの記憶構成を説明する。なお、以下では便宜上、図５Ａから図５Ｇに分けて説明するが、複数の記憶構成が一体であっても、さらに分けられていてもよい。
　（ＳＡＲ画像記憶部）
　図５Ａは、本実施形態に係るＳＡＲ画像記憶部２０４の構成を示す図である。ＳＡＲ画像記憶部２０４は、衛星などで撮像して送信されたＳＡＲ画像を蓄積する。なお、ＳＡＲ画像記憶部２０４の構成は、図５Ａに限定されない。
　ＳＡＲ画像記憶部２０４は、各ＳＡＲ画像ＩＤ５１１のＳＡＲ画像に含まれる反射点ＩＤ（画像ＩＤ）５１２を記憶する。各反射点ＩＤ（画像ＩＤ）５１２に対応つけて、輝度５１３、位置５１４、位相５１５などを記憶する。
　（位相配列蓄積部）
　図５Ｂは、本実施形態に係る位相配列蓄積部２０１の構成を示す図である。位相配列蓄積部２０１は、ＳＡＲ画像に含まれる各反射点（画素）の複数のＳＡＲ画像に渡る位相配列を蓄積する。なお、位相配列蓄積部２０１の構成は、図５Ｂに限定されない。また、位相配列蓄積部２０１は、ＳＡＲ画像記憶部２０４から各反射点（画素）の複数のＳＡＲ画像に渡る位相配列を呼び出す位相配列生成部として実現してもよい。
　位相配列蓄積部２０１は、各反射点ＩＤ（画素ＩＤ）５２１に対応して、各ＳＡＲ画像のタイムスタンプと位相との組からなる位相配列５２２を記憶する。
　（安定反射点候補テーブル）
　図５Ｃは、本実施形態に係る安定反射点候補テーブル５３０の構成を示す図である。安定反射点候補テーブル５３０は、安定反射点候補特定部２０５が複数のＳＡＲ画像における各反射点（画素）の中で、安定した反射点の候補を特定するために使用する。なお、安定反射点候補テーブル５３０の構成は、図５Ｃに限定されない。
　安定反射点候補テーブル５３０は、各反射点ＩＤ（画素ＩＤ）５３１に対応付けて、安定反射点候補を特定するための安定反射点候補条件５３２と、安定反射点候補条件５３２による判定結果の安定反射点候補フラグ５３３とを記憶する。ここで、安定反射点候補フラグ５３３は、“１”の場合が安定反射点候補であり、“０”の場合が安定反射点候補でないとする。また、安定反射点候補条件５３２としては、輝度の安定、位置の安定、位相の安定などがあるが、これらに限定されない少なくとも１つの条件により判定する。
　（クラスタテーブル）
　図５Ｄは、本実施形態に係るクラスタテーブル５４０の構成を示す図である。クラスタテーブル５４０は、クラスタリング部２０２が位相配列に基づいて反射点（画素）をクラスタリングするために使用する。なお、クラスタテーブル５４０の構成は、図５Ｄに限定されない。
　クラスタテーブル５４０は、クラスタＩＤ５４１に対応して、そのクラスタに分類された反射点ＩＤ（画素ＩＤ）５４２と、その位相配列５４３と、安定反射点候補フラグ５４４と、を記憶する。ここで、安定反射点候補のみをクラスタリングする場合には、安定反射点候補フラグ５４４は全て“１”の反射点（画素）がクラスタリングされる。
　（位相統計量テーブル）
　図５Ｅは、本実施形態に係る位相統計量テーブル５５０の構成を示す図である。位相統計量テーブル５５０は、位相統計量演算部４３１が各クラスタの代表的な位相変化を演算するために使用される。なお、位相統計量テーブル５５０の構成は、図５Ｅに限定されない。
　位相統計量テーブル５５０は、クラスタＩＤ５５１に対応して、各クラスタに分類された反射点ＩＤ（画素ＩＤ）５５２と、その位相配列５５３とを記憶する。そして、各クラスタに分類された複数の反射点（画素）に基づいて算出された、各クラスタの位相統計量５５４を記憶する。位相統計量５５４には、例えば、位相配列の平均や位相配列の加重平均などが位相変化として記憶される。
　（位相相関テーブル）
　図５Ｆは、本実施形態に係る位相相関テーブル５６０の構成を示す図である。位相相関テーブル５６０は、位相相関演算部４３２が各反射点（画素）の位相配列と各クラスタの位相統計量（位相変化）との相関を演算するために使用される。なお、位相相関テーブル５６０の構成は、図５Ｆに限定されない。
　位相相関テーブル５６０は、各反射点（画素）の位相配列群５６１と、各クラスタの位相統計量（位相変化）群５６２との間の、全組合せの位相相関値が算出され記憶される。
　（安定反射点テーブル）
　図５Ｇは、本実施形態に係る安定反射点テーブル５７０の構成を示す図である。安定反射点テーブル５７０は、安定反射点特定部２０６が位相相関算出部２０３からの位相相関値に基づいて安定反射点を特定するために使用される。なお、安定反射点テーブル５７０の構成は、図５Ｇに限定されない。
　安定反射点テーブル５７０は、各反射点ＩＤ（画素ＩＤ）５７１に対応付けて、安定反射点候補フラグ５７２と、位相相関算出結果５７３と、安定反射点フラグ５７４と、を記憶する。ここで、安定反射点フラグ５７４は、“１”の場合が安定反射点であり、“０”の場合が安定反射点でないとする。
　《画像処理装置のハードウェア構成》
　図６は、本実施形態に係る画像処理装置２００のハードウェア構成を示すブロック図である。なお、図６においては、安定反射点候補特定部２０５を有し、安定反射点特定部２０６を有しない構成を示すが、これに限定されない。
　図６で、ＣＰＵ(Central Processing Unit)６１０は演算制御用のプロセッサであり、プログラムを実行することで図２および図４の機能構成部を実現する。なお、ＣＰＵ６１０は、それぞれの機能に対応して複数あってもよい。ＲＯＭ(Read Only Memory)６２０は、初期データおよびプログラムなどの固定データおよびプログラムを記憶する。ネットワークインタフェース６３０は、ネットワークを介して、ＳＡＲ画像記憶部２０４などの他の装置との通信を制御する。
　ＲＡＭ(Random Access Memory)６４０は、ＣＰＵ６１０が一時記憶のワークエリアとして使用するランダムアクセスメモリである。ＲＡＭ６４０には、本実施形態の実現に必要なデータを記憶する領域が確保されている。ＳＡＲ画像データ６４１は、ＳＡＲ画像記憶部２０４から取得したＳＡＲ画像のデータである。位相配列データ６４２は、位相配列蓄積部２０１から取得した各反射点（画素）の位相配列のデータである。安定反射点候補データ６４３は、安定反射点候補テーブル５３０から取得した各反射点（画素）の安定反射点候補のデータである。クラスタデータ６４４は、クラスタテーブル５４０から取得した各反射点（画素）が分類されたクラスタのデータである。位相統計量算出データ６４５は、位相統計量テーブル５５０から取得した各クラスタの位相統計量のデータである。位相相関算出データ６４６は、位相相関テーブル５６０から取得した各反射点（画素）の位相配列と各クラスタの位相統計量（位相変化）との相関のデータである。入出力データ６４７は、入出力インタフェース６６０を介して、入出力機器と入出力するデータである。送受信データ６４８は、ネットワークインタフェース６３０を介して、他の装置と送受信を行うデータである。
　ストレージ６５０は、データベースや各種のパラメータ、あるいは本実施形態の実現に必要な以下のデータまたはプログラムが記憶されている。位相配列蓄積部２０１は、図５Ｂに示したように、各反射点（画素）の位相配列を蓄積する。安定反射点候補テーブル５３０は、図５Ｃに示したように、各反射点（画素）が安定反射点候補か否かを記憶する。クラスタテーブル５４０は、図５Ｄに示したように、各クラスタに分類された反射点（画素）を記憶する。位相統計量テーブル５５０は、図５Ｅに示したように、各クラスタの位相統計量（位相変化）を記憶する。位相相関テーブル５６０は、図５Ｆに示したように、各反射点（画素）の位相配列と各クラスタの位相統計量（位相変化）との相関を記憶する。
　ストレージ６５０には、以下のプログラムが格納される。画像処理装置制御プログラム６５１は、本実施形態の画像処理装置２００の全体の処理を制御するプログラムである。位相配列蓄積モジュール６５２は、複数のＳＡＲ画像に渡る各反射点（画素）の位相配列を蓄積するモジュールである。安定反射点候補特定モジュール６５３は、各反射点（画素）のデータの安定らしさにより安定反射点候補を特定するモジュールである。クラスタリングモジュール６５４は、各反射点（画素）の位相配列に基づいて、各反射点（画素）を複数のクラスタにクラスタリングするモジュールである。位相統計量演算モジュール６５５は、各クラスタの代表的な位相特徴量（位相変化）を演算するモジュールである。位相相関演算モジュール６５６は、各反射点（画素）の位相配列と各クラスタの位相統計量（位相変化）との相関を演算するモジュールである。
　入出力インタフェース６６０は、入出力デバイスとのデータ入出力を制御するためのインタフェースを行なう。本実施形態においては、入出力インタフェース６６０には、画像処理装置２００の状態を表示する表示部６６１や画像処理装置２００への操作入力をする操作部６６２が接続される。さらに、ＳＡＲ画像を格納した記憶媒体６６３が接続されてもよい。
　なお、図６のＲＡＭ６４０やストレージ６５０には、画像処理装置２００が有する汎用の機能や他の実現可能な機能に関連するプログラムやデータは図示されていない。
　《画像処理装置の処理手順》
　図７は、本実施形態に係る画像処理装置２００の処理手順を示すフローチャートである。このフローチャートは、図６のＣＰＵ６１０がＲＡＭ６４０を用いて実行し、図２および図４の機能構成部を実現する。
　画像処理装置２００は、ステップＳ７０１において、複数のＳＡＲ画像から各反射点（画素）の位相配列を生成して蓄積する。画像処理装置２００は、ステップＳ７０３において、ＳＡＲ画像の各反射点（画素）を、位相配列に基づいて複数のクラスタにクラスタリングする。画像処理装置２００は、ステップＳ７０５において、各クラスタ内の反射点（画素）の位相配列から各クラスタの代表的な位相変化を表す統計量を算出する。画像処理装置２００は、ステップＳ７０７において、各クラスタの代表的な位相変化を表す統計量と、画像全体の各反射点（画素）の位相配列との間での相関を算出する。画像処理装置２００は、ステップＳ７０９において、算出された相関をそのまま、あるいは、まとめて出力する。相関が強い反射点（画素）ほど、ノイズの影響が少なく、変位解析に適した安定反射点らしいであると判断することができ、安定した評価指標となる。
　本実施形態によれば、複数のクラスタの各クラスタ内での位相変化と各画素における位相配列との相関を安定反射点の評価に使用するので、変位が非線形である対象や、標高が高い対象、変位が大きな対象について、大きな位相変化に影響されずに安定反射点らしさを評価できる。
　すなわち、本実施形態によれば、非線形な変位を起こす対象や、標高が高い対象、変位が大きな対象に対して、変位解析を行うのに必要な安定反射点の特定を行うことができる。というのは、画像全体の各々の画素に対して時間に線形に起こる変位等の仮定を行う必要が無く、モデルの当て嵌めも行う必要がないためである。本実施形態では、一旦、画素を類似した位相変化を持つクラスタに分割し、そのクラスタ各々における代表的な位相変化についての統計量を算出する。これにより、ノイズによる影響を取り除いた位相変化がクラスタの各々に対して求まる。画像全体の画素の各々に対して、クラスタ各々の代表的な位相変化との相関を算出することにより、画像全体の画素の各々におけるノイズが含まれた位相変化と、クラスタの各々におけるノイズが含まれていない位相変化がどの程度相関するかが分かる。そのため、相関が強い画素はノイズの少ない安定反射点であると特定することができる。
　［第３実施形態］
　次に、本発明の第３実施形態に係る画像処理装置について説明する。本実施形態に係る画像処理装置は、上記第２実施形態と比べると、算出された位相統計量からノイズを除去するノイズ除去部を有する点で異なる。その他の構成および動作は、第２実施形態と同様であるため、同じ構成および動作については同じ符号を付してその詳しい説明を省略する。
　（位相相関算出部）
　図８は、本実施形態に係る画像処理装置８００の位相相関算出部８０３の機能構成を示すブロック図である。なお、画像処理装置８００の構成は、図２の位相相関算出部２０３が位相相関算出部８０３に置き換えられたものである。
　本実施形態における位相相関算出部８０３は、位相統計量演算部４３１と、ノイズ除去部８３３と、位相相関演算部４３２とを有する。すなわち、本実施形態における位相相関算出部８０３は、位相統計量演算部４３１と位相相関演算部４３２との間に、ノイズ除去部８３３を有する点のみが、第２実施形態の位相相関算出部２０３と異なる。
　位相統計量演算部４３１は、クラスタの代表的な位相変化を表す統計量を算出する。ノイズ除去部８３３は、位相統計量演算部４３１で算出された統計量から、ノイズの影響を除去する。位相相関演算部４３２は、ノイズ除去部８３３によりノイズの影響が除去された統計量を用いて、各反射点（画素）との相関を算出する。
　（ノイズ除去部）
　続いて、ノイズ除去部８３３について詳細に説明する。ノイズ除去部８３３は、位相統計量演算部４３１で算出された、クラスタの代表的な位相変化に関する統計量からノイズの影響を除去する。位相統計量演算部４３１において、クラスタ各々の代表的な位相変化に関する統計量を算出した際に、ある程度のノイズの影響は除去されている。本実施形態においては、さらに、各クラスタ内の画素の各々において特定されうる位相配列についての事前知識を用い、事前知識に合致しないものをノイズとして除去することにより、ノイズに対する頑健さを増すことができる。
　ノイズの影響を除去するにあたっての事前知識としては、例えば、各クラスタ内の画素各々における位相変化が、各クラスタを代表する位相変化のただ１種に対して、ノイズが加わったものであるという知識を用いる。この事前知識によれば、数式（２）によって算出された統計量であるところの行列ＳC(r)は、固有値分解によってクラスタの代表的な位相変化に対応した固有ベクトル１本と、ノイズに対応する残りの固有ベクトルとに分かれる性質を持つことが分かる。したがって、例えば、数式（２）によって算出された特徴量であるところの行列ＳC(r)を固有値分解し、最大固有ベクトルをクラスタの代表的な位相変化に対応する固有ベクトルであるとし、それ以外の固有ベクトルをノイズに起因するものであるとする。このように、行列ＳC(r)から、最大固有ベクトルに対応した行列と、それ以外の行列とに分離することによって、ノイズの影響が除去された統計量を得ることができる。
　また、クラスタの代表的な位相変化と、各クラスタの画素各々における位相変化の間に、ノイズに起因するものではない変化が少量含まれているという事前知識を用いてもよい。この事前知識によれば、この事前知識によれば、数式（２）によって特定された統計量であるところの行列ＳC(r)は、固有値分解によってクラスタの代表的な位相変化に対応した固有ベクトル数本と、ノイズに対応する残りの固有ベクトルとに分かれる性質を持つことが分かる。したがって、例えば、数式（２）によって算出された特徴量であるところの行列ＳC(r)を固有値分解し、固有値の大きい固有ベクトル数本をクラスタの代表的な位相変化に対応する固有ベクトルであるとし、それ以外の固有ベクトルをノイズに起因するものであるとする。このように、行列ＳC(r)から、固有値の大きい固有ベクトル数本に対応した行列と、それ以外の行列とに分離することによって、ノイズの影響が除去された統計量を得ることができる。
　また、変位について時間に対して緩やかであるといったモデルを仮定して、ノイズの影響を除去してもよい。この場合には、数式（２）によって算出された統計量であるところの行列ＳC(r)上では、画像取得時期としての撮影時期の近い画像に対応する位相配列の要素同士で位相比較を行っている成分については絶対値が大きくなる。一方、撮影時期の離れた画像に対応する位相配列の要素同士で位相比較を行っている成分については絶対値が小さくなる。ただし、ノイズに関しては、行列ＳC(r)の対角成分を除く成分については、時期が離れているかどうか関係なく絶対値を大きくするような性質がある。以上のことから、撮影時期の近いもの同士での位相比較に対応する行列ＳC(r)中の成分はそのまま、撮影時期の離れたもの同士での位相比較に対応した行列ＳC(r)中の成分は強制的に“０（ゼロ）”に近づけることで、ノイズの影響を除去することができる。
　（ノイズ除去テーブル）
　図９は、本実施形態に係るノイズ除去テーブル９１０の構成を示す図である。ノイズ除去テーブル９１０は、ノイズ除去部８３３が位相統計量演算部４３１により算出された位相統計量（位相変化）のノイズを除去するために使用される。なお、ノイズ除去テーブル９１０の構成は、図９に限定されない。
　ノイズ除去テーブル９１０は、クラスＩＤ９１１とその位相統計量９１２に対応して、複数のノイズ除去方法９１３と、各ノイズ除去方法９１３によりノイズ除去された位相統計量９１４と、を記憶する。ここで、ノイズ除去方法９１３としては、固有値分解による最大固有ベクトルを選択、固有値分解による大きい数本の固有ベクトルを選択、撮影時期の離れたもの同士の位相比較を排除、などが含まれる。
　《画像処理装置の処理手順》
　図１０は、本実施形態に係る画像処理装置８００の処理手順を示すフローチャートである。このフローチャートは、ノイズ除去テーブル９１０が追加された図６のＣＰＵ６１０がＲＡＭ６４０を用いて実行し、図２および図８の機能構成部を実現する。なお、図１０において、図７と同様のステップには同じステップ番号を付して、重複する説明を省略する。
　ステップＳ２０７からステップＳ７０５までは、図７と同様に動作する。その後、画像処理装置８００は、ステップＳ１００６において、ステップＳ７０５で算出されたクラスタ内の統計量からノイズの影響を除去する。画像処理装置８００は、ステップＳ７０７では、ノイズの影響を除去された統計量を用いて、図７と同様の方法で、各クラスタの代表的な位相変化と、評価対象の位相変化がどの程度合致するかという相関を算出する。したがって、相関が強い画素ほど、ノイズの影響が少なく、変位解析に適した安定反射点らしい判断することができる。
　本実施形態によれば、非線形な変位を起こす対象や、標高が高い対象、変位が大きな対象に対して、変位解析を行うのに必要な安定反射点の特定を行うことができる。特に、本実施形態には、構成としてノイズの影響を除去する部分が含まれており、より頑健に安定反射点を特定することができる。
　［第４実施形態］
　次に、本発明の第４実施形態に係る画像処理装置について説明する。本実施形態に係る画像処理装置は、上記第２実施形態と比べると、位相相関算出部の相関演算部と位相統計量演算部とが逆に配置される点で異なる。その他の構成および動作は、第２実施形態と同様であるため、同じ構成および動作については同じ符号を付してその詳しい説明を省略する。
　（位相相関算出部）
　図１１は、本実施形態に係る画像処理装置１１００の位相相関算出部１１０３の機能構成を示すブロック図である。なお、画像処理装置１１００の構成は、図２の位相相関算出部２０３が位相相関算出部１１０３に置き換えられたものである。
　本実施形態における位相相関算出部１１０３は、位相相関演算部１１３１と、位相統計量演算部１１３２とを有する。位相相関演算部１１３１は、画像全体の画素の各々において位相配列蓄積部２０１に蓄積された位相配列と、クラスタリング部２０２でクラスタリングされた各クラスタを構成する各反射点（画素）における位相配列との相関を演算する。位相統計量演算部１１３２は、画像全体の各反射点（画素）に対して算出された、各クラスタを構成する各反射点（画素）における位相変化との相関に対して、例えば、クラスタ内における相関の平均値等といった統計量を演算する。
　位相相関算出部１１０３においては、評価対象の各々と、各クラスタを構成する画素の各々との、位相変化の相関の統計量を算出する。
　（位相相関演算部）
　位相相関演算部１１３１は、画像全体の画素の各々と、クラスタを構成する画素の各々との位相変化の相関を算出する。画素の位相変化同士間の相関は、例えば次のように演算する。仮に、相関を特定する対象の画素のインデックスをiとし、クラスタをCとし、クラスタC内の画素のインデックスをj∈Cとする。画素iにおける位相配列の要素を列ベクトルｓiに、画素jにおいて観測された位相配列の要素を列ベクトルｓjに格納するとすれば、画素iと画素jの位相変化同士の相関corr(i,j)は、下記の数式（５）で算出される。
Figure JPOXMLDOC01-appb-M000005
　また、上記の画素iとjの位相変化同士の相関に対して、正の数に対して単調非減少な関数を適用しても構わない。例えば二乗しても構わない。または、ある閾値以下の入力に対しては“０（ゼロ）”となるような関数を適用してもよい。このような関数の適用により、極端に相関が弱い画素の影響を除去して評価を行うことが可能になる。
　（位相統計量演算部）
　位相統計量演算部１１３２は、位相相関演算部１１３１において算出された、画像全体の各画素の位相配列と、各クラスタを構成する各画素における位相配列との相関を入力とする。そして、画像全体の各画素における位相配列と各クラスタを構成する画素における全体としての位相変化との相関を表す統計量を算出する。仮に、相関を演算する対象の画素のi、相関を演算する対象のクラスタをCとすると、位相統計量演算部１１３２は、クラスタ全体との相関を表す統計量corr(i,C)を、数式（５）で求まったcorr(i,j)に対して、クラスタCに含まれる全ての画素j∈Cに関する統計量として算出する。例えば、下記の数式（６）のように、corr(i,j)のj∈Cに対する平均値を演算する。
Figure JPOXMLDOC01-appb-M000006
　あるいは、最大値、最頻値、中央値等の他の統計量を用いてもよい。また、例えば、距離に応じて単調非増加な荷重関数を用いて重み付けした評価を行ってもよい。例えば、荷重平均を用いるならば、w(d)をdに対して単調非増加な関数とし、ｒi,ｒjをそれぞれ画素i,jの位置として、下記の数式（７）のように特定してもよい。
Figure JPOXMLDOC01-appb-M000007
　以上のようにして算出された、画像全体の画素の各々における位相変化と、各クラスタについて構成する画素における位相変化全体としての相関を表す統計量は、第２実施形態における位相相関算出部２０３で特定した画素の各々における位相変化とクラスタ各々の代表的な位相変化との相関と、ほとんど同様の扱いをすることができる。すなわち、画像全体の画素の各々における位相変化と、各クラスタについて構成する画素における位相変化全体としての相関を表す統計量が、強くなるほどより安定反射点らしいと判断することができる。本実施形態における、画素各々に対するクラスタ内全体での相関の統計量は、クラスタごとに求まるため、各画素に対して複数求まるが、第２実施形態と同様にして、画素ごとに１つの値にまとめられる。
　（位相相関テーブル）
　図１２Ａは、本実施形態に係る位相相関テーブル１２１０の構成を示す図である。位相相関テーブル１２１０は、位相相関演算部１１３１が画像全体の各画素の位相配列とクラスタを構成する各画素の位相配列との相関を算出するために使用される。なお、位相相関テーブル１２１０の構成は、図１２Ａに限定されない。
　位相相関テーブル１２１０は、各反射点（画素）の位相配列群５６１と、各クラスタに分類された各反射点（画素）の位相配列群１２６２との間の、全組合せの位相相関値が算出され記憶される。
　（位相統計量テーブル）
　図１２Ｂは、本実施形態に係る位相統計量テーブル１２２０の構成を示す図である。位相統計量テーブル１２２０は、位相統計量演算部１１３２が位相相関演算部１１３１において算出された各画素の位相配列と各クラスタを構成する各画素の位相配列との相関から、各画素における位相配列と各クラスタの位相変化との相関を表す統計量を算出するために使用される。なお、位相統計量テーブル１２２０の構成は、図１２Ｂに限定されない。
　位相統計量テーブル１２２０は、反射点ＩＤ（画素ＩＤ）１２２１に対応して、各クラスタの相関統計量１２２２、１２２３、…を記憶する。各クラスタの相関統計量１２２２、１２２３、…は、各反射点ＩＤ（画素ＩＤ）１２２１と各クラスタ内の各反射点（画素）との相関群と、その相関群から算出された相関統計量として相関平均や相関加重平均を記憶する。
　《画像処理装置の処理手順》
　図１３は、本実施形態に係る画像処理装置１１００の処理手順を示すフローチャートである。このフローチャートは、位相相関テーブル１２１０や位相統計量テーブル１２２０を有する図６のＣＰＵ６１０がＲＡＭ６４０を用いて実行し、図２および図１１の機能構成部を実現する。なお、図１３において、図７と同様のステップには同じステップ番号を付して、重複する説明を省略する。ステップＳ７０１とＳ７０３は、図７と同様に動作する。
　画像処理装置１１００は、ステップＳ１３０５において、画像全体の各画素における位相配列と、各クラスタを構成する各画素における位相配列との相関を算出する。次に、画像処理装置１１００は、ステップＳ１３０７においては、ステップＳ１３０５で算出された、各画素における位相配列と各クラスタを構成する各画素における位相配列との相関に対して、クラスタごとに相関統計量を算出する。画像処理装置１１００は、ステップＳ１３０９において、算出された相関統計量をそのまま、あるいは、まとめて出力する。最終的に算出された、クラスタ内の相関についての統計量が強いほど、ノイズの影響が少なく、変位解析に適した安定反射点であると判断することができる。
　本実施形態によれば、非線形な変位を起こす対象や、標高が高い対象、変位が大きな対象に対して、変位解析を行うのに必要な安定反射点の特定を行うことができる。特に、本実施形態では、画像全体の画素の各々における位相変化と、各クラスタを構成する画素の各々における位相変化との相関を算出するため、画像全体の画素の各々に対して、各クラスタを構成する画素の中から、極端に異なる位相変化を持つ画素の影響を除去して評価を行うことができる。このような評価方法は、例えば、クラスタリングのアルゴリズムによって、クラスタが十分に分かれず、複数のクラスタが混ざったような状態で特定されるような場合に有用である。
　なお、本実施形態において、統計量演算部１１３２で、平均以外の、例えば、最頻値を統計量として算出するなどによって、さらにノイズの影響を除去する効果を得ることができる。
　［第５実施形態］
　次に、本発明の第５実施形態に係る画像処理装置について説明する。本実施形態に係る画像処理装置は、上記第２実施形態乃至第４実施形態と比べると、位相相関算出部の算出結果を安定反射点候補特定部にフィードバックする点で異なる。その他の構成および動作は、第２実施形態と同様であるため、同じ構成および動作については同じ符号を付してその詳しい説明を省略する。
　《画像処理装置の機能構成》
　図１４は、本実施形態に係る画像処理装置１４００の機能構成を示すブロック図である。なお、図１４において、図２と同様の機能構成部には同じ参照番号を付して、重複する説明は省略する。図１４においては、図２の位相相関算出部２０３および安定反射点候補特定部２０５を、位相相関算出部１４０３および安定反射点候補特定部１４０５に置き換え、その間をフィードバック線で結んだものである。
　安定反射点候補特定部１４０５は、安定反射点になりうる画素を特定する。特に、本実施形態における安定反射点候補特定部１４０５は、第２実施形態における安定反射点候補特定部２０５と異なり、位相相関算出部１４０３で算出された相関を元に、安定反射点候補を再度、特定しなおす機能を有する。クラスタリング部２０２は、安定反射点候補特定部１４０５で再度、特定された安定反射候補を入力として、再度、クラスタリングを行う。位相相関算出部１４０３は、画像全体の各画素における位相配列に対して、クラスタリング部２０２で再度、クラスタリングされたクラスタ各々を代表する位相変化との相関を算出する。このように、位相相関算出部１４０３で算出された相関は安定反射点候補特定部１４０５にフィードバックされ、再度、安定反射点候補の特定が行われる。
　なお、位相相関算出部１４０３としては、前述した第２実施形態における位相相関算出部２０３、あるいは、第３実施形態における位相相関算出部８０３、あるいは、第４実施形態における位相相関算出部１１０３、のいずれを用いてもよい。
　（安定反射点候補テーブル）
　図１５は、本実施形態に係る安定反射点候補テーブル１５３０の構成を示す図である。安定反射点候補テーブル１５３０は、安定反射点候補特定部１４０５が複数のＳＡＲ画像における各反射点（画素）の中で、安定した反射点の候補の特定を、相関を基に繰り返すために使用する。なお、図１５において、図５Ｃと同様の構成要素には同じ参照番号を付して、重複する説明は省略する。また、安定反射点候補テーブル１５３０の構成は、図１５に限定されない。
　安定反射点候補テーブル１５３０は、各反射点ＩＤ（画素ＩＤ）５３１に対応付けて、位相相関算出部１４０３からフィードバックされた位相相関算出結果１５３２と、位相相関算出結果１５３２に基づいて変更される安定反射点候補フラグ１５３３と、を記憶する。ここで、安定反射点候補フラグ１５３３は、“１”の場合が安定反射点候補であり、“０”の場合が安定反射点候補でないとする。例えば、反射点ＩＤ（画素ＩＤ）５３１が“A0001-0002”の反射点（画素）は、最初は安定反射点候補であったが位相相関算出結果１５３２に基づいて安定反射点候補から外されている。
　《画像処理装置の処理手順》
　図１６は、本実施形態に係る画像処理装置１４００の処理手順を示すフローチャートである。このフローチャートは、安定反射点候補テーブル１５３０を有する図６のＣＰＵ６１０がＲＡＭ６４０を用いて実行し、図１４の機能構成部を実現する。なお、図１４において、図７と同様のステップには同じステップ番号を付して、重複する説明を省略する。ステップＳ７０１は、図７と同様に動作する。
　画像処理装置１４００は、ステップＳ１６０３において、安定反射点となりうる画素を安定反射点候補として特定する。画像処理装置１４００は、ステップＳ１６０５において、安定反射点候補を位相配列に基づいてクラスタリングする。画像処理装置１４００は、ステップＳ１６０７において、各クラスタにおける代表的な位相変化と画像全体の各画素における位相配列との相関を算出する。なお、ステップＳ１６０７は、図７の（Ｓ７０５、Ｓ７０７）、図１０の（Ｓ７０５、Ｓ１００６、Ｓ７０７）、あるいは、図１３の（Ｓ１３０５、Ｓ１３０７）のいずれかを実行する。
　画像処理装置１４００は、ステップＳ１６０９において、特定されたクラスタ、あるいは、特定されたクラスタとの位相の相関が終了条件を満たしているかを判定する。終了条件としては、例えば、所定回数の繰り返しを行ったかどうかを用いることができる。あるいは、評価対象の画素のうち、安定反射点候補となった画素における相関が一定以上に強かった場合に、全ての安定反射点候補が安定反射点としての条件を満たす、つまり安定反射点候補特定部１４０５が安定反射点を正確に特定したとして、終了条件を満たすと判断してもよい。
　終了条件を満たさない場合、画像処理装置１４００は、ステップＳ１６１１において、現在の位相相関算出結果を安定反射点候補特定部１４０５にフィードバックして、ステップＳ１６０３から再度、再度安定反射点候補の特定を更新あるいはやり直して処理を繰り返す。終了条件を満たした場合、画像処理装置１４００は、ステップＳ１６１３において、絞られた安定反射点候補の位相相関算出結果を出力して処理を終了する。
　最終的に算出されたクラスタ内の相関についての統計量が強いほど、ノイズの影響が少なく、変位解析に適した安定反射点であると判断することができる。
　本実施形態によれば、非線形な変位を起こす対象や、標高が高い対象、変位が大きな対象に対して、変位解析を行うのに必要な安定反射点の特定を行うことができる。特に、本実施形態は、安定反射点候補の特定を、画像全体の画素の各々に対して特定された相関の強さに応じてやり直すことにより、安定反射点らしさをより正確に評価できるような相関を算出できる。というのは、ノイズが少ないと判断された安定反射点候補に対してクラスタを抽出し、その安定反射点候補の特定結果をより正確にしていくように構成したことによって、安定反射点候補のクラスタを代表する位相変化と、画像全体の画素の各々における位相変化との相関をする際のノイズの影響を低減することができるためである。
　［他の実施形態］
　本発明は、高層ビルや橋等非線形な変位や大きな変位を持ちうる対象の変位を、ＳＡＲ画像を用いて解析するようなプログラムに適用できる。また、音波を利用すれば、合成開口ソナーにも適用可能である。本発明によってクラスタ内で安定した位相を持つ画素を抽出すれば、反射強度だけではノイズに埋もれてしまっていた対象も見つけることが可能である。また、Time Of Flightカメラ等で撮影された位相画像にも適用できる。本発明を適用すれば、対象の位置が多少動くものに対しても、観測値の揺れがノイズに起因するものであるかを識別できるようになり、頑健に対象物の立体形状を抽出することができるようになる。また、コヒーレントな光の干渉を利用した位相画像にも適用することができる。
　さらに、撮像素子等のパターンと、撮影対象に付加されたパターンとの間での干渉による位相を利用した精密な歪み計測手法として、サンプリングモアレ法などがあるが、干渉結果を測定する際にノイズで信号が揺らいでいるのか、対象が動いているのかを識別しにくい。本手法により、位相が同期して変化しているかどうかを調べることにより、対象が動いているのかノイズが多いのかを分離することができるようになる。
　なお、実施形態を参照して本願発明を説明したが、本願発明は上記実施形態に限定されるものではない。本願発明の構成や詳細には、本願発明のスコープ内で当業者が理解し得る様々な変更をすることができる。また、それぞれの実施形態に含まれる別々の特徴を如何様に組み合わせたシステムまたは装置も、本発明の範疇に含まれる。
　また、本発明は、複数の機器から構成されるシステムに適用されてもよいし、単体の装置に適用されてもよい。さらに、本発明は、実施形態の機能を実現する画像処理プログラムが、システムあるいは装置に直接あるいは遠隔から供給される場合にも適用可能である。したがって、本発明の機能をコンピュータで実現するために、コンピュータにインストールされるプログラム、あるいはそのプログラムを格納した媒体、そのプログラムをダウンロードさせるＷＷＷ(World Wide Web)サーバも、本発明の範疇に含まれる。特に、少なくとも、上述した実施形態に含まれる処理ステップをコンピュータに実行させるプログラムを格納した非一時的コンピュータ可読媒体（non-transitory computer readable medium）は本発明の範疇に含まれる。
　［実施形態の他の表現］
　上記の実施形態の一部または全部は、以下の付記のようにも記載されうるが、以下には限られない。
（付記１）
　複数の画像に渡る各画素の位相配列を蓄積する位相配列蓄積手段と、
　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリング手段と、
　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出手段と、
　を備える画像処理装置。
（付記２）
　前記位相相関算出手段は、
　　各クラスタ内での代表的な位相変化の統計量を演算する第１位相統計量演算手段と、
　　前記代表的な位相変化と前記各画素における位相配列との相関を演算する第１位相相関演算手段と、
　を有する付記１に記載の画像処理装置。
（付記３）
　前記第１位相統計量演算手段は、各クラスタ内の複数の画素の位相配列の平均を前記代表的な位相変化の統計量とする付記２に記載の画像処理装置。
（付記４）
　前記第１位相統計量演算手段は、前記各画素の位置に対応して重み付けを行う付記３に記載の画像処理装置。
（付記５）
　前記位相相関算出手段は、さらに、クラスタ内での前記代表的な位相変化へのノイズの影響を除去するノイズ除去手段を備える付記２乃至４のいずれか１項に記載の画像処理装置。
（付記６）
　前記ノイズ除去手段は、前記第１位相統計量演算手段が演算した代表的な位相変化の統計量を固有値分解して、その最大固有ベクトルをノイズの影響が除去された前記代表的な位相変化とする付記５に記載の画像処理装置。
（付記７）
　前記ノイズ除去手段は、前記第１位相統計量演算手段が演算した代表的な位相変化の統計量を固有値分解して、その最大固有ベクトルからの数本の固有ベクトルをノイズの影響が除去された前記代表的な位相変化とする付記５に記載の画像処理装置。
（付記８）
　前記ノイズ除去手段は、画像取得時期の離れた位相比較に対応する成分を排除することにより、ノイズの影響が除去された前記代表的な位相変化とする付記５に記載の画像処理装置。
（付記９）
　前記位相相関算出手段は、
　　各クラスタ内の画素の位相配列と前記各画素における位相配列との相関を演算する第２位相相関演算手段と、
　　各クラスタ内での代表的な相関の統計量を演算する第２位相統計量演算手段と、
　を有する付記１に記載の画像処理装置。
（付記１０）
　前記第２位相統計量演算手段は、各クラスタ内の複数の画素における相関の平均を前記代表的な相関の統計量とする付記９に記載の画像処理装置。
（付記１１）
　前記第２位相統計量演算手段は、前記各画素の位置に対応して重み付けを行う付記１０に記載の画像処理装置。
（付記１２）
　前記位相相関算出手段は、前記算出された相関をまとめる処理を含む付記１乃至１１のいずれか１項に記載の画像処理装置。
（付記１３）
　前記複数の画像に渡る各画素から、反射が安定した安定反射点となる候補の画素を特定する安定反射点候補特定手段をさらに備える付記１乃至１２のいずれか１項に記載の画像処理装置。
（付記１４）
　前記安定反射点候補特定手段は、前記位相相関算出手段により算出された相関に基づいて、再度、安定反射点となる候補の画素の特定を更新する付記１３に記載の画像処理装置。
（付記１５）
　前記安定反射点候補特定手段による安定反射点となる候補の画素の特定と、前記位相相関算出手段による相関の算出とを、所定回数の繰り返し、または、安定反射点候補の相関が一定以上となるまで繰り返す付記１４に記載の画像処理装置。
（付記１６）
　前記位相相関算出手段が算出した位相相関を参照して、安定反射点を特定する安定反射点特定手段をさらに備える付記１乃至１５にいずれか１項に記載の画像処理装置。
（付記１７）
　複数の画像に渡る各画素の位相配列を位相配列蓄積手段に蓄積する位相配列蓄積ステップと、
　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリングステップと、
　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出ステップと、
　を含む画像処理方法。
（付記１８）
　複数の画像に渡る各画素の位相配列を位相配列蓄積手段に蓄積する位相配列蓄積ステップと、
　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリングステップと、
　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出ステップと、
　をコンピュータに実行させる画像処理プログラム。
（付記１９）
　複数のＳＡＲ画像に基づいて安定反射点を特定するための画像処理装置と、
　前記安定反射点のデータを用いてＳＡＲ画像を解析するＳＡＲ画像解析装置と、
　を備えるＳＡＲ画像解析システムであって、
　前記画像処理装置は、少なくとも、
　　前記複数のＳＡＲ画像に渡る各画素の位相配列を蓄積する位相配列蓄積手段と、
　　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリング手段と、
　　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出手段と、
　を備えるＳＡＲ画像解析システム。
（付記２０）
　複数のＳＡＲ画像に基づいて安定反射点を特定するための画像処理ステップと、
　前記安定反射点のデータを用いてＳＡＲ画像を解析するＳＡＲ画像解析ステップと、
　を含むＳＡＲ画像解析方法であって、
　前記画像処理ステップは、少なくとも、
　　前記複数のＳＡＲ画像に渡る各画素の位相配列を位相配列蓄積手段に蓄積する位相配列蓄積ステップと、
　　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリングステップと、
　　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出ステップと、
　を含むＳＡＲ画像解析方法。

Claims (20)

  1. 　複数の画像に渡る各画素の位相配列を蓄積する位相配列蓄積手段と、
    　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリング手段と、
    　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出手段と、
    　を備える画像処理装置。
  2. 　前記位相相関算出手段は、
    　　各クラスタ内での代表的な位相変化の統計量を演算する第１位相統計量演算手段と、
    　　前記代表的な位相変化と前記各画素における位相配列との相関を演算する第１位相相関演算手段と、
    　を有する請求項１に記載の画像処理装置。
  3. 　前記第１位相統計量演算手段は、各クラスタ内の複数の画素の位相配列の平均を前記代表的な位相変化の統計量とする請求項２に記載の画像処理装置。
  4. 　前記第１位相統計量演算手段は、前記各画素の位置に対応して重み付けを行う請求項３に記載の画像処理装置。
  5. 　前記位相相関算出手段は、さらに、クラスタ内での前記代表的な位相変化へのノイズの影響を除去するノイズ除去手段を備える請求項２乃至４のいずれか１項に記載の画像処理装置。
  6. 　前記ノイズ除去手段は、前記第１位相統計量演算手段が演算した代表的な位相変化の統計量を固有値分解して、その最大固有ベクトルをノイズの影響が除去された前記代表的な位相変化とする請求項５に記載の画像処理装置。
  7. 　前記ノイズ除去手段は、前記第１位相統計量演算手段が演算した代表的な位相変化の統計量を固有値分解して、その最大固有ベクトルからの数本の固有ベクトルをノイズの影響が除去された前記代表的な位相変化とする請求項５に記載の画像処理装置。
  8. 　前記ノイズ除去手段は、画像取得時期の離れた位相比較に対応する成分を排除することにより、ノイズの影響が除去された前記代表的な位相変化とする請求項５に記載の画像処理装置。
  9. 　前記位相相関算出手段は、
    　　各クラスタ内の画素の位相配列と前記各画素における位相配列との相関を演算する第２位相相関演算手段と、
    　　各クラスタ内での代表的な相関の統計量を演算する第２位相統計量演算手段と、
    　を有する請求項１に記載の画像処理装置。
  10. 　前記第２位相統計量演算手段は、各クラスタ内の複数の画素における相関の平均を前記代表的な相関の統計量とする請求項９に記載の画像処理装置。
  11. 　前記第２位相統計量演算手段は、前記各画素の位置に対応して重み付けを行う請求項１０に記載の画像処理装置。
  12. 　前記位相相関算出手段は、前記算出された相関をまとめる処理を含む請求項１乃至１１のいずれか１項に記載の画像処理装置。
  13. 　前記複数の画像に渡る各画素から、反射が安定した安定反射点となる候補の画素を特定する安定反射点候補特定手段をさらに備える請求項１乃至１２のいずれか１項に記載の画像処理装置。
  14. 　前記安定反射点候補特定手段は、前記位相相関算出手段により算出された相関に基づいて、再度、安定反射点となる候補の画素の特定を更新する請求項１３に記載の画像処理装置。
  15. 　前記安定反射点候補特定手段による安定反射点となる候補の画素の特定と、前記位相相関算出手段による相関の算出とを、所定回数の繰り返し、または、安定反射点候補の相関が一定以上となるまで繰り返す請求項１４に記載の画像処理装置。
  16. 　前記位相相関算出手段が算出した位相相関を参照して、安定反射点を特定する安定反射点特定手段をさらに備える請求項１乃至１５にいずれか１項に記載の画像処理装置。
  17. 　複数の画像に渡る各画素の位相配列を位相配列蓄積手段に蓄積する位相配列蓄積ステップと、
    　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリングステップと、
    　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出ステップと、
    　を含む画像処理方法。
  18. 　複数の画像に渡る各画素の位相配列を位相配列蓄積手段に蓄積する位相配列蓄積ステップと、
    　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリングステップと、
    　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出ステップと、
    　をコンピュータに実行させる画像処理プログラム。
  19. 　複数のＳＡＲ画像に基づいて安定反射点を特定するための画像処理装置と、
    　前記安定反射点のデータを用いてＳＡＲ画像を解析するＳＡＲ画像解析装置と、
    　を備えるＳＡＲ画像解析システムであって、
    　前記画像処理装置は、少なくとも、
    　　前記複数のＳＡＲ画像に渡る各画素の位相配列を蓄積する位相配列蓄積手段と、
    　　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリング手段と、
    　　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出手段と、
    　を備えるＳＡＲ画像解析システム。
  20. 　複数のＳＡＲ画像に基づいて安定反射点を特定するための画像処理ステップと、
    　前記安定反射点のデータを用いてＳＡＲ画像を解析するＳＡＲ画像解析ステップと、
    　を含むＳＡＲ画像解析方法であって、
    　前記画像処理ステップは、少なくとも、
    　　前記複数のＳＡＲ画像に渡る各画素の位相配列を位相配列蓄積手段に蓄積する位相配列蓄積ステップと、
    　　前記位相配列に基づいて、前記各画素を複数のクラスタにクラスタリングするクラスタリングステップと、
    　　前記複数のクラスタの各クラスタ内での位相変化と前記各画素における位相配列との相関を算出する位相相関算出ステップと、
    　を含むＳＡＲ画像解析方法。
PCT/JP2017/043377 2017-12-01 2017-12-01 Ｓａｒ画像解析システム、画像処理装置、画像処理方法および画像処理プログラム WO2019106850A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US16/767,805 US11709254B2 (en) 2017-12-01 2017-12-01 SAR image analysis system, image processing apparatus, image processing method, and image processing program
JP2019556532A JP7020494B2 (ja) 2017-12-01 2017-12-01 Ｓａｒ画像解析システム、画像処理装置、画像処理方法および画像処理プログラム
PCT/JP2017/043377 WO2019106850A1 (ja) 2017-12-01 2017-12-01 Ｓａｒ画像解析システム、画像処理装置、画像処理方法および画像処理プログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2017/043377 WO2019106850A1 (ja) 2017-12-01 2017-12-01 Ｓａｒ画像解析システム、画像処理装置、画像処理方法および画像処理プログラム

Publications (1)

Publication Number Publication Date
WO2019106850A1 true WO2019106850A1 (ja) 2019-06-06

Family

ID=66664839

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2017/043377 WO2019106850A1 (ja) 2017-12-01 2017-12-01 Ｓａｒ画像解析システム、画像処理装置、画像処理方法および画像処理プログラム

Country Status (3)

Country Link
US (1) US11709254B2 (ja)
JP (1) JP7020494B2 (ja)
WO (1) WO2019106850A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021009905A1 (ja) * 2019-07-18 2021-01-21 日本電気株式会社 画像処理装置および画像処理方法
WO2021009904A1 (ja) * 2019-07-18 2021-01-21 日本電気株式会社 画像処理装置および画像処理方法
WO2022070242A1 (ja) * 2020-09-29 2022-04-07 日本電気株式会社 画像解析装置および画像解析方法

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019220574A1 (ja) * 2018-05-16 2019-11-21 日本電気株式会社 合成開口レーダ信号解析装置、合成開口レーダ信号解析方法および合成開口レーダ信号解析プログラム
JP7281660B2 (ja) * 2018-08-30 2023-05-26 パナソニックＩｐマネジメント株式会社 変位計測装置及び変位計測方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0894747A (ja) * 1994-09-28 1996-04-12 Mitsubishi Electric Corp 合成開口レーダ画像による地形抽出装置
JPH09281231A (ja) * 1996-04-15 1997-10-31 Mitsubishi Electric Corp 信号処理装置
JP2004191053A (ja) * 2002-12-06 2004-07-08 Mitsubishi Electric Corp 合成開口レーダ装置及び数値標高モデル作成方法
JP2008014709A (ja) * 2006-07-04 2008-01-24 Mitsubishi Electric Corp 画像処理装置
JP2008157764A (ja) * 2006-12-25 2008-07-10 Mitsubishi Electric Corp レーダ画像処理装置
JP2009289111A (ja) * 2008-05-30 2009-12-10 Mitsubishi Electric Corp 画像処理装置
WO2015008554A1 (ja) * 2013-07-19 2015-01-22 国立大学法人東北大学 合成開口処理を伴うセンサ、そのセンサの処理方法、および、プログラム
JP2017072473A (ja) * 2015-10-07 2017-04-13 三菱電機株式会社 画像処理装置および画像処理方法

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ITMI991154A1 (it) * 1999-05-25 2000-11-25 Milano Politecnico Procedimento per misure radar di spostamento di aere urbane e zone franose
JP2006221520A (ja) * 2005-02-14 2006-08-24 Chuden Gijutsu Consultant Kk 画像特徴量抽出システムおよび方法
EP2053423A4 (en) 2006-08-03 2012-05-09 Pasco Corp DISASTERS-AGAINST ACTION SUPPORT METHOD
EP3876510A1 (en) * 2008-05-20 2021-09-08 FotoNation Limited Capturing and processing of images using monolithic camera array with heterogeneous imagers
US8866920B2 (en) * 2008-05-20 2014-10-21 Pelican Imaging Corporation Capturing and processing of images using monolithic camera array with heterogeneous imagers
ES2566132T3 (es) * 2008-07-04 2016-04-11 Telespazio S.P.A. Identificación y análisis de dispersores persistentes en series de imágenes de SAR
US9289191B2 (en) * 2011-10-12 2016-03-22 Seno Medical Instruments, Inc. System and method for acquiring optoacoustic data and producing parametric maps thereof
US9730587B2 (en) * 2011-11-02 2017-08-15 Seno Medical Instruments, Inc. Diagnostic simulator
US9814394B2 (en) * 2011-11-02 2017-11-14 Seno Medical Instruments, Inc. Noise suppression in an optoacoustic system
US20140005544A1 (en) * 2011-11-02 2014-01-02 Seno Medical Instruments, Inc. System and method for providing selective channel sensitivity in an optoacoustic imaging system
US9445786B2 (en) * 2011-11-02 2016-09-20 Seno Medical Instruments, Inc. Interframe energy normalization in an optoacoustic imaging system
US9743839B2 (en) * 2011-11-02 2017-08-29 Seno Medical Instruments, Inc. Playback mode in an optoacoustic imaging system
US20130116538A1 (en) * 2011-11-02 2013-05-09 Seno Medical Instruments, Inc. Optoacoustic imaging systems and methods with enhanced safety
US9239384B1 (en) * 2014-10-21 2016-01-19 Sandia Corporation Terrain detection and classification using single polarization SAR
US10254399B1 (en) * 2015-12-02 2019-04-09 National Technology & Engineering Solutions Of Sandia, Llc Small mover target detection in synthetic aperture radar imagery
JP7077961B2 (ja) * 2016-12-27 2022-05-31 日本電気株式会社 画像解析装置、画像解析方法及びプログラム
US11397261B2 (en) * 2017-10-30 2022-07-26 Nec Corporation Image processing apparatus, image processing method, image processing program, and image processing system

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0894747A (ja) * 1994-09-28 1996-04-12 Mitsubishi Electric Corp 合成開口レーダ画像による地形抽出装置
JPH09281231A (ja) * 1996-04-15 1997-10-31 Mitsubishi Electric Corp 信号処理装置
JP2004191053A (ja) * 2002-12-06 2004-07-08 Mitsubishi Electric Corp 合成開口レーダ装置及び数値標高モデル作成方法
JP2008014709A (ja) * 2006-07-04 2008-01-24 Mitsubishi Electric Corp 画像処理装置
JP2008157764A (ja) * 2006-12-25 2008-07-10 Mitsubishi Electric Corp レーダ画像処理装置
JP2009289111A (ja) * 2008-05-30 2009-12-10 Mitsubishi Electric Corp 画像処理装置
WO2015008554A1 (ja) * 2013-07-19 2015-01-22 国立大学法人東北大学 合成開口処理を伴うセンサ、そのセンサの処理方法、および、プログラム
JP2017072473A (ja) * 2015-10-07 2017-04-13 三菱電機株式会社 画像処理装置および画像処理方法

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021009905A1 (ja) * 2019-07-18 2021-01-21 日本電気株式会社 画像処理装置および画像処理方法
WO2021009904A1 (ja) * 2019-07-18 2021-01-21 日本電気株式会社 画像処理装置および画像処理方法
JPWO2021009905A1 (ja) * 2019-07-18 2021-01-21
JPWO2021009904A1 (ja) * 2019-07-18 2021-01-21
JP7188594B2 (ja) 2019-07-18 2022-12-13 日本電気株式会社 画像処理装置および画像処理方法
JP7188595B2 (ja) 2019-07-18 2022-12-13 日本電気株式会社 画像処理装置および画像処理方法
US11846702B2 (en) 2019-07-18 2023-12-19 Nec Corporation Image processing device and image processing method
WO2022070242A1 (ja) * 2020-09-29 2022-04-07 日本電気株式会社 画像解析装置および画像解析方法
JP7416277B2 (ja) 2020-09-29 2024-01-17 日本電気株式会社 画像解析装置および画像解析方法

Also Published As

Publication number Publication date
US11709254B2 (en) 2023-07-25
US20200348411A1 (en) 2020-11-05
JPWO2019106850A1 (ja) 2020-12-03
JP7020494B2 (ja) 2022-02-16

Similar Documents

Publication Publication Date Title
WO2019106850A1 (ja) Ｓａｒ画像解析システム、画像処理装置、画像処理方法および画像処理プログラム
CN108596108B (zh) 基于三元组语义关系学习的航拍遥感图像变化检测方法
Tao et al. Dense point cloud extraction from UAV captured images in forest area
CN111797717A (zh) 一种高速高精度的sar图像船只检测方法
JP6397379B2 (ja) 変化領域検出装置、方法、及びプログラム
CN106056070B (zh) 基于低秩矩阵恢复和稀疏表示的sar目标识别方法
CN110826428A (zh) 一种高速的sar图像中船只检测方法
WO2018123748A1 (ja) 画像解析装置、画像解析方法及びコンピュータ読み取り可能記録媒体
US11841421B2 (en) Synthetic aperture radar image analysis system, synthetic aperture radar image analysis method, and synthetic aperture radar image analysis program
CN106772273B (zh) 一种基于动态孔径的sar虚假目标干扰抑制方法及系统
CN108230375A (zh) 基于结构相似性快速鲁棒的可见光图像与sar图像配准方法
WO2019087673A1 (ja) 画像処理装置、画像処理方法、画像処理プログラムおよび画像処理システム
CN103093243B (zh) 高分辨率全色遥感图像云判方法
CN110113560A (zh) 视频智能联动的方法及服务器
JP7152131B2 (ja) 画像処理装置、画像処理方法及びプログラム
Feng et al. A novel saliency detection method for wild animal monitoring images with WMSN
Xia et al. Quality assessment for remote sensing images: approaches and applications
CN111735447B (zh) 一种仿星敏式室内相对位姿测量系统及其工作方法
CN105678716B (zh) 一种地基sar大气干扰相位校正方法及装置
CN110728311A (zh) 一种图像处理方法、装置和存储介质
EP3371543A1 (en) Improvements in and relating to missile targeting
CN111126508A (zh) 一种基于hopc改进的异源图像匹配方法
CN106530326B (zh) 基于影像纹理特征和dsm的变化检测方法
Latif et al. Self-Organizing Maps for processing of data with missing values and outliers: application to remote sensing images
CN110399832B (zh) 基于相干性的TomoSAR植被病虫害监测方法及装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17933746

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2019556532

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 17933746

Country of ref document: EP

Kind code of ref document: A1