Building height change identification method, device and equipment based on remote sensing spatio-temporal information and storage medium

Abstract

The invention discloses a building height change identification method, device and equipment based on remote sensing spatio-temporal information, and a storage medium. The method comprises the following steps: acquiring a building bottom contour by using a vector map; building top contours at different moments are obtained through the multi-temporal high-resolution optical remote sensing images; according to the Euclidean distance between the geometric center point of the bottom contour of the building and the geometric center point of the top contour of the building, obtaining projection displacement of the top and bottom contours of the building; calculating a multi-building relative projection displacement ratio and a single-building absolute projection displacement ratio of different buildings according to the top and bottom contour projection displacements of the buildings at different moments; calculating a satellite inclination angle tangent ratio according to the satellite observation angles at different moments; and according to the multi-building relative projection displacement ratio, the single-building absolute projection displacement ratio and the satellite dip angle tangent ratio, identifying the height change of the building. According to the invention, the accuracy of building height change identification can be improved.

Description

Technical Field

[0001] The present application relates to the field of image data processing, and in particular to a method, device, equipment and storage medium for identifying building height changes based on remote sensing spatiotemporal information. Background Art

[0002] The increase in building height may bring a series of safety hazards, such as uneven foundation settlement, load-bearing structure damage, cracks, deformation, etc., which will directly affect the structural safety and stability of the building. In the method of identifying building height changes, domestic and foreign experts and scholars generally perform shadow inversion based on the building shadow characteristics of remote sensing images to derive the relationship between shadow length and building height. However, in the scene of dense urban buildings, factors such as shadow overlap, confusion, changes in sun angle, and building geometry often cause certain limitations to this method. Summary of the invention

[0003] The main purpose of this application is to provide a method, device, equipment and storage medium for identifying building height changes based on remote sensing spatiotemporal information, aiming to solve the technical problem of how to improve the accuracy of identifying height changes of single buildings in dense urban building complexes.

[0004] To achieve the above objectives, the present application proposes a method for identifying building height changes based on remote sensing spatiotemporal information, the method comprising:

[0005] Use vector maps to obtain the bottom outline of buildings;

[0006] Use multi-temporal high-resolution optical remote sensing images to obtain the top contours of buildings at different times;

[0007] Obtaining the projection displacement of the top and bottom contours of the building according to the Euclidean distance between the geometric center point of the bottom contour of the building and the geometric center point of the top contour of the building;

[0008] Calculate the relative projection displacement ratio of multiple buildings and the absolute projection displacement ratio of a single building of different buildings according to the projection displacement of the top and bottom contours of the buildings at different times;

[0009] Calculating the satellite inclination tangent ratio according to the satellite observation angles at different times;

[0010] The height change of the building is identified according to the relative projection displacement ratio of multiple buildings, the absolute projection displacement ratio of the single building and the satellite inclination tangent ratio.

[0011] In one embodiment, the step of calculating the relative projection displacement ratio of multiple buildings and the absolute projection displacement ratio of a single building of different buildings according to the projection displacement of the top and bottom contours of the buildings at different times comprises:

[0012] Selecting a first building and a second building from the buildings;

[0013] According to the projection displacement of the top and bottom contours of the building, a first projection displacement of the first building and a second projection displacement of the second building at a first moment, a third projection displacement of the first building and a fourth projection displacement of the second building at a second moment, a fifth projection displacement of the first building at a third moment, and a sixth projection displacement of the first building at a fourth moment are obtained;

[0014] Obtaining a relative projection displacement ratio of multiple buildings according to the first projection displacement, the second projection displacement, the third projection displacement and the fourth projection displacement;

[0015] According to the fifth projection displacement and the sixth projection displacement, the absolute projection displacement ratio of a single building is obtained.

[0016] In one embodiment, the step of obtaining the relative projection displacement ratio of multiple buildings according to the first projection displacement, the second projection displacement, the third projection displacement and the fourth projection displacement comprises:

[0017] Obtain a first relative projection displacement ratio according to the first projection displacement and the second projection displacement, and obtain a second relative projection displacement ratio according to the third projection displacement and the fourth projection displacement;

[0018] The first relative projection displacement ratio and the second relative projection displacement ratio are used as the multi-building relative projection displacement ratio.

[0019] In one embodiment, the step of calculating the satellite inclination tangent ratio according to the satellite observation angles at different times comprises:

[0020] Acquire a first observed inclination angle of a satellite relative to a first building at a third moment, and a second observed inclination angle of the satellite relative to the first building at a fourth moment;

[0021] The satellite inclination tangent ratio is obtained according to the first observed inclination angle and the second observed inclination angle.

[0022] In one embodiment, the step of identifying the height change of the building according to the relative projection displacement ratio of multiple buildings, the absolute projection displacement ratio of the single building and the satellite inclination tangent ratio comprises:

[0023] Obtaining a first relative projection displacement ratio and a second relative projection displacement ratio according to the relative projection displacement ratios of the multiple buildings;

[0024] Calculate a first ratio of the first relative projection displacement ratio and the second relative projection displacement ratio;

[0025] When the first ratio satisfies a first criterion and the single building absolute projection displacement ratio is equal to the satellite inclination tangent ratio, determining that the height of the building has not changed;

[0026] When the first ratio does not meet the first standard or the single building absolute projection displacement ratio is not equal to the satellite inclination tangent ratio, it is determined that the height of the building has changed.

[0027] In one embodiment, when the first ratio does not meet the first standard or the single building absolute projection displacement ratio is not equal to the satellite inclination tangent ratio, the step of determining that the height of the building has changed includes:

[0028] Obtaining the building height change difference according to the relative projection displacement ratio of the multiple buildings and / or the absolute projection displacement ratio of the single building;

[0029] When the height change difference of the building is not greater than the height change threshold, the building is included in the observation list;

[0030] When the building height change difference is greater than the height change threshold, the building is included in the danger list.

[0031] In one embodiment, the step of acquiring the top contour of the building at different times using multi-temporal high-resolution optical remote sensing images comprises:

[0032] Using multi-temporal high-resolution optical remote sensing images, we can obtain images of building tops at different times;

[0033] The image of the building top is processed based on an image segmentation algorithm to obtain a building top outline.

[0034] In addition, to achieve the above purpose, the present application also proposes a device for identifying building height changes based on remote sensing spatiotemporal information, the device comprising:

[0035] A data acquisition module, used for acquiring the bottom contour of the building using a vector map;

[0036] The data acquisition module is also used to acquire the top contours of buildings at different times using multi-temporal high-resolution optical remote sensing images;

[0037] A data calculation module, used to obtain the projection displacement of the top and bottom contours of the building according to the Euclidean distance between the geometric center point of the bottom contour of the building and the geometric center point of the top contour of the building;

[0038] The data calculation module is further used to calculate the relative projection displacement ratio of multiple buildings and the absolute projection displacement ratio of a single building of different buildings according to the projection displacement of the top and bottom contours of the buildings at different times;

[0039] The data calculation module is further used to calculate the satellite inclination tangent ratio according to the satellite observation angles at different times;

[0040] The identification and analysis module is used to identify the height change of the building according to the relative projection displacement ratio of multiple buildings, the absolute projection displacement ratio of the single building and the satellite inclination tangent ratio.

[0041] In addition, to achieve the above-mentioned purpose, the present application also proposes a device for identifying building height changes based on remote sensing spatiotemporal information, the device comprising: a memory, a processor, and a computer program stored on the memory and executable on the processor, the computer program being configured to implement the steps of the method for identifying building height changes based on remote sensing spatiotemporal information as described above.

[0042] In addition, to achieve the above-mentioned purpose, the present application also proposes a storage medium, which is a computer-readable storage medium, and a computer program is stored on the storage medium. When the computer program is executed by the processor, the steps of the method for identifying building height changes based on remote sensing spatiotemporal information as described above are implemented.

[0043] In addition, to achieve the above-mentioned purpose, the present application also provides a computer program product, which includes a computer program, and when the computer program is executed by a processor, it implements the steps of the method for identifying building height changes based on remote sensing spatiotemporal information as described above.

[0044] One or more technical solutions proposed in this application have at least the following technical effects:

[0045] By comparing the projection displacement of buildings at different times based on their geometric projection displacement in remote sensing images, the method can realize the identification of building height changes, solve the technical problem that building height identification in dense urban building complexes is affected by factors such as shadow overlap, confusion, changes in sun angle, and building geometric shape. Compared with the existing technology, the method can improve the accuracy of identifying height changes of single buildings in dense urban building complexes. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application.

[0047] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, for ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative labor.

[0048] Figure 1 A flow chart of a first embodiment of a method for identifying building height changes based on remote sensing spatiotemporal information provided by the present application;

[0049] Figure 2 A flow chart of a second embodiment of a method for identifying building height changes based on remote sensing spatiotemporal information provided by the present application;

[0050] Figure 3 A schematic diagram of the projection displacement of two buildings in two remote sensing images of a building height change identification method based on remote sensing spatiotemporal information provided in Example 2 of the present application;

[0051] Figure 4 A schematic diagram of the projection displacement of a single building in two phases of remote sensing images of a building height change identification method based on remote sensing spatiotemporal information provided in Example 2 of the present application;

[0052] Figure 5 A flowchart diagram of a third embodiment of a method for identifying building height changes based on remote sensing spatiotemporal information provided by the present application;

[0053] Figure 6 A flowchart diagram of a fourth embodiment of a method for identifying building height changes based on remote sensing spatiotemporal information provided by the present application;

[0054] Figure 7 A brief flowchart of a building height change identification method based on remote sensing spatiotemporal information in this application;

[0055] Figure 8 A flowchart diagram of a fifth embodiment of a method for identifying building height changes based on remote sensing spatiotemporal information provided by the present application;

[0056] Fig. 9 This is a schematic diagram of the module structure of a device for identifying building height changes based on remote sensing spatiotemporal information according to an embodiment of the present application;

[0057] Fig.10 Schematic diagram of the equipment structure of the hardware operating environment involved in the method for identifying building height changes based on remote sensing spatiotemporal information in an embodiment of the present application.

[0058] The purpose, features and advantages of this application will be further described in conjunction with the embodiments and with reference to the accompanying drawings. DETAILED DESCRIPTION

[0059] It should be understood that the specific embodiments described herein are only used to explain the technical solutions of the present application and are not used to limit the present application.

[0060] In order to better understand the technical solution of the present application, a detailed description will be given below in conjunction with the accompanying drawings and specific implementation methods.

[0061] The main solution of the embodiment of the present application is: using a vector map to obtain the bottom contour of the building; using multi-phase high-resolution optical remote sensing images to obtain the top contour of the building at different times; according to the Euclidean distance between the geometric center point of the bottom contour of the building and the geometric center point of the top contour of the building, the projection displacement of the top and bottom contours of the building is obtained; according to the projection displacement of the top and bottom contours of the building at different times, the relative projection displacement ratio of multiple buildings and the absolute projection displacement ratio of a single building of different buildings are calculated; according to the satellite observation angles at different times, the satellite inclination tangent ratio is calculated; according to the relative projection displacement ratio of multiple buildings, the absolute projection displacement ratio of a single building and the satellite inclination tangent ratio, the height change of the building is identified.

[0062] In this embodiment, for the convenience of description, the following description is made by taking the internal actuator of the height change recognition system as the execution subject.

[0063] In the existing technology, the building height recognition of dense urban building complexes is affected by factors such as shadow overlap, confusion, changes in sun angles, and building geometry.

[0064] The present application provides a solution, which uses the geometric projection displacement of the building in the remote sensing image to compare the projection displacement of the building at different times to realize the identification of the change of building height, thereby improving the accuracy of the identification of the height change of a single building in a dense urban building complex.

[0065] It can be seen from the above embodiments that the present application realizes the recognition of building height changes by comparing the projection displacement of the building at different times according to the geometric projection displacement of the building in the remote sensing image, thereby solving the technical problem that the building height recognition in dense urban building complexes is affected by factors such as shadow overlap, confusion, change in sun angle, and building geometric shape, and achieves the improvement of the accuracy of identifying the height changes of single buildings in dense urban building complexes.

[0066] It should be noted that the execution subject of this embodiment can be a computing service device with data processing, network communication and program running functions, such as a tablet computer, a personal computer, a mobile phone, etc., or an electronic device capable of realizing the above functions, etc. The following takes the internal actuator of the height change recognition system as an example to illustrate this embodiment and the following embodiments.

[0067] Based on this, the embodiment of the present application provides a method for identifying building height changes based on remote sensing spatiotemporal information, referring to Figure 1 , Figure 1 This is a flow chart of the first embodiment of the method for identifying building height changes based on remote sensing spatiotemporal information of the present application.

[0068] In this embodiment, the building height change identification method based on remote sensing spatiotemporal information includes steps S10 to S60:

[0069] Step S10, using the vector map to obtain the bottom outline of the building.

[0070] It should be noted that buildings are fixed structures built by people to meet the needs of activities such as living, working, studying, production, and entertainment, including residential buildings, such as single-story / multi-story houses, public buildings, such as government office buildings, libraries, hospitals, etc., commercial buildings, such as shopping centers, industrial buildings, such as factories, transportation buildings, such as railway stations, agricultural buildings, such as greenhouses, scientific research buildings, such as observatories, historical buildings, such as temples, memorial buildings, such as monuments, and mixed-function buildings, which are a collection of the above-mentioned various types of buildings; a large number of buildings gathered in a specific area of ​​a city form a dense urban building complex, which is closely linked in spatial organization and is an important carrier of urban functions and urban culture. The following takes a dense urban building complex as an example to illustrate this embodiment and the following embodiments.

[0071] In addition, it should be noted that a vector map is a type of map that uses mathematical algorithms to represent geographic spatial data. Compared with bitmaps (such as common JPEG, PNG and other formats), vector maps have the characteristics of infinite scalability, strong interactivity, small data volume, dynamic rendering, etc., and are widely used in navigation software, urban planning, geographic information systems, online map services and other fields. The vector map in this embodiment is vector map data provided by the National Geographic Information Public Service Platform (TianDiTu); the bottom outline of a building refers to the data in which the geometric shape of the bottom boundary of a building is represented in vector form in a computer, which can accurately describe the outline of the building and can be scaled without distortion.

[0072] Step S20, using multi-temporal high-resolution optical remote sensing images to obtain the top contours of the buildings at different times.

[0073] It should be noted that multi-phase high-resolution optical remote sensing images refer to surface images acquired through high-resolution optical remote sensing satellites, which can cover multiple perspectives or multiple scenes and have the characteristics of high resolution, multi-perspective, multi-mode imaging, and stereoscopic imaging; the top outline of a building refers to the data in which the geometric shape of the top boundary of the building is represented in vector form in a computer.

[0074] Step S30, obtaining the projection displacement of the top and bottom contours of the building according to the Euclidean distance between the geometric center point of the bottom contour of the building and the geometric center point of the top contour of the building.

[0075] It should be noted that the projection displacement of the top and bottom contours of the building is the modulus of the projection offset vector of the top and bottom contours of the building on the image plane during satellite imaging, which is obtained based on the geometric projection characteristics of the building in the remote sensing image. The geometric projection characteristics of the building in the remote sensing image, that is, the top geometric contour of the building and the bottom geometric contour of the building have a certain horizontal displacement in the image, that is, the projection displacement. When the satellite sensor shoots the earth's surface in a non-vertical direction, the relative position of the top and bottom contours of the building will be displaced due to the height. Euclidean distance is the straight-line distance between two points in Euclidean space. It is a distance metric in metric space. It is widely used in geometry, physics, engineering, data science and other fields. In a two-dimensional plane, if the coordinates of point A are (x 1 ,y 1 ), the coordinates of point B are (x 2 ,y 2 ), then the Euclidean distance between A and B can be calculated by the following formula:

[0076]

[0077] According to the top contour data and bottom contour data of the building, the vector displacement of the center points of the top contour and the bottom contour is calculated, and the magnitude of the vector displacement is used as the top and bottom contour projection displacement of the building. The top and bottom contour projection displacement of the building is a scalar data, and its value only represents the distance of the top and bottom geometric projections of the building, and does not involve vector meanings such as direction.

[0078] The size of the projection displacement of the top and bottom contours of the building is related to the height of the building and the viewing angle of the satellite. The calculation formula is as follows:

[0079]

[0080] Among them, (x t ,y t ) is the coordinate of the geometric center point of the top outline of the building, (x b ,y b ) is the coordinate of the geometric center point of the bottom contour of the building, H is the building height, D is the plane projection displacement between the top and bottom contours, that is, the projection displacement of the top and bottom contours of the building, and θ is the inclination angle between the satellite and the building, that is, the angle between the line from the satellite to the top and bottom of the building and the horizontal line.

[0081] Step S40, calculating the relative projection displacement ratio of multiple buildings and the absolute projection displacement ratio of a single building according to the projection displacement of the top and bottom contours of the buildings at different times.

[0082] It should be noted that the relative projection displacement ratio of multiple buildings is the data used to compare multiple buildings and identify their height changes, and it is obtained by calculating the ratio of the projection displacements of the top and bottom contours of multiple buildings; the absolute projection displacement ratio of a single building is the data used to compare the projection displacements of the top and bottom contours of a single building at different times, and identify its height changes by calculating the ratio of the projection displacements.

[0083] Step S50, calculating the satellite inclination tangent ratio according to the satellite observation angles at different times.

[0084] It should be noted that the satellite observation angle, that is, the observation inclination of the satellite relative to the building, refers to the viewing angle and orientation of the satellite relative to the target on the earth's surface when performing remote sensing observation, and also refers to the angle between the point directly below the satellite and the observed target point on the earth's surface. This angle measures the angle of the optical remote sensing satellite's tilt shooting, which has an important impact on the quality of the satellite image and the target characteristics. In this embodiment, the satellite observation angle is the viewing angle and orientation of the satellite relative to the building when shooting the remote sensing image of the building; the satellite inclination tangent ratio is the ratio of the tangent of the satellite observation angle of a building at two different times.

[0085] Step S60, identifying the height change of the building according to the relative projection displacement ratio of multiple buildings, the absolute projection displacement ratio of the single building and the satellite inclination tangent ratio.

[0086] The relative projection displacement ratio of multiple buildings is used to compare the projection displacements of the top and bottom contours of multiple buildings, and the absolute projection displacement ratio of a single building and the satellite inclination tangent ratio are used to compare the projection displacements of the top and bottom contours of a single building.

[0087] When identifying the heights of multiple buildings, the identification results include that only one building height has changed, or that multiple buildings heights have changed at the same time, or that multiple buildings heights have not changed; when a building height changes, it is necessary to combine the identification results of the height change of a single building to further determine which one or several buildings have changed in height, and finally obtain the specific results of the building height change.

[0088] This embodiment provides a method for identifying building height changes based on remote sensing spatiotemporal information. It uses the geometric projection displacement of the building in the remote sensing image to compare the projection displacement of the building at different times to realize the identification of building height changes. It solves the technical problem that in dense urban building complexes, building height identification is affected by factors such as shadow overlap, confusion, changes in sun angle, and building geometric shape, and achieves the improvement of the accuracy of identifying height changes of single buildings in dense urban building complexes.

[0089] Based on the first embodiment of the present application, in the second embodiment of the present application, the same or similar contents as those in the above-mentioned embodiment 1 can be referred to the above introduction, and will not be repeated in the following. Figure 2 , step S40 includes steps S41 to S44:

[0090] Step S41, selecting a first building and a second building from the buildings.

[0091] Select a building from the building set as the target body for height change recognition. Two, three or more buildings can be selected for recognition. It is recommended that no more than five buildings should be selected. When recognizing height changes of multiple buildings, the relative projection displacement ratio of multiple buildings can be calculated in groups by grouping two buildings. Based on the selected buildings, the height changes of these buildings are recognized according to their projection displacements. The following takes the selection of two buildings for height change recognition as an example to illustrate this embodiment and the following embodiments.

[0092] From the buildings in the dense urban building complex, two buildings with a moderate distance, for example, a straight-line distance greater than fifty meters, are selected as the first building and the second building, respectively, which can also be called building A and building B. The names are only used to distinguish the two buildings.

[0093] When selecting buildings, buildings that are too close have unclear distance boundaries in the remote sensing image, which may cause their projections to overlap, resulting in calculation errors in the projection displacement. For buildings that are too far away, when the satellite collects remote sensing images, the angle of the satellite relative to the building is quite different. Since the projection displacement of the top and bottom contours of the building is related to the height of the building and the viewing angle of the satellite, it will lead to calculation errors in the projection displacement.

[0094] Step S42, according to the projection displacement of the top and bottom contours of the building, obtain the first projection displacement of the first building and the second projection displacement of the second building at the first moment, the third projection displacement of the first building and the fourth projection displacement of the second building at the second moment, the fifth projection displacement of the first building at the third moment, and the sixth projection displacement of the first building at the fourth moment.

[0095] According to the top and bottom contour projection displacement of the building set, we can obtain the displacement at the first moment, which can also be called T i During the period, the projection displacements of the first building (building A) and the second building (building B) are respectively taken as the first projection displacement D A and the second projection displacement D B , get the second moment, also called T i+1During the period, the projection displacements of the first building (building A) and the second building (building B) are respectively used as the third projection displacement D A ′ and the fourth projection displacement D ′ B , get the value at the third moment T m and the fourth moment T n The projection displacement of the first building (building A) is respectively used as the fifth projection displacement D 1 and the sixth projection displacement D 2 .

[0096] When identifying the height change of two buildings, the third moment can be the same as or different from the first moment, and the fourth moment can be the same as or different from the second moment. If they are the same, then the first projection displacement D A Displacement D with the fifth projection 1 Same, the third projection displacement D A ′ Displacement D with the sixth projection 2 If they are different, the time difference between the first moment and the third moment, and between the second moment and the fourth moment should be controlled within an appropriate range, for example, the time difference is 1 year, and the longest should not exceed 5 years. In practical applications, the second moment and the fourth moment are generally close to the time for height change recognition detection, so as to avoid the data at the moment used for height change recognition losing its reference value due to the large time difference. The following takes the first moment and the third moment, and the second moment and the fourth moment being different as an example to illustrate this embodiment and the following embodiments.

[0097] The schematic diagram of the projection displacement of two buildings in the two phases of remote sensing images is as follows Figure 3 As shown in the figure, the schematic diagram of the projection displacement of a single building in the two phases of remote sensing images is as follows Figure 4 shown.

[0098] Step S43: obtaining a relative projection displacement ratio of multiple buildings according to the first projection displacement, the second projection displacement, the third projection displacement and the fourth projection displacement.

[0099] It should be noted that the relative projection displacement ratio of multiple buildings is the ratio of the projection displacements of two buildings in remote sensing images in the same period.

[0100] According to the remote sensing image at a certain moment, the projection displacement of the top and bottom geometric contours of the two buildings is calculated, and then the projection displacements of the two buildings are compared to obtain the relative projection displacement ratio K of the group of buildings, that is, according to the first projection displacement D A , the second projection displacement D B , the third projection displacement D ′ A and the fourth projection displacement D′ B , the relative projection displacement ratio K of multiple buildings is calculated.

[0101] In a feasible implementation manner, step S43 includes steps S431 to S432:

[0102] Step S431: obtaining a first relative projection displacement ratio according to the first projection displacement and the second projection displacement, and obtaining a second relative projection displacement ratio according to the third projection displacement and the fourth projection displacement.

[0103] In T i At this moment, assuming that the projection displacements between the top and bottom of the first building (building A) and the second building (building B) are D A and D B , then the relative projection displacement ratio of the two is That is, the first relative projection displacement ratio, which can be expressed by the following formula:

[0104]

[0105] In T i+1 At the moment, the projection displacement between the top and bottom of the first building (building A) and the second building (building B) is D′ A and D′ B , then the relative projection displacement ratio of the two is That is, the second relative projection displacement ratio, which can be expressed by the following formula:

[0106]

[0107] Step S432: taking the first relative projection displacement ratio and the second relative projection displacement ratio as a multi-building relative projection displacement ratio.

[0108] After calculating the relative projection displacement ratios of the first building and the second building at two moments, the first relative projection displacement ratio and the second relative projection displacement ratio are used as the relative projection displacement ratios of multiple buildings for subsequent identification of building height changes based on the relative projection displacement ratios of multiple buildings.

[0109] By calculating the relative projection displacement ratio of two buildings, a data basis is provided for height change identification.

[0110] Step S44: obtaining a single building absolute projection displacement ratio according to the fifth projection displacement and the sixth projection displacement.

[0111] It should be noted that the absolute projection displacement ratio of a single building is the ratio of the projection displacement of a single building in remote sensing images in two periods.

[0112] In the two remote sensing images before and after, a building is selected, and the projection displacement of the geometric contours of the top and bottom of the building are calculated respectively. Then, the projection displacement of the building at the two moments is compared to obtain the absolute projection displacement ratio R of the building.

[0113] In T m At the moment, the projection displacement between the top and bottom of the first building (building A) is D 1 ; in T n At the moment, the projection displacement between the top and bottom of the first building (building A) is D 2 , the absolute projection displacement ratio R of the first building (building A) at two moments 2 , that is, the first absolute projection displacement ratio, can be expressed by the following formula:

[0114]

[0115] This embodiment provides a method for identifying building height changes based on remote sensing spatiotemporal information, which includes selecting a first building and a second building from buildings; obtaining a first projection displacement of the first building and a second projection displacement of the second building at a first moment, a third projection displacement of the first building and a fourth projection displacement of the second building at a second moment, a fifth projection displacement of the first building at a third moment, and a sixth projection displacement of the first building at a fourth moment according to the projection displacement of the top and bottom contours of the buildings; obtaining a relative projection displacement ratio of multiple buildings according to the first projection displacement, the second projection displacement, the third projection displacement, and the fourth projection displacement; obtaining an absolute projection displacement ratio of a single building according to the fifth projection displacement and the sixth projection displacement, thereby providing a data basis for height change identification.

[0116] Based on the first embodiment of the present application, in the third embodiment of the present application, the same or similar contents as those in the above-mentioned embodiment 1 can be referred to the above introduction, and will not be repeated in the following. Figure 5 , step S50 includes steps S51-S52:

[0117] Step S51, obtaining a first observation inclination angle of a satellite relative to a first building at a third moment and a second observation inclination angle of the satellite relative to the first building at a fourth moment.

[0118] At the third moment T m , obtain the observation inclination angle of the satellite relative to the first building (building A) at this time as the first observation inclination angle θ m , at the fourth moment T n , obtain the observation inclination angle of the satellite relative to the first building (building A) at this time as the second observation inclination angle θ n .

[0119] Step S52, obtaining the satellite inclination tangent ratio according to the first observed inclination and the second observed inclination.

[0120] In T m At time θ, the satellite's observation angle relative to the first building (building A) is m ; in T n At time θ, the satellite's observation angle relative to the first building (building A) is n , the tangent ratio of the observed inclination angle of the first building (building A) at the two moments is R 2 , which can be expressed by the following formula:

[0121]

[0122] This embodiment provides a method for identifying building height changes based on remote sensing spatiotemporal information, obtaining a first observed inclination angle of a satellite relative to a first building at a third moment and a second observed inclination angle of the satellite relative to the first building at a fourth moment; the satellite inclination tangent ratio is obtained based on the first observed inclination angle and the second observed inclination angle, thereby providing a data basis for height change identification.

[0123] Based on the first embodiment of the present application, in the fourth embodiment of the present application, the same or similar contents as those in the first embodiment can be referred to the above description, and will not be described in detail later. Figure 6 , step S60 includes steps S61 to S64:

[0124] Step S61, obtaining a first relative projection displacement ratio and a second relative projection displacement ratio according to the relative projection displacement ratios of the multiple buildings.

[0125] A first relative projection displacement ratio and a second relative projection displacement ratio are obtained from the relative projection displacement ratios of multiple buildings.

[0126] It should be noted that the relative projection displacement ratio describes the relative change between the projection positions of a building at two different time points. It compares the projection positions of a building at two different time points to determine whether the building has shifted and the direction and degree of the displacement. This comparison is usually made between images at two time points. For example, by comparing two satellite images or aerial photographs of the same building taken at different times, it can be determined whether the building has changed in height.

[0127] Step S62, calculating a first ratio of the first relative projection displacement ratio and the second relative projection displacement ratio.

[0128] Calculate the first ratio X of the first relative projection displacement ratio and the second relative projection displacement ratio 1 , which can be expressed by the following formula:

[0129]

[0130] Since D = H × tanθ, X 1 It can also be expressed as the following formula:

[0131]

[0132] Among them, H 1 T i The height of the first building (building A) at the moment, H 2 T i The height of the second building (building B) at this moment, H ′ 1 T i+1 The height of the first building (building A) at the moment, H ′ 2 T i+1 The height of the second building (building B) at this moment, θ A is the inclination angle between the satellite and the first building (building A), θ 1 is the inclination angle between the satellite and the second building (building B). At different times, the error of the satellite's observed inclination angle relative to the building is very small. Therefore, when identifying height changes in multiple buildings, the inclination error is negligible. The satellite's observed inclination angle relative to the same building at two times is expressed as the same angle.

[0133] Step S63, when the first ratio meets the first standard and the single building absolute projection displacement ratio is equal to the satellite inclination tangent ratio, determining that the height of the building has not changed.

[0134] It should be noted that the first standard is the first ratio X 1 The error threshold standard, for example, |X 1 The value of -1| does not exceed 0.1. This value can be set according to the specific circumstances of building height change identification. When the original height or foundation height of the building whose height change needs to be identified is low, this value can be appropriately increased, for example, increased to 0.2. When the original height or foundation height of the building whose height change needs to be identified is high, or the building is a key observation and protection building, this value can be appropriately reduced, for example, reduced to 0.05.

[0135] For example, take the first standard as the first ratio X 1 The error does not exceed 0.1 in the first ratio X 1 When the first criterion is met, that is, |X 1 When -1|≤0.1, it means X 1 The magnitude of is close to 1, where:

[0136]

[0137] This indicates that the relative heights of the first building (building A) and the second building (building B) at different times have not changed. To determine whether the absolute heights of the first building (building A) and the second building (building B) have changed, we continue to calculate the absolute heights of the first building (building A) and the second building (building B) according to the second ratio X. 2 Make a judgment.

[0138] The ratio of the absolute projection displacement ratio of a single building to the tangent ratio of the satellite inclination angle is X 2 , which can be expressed by the following formula:

[0139]

[0140] Since D = H × tanθ, X 2 It can also be expressed as the following formula:

[0141]

[0142] Among them, H m T m The height of the first building (building A) at the moment, H n T n The height of the first building (Building A) at the moment.

[0143] When the absolute projection displacement ratio of a single building is equal to the tangent ratio of the satellite inclination angle, that is, the second ratio X 2 =H m / H n =1, it means that the absolute height of the first building (building A) has not changed. Since the relative heights of the first building (building A) and the second building (building B) have not changed, it is determined that the heights of both the first building (building A) and the second building (building B) have not changed.

[0144] Step S64, when the first ratio does not meet the first standard or the single building absolute projection displacement ratio is not equal to the satellite inclination tangent ratio, determining that the height of the building has changed.

[0145] When the first ratio does not meet the first standard or the second ratio does not meet the second standard, that is, the first ratio X 1 The magnitude is not close to 1, or the second ratio X 2If it is not equal to 1, it means that the absolute height of at least one of the first building (building A) and the second building (building B) has changed. There are three cases in which the first ratio and the second ratio meet the conditions: the first ratio does not meet the first standard and the absolute projection displacement ratio of the single building is not equal to the satellite inclination tangent ratio; the first ratio meets the first standard and the absolute projection displacement ratio of the single building is not equal to the satellite inclination tangent ratio; the first ratio does not meet the first standard and the absolute projection displacement ratio of the single building is equal to the satellite inclination tangent ratio. The following takes the first standard as the first ratio X 1 Taking the error of no more than 0.1 as an example, the three cases are explained respectively.

[0146] In the first case, when the first ratio does not meet the first standard and the absolute projection displacement ratio of a single building is not equal to the tangent ratio of the satellite inclination, the first ratio X 1 Does not meet the first criterion, that is, |X 1 When -1|>0.1, it means X 1 The magnitude of is not close to 1, where:

[0147]

[0148] This means that the relative height of the first building (building A) and the second building (building B) has changed, that is, the absolute height of at least one of the first building (building A) and the second building (building B) has changed; and because the absolute projection displacement ratio of a single building is not equal to the tangent ratio of the satellite inclination, that is, the second ratio X 2 =H 1 / H ′ 1 ≠1, indicating that the absolute height of the first building (building A) has changed; to further determine whether the absolute height of the second building (building B) has changed, it is necessary to return to the step of obtaining the relative projection displacement ratio of multiple buildings and the absolute projection displacement ratio of a single building according to the top and bottom contour projection displacements, and continue to obtain the single building absolute projection displacement ratio of the second building (building B), and then determine whether the absolute height of the second building (building B) has changed.

[0149] In the second case, the first ratio meets the first criterion and the absolute projection displacement ratio of a single building is not equal to the tangent ratio of the satellite inclination angle, because the first ratio X 1 Satisfies the first criterion, that is, |X 1 When -1|≤0.1, it means X 1 The magnitude of is close to 1, where:

[0150]

[0151] This means that the relative heights of the first building (building A) and the second building (building B) have not changed; and because the absolute projection displacement ratio of a single building is not equal to the tangent ratio of the satellite inclination, that is, the second ratio X 2 =H 1 / H ′ 1 ≠1, indicating that the absolute height of the first building (building A) has changed, it can be determined that the absolute height of the second building (building B) has not changed; therefore, when the first ratio meets the first standard and the single building absolute projection displacement ratio is not equal to the satellite inclination tangent ratio, the absolute height of the first building (building A) has changed, and the absolute height of the second building (building B) has not changed.

[0152] In the third case, when the first ratio does not meet the first criterion and the absolute projection displacement ratio of a single building is equal to the tangent ratio of the satellite inclination, because the first ratio X 1 Does not meet the first criterion, that is, |X 1 When -1|>0.1, it means X 1 The magnitude of is not close to 1, where:

[0153]

[0154] This means that the relative height of the first building (building A) and the second building (building B) has changed, that is, the absolute height of at least one of the first building (building A) and the second building (building B) has changed; and because the absolute projection displacement ratio of a single building is equal to the tangent ratio of the satellite inclination, that is, the second ratio X 2 =H 1 / H ′ 1 =1, indicating that the absolute height of the first building (building A) has not changed, it can be determined that the absolute height of the second building (building B) has changed; therefore, when the first ratio does not meet the first standard and the single building absolute projection displacement ratio is equal to the satellite inclination tangent ratio, the absolute height of the first building (building A) has not changed, and the absolute height of the second building (building B) has changed.

[0155] In a feasible implementation manner, step S64 includes steps S01 to S03:

[0156] Step S01, obtaining building height change differences according to the relative projection displacement ratios of multiple buildings and / or the absolute projection displacement ratios of a single building.

[0157] It should be noted that the building height change difference is data obtained based on the absolute value of the difference between the relative projection displacement ratio of multiple buildings and / or the absolute projection displacement ratio of the single building and 1, and is represented by the symbol Δ.

[0158] When the first ratio satisfies the first standard and the absolute projection displacement ratio of a single building is equal to the tangent ratio of the satellite inclination angle, the absolute heights of the first building (building A) and the second building (building B) have not changed, and there is no need to calculate the difference in building height changes.

[0159] When the first ratio does not meet the first standard and the absolute projection displacement ratio of the single building is not equal to the satellite inclination tangent ratio, the building height change difference Δ of the first building (building A) 1 =|X 2 -1|, the building height change difference of the second building (building B) Δ 2 =|X 3 -1|, where X 2 is the ratio of the absolute projection displacement ratio of the first building (building A) to the tangent ratio of the satellite inclination angle, X 3 It is the ratio of the single building absolute projection displacement ratio of the second building (building B) and the satellite inclination tangent ratio.

[0160] When the first ratio satisfies the first standard and the absolute projection displacement ratio of the single building is not equal to the satellite inclination tangent ratio, the building height change difference Δ of the first building (building A) 1 =|X 1 -1|, since the absolute height of the second building (building B) has not changed at this time, the building height change difference of the second building (building B) is not calculated.

[0161] When the first ratio does not meet the first standard and the absolute projection displacement ratio of the single building is equal to the satellite inclination tangent ratio, the building height change difference Δ of the second building (building B) 2 =|X 3 -1|, since the absolute height of the first building (building A) has not changed at this time, the building height change difference of the first building (building A) is not calculated.

[0162] Step S02: When the height change difference of the building is not greater than the height change threshold, the building is included in the observation list.

[0163] It should be noted that the observation list is a list of buildings in densely populated urban building complexes whose building heights have changed, but the changes are small and the range of changes is within a safe range. The list includes the building's name, location information, height information, building type, pedestrian flow and other information.

[0164] In addition, it should be noted that the height change threshold is a critical value that limits the numerical value of the building height change difference. For example, the height change threshold is 0.1, which means that the building height change difference is limited to no less than 0.1, that is, Δ≤0.1. The size of the height change threshold can be changed according to the specific circumstances of the building height change identification, and can take any value from 0 to 1. In practical applications, it is generally recommended that the height change threshold should be less than 0.2. The following takes the height change threshold of 0.1 as an example to illustrate this embodiment and the following embodiments.

[0165] When the difference in building height change is not greater than the height change threshold, that is, Δ≤0.1, it means that the height change of the building is small and is temporarily within the safe range. The building will be included in the observation list for relevant units to conduct regular observations on the height, structure, etc. of the building, and promptly discover safety and stability issues of the building, provide important decision-making support for housing administrative departments, and better monitor and manage housing safety risk issues.

[0166] Step S03, when the building height change difference is greater than the height change threshold, the building is included in the danger list.

[0167] It should be noted that the danger list is a list of buildings in densely populated urban building complexes whose heights have changed and the range of change exceeds the safe range. The list includes the building's name, location information, height information, building type, pedestrian flow and other information.

[0168] When the height change difference of the building is not greater than the height change threshold, that is, Δ>0.1, it means that the height change of the building is large and the change range has exceeded the safe range. The building will be included in the dangerous list for relevant units to inspect the height, structure, etc. of the building and determine whether it is necessary to take measures such as demolition and warning to avoid safety accidents such as building collapse, structural damage, and material falling.

[0169] By judging the risk of buildings through the difference in building height changes, it is possible to monitor and manage housing safety risks.

[0170] This embodiment provides a method for identifying building height changes based on remote sensing spatiotemporal information, wherein a first relative projection displacement ratio and a second relative projection displacement ratio are obtained according to the relative projection displacement ratios of multiple buildings; a first ratio of the first relative projection displacement ratio and the second relative projection displacement ratio is calculated; when the first ratio meets a first standard and the absolute projection displacement ratio of the single building is equal to the satellite inclination tangent ratio, it is determined that the height of the building has not changed; when the first ratio does not meet the first standard or the absolute projection displacement ratio of the single building is not equal to the satellite inclination tangent ratio, it is determined that the height of the building has changed, thereby providing a monitoring basis for housing safety risk issues.

[0171] For example, to help understand the implementation process of the building height change identification method based on remote sensing spatiotemporal information of the first to fourth embodiments, please refer to Figure 7 , Figure 7 A brief flow chart of a building height change recognition method based on remote sensing spatiotemporal information is provided. Specifically, the top contour of the building is obtained through time-series high-resolution remote sensing images; the bottom contour of the building is obtained through vector maps; the top and bottom projection displacements of the building are calculated; when comparing multiple buildings, T i Period and T i+1 Relative projection displacement ratio of a group of buildings in the period and When comparing single buildings, obtain T i Period and T i+1 The absolute projection displacement ratio R of a single building during the period; and When they are approximately equal and R is equal to the tangent ratio of the two satellite inclination angles, it is determined that the building height has not changed; and When they are not approximately equal or R is not equal to the tangent ratio of the two satellite inclination angles, the change in building height is determined.

[0172] Based on the first embodiment of the present application, in the fifth embodiment of the present application, the same or similar contents as those in the first embodiment can be referred to the above description, and will not be described in detail later. Figure 8 , step S20 includes steps S21-S22:

[0173] Step S21, using multi-temporal high-resolution optical remote sensing images to obtain images of the top of the building at different times.

[0174] It should be noted that the building top image refers to the image obtained through remote sensing imaging technology, which shows the edge lines of the building top. This image can clearly depict the shape and structure of the building top, and has important application value in urban planning, geographic information system (GIS), disaster assessment and building monitoring. The building top outline image refers to the representation of the outline of the building top on a two-dimensional plane, which is usually extracted from aerial photographs, satellite images or other remote sensing images through image processing technology, and can be used in various application scenarios such as urban planning and building analysis.

[0175] Step S22, processing the building top image based on an image segmentation algorithm to obtain a building top outline.

[0176] It should be noted that the image segmentation algorithm is a technology used in the field of computer vision and image processing to divide an image into multiple regions or objects. The purpose is to organize the pixels in the image into meaningful parts in order to better understand and analyze the image content. There are various types of image segmentation algorithms, including threshold segmentation, edge detection algorithms, graph cut algorithms, clustering algorithms, deep learning methods, etc. The selection of image segmentation algorithms is not limited in this embodiment.

[0177] In addition, it should be noted that the building top outline refers to data in which the geometric shape of the building top boundary is represented in a vector form in a computer.

[0178] Based on the image segmentation algorithm, the top outline is extracted from the image of the building top. The process of extracting the top outline data is iterative and may require multiple adjustments and optimizations to achieve the best results. In actual applications, this process may be more complicated and involve advanced technologies such as multi-scale analysis and multi-modal data fusion.

[0179] This embodiment provides a method for identifying changes in building height based on remote sensing spatiotemporal information. It uses multi-phase high-resolution optical remote sensing images to obtain images of building tops at different times. The images of building tops are processed based on an image segmentation algorithm to obtain the outline of the building top, thereby providing a data basis for identifying changes in building height.

[0180] It should be noted that the above examples are only used to understand the present application and do not constitute a limitation on the method for identifying building height changes based on remote sensing spatiotemporal information of the present application. More simple transformations based on this technical concept are all within the scope of protection of the present application.

[0181] This application also provides a building height change recognition device based on remote sensing spatiotemporal information, please refer to Fig. 9 , the device comprises:

[0182] The data acquisition module 10 is used to acquire the bottom contour of the building using the vector map;

[0183] The data acquisition module 10 is also used to acquire the top contours of buildings at different times using multi-temporal high-resolution optical remote sensing images;

[0184] A data calculation module 20, for obtaining the projection displacement of the top and bottom contours of the building according to the Euclidean distance between the geometric center point of the bottom contour of the building and the geometric center point of the top contour of the building;

[0185] The data calculation module 20 is further used to calculate the relative projection displacement ratio of multiple buildings and the absolute projection displacement ratio of a single building of different buildings according to the projection displacement of the top and bottom contours of the buildings at different times;

[0186] The data calculation module 20 is further used to calculate the satellite inclination tangent ratio according to the satellite observation angles at different times;

[0187] The identification and analysis module 30 is used to identify the height change of the building according to the relative projection displacement ratio of multiple buildings, the absolute projection displacement ratio of a single building and the satellite inclination tangent ratio.

[0188] In one embodiment, the data calculation module 20 is further used to select a first building and a second building from the buildings; obtain a first projection displacement of the first building and a second projection displacement of the second building at a first moment, a third projection displacement of the first building and a fourth projection displacement of the second building at a second moment, a fifth projection displacement of the first building at a third moment, and a sixth projection displacement of the first building at a fourth moment according to the projection displacement of the top and bottom contours of the buildings; obtain a relative projection displacement ratio of multiple buildings according to the first projection displacement, the second projection displacement, the third projection displacement and the fourth projection displacement; obtain an absolute projection displacement ratio of a single building according to the fifth projection displacement and the sixth projection displacement.

[0189] In one embodiment, the data calculation module 20 is further used to obtain a first relative projection displacement ratio according to the first projection displacement and the second projection displacement, and to obtain a second relative projection displacement ratio according to the third projection displacement and the fourth projection displacement; and the first relative projection displacement ratio and the second relative projection displacement ratio are used as the relative projection displacement ratio of multiple buildings.

[0190] In one embodiment, the data calculation module 20 is also used to obtain a first observation inclination angle of the satellite relative to the first building at a third moment and a second observation inclination angle of the satellite relative to the first building at a fourth moment; and obtain the satellite inclination tangent ratio based on the first observation inclination angle and the second observation inclination angle.

[0191] In one embodiment, the identification and analysis module 30 is also used to obtain a first relative projection displacement ratio and a second relative projection displacement ratio based on the relative projection displacement ratios of the multiple buildings; calculate a first ratio of the first relative projection displacement ratio and the second relative projection displacement ratio; when the first ratio meets a first standard and the absolute projection displacement ratio of the single building is equal to the satellite inclination tangent ratio, determine that the height of the building has not changed; when the first ratio does not meet the first standard or the absolute projection displacement ratio of the single building is not equal to the satellite inclination tangent ratio, determine that the height of the building has changed.

[0192] In one embodiment, the identification and analysis module 30 is also used to obtain the building height change difference based on the relative projection displacement ratio of the multiple buildings and / or the absolute projection displacement ratio of the single building; when the building height change difference is not greater than the height change threshold, the building is included in the observation list; when the building height change difference is greater than the height change threshold, the building is included in the danger list.

[0193] In one embodiment, the data acquisition module 10 is further used to obtain images of the top of a building at different times using multi-temporal high-resolution optical remote sensing images; and to process the images of the top of the building based on an image segmentation algorithm to obtain the outline of the top of the building.

[0194] The building height change identification device based on remote sensing spatiotemporal information provided by the present application adopts the building height change identification method based on remote sensing spatiotemporal information in the above-mentioned embodiment, which can solve the technical problem of how to improve the accuracy of identifying the height change of a single building in a dense urban building complex. Compared with the prior art, the beneficial effects of the building height change identification device based on remote sensing spatiotemporal information provided by the present application are the same as the beneficial effects of the building height change identification method based on remote sensing spatiotemporal information provided by the above-mentioned embodiment, and the other technical features of the building height change identification device based on remote sensing spatiotemporal information are the same as the features disclosed in the above-mentioned embodiment method, which will not be repeated here.

[0195] The present application provides a device for identifying building height changes based on remote sensing spatiotemporal information. The device for identifying building height changes based on remote sensing spatiotemporal information includes: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor can execute the method for identifying building height changes based on remote sensing spatiotemporal information in the above-mentioned embodiment one.

[0196] Reference below Fig.10, which shows a schematic diagram of the structure of a building height change recognition device based on remote sensing spatiotemporal information suitable for implementing the embodiment of the present application. The building height change recognition device based on remote sensing spatiotemporal information in the embodiment of the present application may include but is not limited to mobile terminals such as mobile phones, laptop computers, digital broadcast receivers, PDAs (Personal Digital Assistants), PADs (Portable Application Descriptions), PMPs (Portable Media Players), vehicle-mounted terminals (such as vehicle-mounted navigation terminals), etc., and fixed terminals such as digital TVs, desktop computers, etc. Fig.10 The building height change identification device based on remote sensing spatiotemporal information shown is merely an example and should not bring any limitation to the functions and scope of use of the embodiments of the present application.

[0197] like Fig.10 As shown, the building height change recognition device based on remote sensing spatiotemporal information may include a processing device 1001 (e.g., a central processing unit, a graphics processing unit, etc.), which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM: Read Only Memory) 1002 or a program loaded from a storage device 1003 to a random access memory (RAM: Random Access Memory) 1004. Various programs and data required for the operation of the building height change recognition device based on remote sensing spatiotemporal information are also stored in RAM1004. The processing device 1001, ROM1002, and RAM1004 are connected to each other via a bus 1005. An input / output (I / O) interface 1006 is also connected to the bus. Typically, the following systems can be connected to the I / O interface 1006: input devices 1007 including, for example, a touch screen, a touchpad, a keyboard, a mouse, an image sensor, a microphone, an accelerometer, a gyroscope, etc.; output devices 1008 including, for example, a liquid crystal display (LCD), a speaker, a vibrator, etc.; storage devices 1003 including, for example, a magnetic tape, a hard disk, etc.; and communication devices 1009. The communication device 1009 can allow the building height change identification device based on remote sensing spatiotemporal information to communicate wirelessly or wired with other devices to exchange data. Although the figure shows a building height change identification device based on remote sensing spatiotemporal information with various systems, it should be understood that it is not required to implement or have all the systems shown. More or fewer systems may be implemented or have instead.

[0198] In particular, according to the embodiments disclosed in the present application, the process described above with reference to the flowchart can be implemented as a computer software program. For example, the embodiments disclosed in the present application include a computer program product, which includes a computer program carried on a computer-readable medium, and the computer program includes a program code for executing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from a network through a communication device, or installed from a storage device 1003, or installed from a ROM 1002. When the computer program is executed by the processing device 1001, the above-mentioned functions defined in the method of the embodiment disclosed in the present application are executed.

[0199] The building height change recognition device based on remote sensing spatiotemporal information provided by the present application adopts the building height change recognition method based on remote sensing spatiotemporal information in the above-mentioned embodiment, which can solve the technical problem of how to improve the accuracy of identifying the height change of a single building in a dense urban building complex. Compared with the prior art, the beneficial effects of the building height change recognition device based on remote sensing spatiotemporal information provided by the present application are the same as the beneficial effects of the building height change recognition method based on remote sensing spatiotemporal information provided by the above-mentioned embodiment, and the other technical features of the building height change recognition device based on remote sensing spatiotemporal information are the same as the features disclosed in the method of the previous embodiment, which will not be repeated here.

[0200] It should be understood that the various parts disclosed in this application can be implemented by hardware, software, firmware or a combination thereof. In the description of the above embodiments, specific features, structures, materials or characteristics can be combined in any one or more embodiments or examples in a suitable manner.

[0201] The above is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any person skilled in the art who is familiar with the present technical field can easily think of changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application. Therefore, the protection scope of the present application should be based on the protection scope of the claims.

[0202] The present application provides a computer-readable storage medium having computer-readable program instructions (ie, computer programs) stored thereon, and the computer-readable program instructions are used to execute the building height change identification method based on remote sensing spatiotemporal information in the above-mentioned embodiment.

[0203] The computer-readable storage medium provided in the present application may be, for example, a USB flash drive, but is not limited to electrical, magnetic, optical, electromagnetic, infrared, or semiconductor systems, systems or devices, or any combination of the above. More specific examples of computer-readable storage media may include, but are not limited to: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In this embodiment, the computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in combination with an instruction execution system, system or device. The program code contained on the computer-readable storage medium may be transmitted using any appropriate medium, including but not limited to: wires, optical cables, RF (Radio Frequency), etc., or any suitable combination of the above.

[0204] The computer-readable storage medium may be included in the device for identifying building height changes based on remote sensing spatiotemporal information; or may exist independently without being assembled into the device for identifying building height changes based on remote sensing spatiotemporal information.

[0205] The computer-readable storage medium carries one or more programs. When the one or more programs are executed by a building height change identification device based on remote sensing spatiotemporal information, the building height change identification device based on remote sensing spatiotemporal information: obtains the bottom contour of the building using a vector map; obtains the top contour of the building at different times using a multi-phase high-resolution optical remote sensing image; obtains the projection displacement of the top and bottom contours of the building based on the Euclidean distance between the geometric center point of the bottom contour of the building and the geometric center point of the top contour of the building; calculates the relative projection displacement ratio of multiple buildings and the absolute projection displacement ratio of a single building of different buildings based on the projection displacement of the top and bottom contours of the building at different times; calculates the satellite inclination tangent ratio based on the satellite observation angles at the different times; and identifies the height change of the building based on the relative projection displacement ratio of multiple buildings, the absolute projection displacement ratio of a single building and the satellite inclination tangent ratio.

[0206] Computer program code for performing the operations of the present application may be written in one or more programming languages ​​or a combination thereof, including object-oriented programming languages ​​such as Java, Smalltalk, C++, and conventional procedural programming languages ​​such as "C" or similar programming languages. The program code may be executed entirely on the user's computer, partially on the user's computer, as a separate software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0207] The flow chart and block diagram in the accompanying drawings illustrate the possible architecture, function and operation of the system, method and computer program product according to various embodiments of the present application. In this regard, each square box in the flow chart or block diagram can represent a module, a program segment or a part of a code, and the module, the program segment or a part of the code contains one or more executable instructions for realizing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the square box can also occur in a sequence different from that marked in the accompanying drawings. For example, two square boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each square box in the block diagram and / or flow chart, and the combination of the square boxes in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs a specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0208] The modules involved in the embodiments described in this application may be implemented by software or hardware, wherein the name of the module does not constitute a limitation on the unit itself in some cases.

[0209] The readable storage medium provided by the present application is a computer-readable storage medium, which stores computer-readable program instructions (i.e., computer programs) for executing the above-mentioned building height change identification method based on remote sensing spatiotemporal information, and can solve the technical problem of how to improve the accuracy of identifying the height change of a single building in a dense urban building complex. Compared with the prior art, the beneficial effects of the computer-readable storage medium provided by the present application are the same as the beneficial effects of the building height change identification method based on remote sensing spatiotemporal information provided by the above-mentioned embodiment, and will not be repeated here.

[0210] The present application also provides a computer program product, including a computer program, which, when executed by a processor, implements the steps of the above-mentioned method for identifying building height changes based on remote sensing spatiotemporal information.

[0211] The computer program product provided by the present application can solve the technical problem of how to improve the accuracy of identifying the height changes of individual buildings in dense urban building complexes. Compared with the prior art, the beneficial effects of the computer program product provided by the present application are the same as the beneficial effects of the building height change identification method based on remote sensing spatiotemporal information provided by the above embodiment, and will not be repeated here.

[0212] The above descriptions are only some embodiments of the present application, and are not intended to limit the patent scope of the present application. All equivalent structural changes made using the contents of the present application specification and drawings under the technical concept of the present application, or direct / indirect applications in other related technical fields are included in the patent protection scope of the present application.

Claims

1. A method for identifying building height changes based on remote sensing spatiotemporal information, characterized in that: The method comprises: Use vector maps to obtain the bottom outline of buildings; Use multi-temporal high-resolution optical remote sensing images to obtain the top contours of buildings at different times; Obtaining the projection displacement of the top and bottom contours of the building according to the Euclidean distance between the geometric center point of the bottom contour of the building and the geometric center point of the top contour of the building; Calculate the relative projection displacement ratio of multiple buildings and the absolute projection displacement ratio of a single building of different buildings according to the projection displacement of the top and bottom contours of the buildings at different times; Calculating the satellite inclination tangent ratio according to the satellite observation angles at different times; The height change of the building is identified according to the relative projection displacement ratio of multiple buildings, the absolute projection displacement ratio of the single building and the satellite inclination tangent ratio.

2. The method according to claim 1, characterized in that The step of calculating the relative projection displacement ratio of multiple buildings and the absolute projection displacement ratio of a single building of different buildings according to the projection displacement of the top and bottom contours of the buildings at different times comprises: Selecting a first building and a second building from the buildings; According to the projection displacement of the top and bottom contours of the building, a first projection displacement of the first building and a second projection displacement of the second building at a first moment, a third projection displacement of the first building and a fourth projection displacement of the second building at a second moment, a fifth projection displacement of the first building at a third moment, and a sixth projection displacement of the first building at a fourth moment are obtained; Obtaining a relative projection displacement ratio of multiple buildings according to the first projection displacement, the second projection displacement, the third projection displacement and the fourth projection displacement; According to the fifth projection displacement and the sixth projection displacement, the absolute projection displacement ratio of a single building is obtained.

3. The method according to claim 2, characterized in that The step of obtaining the relative projection displacement ratio of multiple buildings according to the first projection displacement, the second projection displacement, the third projection displacement and the fourth projection displacement comprises: Obtain a first relative projection displacement ratio according to the first projection displacement and the second projection displacement, and obtain a second relative projection displacement ratio according to the third projection displacement and the fourth projection displacement; The first relative projection displacement ratio and the second relative projection displacement ratio are used as the multi-building relative projection displacement ratio.

4. The method according to claim 1, characterized in that The step of calculating the satellite inclination tangent ratio according to the satellite observation angles at different times comprises: Acquire a first observed inclination angle of a satellite relative to a first building at a third moment, and a second observed inclination angle of the satellite relative to the first building at a fourth moment; The satellite inclination tangent ratio is obtained according to the first observed inclination angle and the second observed inclination angle.

5. The method according to claim 1, characterized in that The step of identifying the height change of the building according to the relative projection displacement ratio of multiple buildings, the absolute projection displacement ratio of the single building and the satellite inclination tangent ratio comprises: Obtaining a first relative projection displacement ratio and a second relative projection displacement ratio according to the relative projection displacement ratios of the multiple buildings; Calculate a first ratio of the first relative projection displacement ratio and the second relative projection displacement ratio; When the first ratio satisfies a first criterion and the single building absolute projection displacement ratio is equal to the satellite inclination tangent ratio, determining that the height of the building has not changed; When the first ratio does not meet the first standard or the single building absolute projection displacement ratio is not equal to the satellite inclination tangent ratio, it is determined that the height of the building has changed.

6. The method according to claim 5, characterized in that The step of determining that the height of the building has changed when the first ratio does not meet the first standard or the single building absolute projection displacement ratio is not equal to the satellite inclination tangent ratio comprises: Obtaining the building height change difference according to the relative projection displacement ratio of the multiple buildings and / or the absolute projection displacement ratio of the single building; When the height change difference of the building is not greater than the height change threshold, the building is included in the observation list; When the building height change difference is greater than the height change threshold, the building is included in the danger list.

7. The method according to any one of claims 1 to 6, characterized in that The step of obtaining the top contours of buildings at different times using multi-temporal high-resolution optical remote sensing images comprises: Using multi-temporal high-resolution optical remote sensing images, we can obtain images of building tops at different times; The image of the building top is processed based on an image segmentation algorithm to obtain a building top outline.

8. A device for identifying building height changes based on remote sensing spatiotemporal information, characterized in that: The device comprises: A data acquisition module, used for acquiring the bottom contour of the building using a vector map; The data acquisition module is also used to acquire the top contours of buildings at different times using multi-temporal high-resolution optical remote sensing images; A data calculation module, used to obtain the projection displacement of the top and bottom contours of the building according to the Euclidean distance between the geometric center point of the bottom contour of the building and the geometric center point of the top contour of the building; The data calculation module is further used to calculate the relative projection displacement ratio of multiple buildings and the absolute projection displacement ratio of a single building of different buildings according to the projection displacement of the top and bottom contours of the buildings at different times; The data calculation module is further used to calculate the satellite inclination tangent ratio according to the satellite observation angles at different times; The identification and analysis module is used to identify the height change of the building according to the relative projection displacement ratio of multiple buildings, the absolute projection displacement ratio of the single building and the satellite inclination tangent ratio.

9. A device for identifying building height changes based on remote sensing spatiotemporal information, characterized in that: The device comprises: a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the computer program is configured to implement the steps of the method for identifying building height changes based on remote sensing spatiotemporal information as described in any one of claims 1 to 7.

10. A storage medium, characterized in that: The storage medium is a computer-readable storage medium, and a computer program is stored on the storage medium. When the computer program is executed by a processor, the steps of the method for identifying building height changes based on remote sensing spatiotemporal information as described in any one of claims 1 to 7 are implemented.

Images