JP4675868B2 - Ground surface change discrimination device and ground surface change discrimination program - Google Patents

Description

  The present invention relates to a ground surface change discriminating apparatus and a ground surface change discriminating program for discriminating changes in the ground surface of a corresponding area such as topography and structures.

  In the conventional surface change detection device, the altitude information of the surface newly obtained by scanning the surface of the area with radar from an aircraft etc., the color spectrum information newly obtained by photographing with the camera, and the surface information obtained in advance From the difference between the altitude information of the ground surface and the color spectrum information of the ground surface, the change of the ground surface of the corresponding area is discriminated. Examples of changes in the surface of the earth include, for example, generation or disappearance of a building, generation or disappearance of a depression including flooding of a river, generation or disappearance of a constant slope including landslide.

  As described above, in the conventional surface change determination device, the surface of the corresponding area is photographed with visible light, infrared rays, radio waves, etc. from an aircraft or an artificial satellite, and the altitude information of the surface and the color spectrum information of the surface are obtained. Further, it is performed to determine a change in the ground surface including a building with the passage of time from a plurality of ground surface images taken or scanned with a time difference.

  Conventionally, an invention for automatically performing the process of discriminating the change in the ground surface without depending on human visual ability has been made. For example, in Patent Document 1, an object such as a building or a natural object exists in an area. Edge image generation means for generating an edge image from the photographed image, and for each edge image generated by the edge image generation means, the existence accuracy for the entire photographed image is calculated, and an edge image whose existence accuracy is a predetermined value or less And an image processing apparatus that generates a preprocessed image from which the shadow component of the object is deleted.

  Further, in Patent Document 2, new and old images obtained by shooting the same area from the sky at a time are prepared, and first and second converted images are obtained from each of the old and new images by a first image conversion method. The first and second converted images are compared to obtain a first binarized image having data estimated as a change region, and the third and fourth conversions are performed from the old and new images by a second image conversion method. An image is obtained, the third and fourth transformed images are compared, a second binarized image having data estimated as a change area is obtained, and the first binarized image is composed of a predetermined number of pixels. Dividing into a plurality of pixel areas, obtaining the number of pixels indicating the area change in each pixel area, comparing the number of pixels with a threshold value to obtain a third binary image, and obtaining the second binary value The divided image is divided into a plurality of pixel areas each having a predetermined number of pixels, and the area change in each pixel area The number of pixels to be shown is obtained, the number of pixels is compared with a threshold value to obtain a fourth binarized image, and the change area between the old and new images is specified based on the third and fourth binarized images. An image data processing method is disclosed.

  In the conventional ground surface change detection method of the above-mentioned patent document, it is necessary to accurately align the images obtained by photographing the old and new areas, and then detect the ground surface change. Is described on the premise that alignment is performed only by parallel movement of the same scale, or manual alignment is performed which is not disclosed.

  On the other hand, at the time of surveying, it is described in Non-Patent Document 1 that an anti-air sign is installed on the ground and used as a reference for alignment. However, in the case of measurements other than normal surveying, a disaster occurred in particular. In the subsequent shooting on the ground, such an anti-air marking for positioning is not set on the ground as the subject.

  In addition, the position and orientation of the camera and laser irradiation device to be photographed can be obtained with a GPS (Global Positioning System) device and posture control device installed on the aircraft, and this can be used to align the image of the old and new areas. Although it can be used, there are many errors, and relying only on the alignment is not sufficient as an alignment for capturing changes in the ground surface.

Japanese Patent Laid-Open No. 10-312466 (paragraph 0006) JP 2001-109872 (paragraph 0006) Ministry of Land, Infrastructure, Transport and Tourism public survey work regulations, supervised by the Ministry of Land, Infrastructure, Transport and Tourism Minister's Secretariat Technology Research Division, Japan Surveyor Association, June 6, 2002, Chapter 3 Aerial Photogrammetry, Section 1, Summary, Article 98, p. 43

  Since the conventional surface change determination device is configured as described above, the accuracy of the alignment of the captured images of the old and new applicable areas is low, and the change in the surface of the applicable area cannot be analyzed accurately at high speed. was there.

  The present invention has been made to solve the above-described problems, and an object of the present invention is to obtain a ground surface change discriminating apparatus and a ground surface change discriminating program capable of analyzing a ground surface change in a corresponding area at high speed and with high accuracy. .

The ground surface change determination device according to the present invention detects a specific distribution of the surface information of the corresponding area newly obtained, calculates the singular area of the corresponding area newly obtained, and obtains the surface information of the corresponding area obtained in advance. Singular surface information detection means for calculating the singular area of the corresponding area obtained in advance by detecting the singular distribution, and the singular area of the corresponding area newly obtained by the singular ground information detection means and the obtained in advance Based on the specific area of the corresponding area, the position correction means for aligning the newly obtained surface information of the corresponding area and the surface information of the corresponding area obtained in advance, and the position correction means The surface information difference detecting means for detecting the difference between the newly obtained surface information of the corresponding area and the surface information of the corresponding area obtained in advance, and the ground surface set in advance for the difference detected by the surface information difference detecting means. Of information And a surface information change evaluation means for evaluating whether or not the change in the surface information actually occurs compared to a parameter for evaluating the conversion. An area dividing means for dividing the medium area and the small area, an area statistics calculating means for calculating area statistics of the ground information for each predetermined medium area and small area divided by the area dividing means, and the area statistics calculating means. Based on the calculated area statistics of the surface information for each of the medium area and the small area, a specific area calculating means for detecting a specific distribution of the surface information of the corresponding area and calculating a specific area of the corresponding area is provided. .

  By this invention, the effect that the change of the ground surface of an applicable area can be analyzed accurately at high speed is acquired.

Embodiment 1 FIG.
Embodiments of the present invention will be described below with reference to the drawings.
1 is a block diagram showing the structure of a ground surface change determination device according to Embodiment 1 of the present invention. This ground surface change discriminating device discriminates the topography of the corresponding area and the ground surface change of the structure. New arrival information input means 101, new arrival information memory 102, peculiar ground information detection means 103, new arrival peculiar information memory 104, position correction Means 105, new arrival position correction information memory 106, ground surface information difference detection means 107, new arrival information difference memory 108, ground surface information change evaluation means 109, change information memory 110, change information output means 111, old information memory 202, singular ground surface information detection Means 203, a conventional singular information memory 204, and a ground surface information change model memory 208 are provided.

  In FIG. 1, new arrival information input means 101 is, for example, a disk device, and newly obtains ground surface information that is altitude information on the surface of the corresponding area or color spectrum information of the ground surface and stores it in the new arrival information memory 102. The peculiar ground surface information detecting means 103 detects a peculiar distribution of the ground surface information of the newly obtained corresponding area stored in the newly arrived information memory 102, calculates a singular area of the newly obtained corresponding area, and newly arrived peculiar information Store in the memory 104.

  The conventional information memory 202 stores the ground surface information which is the altitude information of the surface of the corresponding area or the color spectrum information of the ground surface obtained in advance. The peculiar land information detecting means 203 detects a peculiar distribution of the surface information of the corresponding area obtained in advance stored in the traditional information memory 202, calculates the peculiar area of the corresponding area obtained in advance, and sets the peculiar singular information. Store in the memory 204.

The position correcting means 105 is newly created based on the newly obtained unique area of the relevant area stored in the new arrival unique information memory 104 and the unique area of the relevant area previously obtained stored in the old unique information memory 204. Obtain the vector amount of positional deviation for each middle region in the obtained ground surface information and the ground surface information of the relevant region obtained in advance, and newly obtain the relevant region based on the obtained positional displacement vector amount for each middle region The position of the ground surface information of the corresponding area newly obtained is corrected by shifting the ground surface information for each small area and stored in the newly arrived position correction information memory 106.
In this way, the position correcting unit 105 is configured so that the singular region of the corresponding area newly calculated by the singular ground information detecting unit 103 and the singular region of the corresponding region obtained in advance calculated by the singular ground information detecting unit 203 are obtained. Based on the above, the ground information of the corresponding area newly obtained is aligned with the surface information of the corresponding area obtained in advance.

  The ground surface information difference detection means 107 includes the ground surface information of the newly obtained corresponding region whose position stored in the newly arrived position correction information memory 106 is corrected, and the corresponding region obtained in advance stored in the old information memory 202. The difference between the newly obtained surface information of the corresponding area and the previously obtained surface information of the corresponding area is stored in the new arrival information difference memory 108.

  The surface information change model memory 208 stores parameters for evaluating changes in the surface information of the corresponding area in advance. The surface information change evaluation means 109 stores, in the surface information change model memory 208, the difference between the newly obtained surface information of the corresponding area stored in the new arrival information difference memory 108 and the surface information of the corresponding area obtained in advance. Compared with the preset parameter for evaluating the change in the surface information, it is evaluated whether or not the change has actually occurred in the surface information, and is stored in the change information memory 110.

  The change information output unit 111 outputs a change in the ground information actually generated stored in the change information memory 110.

  FIG. 2 is a block diagram showing a configuration of the peculiar ground surface information detecting means 103. This peculiar ground surface information detection means 103 includes a region division control unit 301, a region division unit 302, a region statistics calculation unit 303, a region statistics memory 304, and a unique region calculation unit 305. Moreover, although not shown in figure, the specific surface information detection means 203 is also provided with the same structure as the structure of the specific surface information detection means 103 shown in FIG.

  In FIG. 2, the region division control unit 301 is preset with a plurality of region division methods for the corresponding region that divide the entire region of the corresponding region into predetermined middle regions and divide the divided middle region into predetermined small regions. Yes. The area dividing unit 302 divides the corresponding area into predetermined medium areas and small areas based on the area dividing method of the corresponding area set in the area dividing control unit 301.

  The area statistics calculation unit 303 calculates the area statistics of the ground information for each predetermined middle area and small area divided by the area dividing unit 302 and stores the area statistics in the area statistics memory 304. Based on the area statistics of the surface information stored in the area statistics memory 304, the singular area calculation unit 305 detects a singular distribution of the surface information of the corresponding area, and calculates the singular area of the corresponding area.

Next, the operation will be described.
The new arrival information input means 101 obtains newly the ground surface information which is the altitude information of the surface of the corresponding area or the color spectrum information of the ground surface and stores it in the new arrival information memory 102. Here, the surface information is scanned with a laser or the like. The altitude data of the ground surface shall be obtained. The surface information newly obtained by scanning with this laser, etc. is not affected by the shadows of clouds and sunlight when compared with images of optical sensors including infrared rays, but it is possible to obtain high altitude information. It is data and has a feature that a difference in color and heat due to the reflection spectrum of the ground surface is not detected.

  Altitude data acquired by a radar or the like is basically altitude data for each discrete point. The altitude data for each point is handled in a format of discrete altitude data (for example, data format: x, y, z, where x is latitude, y is longitude, and z is altitude). The ground surface information stored in the new arrival information memory 102 is transferred to the specific ground surface information detecting means 103 and the position correcting means 105.

In the region division control unit 301 in the unique ground surface information detection unit 103, a plurality of region division methods for the corresponding region are set in advance.
FIG. 3 is a diagram showing an example of the region dividing method for the corresponding area preset in the region dividing control unit 301 in the specific ground surface information detecting unit 103. In the area dividing method A shown in FIG. 3, the corresponding area is divided into 300 m × 200 m medium areas 401, and the divided medium area 401 is divided into 60 m × 40 m small areas 402. In the area dividing method B, the corresponding area is divided into 300 m × 200 m medium areas 401, and the divided medium area 401 is divided into 50 m × 50 m small areas 402.

  The newly arrived information input means 101 normally holds a region of the corresponding area whose one side is about 3 km to 10 km, and the average resolution (minimum interval of x values or minimum interval of y values) is 0.3 m or less. It is about 10m. In the area division method of the corresponding area set in advance in the area division control means 301, this whole area is divided into a middle area 401 having one side of about 100 m to 300 m, and this middle area 401 is further divided by 30 m on one side. Or it is divided into small areas 402 of about 100 m.

  As the area dividing method, the area dividing methods A and B are illustrated in FIG. 3, but other than this, a plurality of types of area dividing methods for dividing vertically and obliquely are prepared and set. Note that the length of one side of the small area takes into account the magnitude of events (new construction or removal of buildings, changes in riverbed, etc.) that should be detected as ground changes. In the case of crustal deformation due to a large-scale earthquake or the like, it can be changed and set in the area division control means 301.

  The area dividing unit 302 is based on the area dividing method A of the corresponding area set in the area dividing control unit 301, for example, with the corresponding area as a predetermined middle area and small area with the lowest end as the base point. Divide into

  The area statistics calculation unit 303 calculates the area statistics of the ground information for each predetermined middle area and small area divided by the area dividing unit 302 and stores the area statistics in the area statistics memory 304. Here, when the ground surface information is x, y, z, x is latitude, y is longitude, and z is altitude, the area statistics calculation unit 303 determines whether the middle area and each small area have the latitude x and longitude y. Based on the altitude z at each point, the arithmetic average value (hereinafter referred to as the average value) of the altitude z in each middle area and each small area is calculated as area statistics, and the result is stored in the area statistics memory 304.

  Next, the area dividing unit 302 shifts the latitude x, which is half the size of the small area, by 30 m from the base point of the corresponding area based on the area dividing method A of the corresponding area set in the area dividing control unit 301. However, the area statistics calculation means 303 divides the corresponding area into predetermined middle areas and small areas with the origin y at the place where the longitude y is shifted by 20 m and the position at which the latitude x is shifted by 30 m and the longitude y by 20 m. Similarly, for each predetermined medium area and small area divided by the area dividing means 302, based on the altitude z of each point in each latitude x and longitude y of each medium area and each small area, The average value of the altitude z in each small area is calculated as area statistics, and the result is stored in the area statistics memory 304.

  In addition, when the area dividing unit 302 divides the corresponding area from the base point of the corresponding area into a predetermined middle area and small area based on the area dividing method B of the corresponding area set in the area dividing control unit 301, and When the corresponding area is divided into a predetermined medium area and a small area with a point shifted by half the size of the small area from the base point of the area as a starting point, the area statistics calculation unit 303 performs the predetermined medium divided by the area dividing unit 302. Similarly, for each region and small region, based on the altitude z of each point in each latitude x and longitude y of each middle region and each small region, the average value of the height z in each middle region and each small region is region statistics. And the result is stored in the area statistics memory 304.

  Further, with respect to other area dividing methods set in the area division control unit 301, similarly, the area dividing unit 302 divides the area into predetermined medium areas and small areas, and calculates area statistics. In step 303, the area statistics of the ground surface information are calculated and stored in the area statistics memory 304.

The area statistics calculated by the area statistics calculation means 303 and stored in the area memory 304 will be described in relation to the surface information.
FIG. 4 shows the ground surface information as a schematic map, and shows a small area 402a and a small area 402b in which the longitude y of the small area 402a is shifted by half the size of the small area. There are a building 501 and a riverbed 502 in the surface information of FIG. The altitude z inside the building 501 is higher than the surrounding land, and the altitude z inside the riverbed 502 is lower than the surrounding land. The average value of the altitude z is higher in the small area 402b than in the adjacent small area where the building 501 and the river bed 502 are not present, if this vicinity is flat. In the small area 402a, since the average value of the altitude z is affected by the building 501 and the river bed 502, there is no particular feature when compared with a neighboring small area without the building 501 or the river bed 502, but the variance value is large.

  Therefore, the singular region calculation means 305 compares the region statistic of the middle region stored in the region statistic memory 304 with the region statistic for each small region inside the region statistic, and compared with a preset threshold value, A small area having the following characteristics is regarded as a unique area, and a unique reason code is set.

  FIG. 5 is a diagram showing an example of a singular reason code table preset in the singular region calculation means 305. This singular reason code table is set as a pair of a small region's region statistics state and a singular reason code, and the singular region calculation means 305 for the small region's region statistics state stored in the region statistics memory 304. Indicates a singular reason code given to the center coordinates (x, y) of the corresponding small area. FIG. 5 shows the detection purpose for explanation for each unique reason code.

  As shown in FIG. 5, there are singular reason codes 01, 02, 03, 04,... Based on average values as singular reason codes. For example, it is indicated that the singular reason code 01 is lower than the surroundings. The singular area calculation means 305 stores (x, y, singular reason code) as a singular area in the newly arrived singular information memory 104 together with the center coordinates (x, y) of the corresponding small area.

  That is, the singular area calculation means 305 compares the average value of the ground information of the small area calculated by the area statistics calculation means 303 with the average value of the ground information of the middle area including the small area, or the area statistics calculation means. If the average value of the ground surface information of the small area calculated in step 303 is compared with the average value of the ground information of the small area adjacent to the small area, and if the predetermined condition is set in advance, the small area is determined to be a singular area. Calculate as

  For example, as shown in FIG. 5, when the average value of the ground information of the small area is lower than the average value of the ground information of the middle area including the small area, the singular area calculation unit 305, Judgment reason code 01 is given to the small area judging that it is a depression or a riverbed, and the average value of the ground information of the small area is higher than the average value of the ground information of the middle area including the small area by a predetermined threshold or more. In this case, it is determined as a building or a mountain peak, and a special reason code 02 is given to the small area, and the average value of the ground information of the small area is lower than the average value of the ground information of the small area above (north) the small area. If it is higher than the average value of the ground information of the lower (south) small area, it is determined to be a slope facing south, and the special reason code 03 is given to the small area, and the average value of the ground information of the small area is Lower than the average value of the ground information of the left (west) small area, Is higher than the average value of the table information, gives a specific reason code 04 to the small region is determined that the East facing slopes.

  On the other hand, for the ground surface information stored in the traditional information memory 202, the unique local information detecting means 203 gives a unique reason code to the small area in the same manner as the unique local information detecting means 103, and the traditional unique information memory 204 Is stored as a singular region. Needless to say, this processing can be performed before newly obtained ground surface information is input to the new arrival information input means 101.

  After the singular reason code for each small area as the singular area is stored in the new arrival singular information memory 104 and the old singular information memory 204, the position correction means 105 is based on the singular reason code for each small area as the singular area. The newly obtained ground information stored in the new arrival information memory 102 and the previously obtained ground information stored in the old information memory 202 are aligned. If the position of the newly obtained ground information and the previously obtained ground information are not misaligned, the unique reason code in the newly arrived unique information memory 104 and the unique free code in the old unique information memory 204 are completely identical. However, since this is unlikely to occur based on actual observation data, alignment is performed as follows.

  The position correcting means 105 refers to (x1, y1, singular reason code) as the singular area stored in the newly arrived singular information memory 104, and (x2 as the singular area stored in the old singular information memory 204) , Y2, peculiar reason code), compare the singular reason code of each small region in the newly obtained middle region with the singular reason code of each small region in the middle region obtained in advance, and the same singular cause code When the number of small areas having a distance within a predetermined threshold is equal to or greater than the predetermined threshold, the positional deviation vector amount (x1-x2, y1-y2) in the middle area is As a candidate.

  Further, the position correcting means 105 performs the same processing by another region dividing method, for example, the region dividing method B shown in FIG. At this time, in order to surely obtain a positional deviation candidate, similar processing for obtaining a positional deviation candidate is performed for the middle area obtained by dividing the corresponding area by shifting the origin.

  In this way, the position correcting unit 105 obtains the positional deviation vector amount for each middle region as a positional deviation candidate in the middle region. At this time, if there is a positional deviation within 5 degrees in the scanning direction, the entire area has a different vector amount of positional deviation. For example, sin5 degrees = 0.087, and if one side is a middle area of 300 m, the positional deviation is different. Since the maximum of is about 2.6 m, the influence of the angle can be ignored if the shift is obtained for each middle region.

  Then, the position correction unit 105 selects the region dividing method for the corresponding region for which the maximum number of vector displacement amounts for the middle region is calculated, and the calculated position displacement is calculated for the middle region for which the amount of vector displacement is calculated. In the middle region where the vector amount is employed and the positional deviation vector amount is not calculated, the positional displacement vector amount of the middle region having the closest center coordinate is employed. Further, the position correction means 105 corrects the position by shifting newly obtained ground information stored in the newly arrived information memory 102 for each small area based on the adopted vector amount of the positional deviation for each middle area. The newly obtained ground surface information is stored in the newly arrived position correction information memory 106.

  Next, the ground surface information difference detection means 107 obtains in advance the ground surface information obtained by correcting the position stored in the newly arrived position correction information memory 106 and the corresponding position stored in the old information memory 202 in advance. Compared with the ground surface information, a significant difference in ground surface information is detected by the following method. Since there is a difference between sensors and a difference in height detection point (x, y), an altitude difference occurs even on a ground surface where there is essentially no difference. It's not a difference. Further, when calculating the positional deviation by the position correcting means 105, even if the corresponding area has the same singular reason code, since this is a comparison based on statistics, there can be a significant difference in ground information.

  The ground surface information difference detecting means 107 is the ground surface closest to the distance stored in the newly arrived position correction information memory 106 with respect to the values (x2, y2, z2) of the ground surface information in units of small areas of the ground surface information obtained in advance. The value (x1, y1, z1) of information is calculated | required and the difference of the height is calculated as z3 = z2-z1. The value (x2, y2, z3) is stored in the newly arrived information difference memory 108 for the altitude difference z3 in which the absolute value | z2-z1 | of the altitude difference exceeds a predetermined threshold value.

  What is stored in the new arrival information difference memory 108 is a list of points that have a highly significant difference for each point. Whether or not this is an actual ground surface change that is a detection purpose is determined by the following method. Evaluation is performed by the change evaluation means 109.

  The surface information change evaluation means 109 uses a value (x2, y2, z3) stored in the newly arrived information difference memory 108 within a predetermined distance (x2) having a positive altitude difference z3 larger than a predetermined threshold. , Y2) are traced from a small area near the base point of the corresponding area to form a counterclockwise chain, thereby forming an ascending area contour line and obtaining an area S1 of the outline area surrounded by the ascending area outline line. At the same time, a point (x2, y2) within a predetermined distance having a negative altitude difference z3 smaller than a predetermined threshold is traced from a small area near the base point of the corresponding area to form a counterclockwise chain. To form a descending area contour line, and obtain an area S2 of the contour region surrounded by the descending area contour line.

  Then, the surface information change evaluation means 109 evaluates that the rising area has actually occurred when the area S1 of the contour area surrounded by the rising area outline is larger than the predetermined parameter 1, and is surrounded by the falling area outline. If the area S2 of the contour region is larger than the predetermined parameter 1, it is evaluated that a descending area has actually occurred. Here, the parameter 1 is a parameter for eliminating an error in evaluating that an ascending zone or a descending zone has occurred.

  The surface information change evaluation means 109 obtains the circumscribed rectangle area S3 of the rising area outline, the circumscribed rectangle area S4 of the falling area outline, and the overlapping area S5 of both circumscribed rectangles, and S5 / S3> parameter 2 or S5. When / S4> parameter 2, it is evaluated that the rising area and the falling area are mixed. Here, the parameter 2 is a parameter for eliminating an error when evaluating that the rising area and the falling area are mixed.

  Then, the surface information change evaluation unit 109 stores in the change information memory 110 the values (x2, y2, z3) evaluated as actual changes in the surface information, or the ascending area outline and the descending area outline.

  The change information output unit 111 outputs the values (x2, y2, z3) stored in the change information memory 110, or the rising area outline and the falling area outline as change information.

  As described above, according to the first embodiment, the peculiar surface information detection unit 103 detects the peculiar distribution of the surface information of the corresponding area newly obtained and calculates the singular area of the corresponding area newly obtained. The singular ground information detection means 203 detects the singular distribution of the ground surface information of the corresponding area obtained in advance, calculates the singular area of the corresponding area obtained in advance, and the position correction means 105 obtains the newly obtained Based on the singular area of the corresponding area and the singular area of the corresponding area obtained in advance, the ground information of the corresponding area newly obtained is aligned with the ground information of the corresponding area obtained in advance, and the ground surface information difference detecting means 107 Detects the difference between the newly obtained ground information of the corresponding area and the ground information of the relevant area obtained in advance, and the ground information change evaluation means 109 is preset for the detected difference. The change in surface information Compared to the parameters for, by evaluating whether the change in actual generated surface information, the effect of the surface changes in those areas can be accurately analyzed quickly is obtained.

  In Embodiment 1 described above, altitude information is used as the z information of (x, y, z) as the ground surface information. Instead, the reflection value in the infrared sensor is a discrete value (for example, 0 to 255). This process can also be applied to In that case, instead of altitude, changes are evaluated in areas with high and low temperature radiation. Further, this processing can also be applied when color spectrum information based on the three primary colors red, blue, and green is obtained as discrete values instead of altitude information. In that case, the change is first evaluated by reflection by each primary color, and the change is determined by taking the logical sum of the evaluation results for each primary color by the ground surface information change evaluation means 109.

  The ground surface change determination apparatus according to the first embodiment can be realized by a computer and a ground surface change determination program for causing the computer to function as each unit according to the first embodiment.

It is a block diagram which shows the structure of the surface information change evaluation apparatus by Embodiment 1 of this invention. It is a block diagram which shows the structure of the specific surface information detection means of the surface information change evaluation apparatus by Embodiment 1 of this invention. It is a figure which shows the example of the area division method of the applicable area preset by the area division control means in the specific surface information detection means of the surface information change evaluation apparatus by Embodiment 1 of this invention. In the ground surface information change evaluation apparatus according to Embodiment 1 of the present invention, the ground surface information is represented as a schematic map, and a small region and a small region obtained by shifting the longitude of the small region by half the size of the small region. It is a figure which shows the example of the peculiar reason code table preset by the peculiar area calculation means in the peculiar surface information detection means of the surface information change evaluation apparatus which concerns on Embodiment 1 of this invention.

Explanation of symbols

  101 newly arrived information input means, 102 newly arrived information memory, 103 unique land surface information detecting means, 104 newly arrived unique information memory, 105 position correcting means, 106 newly arrived position correction information memory, 107 ground surface information difference detecting means, 108 newly arrived information difference memory, 109 Land surface information change evaluation means, 110 Change information memory, 111 Change information output means, 202 Conventional information memory, 203 Singular ground surface information detection means, 204 Conventional singular information memory, 208 Ground surface information change model memory, 301 Area division control means, 302 area Dividing means, 303 area statistics calculating means, 304 area statistics memory, 305 singular area calculating means, 401 middle area, 402 small area, 402a small area, 402b small area, 501 building, 502 riverbed.

Claims (2)

Detecting the unique distribution of the surface information of the corresponding area newly obtained and calculating the unique area of the corresponding area newly obtained, and detecting the unique distribution of the surface information of the corresponding area obtained in advance A peculiar surface information detecting means for calculating a peculiar region of the obtained area;
Based on the singular area of the corresponding area newly obtained and the singular area of the corresponding area obtained in advance calculated by the above-mentioned singular ground information detection means, the ground information of the corresponding area newly obtained and the surface of the corresponding area obtained in advance Position correction means for aligning with information;
A ground surface information difference detecting means for detecting a difference between the ground information of the corresponding area newly obtained by the position correction means and the ground surface information of the corresponding area obtained in advance;
The ground surface information for evaluating whether or not the difference detected by the ground surface information difference detecting means is a change in the ground surface information actually generated by comparing with a preset parameter for evaluating the ground surface information change. Bei example and change evaluation means,
The singular ground surface information detecting means is
Area dividing means for dividing the area into predetermined medium and small areas;
Area statistics calculating means for calculating area statistics of the ground information for each predetermined medium area and small area divided by the area dividing means;
Singular area calculation means for detecting a specific distribution of ground information of the corresponding area and calculating a singular area of the corresponding area on the basis of the area statistics of the ground information for each medium area and small area calculated by the area statistics calculating means; A ground surface change determination device characterized by comprising: The area statistics calculation means calculates the average value of the surface information of each middle area and each small area as area statistics based on the ground surface information of each point of each middle area and each small area,
The singular area calculation means compares the average value of the ground information of the small area calculated by the area statistics calculation means with the average value of the ground information of the middle area including the small area, or by the area statistics calculation means. Compare the calculated average value of the ground surface information of the small area with the average value of the ground surface information of the small area adjacent to the small area, and calculate the small area as a singular area when a predetermined condition is set in advance. The ground surface change determination device according to claim 1 .

Images