CN116309915A - Remote sensing image real-time rendering method, device and storage medium based on B/S architecture - Google Patents

Remote sensing image real-time rendering method, device and storage medium based on B/S architecture Download PDF

Info

Publication number
CN116309915A
CN116309915A CN202310285692.8A CN202310285692A CN116309915A CN 116309915 A CN116309915 A CN 116309915A CN 202310285692 A CN202310285692 A CN 202310285692A CN 116309915 A CN116309915 A CN 116309915A
Authority
CN
China
Prior art keywords
remote sensing
sensing image
data
optical remote
image data
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.)
Pending
Application number
CN202310285692.8A
Other languages
Chinese (zh)
Inventor
张龙星
张奇
孟祥国
吴洋
卢传芳
刘斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Space Star Technology Co Ltd
Original Assignee
Space Star Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Space Star Technology Co Ltd filed Critical Space Star Technology Co Ltd
Priority to CN202310285692.8A priority Critical patent/CN116309915A/en
Publication of CN116309915A publication Critical patent/CN116309915A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D10/00—Energy efficient computing, e.g. low power processors, power management or thermal management

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)

Abstract

The invention relates to a remote sensing image real-time rendering method, equipment and storage medium based on a B/S architecture, wherein the remote sensing image real-time rendering method comprises the following steps: requesting tile data by a foreground; reading optical remote sensing image data attribute information; establishing a conversion geometric relationship between the optical remote sensing image data and the PNG tile data; calculating the range intersection of the optical remote sensing image data and the PNG tile data; performing data processing on the read optical remote sensing image data block data; and constructing PNG tile data, and returning the generated tile data to the foreground. The invention solves the problem that the B/S architecture can not render and display the standard scene optical remote sensing images of high-resolution series, resource series and the like of level 1 and more than level 1 in real time, and achieves the purposes of needing no tile cutting of the optical remote sensing image data under the emergency condition and real-time interpretation condition, increasing the display mode of the optical remote sensing image data under the B/S architecture, realizing the optical remote sensing image real-time rendering and displaying method meeting the rendering and displaying rules of the B/S architecture, and improving the working efficiency of operators.

Description

Remote sensing image real-time rendering method, device and storage medium based on B/S architecture
Technical Field
The invention relates to the technical fields of remote sensing image interpretation and intelligent image interpretation, in particular to a remote sensing image real-time rendering method system, equipment and storage medium based on a B/S architecture.
Background
The remote sensing image real-time rendering technology is widely applied to remote sensing GIS software and has become a basic function in the software. The remote sensing GIS software is generally divided into desktop software, webGIS software and mobile software; the desktop software has the most powerful function and can process, analyze and display various data; the WebGIS terminal software and functions are weaker, and can display and simplify analysis processing on part of standardized data; the mobile terminal software is similar to WebGIS software, has weaker functions, can display partial standardized data and simplify analysis processing, and is more suitable for mobile devices such as mobile phones, PADs and the like.
With the development of WebGIS, more and more remote sensing GIS applications are developed based on WebGIS, so that the WebGIS applications are all put together. Meanwhile, image interpretation and image interpretation software developed and designed based on desktop software architecture also tend to develop to WebGIS application. Before the existing WebGIS platform loads the optical remote sensing image, the optical remote sensing image needs to be subjected to tile generation processing in advance, then tile data generated by processing is loaded, the purpose of loading the optical remote sensing image data is achieved, the optical remote sensing image data cannot be directly loaded, and real-time image stretching and tone adjustment cannot be carried out on the loaded optical remote sensing image data, so that the WebGIS platform cannot meet the development of image interpretation software to the WebGIS direction.
With the continuous progress of the remote sensing GIS field and the computer technology, the WebGIS platform can bear more data display and processing tasks, and compared with the existing WebGIS platform, the data display and processing processes of the WebGIS platform can be respectively placed on a foreground and a background so as to achieve the improvement of the image display capability of the WebGIS platform, which is one of the development directions of the next generation WebGIS platform. With the new demands and new requirements of image interpretation and image interpretation users, the development of image interpretation and image interpretation software by using a WebGIS platform has become an important research direction. Therefore, the research of the remote sensing image real-time rendering technology based on the B/S architecture becomes an important way for improving the capacity of the WebGIS platform, and has become the focus of the research of various manufacturers.
Disclosure of Invention
In order to solve the technical problems in the prior art, the invention aims to provide a remote sensing image real-time rendering method, equipment and a storage medium based on a B/S architecture, which can solve the problem that an optical remote sensing image cannot be rendered and displayed in real time in the display process of the optical remote sensing image under the B/S architecture, and an algorithm process for generating PNG tile data in real time is provided by establishing a geometric relationship between the optical remote sensing image data and the PNG tile data, and can meet the requirements of real-time display, image stretching and color adjustment of the optical remote sensing image data under the B/S architecture, effectively utilize computer resources, improve the capability of a WebGIS platform, and meet the requirements of image interpretation and image interpretation users on development of the WebGIS platform.
In order to achieve the above object, the present invention provides a remote sensing image real-time rendering method based on a B/S architecture, comprising the following steps:
step S1, requesting tile data by a foreground;
s2, reading attribute information of the optical remote sensing image data;
s3, establishing a conversion geometric relationship between the optical remote sensing image data and the PNG tile data;
s4, calculating a range intersection of the optical remote sensing image data and the PNG tile data;
s5, performing data processing on the read optical remote sensing image data block data;
and S6, constructing PNG tile data, and returning the generated tile data to the foreground.
In accordance with one aspect of the present invention, in the step S1, the foreground requests the required tile data in the form of an Http request,the request information comprises an optical remote sensing image data path and a tile data longitude and latitude range E w 。
According to one aspect of the present invention, in the step S2, the optical remote sensing image data is opened, and the width W of the optical remote sensing image data is read s High H s Band number B s Geographic reference information, RPC parameters.
According to one aspect of the present invention, in the step S3, if the optical remote sensing image data has the geographic reference information, a coordinate projection conversion geometric relationship between the optical remote sensing image data and the tile data is established; if the RPC parameters exist in the optical remote sensing image data, establishing an RPC conversion geometric relationship between the optical remote sensing image data and the tile data;
establishing a coordinate projection conversion geometric relation: setting a pixel coordinate (i) of a certain point A of the optical remote sensing image data s ，j s ) Corresponding projection coordinates (X s ，Y s ) Corresponding geographic coordinates (L s ，B s ) Affine transformation six parameters (a s0 ，a s1 ，a s2 ，b s0 ，b s1 ，b s2 ) The coordinate reference is in WKT form, and the conversion relationship between the pixel coordinate and the projection coordinate is as follows:
Figure BDA0004139756280000031
the conversion relation between the projection coordinates and the geographic coordinates can be created by using the OGRCoordinates conversion information in the GDAL in a coordinate reference WKT form, and the conversion of the projection coordinates and the geographic coordinates is performed by using a function conversion form;
establishing an RPC conversion geometric relation: setting a pixel coordinate (i) of a certain point A of the optical remote sensing image data s ，j s ) Corresponding projection coordinates (X s ，Y s ) Corresponding geographic coordinates (L s ，B s ) The conversion relation between the pixel coordinates and the geographic coordinates uses GDALCREATE RPCTransformer in GDAL to construct an RPC parameter conversion model, and the pixel coordinates and the geographic coordinates are converted by the function GDALRPCTransformConversion is performed between geographic coordinates.
According to an aspect of the present invention, in the step S4, specifically includes:
calculating the geographic range E of the optical remote sensing image data through the conversion geometric relationship established in the step S3 s Intersection E of tile data and optical remote sensing image data geographic range c The calculation formula of (2) is as follows:
E c ＝E w ∩E s
wherein E is w Represented as tile data latitude and longitude ranges.
According to an aspect of the present invention, in the step S5, specifically includes:
step S51, back-calculating the intersection E of the geographic ranges according to the conversion geometric relationship established in the step S3 c Pixel coordinates on the optical remote sensing image;
step S52, reading optical remote sensing image block data according to the pixel coordinate range on the optical remote sensing image;
step S53, performing data value-free processing on the read optical remote sensing image block data;
and step S54, performing one mode of maximum and minimum stretching, percentage truncation, histogram equalization, gama enhancement or standard deviation stretching on the data processed without the data value.
According to an aspect of the present invention, in the step S53, the data value processing includes: if the optical remote sensing image data does not have a data-free value, the pixel value 0 is regarded as the data-free value; if no data value exists, the output pixel value is assigned to 0, and the output tile pixel value formula is as follows:
Figure BDA0004139756280000041
wherein, the original pixel value of the input optical remote sensing image data is P s The output tile pixel value is P d 。
According to one aspect of the invention, in said step S54,
the maximum minimum stretching includes: maximum and minimum stretching and 16-bit to 8-bit processing, wherein the formula is as follows:
Figure BDA0004139756280000042
the percentage cutoff includes: s is intercepted by the left end and the right end of the histogram frequency N respectively by using a percentage cut-off value p Calculating a new pixel maximum p' m And pixel minimum value P' n The output tile pixel value is calculated using the following formula:
Figure BDA0004139756280000043
the histogram equalization includes: calculating the gray cumulative distribution frequency S of each level of the histogram according to the histogram level and the histogram frequency L The calculation formula is as follows:
Figure BDA0004139756280000044
calculating histogram equalized pixel value P using gray scale cumulative distribution frequency z And finally, processing by using maximum and minimum stretching, wherein the calculation formula is as follows:
P z ＝S L *(L-1)+P n
Figure BDA0004139756280000045
the Gama enhancement includes: the original pixel value is normalized, gamma processing and inverse normalization processing are carried out, and finally maximum and minimum stretching processing is carried out, wherein the calculation formula is as follows:
Figure BDA0004139756280000051
wherein t is a normalized value, t' is a Gamma-processed value, P z Is an inverse normalized value;
the standard deviation stretching includes: calculating new pixel maximum and minimum values by using the standard deviation and the average value of the pixels of the image, and finally performing maximum and minimum stretching treatment, wherein the calculation formula is as follows:
Figure BDA0004139756280000052
wherein the maximum value of the pixel value of the optical remote sensing image data is P m The minimum value of the pixel value is P n The average value of the pixel values is P p The mean square error of the pixel value is P σ The pixel has no data value P o The histogram level is P L N is the total number of pixels, and the histogram frequency is [ N ] 1 ,N 2 ......N L ]Gama value of G a The percentage cut-off value is S p Standard deviation ratio S s The method comprises the steps of carrying out a first treatment on the surface of the Let the original pixel value of the input optical remote sensing image data be P s The output tile pixel value is P d 。
According to an aspect of the present invention, there is provided an electronic apparatus including: one or more processors, one or more memories, and one or more computer programs; the processor is connected with the memory, the one or more computer programs are stored in the memory, and when the electronic device runs, the processor executes the one or more computer programs stored in the memory, so that the electronic device executes a remote sensing image real-time rendering method based on the B/S architecture according to any one of the technical schemes.
According to an aspect of the present invention, there is provided a computer readable storage medium storing computer instructions that, when executed by a processor, implement a B/S architecture-based remote sensing image real-time rendering method according to any one of the above technical solutions.
Compared with the prior art, the invention has the following beneficial effects:
according to the remote sensing image real-time rendering method system, equipment and storage medium based on the B/S architecture, real-time generation of optical remote sensing image tile data is achieved, a conversion geometric relation between the optical remote sensing image data and PNG tile data is established in the background, the range of an optical remote sensing image data block to be read is calculated according to the PNG tile data range and the conversion geometric relation, reading and processing are carried out on the optical remote sensing image data block, conversion geometric relation conversion at the pixel level is achieved, PNG tile data is built and generated, the problem that the B/S architecture cannot conduct real-time rendering display on standard scene optical remote sensing images such as high-resolution series and resource series of 1 level and above is solved, the problem that the optical remote sensing image data are not required to be cut into tiles under the emergency condition and real-time interpretation condition is achieved, the display mode of the optical remote sensing image data under the B/S architecture is increased, the optical remote sensing image real-time rendering display method meeting the rendering display rule of the B/S architecture is achieved, and the working efficiency of operators is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required to be used in the embodiments will be briefly described below. It is apparent that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained from these drawings without inventive effort for a person of ordinary skill in the art.
FIG. 1 schematically shows a flowchart of a remote sensing image real-time rendering method based on a B/S architecture according to an embodiment of the present invention;
FIG. 2 schematically shows that the percentage cut-off value cuts S at the left and right ends of the histogram frequency N, respectively, in the embodiment of the invention p Schematic of (2);
FIG. 3 is a flow chart schematically illustrating a method for rendering remote sensing images in real time based on a B/S architecture according to another embodiment of the present invention;
fig. 4 schematically shows an effect diagram of a remote sensing image real-time rendering method based on a B/S architecture according to an embodiment of the present invention.
Detailed Description
The description of the embodiments of this specification should be taken in conjunction with the accompanying drawings, which are a complete description of the embodiments. In the drawings, the shape or thickness of the embodiments may be enlarged and indicated simply or conveniently. Furthermore, portions of the structures in the drawings will be described in terms of separate descriptions, and it should be noted that elements not shown or described in the drawings are in a form known to those of ordinary skill in the art.
Any references to directions and orientations in the description of the embodiments herein are for convenience only and should not be construed as limiting the scope of the invention in any way. The following description of the preferred embodiments will refer to combinations of features, which may be present alone or in combination, and the invention is not particularly limited to the preferred embodiments. The scope of the invention is defined by the claims.
As shown in fig. 1 and 4, the remote sensing image real-time rendering method based on the B/S architecture of the present invention includes the following steps:
step S1, requesting tile data by a foreground;
s2, reading attribute information of the optical remote sensing image data;
s3, establishing a conversion geometric relationship between the optical remote sensing image data and the PNG tile data;
s4, calculating a range intersection of the optical remote sensing image data and the PNG tile data;
s5, performing data processing on the read optical remote sensing image data block data;
and S6, constructing PNG tile data, and returning the generated tile data to the foreground.
In this embodiment, according to the B/S architecture tile display rule, the optical remote sensing image data is rendered in real time into 256-pixel PNG tile data in the background, and the rendered PNG tiles are transmitted to the foreground view for display as required. In order to realize real-time generation of the optical remote sensing image tile data, a conversion geometric relation between the optical remote sensing image data and the PNG tile data is established in the background, the range of an optical remote sensing image data block to be read is calculated according to the PNG tile data range and the conversion geometric relation, the optical remote sensing image data block is read and processed, pixel-level conversion geometric relation conversion is realized, PNG tile data is constructed and generated, the problem that a B/S architecture cannot render and display standard scene optical remote sensing images such as 1-level and more than 1-level high-resolution series, resource series and the like in real time is solved, the situation that the optical remote sensing image data is cut into tiles under emergency conditions and real-time interpretation is achieved, the display mode of the optical remote sensing image data under a B/S architecture is increased, and the real-time rendering and displaying method of the optical remote sensing image meeting the rendering and displaying rule of the B/S architecture is realized, and the working efficiency of operators is improved.
In one embodiment of the present invention, preferably, in step S1, the foreground requests the required tile data in an Http request manner, and the request information includes an optical remote sensing image data path, and a latitude and longitude range E of the tile data w 。
In one embodiment of the present invention, preferably, in step S2, the optical remote sensing image data is opened, and the width W of the optical remote sensing image data is read s High H s Band number B s Geographic reference information, RPC parameters.
In one embodiment of the present invention, preferably, in step S3, if the optical remote sensing image data has the geographic reference information, a coordinate projection transformation geometric relationship between the optical remote sensing image data and the tile data is established; if the RPC parameters exist in the optical remote sensing image data, establishing an RPC conversion geometric relationship between the optical remote sensing image data and the tile data;
establishing a coordinate projection conversion geometric relation: setting a pixel coordinate (i) of a certain point A of the optical remote sensing image data s ，j s ) Corresponding projection coordinates (X s ，Y s ) Corresponding geographic coordinates (L s ，B s ) Affine transformation six parameters (a s0 ，a s1 ，a s2 ，b s0 ，b s1 ，b s2 ) The coordinate reference is in WKT form, and the conversion relationship between the pixel coordinate and the projection coordinate is as follows:
Figure BDA0004139756280000081
the conversion relation between the projection coordinates and the geographic coordinates can be created by using the OGRCoordinates conversion information in the GDAL in a coordinate reference WKT form, and the conversion of the projection coordinates and the geographic coordinates is performed by using a function conversion form;
establishing an RPC conversion geometric relation: setting a pixel coordinate (i) of a certain point A of the optical remote sensing image data s ，j s ) Corresponding projection coordinates (X s ，Y s ) Corresponding geographic coordinates (L s ，B s ) The conversion relation between the pixel coordinates and the geographic coordinates uses GDALCREATE RPCTransformer in GDAL to construct an RPC parameter conversion model, and the conversion between the pixel coordinates and the geographic coordinates is carried out through a function GDALRPCTransform.
In the embodiment, the optical remote sensing image data block is read, rendered and displayed by constructing the conversion geometric relation between the optical remote sensing image data and the PNG tile data, so that the optical remote sensing image data can be directly rendered and displayed on the existing B/S architecture software.
In one embodiment of the present invention, preferably, in step S4, the method specifically includes:
calculating the geographic range E of the optical remote sensing image data through the conversion geometric relationship established in the step S3 s Intersection E of tile data and optical remote sensing image data geographic range c The calculation formula of (2) is as follows:
E c ＝E w ∩E s
wherein E is w Represented as tile data latitude and longitude ranges.
In one embodiment of the present invention, preferably, in step S5, the method specifically includes:
step S51, back-calculating the intersection E of the geographic ranges according to the conversion geometric relationship established in the step S3 c Pixel coordinates on the optical remote sensing image;
step S52, reading optical remote sensing image block data according to the pixel coordinate range on the optical remote sensing image;
step S53, performing data value-free processing on the read optical remote sensing image block data;
and step S54, performing one mode of maximum and minimum stretching, percentage truncation, histogram equalization, gama enhancement or standard deviation stretching on the data processed without the data value.
In the embodiment, the real-time image stretching of the optical remote sensing image on the B/S architecture software is realized by performing data value-free processing on the optical remote sensing image data block and performing one of maximum and minimum stretching, percentage truncation, histogram equalization, gama enhancement and standard deviation stretching.
In one embodiment of the present invention, preferably, in step S53, the data value processing includes: if the optical remote sensing image data does not have a data-free value, the pixel value 0 is regarded as the data-free value; if no data value exists, the output pixel value is assigned to 0, and the output tile pixel value formula is as follows:
Figure BDA0004139756280000091
wherein, the original pixel value of the input optical remote sensing image data is P s The output tile pixel value is P d 。
In one embodiment of the present invention, preferably, in step S54,
maximum minimum stretching includes: maximum and minimum stretching and 16-bit to 8-bit processing, wherein the formula is as follows:
Figure BDA0004139756280000092
as shown in fig. 2, the percentage cutoff includes: s is intercepted by the left end and the right end of the histogram frequency N respectively by using a percentage cut-off value p Calculate a new pixel maximum P' m And pixel minimum value P' n The output tile pixel value is calculated using the following formula:
Figure BDA0004139756280000093
histogram equalization includes: calculating the gray cumulative distribution frequency S of each level of the histogram according to the histogram level and the histogram frequency L The calculation formula is as follows:
Figure BDA0004139756280000094
calculating histogram equalized pixel value P using gray scale cumulative distribution frequency z And finally, processing by using maximum and minimum stretching, wherein the calculation formula is as follows:
P z ＝S L *(L-1)+P n
Figure BDA0004139756280000101
gama enhancement includes: the original pixel value is normalized, gamma processing and inverse normalization processing are carried out, and finally maximum and minimum stretching processing is carried out, wherein the calculation formula is as follows:
Figure BDA0004139756280000102
wherein t is a normalized value, t' is a Gamma-processed value, P z Is an inverse normalized value;
standard deviation stretching includes: calculating new pixel maximum and minimum values by using the standard deviation and the average value of the pixels of the image, and finally performing maximum and minimum stretching treatment, wherein the calculation formula is as follows:
Figure BDA0004139756280000103
wherein the maximum value of the pixel value of the optical remote sensing image data is P m A pixelMinimum value of P n The average value of the pixel values is P p The mean square error of the pixel value is P σ The pixel has no data value P o The histogram level is P L N is the total number of pixels, and the histogram frequency is [ N ] 1 ，N 2 ......N L ]Gama value of G a The percentage cut-off value is S p Standard deviation ratio S s The method comprises the steps of carrying out a first treatment on the surface of the Let the original pixel value of the input optical remote sensing image data be P s The output tile pixel value is P d 。
As shown in fig. 3 and 4, the remote sensing image real-time rendering method based on the B/S architecture of the present invention includes the following steps:
(1) Requesting tile data by a foreground;
(2) Judging the legality of rendering parameters, wherein the rendering parameters comprise the range of PNG tile data, an optical remote sensing image data path, a Gamma value, a percentage cut-off value and a standard deviation value;
(3) If yes, opening the optical remote sensing image data to acquire the information of the width, height, wave band number and data type of the optical remote sensing image; if not, directly ending;
(4) Reading the data attribute of the optical remote sensing image, acquiring coordinate reference information (WKT, affine transformation six parameters), and shortening the validity of the coordinate reference information, if so, using the coordinate projection conversion geometric relationship between the optical remote sensing image data and the tile data; if not, acquiring RPC parameters of the optical remote sensing image;
(5) If the RPC parameter is invalid, acquiring the RPC parameter of the optical remote sensing image, and judging the validity of the RPC parameter; if so, using an RPC conversion geometric relation between the optical remote sensing image data and the tile data;
(6) If so, using an RPC conversion geometric relation between the optical remote sensing image data and the tile data; if not, directly ending;
(7) Assigning a memory file attribute, wherein the memory file attribute comprises a width, a height, a wave band number, a data type and a wave band index value;
(8) Calculating affine transformation six parameters of PNG tile data according to the PNG tile data range and the tile size;
(9) Calculating four corner geographic coordinates of the optical remote sensing image data according to the width, height and geometric conversion relation of the optical remote sensing image data;
(10) Calculating the intersection of the geographical range of the optical remote sensing image data and the geographical range of the PNG tile data;
(11) Reversely calculating a start row number index value and an end row number index value of the optical remote sensing image data according to the intersection geographic range and the geometric conversion relation; similarly, calculating a start line index value and an end line index value of the intersection geographic range on the PNG memory file according to the intersection geographic range and affine transformation six parameters of the PNG tile data, and finally calculating the number of lines and the number of columns of the intersection geographic range on the PNG memory file;
(12) Calculating the width and height of an optical remote sensing image data block according to the number of rows and columns of the step intersection geographic range on the PNG memory file;
(13) Back-calculating affine transformation six parameters or RPC parameters of the optical remote sensing image data block according to the conversion geometric relationship;
(14) Reading an optical remote sensing image data block;
(15) Taking the PNG memory file as a processing object, and circularly processing each pixel value in the PNG memory file from left to right and from top to bottom; calculating the geographic coordinates of pixels on the PNG memory file according to the line number, the column number and the conversion geometric relation of the intersection geographic range on the PNG memory file;
(13) If the optical remote sensing image data block is in a coordinate projection relationship, converting the geographic coordinates into projection coordinates according to the conversion geometric relationship, and converting the projection coordinates into pixel coordinates of the optical remote sensing image data block;
(17) If the optical remote sensing image data block is in the RPC relationship, converting the geographic coordinate into the pixel coordinate of the optical remote sensing image data block according to the conversion geometric relationship;
(18) Reading a specified pixel coordinate pixel value from an optical remote sensing image data block;
(19) Carrying out data value-free processing on the pixel values, and carrying out one mode of maximum and minimum value stretching, percentage truncation, histogram equalization, gama enhancement or standard deviation stretching on the data which are subjected to data value-free processing;
(20) Sequentially cycling all pixel data in the PNG memory file;
(21) And constructing PNG memory data by using the output tile data, setting a data-free value of the memory file, and returning the generated tile data to the foreground.
According to an aspect of the present invention, there is provided an electronic apparatus including: one or more processors, one or more memories, and one or more computer programs; the processor is connected with the memory, the one or more computer programs are stored in the memory, and when the electronic device runs, the processor executes the one or more computer programs stored in the memory, so that the electronic device executes a remote sensing image real-time rendering method based on the B/S architecture according to any one of the technical schemes.
According to one aspect of the present invention, a computer readable storage medium is provided, configured to store computer instructions, where the computer instructions, when executed by a processor, implement a B/S architecture-based remote sensing image real-time rendering method according to any one of the above technical solutions.
According to the remote sensing image real-time rendering method, equipment and storage medium based on the B/S architecture, according to the tile display rule of the B/S architecture, optical remote sensing image data are rendered into PNG tile data with the size of 256 pixels in real time in the background, and the rendered PNG tiles are transmitted to a foreground view for display according to requirements. In order to realize real-time generation of the optical remote sensing image tile data, a geometric relation between the optical remote sensing image data and the PNG tile data is established in the background, an optical remote sensing image data block range to be read is calculated according to the PNG tile data range and the geometric relation, the optical remote sensing image data block is read and processed, pixel-level geometric relation conversion is realized, PNG tile data is constructed and generated, and finally an algorithm scheme meeting B/S architecture tile data is realized for real-time generation of the optical remote sensing image data.
Furthermore, it should be noted that the present invention can be provided as a method, an apparatus, or a computer program product. Accordingly, embodiments of the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, embodiments of the invention may take the form of a computer program product on one or more computer-usable storage media having computer-usable program code embodied therein.
Embodiments of the present invention are described with reference to flowchart illustrations and/or block diagrams of methods, terminal devices (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, embedded processor, or other programmable data processing terminal device to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing terminal device, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks. These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It should also be noted that, in this document, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or terminal that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or terminal. Without further limitation, an element defined by the phrase "comprising one does not exclude the presence of other like elements in a process, method, article or terminal device comprising the element.
It is finally pointed out that the above description of the preferred embodiments of the invention, it being understood that although preferred embodiments of the invention have been described, it will be obvious to those skilled in the art that, once the basic inventive concepts of the invention are known, several modifications and adaptations can be made without departing from the principles of the invention, and these modifications and adaptations are intended to be within the scope of the invention. It is therefore intended that the following claims be interpreted as including the preferred embodiment and all such alterations and modifications as fall within the scope of the embodiments of the invention.

Claims (10)

1. A remote sensing image real-time rendering method based on a B/S architecture is characterized by comprising the following steps:
step S1, requesting tile data by a foreground;
s2, reading attribute information of the optical remote sensing image data;
s3, establishing a conversion geometric relationship between the optical remote sensing image data and the PNG tile data;
s4, calculating a range intersection of the optical remote sensing image data and the PNG tile data;
s5, performing data processing on the read optical remote sensing image data block data;
and S6, constructing PNG tile data, and returning the generated tile data to the foreground.
2. The method according to claim 1, wherein in the step S1, the foreground requests the required tile data in the Http request manner, and the request information includes the optical informationRemote sensing image data path and tile data longitude and latitude range E w 。
3. The B/S architecture-based remote sensing image real-time rendering method according to claim 1, wherein in the step S2, the optical remote sensing image data is opened, and the width W of the optical remote sensing image data is read s High H s Band number B s Geographic reference information, RPC parameters.
4. The B/S architecture-based remote sensing image real-time rendering method according to claim 1, wherein in the step S3, if the optical remote sensing image data has the geographic reference information, a coordinate projection conversion geometric relationship between the optical remote sensing image data and the tile data is established; if the RPC parameters exist in the optical remote sensing image data, establishing an RPC conversion geometric relationship between the optical remote sensing image data and the tile data;
establishing a coordinate projection conversion geometric relation: setting a pixel coordinate (i) of a certain point A of the optical remote sensing image data s ,j s ) Corresponding projection coordinates (X s ,Y s ) Corresponding geographic coordinates (L s ,B s ) Affine transformation six parameters (a s0 ,a s1 ,a s2 ,b s0 ,b s1 ,b s2 ) The coordinate reference is in WKT form, and the conversion relationship between the pixel coordinate and the projection coordinate is as follows:
Figure FDA0004139756270000011
the conversion relation between the projection coordinates and the geographic coordinates can be created by using the OGRCoordinates conversion information in the GDAL in a coordinate reference WKT form, and the conversion of the projection coordinates and the geographic coordinates is performed by using a function conversion form;
establishing an RPC conversion geometric relation: setting a pixel coordinate (i) of a certain point A of the optical remote sensing image data s ,j s ) Corresponding projection coordinates (X s ,Y s ) Corresponding geographic coordinates(L s ,B s ) The conversion relation between the pixel coordinates and the geographic coordinates uses GDALCREATE RPCTransformer in GDAL to construct an RPC parameter conversion model, and the conversion between the pixel coordinates and the geographic coordinates is carried out through a function GDALRPCTransform.
5. The B/S architecture-based remote sensing image real-time rendering method according to claim 4, wherein in the step S4, specifically comprising:
calculating the geographic range E of the optical remote sensing image data through the conversion geometric relationship established in the step S3 s Intersection E of tile data and optical remote sensing image data geographic range c The calculation formula of (2) is as follows:
E c ＝E w ∩E s
wherein E is w Represented as tile data latitude and longitude ranges.
6. The B/S architecture-based remote sensing image real-time rendering method according to claim 5, wherein in step S5, specifically comprising:
step S51, back-calculating the intersection E of the geographic ranges according to the conversion geometric relationship established in the step S3 c Pixel coordinates on the optical remote sensing image;
step S52, reading optical remote sensing image block data according to the pixel coordinate range on the optical remote sensing image;
step S53, performing data value-free processing on the read optical remote sensing image block data;
and step S54, performing one mode of maximum and minimum stretching, percentage truncation, histogram equalization, gama enhancement or standard deviation stretching on the data processed without the data value.
7. The B/S architecture-based remote sensing image real-time rendering method according to claim 6, wherein in the step S53, the data value processing includes: if the optical remote sensing image data does not have a data-free value, the pixel value 0 is regarded as the data-free value; if no data value exists, the output pixel value is assigned to 0, and the output tile pixel value formula is as follows:
Figure FDA0004139756270000021
wherein, the original pixel value of the input optical remote sensing image data is P s The output tile pixel value is P d 。
8. The method of claim 7, wherein in the step S54,
the maximum minimum stretching includes: maximum and minimum stretching and 16-bit to 8-bit processing, wherein the formula is as follows:
Figure FDA0004139756270000031
the percentage cutoff includes: s is intercepted by the left end and the right end of the histogram frequency N respectively by using a percentage cut-off value p Calculate a new pixel maximum P' m And pixel minimum value P' n The output tile pixel value is calculated using the following formula:
Figure FDA0004139756270000032
the histogram equalization includes: calculating the gray cumulative distribution frequency S of each level of the histogram according to the histogram level and the histogram frequency L The calculation formula is as follows:
Figure FDA0004139756270000033
calculating histogram equalized pixel value P using gray scale cumulative distribution frequency z Finally, using the maximum and minimumThe value stretching is processed, and the calculation formula is as follows:
P z ＝S L *(L-1)+P n
Figure FDA0004139756270000034
the Gama enhancement includes: the original pixel value is normalized, gamma processing and inverse normalization processing are carried out, and finally maximum and minimum stretching processing is carried out, wherein the calculation formula is as follows:
Figure FDA0004139756270000035
wherein t is a normalized value, t' is a Gamma-processed value, P z Is an inverse normalized value;
the standard deviation stretching includes: calculating new pixel maximum and minimum values by using the standard deviation and the average value of the pixels of the image, and finally performing maximum and minimum stretching treatment, wherein the calculation formula is as follows:
Figure FDA0004139756270000041
wherein the maximum value of the pixel value of the optical remote sensing image data is P m The minimum value of the pixel value is P n The average value of the pixel values is P p The mean square error of the pixel value is P σ The pixel has no data value P o The histogram level is P L N is the total number of pixels, and the histogram frequency is [ N ] 1 ,N 2 ……N L ]Gama value of G a The percentage cut-off value is S p Standard deviation ratio S s The method comprises the steps of carrying out a first treatment on the surface of the Let the original pixel value of the input optical remote sensing image data be P s The output tile pixel value is P d 。
9. An electronic device, comprising: one or more processors, one or more memories, and one or more computer programs; wherein the processor is connected to the memory, the one or more computer programs are stored in the memory, and when the electronic device is running, the processor executes the one or more computer programs stored in the memory, so that the electronic device executes the remote sensing image real-time rendering method based on the B/S architecture as claimed in any one of claims 1 to 8.
10. A computer readable storage medium storing computer instructions which, when executed by a processor, implement the B/S architecture-based remote sensing image real-time rendering method of any one of claims 1 to 8.
CN202310285692.8A 2023-03-22 2023-03-22 Remote sensing image real-time rendering method, device and storage medium based on B/S architecture Pending CN116309915A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310285692.8A CN116309915A (en) 2023-03-22 2023-03-22 Remote sensing image real-time rendering method, device and storage medium based on B/S architecture

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310285692.8A CN116309915A (en) 2023-03-22 2023-03-22 Remote sensing image real-time rendering method, device and storage medium based on B/S architecture

Publications (1)

Publication Number Publication Date
CN116309915A true CN116309915A (en) 2023-06-23

Family

ID=86788394

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310285692.8A Pending CN116309915A (en) 2023-03-22 2023-03-22 Remote sensing image real-time rendering method, device and storage medium based on B/S architecture

Country Status (1)

Country Link
CN (1) CN116309915A (en)

Similar Documents

Publication Publication Date Title
US8411954B2 (en) Image processing system, image processing program, and image processing method
US10402941B2 (en) Guided image upsampling using bitmap tracing
CN110917617B (en) Method, device, equipment and storage medium for generating water ripple image
CN111369657A (en) Three-dimensional thermodynamic diagram generation method and device, computer equipment and storage medium
CN108242063B (en) Light field image depth estimation method based on GPU acceleration
US11205290B2 (en) Method and device for inserting an image into a determined region of a target eye image
CN111476096A (en) Face image key point method and device, computer equipment and storage medium
CN115082322B (en) Image processing method and device, and training method and device of image reconstruction model
CN110443252A (en) A kind of character detecting method, device and equipment
CN113222921A (en) Image processing method and system
CN116109765A (en) Three-dimensional rendering method and device for labeling objects, computer equipment and storage medium
CN112651399B (en) Method for detecting same-line characters in inclined image and related equipment thereof
CN117611703A (en) Barrage character rendering method, barrage character rendering device, barrage character rendering equipment, storage medium and program product
CN112927163A (en) Image data enhancement method and device, electronic equipment and storage medium
CN113538502A (en) Picture clipping method and device, electronic equipment and storage medium
CN116309915A (en) Remote sensing image real-time rendering method, device and storage medium based on B/S architecture
CN112215308A (en) Single-order detection method and device for hoisted object, electronic equipment and storage medium
US11450005B2 (en) Image processing method, apparatus, terminal device and storage medium
US20220318950A1 (en) Video enhancement method and apparatus, and electronic device and storage medium
CN114996374A (en) Online data visualization implementation method, system, device and medium
CN113034368A (en) Image super-resolution model training method and device, computer equipment and medium
CN117314756B (en) Verification and protection method and device based on remote sensing image, computer equipment and storage medium
CN115526775B (en) Image data processing method and device
CN110942500B (en) Method and device for converting static diagram into dynamic diagram
CN114419322B (en) Image instance segmentation method and device, electronic equipment and storage medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination