CN114972482B - Remote sensing image application area statistical method - Google Patents
Remote sensing image application area statistical method Download PDFInfo
- Publication number
- CN114972482B CN114972482B CN202210404351.3A CN202210404351A CN114972482B CN 114972482 B CN114972482 B CN 114972482B CN 202210404351 A CN202210404351 A CN 202210404351A CN 114972482 B CN114972482 B CN 114972482B
- Authority
- CN
- China
- Prior art keywords
- remote sensing
- sensing image
- area
- data
- appointed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000007619 statistical method Methods 0.000 title abstract description 7
- 238000000034 method Methods 0.000 claims description 25
- 238000012545 processing Methods 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000013589 supplement Substances 0.000 claims description 4
- 230000001502 supplementing effect Effects 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 3
- 230000001788 irregular Effects 0.000 description 13
- 238000010586 diagram Methods 0.000 description 11
- 238000011835 investigation Methods 0.000 description 7
- 238000012544 monitoring process Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 238000007726 management method Methods 0.000 description 3
- 239000002689 soil Substances 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Abstract
The invention provides a remote sensing image application area statistical method, which is used for obtaining accurate statistical data of remote sensing image application areas through analysis of available remote sensing image data in a plurality of appointed geographic areas with appointed time. The acquisition capability and coverage rule of the domestic public welfare satellite remote sensing image are fully mastered, the domestic public welfare satellite remote sensing image is utilized to the maximum extent, and data reference can be provided for making and supplementing domestic and foreign commercial satellite remote sensing image plans in advance.
Description
Technical Field
The invention relates to the technical field of satellite remote sensing image processing, in particular to a method for counting application conditions of data of remote sensing images from different sources.
Background
Along with the rapid development of remote sensing technology, information technology and IT technology, the remote sensing image acquisition capability is obviously improved, the engineering application field of the remote sensing image is continuously expanded, the remote sensing image acquisition method is applied to the aspects of natural resource investigation and monitoring, supervision and law enforcement, homeland space planning, ecological environment protection and the like, and good social benefit and economic benefit are obtained.
The data source of the remote sensing image relates to different commercial or government institutions at home and abroad, and besides satellite data, the remote sensing image also comprises data acquired by utilizing an unmanned plane and other aviation platforms, wherein the satellite data comprises WorldView Gao Fenji columns in the United states, SPOT high-score series in France, BJ, GF high-score series in China and the like. Under the condition that the application of the remote sensing image is wider and wider, if the same remote sensing image data is repeatedly purchased by different department units, the overall expenditure of the country is greatly increased, so that the country gradually advances the overall acquisition and sharing application of the remote sensing image from the top-level designed height, the 2021 7 month and 28 days later, and the 'natural resource office remote sensing image overall sharing management method' printed by the natural resource office is a substantial measure, and aims to realize 'one-purpose' of the remote sensing image, promote the overall planning and sharing of the remote sensing image, fully exert the benefits of the remote sensing image, avoid repeated purchase and processing, reduce the waste of resources and greatly improve the application depth and value of the remote sensing image.
The existing remote sensing image application area statistics generally uses scenes or frames as units for statistics according to the habits of image management and processing, calculates the coverage area of each scene or frame of remote sensing image, and then performs accumulated statistics. There are two statistical approaches to this,
The statistics method is carried out by a remote sensing image acquisition management provider in a dominant mode, statistics is carried out from the angle of providing services for the remote sensing images, and is used for grasping the service range and the service capability of the remote sensing images.
The other is conducted by remote sensing image users, statistics is conducted from the application angle of the remote sensing images, and the statistics is used for grasping the condition of covering an application area of the remote sensing images so as to ensure that the seamless covering of the remote sensing images is oriented, and the statistics method generally does not consider the overlapping covering condition.
Redundancy exists in the application areas of the remote sensing images obtained by the existing two statistical modes, and more accurate data cannot be provided.
In major national condition national force investigation projects such as national soil investigation and national geographic national condition monitoring implemented by the government of China, the remote sensing images are required to be applied to carry out spatial geographic information acquisition, and the requirements of unified time points are met, namely, the investigation result reflects the national condition at which time point. Thus, for a spatial location, the applied remote sensing images need to have explicit time information to capture, rather than overlapping overlays blended together to make the result time-uncertainty. This design is significant for monitoring outcome popularization and application, and the user can determine which time point the spatial position information is. Thus, a full-coverage, non-overlapping and seamless one-edition remote sensing image and corresponding homeland space information result can be formed on the national scale. Meanwhile, the government projects are national overall guarantee image data sources, the domestic public welfare satellite remote sensing images are fully utilized, and areas which cannot meet the demands are supplemented and purchased, so that the acquisition cost of the remote sensing images can be reduced to the greatest extent. For the application requirements of the remote sensing images, a statistical method can be provided, and the coverage range of the existing remote sensing image resources can be obtained only by using less time in each project, so that the domestic public welfare satellite remote sensing images can be utilized to the maximum extent, and data reference can be provided for making and supplementing domestic and foreign commercial satellite remote sensing image plans in advance.
Disclosure of Invention
The present invention is directed to a method for applying area statistics to remote sensing images, so as to reduce or avoid the above-mentioned problems.
The invention provides a remote sensing image application area statistical method, which obtains accurate statistical data of remote sensing image application areas through analyzing available remote sensing image data in a plurality of appointed geographic areas at appointed time, and comprises the following steps:
And step A, selecting a plurality of standby remote sensing image data for a designated geographical area at a designated time, and respectively carrying out orthographic correction according to a unified coordinate system. The selected remote sensing image data covers the appointed geographical area, and the shooting time can be within a range of 90 days before and after the appointed time.
And B, selecting one digital orthographic image from the plurality of digital orthographic images obtained in the step A as a first reference image according to the priority of resolution, shooting time and quality, if the first reference image cannot fully cover the appointed geographic area, selecting a second reference image to supplement, and so on until the appointed geographic area is fully covered, generating a fully covered image, selecting the content of the reference image with high priority as the content of an overlapping area when the areas of the reference images overlap each other in the fully covered image, carrying out a matting processing on the areas with invalid value coverage and cloud coverage in the fully covered image, and sequentially selecting valid and cloud-free covered data in the reference images with low priority to replace.
And C, creating metadata of a vector data structure for polygons formed by boundaries of the digital orthographic images according to the full-coverage digital orthographic images of the specified geographic areas obtained in the step B, and calculating the areas of the polygons in the metadata to obtain the plane areas of the polygons under a unified coordinate system. Thus, the accurate data of the application area of the remote sensing image data in the step B is obtained.
And D, processing the available remote sensing image data in a plurality of appointed geographical areas at appointed time according to the methods from the step A to the step C, so that the accurate statistical data of the application area of the remote sensing image in a large range can be obtained.
Preferably, in step a, the coordinate system is a 2000 national earth coordinate system.
Preferably, in step a, the number of standby remote sensing image data is 2-5.
Preferably, in step B, the calculation formula of the priority is that,
Priority=k1/resolution+k2 (100-capturing time difference)/100+k3 mass
In the above, the resolution is a parameter of the data source of each standby remote sensing image data,
The shooting time difference is the date difference between the shooting time of each standby remote sensing image data and the appointed time,
And in the quality of each standby remote sensing image data, taking the area of the range covering the designated geographic area as a denominator, and taking the area of the part covered by the invalid value and covered by the cloud in the range covering the designated geographic area as a numerator, so as to obtain the percentage data.
K1, k2, k3 are adjustment constants for balancing the parameters, k1=0.4, k2=0.4, k3=0.2.
Preferably, in the step B, when the region having invalid value coverage and cloud coverage is processed, a straight line segment polygon mode is adopted. And carrying out line segmentation processing on boundary parts between different reference images.
Preferably, in step C, the metadata is in shapefile format.
Preferably, in step C, the metadata records the save resolution and the shooting time in addition to the application area.
The application area statistical method of the remote sensing image provided by the invention can accurately grasp the application situation of the remote sensing image, especially on the national scale, the actual effective coverage capacity and the actual application situation of the remote sensing image comprise key parameter information such as the resolution of the remote sensing image, shooting time, space coverage, application area and the like, and support overall acquisition, reasonable layout and scientific decision of the remote sensing image. The acquisition capability and coverage rule of the domestic public welfare satellite remote sensing image are fully mastered, the domestic public welfare satellite remote sensing image is utilized to the maximum extent, and data reference can be provided for making and supplementing domestic and foreign commercial satellite remote sensing image plans in advance.
Drawings
The following drawings are only for purposes of illustration and explanation of the present invention and are not intended to limit the scope of the invention. Wherein,
FIG. 1 is a schematic diagram of a designated geographic area of a remote sensing image using an area statistics method according to an embodiment of the present invention;
FIG. 2 is a digital orthophoto schematic of a remote sensing image covering the designated geographic area of FIG. 1;
FIG. 3 is a digital orthophoto schematic of another remote sensing image covering the designated geographic area of FIG. 1;
FIG. 4 is a schematic diagram showing the overlapping state of the digital orthographic images of FIGS. 2 and 3;
FIG. 5 is a schematic diagram of a full-coverage digital orthographic image of the designated geographic area of FIG. 1 obtained by the digital orthographic image adjustment of FIGS. 2 and 3;
FIG. 6 is a schematic diagram of the full-coverage digital orthophoto map of FIG. 5 after correction;
fig. 7 is a schematic diagram after the full-coverage digital orthophoto boundary vectorization of fig. 6.
Detailed Description
For a clearer understanding of technical features, objects, and effects of the present invention, a specific embodiment of the present invention will be described with reference to the accompanying drawings.
The invention provides a remote sensing image application area statistical method, which obtains accurate statistical data of remote sensing image application areas through analyzing available remote sensing image data in a plurality of appointed geographic areas at appointed time, and comprises the following steps:
And step A, selecting a plurality of standby remote sensing image data for a designated geographical area at a designated time, and respectively carrying out orthographic correction according to a unified coordinate system. The selected remote sensing image data covers the appointed geographical area, and the shooting time can be within a range of 90 days before and after the appointed time.
The main source of the remote sensing image data is satellite data, and different satellites have respective operation rules, so that the satellite remote sensing image data cannot be guaranteed to exist in a specific designated geographic area, and therefore, in practical project application, the shooting time of the satellite data can be widened to a time period of the first 90 days and the last 90 days of the designated geographic area, and the satellite data covering the designated geographic area can be guaranteed. In addition, since different satellite data may not necessarily ensure complete coverage for the specified geographic area, it is possible to cover only a portion of the specified geographic area, and in order to ensure operational efficiency, satellite data is selected to cover at least a range of 5% of the specified geographic area.
The mathematical basis of the remote sensing image data obtained by different satellites is different, and the major engineering has uniform requirements on a coordinate system, so that after a plurality of standby remote sensing image data are selected, a uniform coordinate system needs to be defined, and different satellite data related to the standby remote sensing image data can be subjected to orthographic correction according to the uniform coordinate system, so that a digital orthographic image taking a scene or a picture as a unit is generated, and the national common use is a 2000-country geodetic coordinate system.
With the development of the aerospace technology in China, satellite image data which can be provided by China are more and more abundant, in the major national condition national force investigation engineering such as national soil investigation and national geographic national condition monitoring, in order to facilitate visual analysis, high-resolution satellite image data with resolution not more than 2.5 meters and with focal area not more than 1 meter are mainly used as data supplement, and when the quantity of the high-resolution satellite image data is insufficient to cover a designated geographic area, satellite image data with resolution lower than 2.5 meters can also be used for standby remote sensing image data.
In order to ensure the acquisition integrity of the information of the appointed geographical area, theoretically, the more standby remote sensing image data are, the better, but the more data are, the more complicated operation is increased, and the inventor team can better consider the problems of data coverage and calculation balance by researching the invention, namely 2-5 standby remote sensing image data.
FIG. 1 is a schematic diagram of a designated geographic area of a remote sensing image using an area statistics method according to an embodiment of the present invention; FIG. 2 is a digital orthophoto schematic of a remote sensing image covering the designated geographic area of FIG. 1; FIG. 3 is a digital orthophoto schematic of another remote sensing image covering the designated geographic area of FIG. 1; referring to fig. 1-3, in a specific embodiment, a specific remote sensing image project needs to process data of a specific geographical area shown in fig. 1, a specific time point is 2021, 6 and 30 days, a team of inventors finds out two remote sensing images covering the specific geographical area shown in fig. 1 in remote sensing image data of 90 days around 6 and 30 days, the resolution of the two remote sensing images is 1 meter, the shooting time is 2021, 6 and 28 days, 2021, 4 and 15 days, and after orthorectified according to a 2000-country geodetic coordinate system, digital orthographic images shown in fig. 2 and 3 are obtained, respectively, the shooting time of fig. 2 is 2021, 6 and 28 days, and the shooting time of fig. 3 is 2021, 4 and 15 days.
And B, selecting one digital orthographic image from the plurality of digital orthographic images obtained in the step A as a first reference image according to the priority of resolution, shooting time and quality, if the first reference image cannot fully cover the appointed geographic area, selecting a second reference image to supplement, and so on until the appointed geographic area is fully covered, generating a fully covered image, selecting the content of the reference image with high priority as the content of an overlapping area when the areas of the reference images overlap each other in the fully covered image, carrying out a matting processing on the areas with invalid value coverage and cloud coverage in the fully covered image, and sequentially selecting valid and cloud-free covered data in the reference images with low priority to replace.
As described above, since different satellites have respective operation rules, the optimally selected satellite image does not necessarily cover the designated geographical area, and the remote sensing image data obtained by different satellites may have different parameters such as shooting angle, resolution, coverage band, etc., therefore, among the plurality of digital orthographic images obtained in step a, it is necessary to select the sequentially numbered reference images according to the priority of resolution, shooting time, quality, specifically, the priority of the standby remote sensing image data may be set according to the following formula,
Priority=k1/resolution+k2 (100-capturing time difference)/100+k3 mass
In the above description, the resolution is the parameter of the data source (satellite, aerial) of each standby remote sensing image data,
The photographing time difference is a date difference between the photographing time of each standby remote sensing image data and the specified time, for example, if the photographing time is the same day as the specified time, the photographing time difference is 0, if the photographing time is the first 30 days or the last 30 days of the specified time, the photographing time difference is 30,
And in the quality of each standby remote sensing image data, taking the area of the range covering the designated geographic area as a denominator, and taking the area of the part covered by the invalid value and covered by the cloud in the range covering the designated geographic area as a numerator, so as to obtain the percentage data.
K1, k2 and k3 are adjustment constants set by the inventor and used for balancing parameters, and in the research process of the inventor, the inventors find that when the following values are respectively set, the standby remote sensing image data can be balanced better, k1=0.4, k2=0.4 and k3=0.2.
Taking the digital orthographic images shown in fig. 2 and 3 as an example, the resolution of the two standby remote sensing image data shown in fig. 2 and 3 is 1 meter, the photographing time of fig. 2 is 2021, 6, 28, and the photographing time of fig. 3 is 2021, 4, 15. Through the calculation of the computer, the computer program,
The resolution of fig. 2 is 1 meter, the photographing time difference is 2 days, and the mass is about 92%, so the priority of fig. 2=0.4/1+0.4 (100-2)/100+0.2×92% =0.4+0.392+0.184=0.976.
The resolution of fig. 3 is 1 meter, the photographing time difference is 76 days, and the mass is about 99%, so the priority of fig. 3=0.4/1+0.4 (100-76)/100+0.2×99% =0.4+0.096+0.198=0.694.
Thus, for the designated geographic area shown in fig. 1, the digital orthographic image shown in fig. 2 may be selected as the first reference image and the digital orthographic image shown in fig. 3 may be selected as the second reference image.
Fig. 4 is a schematic diagram showing an overlapping state of the digital orthographic images of fig. 2 and 3, and fig. 5 is a schematic diagram showing a full coverage digital orthographic image of the designated geographic area of fig. 1 obtained by adjusting the digital orthographic images of fig. 2 and 3;
Referring to fig. 4, although the digital orthographic image shown in fig. 3 covers the designated geographical area more, the digital orthographic image shown in fig. 2 has higher priority, so that the digital orthographic image shown in fig. 2 is selected as the first reference image and the digital orthographic image shown in fig. 3 is the second reference image among the full-coverage images of the designated geographical area, which is mainly the digital orthographic image shown in fig. 5.
Referring to fig. 5, the white area in the middle of fig. 5 is an area covered by invalid values or cloud in the digital orthographic image shown in fig. 2, and for this area, a matting process is required, and the corresponding area data in the digital orthographic image of fig. 3 is used for replacement. The white area in fig. 5 has irregular boundaries, and if the white area is scratched out according to the irregular boundaries, the operation efficiency is greatly affected.
FIG. 6 is a schematic diagram of the full-coverage digital orthophoto map of FIG. 5 after correction; for the irregular area to be scratched out in fig. 5, the inventor has selected to process in a polygonal manner, specifically, the boundary end point of the irregular area may be first marked on a computer, and as shown in fig. 5 and 6, the upper right corner end point and the lower right corner end point may be first marked, and these two end points may be connected, if all the right side boundary of the irregular area is covered, the connection is selected as a right side straight line boundary, then the lower left corner end point of the irregular area is selected, the lower left corner end point and the lower right corner end point may be connected, if all the lower boundary of the irregular area is covered, then the connection is selected as a lower straight line boundary, then the upper left corner end point of the irregular area is selected, and the upper left corner end point and the upper right corner end point may be connected, if all the upper boundary of the irregular area is covered, then the connection line is selected as an upper straight line boundary, then the upper left corner end point and the lower left corner end point are connected, and as the connection line can not cover all left side boundaries of the irregular area, the leftmost end point of the irregular area is marked again, then the leftmost end point is respectively connected with the upper left corner end point and the lower left corner end point, the coverage condition is detected, the connection line of the leftmost end point and the lower left corner end point can cover the boundary of the irregular area, so that the connection line is reserved as the lower left straight line boundary, and the connection line of the leftmost end point and the upper left corner end point can not cover the boundary of the irregular area, therefore, the point with the farthest connection line distance can be selected as a correction end point in the vertical direction of the connection line of the leftmost end point and the upper left corner end point, and the upper left corner end point are respectively connected to obtain the upper left straight line boundary.
As described above, for the region covered by the invalid value or the cloud in the first reference image, after analysis is performed according to the boundary nodes of the regions, a straight-line segment polygon for cutting out may be obtained, and according to the above method, generally, the sides of the polygon will not exceed eight sides. Therefore, the operation efficiency of the computer can be greatly improved.
In addition, the boundary part between different reference images can be subjected to line segmentation processing, so that the processing efficiency of a subsequent computer can be improved. The operation of the line segmentation process may be performed by using computer software, specifically, referring to fig. 6, for the boundary between the first reference image shown in fig. 2 and the second reference image shown in fig. 3, the automatic capturing function of ArcGIS software may be used, and the boundary may be similarly adjusted to a line segment boundary composed of a plurality of straight lines.
And C, creating metadata of a vector data structure for polygons formed by boundaries of the digital orthographic images according to the full-coverage digital orthographic images of the specified geographic areas obtained in the step B, and calculating the areas of the polygons in the metadata to obtain the plane areas of the polygons under a unified coordinate system. Thus, the accurate data of the application area of the remote sensing image data in the step B is obtained.
Fig. 7 is a schematic diagram after the full-coverage digital orthophoto boundary vectorization of fig. 6. Fig. 7 is a schematic diagram, in which boundary points of polygons of right images are not all shown, and as shown in fig. 6 and 7, as described above, the boundaries between different reference images are also line segment boundaries composed of a plurality of straight lines, except that the cutout portions are polygons, so that the whole full-coverage digital orthographic image is composed of three polygons, which is very beneficial to creating metadata of a vector data structure. Specifically, metadata created in shape format may be used. After the metadata of the vector data structure is created, the area of each polygon in the metadata can be calculated by a computer easily, and the plane area of each polygon under a unified coordinate system can be obtained. Thus, each reference image, namely the application area of each standby remote sensing image data, is obtained.
Since the vector data structure is adopted, basic information such as resolution, shooting time and the like can be recorded besides the application area.
And D, processing available remote sensing image data in a plurality of appointed geographical areas at appointed time according to the methods from the step A to the step C to obtain a large amount of metadata, and thus obtaining accurate statistical data of a large-scale remote sensing image application area.
After obtaining a large amount of metadata in the step C, the resolution can be classified (for example, the resolution is classified into 'better than 1 meter', '1 meter to 2.5 meters', 'lower than 2.5 meters') according to the resolution of the remote sensing image block recorded by the metadata, and the classification result is recorded into a resolution level attribute item, namely, each polygon is marked with a resolution level, and according to the mark, the application area cumulative value of the remote sensing image of each resolution level in the whole country can be obtained through statistics.
In addition, the photographing time may be segmented (for example, divided into "2021 year 1-3 months", "2021 year 4-6 months", "2021 year 7-9 months") according to the photographing time of the remote sensing image block recorded by the metadata, and the segmentation result may be recorded in the photographing time period attribute item, that is, the photographing time period is marked for each polygon, and the remote sensing image application area cumulative value of each photographing time period in the whole country is obtained by statistics according to the mark.
The area accumulated value of the corresponding attribute in the whole country can be obtained by statistics from the angles of other attributes according to the needs.
The metadata of the vector data structure obtained in the step C can easily store the accurate data of the parameters of the application area, the resolution and the shooting time. Therefore, as described in the background art, for massive remote sensing image application scenes in projects such as national soil investigation and national geographic state monitoring which are implemented by the government, in different project projects of the national scale, the application situation of the remote sensing image can be accurately mastered by adopting the method of the application, especially on the national scale, the actual effective coverage capability and the actual application situation of the remote sensing image, including key parameter information such as remote sensing image resolution, shooting time, space coverage, application area and the like, support overall acquisition, reasonable layout and scientific decision of the remote sensing image. The acquisition capability and coverage rule of the domestic public welfare satellite remote sensing image are fully mastered, the domestic public welfare satellite remote sensing image is utilized to the maximum extent, and data reference can be provided for making and supplementing domestic and foreign commercial satellite remote sensing image plans in advance.
It should be understood by those skilled in the art that while the present invention has been described in terms of several embodiments, not every embodiment contains only one independent technical solution. The description is given for clearness of understanding only, and those skilled in the art will understand the description as a whole and will recognize that the technical solutions described in the various embodiments may be combined with one another to understand the scope of the present invention.
The foregoing is illustrative of the present invention and is not to be construed as limiting the scope of the invention. Any equivalent alterations, modifications and combinations thereof will be effected by those skilled in the art without departing from the spirit and principles of this invention, and it is intended to be within the scope of the invention.
Claims (6)
1. The method is characterized in that the accurate statistical data of the application area of the remote sensing image is obtained by analyzing the available remote sensing image data in a specified geographic area at a plurality of specified times, and the method comprises the following steps:
Step A, selecting a plurality of standby remote sensing image data for a designated geographic area at a designated time, and respectively carrying out orthographic correction according to a unified coordinate system; the selected remote sensing image data cover the appointed geographical area, and the shooting time is within the range of 90 days before and after the appointed time;
B, selecting one digital orthographic image from the plurality of digital orthographic images obtained in the step A as a first reference image according to the priority of resolution, shooting time and quality, if the first reference image cannot fully cover the appointed geographic area, selecting a second reference image to supplement, and so on until the appointed geographic area is fully covered, generating a fully covered image, selecting the content of the reference image with high priority as the content of an overlapping area when the areas of the reference images overlap each other in the fully covered image, carrying out a matting processing on the areas with invalid value coverage and cloud coverage in the fully covered image, and sequentially selecting valid and cloud-free covered data in the reference images with low priority to replace;
The calculation formula of the priority is that,
Priority=k1/resolution+k2 (100-capturing time difference)/100+k3 mass
In the above, the resolution is a parameter of the data source of each standby remote sensing image data,
The shooting time difference is the date difference between the shooting time of each standby remote sensing image data and the appointed time,
In the quality of each standby remote sensing image data, taking the area of the range covering the appointed geographic area as a denominator, and taking the area of the part covered by the invalid value and covered by the cloud in the range covering the appointed geographic area as a numerator, so as to obtain percentage data;
k1, k2, k3 are tuning constants for balancing the parameters, k1=0.4, k2=0.4, k3=0.2;
Step C, creating metadata of a vector data structure for polygons formed by boundaries of the digital orthographic images according to the full-coverage digital orthographic images of the specified geographic areas obtained in the step B, and calculating the areas of the polygons in the metadata to obtain the plane areas of the polygons under a unified coordinate system; thereby obtaining the accurate data of the application area of the remote sensing image data in the step B;
And D, processing the available remote sensing image data in a plurality of appointed geographical areas at appointed time according to the methods from the step A to the step C, so that the accurate statistical data of the application area of the remote sensing image in a large range can be obtained.
2. The method according to claim 1, wherein in step a, the coordinate system is a 2000-country geodetic coordinate system.
3. The method of claim 1, wherein in step a, the number of backup remote sensing image data is 2-5.
4. The method according to claim 1, wherein in step B, when the region having invalid value coverage and cloud coverage is processed, a linear segment polygon method is used to segment the boundary between different reference images.
5. The method according to claim 1, wherein in step C, the metadata is in shape format.
6. The method according to claim 1, wherein in step C, the metadata is recorded with a storage resolution and a photographing time in addition to an application area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210404351.3A CN114972482B (en) | 2022-04-18 | Remote sensing image application area statistical method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210404351.3A CN114972482B (en) | 2022-04-18 | Remote sensing image application area statistical method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114972482A CN114972482A (en) | 2022-08-30 |
| CN114972482B true CN114972482B (en) | 2024-06-21 |
Family
ID=
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106327452A (en) * | 2016-08-14 | 2017-01-11 | 曾志康 | Fragmented remote sensing image synthesis method and device for cloudy and rainy region |
| CN107527014A (en) * | 2017-07-20 | 2017-12-29 | 武汉珈和科技有限公司 | Crops planting area RS statistics scheme of sample survey design method at county level |
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106327452A (en) * | 2016-08-14 | 2017-01-11 | 曾志康 | Fragmented remote sensing image synthesis method and device for cloudy and rainy region |
| CN107527014A (en) * | 2017-07-20 | 2017-12-29 | 武汉珈和科技有限公司 | Crops planting area RS statistics scheme of sample survey design method at county level |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109949399B (en) | Scene three-dimensional reconstruction method based on unmanned aerial vehicle aerial image | |
| US8260006B1 (en) | System and method of aligning images | |
| US11682170B2 (en) | Generating three-dimensional geo-registered maps from image data | |
| US20070269102A1 (en) | Method and System of Generating 3D Images with Airborne Oblique/Vertical Imagery, GPS/IMU Data, and LIDAR Elevation Data | |
| CN112559534B (en) | Remote sensing image data filing management system and method | |
| CN112765095B (en) | Method and system for filing image data of stereo mapping satellite | |
| CN112113544B (en) | Remote sensing data processing method and system based on unmanned aerial vehicle image | |
| CN112561832B (en) | Remote sensing image data storage method and system | |
| Gruen | Reality-based generation of virtual environments for digital earth | |
| CN103606188A (en) | Geographical information on-demand acquisition method based on image point cloud | |
| CN113837892A (en) | Slope farmland dividing method based on 3S technology | |
| CN110660125B (en) | Three-dimensional modeling device for power distribution network system | |
| CN110659369B (en) | On-orbit high-precision lightweight global image control point database construction method and system | |
| CN104657968A (en) | Automatic vehicle-mounted three-dimensional laser point cloud facade classification and outline extraction method | |
| CN105005580A (en) | Method and device for displaying reservoir terrain | |
| CN112527916A (en) | Grid visual definition and intelligent feature analysis method based on remote sensing image | |
| Qiu et al. | Design and development of a web‐based interactive twin platform for watershed management | |
| CN111986074A (en) | Real projective image manufacturing method, device, equipment and storage medium | |
| CN109035365B (en) | Mosaic processing method of high-resolution image | |
| Li et al. | OmniCity: Omnipotent city understanding with multi-level and multi-view images | |
| CN113268533A (en) | Method for making two-label data in two-label four-real data platform | |
| CN114972482B (en) | Remote sensing image application area statistical method | |
| CN115357675B (en) | Method and system for establishing image control point database through standardized processing of image control points | |
| Gu et al. | Surveying and mapping of large-scale 3D digital topographic map based on oblique photography technology | |
| CN114972482A (en) | Application area statistical method for remote sensing image |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |