WO2022054186A1 - 推定装置、推定方法及び記憶媒体 - Google Patents

推定装置、推定方法及び記憶媒体 Download PDF

Info

Publication number
WO2022054186A1
WO2022054186A1 PCT/JP2020/034216 JP2020034216W WO2022054186A1 WO 2022054186 A1 WO2022054186 A1 WO 2022054186A1 JP 2020034216 W JP2020034216 W JP 2020034216W WO 2022054186 A1 WO2022054186 A1 WO 2022054186A1
Authority
WO
WIPO (PCT)
Prior art keywords
point
target point
measurement
ground surface
estimation
Prior art date
Application number
PCT/JP2020/034216
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 JP2022548307A priority Critical patent/JPWO2022054186A5/ja
Priority to PCT/JP2020/034216 priority patent/WO2022054186A1/ja
Publication of WO2022054186A1 publication Critical patent/WO2022054186A1/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

Definitions

  • the present disclosure relates to a technique for estimating an estimated value from which noise components have been removed from the measurement results, and more particularly to a technique for estimating an estimated value from the results measured by a synthetic aperture radar.
  • Measurements measured by sensors such as changes in ground surface height measured by synthetic aperture radar, generally include noise.
  • Patent Document 1 discloses one of the techniques for obtaining an estimated value from which noise is removed from the measurement result.
  • a device for adding and averaging the signals of pixels at the same position of image signals obtained by a plurality of infrared image pickup elements arranged so as to obtain the same image is described.
  • Patent Document 2 describes noise from three-dimensional point cloud data, which is a set of points for which three-dimensional coordinates are obtained, which is obtained as a result of measurement of a laser scanning device that measures the distance to a measurement point that reflects laser light.
  • a method of performing the removal is disclosed. The method divides a curved surface formed by measurement points corresponding to a two-dimensional rectangular region set for three-dimensional point cloud data into a set of triangles having measurement points as vertices. The method is to obtain a composite normal vector that combines the normal vectors of all the triangles into one. Then, in that method, when the standard deviation of the measured values in the coordinate axis direction when the coordinates are converted so that the composite normal vector matches the coordinate axes exceeds a predetermined value, one of those measurement points is treated as noise.
  • the techniques of Patent Document 1 the signals of pixels at the same position of a plurality of images captured by a plurality of sensors at the same time are added and averaged.
  • the technique of Patent Document 2 is a technique of specifying noise at a plurality of measurement points where distances are measured at the same timing.
  • the technique of Patent Document 2 when noise is included in the measured value of a predetermined measurement point, the three-dimensional coordinates of the measurement point cannot be estimated.
  • the techniques of Patent Documents 1 and 2 cannot estimate changes in measured values and the like. Therefore, the techniques of Patent Documents 1 and 2 cannot be applied to, for example, the estimation of the height of the ground surface. In the techniques of Patent Documents 1 and 2, it is not possible to accurately estimate the change in the measured value when the change in the measured value of the target varies depending on the difference in the properties of the area to be measured.
  • One of the purposes of the present disclosure is to provide an estimation device and the like that can improve the accuracy of estimation of changes even when changes in measured values vary depending on a target area.
  • the estimation device is a determination means for determining selection conditions based on the geospatial information of the target point among a plurality of measurement points whose changes in the height of the ground surface are measured by the synthetic aperture radar. , The height of the ground surface at the target point based on the measured value of the transition of the ground surface height at the target point and the measured value of the transition of the ground surface height at the measurement point satisfying the selection condition. It is provided with an estimation means for estimating an estimated value of change in.
  • selection conditions are determined based on the geospatial information of the target point among a plurality of measurement points whose changes in the height of the ground surface are measured by the synthetic aperture radar, and the target is described above.
  • the change in the height of the ground surface at the target point is based on the measured value of the transition of the height of the ground surface at the point and the measured value of the transition of the height of the ground surface at the measurement point satisfying the selection condition. Estimate the estimated value.
  • the storage medium is a determination process for determining selection conditions based on the geospatial information of the target point among a plurality of measurement points whose changes in the height of the ground surface are measured by the synthetic aperture radar. , The height of the ground surface at the target point based on the measured value of the transition of the ground surface height at the target point and the measured value of the transition of the ground surface height at the measurement point satisfying the selection condition. It stores an estimation process that estimates the estimated value of the change in and a program that causes the computer to execute.
  • One aspect of the present disclosure is also realized by the program described above.
  • the present disclosure has an effect that the accuracy of estimation of the change can be improved even if the change of the measured value varies depending on the target area.
  • FIG. 1 is a block diagram showing an example of the configuration of the estimation system according to the first embodiment of the present disclosure.
  • FIG. 2 is a flowchart showing an example of the operation of the estimation device according to the first embodiment of the present disclosure.
  • FIG. 3 is a block diagram showing an example of the configuration of the estimation system according to the second embodiment of the present disclosure.
  • FIG. 4 is a flowchart showing an example of the operation of the estimation device according to the second embodiment of the present disclosure.
  • FIG. 5 is a block diagram showing an example of the configuration of the estimation device according to the third embodiment of the present disclosure.
  • FIG. 6 is a flowchart showing an example of the operation of the estimation device according to the third embodiment of the present disclosure.
  • FIG. 7 is a block diagram showing an example of a computer configuration capable of realizing an apparatus such as an estimation apparatus according to the embodiment of the present disclosure.
  • FIG. 1 is a block diagram showing an example of the configuration of the estimation system 1 according to the first embodiment of the present disclosure.
  • the estimation system 1 includes an estimation device 100, a measurement data storage device 200, a geospatial information storage device 300, and an output destination device 400.
  • the estimation device 100 includes a determination unit 110, an estimation unit 120, an output unit 140, a measurement data acquisition unit 150, and a geospatial information acquisition unit 160.
  • the measurement data storage device 200 includes a measurement data storage unit 210 and a first input / output unit 220.
  • the geospatial information storage device 300 includes a geospatial information storage unit 310 and a second input / output unit 320.
  • the determination unit 110 of the estimation device 100 is based on the geospatial information stored in the geospatial information storage device 300 of the target point among the plurality of measurement points whose changes in the height of the ground surface are measured by the synthetic aperture radar. To determine the selection conditions.
  • the geospatial information is, for example, information representing at least one of the state of the ground surface and the underground state of the ground surface. Geospatial information is related, for example, to the geology of the measurement point and the presence or absence of buildings.
  • the estimation unit 120 determines the height of the ground surface at the target point based on the measured value of the transition of the height of the ground surface at the target point and the measured value of the transition of the height of the ground surface at the measurement point satisfying the selection condition. Estimate the estimated value of the change in. Therefore, even if the change in the measured value of the target varies depending on the geology of the area to be measured and the presence or absence of a building, it is possible to estimate the change in the measured value with high accuracy.
  • the measured value of the transition of the height of the ground surface at the measurement point is stored in the measurement data storage device 200 as measurement data.
  • the output unit 140 outputs an estimated value of the change in the height of the ground surface at the target point to the output destination device 400.
  • Measurement data storage unit 210 stores the measurement data.
  • the measurement data is, for example, data showing the transition of the result of the measurement of the height of the ground surface at a plurality of points by the synthetic aperture radar (hereinafter referred to as the transition of the height of the ground surface).
  • the height of the ground surface may be the height of the ground surface from a predetermined reference plane.
  • the height of the ground surface of the measured data is multiple on the ground surface obtained by observation using, for example, a radar mounted on a flying object such as an artificial satellite or an aircraft as a synthetic aperture radar (SAR). It is the height at the measurement point of.
  • SAR synthetic aperture radar
  • the data representing the transition of the height of the ground surface may be, for example, data including a plurality of combinations of a value representing the height and data representing the time point at which the height is obtained by observation.
  • the unit of data representing the time point may be appropriately determined.
  • the data representing a time point may represent a date, or may represent a date and a time.
  • the unit of time may also be appropriately determined.
  • the measurement data may be data including a combination of position data representing the position of each measurement point of a plurality of measurement points and data representing the transition of the height of the ground surface at the measurement points.
  • the information representing the position of a point (that is, a measurement point) on the ground surface may be, for example, latitude and longitude information.
  • the information indicating the position may be other information that can identify the position on the ground surface. In the following, the information indicating the position of the point is referred to as the point information.
  • the data showing the transition of the height of the ground surface at the measurement point on the ground surface will be referred to as the height transition data.
  • the height transition data included in the measurement data is data in which noise is not removed.
  • the height transition data of the measurement data includes noise and the like.
  • the height transition data may be associated with the point information of the measurement point where the height transition represented by the height transition data is measured.
  • the measurement data includes point information of a plurality of measurement points and height transition data at each of the plurality of measurement points.
  • the first input / output unit 220 reads out the measurement data stored in the measurement data storage unit 210 and supplies it to the estimation device 100.
  • the first input / output unit 220 may receive, for example, a request for measurement data from the estimation device 100.
  • the request for measurement data may include information that identifies a range on the ground surface.
  • the first input / output unit 220 may read the measurement data of the measurement point included in the received request from the measurement data storage unit 210, and supply the read measurement data to the estimation device 100.
  • Geospatial Information Storage Unit 310 stores geospatial information.
  • the geospatial information is, for example, information representing at least one of the state of the ground surface and the underground state of the ground surface.
  • the geospatial information may be at least one of the information obtained from a so-called geographic information system (Geographic Information System).
  • Geographic Information System Geographic Information System
  • the geospatial information may be data obtained by observation from an artificial satellite, an aircraft, or the like.
  • the geospatial information may be data obtained by a field survey.
  • the geospatial information may be information representing the result of analysis based on the data obtained by measurement or survey.
  • the geospatial information may be artificially determined information based on the data obtained by measurement or survey.
  • Geographical space information may be referred to as GIS (Geographic Information System) data.
  • the geospatial information may be acquired from the geospatial information system in advance and stored in the geospatial information storage unit 310.
  • a plurality of types of geospatial information may be stored in the geospatial information storage unit 310.
  • the geospatial information may be expressed in a format that can specify the value of the geospatial information at the point specified by the point information (for example, latitude and longitude). Specific examples of geospatial information will be described in detail later.
  • Second input / output unit 320 extracts the value of the geospatial information from the geospatial information stored in the geospatial information storage unit 310, and supplies the extracted value of the geospatial information to the estimation device 100.
  • the second input / output unit 320 may receive, for example, a request for a value of geospatial information including point information from the estimation device 100.
  • the second input / output unit 320 may extract the value of the geospatial information of the point specified by the point information of the received request from the geospatial information stored in the geospatial information storage unit 310.
  • the request for the value of the geospatial information may include the point information of a plurality of points.
  • the second input / output unit 320 stores the value of the geospatial information of the point specified by the point information of a plurality of points included in the request for the value of the geospatial information in the geospatial information storage unit 310. It may be extracted from the spatial information. The second input / output unit 320 may supply the extracted geospatial information values of a plurality of points to the estimation device 100.
  • the second input / output unit 320 may extract the value of the predetermined type of geospatial information at the point specified by the point information.
  • the second input / output unit 320 may extract the values of all types of geospatial information stored in the geospatial information storage unit 310 at the points specified by the point information. If there is geospatial information that does not have a value at the point specified by the point information, the second input / output unit 320 does not have to extract the value of the geospatial information.
  • the second input / output unit 320 sets the value of the geospatial information to a value indicating that the value does not exist (for example, 0, etc.). May be set to.
  • the geospatial information may be represented by, for example, a state-representing value for each mesh on which the ground surface is delimited.
  • the second input / output unit 320 uses the value of the geospatial information representing the state in the mesh including the position specified by the point information as the value of the geospatial information of the point specified by the point information. Extract.
  • the size and shape of the mesh may be determined for each type of geospatial information.
  • the geospatial information may be expressed in other formats.
  • the geospatial information may be represented, for example, by a boundary line between regions having different states and a value representing a state within the region separated by the boundary line.
  • the second input / output unit 320 extracts a value representing the state in the area including the position specified by the point information as the value of the geospatial information of the point specified by the point information.
  • the format of the geospatial information may be defined for each type of geospatial information.
  • Specific geospatial information includes, for example, type of embankment site, average slope angle, average precipitation (for example, average annual precipitation), surface geology, steep slope designation, sediment disaster warning area designation, liquefaction risk, etc. It may be whether or not a rainwater infiltration basin can be installed, easiness of shaking during an earthquake, lowland where drainage is difficult, land use in urban areas, natural terrain classification, artificial terrain classification, surface geology, riverbed, facility information (presence or absence of construction, etc.).
  • the geospatial information of this embodiment is not necessarily limited to these examples.
  • the geospatial information of this embodiment may not include some of these examples.
  • the type of embankment construction site may represent the method of embankment, which is determined by the shape of the ground surface on which the embankment was made.
  • the types of embankment construction sites are, for example, “valley-filled embankment”, which is an embankment in which valleys and swamps are filled with embankment, and "belly-type embankment”, which is an embankment made on slopes.
  • the type of embankment site may further represent the scale of the embankment.
  • the type of embankment satisfying the standard for example, embankment having an area of 3000 square meters or more
  • the type of embankment that does not meet the criteria may be "valley-filled embankment”.
  • the angle of the embankment that meets the criteria is 20 degrees or more with respect to the horizontal plane, and the height of the embankment is 5.
  • the type of embankment) that is metric or larger may be a large-scale embankment.
  • the type of embankment that does not meet the criteria may be "belly-type embankment".
  • the value of the type of embankment construction site may be, for example, one of different numerical values assigned in advance to "valley-filled embankment" or "belly-filled embankment”.
  • the average inclination angle may be, for example, data of the average inclination angle of the ground surface calculated in mesh units.
  • the value of the average inclination angle may be the calculated average inclination angle of the ground surface.
  • the average precipitation may be, for example, data of the average precipitation on the ground surface calculated in mesh units.
  • the average precipitation value may be the calculated average precipitation on the ground surface.
  • the surface geology may be data representing the surface geology (in other words, the type of geology) on the surface of the earth.
  • the type of geology may be predetermined. Different numerical values may be assigned to each type of geology in advance.
  • the value of the surface geology may be any one of the numerical values assigned to the geology in advance.
  • the steep slope designation may be data indicating whether or not it is designated as a steep slope by, for example, a local government.
  • the value of the steep slope designation may be, for example, a numerical value indicating that it is designated as a steep slope, or a numerical value indicating that it is not designated as a steep slope. As these numerical values, numerical values different from each other may be appropriately determined in advance.
  • the sediment-related disaster warning area designation may indicate whether or not it has been designated as a sediment-related disaster warning area by, for example, a local government.
  • the value of the sediment-related disaster warning area designation may be, for example, a numerical value indicating that it is designated as a sediment-related disaster warning area or a numerical value indicating that it is not designated as a sediment-related disaster warning area. As these numerical values, numerical values different from each other may be predetermined.
  • the liquefaction risk may be, for example, data indicating the degree of risk of land liquefaction.
  • the liquefaction risk value may be a numerical value indicating the degree of liquefaction risk of the land.
  • the value of the liquefaction risk may be any one of a plurality of different numerical values representing different degrees.
  • the numerical value indicating the degree of danger may be appropriately determined in advance.
  • Whether or not a rainwater infiltration basin can be installed is information indicating whether or not an infiltration facility can be installed based on the "infiltration facility installation judgment map" that indicates the result of determining whether or not an infiltration facility can be installed based on, for example, topography, soil quality, and groundwater level. It's okay.
  • the value of whether or not the rainwater infiltration basin can be installed may be a numerical value indicating that the installation is possible or a numerical value indicating that the installation is not possible. As these numerical values, numerical values different from each other may be appropriately determined in advance.
  • the easiness of shaking at the time of an earthquake may be, for example, data indicating the degree of easiness of shaking of the ground surface when an earthquake occurs.
  • the value of easiness of shaking at the time of an earthquake may be a numerical value indicating the degree of easiness of shaking of the ground surface in the event of an earthquake.
  • the value of easiness of shaking at the time of an earthquake may be any one of a plurality of numerical values indicating the degree of easiness of shaking of the ground surface in the event of an earthquake. Numerical values indicating the degree of easiness of shaking of the ground surface may be appropriately determined in advance.
  • the difficult-to-drain lowland may represent, for example, whether or not the land is a difficult-to-drain lowland estimated from the altitude of the land or the difference in altitude from the surrounding area.
  • the value of the lowland where drainage is difficult may be a numerical value indicating that it is a lowland where drainage is difficult, or a numerical value indicating that it is not a lowland where drainage is difficult. These numerical values may be appropriately determined in advance.
  • Urban land use may be a type of land use in the area designated as a city.
  • the type of land use in urban land use may be read, for example, from satellite images.
  • a land use type selected from a plurality of predetermined types may be set for an area included in a city area. Different numerical values may be appropriately assigned to each of the plurality of predetermined types.
  • the value of the land use type set for the area may be a numerical value assigned to the type.
  • the natural terrain classification may be, for example, a type of terrain in a place that is not a building built by humans.
  • a plurality of terrain types that can be set as natural terrain classification may be appropriately determined in advance.
  • the terrain type selected from a plurality of terrain types predetermined as the types that can be set as the natural terrain classification may be set.
  • Different numerical values may be assigned to each of the plurality of types.
  • the terrain value in the natural terrain classification set in the area may be a numerical value assigned to the type of terrain set in the area.
  • the artificial terrain classification may be, for example, a type of terrain in a place where a human has modified the terrain or a building constructed by a human.
  • a plurality of terrain types that can be set as artificial terrain classification may be appropriately determined in advance.
  • the terrain type selected from a plurality of terrain types predetermined as the types that can be set as the artificial terrain classification may be set.
  • Different numerical values may be assigned to each of the plurality of types.
  • the terrain value in the artificial terrain classification set in the area may be a numerical value assigned to the type of terrain set in the area.
  • the surface geology may represent, for example, the type of soil on the surface of the earth.
  • a plurality of soil types may be appropriately determined in advance.
  • different numerical values which are appropriately determined in advance, may be assigned to each of the plurality of soil types.
  • the type of soil based on the results of the survey may be set for the area.
  • the surface geological value of the area may be a numerical value assigned to the soil type set in the area.
  • the specific geospatial information may be building information.
  • the building information may represent, for example, the range of individual buildings.
  • the value of the building information may be, for example, assigned different identification numbers to each building.
  • the value of the building information at a point not included in the range of the building may be a value indicating something that is not a building (for example, 0 or the like).
  • Measurement data acquisition unit 150 acquires measurement data from the measurement data storage device 200.
  • the measurement data acquisition unit 150 may, for example, transmit a request for measurement data in a designated range (hereinafter referred to as a designated range) to the measurement data storage device 200.
  • the designated range may be predetermined.
  • the designated range may be designated by the user of the estimation device 100, for example, by an input unit such as a keyboard or a mouse of the estimation device 100. The input unit is omitted in FIG.
  • the measurement data acquisition unit 150 receives the measurement data at the measurement points in the designated range from the measurement data storage device 200.
  • the measurement data acquisition unit 150 sends the received measurement data to the estimation unit 120.
  • Estimating unit 120 receives the measurement data in the designated area from the measurement data acquisition unit 150.
  • the estimation unit 120 selects one measurement point as the target point from the measurement points for which the height transition data is obtained by the measurement in the received measurement data.
  • the estimation unit 120 sends the point information of the target point to the determination unit 110.
  • the determination unit 110 will be described in detail later.
  • the estimation unit 120 receives from the determination unit 110 a selection condition for selecting a measurement point, which is determined based on the point information of the target point.
  • the estimation unit 120 receives information representing the determined selection condition from the determination unit 110.
  • the selection condition is determined by the determination unit 110.
  • the selection condition is determined so that the target point satisfies the selection condition, for example.
  • the selection criteria are further determined based on the value of the geospatial information at the point of interest.
  • the type of geospatial information associated with the value of the geospatial information may be predetermined.
  • the estimation unit 120 sends the value of the geospatial information to the determination unit 110 in order to cause the determination unit 110 to determine the selection condition.
  • the estimation unit 120 acquires the value of the geospatial information at each of the measurement points included in the target range from the geospatial information storage device 300 via the geospatial information acquisition unit 160.
  • the estimation unit 120 requests, for example, a geospatial information request in which information for specifying the type of geospatial information and the point information of the measurement point are specified, via the geospatial information acquisition unit 160, as a geospatial information storage device. Send to 300.
  • the estimation unit 120 requests the geospatial information including the information for specifying the type of the geospatial information and the point information of the measurement point to the geospatial information storage device 300 via the geospatial information acquisition unit 160.
  • the estimation unit 120 receives the value of the geospatial information at the measurement point from the geospatial information storage device 300 via the geospatial information acquisition unit 160.
  • the estimation unit 120 may send a request for geospatial information to the geospatial information acquisition unit 160.
  • the geospatial information acquisition unit 160 receives a request for geospatial information, and transmits the received request for geospatial information to the geospatial information storage device 300.
  • the geospatial information acquisition unit 160 receives the value of the geospatial information from the geospatial information storage device 300, and sends the received value of the geospatial information to the estimation unit 120.
  • the estimation unit 120 receives the value of the geospatial information from the geospatial information acquisition unit 160 and outputs it to the determination unit 110.
  • the estimation unit 120 may send the value of the geospatial information of the target point to the determination unit 110 together with the point information of the target point.
  • the determination unit 110 may receive the point information of the target point and the value of the geospatial information of the target point from the estimation unit 120.
  • the determination unit 110 acquires the value of the geospatial information at the point indicated by the geospatial information of the received target point. It may be acquired from the geospatial information storage device 300 via the unit 160.
  • the determination unit 110 sends a request for geospatial information to the geospatial information acquisition unit 160, which specifies information for specifying the type of geospatial information and the point information of the target point.
  • the determination unit 110 transmits a request for geospatial information including information for specifying the type of geospatial information and location information of the target point to the geospatial information acquisition unit 160.
  • the geospatial information acquisition unit 160 transmits the received geospatial information request to the geospatial information storage device 300.
  • the geospatial information acquisition unit 160 receives the value of the geospatial information at the target point from the geospatial information storage device 300, and sends the received value of the geospatial information at the target point to the determination unit 110. Then, the determination unit 110 receives the value of the geospatial information at the target point from the geospatial information acquisition unit 160.
  • the determination unit 110 determines the selection condition based on the value of the geospatial information at the target point.
  • the estimation unit 120 receives information representing the determined selection condition from the determination unit 110.
  • the estimation unit 120 selects a measurement point satisfying the selection condition from the measurement points from which the height transition data is obtained in the received measurement data.
  • the estimation unit 120 estimates the estimated change transition data of the target point based on the height transition data of the selected measurement point and the height transition data of the target point. The estimated change transition data will be described in detail later. If the selection condition is determined so that the target point satisfies the selection condition, the target point is included in the selected measurement point.
  • the estimation unit 120 estimates the estimated change transition data of the target point based on the height transition data of the selected measurement point.
  • the estimation unit 120 may calculate the transition of the height change at the plurality of measurement points from the height transition data at the plurality of measurement points included in the received measurement data.
  • the estimation unit 120 determines the height of the ground surface at each time point where the height of the ground surface at the previous time point is obtained as the transition of the change in the height of the ground surface at a certain measurement point, and the ground surface at the immediately preceding time point.
  • the change from the height of may be calculated.
  • the time point at which the measurement is performed is referred to as the measurement time point.
  • the change in height at a certain measurement time point may be, for example, a value obtained by subtracting the height of the ground surface at the measurement time point immediately before the measurement time point from the height of the ground surface at the measurement time point.
  • the estimation unit 120 calculates, as a change in height at a certain measurement time point, a change in the height of the ground surface at the time of the measurement from the height of the ground surface at the time of the measurement immediately before the time of the measurement, per predetermined period.
  • the predetermined period may be appropriately set.
  • the predetermined period may be, for example, one month.
  • the predetermined period may be, for example, one year.
  • the estimation unit 120 multiplies, for example, the value obtained by subtracting the height of the ground surface at the measurement time immediately before the measurement time from the height of the ground surface at the measurement time, and multiplying the value representing a predetermined time interval. Further, it may be divided by a value representing the period from the measurement time point and the immediate time point immediately before the measurement time point.
  • the estimation unit 120 may use the obtained value as a change in height at a certain measurement time point.
  • the difference between the measured values of the height of the ground surface measured at two time points, or the change per predetermined time interval calculated from the difference is also referred to as the measured value of the change in height. ..
  • the change data includes height transition data of a plurality of measurement points.
  • the change transition data at a certain measurement point includes data showing the change in the height of the ground surface at each of the measurement points at which the previous measurement point exists.
  • the change transition data may be a plurality of combinations of data representing a measurement time point and data representing a change in the height of the ground surface at the measurement time point.
  • the change data may include, for example, a combination of point information and change transition data of a plurality of measurement points.
  • the change transition data represents the transition of the change in the height of the ground surface, that is, the time series of the change in the height of the ground surface calculated from the measured value of the height of the ground surface.
  • the estimation unit 120 includes a change in the height of the ground surface at a certain measurement time point and a change in the height of the ground surface at the selected measurement point at the same measurement time point as the measurement time point, which is included in the change transition data at the target point. To identify.
  • the estimation unit 120 determines, for example, a change in the height of the ground surface at a certain measurement time point and a change in the height of the ground surface specified as a change in the height of the ground surface at the same measurement time point as the measurement time point. And, the statistical value may be calculated.
  • the statistical value may be, for example, a mean value, a median value, an intermediate value, a mode value, or the like.
  • the estimation unit 120 estimates that the calculated statistical value is an estimated value of the change in the height of the ground surface at the target point. If there is no measurement point other than the target point that meets the target conditions at the target point included in the designated area, the change in the height of the ground surface at the target point is the change in the height of the ground surface at the target point. It may be an estimated value.
  • the estimation unit 120 similarly estimates each estimated value of the change in the height of the ground surface at the target point.
  • the estimation unit 120 estimates that the time series of the estimated value of the change in the height of the ground surface at the target point is the estimated value of the change in the height of the ground surface at the target point.
  • the estimation unit 120 may estimate the estimated value of the change in the height of the ground surface at the measurement point included in the designated area by selecting the measurement point as the target point and performing the above-mentioned processing.
  • the estimation unit 120 may similarly estimate the estimated value of the transition of the change in the height of the ground surface at each of the measurement points included in the designated area.
  • the estimation unit 120 generates estimated change transition data, which is data on the transition of the estimated value of the change in the height of the ground surface at each of the measurement points within the designated range.
  • the estimated change transition data may be data including a combination of the data representing the measurement time point at each measurement time point in which the estimated value of the height change is obtained and the estimated value of the change in the height of the ground surface.
  • the estimation unit 120 includes estimation of each of the measurement points within the designated range, including point information indicating the position of the measurement point and estimated change transition data which is data of the transition of the estimated value of the change in the height of the ground surface. Generate data.
  • the estimation unit 120 is a combination of the point information indicating the position of the measurement point and the estimated change transition data representing the transition of the estimated value of the change in the height of the ground surface of each of the measurement points included in the designated area (that is, the above-mentioned). Estimated data) is sent to the output unit 140.
  • the determination unit 110 receives the point information of the target point from the estimation unit 120.
  • the determination unit 110 may further receive the value of the geospatial information of the target point from the estimation unit 120.
  • the determination unit 110 may receive the value of the geospatial information of the target point from the geospatial information storage device 300 via the geospatial information acquisition unit 160.
  • the estimation unit 120 is configured to send the value of the geodata information of the target point to the determination unit 110
  • the determination unit 110 is configured to receive the value of the geodata information of the target point from the estimation unit 120. Will be explained.
  • the determination unit 110 may determine the selection condition based on the point information of the target point.
  • the determination unit 110 may determine the selection condition based on the value of the geospatial information of the target point.
  • the determination unit 110 may determine the selection condition based on the point information of the target point and the value of the geospatial information of the target point.
  • the determination unit 110 sends information representing the determined selection condition, in other words, the determined selection condition, to the estimation unit 120.
  • the selection condition based on the point information of the target point (hereinafter referred to as the first selection condition) may be, for example, that the distance from the position of the target point is a predetermined distance or less.
  • the selection condition based on the point information of the target point and the geospatial information of the target point satisfies the first selection condition and the selection condition based on the geospatial information of the target point (hereinafter referred to as the second selection condition). May be.
  • the second selection condition may be, for example, that the value of the predetermined type of geospatial information is the same as the value of the predetermined type of geospatial information at the target point.
  • the second selection condition may be, for example, that the value of the predetermined type of geospatial information is closer to the value of the predetermined type of geospatial information at the target point in comparison with the predetermined standard. In this case, if the distance can be defined for the value of the predetermined type of geodata information, the second selection condition is from the value of the predetermined type of geodata information to the value of the predetermined type of geodata information of the target point. The distance to the value may be less than the predetermined value.
  • the values of geospatial information may be divided into a plurality of categories. In that case, the second selection condition may be, for example, that the value of the geospatial information is included in the category including the value of the geospatial information of the target point.
  • the second selection condition may be a condition relating to the values of a plurality of types of geospatial information.
  • the second selection condition may be that, for example, the value of the urban land use in the geospatial information is the same as the value of the urban land use of the target point.
  • the second selection condition may be, for example, that the value of the artificial terrain classification is the same as the value of the artificial terrain classification of the target point.
  • the second selection condition may be, for example, that the value of the surface geology is the same as the value of the surface geology at the target point.
  • the second selection condition may be, for example, that at least one of the value of urban land use, the value of artificial topography classification, and the value of surface geology is common to the target point.
  • the second selection condition may be, for example, that the value of the building information in the geospatial information is the same as the value of the building information at the target point.
  • the selection condition may satisfy only the second selection condition.
  • the geospatial information acquisition unit 160 acquires the value of the geospatial information from the geospatial information storage device 300.
  • the geospatial information acquisition unit 160 receives, for example, a request for a value of the geospatial information from the estimation unit 120 for which the point information of the measurement point is specified (that is, the point information of the measurement point is included).
  • the geospatial information acquisition unit 160 may receive a plurality of requests for the values of the geospatial information including the point information of the plurality of measurement points from the estimation unit 120.
  • the geospatial information acquisition unit 160 reads out a predetermined type of geospatial information at a point represented by the point information specified in the received request from the geospatial information storage device 300.
  • the geospatial information acquisition unit 160 transmits the request for the value of the geospatial information requested by the received request to the geospatial information storage device 300. Then, the geospatial information acquisition unit 160 receives from the geospatial information storage device 300 a predetermined type of geospatial information at the point represented by the point information specified in the received request.
  • the predetermined type may be a predetermined type of all types of geospatial information.
  • the predetermined type may be all types of geospatial information.
  • the geospatial information acquisition unit 160 may similarly receive a request for a value of geospatial information at a target point (for example, a request for a value of geospatial information including the point information of a target point) from the determination unit 110. Similarly, the geospatial information acquisition unit 160 may read the value of the geospatial information at the point represented by the point information of the target point from the geospatial information storage device 300. The geospatial information acquisition unit 160 causes the determination unit 110 to determine the value of the geospatial information at the point represented by the point information of the target point, which is read from the geospatial information storage device 300, that is, received from the geospatial information storage device 300. Send out.
  • the output unit 140 receives an estimated value of the change in the height of the ground surface at the measurement point in the designated area from the estimation unit 120. Specifically, the output unit 140 receives from the estimation unit 120 the above-mentioned estimation data including the estimation value of the change in the height of the ground surface at the measurement point in the designated region. The output unit 140 outputs the received estimated value of the change in the height of the ground surface to the output destination device 400. The output unit 140 may output an estimated value of the change in the height of the ground surface at a predetermined measurement time point to the output destination device 400. The output unit 140 may output, for example, an estimated value of the change in the height of the ground surface at the time of measurement instructed by the user of the estimation device 100 to the output destination device 400.
  • the output unit 140 Even if the output unit 140 outputs the estimated value of the change in the height of the ground surface at each measurement time point from the oldest measurement time point to the newest measurement point point included in the estimation data to the output destination device 400 in order. good.
  • the output unit 140 may output the estimated data to the output destination device 400.
  • the output destination device 400 is a computer as described later, the output unit 140 outputs an estimated value of the change in the height of the ground surface at the measurement time point instructed by the output destination device 400 to the output destination device 400. May be good.
  • the output destination device 400 may be, for example, a display. In that case, the output unit 140 may display graphic or character information such as a graph showing an estimated value of the change in the height of the ground surface on the output destination device 400.
  • the output destination device 400 may be, for example, a storage device. In that case, the output unit 140 may store information representing an estimated value of the change in the height of the ground surface in the output destination device 400.
  • the output destination device 400 may be a computer such as a terminal device or a server device. In that case, the output unit 140 may transmit information representing an estimated value of the change in the height of the ground surface to the output destination device 400.
  • FIG. 2 is a flowchart showing an example of the operation of the estimation device 100 of the present embodiment.
  • the estimation unit 120 receives the measurement data acquired by the measurement data acquisition unit 150 from the measurement data storage device 200 (step S101).
  • the estimation unit 120 acquires the value of the geospatial information at the measurement point included in the designated range of the measurement data from the geospatial information storage device 300 via the geospatial information acquisition unit 160 (step S102).
  • the estimation unit 120 selects a target point from the measurement points included in the designated range of the measurement data (step S103).
  • the estimation unit 120 may select one of the measurement points that has not yet been selected as the target point as the target point.
  • the determination unit 110 determines the selection condition as described above (step S104).
  • the estimation unit 120 selects a measurement point that satisfies the selection condition (step S105).
  • the estimation unit 120 estimates the estimated value of the change transition data at the target point (that is, the above-mentioned estimated change transition data) from the measurement data at the selected measurement point (that is, the measurement point satisfying the selection condition) (step S106). ..
  • the target point is also included in the measurement point satisfying the selection condition. Therefore, in other words, the estimation unit 120 estimates the estimated value of the change transition data of the target point based on the measurement data of the target point and the measurement data of the measurement point satisfying the selection condition.
  • step S107 If there is a measurement point not selected as the target point in the measurement point included in the target range (NO in step S107), the operation of the estimation device 100 returns to step S103.
  • the output unit 140 transfers the estimated value of the change transition data (that is, the above-mentioned estimated data) to the output destination device 400. Output (step S108).
  • the estimation unit 120 is based on the measured value of the transition of the ground surface height of the target point among the measurement points and the measured value of the transition of the ground surface height of the measurement point satisfying the selection condition. This is because the estimated value of the change in the height of the ground surface at the target point is estimated. Then, the determination unit 110 determines the selection condition based on the geospatial information of the target point (for example, the value of the geospatial information of a predetermined type).
  • the estimation unit 120 may estimate the estimated value of the change in the height of the ground surface at the target point based on the weight based on the similarity of the geospatial information between the measurement point and the target point. For example, the estimation unit 120 may determine the weight according to the distance of the value of the geospatial information between the target point and the measurement point. The estimation unit 120 may determine the weight so that the shorter the distance between the values of the geospatial information, the smaller the weight value.
  • the estimation unit 120 sets the weight value when the value of the geospatial information of the measurement point is the same as the value of the geospatial information of the target point. It may be decided to 1.
  • the estimation unit 120 May determine the weight value to be greater than 0 and less than 1. The estimation unit 120 does not have to determine the weight of the measurement point that does not satisfy the selection condition.
  • the estimation unit 120 may determine the weight of the measurement point, which has a size corresponding to the distance between the measurement point and the target point. For example, the estimation unit 120 determines the weight for the measurement point from the distance between the measurement point and the target point by a function that converts the distance into a weight so that the magnitude relation of the distance and the magnitude relation of the weight are not reversed. You may.
  • the estimation unit 120 may determine the weight based on the similarity of the geospatial information between the measurement point and the target point and the distance between the measurement point and the target point. In this case, the estimation unit 120 sets the product of the weight based on the similarity of the geospatial information between the measurement point and the target point and the weight based on the distance between the measurement point and the target point at the measurement point. May be determined by the weight of.
  • the estimation unit 120 may calculate the product of the measured value of the change in the height of the ground surface at the measurement point satisfying the selection condition and the weight determined for the measurement point.
  • the estimation unit 120 may calculate the sum of the above-mentioned products calculated for all the measurement points satisfying the selection condition.
  • the estimation unit 120 may calculate the total weight of the measurement points that satisfy the selected scene. Then, the estimation unit 120 may use a value obtained by dividing the sum of the above-mentioned products by the sum of the above-mentioned weights as an estimated value of the change in the height of the ground surface at the target point.
  • the electromagnetic wave emitted from the antenna is specularly reflected by the rippling water surface and happens to be normal at the time of measurement.
  • the electromagnetic wave reflected at the point on the water surface where the direction and the direction of the electromagnetic wave match returns to the direction of the antenna. Since the shape of the water surface changes variously depending on the size and shape of the wave, it can be considered that the moving antenna does not receive electromagnetic waves from the same point at the same height at a plurality of time points. Therefore, when the measurement point is the water surface, it is difficult to obtain an accurate measurement value of the transition of the water surface height at the measurement point by the synthetic aperture radar.
  • the measured value of such a transition of the water surface height can be regarded as noise. Furthermore, in this disclosure, changes in the height of buildings and the ground are subject to measurement, but water surfaces are not.
  • the estimation unit 120 determines the height of the ground surface of the target point. It is not necessary to estimate the estimated value of change. In this case, the estimation unit 120 may set all the estimated values of changes in the height of the ground surface of the target point to 0. Further, it is not necessary to use the measured value of the height of the ground surface at the measurement point where the ground surface is the water surface to calculate the estimated value of the height of the ground surface at the target point.
  • the geospatial information may include a history of land changes.
  • the land change history may be represented by the value of at least one of the above geodata information at a plurality of past time points.
  • the estimation unit 120 may select a measurement point from the measurement points satisfying the selection condition based on the change history of the land. Specifically, the estimation unit 120 may select a measurement point whose land change history is the same as that of the land change history of the target point from the measurement points satisfying the selection condition.
  • the estimation unit 120 may estimate the estimated value of the change in the height of the ground surface at the target point based on the measured value of the change in the height of the ground surface at the measurement point selected in this way.
  • the estimation device 100 may include a measurement data storage device 200.
  • the estimation device 100 may include a geospatial information storage device 300.
  • the estimation device 100 may include a measurement data storage device 200 and a geospatial information storage device 300.
  • the measurement data stored in the measurement data storage device 200 may further include change transition data calculated in advance based on the height transition data. In that case, the estimation unit 120 does not have to calculate the change transition data.
  • a sixth modification of the first embodiment >> It is also possible to carry out a modification in which two or more of the above-mentioned first to fifth modifications are combined with the first embodiment.
  • FIG. 3 is a block diagram showing an example of the configuration of the estimation system 2 according to the present embodiment.
  • the estimation system 2 shown in FIG. 3 is the same as the estimation system 1 of the first embodiment, except for the differences described below.
  • the measurement data storage device 200 of the present embodiment is the same as the measurement data storage device 200 of the first embodiment.
  • the geospatial information storage device 300 of the present embodiment is the same as the geospatial information storage device 300 of the first embodiment.
  • the estimation system 2 includes the estimation device 101 instead of the estimation device 100.
  • the estimation device 101 includes a determination unit 130 in addition to each element of the estimation device 100 of the first embodiment.
  • the elements other than the determination unit 130 of the estimation device 101 of the present embodiment are the same as the elements of the estimation device 100 of the first embodiment to which the same name and reference numeral are given, except for the differences described below. be.
  • the selection condition determined by the determination unit 110 of the present embodiment may be that the value of the building information, which is a kind of geospatial information, is the same as the value of the building information at the target point.
  • the determination unit 110 of the present embodiment may determine the above-mentioned selection conditions.
  • the estimation unit 120 of the present embodiment operates in the same manner as the estimation unit 120 of the first embodiment.
  • the estimation unit 120 of the present embodiment further sends the above-mentioned estimation data to the determination unit 130.
  • the estimation unit 120 of the present embodiment may further send the value of the building information, which is a kind of geospatial information, to the determination unit 130 at the measurement point within the designated range.
  • the determination unit 130 receives the estimation data from the estimation unit 120.
  • the determination unit 130 may further receive the value of the building information at the measurement point within the designated range from the estimation unit 120.
  • the determination unit 130 receives the value of the building information, which is a kind of geospatial information, at the point indicated by the point information included in the estimation data from the geospatial information storage device 300 via the geospatial information acquisition unit 160. May be good.
  • the determination unit 130 may estimate the type of the building including the measurement point based on the estimated value of the transition of the change in the height of the ground surface.
  • the determination unit 130 may determine the type of building by, for example, a determination device that determines the magnitude of seasonal variation in height depending on the type of building, which has been learned in advance by machine learning. As the machine learning method in this case, any of the existing machine learning methods may be used.
  • the height fluctuation is not necessarily limited to the seasonal height fluctuation.
  • the determination unit 130 may determine the type of the building at each of the measurement points that are the buildings included in the designated area.
  • the determination unit 130 may set information indicating that the building is not a building as the type of the building at the measurement point that is not the building.
  • the determination unit 130 may acquire the building map of the designated area from the geospatial information storage device 300 via the geospatial information acquisition unit 160.
  • a building map is stored in advance in the geospatial information storage device 300.
  • the determination unit 130 may send the determined information indicating the type of the building at the measurement point included in the designated area to the output unit 140.
  • the determination unit 130 may further send information representing the acquired building map to the output unit 140.
  • the output unit 140 of the present embodiment has the same function as the output unit 140 of the first embodiment.
  • the output unit 140 of the present embodiment further receives information from the determination unit 130 indicating the type of the building at the measurement point included in the designated area.
  • the output unit 140 of the present embodiment may further receive information representing a building map from the determination unit 130.
  • the output unit 140 outputs the received information indicating the type of the building at the measurement point included in the designated area to the output destination device 400.
  • the output unit 140 may superimpose a mark indicating the type of the building at the measurement point on the position of the building map corresponding to the position of the measurement point.
  • the output unit 140 may output a building map on which a mark indicating the type of the building is superimposed to the output destination device 400.
  • the output unit 140 may determine the type of building included in the building map based on the estimated type of building at the measurement point included in the area of the building.
  • the output unit 140 may change at least one of the patterns and colors of the area of the building included in the building map according to the determined type of building. In this case, at least one of the patterns and colors of the area of the building may be predetermined for each type of building.
  • the method for determining the type of building may be a predetermined method. For example, when the building contains only one measurement point, the output unit 140 may determine the type of the building as the type of the building determined for the measurement point. For example, when a building contains a plurality of measurement points, the output unit 140 may determine the type of the building as the type of the building determined for the most measurement points.
  • the output unit 140 may be used. It is not necessary to change the pattern and color of the area of the building.
  • FIG. 4 is a flowchart showing an example of the operation of the estimation device 101 of the present embodiment.
  • the operation of the estimation device 101 of the present embodiment from step S101 to step S107 is the same as the operation of the estimation device 100 of the first embodiment from step S101 to step S107.
  • the estimation device 101 of the present embodiment performs the operation described below after the operation of step S107.
  • the determination unit 130 determines the type of the building including the measurement point in the area based on the above-mentioned estimation data of the measurement point (step S208).
  • the determination unit 130 may determine the type of the building for each of the measurement points included in the area of any of the buildings.
  • the output unit 140 outputs the estimated value of the measurement data (that is, the above-mentioned estimated data) and the type of the building (step S209).
  • This embodiment has the effect of being able to know the type of building.
  • the reason is that the determination unit 130 determines the type of the building based on the estimation data of the measurement points included in the area of the building.
  • the second embodiment can also be modified as in the first to fifth modifications of the first embodiment.
  • a second modification of the second embodiment >> The output unit 140 of this modification outputs the received information indicating the type of the building at the measurement point included in the designated area to the output destination device 400, and does not have to output the estimation data.
  • the estimation unit 120 of this modification does not have to send the estimation data to the output unit 140.
  • the output unit 140 of this modification does not have to receive the estimation data from the estimation unit 120.
  • FIG. 5 is a block diagram showing an example of the configuration of the estimation device 102 according to the present embodiment.
  • the estimation device 102 includes a determination unit 110 and an estimation unit 120.
  • the determination unit 110 determines the selection condition based on the geospatial information of the target point among the plurality of measurement points whose transition of the height of the ground surface is measured by the synthetic aperture radar.
  • the determination unit 110 may determine the selection conditions in the same manner as the determination unit 110 of the first embodiment.
  • the estimation unit 120 determines the ground at the target point based on the measured value of the transition of the ground surface height at the target point and the measured value of the transition of the ground surface height at the measurement point satisfying the selection condition. Estimate an estimate of changes in surface height.
  • the estimation unit 120 may estimate the estimated value of the change in the height of the ground surface at the target point, as in the estimation unit 120 of the first embodiment.
  • FIG. 6 is a flowchart showing an example of the operation of the estimation device 102 according to the present embodiment.
  • the determination unit 110 determines the selection condition based on the geospatial information of the target point (step S301).
  • the estimation unit 120 determines the ground surface of the target point based on the measured value of the transition of the ground surface height at the target point and the measured value of the transition of the ground surface height at the measurement point satisfying the selection condition. Estimate the estimated value of the change in height of (step S302).
  • the above-mentioned estimation device 100, estimation device 101, estimation device 102, measurement data storage device 200, and geospatial information storage device 300 can be realized by a computer, respectively.
  • the computer includes a memory loaded with a program read from the storage medium and a processor executing the program.
  • the computer may be realized by a combination of a plurality of computers connected so as to be able to communicate with each other.
  • the estimation device 100, the estimation device 101, the estimation device 102, the measurement data storage device 200, and the geospatial information storage device 300 can each be realized by dedicated hardware. Dedicated hardware may be implemented, for example, by a single circuit or by a plurality of circuits communicably connected to each other.
  • the estimation device 100, the estimation device 101, the estimation device 102, the measurement data storage device 200, and the geospatial information storage device 300 can also be realized by the combination of the above-mentioned computer and dedicated hardware.
  • FIG. 7 is a diagram showing an example of a hardware configuration of a computer 1000 that can realize an apparatus such as an estimation apparatus according to the embodiment of the present disclosure.
  • the computer 1000 includes a processor 1001, a memory 1002, a storage device 1003, and an I / O (Output / Output) interface 1004.
  • the computer 1000 can access the storage medium 1005.
  • the memory 1002 and the storage device 1003 are storage devices such as a RAM (Random Access Memory) and a hard disk, for example.
  • the storage medium 1005 is, for example, a storage device such as a RAM or a hard disk, a ROM (Read Only Memory), or a portable storage medium.
  • the storage device 1003 may be a storage medium 1005.
  • the processor 1001 can read and write data and programs to the memory 1002 and the storage device 1003. Processor 1001 can access other devices, for example, via the I / O interface 1004. The processor 1001 can access the storage medium 1005.
  • the storage medium 1005 stores a program for operating the computer 1000 as the device according to the above-described embodiment.
  • the processor 1001 loads the memory 1002 with a program stored in the storage medium 1005 that causes the computer 1000 to operate as the estimation device according to the above-described embodiment. Then, the processor 1001 executes the program loaded in the memory 1002, so that the computer 1000 operates as the estimation device according to the above-described embodiment.
  • the processor 1001 loads the memory 1002 with a program stored in the storage medium 1005 for operating the computer 1000 as the measurement data storage device according to the above-described embodiment. Then, when the processor 1001 executes the program loaded in the memory 1002, the computer 1000 operates as the measurement data storage device according to the above-described embodiment.
  • the processor 1001 loads the memory 1002 with a program stored in the storage medium 1005 that causes the computer 1000 to operate as the geospatial information storage device according to the above-described embodiment. Then, when the processor 1001 executes the program loaded in the memory 1002, the computer 1000 operates as the geospatial information storage device according to the above-described embodiment.
  • the determination unit 110, the estimation unit 120, the determination unit 130, the output unit 140, the measurement data acquisition unit 150, and the geospatial information acquisition unit 160 are realized by, for example, a processor 1001 that executes a program loaded in the memory 1002. Can be done.
  • the first input / output unit 220 and the second input / output unit 320 can be realized by, for example, a processor 1001 that executes a program loaded in the memory 1002.
  • the measurement data storage unit 210 and the geospatial information storage unit 310 can be realized by a storage device 1003 such as a memory 1002 or a hard disk device included in the computer 1000.
  • the height of the ground surface at the target point is based on the measured value of the transition of the height of the ground surface at the target point and the measured value of the transition of the height of the ground surface at the measurement point satisfying the selection condition.
  • Appendix 2 The estimation device according to Appendix 1, wherein the determination means determines the selection condition based on the type of the ground surface based on the geospatial information of the target point and the distance from the target point.
  • Appendix 4 The estimation device according to Appendix 3, wherein the estimation means estimates the estimated value based on the weight further based on the distance between the measurement point and the target point.
  • the estimation means further selects the measurement point from the measurement point satisfying the selection condition based on the change history of the land at the measurement point, and estimates the estimated value based on the measurement value of the measurement point.
  • the estimation device according to any one of Supplementary note 1 to 4.
  • the estimation means does not estimate the estimated value at the target point when the ground surface of the target point is the water surface, and when the ground surface of the measurement point is the water surface, the height fluctuation of the measurement point changes.
  • the estimation device according to any one of Supplementary note 1 to 5, wherein the measured value is not used for estimating the estimated value at the target point.
  • the determination means determines the selection condition so as to select the measurement point included in the area of the building when the target point is included in the area of the building in the geospatial information.
  • the estimation device according to any one item.
  • any of the appendices 1 to 7 further provided with a determination means for determining the type of the building including the target point based on the characteristics of the estimated value.
  • the estimation device according to item 1.
  • the selection conditions are determined based on the geospatial information of the target point among multiple measurement points where the transition of the height of the ground surface is measured by the synthetic aperture radar.
  • the height of the ground surface at the target point is based on the measured value of the transition of the height of the ground surface at the target point and the measured value of the transition of the height of the ground surface at the measurement point satisfying the selection condition. Estimate the estimated value of change, Estimating method.
  • Appendix 10 The estimation method according to Appendix 9, which determines the selection condition based on the type of the ground surface based on the geospatial information of the target point and the distance from the target point.
  • Appendix 12 The estimation method according to Appendix 11, which estimates the estimated value based on the weight further based on the distance between the measurement point and the target point.
  • the measurement points are further selected based on the change history of the land at the measurement points, and the estimated values are estimated based on the measurement values at the measurement points.
  • the estimation method according to any one item.
  • (Appendix 17) A decision process that determines selection conditions based on the geospatial information of the target point among multiple measurement points where the transition of the height of the ground surface was measured by the synthetic aperture radar.
  • the height of the ground surface at the target point is based on the measured value of the transition of the height of the ground surface at the target point and the measured value of the transition of the height of the ground surface at the measurement point satisfying the selection condition.
  • Estimating processing to estimate the estimated value of change and A storage medium that stores a program that causes a computer to execute a program.
  • Appendix 18 The storage medium according to Appendix 17, wherein the determination process determines the selection condition based on the type of the ground surface based on the geospatial information of the target point and the distance from the target point.
  • Appendix 20 The storage medium according to Appendix 19, wherein the estimation process estimates the estimated value based on the weight further based on the distance between the measurement point and the target point.
  • the measurement point is further selected from the measurement point satisfying the selection condition based on the change history of the land at the measurement point, and the estimated value is estimated based on the measurement value of the measurement point.
  • the storage medium according to any one of Supplementary note 17 to 20.
  • the estimation process does not estimate the estimated value at the target point when the ground surface of the target point is the water surface, and when the ground surface of the measurement point is the water surface, the height fluctuation of the measurement point changes.
  • the storage medium according to any one of Supplementary note 17 to 21, wherein the measured value is not used for estimating the estimated value at the target point.
  • Estimating system 2 Estimating system 100 Estimating device 101 Estimating device 102 Estimating device 110 Determining unit 120 Estimating unit 130 Judgment unit 140 Output unit 150 Measurement data acquisition unit 160 Geospatial information acquisition unit 200 Measurement data storage device 210 Measurement data storage unit 220 1 Input / output unit 300 Geospatial information storage device 310 Geospatial information storage unit 320 Second input / output unit 400 Output destination device 1000 Computer 1001 Processor 1002 Memory 1003 Storage device 1004 I / O interface 1005 Storage medium

Landscapes

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

Abstract

測定値の変化に対象の領域に応じたばらつきがある場合であっても、変化の推定の精度を向上させることができる推定装置等を提供する。本開示の一態様に係る推定方法は、合成開口レーダーによって高さの推移が測定された複数の測定地点のうち対象地点の地理空間情報に基づいて、選択条件を決定し、前記対象地点における高さの推移の測定値と、前記選択条件を満たす前記測定地点における高さの推移の測定値とに基づいて、前記対象地点における高さの変化の推定値を推定する。

Description

推定装置、推定方法及び記憶媒体
　本開示は、測定の結果からノイズ成分が除去された推定値を推定する技術に関し、特に、合成開口レーダーによって測定された結果から推定値を推定する技術に関する。
　例えば合成開口レーダーによって測定された地表面の高さの変動などの、センサによって測定された測定値には、一般に、ノイズが含まれる。
　測定の結果からノイズが除去された推定値を得る技術の１つが、特許文献１によって開示されている。特許文献１の技術では、同一の画像が得られるように配置された複数の赤外線撮像素子によって得られた画像信号の、同一の位置の画素の信号を加算平均する装置が記載されている。
　特許文献２には、レーザ光を反射する計測点までの距離を計測するレーザスキャニング装置の計測の結果として得られた、三次元座標が得られた点の集合である三次元点群データからノイズ除去を行う方法が開示されている。その方法は、三次元点群データに対して設定された二次元の矩形領域に対応する計測点により形成される曲面を、計測点を頂点とする三角形の集合に分割する。その方法は、全ての三角形の法線ベクトルを一つにまとめた合成法線ベクトルを求める。そして、その方法は、合成法線ベクトルが座標軸に一致するよう座標を変換した場合の座標軸方向の測定値の標準偏差が所定値を超える場合、それらの測定点の１つをノイズとして取り扱う。
特開平９－２３３３９０号公報 特開２００５－０２４３７０号公報
　特許文献１の技術では、同時に複数のセンサにより撮像された複数の画像の同一の位置の画素の信号が加算平均される。また、特許文献２の技術は、同じタイミングで距離が計測された複数の計測点において、ノイズを特定する技術である。特許文献２の技術では、所定の計測点の計測値にノイズが含まれている場合に、その計測点の三次元座標を推定することはできない。さらに、特許文献１及び２の技術では、計測値等の変化を推定することはできない。従って、特許文献１及び２の技術を、例えば地表面の高さの推定に応用することはできない。特許文献１及び２の技術では、測定の対象の領域の性質の違いに応じて、対象の測定値の変化にばらつきが生じる場合に、精度のよい測定値の変化の推定ができない。
　本開示の目的の１つは、測定値の変化に対象の領域に応じたばらつきがある場合であっても、変化の推定の精度を向上させることができる推定装置等を提供することにある。
　本開示の一態様に係る推定装置は、合成開口レーダーによって地表面の高さの推移が測定された複数の測定地点のうち対象地点の地理空間情報に基づいて、選択条件を決定する決定手段と、前記対象地点における地表面の高さの推移の測定値と、前記選択条件を満たす前記測定地点における地表面の高さの推移の測定値とに基づいて、前記対象地点における地表面の高さの変化の推定値を推定する推定手段と、を備える。
　本開示の一態様に係る推定方法は、合成開口レーダーによって地表面の高さの推移が測定された複数の測定地点のうち対象地点の地理空間情報に基づいて、選択条件を決定し、前記対象地点における地表面の高さの推移の測定値と、前記選択条件を満たす前記測定地点における地表面の高さの推移の測定値とに基づいて、前記対象地点における地表面の高さの変化の推定値を推定する。
　本開示の一態様に係る記憶媒体は、合成開口レーダーによって地表面の高さの推移が測定された複数の測定地点のうち対象地点の地理空間情報に基づいて、選択条件を決定する決定処理と、前記対象地点における地表面の高さの推移の測定値と、前記選択条件を満たす前記測定地点における地表面の高さの推移の測定値とに基づいて、前記対象地点における地表面の高さの変化の推定値を推定する推定処理と、をコンピュータに実行させるプログラムを記憶する。本開示の一態様は、上述のプログラムによっても実現される。
　本開示には、測定値の変化に対象の領域に応じたばらつきがある場合であっても、変化の推定の精度を向上させることができるという効果がある。
図１は、本開示の第１の実施形態に係る推定システムの構成の例を表すブロック図である。 図２は、本開示の第１の実施形態に係る推定装置の動作の例を表すフローチャートである。 図３は、本開示の第２の実施形態に係る推定システムの構成の例を表すブロック図である。 図４は、本開示の第２の実施形態に係る推定装置の動作の例を表すフローチャートである。 図５は、本開示の第３の実施形態に係る推定装置の構成の例を表すブロック図である。 図６は、本開示の第３の実施形態に係る推定装置の動作の例を表すフローチャートである。 図７は、本開示の実施形態に係る推定装置などの装置を実現できるコンピュータの構成の例を表すブロック図である。
　次に、本開示の実施の形態を、図面を参照しながら詳細に説明する。
　＜第１の実施形態＞
　＜＜構成＞＞
　図１は、本開示の第１の実施形態に係る推定システム１の構成の例を表すブロック図である。図１に示す例では、推定システム１は、推定装置１００と、測定データ記憶装置２００と、地理空間情報記憶装置３００と、出力先装置４００とを含む。
　推定装置１００は、決定部１１０と、推定部１２０と、出力部１４０と、測定データ取得部１５０と、地理空間情報取得部１６０とを含む。測定データ記憶装置２００は、測定データ記憶部２１０と、第１入出力部２２０とを含む。地理空間情報記憶装置３００は、地理空間情報記憶部３１０と、第２入出力部３２０を含む。
　推定装置１００の決定部１１０は、合成開口レーダーによって地表面の高さの推移が測定された複数の測定地点のうち対象地点の、地理空間情報記憶装置３００に格納されている地理空間情報に基づいて、選択条件を決定する。地理空間情報は、例えば、地表面の状態及びその地表面の地下の状態の少なくとも一方を表す情報である。地理空間情報は、例えば、測定地点の地質や建造物の有無に関連する。推定部１２０は、対象地点における地表面の高さの推移の測定値と、選択条件を満たす測定地点における地表面の高さの推移の測定値とに基づいて、対象地点における地表面の高さの変化の推定値を推定する。したがって、測定の対象の領域の地質や建造物の有無の違いに応じて、対象の測定値の変化にばらつきが生じる場合であっても、精度のよい測定値の変化の推定ができる。
　測定地点における地表面の高さの推移の測定値は、測定データとして、測定データ記憶装置２００に格納されている。出力部１４０は、対象地点における地表面の高さの変化の推定値を、出力先装置４００に出力する。
　以下、詳細に説明する。
　＜＜測定データ記憶部２１０＞＞
　測定データ記憶部２１０は、測定データを記憶する。
　測定データは、例えば、合成開口レーダーによる、複数の地点における、地表面の高さの測定の結果の推移（以下、地表面の高さの推移と表記）を表すデータである。地表面の高さは、地表面の、あらかじめ定められた基準面からの高さであってよい。測定データの地表面の高さは、例えば、人工衛星や航空機などの飛翔体に搭載されたレーダーを合成開口レーダー（ＳＡＲ；Ｓｙｎｔｈｅｔｉｃ　Ａｐｅｒｔｕｒｅ　Ｒａｄａｒ）として使用する観測によって得られた、地表面上の複数の測定地点における高さである。以下の説明では、このような観測を、合成開口レーダーによる観測と表記する。地表面の高さの推移を表すデータは、例えば、高さを表す値と、その高さが観測によって得られた時点を表すデータとの組み合わせを複数含むデータであってよい。時点を表すデータの単位は、適宜定められていてよい。例えば、時点を表すデータは、日付を表していてもよく、日付及び時刻を表していてもよい。時刻の単位も、適宜定められていてよい。
　地表面において、地表面の高さの推移が得られている地点を、測定地点と表記する。測定データは、複数の測定地点の、おのおのの測定地点の位置を表す位置データと、その測定地点における地表面の高さの推移を表すデータと、の組み合わせを含むデータであってよい。地表面上の地点（すなわち測定地点）の位置を表す情報は、例えば、緯度及び経度の情報であってよい。位置を表す情報は、地表面上の位置を特定できる他の情報であってもよい。以下では、地点の位置を表す情報を、地点情報と表記する。
　地表面上の測定地点における、地表面の高さの推移を表すデータを、以下では、高さ推移データと表記する。測定データに含まれる高さ推移データは、ノイズの除去などが行われていないデータである。測定データの高さ推移データには、ノイズなどが含まれている。
　言い換えると、測定データにおいて、高さ推移データは、その高さ推移データが表す高さの推移が測定された測定地点の地点情報と関連付けられていてよい。さらに言い換えると、測定データは、複数の測定地点の地点情報と、複数の測定地点の各々における高さ推移データとを含む。
　＜＜第１入出力部２２０＞＞
　第１入出力部２２０は、測定データ記憶部２１０に格納されている測定データを読み出し、推定装置１００に供給する。第１入出力部２２０は、例えば、推定装置１００から測定データの要求を受け取ってよい。測定データの要求は、地表面における範囲を特定する情報を含んでいてよい。第１入出力部２２０は、受け取った要求に含まれる測定地点の測定データを、測定データ記憶部２１０から読み出し、読み出した測定データを推定装置１００に供給してよい。
　＜＜地理空間情報記憶部３１０＞＞
　地理空間情報記憶部３１０は、地理空間情報を記憶する。
　前述したように、本実施形態では、地理空間情報は、例えば、地表面の状態及びその地表面の地下の状態の少なくとも一方を表す情報である。地理空間情報は、いわゆる地理情報システム（Ｇｅｏｇｒａｐｈｉｃ　Ｉｎｆｏｒｍａｔｉｏｎ　Ｓｙｓｔｅｍ）から得られる情報の少なくともいずれかであってもよい。地理空間情報は、人工衛星や航空機などから観測によって得られたデータであってもよい。地理空間情報は、現地調査によって得られたデータであってもよい。地理空間情報は、計測や調査によって得られたデータに基づく分析の結果を表す情報であってもよい。地理空間情報は、計測や調査によって得られたデータに基づいて、人為的に定められた情報であってもよい。地理空間情報は、ＧＩＳ（Ｇｅｏｇｒａｐｈｉｃ　Ｉｎｆｏｒｍａｔｉｏｎ　Ｓｙｓｔｅｍ）データと表記されることもある。
　地理空間情報は、あらかじめ、地理情報システムから取得され、地理空間情報記憶部３１０に格納されていてよい。複数の種類の地理空間情報が、地理空間情報記憶部３１０に格納されていてよい。地理空間情報は、地点情報（例えば緯度及び経度）によって特定される地点の、地理空間情報の値を特定できる形式で表されていてよい。地理空間情報の具体例については、後で詳細に説明する。
　＜＜第２入出力部３２０＞＞
　第２入出力部３２０は、地理空間情報記憶部３１０に格納されている地理空間情報から地理空間情報の値を抽出し、抽出した地理空間情報の値を推定装置１００に供給する。第２入出力部３２０は、例えば、地点情報を含む、地理空間情報の値の要求を推定装置１００から受け取ってよい。第２入出力部３２０は、受け取った要求の地点情報によって特定される地点の地理空間情報の値を、地理空間情報記憶部３１０に格納されている地理空間情報から抽出してよい。地理空間情報の値の要求は、複数の地点の地点情報を含んでいてよい。第２入出力部３２０は、地理空間情報の値の要求に含まれる、複数の地点の地点情報によって特定される地点の地理空間情報の値を、地理空間情報記憶部３１０に格納されている地理空間情報から抽出してよい。第２入出力部３２０は、抽出した、複数の地点の地理空間情報の値を、推定装置１００に供給してよい。
　第２入出力部３２０は、あらかじめ定められている種類の地理空間情報の、地点情報によって特定される地点における値を抽出してよい。第２入出力部３２０は、地理空間情報記憶部３１０に格納されている全ての種類の地理空間情報の、地点情報によって特定される地点における値を抽出してもよい。地点情報によって特定される地点における値が存在しない地理空間情報が存在する場合、第２入出力部３２０は、その地理空間情報の値を抽出しなくてもよい。地点情報によって特定される地点における値が存在しない地理空間情報が存在する場合、第２入出力部３２０は、その地理空間情報の値を、値が存在しないことを表す値（例えば、０等）に設定してもよい。
　＜＜地理空間情報の詳細＞＞
　地理空間情報は、例えば、地表面が区切られたメッシュごとの、状態を表す値によって表されていてもよい。この場合、第２入出力部３２０は、地点情報によって特定される位置が含まれるメッシュの中における状態を表す、地理空間情報の値を、地点情報によって特定される地点の地理空間情報の値として抽出する。メッシュのサイズや形状は、地理空間情報の種類ごとに定められていてよい。
　地理空間情報は、他の形式で表されていてもよい。地理空間情報は、例えば、状態が異なる領域の間の境界線と、境界線によって区切られる領域内における状態を表す値とによって表されていてよい。この場合、第２入出力部３２０は、地点情報によって特定される位置が含まれる領域内における状態を表す値を、地点情報によって特定される地点の地理空間情報の値として抽出する。地理空間情報の形式は、地理空間情報の種類ごとに定められていてよい。
　具体的な地理空間情報は、例えば、盛土造成地のタイプ、平均傾斜角度、平均降水量（例えば、平均年間降水量）、表層地質、急傾斜地指定、土砂災害警戒区域指定、液状化危険度、雨水浸透桝設置可否、地震時の揺れ易さ、排水困難低地、都市域土地利用、自然地形分類、人工地形分類、表層地質、河川敷、施設情報（工事の有無等）であってよい。本実施形態の地理空間情報は、必ずしもこれらの例に限られない。本実施形態の地理空間情報は、これらの例の一部を含んでいなくてもよい。
　盛土造成地のタイプは、盛土がなされた地表面の形状によって定まる、土の盛り方を表していてよい。盛土造成地のタイプは、例えば、谷や沢を盛土によって埋めた盛土である「谷埋め型盛土」や、傾斜地に行われた盛土である「腹付け型盛土」である。
　盛土造成地のタイプは、さらに、盛土の規模を表していてもよい。この場合、例えば、谷や沢を盛土によって埋めた盛土のうち、基準を満たす盛土（例えば、面積が３０００平方メートル以上の盛土）のタイプは、「大規模谷埋め型盛土」であってよい。この場合、谷や沢を盛土によって埋めた盛土のうち、基準を満たさない盛土のタイプが、「谷埋め型盛土」であってよい。また、傾斜地に行われた盛土のうち、基準を満たす盛土（例えば、盛土をする前の地盤面（原地盤面とも呼ばれる）の水平面に対する角度が２０度以上で、かつ、盛土の高さが５メートル以上である盛土）のタイプは、大規模腹付け型盛土であってよい。この場合、傾斜地に行われた盛土のうち、基準を満たさない盛土のタイプが、「腹付け型盛土」であってよい。
　盛土造成地のタイプの値は、例えば、「谷埋め型盛土」や「腹付け型盛土」などにあらかじめそれぞれ適宜割り振られた、互いに異なる数値のいずれか１つであってもよい。
　平均傾斜角度は、例えば、メッシュ単位で算出された、地表面の平均の傾斜角度のデータであってよい。平均傾斜角度の値は、算出された地表面の平均の傾斜角度であってよい。
　平均降水量は、例えば、メッシュ単位で算出された、地表面の平均の降水量のデータであってよい。平均降水量の値は、算出された地表面の平均の降水量であってよい。
　表層地質は、地表面の表層の地質（言い換えると、地質の種類）を表すデータであってよい。地質の種類は、あらかじめ定められていてよい。地質の種類の各々には、あらかじめ互いに異なる数値が割り当てられていてよい。表層地質の値は、地質にあらかじめそれぞれ適宜割り当てられている数値のいずれか１つであってよい。
　急傾斜地指定は、例えば自治体などによって、急傾斜地として指定されているか否かを表すデータであってよい。急傾斜地指定の値は、例えば、急傾斜地として指定されていることを表す数値、または、急傾斜地として指定されていないことを表す数値であってよい。これらの数値として、互いに異なる数値があらかじめ適宜決められていてよい。
　土砂災害警戒区域指定は、例えば自治体などによって、土砂災害警戒区域として指定されているか否かを表していてよい。土砂災害警戒区域指定の値は、例えば、土砂災害警戒区域として指定されていることを表す数値、又は、土砂災害警戒区域として指定されていないことを表す数値であってよい。これらの数値として、互いに異なる数値があらかじめ決められていてよい。
　液状化危険度は、例えば、土地が液状化する危険性の程度を表すデータであってよい。液状化危険度の値は、土地が液状化する危険性の程度を表す数値であってよい。液状化危険度の値は、それぞれ異なる程度を表す異なる複数の数値のいずれか１つであってよい。危険性の程度を表す数値は、あらかじめ適宜定められていてよい。
　雨水浸透桝設置可否は、例えば、地形、土質、地下水位に基づいて浸透施設の設置が可能か否かについて判断された結果を表す「浸透施設設置判断マップ」に基づく設置可否を表す情報であってよい。雨水浸透桝設置可否の値は、設置が可能であることを表す数値、又は、設置が可能でないことを表す数値であってよい。これらの数値として、互いに異なる数値があらかじめ適宜決められていてよい。
　地震時の揺れ易さは、例えば、地震が発生した場合の地表面の揺れ易さの程度を表すデータであってよい。地震時の揺れ易さの値は、地震が発生した場合の地表面の揺れ易さの程度を表す数値であってよい。地震時の揺れ易さの値は、地震が発生した場合の地表面の揺れ易さの程度を表す、複数の数値のいずれか１つであってよい。地表面の揺れ易さの程度を表す数値は、あらかじめ適宜定められていてよい。
　排水困難低地は、例えば、土地の標高や周囲との標高差等によって推定された、その土地が排水困難低地であるか否かを表していてよい。排水困難低地の値は、排水困難低地であることを表す数値、又は、排水困難低地でないことを表す数値であってよい。これらの数値は、あらかじめ適宜定められていてよい。
　都市域土地利用は、都市として指定されている地域における土地利用の種別であってよい。都市域土地利用における、土地利用の種別は、例えば衛星写真から判読されてもよい。あらかじめ定められた複数の種別から選択された土地利用の種別が、都市の地域に含まれる領域に対して設定されてよい。あらかじめ定められた複数の種別には、それぞれ、異なる数値があらかじめ適宜割り当てられていてよい。領域に対して設定される土地利用の種別の値は、その種別に割り当てられている数値であってよい。
　自然地形分類は、例えば、人間によって建造された建造物ではない場所における、地形の種別であってよい。自然地形分類として設定可能な、複数の地形の種別が、あらかじめ適宜定められていてよい。自然地形分類が設定される領域に対して、自然地形分類として設定可能な種別としてあらかじめ定められている複数の地形の種別から選択された地形の種別が設定されてよい。それらの複数の種別に対して、それぞれ異なる数値が割り当てられていてよい。領域に設定される自然地形分類における地形の値は、その領域に設定されている地形の種別に割り当てられている数値であってよい。
　人工地形分類は、例えば、人間が地形を改変した場所や人間によって建造された建造物である場所における、地形の種別であってよい。人工地形分類として設定可能な、複数の地形の種別が、あらかじめ適宜定められていてよい。人工地形分類が設定される領域に対して、人工地形分類として設定可能な種別としてあらかじめ定められている複数の地形の種別から選択された地形の種別が設定されてよい。それらの複数の種別に対して、それぞれ異なる数値が割り当てられていてよい。領域に設定される人工地形分類における地形の値は、その領域に設定されている地形の種別に割り当てられている数値であってよい。
　表層地質は、例えば、地表における土壌の種類を表していてよい。複数の土壌の種類が、あらかじめ適宜定められていてよい。そして、あらかじめ適宜定められている異なる数値が、それらの複数の土壌の種類に、それぞれ割り当てられていてよい。例えば調査の結果に基づく土壌の種類が、領域に対して設定されてよい。領域の表層地質の値は、その領域に設定されている土壌の種類に割り当てられている数値であってよい。
　具体的な地理空間情報は、建造物情報であってもよい。建造物情報は、例えば、個々の建造物の範囲を表していてよい。建造物情報の値は、例えば、個々の建造物に互いに異なる識別番号が割り当てられていてよい。建造物の範囲に含まれない地点の、建造物情報の値は、建造物ではないものを表す値（例えば、０等）であってもよい。
　＜＜測定データ取得部１５０＞＞
　測定データ取得部１５０は、測定データ記憶装置２００から、測定データを取得する。測定データ取得部１５０は、例えば、指定された範囲（以下、指定範囲と表記）の測定データの要求を、測定データ記憶装置２００に送信してよい。指定範囲は、あらかじめ定められていてもよい。指定範囲は、推定装置１００のユーザによって、例えば推定装置１００のキーボードやマウスなどの入力部によって、指定されてもよい。入力部は、図１では省略されている。測定データ取得部１５０は、指定された範囲の測定地点における測定データを、測定データ記憶装置２００から受け取る。
　測定データ取得部１５０は、受け取った測定データを、推定部１２０に送出する。
　＜＜推定部１２０＞＞
　推定部１２０は、指定領域の測定データを測定データ取得部１５０から受け取る。
　推定部１２０は、受け取った測定データにおいて、測定によって高さ推移データが得られている測定地点から、１つの測定地点を、対象地点として選択する。
　推定部１２０は、対象地点の地点情報を、決定部１１０に送出する。決定部１１０については後で詳細に説明する。推定部１２０は、対象地点の地点情報に基づいて決定された、測定地点を選択するための選択条件を、決定部１１０から受け取る。言い換えると、推定部１２０は、決定された選択条件を表す情報を、決定部１１０から受け取る。選択条件は、決定部１１０によって決定される。選択条件は、例えば、対象地点が選択条件を満たすように決定される。選択条件は、さらに、対象地点における地理空間情報の値に基づいて決定される。地理空間情報の値に関連する地理空間情報の種類は、あらかじめ定められていてよい。
　推定部１２０は、決定部１１０に選択条件を決定させるために地理空間情報の値を決定部１１０に送出する。この場合、推定部１２０は、対象範囲に含まれる測定地点の各々における地理空間情報の値を、地理空間情報取得部１６０を介して、地理空間情報記憶装置３００から取得する。推定部１２０は、例えば、地理空間情報の種類を特定する情報と測定地点の地点情報とが指定された、地理空間情報の要求を、地理空間情報取得部１６０を介して、地理空間情報記憶装置３００に送信する。言い換えると、推定部１２０は、地理空間情報の種類を特定する情報と測定地点の地点情報とを含む、地理空間情報の要求を、地理空間情報取得部１６０を介して地理空間情報記憶装置３００に送信する。推定部１２０は、地理空間情報記憶装置３００から、地理空間情報取得部１６０を介して、測定地点における地理空間情報の値を受け取る。具体的には、推定部１２０は、地理空間情報の要求を地理空間情報取得部１６０に送出してよい。地理空間情報取得部１６０は、地理空間情報の要求を受け取り、受け取った地理空間情報の要求を、地理空間情報記憶装置３００に送信する。地理空間情報取得部１６０は、地理空間情報記憶装置３００から、地理空間情報の値を受け取り、受け取った地理空間情報の値を推定部１２０に送出する。推定部１２０は、地理空間情報の値を地理空間情報取得部１６０から受信し、決定部１１０に出力する。
　推定部１２０が、対象地点の地点情報と共に、対象地点の地理空間情報の値を、決定部１１０に送出してもよい。決定部１１０は、対象地点の地点情報と、対象地点の地理空間情報の値とを、推定部１２０から受け取ってもよい。
　また、推定部１２０が、決定部１１０に地理空間情報の値を送出する代わりに、決定部１１０が、受け取った対象地点の地理空間情報が示す地点における地理空間情報の値を、地理空間情報取得部１６０を介して、地理空間情報記憶装置３００から取得してよい。
　具体的には、決定部１１０が、地理空間情報の種類を特定する情報と対象地点の地点情報とが指定された、地理空間情報の要求を、地理空間情報取得部１６０に送信する。言い換えると、決定部１１０が、地理空間情報の種類を特定する情報と対象地点の地点情報とを含む、地理空間情報の要求を、地理空間情報取得部１６０に送信する。そして、地理空間情報取得部１６０が、受け取った地理空間情報の要求を、地理空間情報記憶装置３００に送信する。地理空間情報取得部１６０は、対象地点における地理空間情報の値を、地理空間情報記憶装置３００から受け取り、受け取った、対象地点における地理空間情報の値を、決定部１１０に送出する。そして、決定部１１０は、対象地点における地理空間情報の値を地理空間情報取得部１６０から受け取る。
　決定部１１０は、対象地点における地理空間情報の値に基づく選択条件を決定する。推定部１２０は、決定された選択条件を表す情報を、決定部１１０から受け取る。
　選択条件については、後で詳細に説明する。
　推定部１２０は、受け取った測定データにおいて、高さ推移データが得られている測定地点から、選択条件を満たす測定地点を選択する。推定部１２０は、選択した測定地点の高さ推移データと対象地点の高さ推移データとに基づいて、対象地点の推定変化推移データを推定する。推定変化推移データについては、後で詳細に説明する。対象地点が選択条件を満たすように、選択条件が決定されている場合、選択された測定地点に対象地点が含まれている。推定部１２０は、選択した測定地点の高さ推移データに基づく、対象地点の推定変化推移データを推定する。
　次に、対象地点の推定変化推移データの推定について、詳述する。推定部１２０は、まず、受け取った測定データに含まれる、複数の測定地点における高さ推移データから、それらの複数の測定地点における、高さの変化の推移を算出してよい。
　推定部１２０は、ある測定地点における地表面の高さの変化の推移として、前の時点の地表面の高さが得られている各時点における地表面の高さの、直前の時点における地表面の高さからの変化を算出してよい。以下、測定が行われた時点を、測定時点と表記する。この場合、ある測定時点における高さの変化は、例えば、その測定時点における地表面の高さから、その測定時点の直前の測定時点における地表面の高さを引いた値であってよい。推定部１２０は、ある測定時点における高さの変化として、その測定時点における地表面の高さの、その測定時点の直前の測定時点における地表面の高さからの、所定期間当たりの変化を算出してもよい。所定期間は、適宜定められていてよい。所定期間は、例えば、１か月であってもよい。所定期間は、例えば、１年であってもよい。この場合、推定部１２０は、例えば、その測定時点における地表面の高さから、その測定時点の直前の測定時点における地表面の高さを引いた値に、所定の時間間隔を表す値を掛け、さらに、その測定時点とその測定時点の直前の即手時点から期間を表す値で割ってよい。推定部１２０は、得られた値を、ある測定時点における高さの変化としてよい。以下の説明において、２つの時点において測定された地表目の高さの測定値の差、又は、その差から算出された所定の時間間隔当たりの変化を、高さの変化の測定値とも表記する。
　以下では、複数の測定地点における地表面の高さの変化の推移を表すデータを、変化データと表記する。ある測定地点における、地表面の高さの変化の推移を表すデータを、変化推移データと表記する。変化データは、複数の測定地点の高さ推移データを含む。
　ある測定地点における変化推移データは、その測定地点の、前の測定時点が存在する測定時点の各々における、地表面の高さの変化を表すデータを含む。変化推移データは、測定時点を表すデータと、その測定時点における地表面の高さの変化を表すデータとの、複数の組み合わせであってもよい。変化データは、例えば、複数の測定地点の、地点情報と変化推移データとの組み合わせを含んでいてよい。
　上述のように、変化推移データは、地表面の高さの変化の推移、すなわち、地表面の高さの測定値から算出された、地表面の高さの変化の時系列を表す。推定部１２０は、対象地点における変化推移データに含まれる、ある測定時点における地表面の高さの変化と、その測定時点と同じ測定時点における、選択された測定地点における地表面の高さの変化を特定する。推定部１２０は、例えば、対象地点における、ある測定時点における地表面の高さの変化と、その測定時点と同じ測定時点における地表面の高さの変化として特定された地表面の高さの変化と、の統計値を算出してよい。統計値は、例えば、平均値、中央値、中間値、又は、最頻値等であってよい。推定部１２０は、算出した統計値が、対象地点における地表面の高さの変化の推定値であると推定する。なお、指定領域に含まれる対象地点に、対象地点以外の、対象条件を満たす測定地点が存在しない場合、対象地点における地表面の高さの変化を、対象地点における地表面の高さの変化の推定値としてよい。推定部１２０は、対象地点における地表面の高さの変化の各々の推定値を、同様に推定する。推定部１２０は、対象地点における地表面の高さの変化の推定値の時系列が、対象地点における地表面の高さの変化の推移の推定値であると推定する。
　推定部１２０は、指定領域に含まれる測定地点における地表面の高さの変化の推定値を、その測定地点を対象地点として選択して上述の処理を行うことによって推定してよい。推定部１２０は、指定領域に含まれる測定地点の各々における地表面の高さの変化の推移の推定値を、同様に推定してよい。
　推定部１２０は、指定範囲内の測定地点の各々の、地表面の高さの変化の推定値の推移のデータである、推定変化推移データを生成する。推定変化推移データは、高さの変化の推定値が得られた各測定時点の、測定時点を表すデータと、地表面の高さの変化の推定値との組み合わせを含むデータであってよい。
　推定部１２０は、指定範囲内の測定地点の各々の、測定地点の位置を表す地点情報と地表面の高さの変化の推定値の推移のデータである推定変化推移データと、を含む、推定データを生成する。
　推定部１２０は、指定領域に含まれる測定地点の各々の、測定地点の位置を表す地点情報と地表面の高さの変化の推定値の推移を表す推定変化推移データとの組み合わせ（すなわち、上述の推定データ）を、出力部１４０に送出する。
　＜＜決定部１１０＞＞
　決定部１１０は、推定部１２０から、対象地点の地点情報を受け取る。決定部１１０は、推定部１２０から、さらに、対象地点の地理空間情報の値を受け取ってもよい。決定部１１０は、地理空間情報取得部１６０を介して、地理空間情報記憶装置３００から対象地点の地理空間情報の値を受け取ってもよい。以下では、推定部１２０は対象地点の地理空間情報の値を決定部１１０に送出するよう構成され、決定部１１０は対象地点の地理空間情報の値を推定部１２０から受け取るよう構成されている場合について説明する。決定部１１０は、対象地点の地点情報に基づく選択条件を決定してもよい。決定部１１０は、対象地点の地理空間情報の値に基づく、選択条件を決定してもよい。決定部１１０は、対象地点の地点情報と、対象地点の地理空間情報の値とに基づく、選択条件を決定してもよい。決定部１１０は、決定した選択条件、言い換えると、決定した選択条件を表す情報を、推定部１２０に送出する。
　＜＜選択条件の詳細＞＞
　対象地点の地点情報に基づく選択条件（以下、第１選択条件と表記）は、例えば、対象地点の位置からの距離が所定距離以下であることであってよい。
　対象地点の地点情報と対象地点の地理空間情報とに基づく選択条件は、第１選択条件と、対象地点の地理空間情報とに基づく選択条件（以下、第２選択条件と表記）とを満たすことであってもよい。
　第２選択条件は、例えば、所定の種類の地理空間情報の値が、対象地点のその所定の種類の地理空間情報の値と同じであることであってもよい。
　第２選択条件は、例えば、所定の種類の地理空間情報の値が、対象地点のその所定の種類の地理空間情報の値と、所定の基準と比較して近いことであってもよい。この場合、所定の種類の地理空間情報の値に距離が定義できるのであれば、第２選択条件は、所定の種類の地理空間情報の値から、対象地点のその所定の種類の地理空間情報の値までの距離が、所定値よりも小さいことであってもよい。地理空間情報の値は、複数のカテゴリに分けられていてもよい。その場合、第２選択条件は、例えば、その地理空間情報の値が、対象地点のその地理空間情報の値が含まれるカテゴリに含まれることであってもよい。第２選択条件は、複数の種類の地理空間情報の値に関する条件であってもよい。
　具体的には、第２選択条件は、地理空間情報のうち、例えば、都市域土地利用の値が、対象地点の都市域土地利用の値と同じであることであってもよい。第２選択条件は、例えば、人工地形分類の値が、対象地点の人工地形分類の値と同じであることであってもよい。第２選択条件は、例えば、表層地質の値が、対象地点の表層地質の値と同じであることであってもよい。第２選択条件は、例えば、都市域土地利用の値、人工地形分類の値、及び、表層地質の値の少なくとも１つが、対象地点と共通であることであってもよい。
　第２選択条件は、例えば、地理空間情報のうち建造物情報の値が、対象地点の建造物情報の値と同じであることであってもよい。この場合、選択条件は、第２選択条件のみを満たすことであってよい。
　＜＜地理空間情報取得部１６０＞＞
　地理空間情報取得部１６０は、地理空間情報記憶装置３００から、地理空間情報の値を取得する。地理空間情報取得部１６０は、例えば、推定部１２０から、測定地点の地点情報が指定された（すなわち、測定地点の地点情報を含む）、地理空間情報の値の要求を受け取る。地理空間情報取得部１６０は、複数の測定地点の、地点情報を含む地理空間情報の値の複数の要求を、推定部１２０から受け取ってもよい。地理空間情報取得部１６０は、受け取った要求に指定された地点情報が表す地点における、所定の種類の地理空間情報を、地理空間情報記憶装置３００から読み出す。具体的には、地理空間情報取得部１６０は、受け取った要求によって要求された、地理空間情報の値の要求を、地理空間情報記憶装置３００に送信する。そして、地理空間情報取得部１６０は、受け取った要求に指定された地点情報が表す地点における、所定の種類の地理空間情報を、地理空間情報記憶装置３００から受け取る。所定の種類は、地理空間情報の全ての種類のうち、あらかじめ定められている一部の種類であってもよい。所定の種類は、地理空間情報の全ての種類であってもよい。
　地理空間情報取得部１６０は、決定部１１０から、同様に、対象地点の地理空間情報の値の要求（例えば、対象地点の地点情報を含む、地理空間情報の値の要求）を受け取ってよい。地理空間情報取得部１６０は、同様に、対象地点の地点情報が表す地点における地理空間情報の値を、地理空間情報記憶装置３００から読み出してよい。地理空間情報取得部１６０は、地理空間情報記憶装置３００から読み出した、すなわち、地理空間情報記憶装置３００から受け取った、対象地点の地点情報が表す地点における地理空間情報の値を、決定部１１０に送出する。
　＜＜出力部１４０＞＞
　出力部１４０は、推定部１２０から、指定領域内の測定地点における地表面の高さの変化の推定値を受け取る。具体的には、出力部１４０は、推定部１２０から、指定領域内の測定地点における地表面の高さの変化の推定値を含む、上述の推定データを受け取る。出力部１４０は、受け取った、地表面の高さの変化の推定値を、出力先装置４００に出力する。出力部１４０は、所定の測定時点における、地表面の高さの変化の推定値を、出力先装置４００に出力してもよい。出力部１４０は、例えば、推定装置１００のユーザによって指示された測定時点における、地表面の高さの変化の推定値を、出力先装置４００に出力してもよい。出力部１４０は、推定データに含まれる、最も古い測定時点から最も新しい測定時点までの、各測定時点における地表面の高さの変化の推定値を、順に、出力先装置４００に出力してもよい。出力部１４０は、推定データを出力先装置４００に出力してもよい。後述のように出力先装置４００がコンピュータである場合、出力部１４０は、出力先装置４００から指示された測定時点における地表面の高さの変化の推定値を、出力先装置４００に出力してもよい。
　＜＜出力先装置４００＞＞
　出力先装置４００は、例えば、ディスプレイであってもよい。その場合、出力部１４０は、地表面の高さの変化の推定値を表すグラフなどの図形や文字情報を、出力先装置４００に表示してよい。出力先装置４００は、例えば、記憶装置であってもよい。その場合、出力部１４０は、地表面の高さの変化の推定値を表す情報を、出力先装置４００に格納してよい。出力先装置４００は、端末装置又はサーバ装置などのコンピュータであってもよい。その場合、出力部１４０は、地表面の高さの変化の推定値を表す情報を、出力先装置４００に送信してよい。
　＜＜動作＞＞
　次に、本実施形態の動作について図面を参照して詳細に説明する。
　図２は、本実施形態の推定装置１００の動作の一例を表すフローチャートである。
　図２に示す例では、まず、推定部１２０が、測定データ取得部１５０が測定データ記憶装置２００から取得した、測定データを受け取る（ステップＳ１０１）。推定部１２０は、測定データの指定範囲に含まれる測定地点における地理空間情報の値を、地理空間情報取得部１６０を介して、地理空間情報記憶装置３００から取得する（ステップＳ１０２）。推定部１２０は、測定データの指定範囲に含まれる測定地点から、対象地点を選択する（ステップＳ１０３）。推定部１２０は、まだ対象地点として選択されていない測定地点の１つを、対象地点として選択してよい。
　次に、決定部１１０が、上述のように、選択条件を決定する（ステップＳ１０４）。
　推定部１２０は、選択条件を満たす測定地点を選択する（ステップＳ１０５）。推定部１２０は、選択した測定地点（すなわち、選択条件を満たす測定地点）における測定データから、対象地点における変化推移データの推定値（すなわち、上述の推定変化推移データ）を推定する（ステップＳ１０６）。なお、この動作の説明では、対象地点も、選択条件を満たす測定地点に含まれる。従って、言い換えると、推定部１２０は、対象地点の測定データと、選択条件を満たす測定地点の測定データとに基づく、対象地点の変化推移データの推定値を推定する。
　対象範囲に含まれる測定地点に対象地点として選択されていない測定地点が存在する場合（ステップＳ１０７においてＮＯ）、推定装置１００の動作は、ステップＳ１０３に戻る。対象範囲に含まれるすべての測定地点が対象地点として選択された場合（ステップＳ１０６においてＹＥＳ）、出力部１４０は、変化推移データの推定値（すなわち、上述の推定データ）を、出力先装置４００に出力する（ステップＳ１０８）。
　＜＜効果＞＞
　測定値の変化（すなわち、地表面の高さの変化の測定値）に対象に応じたばらつきがある場合であっても、変化（すなわち、地表面の高さの変化の推定値）の推定の精度を向上させることができるという効果がある。その理由は、推定部１２０が、測定地点のうち対象地点の地表面の高さの推移の測定値と、選択条件を満たす測定地点の地表面の高さの推移の測定値とに基づいて、対象地点の地表面の高さの変化の推定値を推定するからである。そして、決定部１１０が、対象地点の地理空間情報（例えば、所定の種類の地理空間情報の値）に基づいて、選択条件を決定するからである。
　＜＜第１の実施形態の第１の変形例＞＞
　次に、第１の実施形態の第１の変形例について説明する。以下で説明する第１から第３の変形例に係る推定システム１は、第１の実施形態の推定システム１の構成と同じ構成を備える。以下では、第１の実施形態との相違点を説明する。
　本変形例では、推定部１２０は、測定地点と対象地点における、地理空間情報の類似性に基づく重みに基づいて、対象地点における地表面の高さの変化の推定値を推定してもよい。例えば、推定部１２０は、対象地点と測定地点との間の、地理空間情報の値の距離に応じた重みを決定してもよい。推定部１２０は、地理空間情報の値の距離が短いほど重みの値が小さくなるように、重みを決定してもよい。
　例えば、地理空間情報の値が複数のカテゴリに分類される場合、推定部１２０は、測定地点の地理空間情報の値が、対象地点の地理空間情報の値と同じである場合、重みの値を１に決定してよい。測定地点の地理空間情報の値が対象地点の地理空間情報の値と異なり、測定地点の地理空間情報の値のカテゴリが対象地点の地理空間情報の値のカテゴリと同じである場合、推定部１２０は、重みの値を０よりも大きく１よりも小さい値に決定してもよい。推定部１２０は、選択条件を満たさない測定地点の重みを決定しなくてよい。
　推定部１２０は、測定地点と対象地点との間の距離に応じた大きさの、測定地点についての重みを決定してもよい。推定部１２０は、例えば、距離の大小関係と重みの大小関係が逆転しないように、距離を重みに変換する関数によって、測定地点と対象地点との間の距離から、測定地点についての重みを決定してもよい。
　推定部１２０は、測定地点と対象地点との間の地理空間情報の類似性と、測定地点と対象地点との間の距離とに基づく重みを決定してもよい。この場合、推定部１２０は、測定地点と対象地点との間の地理空間情報の類似性に基づく重みと、測定地点と対象地点との間の距離とに基づく重みと、の積を、測定地点についての重みに決定してもよい。
　例えば、推定部１２０は、選択条件を満たす測定地点の地表面の高さの変化の測定値とその測定地点について決定した重みとの積を算出してよい。推定部１２０は、選択条件を満たす全ての測定地点について算出した上述の積の総和を算出してよい。推定部１２０は、選択情景を満たす測定地点の重みの総和を算出してよい。そして、推定部１２０は、上述の積の総和を上述の重みの総和で割った値を、対象地点における地表面の高さの変化の推定値としてもよい。
　＜＜第１の実施形態の第２の変形例＞＞
　次に、第１の実施形態の第２の変形例について説明する。
　一般的に、人工衛星等の飛翔体に搭載したレーダーによって観測を行う場合、対象地点が水面であれば、アンテナから放射された電磁波は、波立った水面で鏡面反射し、たまたま測定時に法線方向と電磁波の方向が一致した水面上の点で反射した電磁波がアンテナの方向に戻る。水面の形状は、波の大きさや形によって様々に変化するので、移動するアンテナが、複数の時点において、同一の高さの同一の地点からの電磁波を受信することは無いとみなせる。したがって、測定地点が水面である場合、合成開口レーダーによって測定地点における水面の高さの推移の正確な測定値を得ることは困難である。そのような水面の高さの推移の測定値は、ノイズとみなせる。さらに、本開示では、建造物や地面の高さの推移は測定の対象であるが、水面は測定対象ではない。対象地点の地表面が水面である場合（例えば、対象地点の、地理空間情報の一種である自然地形分類の値が水面を表す場合）、推定部１２０は、対象地点の地表面の高さの変化の推定値を推定しなくてもよい。この場合、推定部１２０は、対象地点の地表面の高さの変化の推定値を、全て０に設定してもよい。また、対象地点の地表面の高さの推定値の算出に、地表面が水面である測定地点の地表面の高さの測定値を使用しなくてよい。
　＜＜第１の実施形態の第３の変形例＞＞
　次に、第１の実施形態の第３の変形例について説明する。
　地理空間情報は、土地の変更履歴を含んでいてもよい。土地の変更履歴は、過去の複数の時点における、上述の地理空間情報のうち少なくともいずれかの地理空間情報の値によって表されていてよい。推定部１２０は、選択条件を満たす測定地点から、土地の変更履歴に基づいて測定地点を選択してもよい。具体的には、推定部１２０は、選択条件を満たす測定地点から、さらに、土地の変更履歴が対象地点の土地の変更履歴と同じ測定地点を選択してもよい。推定部１２０は、このように選択された測定地点の地表面の高さの推移の測定値に基づいて、対象地点の地表面の高さの変化の推移の推定値を推定してもよい。
　＜＜第１の実施形態の第４の変形例＞＞
　推定装置１００は、測定データ記憶装置２００を含んでいてもよい。推定装置１００は、地理空間情報記憶装置３００を含んでいてもよい。推定装置１００は、測定データ記憶装置２００と地理空間情報記憶装置３００とを含んでいてもよい。
　＜＜第１の実施形態の第５の変形例＞＞
　測定データ記憶装置２００に格納されている測定データは、高さ推移データに基づいてあらかじめ算出された変化推移データを更に含んでいてもよい。その場合、推定部１２０は、変化推移データを算出しなくてもよい。
　＜＜第１の実施形態の第６の変形例＞＞
　第１の実施形態に、上述の第１から第５の変形例のいずれか２つ以上を組み合わせた変形を行うこともできる。
　＜第２の実施形態＞
　次に、本開示の第２の実施形態について図面を参照して詳細に説明する。
　＜＜構成＞＞
　図３は、本実施形態に係る推定システム２の構成の例を表すブロック図である。
　図３に示す推定システム２は、以下で説明する相違を除いて、第１の実施形態の推定システム１と同じである。本実施形態の測定データ記憶装置２００は、第１の実施形態の測定データ記憶装置２００と同じである。本実施形態の地理空間情報記憶装置３００は、第１の実施形態の地理空間情報記憶装置３００と同じである。第１の実施形態の推定システム１と比較すると、推定システム２は、推定装置１００の代わりに推定装置１０１を含む。推定装置１０１は、第１の実施形態の推定装置１００の各要素に加えて、判定部１３０を含む。本実施形態の推定装置１０１の判定部１３０以外の要素は、以下で説明する相違点を除いて、第１の実施形態の推定装置１００の、同じ名称及び符号が付与されている要素と同じである。
　＜＜決定部１１０＞＞
　本実施形態の決定部１１０によって決定される選択条件は、地理空間情報の一種である建造物情報の値が、対象地点の建造物情報の値と同じであることであってよい。本実施形態の決定部１１０は、上述の選択条件を決定してよい。
　＜＜推定部１２０＞＞
　本実施形態の推定部１２０は、第１の実施形態の推定部１２０と同様に動作する。本実施形態の推定部１２０は、さらに、上述の推定データを、判定部１３０に送出する。本実施形態の推定部１２０は、更に、指定範囲内の測定地点における、地理空間情報の一種である建造物情報の値を、判定部１３０に送出してもよい。
　＜＜判定部１３０＞＞
　判定部１３０は、推定データを推定部１２０から受け取る。判定部１３０は、さらに、指定範囲内の測定地点における建造物情報の値を、推定部１２０から受け取ってもよい。判定部１３０は、推定データに含まれる地点情報が示す地点における、地理空間情報の一種である建造物情報の値を、地理空間情報取得部１６０を介して、地理空間情報記憶装置３００から受け取ってもよい。
　一般に、例えば航空写真によって地表面が建造物であるか否かを判定することはできるが、建物の種類を航空写真のみによって判定することは必ずしも容易ではない。建造物情報として範囲が得られている個々の建物の種類を得たい場合、例えば調査員が現地に赴いて建造物の種類を判定することは可能ではあるが、時間やコストが必要である。また、建造物の外見から建造物の種類を判断できない場合もある。判定部１３０は、測定地点が建造物である場合、地表面の高さの変化の推移の推定値に基づいて、測定地点を含む建造物の種類を推定してもよい。例えば、建造物が鉄筋コンクリートである場合、建造物が木造建築である場合と比較して、夏季と冬季とで温度の差による高さの変動（高さの季節変動と表記）が大きい。建造物が、鉄骨の骨組みとパネルの壁とによって建造された建造物である場合、木造建築や鉄筋コンクリートの建造物と比較して、更に高さの季節変動が大きい傾向がある。判定部１３０は、例えば、あらかじめ機械学習によって学習した、建造物の種類による高さの季節変動の大きさを判定する判定器によって、建造物の種類を判定してよい。この場合の機械学習の方法として、既存の機械学習の方法のいずれかが使用されてよい。なお、高さの変動は、必ずしも高さの季節変動に限られない。判定部１３０は、指定領域に含まれる、建造物である測定地点の各々の、建造物の種類を判定してよい。判定部１３０は、建造物ではない測定地点の建造物の種類として、建造物ではないことを表す情報を設定してもよい。
　判定部１３０は、指定領域の建造物地図を、地理空間情報取得部１６０を介して、地理空間情報記憶装置３００から取得してもよい。この場合、地理空間情報の１種類として、建造物地図が地理空間情報記憶装置３００にあらかじめ格納されている。
　判定部１３０は、判定した、指定領域に含まれる測定地点の建造物の種類を表す情報を、出力部１４０に送出してよい。判定部１３０は、さらに、取得した建造物地図を表す情報を、出力部１４０に送出してもよい。
　＜＜出力部１４０＞＞
　本実施形態の出力部１４０は、第１の実施形態の出力部１４０と同じ機能を備える。本実施形態の出力部１４０は、さらに、判定部１３０から、指定領域に含まれる測定地点の建造物の種類を表す情報を受け取る。本実施形態の出力部１４０は、さらに、判定部１３０から、建造物地図を表す情報を受け取ってもよい。出力部１４０は、受け取った、指定領域に含まれる測定地点の建造物の種類を表す情報を、出力先装置４００に出力する。
　出力部１４０は、測定地点の位置に対応する建造物地図の位置に、測定地点の建造物の種類を表すマークを重畳してもよい。出力部１４０は、建造物の種類を表すマークが重畳された建造物地図を、出力先装置４００に出力してもよい。
　出力部１４０は、建造物地図に含まれる建造物の種類を、建造物の領域に含まれる測定地点の、推定された建造物の種類に基づいて決定してもよい。出力部１４０は、建造物地図に含まれる建造物の領域の模様及び色の少なくとも何れかを、決定された建造物の種類に応じて変更してもよい。この場合、建造物の領域の模様及び色の少なくともいずれかは、建造物の種類ごとにあらかじめ定められていてよい。建造物の種類の決定方法は、あらかじめ適宜定められた方法であってよい。出力部１４０は、例えば、建造物に１つの測定地点のみが含まれる場合、その建造物の種類を、その測定地点について判定された建造物の種類に決定してよい。出力部１４０は、例えば、建造物に複数の測定地点が含まれる場合、その建造物の種類を、最も多くの測定地点について判定された建造物の種類に決定してもよい。建造物の領域に含まれる測定地点が存在しない場合、及び、最も多くの測定地点について判定された建造物の種類が複数存在する場合など、建造物の種類を決定できない場合、出力部１４０は、その建造物の領域の模様及び色を変更しなくてよい。
　＜＜動作＞＞
　次に、本実施形態の推定装置１０１の動作について、図面を参照して詳細に説明する。
　図４は、本実施形態の推定装置１０１の動作の例を表すフローチャートである。本実施形態の推定装置１０１の、ステップＳ１０１からステップＳ１０７までの動作は、第１の実施形態の推定装置１００の、ステップＳ１０１からステップＳ１０７までの動作と同じである。本実施形態の推定装置１０１は、ステップＳ１０７の動作の後、以下で説明する動作を行う。
　判定部１３０は、測定地点が建造物の領域に含まれる場合、測定地点の、上述の推定データに基づいて、測定地点を領域に含む建造物の種類を判定する（ステップＳ２０８）。判定部１３０は、いずれかの建造物の領域に含まれる測定地点の各々について、建造物の種類を判定してよい。
　出力部１４０は、測定データの推定値（すなわち、上述の推定データ）と、建造物の種類とを出力する（ステップＳ２０９）。
　＜＜効果＞＞
　本実施形態には、第１の実施形態と同じ効果がある。その理由は、第１の実施形態の効果が生じる理由と同じである。
　本実施形態には、さらに、建造物の種類を知ることができるという効果がある。その理由は、判定部１３０が、建造物の領域に含まれる測定地点の推定データに基づいて、その建造物の種類を判定するからである。
　＜＜第２の実施形態の第１の変形例＞＞
　第２の実施形態を、第１の実施形態の第１から第５の変形例のように変形することもできる。
　＜＜第２の実施形態の第２の変形例＞＞
　本変形例の出力部１４０は、受け取った、指定領域に含まれる測定地点の建造物の種類を表す情報を、出力先装置４００に出力し、推定データを出力しなくてもよい。本変形例の推定部１２０は、推定データを出力部１４０に送出しなくてよい。本変形例の出力部１４０は、推定データを推定部１２０から受け取らなくてよい。
　＜第３の実施形態＞
　次に、本開示の第３の実施形態について図面を参照して詳細に説明する。
　＜＜構成＞＞
　図５は、本実施形態に係る推定装置１０２の構成の例を表すブロック図である。
　図５に示す例では、本開示の一態様に係る推定装置１０２は、決定部１１０と、推定部１２０と、を備える。決定部１１０は、合成開口レーダーによって地表面の高さの推移が測定された複数の測定地点のうち対象地点の地理空間情報に基づいて、選択条件を決定する。決定部１１０は、第１の実施形態の決定部１１０と同様に、選択条件を決定してよい。推定部１２０は、前記対象地点における地表面の高さの推移の測定値と、前記選択条件を満たす前記測定地点における地表面の高さの推移の測定値とに基づいて、前記対象地点における地表面の高さの変化の推定値を推定する。推定部１２０は、第１の実施形態の推定部１２０と同様に、対象地点における地表面の高さの変化の推定値を推定してよい。
　＜＜動作＞＞
　図６は、本実施形態に係る推定装置１０２の動作の例を表すフローチャートである。
　図６に示す例では、まず、決定部１１０が、対象地点の地理空間情報に基づいて、選択条件を決定する（ステップＳ３０１）。次に、推定部１２０が、対象地点における地表面の高さの推移の測定値と、選択条件を満たす測定地点における地表面の高さの推移の測定値とに基づいて、対象地点の地表面の高さの変化の推定値を推定する（ステップＳ３０２）。
　＜＜効果＞＞
　本実施形態には、第１の実施形態と同じ効果がある。その理由は、第１の実施形態の効果が生じる理由と同じである。
　＜＜他の実施形態＞＞
　上述の推定装置１００、推定装置１０１、推定装置１０２、測定データ記憶装置２００、及び、地理空間情報記憶装置３００は、それぞれ、コンピュータによって実現することができる。そのコンピュータは、記憶媒体から読み出されたプログラムがロードされたメモリと、そのプログラムを実行するプロセッサとを含む。コンピュータは、互いに通信可能に接続された複数のコンピュータの組み合わせによって実現されていてもよい。推定装置１００、推定装置１０１、推定装置１０２、測定データ記憶装置２００、及び、地理空間情報記憶装置３００は、それぞれ、専用のハードウェアによって実現することもできる。専用のハードウェアは、例えば、１つの回路、又は、互いに通信可能に接続された複数の回路によって実現されていてよい。推定装置１００、推定装置１０１、推定装置１０２、測定データ記憶装置２００、及び、地理空間情報記憶装置３００は、前述のコンピュータと専用のハードウェアとの組み合わせによって実現することもできる。
　図７は、本開示の実施形態に係る推定装置などの装置を実現することができる、コンピュータ１０００のハードウェア構成の一例を表す図である。図７を参照すると、コンピュータ１０００は、プロセッサ１００１と、メモリ１００２と、記憶装置１００３と、Ｉ／Ｏ（Ｉｎｐｕｔ／Ｏｕｔｐｕｔ）インタフェース１００４とを含む。また、コンピュータ１０００は、記憶媒体１００５にアクセスすることができる。メモリ１００２と記憶装置１００３は、例えば、ＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）、ハードディスクなどの記憶装置である。記憶媒体１００５は、例えば、ＲＡＭ、ハードディスクなどの記憶装置、ＲＯＭ（Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）、可搬記憶媒体である。記憶装置１００３が記憶媒体１００５であってもよい。プロセッサ１００１は、メモリ１００２と、記憶装置１００３に対して、データやプログラムの読み出しと書き込みを行うことができる。プロセッサ１００１は、Ｉ／Ｏインタフェース１００４を介して、例えば、他の装置にアクセスすることができる。プロセッサ１００１は、記憶媒体１００５にアクセスすることができる。記憶媒体１００５には、コンピュータ１０００を、上述の実施形態に係る装置として動作させるプログラムが格納されている。
　プロセッサ１００１は、記憶媒体１００５に格納されている、コンピュータ１０００を、上述の実施形態に係る推定装置として動作させるプログラムを、メモリ１００２にロードする。そして、プロセッサ１００１が、メモリ１００２にロードされたプログラムを実行することにより、コンピュータ１０００は、上述の実施形態に係る推定装置として動作する。
　プロセッサ１００１は、記憶媒体１００５に格納されている、コンピュータ１０００を、上述の実施形態に係る測定データ記憶装置として動作させるプログラムを、メモリ１００２にロードする。そして、プロセッサ１００１が、メモリ１００２にロードされたプログラムを実行することにより、コンピュータ１０００は、上述の実施形態に係る測定データ記憶装置として動作する。
　プロセッサ１００１は、記憶媒体１００５に格納されている、コンピュータ１０００を、上述の実施形態に係る地理空間情報記憶装置として動作させるプログラムを、メモリ１００２にロードする。そして、プロセッサ１００１が、メモリ１００２にロードされたプログラムを実行することにより、コンピュータ１０００は、上述の実施形態に係る地理空間情報記憶装置として動作する。
　決定部１１０、推定部１２０、判定部１３０、出力部１４０、測定データ取得部１５０、及び、地理空間情報取得部１６０は、例えば、メモリ１００２にロードされたプログラムを実行するプロセッサ１００１により実現することができる。第１入出力部２２０、及び、第２入出力部３２０は、例えば、メモリ１００２にロードされたプログラムを実行するプロセッサ１００１により実現することができる。また、測定データ記憶部２１０、及び、地理空間情報記憶部３１０は、コンピュータ１０００が含むメモリ１００２やハードディスク装置等の記憶装置１００３により実現することができる。決定部１１０、推定部１２０、判定部１３０、出力部１４０、測定データ取得部１５０、地理空間情報取得部１６０、第１入出力部２２０、及び、第２入出力部３２０の一部又は全部を、各部の機能を実現する専用の回路によって実現することもできる。測定データ記憶部２１０、及び、地理空間情報記憶部３１０の一部又は全部を、各部の機能を実現する専用の回路によって実現することもできる。
　また、上記の実施形態の一部又は全部は、以下の付記のようにも記載されうるが、以下には限られない。
　（付記１）
　合成開口レーダーによって地表面の高さの推移が測定された複数の測定地点のうち対象地点の地理空間情報に基づいて、選択条件を決定する決定手段と、
　前記対象地点における地表面の高さの推移の測定値と、前記選択条件を満たす前記測定地点における地表面の高さの推移の測定値とに基づいて、前記対象地点における地表面の高さの変化の推定値を推定する推定手段と、
　を備える推定装置。
　（付記２）
　前記決定手段は、前記対象地点の地理空間情報に基づく地表面の種類と、前記対象地点からの距離と、に基づく前記選択条件を決定する
　付記１に記載の推定装置。
　（付記３）
　前記推定手段は、前記測定地点と前記対象地点とにおける、地表面の種類の類似性に基づく重みに基づいて、前記推定値を推定する
　付記１又は２のいずれか１項に記載の推定装置。
　（付記４）
　前記推定手段は、前記測定地点と前記対象地点との距離にさらに基づく前記重みに基づいて、前記推定値を推定する
　付記３に記載の推定装置。
　（付記５）
　前記推定手段は、前記選択条件を満たす前記測定地点から、当該測定地点における土地の変更履歴に基づいて前記測定地点を更に選択し、当該測定地点の前記測定値に基づいて前記推定値を推定する
　付記１乃至４のいずれか１項に記載の推定装置。
　（付記６）
　前記推定手段は、前記対象地点の地表面が水面である場合、当該対象地点における前記推定値を推定せず、前記測定地点の地表面が水面である場合、当該測定地点の高さの変動の測定値を、前記対象地点における推定値の推定に使用しない
　付記１乃至５のいずれか１項に記載の推定装置。
　（付記７）
　前記決定手段は、前記地理空間情報において前記対象地点が建造物の領域に含まれる場合、前記建造物の領域に含まれる前記測定地点を選択するように前記選択条件を決定する
　付記１乃至６のいずれか１項に記載の推定装置。
　（付記８）
　前記地理空間情報において前記対象地点が建造物の領域に含まれる場合、前記推定値の特徴に基づいて、前記対象地点を含む建造物の種類を判定する判定手段
　を更に備える付記１乃至７のいずれか１項に記載の推定装置。
　（付記９）
　合成開口レーダーによって地表面の高さの推移が測定された複数の測定地点のうち対象地点の地理空間情報に基づいて、選択条件を決定し、
　前記対象地点における地表面の高さの推移の測定値と、前記選択条件を満たす前記測定地点における地表面の高さの推移の測定値とに基づいて、前記対象地点における地表面の高さの変化の推定値を推定する、
　推定方法。
　（付記１０）
　前記対象地点の地理空間情報に基づく地表面の種類と、前記対象地点からの距離と、に基づく前記選択条件を決定する
　付記９に記載の推定方法。
　（付記１１）
　前記測定地点と前記対象地点とにおける、地表面の種類の類似性に基づく重みに基づいて、前記推定値を推定する
　付記９又は１０のいずれか１項に記載の推定方法。
　（付記１２）
　前記測定地点と前記対象地点との距離にさらに基づく前記重みに基づいて、前記推定値を推定する
　付記１１に記載の推定方法。
　（付記１３）
　前記選択条件を満たす前記測定地点から、当該測定地点における土地の変更履歴に基づいて前記測定地点を更に選択し、当該測定地点の前記測定値に基づいて前記推定値を推定する
　付記９乃至１２のいずれか１項に記載の推定方法。
　（付記１４）
　前記対象地点の地表面が水面である場合、当該対象地点における前記推定値を推定せず、前記測定地点の地表面が水面である場合、当該測定地点の高さの変動の測定値を、前記対象地点における推定値の推定に使用しない
　付記９乃至１３のいずれか１項に記載の推定方法。
　（付記１５）
　前記地理空間情報において前記対象地点が建造物の領域に含まれる場合、前記建造物の領域に含まれる前記測定地点を選択するように前記選択条件を決定する
　付記９乃至１４のいずれか１項に記載の推定方法。
　（付記１６）
　前記地理空間情報において前記対象地点が建造物の領域に含まれる場合、前記推定値の特徴に基づいて、前記対象地点を含む建造物の種類を判定する
　付記９乃至１５のいずれか１項に記載の推定方法。
　（付記１７）
　合成開口レーダーによって地表面の高さの推移が測定された複数の測定地点のうち対象地点の地理空間情報に基づいて、選択条件を決定する決定処理と、
　前記対象地点における地表面の高さの推移の測定値と、前記選択条件を満たす前記測定地点における地表面の高さの推移の測定値とに基づいて、前記対象地点における地表面の高さの変化の推定値を推定する推定処理と、
　をコンピュータに実行させるプログラムを記憶する記憶媒体。
　（付記１８）
　前記決定処理は、前記対象地点の地理空間情報に基づく地表面の種類と、前記対象地点からの距離と、に基づく前記選択条件を決定する
　付記１７に記載の記憶媒体。
　（付記１９）
　前記推定処理は、前記測定地点と前記対象地点とにおける、地表面の種類の類似性に基づく重みに基づいて、前記推定値を推定する
　付記１７又は１８のいずれか１項に記載の記憶媒体。
　（付記２０）
　前記推定処理は、前記測定地点と前記対象地点との距離にさらに基づく前記重みに基づいて、前記推定値を推定する
　付記１９に記載の記憶媒体。
　（付記２１）
　前記推定処理は、前記選択条件を満たす前記測定地点から、当該測定地点における土地の変更履歴に基づいて前記測定地点を更に選択し、当該測定地点の前記測定値に基づいて前記推定値を推定する
　付記１７乃至２０のいずれか１項に記載の記憶媒体。
　（付記２２）
　前記推定処理は、前記対象地点の地表面が水面である場合、当該対象地点における前記推定値を推定せず、前記測定地点の地表面が水面である場合、当該測定地点の高さの変動の測定値を、前記対象地点における推定値の推定に使用しない
　付記１７乃至２１のいずれか１項に記載の記憶媒体。
　（付記２３）
　前記決定処理は、前記地理空間情報において前記対象地点が建造物の領域に含まれる場合、前記建造物の領域に含まれる前記測定地点を選択するように前記選択条件を決定する
　付記１７乃至２２のいずれか１項に記載の記憶媒体。
　（付記２４）
　前記プログラムは、
　前記地理空間情報において前記対象地点が建造物の領域に含まれる場合、前記推定値の特徴に基づいて、前記対象地点を含む建造物の種類を判定する判定処理
　を更にコンピュータに実行させる付記１７乃至２３のいずれか１項に記載の記憶媒体。
　以上、実施形態を参照して本発明を説明したが、本発明は上記実施形態に限定されるものではない。本発明の構成や詳細には、本発明のスコープ内で当業者が理解し得る様々な変更をすることができる。
　１　　推定システム
　２　　推定システム
　１００　　推定装置
　１０１　　推定装置
　１０２　　推定装置
　１１０　　決定部
　１２０　　推定部
　１３０　　判定部
　１４０　　出力部
　１５０　　測定データ取得部
　１６０　　地理空間情報取得部
　２００　　測定データ記憶装置
　２１０　　測定データ記憶部
　２２０　　第１入出力部
　３００　　地理空間情報記憶装置
　３１０　　地理空間情報記憶部
　３２０　　第２入出力部
　４００　　出力先装置
　１０００　　コンピュータ
　１００１　　プロセッサ
　１００２　　メモリ
　１００３　　記憶装置
　１００４　　Ｉ／Ｏインタフェース
　１００５　　記憶媒体

Claims (24)

  1. 　合成開口レーダーによって地表面の高さの推移が測定された複数の測定地点のうち対象地点の地理空間情報に基づいて、選択条件を決定する決定手段と、
    　前記対象地点における地表面の高さの推移の測定値と、前記選択条件を満たす前記測定地点における地表面の高さの推移の測定値とに基づいて、前記対象地点における地表面の高さの変化の推定値を推定する推定手段と、
    　を備える推定装置。
  2. 　前記決定手段は、前記対象地点の地理空間情報に基づく地表面の種類と、前記対象地点からの距離と、に基づく前記選択条件を決定する
    　請求項１に記載の推定装置。
  3. 　前記推定手段は、前記測定地点と前記対象地点とにおける、地表面の種類の類似性に基づく重みに基づいて、前記推定値を推定する
    　請求項１又は２のいずれか１項に記載の推定装置。
  4. 　前記推定手段は、前記測定地点と前記対象地点との距離にさらに基づく前記重みに基づいて、前記推定値を推定する
    　請求項３に記載の推定装置。
  5. 　前記推定手段は、前記選択条件を満たす前記測定地点から、当該測定地点における土地の変更履歴に基づいて前記測定地点を更に選択し、当該測定地点の前記測定値に基づいて前記推定値を推定する
    　請求項１乃至４のいずれか１項に記載の推定装置。
  6. 　前記推定手段は、前記対象地点の地表面が水面である場合、当該対象地点における前記推定値を推定せず、前記測定地点の地表面が水面である場合、当該測定地点の高さの変動の測定値を、前記対象地点における推定値の推定に使用しない
    　請求項１乃至５のいずれか１項に記載の推定装置。
  7. 　前記決定手段は、前記地理空間情報において前記対象地点が建造物の領域に含まれる場合、前記建造物の領域に含まれる前記測定地点を選択するように前記選択条件を決定する
    　請求項１乃至６のいずれか１項に記載の推定装置。
  8. 　前記地理空間情報において前記対象地点が建造物の領域に含まれる場合、前記推定値の特徴に基づいて、前記対象地点を含む建造物の種類を判定する判定手段
    　を更に備える請求項１乃至７のいずれか１項に記載の推定装置。
  9. 　合成開口レーダーによって地表面の高さの推移が測定された複数の測定地点のうち対象地点の地理空間情報に基づいて、選択条件を決定し、
    　前記対象地点における地表面の高さの推移の測定値と、前記選択条件を満たす前記測定地点における地表面の高さの推移の測定値とに基づいて、前記対象地点における地表面の高さの変化の推定値を推定する、
    　推定方法。
  10. 　前記対象地点の地理空間情報に基づく地表面の種類と、前記対象地点からの距離と、に基づく前記選択条件を決定する
    　請求項９に記載の推定方法。
  11. 　前記測定地点と前記対象地点とにおける、地表面の種類の類似性に基づく重みに基づいて、前記推定値を推定する
    　請求項９又は１０のいずれか１項に記載の推定方法。
  12. 　前記測定地点と前記対象地点との距離にさらに基づく前記重みに基づいて、前記推定値を推定する
    　請求項１１に記載の推定方法。
  13. 　前記選択条件を満たす前記測定地点から、当該測定地点における土地の変更履歴に基づいて前記測定地点を更に選択し、当該測定地点の前記測定値に基づいて前記推定値を推定する
    　請求項９乃至１２のいずれか１項に記載の推定方法。
  14. 　前記対象地点の地表面が水面である場合、当該対象地点における前記推定値を推定せず、前記測定地点の地表面が水面である場合、当該測定地点の高さの変動の測定値を、前記対象地点における推定値の推定に使用しない
    　請求項９乃至１３のいずれか１項に記載の推定方法。
  15. 　前記地理空間情報において前記対象地点が建造物の領域に含まれる場合、前記建造物の領域に含まれる前記測定地点を選択するように前記選択条件を決定する
    　請求項９乃至１４のいずれか１項に記載の推定方法。
  16. 　前記地理空間情報において前記対象地点が建造物の領域に含まれる場合、前記推定値の特徴に基づいて、前記対象地点を含む建造物の種類を判定する
    　請求項９乃至１５のいずれか１項に記載の推定方法。
  17. 　合成開口レーダーによって地表面の高さの推移が測定された複数の測定地点のうち対象地点の地理空間情報に基づいて、選択条件を決定する決定処理と、
    　前記対象地点における地表面の高さの推移の測定値と、前記選択条件を満たす前記測定地点における地表面の高さの推移の測定値とに基づいて、前記対象地点における地表面の高さの変化の推定値を推定する推定処理と、
    　をコンピュータに実行させるプログラムを記憶する記憶媒体。
  18. 　前記決定処理は、前記対象地点の地理空間情報に基づく地表面の種類と、前記対象地点からの距離と、に基づく前記選択条件を決定する
    　請求項１７に記載の記憶媒体。
  19. 　前記推定処理は、前記測定地点と前記対象地点とにおける、地表面の種類の類似性に基づく重みに基づいて、前記推定値を推定する
    　請求項１７又は１８のいずれか１項に記載の記憶媒体。
  20. 　前記推定処理は、前記測定地点と前記対象地点との距離にさらに基づく前記重みに基づいて、前記推定値を推定する
    　請求項１９に記載の記憶媒体。
  21. 　前記推定処理は、前記選択条件を満たす前記測定地点から、当該測定地点における土地の変更履歴に基づいて前記測定地点を更に選択し、当該測定地点の前記測定値に基づいて前記推定値を推定する
    　請求項１７乃至２０のいずれか１項に記載の記憶媒体。
  22. 　前記推定処理は、前記対象地点の地表面が水面である場合、当該対象地点における前記推定値を推定せず、前記測定地点の地表面が水面である場合、当該測定地点の高さの変動の測定値を、前記対象地点における推定値の推定に使用しない
    　請求項１７乃至２１のいずれか１項に記載の記憶媒体。
  23. 　前記決定処理は、前記地理空間情報において前記対象地点が建造物の領域に含まれる場合、前記建造物の領域に含まれる前記測定地点を選択するように前記選択条件を決定する
    　請求項１７乃至２２のいずれか１項に記載の記憶媒体。
  24. 　前記プログラムは、
    　前記地理空間情報において前記対象地点が建造物の領域に含まれる場合、前記推定値の特徴に基づいて、前記対象地点を含む建造物の種類を判定する判定処理
    　を更にコンピュータに実行させる請求項１７乃至２３のいずれか１項に記載の記憶媒体。
PCT/JP2020/034216 2020-09-10 2020-09-10 推定装置、推定方法及び記憶媒体 WO2022054186A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2022548307A JPWO2022054186A5 (ja) 2020-09-10 推定装置、推定方法及びプログラム
PCT/JP2020/034216 WO2022054186A1 (ja) 2020-09-10 2020-09-10 推定装置、推定方法及び記憶媒体

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2020/034216 WO2022054186A1 (ja) 2020-09-10 2020-09-10 推定装置、推定方法及び記憶媒体

Publications (1)

Publication Number Publication Date
WO2022054186A1 true WO2022054186A1 (ja) 2022-03-17

Family

ID=80631814

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/034216 WO2022054186A1 (ja) 2020-09-10 2020-09-10 推定装置、推定方法及び記憶媒体

Country Status (1)

Country Link
WO (1) WO2022054186A1 (ja)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011209780A (ja) * 2010-03-29 2011-10-20 Mitsubishi Space Software Kk 変化領域特定装置および変化領域特定プログラム。
WO2019215819A1 (ja) * 2018-05-08 2019-11-14 日本電気株式会社 合成開口レーダ画像解析システム、合成開口レーダ画像解析方法および合成開口レーダ画像解析プログラム
WO2019220574A1 (ja) * 2018-05-16 2019-11-21 日本電気株式会社 合成開口レーダ信号解析装置、合成開口レーダ信号解析方法および合成開口レーダ信号解析プログラム
WO2020202259A1 (ja) * 2019-03-29 2020-10-08 日本電気株式会社 合成開口レーダの画像処理装置及び画像処理方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011209780A (ja) * 2010-03-29 2011-10-20 Mitsubishi Space Software Kk 変化領域特定装置および変化領域特定プログラム。
WO2019215819A1 (ja) * 2018-05-08 2019-11-14 日本電気株式会社 合成開口レーダ画像解析システム、合成開口レーダ画像解析方法および合成開口レーダ画像解析プログラム
WO2019220574A1 (ja) * 2018-05-16 2019-11-21 日本電気株式会社 合成開口レーダ信号解析装置、合成開口レーダ信号解析方法および合成開口レーダ信号解析プログラム
WO2020202259A1 (ja) * 2019-03-29 2020-10-08 日本電気株式会社 合成開口レーダの画像処理装置及び画像処理方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ANONYMOUS: "Two-dimensional micro-displacement analysis technology utilizing satellite radar that enables inspection of aging of multiple roads and buildings in urban areas", NEC, 20 July 2018 (2018-07-20), XP055913676, Retrieved from the Internet <URL:https://jpn.nec.com/rd/technologies/201802/index.html> [retrieved on 20220420] *
TANAKA DAICHI: "Development of automation technology for microdisplacement analysis of infrastructure structures using satellite synthetic aperture radar images", 3.2 TECHNOLOGY FOR ASSOCIATING MAP INFORMATION WITH REFLECTION POINT CLUSTERS, 14 January 2020 (2020-01-14), pages 78 - 91, XP055913671, Retrieved from the Internet <URL:http://web.archive.org/web/20200114013956if_/http://www.fbi-award.jp:80/sentan/jusyou/2019/06.pdf> [retrieved on 20220420] *
TANAKA TAICHI; IKEFUJI DAISUKE; HOSHUYAMA OSAMU: "Multitemporal Sar And Map Fusion For Extracting Persitent Scatterers On Roads", IGARSS 2019 - 2019 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM, IEEE, 28 July 2019 (2019-07-28), pages 182 - 185, XP033656451, DOI: 10.1109/IGARSS.2019.8899061 *
TANAKA, TAICHI; HOSHUYAMA, OSAMU: "B-2-17 Multi-temporal synthetic aperture radar interferometry for nonlinear displacement analysis using cluster based persistent scatterer extraction and gaussian process regression", PROCEEDINGS OF THE IEICE GENERAL CONFERENCE 2020: COMMUNICATION 1; HIGASHIHIROSHIMA, JAPAN; MARCH 17-20, 2020, 3 March 2020 (2020-03-03) - 20 March 2020 (2020-03-20), pages 175, XP009535298 *

Also Published As

Publication number Publication date
JPWO2022054186A1 (ja) 2022-03-17

Similar Documents

Publication Publication Date Title
Patel et al. Performances evaluation of different open source DEM using Differential Global Positioning System (DGPS)
Brasington et al. Monitoring and modelling morphological change in a braided gravel‐bed river using high resolution GPS‐based survey
Liu et al. LiDAR-derived high quality ground control information and DEM for image orthorectification
Nelson et al. DEM production methods and sources
EP1709396B1 (en) System,computer program and method for 3d object measurement, modeling and mapping from single imagery
Gong et al. ICEsat GLAS data for urban environment monitoring
Pipaud et al. Evaluation of TanDEM-X elevation data for geomorphological mapping and interpretation in high mountain environments—A case study from SE Tibet, China
Oksanen Digital elevation model error in terrain analysis
Gallay Direct acquisition of data: airborne laser scanning
Mason et al. Data utilization in flood inundation modelling
Khalid et al. Open-source digital elevation model (DEMs) evaluation with GPS and LiDAR data
Karan et al. Digital modeling of construction site terrain using remotely sensed data and geographic information systems analyses
Gonçalves et al. Monitoring Local Shoreline Changes by Integrating UASs, Airborne LiDAR, Historical Images and Orthophotos.
Leitão et al. Improving merge methods for grid-based digital elevation models
Miller et al. A robust surface matching technique for coastal geohazard assessment and management
Fraser et al. Towards precise drone-based measurement of elevation change in permafrost terrain experiencing thaw and thermokarst
González et al. Relative height accuracy estimation method for InSAR-based DEMs
JP7294529B2 (ja) 分析装置、分析方法及びプログラム
Mallet et al. Digital terrain models derived from airborne lidar data
Winiwarter et al. Full 4D Change Analysis of Topographic Point Cloud Time Series using Kalman Filtering
Liu et al. Development of building height data in Peru from high-resolution SAR imagery
Winiwarter et al. Full four-dimensional change analysis of topographic point cloud time series using Kalman filtering
WO2022054186A1 (ja) 推定装置、推定方法及び記憶媒体
Okojie et al. Relative canopy height modelling precision from UAV and ALS datasets for forest tree height estimation
Sharma et al. A method for extracting deformation features from terrestrial laser scanner 3d point clouds data in rgipt building

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: 20953252

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2022548307

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: 20953252

Country of ref document: EP

Kind code of ref document: A1