A method and device for determining image coverage area of ​​remote sensing image

Abstract

The present application provides a method and device for determining the image coverage area of ​​a remote sensing image, including: obtaining associated remote sensing images that are spatially associated with a target statistical area; for each associated remote sensing image, obtaining coverage status information corresponding to the associated remote sensing image in the associated area; storing the coverage status information corresponding to the associated remote sensing image in a blank bitmap corresponding to the associated remote sensing image in a preset storage order; obtaining the bitmap corresponding to each associated remote sensing image that is spatially associated with the target statistical area in the predetermined time period; based on the coverage status information on the bitmap corresponding to each associated remote sensing image, determining the real image coverage area of ​​the remote sensing image of the target statistical area in the predetermined time period. Through the method and device, the real image coverage area of ​​the remote sensing image of the target statistical area in any predetermined time period can be calculated immediately when a statistical request is obtained.

Description

Technical Field

[0001] The present application relates to the field of remote sensing image technology, and in particular to a method and device for determining an image coverage area of ​​a remote sensing image. Background Art

[0002] In recent years, with the increasing number of domestic remote sensing satellites, the coverage of remote sensing images has become stronger and stronger. Due to certain business needs, it is usually necessary to count the image coverage area of ​​the remote sensing images of the target statistical area within any predetermined time period, so as to evaluate the image coverage ability of the satellite for the target statistical area within any predetermined time period. However, due to the differences in parameters such as the revisit period and width of different satellites, different satellites can capture different areas within the target statistical area within the predetermined time period for the target statistical area within the predetermined time period; in addition, due to the influence of factors such as cloud and rain weather, sensors and data transmission, the same satellite may collect some areas of the single-phase remote sensing images for the target statistical area within the predetermined time period without being covered, and there will be multiple overlapping areas in the collected multi-phase remote sensing images. Based on the complex situation of the remote sensing images collected by the above satellites, it becomes very difficult to determine the true image coverage area of ​​the remote sensing images of the target statistical area within any predetermined time period.

[0003] Currently, in the known solutions for determining the image coverage area of ​​remote sensing images of a target statistical area within any predetermined time period, when a large number of remote sensing images of the target statistical area within the predetermined time period are collected by the satellite, it will take a lot of time to do an area statistics. Summary of the invention

[0004] In view of this, the purpose of the present application is to provide a method and a device for determining the image coverage area of ​​a remote sensing image, which can overcome the defect in the existing scheme that when determining the image coverage area of ​​a remote sensing image of a target statistical area within any predetermined time period, when there are many remote sensing images of the target statistical area collected by the satellite within the predetermined time period, it will take a lot of time to make an area statistics.

[0005] In a first aspect, an embodiment of the present application provides a method for determining an image coverage area of ​​a remote sensing image, the determination method comprising:

[0006] Acquire an associated remote sensing image that is spatially associated with a target statistical area; wherein the target statistical area corresponds to an associated area, the target statistical area is smaller than the associated area, the associated area is pre-divided into a predetermined number of statistical grids, each associated remote sensing image corresponds to a blank bitmap for storing information, and the number of bits of each blank bitmap is the same as the number of grids in the associated area;

[0007] For each associated remote sensing image, preprocessing the associated remote sensing image is performed to obtain coverage status information corresponding to the associated remote sensing image in the associated area;

[0008] The coverage status information corresponding to the associated remote sensing image is stored in a blank bitmap corresponding to the associated remote sensing image according to a preset storage order;

[0009] In response to a request for counting the image coverage area of ​​a target statistical area in a predetermined time period, obtaining a bitmap corresponding to each associated remote sensing image that is spatially associated with the target statistical area in the predetermined time period;

[0010] Based on the coverage status information on the bitmap corresponding to each associated remote sensing image, the actual image coverage area of ​​the remote sensing image of the target statistical area in the predetermined time period is determined.

[0011] Optionally, for each associated remote sensing image, preprocessing the associated remote sensing image to obtain coverage status information corresponding to the associated remote sensing image in the associated area includes:

[0012] For each associated remote sensing image, determine the image range block of the overlapping area between the associated remote sensing image and the target statistical area;

[0013] Overlay the image range block to the corresponding position on the statistical grid of the associated area to obtain an overlay grid;

[0014] Based on the center line of each row of grids in the superimposed grid, coverage status information of the associated area corresponding to the associated remote sensing image is obtained.

[0015] Optionally, obtaining coverage status information of the associated area corresponding to the associated remote sensing image based on the center line of each row of grids in the superimposed grid includes:

[0016] For each row of grids in the superimposed grid, pre-covered boundary grids are selected from the row of grids, where the boundary of the image range block intersects with the center line of the row of grids;

[0017] Based on the position serial number of each pre-covered boundary grid, all grids except the pre-covered boundary grid in the row of grids are screened to obtain covered grids among all grids except the pre-covered boundary grid; the position serial number indicates the order of each pre-covered boundary grid in the pre-covered boundary grids in the row of grids;

[0018] Based on the position sequence number of each pre-covered boundary grid and the center point of each pre-covered boundary grid, all pre-covered boundary grids are screened to obtain the covering grids in the pre-covered boundary grids;

[0019] Determine the other grids except the covered grids screened out from the row of grids as uncovered grids;

[0020] The first coverage status information is determined as the coverage status information of the covered grid, and the second coverage status information is determined as the coverage status information of the uncovered grid.

[0021] Optionally, the step of screening all grids in the row of grids except the pre-covered boundary grid based on the position serial number of each pre-covered boundary grid to obtain the covering grids in all grids except the pre-covered boundary grid includes:

[0022] According to the order indicated by the position numbers of the pre-covered boundary grids, the pre-covered boundary grids corresponding to the odd position numbers and the pre-covered boundary grids corresponding to the even position numbers are searched from left to right in sequence;

[0023] For each pre-covered boundary grid corresponding to an odd-numbered position serial number, all grids between the pre-covered boundary grid and the pre-covered boundary grid corresponding to an even-numbered position serial number adjacent to the pre-covered boundary grid to the right are determined as covered grids.

[0024] Optionally, the step of screening all pre-covered boundary grids based on the position serial number of each pre-covered boundary grid and the center point of each pre-covered boundary grid to obtain the covering grid in the pre-covered boundary grid includes:

[0025] For each pre-covered boundary grid in the row of grids, determining the parity of the position sequence number of the pre-covered boundary grid;

[0026] When the position sequence number of the pre-covering boundary grid is an odd number, determining whether the intersection point of the pre-covering boundary grid is on the left side of the center point of the pre-covering boundary grid;

[0027] If the intersection point of the pre-covering boundary grid is on the left side of the center point of the pre-covering boundary grid, the pre-covering boundary grid is determined as a covering grid;

[0028] When the position sequence number of the pre-covering boundary grid is an even number, determining whether the intersection point of the pre-covering boundary grid is on the right side of the center point of the pre-covering boundary grid;

[0029] If the intersection point of the pre-covering boundary grid is on the right side of the center point of the pre-covering boundary grid, the pre-covering boundary grid is determined as a covering grid.

[0030] Optionally, determining the real image coverage area of ​​the remote sensing image of the target statistical area in the predetermined time period based on the coverage status information on the bitmap corresponding to each associated remote sensing image includes:

[0031] According to the arrangement order of the position indexes corresponding to the bitmaps, the coverage status information of each bit on each bitmap corresponding to each position index is obtained in sequence;

[0032] For each position index, an OR operation is performed on the coverage status information of the corresponding bit on each bitmap corresponding to the position index to obtain the merged coverage status information corresponding to the position index; the merged coverage status information includes merged coverage and merged non-coverage;

[0033] Filter out the merged coverage status information of the merged coverage from all the merged coverage status information;

[0034] Based on the combined coverage state information of the combined coverage, the actual image coverage area of ​​the remote sensing image of the target statistical area within the predetermined time period is determined.

[0035] Optionally, the step of determining the real image coverage area of ​​the remote sensing image of the target statistical area within the predetermined time period based on the combined coverage state information of the combined coverage includes:

[0036] Determine a target position index corresponding to the merged coverage state information of the merged coverage;

[0037] Based on the correspondence between the preset position index and the statistical grid, a target statistical grid corresponding to the target position index is obtained;

[0038] Based on the correspondence between the preset statistical grids and the real geographical areas, the real geographical area corresponding to each target statistical grid is obtained; wherein the real geographical area corresponding to each statistical grid represents the area of ​​the real geographical region corresponding to the statistical grid;

[0039] The sum of the real geographic areas corresponding to each target statistical grid is determined as the real image coverage area of ​​the remote sensing image of the target statistical area in the predetermined time period.

[0040] In a second aspect, an embodiment of the present application provides a device for determining an image coverage area of ​​a remote sensing image, the device comprising:

[0041] The first acquisition module is used to acquire an associated remote sensing image that is spatially associated with the target statistical area; wherein the target statistical area corresponds to an associated area, the target statistical area is smaller than the associated area, the associated area is pre-divided into a predetermined number of statistical grids, each associated remote sensing image corresponds to a blank bitmap for storing information, and the number of bits of each blank bitmap is the same as the number of grids in the associated area;

[0042] A preprocessing module is used to preprocess each associated remote sensing image to obtain coverage status information of the associated remote sensing image in the associated area;

[0043] A storage module, used for storing the coverage status information corresponding to the associated remote sensing image in a blank bitmap corresponding to the associated remote sensing image according to a preset storage order;

[0044] A second acquisition module is used to respond to a request for counting the image coverage area of ​​a target statistical area in a predetermined time period, and to acquire a bitmap corresponding to each associated remote sensing image that is spatially associated with the target statistical area in the predetermined time period;

[0045] The statistical module is used to determine the real image coverage area of ​​the remote sensing image of the target statistical area within the predetermined time period based on the coverage status information on the bitmap corresponding to each associated remote sensing image.

[0046] In a third aspect, an embodiment of the present application provides an electronic device, comprising: a processor, a memory and a bus, wherein the memory stores machine-readable instructions executable by the processor, and when the electronic device is running, the processor and the memory communicate via the bus, and when the machine-readable instructions are executed by the processor, the steps of the method for determining the image coverage area of ​​the remote sensing image as described above are performed.

[0047] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the steps of the method for determining the image coverage area of ​​the remote sensing image as described above are executed.

[0048] The method and device for determining the image coverage area of ​​a remote sensing image provided in the embodiments of the present application can store the coverage status information of the associated remote sensing image in the associated area in the form of a bitmap when receiving an associated remote sensing image that is spatially associated with a target statistical area, so that when counting the actual image coverage area of ​​the remote sensing image of the target statistical area within any predetermined time period, area statistics can be directly performed based on the bitmap, so that when counting the actual image coverage area of ​​the remote sensing image of the target statistical area within any predetermined time period, the actual image coverage area of ​​the remote sensing image of the target statistical area within the predetermined time period can be calculated immediately when a statistical request is obtained.

[0049] In order to make the above-mentioned objects, features and advantages of the present application more obvious and easy to understand, preferred embodiments are specifically cited below and described in detail with reference to the attached drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the embodiments will be briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present application and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

[0051] Figure 1 A flow chart showing a method for determining an image coverage area of ​​a remote sensing image provided by an exemplary embodiment of the present application is shown;

[0052] Figure 2 A schematic diagram showing a superimposed grid and a corresponding bitmap obtained according to an exemplary embodiment of the present application is shown;

[0053] Figure 3 A schematic diagram showing a bitmap corresponding to a plurality of associated remote sensing images within a predetermined time period provided by an exemplary embodiment of the present application;

[0054] Figure 4 A schematic structural diagram of a device for determining image coverage area of ​​a remote sensing image provided by an exemplary embodiment of the present application is shown;

[0055] Figure 5 A schematic structural diagram of an electronic device provided by an exemplary embodiment of the present application is shown. DETAILED DESCRIPTION

[0056] To make the purpose, technical scheme and advantages of the embodiments of the present application clearer, the technical scheme in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all of the embodiments. The components of the embodiments of the present application usually described and shown in the drawings here can be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present application provided in the drawings is not intended to limit the scope of the application claimed for protection, but merely represents the selected embodiments of the present application. Based on the embodiments of the present application, each other embodiment obtained by those skilled in the art without making creative work belongs to the scope of protection of the present application.

[0057] Currently, in the known solutions for determining the image coverage area of ​​remote sensing images of a target statistical area within any predetermined time period, when a large number of remote sensing images of the target statistical area within the predetermined time period are collected by the satellite, it will take a lot of time to do an area statistics.

[0058] Based on this, the embodiments of the present application provide a method and device for determining the image coverage area of ​​a remote sensing image, which can overcome the defect of the existing scheme that it takes a long time to determine the image coverage area of ​​a remote sensing image of a target statistical area within any predetermined time period.

[0059] See also Figure 1 , Figure 1 A flow chart of a method for determining an image coverage area of ​​a remote sensing image provided by an exemplary embodiment of the present application is shown.

[0060] like Figure 1 As shown in , an exemplary embodiment of the present application provides a method for determining the image coverage area of ​​a remote sensing image, comprising the following steps:

[0061] S1. Acquire an associated remote sensing image that is spatially associated with a target statistical area; wherein the target statistical area corresponds to an associated area, the target statistical area is smaller than the associated area, the associated area is pre-divided into a predetermined number of statistical grids, each associated remote sensing image corresponds to a blank bitmap for storing information, and the number of bits of each blank bitmap is the same as the number of statistical grids of the associated area;

[0062] Here, the target statistical area is an area of ​​expected statistical area, for example, the target statistical area may be the geographical area where China is located.

[0063] Here, the acquired remote sensing image spatially associated with the target statistical area refers to a remote sensing image in which a shooting area indicated by the spatial information of the remote sensing image overlaps with the target statistical area.

[0064] Here, the associated area is an area related to subsequent storage. For example, when the target statistical area is the geographical area where China is located, the associated area is the area corresponding to the outer rectangular frame containing the geographical area where China is located.

[0065] The number of statistical grids into which the associated area is divided may be determined based on statistical accuracy. When the statistical accuracy is high, the number of statistical grids into which the associated area is divided is large, and when the statistical accuracy is low, the number of statistical grids into which the associated area is divided is small.

[0066] A bitmap is a storage format for storing information. Each associated remote sensing image corresponds to a blank bitmap for storing information of the associated remote sensing image.

[0067] S2. For each associated remote sensing image, preprocess the associated remote sensing image to obtain coverage status information corresponding to the associated remote sensing image in the associated area;

[0068] As an example, the steps may include:

[0069] S21, for each associated remote sensing image, determining an image range block of an overlapping area between the associated remote sensing image and the target statistical area;

[0070] As an example, in this step, firstly, for each associated remote sensing image, an image coverage area where the associated remote sensing image overlaps with the target statistical area can be determined;

[0071] Then, based on the image coverage area, an image range block having the same size as the image coverage area is obtained;

[0072] Here, the image range data of the image coverage area can be determined first, then the area surrounded by the image range data can be determined as the area block corresponding to the image range data, and finally the area block corresponding to the image range data can be determined as the image range block having the same size as the image coverage area.

[0073] Here, the image range data indicates the position of the boundary of the image coverage area. As an example, the image range data may be coordinate data.

[0074] Here, in the area block corresponding to the image range data, there is only the position information of the most edge of the area block, and no information of the middle part of the area block. In this way, the operation speed in subsequent calculation and storage can be faster.

[0075] S22, superimposing the image range block to a corresponding position on the statistical grid of the associated area to obtain a superimposed grid;

[0076] Here, the corresponding position refers to the position of the image range block in the target statistical area.

[0077] S23. Based on the center line of each row of grids in the superimposed grid, obtain coverage status information of the associated area corresponding to the associated remote sensing image.

[0078] As an example, the steps may include:

[0079] S231, for each row of grids in the superimposed grid, filter out pre-covered boundary grids from the row of grids, where the boundary of the image range block intersects with the center line of the row of grids;

[0080] S232, based on the position serial number of each pre-covered boundary grid, all grids except the pre-covered boundary grid in the row of grids are screened to obtain covered grids among all grids except the pre-covered boundary grid; the position serial number indicates the order of each pre-covered boundary grid in the pre-covered boundary grids in the row of grids;

[0081] As an example, in this step, the pre-covered boundary grids corresponding to odd-numbered position numbers and the pre-covered boundary grids corresponding to even-numbered position numbers may be searched from left to right in sequence according to the order indicated by the position numbers of the pre-covered boundary grids;

[0082] Then, for each pre-covered boundary grid corresponding to an odd-numbered position serial number, all grids between the pre-covered boundary grid and the pre-covered boundary grid corresponding to an even-numbered position serial number adjacent to the pre-covered boundary grid to the right are determined as covered grids.

[0083] S233, based on the position serial number of each pre-covered boundary grid and the center point of each pre-covered boundary grid, all pre-covered boundary grids are screened to obtain a covering grid in the pre-covered boundary grids;

[0084] For each pre-covered boundary grid in the row of grids, determining the parity of the position sequence number of the pre-covered boundary grid;

[0085] When the position sequence number of the pre-covering boundary grid is an odd number, determining whether the intersection point of the pre-covering boundary grid is on the left side of the center point of the pre-covering boundary grid;

[0086] If the intersection point of the pre-covering boundary grid is on the left side of the center point of the pre-covering boundary grid, the pre-covering boundary grid is determined as a covering grid;

[0087] When the position sequence number of the pre-covering boundary grid is an even number, determining whether the intersection point of the pre-covering boundary grid is on the right side of the center point of the pre-covering boundary grid;

[0088] If the intersection point of the pre-covering boundary grid is on the right side of the center point of the pre-covering boundary grid, the pre-covering boundary grid is determined as a covering grid.

[0089] S234, determining the other grids except the covered grids screened out from the row of grids as uncovered grids;

[0090] S235: Determine the first coverage status information as the coverage status information of the covered grid, and determine the second coverage status information as the coverage status information of the uncovered grid.

[0091] As an example, the first coverage status information may be "1" and the second coverage status information may be "0".

[0092] S3, storing the coverage status information corresponding to the associated remote sensing image into the corresponding blank bitmap according to a preset storage order;

[0093] Here, the number of bits in each blank bitmap is the same as the number of grids in the associated area, so it can be understood that each grid in the associated area corresponds to one bit in the blank bitmap. For example, when the number of grids in the associated area is (m rows, n columns), the number of bits in each blank bitmap may be (1 row, m×n columns).

[0094] Therefore, in this step, the coverage status information of the remote sensing image (m rows, n columns) corresponding to the remote sensing image can be written into the bits in the corresponding blank bitmap in sequence according to the storage order from top to bottom and from left to right to obtain the bitmap corresponding to the associated remote sensing image.

[0095] As an example, after obtaining the bitmap corresponding to the associated remote sensing image, the bitmap can be stored in Redis or memory.

[0096] Next, we will combine Figure 2 In a specific example, it is introduced how to obtain the coverage status information and the corresponding bitmap of the associated area corresponding to the associated remote sensing image based on the center line of each row of grids in the superimposed grid.

[0097] See also Figure 2 , Figure 2 A schematic diagram of a superimposed grid and a corresponding bitmap obtained according to an exemplary embodiment of the present application is shown, wherein m represents the number of grids in each column of the superimposed grid, and n represents the number of grids in each row of the superimposed grid.

[0098] like Figure 3As shown, taking the second row of grids as an example, the pre-covered boundary grids "(3,4)", "(5,4)", "(9,4)", "(13,4)", "(19,4)", and "(22,4)" whose boundaries of the image range block and the center line "a" of the row of grids have intersection points "A", "B", "C", "D", "E", and "F" are screened out from the row of grids; then, according to the order indicated by the position numbers "1", "2", "3", "4", "5", and "6" of the pre-covered boundary grids, the pre-covered boundary grids "(3,4)", "(9,4)", and "19,4)" corresponding to the odd position numbers and the pre-covered boundary grids "(5,4)", "(13,4)", and "22,4)" corresponding to the even position numbers are searched from left to right in sequence; then, for the pre-covered boundary grid "(3,4)" corresponding to the odd position number "1", the pre-covered boundary grid "(3,4)" is removed. The grid (4,4) (not shown in the figure) between the pre-covered boundary grid (5,4) corresponding to the even position number "2" adjacent to the pre-covered boundary grid to the right is determined as the covering grid; for the pre-covered boundary grid "(9,4)" corresponding to the odd position number "3", the grids (10,4), (11,4), and (12,4) (not shown in the figure) between the pre-covered boundary grid and the pre-covered boundary grid (13,4) corresponding to the even position number "4" adjacent to the pre-covered boundary grid to the right are determined as the covering grids; for the pre-covered boundary grid "(19,4)" corresponding to the odd position number "5", the grids (20,4), (21,4) (not shown in the figure) between the pre-covered boundary grid and the pre-covered boundary grid (22,4) corresponding to the even position number "6" adjacent to the pre-covered boundary grid to the right are determined as the covering grids.Then, for the pre-covered boundary grid (3, 4) in the row of grids, first determine that the position serial number "1" of the pre-covered boundary grid is an odd number, and then determine whether the intersection point A of the pre-covered boundary grid is on the left side of the center point of the pre-covered boundary grid; because the intersection point of the pre-covered boundary grid is not on the left side of the center point of the pre-covered boundary grid, the pre-covered boundary grid is not a covered grid; for the pre-covered boundary grid (5, 4) in the row of grids, first determine that the position serial number "2" of the pre-covered boundary grid is an even number, and then determine whether the intersection point B of the pre-covered boundary grid is on the right side of the center point of the pre-covered boundary grid; because the intersection point of the pre-covered boundary grid is not on the right side of the center point of the pre-covered boundary grid, the pre-covered boundary grid is not a covered grid; For the pre-covered boundary grid (9,4) in the row of grids, first determine that the position number "3" of the pre-covered boundary grid is an odd number, and then determine whether the intersection point C of the pre-covered boundary grid is on the left side of the center point of the pre-covered boundary grid; because the intersection point of the pre-covered boundary grid is on the left side of the center point of the pre-covered boundary grid, the pre-covered boundary grid is determined to be a covered grid; for the pre-covered boundary grid (13,4) in the row of grids, first determine that the position number "4" of the pre-covered boundary grid is an odd number, and then determine whether the intersection point D of the pre-covered boundary grid is on the right side of the center point of the pre-covered boundary grid; because the intersection point of the pre-covered boundary grid is on the right side of the center point of the pre-covered boundary grid, the pre-covered boundary grid is determined to be a covered grid. Similarly, it can be obtained that the pre-covered boundary grid (19,4) is a covered grid, and the pre-covered boundary grid (22,4) is not a covered grid; then, the grids in the row of grids other than the covered grids that are screened out are determined to be uncovered grids. Finally, the first coverage status information "1" is determined as the coverage status information of the covered grid, and the second coverage status information "0" is determined as the coverage status information of the uncovered grid.

[0099] After obtaining the coverage status information of the associated area corresponding to the associated remote sensing image according to the above method, the coverage status information of (m rows, n columns) corresponding to the remote sensing image is written into the bits (0, m×n-1) in the blank bitmap corresponding to the associated remote sensing image in order from top to bottom and from left to right to obtain the bitmap corresponding to the associated remote sensing image.

[0100] According to the above method, when a remote sensing image spatially associated with a target statistical area is received, the coverage status information of the associated remote sensing image in the associated area can be stored in the form of a bitmap, so that when counting the actual image coverage area of ​​the remote sensing image of the target statistical area within any predetermined time period, the area statistics can be directly performed based on the bitmap, so that when counting the actual image coverage area of ​​the remote sensing image of the target statistical area within any predetermined time period, the actual image coverage area of ​​the remote sensing image of the target statistical area within any predetermined time period can be calculated immediately when the statistical request is obtained.

[0101] S4, in response to a request to count the image coverage area of ​​the target statistical area in a predetermined time period, obtaining a bitmap corresponding to each associated remote sensing image that is spatially associated with the target statistical area in the predetermined time period;

[0102] Here, as an example, the predetermined time period may be a time period before the moment of obtaining the statistical request, for example, the predetermined time period may be from 2020 to 2021;

[0103] S5. Determine the actual image coverage area of ​​the remote sensing image of the target statistical area in the predetermined time period based on the coverage status information on the bitmap corresponding to each associated remote sensing image.

[0104] As an example, the steps may include:

[0105] S51, according to the arrangement order of the position indexes corresponding to the bitmaps corresponding to each associated remote sensing image, sequentially obtain the coverage status information of each bit on each bitmap corresponding to each position index;

[0106] Here, the coverage status information of each bit includes coverage and non-coverage. For example, the coverage status information of coverage can be represented by "1", and the coverage status information of non-coverage can be represented by "0".

[0107] S52, for each position index, performing an OR operation on the coverage status information of the corresponding bit on each bitmap corresponding to the position index to obtain the merged coverage status information corresponding to the position index; the merged coverage status information includes merged coverage and merged non-coverage;

[0108] Here, the merged coverage status information of the merged coverage can be represented by "1", and the merged coverage status information of the merged non-coverage can be represented by "0".

[0109] See also Figure 3 , Figure 3A schematic diagram of bitmaps corresponding to a plurality of associated remote sensing images within a predetermined time period provided by an exemplary embodiment of the present application is shown, wherein k represents the number of associated remote sensing images within the predetermined time period, and the bitmap corresponding to k represents the bitmap corresponding to the kth associated remote sensing image.

[0110] Here, the OR operation is an existing operation method. For example, when the coverage status information of the corresponding bits on each bitmap corresponding to the first position index is "0", "0", "0", "0", "0", "0", respectively, the merged coverage status information corresponding to the position index obtained through the OR operation is "0", and when the coverage status information of the corresponding bits on each bitmap corresponding to the third position index is "0", "1", "1", "0", "1", respectively, the merged coverage status information corresponding to the position index obtained through the OR operation is "1".

[0111] S53, filtering out the merged coverage status information of the merged coverage from all the merged coverage status information;

[0112] S54: Determine the actual image coverage area of ​​the remote sensing image of the target statistical area within the predetermined time period based on the combined coverage state information of the combined coverage.

[0113] As an example, the steps may include:

[0114] S541, determining a target position index corresponding to the merged coverage state information of the merged coverage;

[0115] S542, based on the correspondence between the preset position index and the statistical grid, obtaining the target statistical grid corresponding to the target position index;

[0116] Here, since each statistical grid in the statistical grid of the associated area corresponds one-to-one to a bit in the bitmap, each position index actually corresponds to one of the grids in the statistical grid.

[0117] S543, based on the correspondence between the preset statistical grids and the real geographical areas, obtaining the real geographical area corresponding to each target statistical grid; the real geographical area corresponding to each statistical grid represents the area of ​​the real geographical region corresponding to the statistical grid;

[0118] Here, since the statistical grids are divided based on real geographical areas, each statistical grid actually corresponds to a real geographical area, and the real geographical area corresponds to a real geographical area, that is, each statistical grid corresponds to a real geographical area.

[0119] S544: Determine the sum of the real geographic areas corresponding to each target statistical grid as the real image coverage area of ​​the remote sensing image of the target statistical area in the predetermined time period.

[0120] A method for determining the image coverage area of ​​a remote sensing image provided in an embodiment of the present application can store the coverage status information of the associated remote sensing image in the associated area in the form of a bitmap when receiving a remote sensing image that is spatially associated with a target statistical area, so that when counting the actual image coverage area of ​​the remote sensing image of the target statistical area within any predetermined time period, area statistics can be directly performed based on the bitmap, so that when counting the actual image coverage area of ​​the remote sensing image of the target statistical area within any predetermined time period, the actual image coverage area of ​​the remote sensing image of the target statistical area within any predetermined time period can be calculated immediately when a statistical request is obtained.

[0121] See also Figure 4 , Figure 4 The following is a schematic diagram of a device for determining the image coverage area of ​​a remote sensing image provided by an exemplary embodiment of the present application. Figure 4 As shown in , the determining device 400 includes:

[0122] The first acquisition module 410 is used to acquire an associated remote sensing image that is spatially associated with a target statistical area; wherein the target statistical area corresponds to an associated area, the target statistical area is smaller than the associated area, the associated area is pre-divided into a predetermined number of statistical grids, each associated remote sensing image corresponds to a blank bitmap for storing information, and the number of bits of each blank bitmap is the same as the number of grids in the associated area;

[0123] The preprocessing module 420 is used to preprocess each associated remote sensing image to obtain coverage status information of the associated remote sensing image in the associated area.

[0124] The storage module 430 is used to store the coverage status information corresponding to the associated remote sensing image in a blank bitmap corresponding to the associated remote sensing image according to a preset storage order;

[0125] The second acquisition module 440 is used to obtain a bitmap corresponding to each associated remote sensing image that is spatially associated with the target statistical area in the predetermined time period in response to a request to count the image coverage area of ​​the target statistical area in the predetermined time period;

[0126] The statistical module 450 is used to determine the actual image coverage area of ​​the remote sensing image of the target statistical area within the predetermined time period based on the coverage status information on the bitmap corresponding to each associated remote sensing image.

[0127] In a possible implementation manner, the preprocessing module 420 is specifically used for:

[0128] For each associated remote sensing image, determine the image range block of the overlapping area between the associated remote sensing image and the target statistical area;

[0129] Overlay the image range block to the corresponding position on the statistical grid of the associated area to obtain an overlay grid;

[0130] Based on the center line of each row of grids in the superimposed grid, coverage status information of the associated area corresponding to the associated remote sensing image is obtained.

[0131] In a possible implementation manner, the preprocessing module 420 is specifically used for:

[0132] For each row of grids in the superimposed grid, pre-covered boundary grids are selected from the row of grids, where the boundary of the image range block intersects with the center line of the row of grids;

[0133] Based on the position serial number of each pre-covered boundary grid, all grids except the pre-covered boundary grid in the row of grids are screened to obtain covered grids among all grids except the pre-covered boundary grid; the position serial number indicates the order of each pre-covered boundary grid in the pre-covered boundary grids in the row of grids;

[0134] Based on the position sequence number of each pre-covered boundary grid and the center point of each pre-covered boundary grid, all pre-covered boundary grids are screened to obtain the covering grids in the pre-covered boundary grids;

[0135] Determine the other grids except the covered grids screened out from the row of grids as uncovered grids;

[0136] The first coverage status information is determined as the coverage status information of the covered grid, and the second coverage status information is determined as the coverage status information of the uncovered grid.

[0137] In a possible implementation manner, the preprocessing module 420 is specifically used for:

[0138] According to the order indicated by the position numbers of the pre-covered boundary grids, the pre-covered boundary grids corresponding to the odd position numbers and the pre-covered boundary grids corresponding to the even position numbers are searched from left to right in sequence;

[0139] For each pre-covered boundary grid corresponding to an odd-numbered position serial number, all grids between the pre-covered boundary grid and the pre-covered boundary grid corresponding to an even-numbered position serial number adjacent to the pre-covered boundary grid to the right are determined as covered grids.

[0140] In a possible implementation manner, the preprocessing module 420 is specifically used for:

[0141] For each pre-covered boundary grid in the row of grids, determining the parity of the position sequence number of the pre-covered boundary grid;

[0142] When the position sequence number of the pre-covering boundary grid is an odd number, determining whether the intersection point of the pre-covering boundary grid is on the left side of the center point of the pre-covering boundary grid;

[0143] If the intersection point of the pre-covering boundary grid is on the left side of the center point of the pre-covering boundary grid, the pre-covering boundary grid is determined as a covering grid;

[0144] When the position sequence number of the pre-covering boundary grid is an even number, determining whether the intersection point of the pre-covering boundary grid is on the right side of the center point of the pre-covering boundary grid;

[0145] If the intersection point of the pre-covering boundary grid is on the right side of the center point of the pre-covering boundary grid, the pre-covering boundary grid is determined as a covering grid.

[0146] In a possible implementation, the statistics module 450 is specifically used for:

[0147] According to the arrangement order of the position indexes corresponding to the bitmaps, the coverage status information of each bit on each bitmap corresponding to each position index is obtained in sequence;

[0148] For each position index, an OR operation is performed on the coverage status information of the corresponding bit on each bitmap corresponding to the position index to obtain the merged coverage status information corresponding to the position index; the merged coverage status information includes merged coverage and merged non-coverage;

[0149] Filter out the merged coverage status information of the merged coverage from all the merged coverage status information;

[0150] Based on the combined coverage state information of the combined coverage, the actual image coverage area of ​​the remote sensing image of the target statistical area within the predetermined time period is determined.

[0151] In a possible implementation, the statistics module 450 is specifically used for:

[0152] Determine a target position index corresponding to the merged coverage state information of the merged coverage;

[0153] Based on the correspondence between the preset position index and the statistical grid, a target statistical grid corresponding to the target position index is obtained;

[0154] Based on the correspondence between the preset statistical grids and the real geographical areas, the real geographical area corresponding to each target statistical grid is obtained; wherein the real geographical area corresponding to each statistical grid represents the area of ​​the real geographical region corresponding to the statistical grid;

[0155] The sum of the real geographic areas corresponding to each target statistical grid is determined as the real image coverage area of ​​the remote sensing image of the target statistical area in the predetermined time period.

[0156] An exemplary embodiment of the present application provides a device for determining the image coverage area of ​​a remote sensing image. When a remote sensing image spatially associated with a target statistical area is received, coverage status information corresponding to the remote sensing image is stored in the form of a bitmap. When counting the actual image coverage area of ​​the remote sensing image of the target statistical area within any predetermined time period, area statistics are directly performed based on the coverage status information in the pre-stored bitmap, so that when counting the actual image coverage area of ​​the target statistical area within any predetermined time period, the actual image coverage area of ​​the remote sensing image of the target statistical area within the predetermined time period can be calculated immediately when a statistical request is obtained.

[0157] See also Figure 5 , Figure 5 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. Figure 5 As shown in , the electronic device 500 includes a processor 510 , a memory 520 and a bus 530 .

[0158] The memory 520 stores machine-readable instructions executable by the processor 510. When the electronic device 500 is running, the processor 510 communicates with the memory 520 through the bus 530. When the machine-readable instructions are executed by the processor 510, the steps of the method for determining the image coverage area of ​​the remote sensing image in the above-mentioned method embodiment can be executed. The specific implementation method can be found in the method embodiment, which will not be repeated here.

[0159] An embodiment of the present application also provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the steps of the method for determining the image coverage area of ​​the remote sensing image in the above method embodiment can be executed. The specific implementation method can be found in the method embodiment and will not be repeated here.

[0160] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0161] In the several embodiments provided in the present application, it should be understood that the disclosed systems, devices and methods can be implemented in other ways. The device embodiments described above are merely schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some communication interfaces, and the indirect coupling or communication connection of devices or units can be electrical, mechanical or other forms.

[0162] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0163] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit.

[0164] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a non-volatile computer-readable storage medium that is executable by a processor. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions for a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present application. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM), random access memory (RAM), disk or optical disk, and other media that can store program codes.

[0165] Finally, it should be noted that the above-described embodiments are only specific implementation methods of the present application, which are used to illustrate the technical solutions of the present application, rather than to limit them. The protection scope of the present application is not limited thereto. Although the present application is described in detail with reference to the above-mentioned embodiments, ordinary technicians in the field should understand that any technician familiar with the technical field can still modify the technical solutions recorded in the above-mentioned embodiments within the technical scope disclosed in the present application, or can easily think of changes, or make equivalent replacements for some of the technical features therein; and these modifications, changes or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application, and should be included in the protection scope of the present application. Therefore, the protection scope of the present application shall be based on the protection scope of the claims.

Claims

1. A method for determining the image coverage area of ​​a remote sensing image, characterized in that: The determination method comprises: Acquire an associated remote sensing image that is spatially associated with a target statistical area; wherein the target statistical area corresponds to an associated area, the target statistical area is smaller than the associated area, the associated area is pre-divided into a predetermined number of statistical grids, each associated remote sensing image corresponds to a blank bitmap for storing information, and the number of bits of each blank bitmap is the same as the number of grids in the associated area; For each associated remote sensing image, preprocessing the associated remote sensing image is performed to obtain coverage status information corresponding to the associated remote sensing image in the associated area; The coverage status information corresponding to the associated remote sensing image is stored in a blank bitmap corresponding to the associated remote sensing image according to a preset storage order; In response to a request for counting the image coverage area of ​​a target statistical area in a predetermined time period, obtaining a bitmap corresponding to each associated remote sensing image that is spatially associated with the target statistical area in the predetermined time period; Determine the real image coverage area of ​​the remote sensing image of the target statistical area in the predetermined time period based on the coverage status information on the bitmap corresponding to each associated remote sensing image; The method of preprocessing each associated remote sensing image to obtain coverage status information of the associated remote sensing image in the associated area includes: For each associated remote sensing image, determine the image range block of the overlapping area between the associated remote sensing image and the target statistical area; Overlay the image range block to the corresponding position on the statistical grid of the associated area to obtain an overlay grid; Based on the center line of each row of grids in the superimposed grid, obtaining coverage status information of the associated remote sensing image in the associated area; The obtaining, based on the center line of each row of grids in the superimposed grid, coverage status information corresponding to the associated remote sensing image in the associated area includes: For each row of grids in the superimposed grid, pre-covered boundary grids are selected from the row of grids, where the boundary of the image range block intersects with the center line of the row of grids; Based on the position serial number of each pre-covered boundary grid, all grids except the pre-covered boundary grid in the row of grids are screened to obtain covered grids among all grids except the pre-covered boundary grid; the position serial number indicates the order of each pre-covered boundary grid in the pre-covered boundary grids in the row of grids; Based on the position sequence number of each pre-covered boundary grid and the center point of each pre-covered boundary grid, all pre-covered boundary grids are screened to obtain the covering grids in the pre-covered boundary grids; Determine the other grids except the covered grids screened out from the row of grids as uncovered grids; The first coverage status information is determined as the coverage status information of the covered grid, and the second coverage status information is determined as the coverage status information of the uncovered grid.

2. The determination method according to claim 1, characterized in that: The step of screening all grids in the row of grids except the pre-covered boundary grid based on the position serial number of each pre-covered boundary grid to obtain the covering grids in all grids except the pre-covered boundary grid comprises: According to the order indicated by the position numbers of the pre-covered boundary grids, the pre-covered boundary grids corresponding to the odd position numbers and the pre-covered boundary grids corresponding to the even position numbers are searched from left to right in sequence; For each pre-covered boundary grid corresponding to an odd-numbered position serial number, all grids between the pre-covered boundary grid and the pre-covered boundary grid corresponding to an even-numbered position serial number adjacent to the pre-covered boundary grid to the right are determined as covered grids.

3. The determination method according to claim 1, characterized in that: The step of screening all pre-covered boundary grids based on the position serial number of each pre-covered boundary grid and the center point of each pre-covered boundary grid to obtain the covering grid in the pre-covered boundary grid includes: For each pre-covered boundary grid in the row of grids, determining the parity of the position sequence number of the pre-covered boundary grid; When the position sequence number of the pre-covering boundary grid is an odd number, determining whether the intersection point of the pre-covering boundary grid is on the left side of the center point of the pre-covering boundary grid; If the intersection point of the pre-covering boundary grid is on the left side of the center point of the pre-covering boundary grid, the pre-covering boundary grid is determined as a covering grid; When the position sequence number of the pre-covering boundary grid is an even number, determining whether the intersection point of the pre-covering boundary grid is on the right side of the center point of the pre-covering boundary grid; If the intersection point of the pre-covering boundary grid is on the right side of the center point of the pre-covering boundary grid, the pre-covering boundary grid is determined as a covering grid.

4. The determination method according to claim 1, characterized in that: The method of determining the real image coverage area of ​​the remote sensing image of the target statistical area in the predetermined time period based on the coverage status information on the bitmap corresponding to each associated remote sensing image includes: According to the arrangement order of the position indexes corresponding to the bitmaps, the coverage status information of each bit on each bitmap corresponding to each position index is obtained in sequence; For each position index, an OR operation is performed on the coverage status information of the corresponding bit on each bitmap corresponding to the position index to obtain the merged coverage status information corresponding to the position index; the merged coverage status information includes merged coverage and merged non-coverage; Filter out the merged coverage status information of the merged coverage from all the merged coverage status information; Based on the combined coverage state information of the combined coverage, the actual image coverage area of ​​the remote sensing image of the target statistical area within the predetermined time period is determined.

5. The determination method according to claim 4, characterized in that: The step of determining the real image coverage area of ​​the remote sensing image of the target statistical area within the predetermined time period based on the combined coverage state information of the combined coverage includes: Determine a target position index corresponding to the merged coverage state information of the merged coverage; Based on the correspondence between the preset position index and the statistical grid, a target statistical grid corresponding to the target position index is obtained; Based on the correspondence between the preset statistical grids and the real geographical areas, the real geographical area corresponding to each target statistical grid is obtained; wherein the real geographical area corresponding to each statistical grid represents the area of ​​the real geographical region corresponding to the statistical grid; The sum of the real geographic areas corresponding to each target statistical grid is determined as the real image coverage area of ​​the remote sensing image of the target statistical area in the predetermined time period.

6. A device for determining image coverage area of ​​a remote sensing image, characterized in that: The determining device comprises: The first acquisition module is used to acquire an associated remote sensing image that is spatially associated with the target statistical area; wherein the target statistical area corresponds to an associated area, the target statistical area is smaller than the associated area, the associated area is pre-divided into a predetermined number of statistical grids, each associated remote sensing image corresponds to a blank bitmap for storing information, and the number of bits of each blank bitmap is the same as the number of grids in the associated area; A preprocessing module is used to preprocess each associated remote sensing image to obtain coverage status information of the associated remote sensing image in the associated area; A storage module, used for storing the coverage status information corresponding to the associated remote sensing image in a blank bitmap corresponding to the associated remote sensing image according to a preset storage order; A second acquisition module is used to respond to a request for counting the image coverage area of ​​a target statistical area in a predetermined time period, and to acquire a bitmap corresponding to each associated remote sensing image that is spatially associated with the target statistical area in the predetermined time period; A statistical module, for determining the real image coverage area of ​​the remote sensing image of the target statistical area within the predetermined time period based on the coverage status information on the bitmap corresponding to each associated remote sensing image; The preprocessing module is specifically used for: For each associated remote sensing image, determine the image range block of the overlapping area between the associated remote sensing image and the target statistical area; Overlay the image range block to the corresponding position on the statistical grid of the associated area to obtain an overlay grid; Based on the center line of each row of grids in the superimposed grid, obtaining coverage status information of the associated remote sensing image in the associated area; The preprocessing module is also specifically used for: For each row of grids in the superimposed grid, pre-covered boundary grids are selected from the row of grids, where the boundary of the image range block intersects with the center line of the row of grids; Based on the position serial number of each pre-covered boundary grid, all grids except the pre-covered boundary grid in the row of grids are screened to obtain covered grids among all grids except the pre-covered boundary grid; the position serial number indicates the order of each pre-covered boundary grid in the pre-covered boundary grids in the row of grids; Based on the position sequence number of each pre-covered boundary grid and the center point of each pre-covered boundary grid, all pre-covered boundary grids are screened to obtain the covering grids in the pre-covered boundary grids; Determine the other grids except the covered grids screened out from the row of grids as uncovered grids; The first coverage status information is determined as the coverage status information of the covered grid, and the second coverage status information is determined as the coverage status information of the uncovered grid.

7. An electronic device, characterized in that: include: A processor, a memory and a bus, wherein the memory stores machine-readable instructions executable by the processor, and when the electronic device is running, the processor and the memory communicate via the bus, and when the machine-readable instructions are run by the processor, the steps of the method for determining the image coverage area of ​​a remote sensing image as described in any one of claims 1 to 5 are executed.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method for determining the image coverage area of ​​a remote sensing image as claimed in any one of claims 1 to 5 are executed.

Images