CN112598801A - Three-dimensional visualization method and device for environmental data - Google Patents

Three-dimensional visualization method and device for environmental data Download PDF

Info

Publication number
CN112598801A
CN112598801A CN202011508725.3A CN202011508725A CN112598801A CN 112598801 A CN112598801 A CN 112598801A CN 202011508725 A CN202011508725 A CN 202011508725A CN 112598801 A CN112598801 A CN 112598801A
Authority
CN
China
Prior art keywords
grid
data
grid point
point
dimensional
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
CN202011508725.3A
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.)
3Clear Technology Co Ltd
Original Assignee
3Clear 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 3Clear Technology Co Ltd filed Critical 3Clear Technology Co Ltd
Priority to CN202011508725.3A priority Critical patent/CN112598801A/en
Publication of CN112598801A publication Critical patent/CN112598801A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/10Geometric effects
    • G06T15/20Perspective computation
    • G06T15/205Image-based rendering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/50Lighting effects
    • G06T15/80Shading
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05Geographic models

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Computer Graphics (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Remote Sensing (AREA)
  • Computing Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a three-dimensional visualization method and a three-dimensional visualization device for environmental data, wherein a three-dimensional space coordinate system corresponding to a designated area is established; constructing a three-dimensional grid model according to a NAQPMS system based on a three-dimensional space coordinate system; wherein, the range of coordinate axes of the three-dimensional grid model and the grid size correspond to the NAQPMS system; acquiring environmental data of each grid in the three-dimensional grid model through a NAQPMS system; acquiring first environment data of a first designated grid at a first moment and second environment data of the first designated grid at a second moment; acquiring environmental data of a third moment according to the first environmental data and the second environmental data through the following formula; the third moment is a moment between the first moment and the second moment;
Figure DDA0002845693350000011
Figure DDA0002845693350000012
the method solves the problems that pollutant concentration distribution cannot be embodied through three-dimensional visualization in the prior art, and air quality prediction data of any grid point and any moment cannot be acquired, and realizes three-dimensional visualization of air quality mode prediction data at any moment.

Description

Three-dimensional visualization method and device for environmental data
Technical Field
The invention relates to the technical field of drawing and environmental protection, in particular to a three-dimensional visualization method and device for environmental data.
Background
With the improvement of the quality of life of people, people pay more and more attention to the environmental problem, and the environmental air quality forecast is gradually the same as the weather forecast, so that the environmental air quality forecast becomes the information concerned by people.
The current air quality forecast data is usually displayed on a two-dimensional map, that is, conventional terrain visualization technologies (such as contour topographic maps, section maps, scenic drawings and the like) are based on two dimensions, and the pollutant concentration forecast data is displayed in a point and surface form.
The two-dimensional display mode cannot reflect the three-dimensional distribution characteristics of pollutant concentration forecast, however, the atmospheric pollution process occurs in a three-dimensional space, so that the three-dimensional overall characteristics of pollutant concentration distribution can be known, and the improvement of the environmental air quality is greatly facilitated.
Aiming at the problems that the existing air quality forecast data is usually displayed on a two-dimensional map, the pollutant concentration distribution cannot be embodied through three-dimensional visualization, and the air quality forecast data of any grid point and any moment cannot be acquired in the prior art, an effective solution is not provided.
Disclosure of Invention
In view of this, embodiments of the present invention provide a method and an apparatus for three-dimensional visualization of environmental data, so as to solve the problems that air quality prediction data is usually displayed on a two-dimensional map, pollutant concentration distribution cannot be reflected through three-dimensional visualization, and air quality prediction data at any grid point and at any time cannot be acquired.
Therefore, the embodiment of the invention provides the following technical scheme:
in a first aspect of the present invention, a method for three-dimensional visualization of environmental data is provided, including:
constructing a three-dimensional space coordinate system corresponding to the designated area;
constructing a three-dimensional grid model according to a NAQPMS system based on the three-dimensional space coordinate system; wherein the range of coordinate axes and the mesh size of the three-dimensional mesh model correspond to the NAQPMS system;
acquiring environmental data of each grid in the three-dimensional grid model through the NAQPMS system;
acquiring first environment data of a first designated grid at a first moment and second environment data of the first designated grid at a second moment;
acquiring environmental data of a third moment according to the first environmental data and the second environmental data through the following formula; the third moment is a moment between the first moment and the second moment;
Figure BDA0002845693330000021
a represents the first time, B represents the second time, T represents the third time, Tn represents the environment data at the third time, An represents the first environment data, and Bn represents the second environment data;
and displaying the environment data of each grid through the three-dimensional grid model.
Optionally, the method further comprises:
acquiring a designated three-dimensional grid and a designated grid point P in the designated three-dimensional grid; wherein, the four grid points of one face of the specified three-dimensional grid are A, B, C, D, and the four grid points of the corresponding face of the one face are E, F, G, H;
acquiring environment data of a grid point A, a grid point B, a grid point C, a grid point D, a grid point E, a grid point F, a grid point G and a grid point H through the NAQPMS system;
the environment data of the specified grid point P is calculated by the following formula:
P.v=A.v*Aq+B.v*Bq+C.v*Cq+D.v*Dq+E.v*Eq+F.v*Fq+G.v*Gq+H.v*Hq;
where P.v denotes the environmental data of the specified grid point P, A.v denotes the environmental data of grid point A, B.v denotes the environmental data of grid point B, C.v denotes the environmental data of grid point C, D.v denotes the environmental data of grid point D, E.v denotes the environmental data of grid point E, F.v denotes the environmental data of grid point F, G.v denotes the environmental data of grid point G, H.v denotes the environmental data of grid point H, Aq denotes the pollution weight of grid point A to the specified grid point P, Bq denotes the pollution weight of grid point B to the specified grid point P, Cq denotes the pollution weight of grid point C to the specified grid point P, Dq denotes the pollution weight of grid point D to the specified grid point P, Eq denotes the pollution weight of grid point E to the specified grid point P, Fq denotes the pollution weight of grid point F to the specified grid point P, Gq denotes the pollution weight of grid point G to the specified grid point P, hq represents the pollution weight of grid point H to the designated grid point P;
wherein the content of the first and second substances,
Aq=(1–Qx)*(1–Qy)*(1–Qz);
Bq=Qx*(1–Qy)*(1–Qz);
Cq=Qx*Qy*(1–Qz);
Dq=(1–Qx)*Qy*(1–Qz);
Eq=(1–Qx)*(1–Qy)*Qz;
Fq=Qx*(1–Qy)*Qz;
Gq=Qx*Qy*Qz;
Hq=(1–Qx)*Qy*Qz;
wherein the content of the first and second substances,
Qx=(x–A.x)/Lx;
Qy=(y–A.y)/Ly;
Qz=(z–A.z)/Lz;
wherein the content of the first and second substances,
Lx=G.x–A.x;
Ly=G.y–A.y;
Lz=G.z–A.z;
lx, Ly and Lz are respectively the length, width and height of the designated three-dimensional grid, A.x, A.y and A.z are respectively the X-axis coordinate of the point A, the Y-axis coordinate of the point A and the Z-axis coordinate of the point A, G.x, G.y and G.z are respectively the X-axis coordinate of the point G, the Y-axis coordinate of the point G and the Z-axis coordinate of the point G, and X, Y and Z are respectively the X, Y and Z coordinate values of the designated grid point P in the three-dimensional coordinates.
Optionally, the method further comprises:
acquiring a first data display area parallel to the ground and a second data display area correspondingly perpendicular to the first data display area according to the longitude and latitude coordinates of the designated area;
correspondingly drawing the first data display area and the second data display area to a virtual earth;
acquiring environment data of each grid in the first data presentation area and the second data presentation area through the NAQPMS system;
and displaying the environment data of each grid in the first data display area and the second data display area in the virtual earth.
Optionally, the method further comprises:
obtaining environmental data for the second designated grid by:
Mvalue=Tvalue*QT+Bvalue*QB;
QT=(dataHgt-BottomHgt)/(TopHgt–BottomHgt);
QB=(TopHgt-dataHgt)/(TopHgt–BottomHgt);
wherein Mvalue represents environment data of the second specified grid, Tvalue represents environment data of a first grid, QT represents concentration contribution factor of the first grid, Bvalue represents environment data of a second grid, QB represents concentration contribution factor of the second grid, dataHgt represents height of the second specified grid, TopHgt represents height of the first grid, and BottomHgt represents height of the second grid; wherein the layer height of the second designated mesh is between the layer heights of the first mesh and the second mesh, and the layer height of the first mesh is higher than the layer height of the second mesh.
Optionally, the displaying the environmental data of each mesh by the three-dimensional mesh model includes:
acquiring color indexes corresponding to the environment data of each grid;
acquiring colors corresponding to the environment data of each grid according to the color indexes;
wherein the color index corresponding to the environment data of each grid is obtained by the following formula:
cIndex=((V–Nmin)/(Nmax–Nmin))*CLength;
wherein, cndex represents a color index, V represents environment data corresponding to each grid, Nmin represents a minimum value of the environment data, Nmax represents a maximum value of the environment data, and CLength represents the number of colors of the rendering legend of the environmental pollution.
In a second aspect of the present invention, there is provided an apparatus for three-dimensional visualization of environmental data, comprising:
the first construction module is used for constructing a three-dimensional space coordinate system corresponding to the designated area;
the second construction module is used for constructing a three-dimensional grid model based on the three-dimensional space coordinate system according to a NAQPMS system; wherein the range of coordinate axes and the mesh size of the three-dimensional mesh model correspond to the NAQPMS system;
a first obtaining module, configured to obtain, by the NAQPMS system, environment data of each grid in the three-dimensional grid model;
the second acquisition module is used for acquiring first environmental data of the first designated grid at a first moment and second environmental data of the first designated grid at a second moment;
the third obtaining module is used for obtaining the environmental data of a third moment according to the first environmental data and the second environmental data through the following formula; the third moment is a moment between the first moment and the second moment;
Figure BDA0002845693330000051
a represents the first time, B represents the second time, T represents the third time, Tn represents the environment data at the third time, An represents the first environment data, and Bn represents the second environment data;
and the first display module is used for displaying the environment data of each grid through the three-dimensional grid model.
Optionally, the apparatus further comprises:
the fourth acquisition module is used for acquiring a specified three-dimensional grid and a specified grid point P in the specified three-dimensional grid; wherein, the four grid points of one face of the specified three-dimensional grid are A, B, C, D, and the four grid points of the corresponding face of the one face are E, F, G, H;
a fifth obtaining module, configured to obtain, by the NAQPMS system, environment data of the grid point a, the grid point B, the grid point C, the grid point D, the grid point E, the grid point F, the grid point G, and the grid point H;
a calculating module, configured to calculate the environment data of the specified grid point P by using the following formula:
P.v=A.v*Aq+B.v*Bq+C.v*Cq+D.v*Dq+E.v*Eq+F.v*Fq+G.v*Gq+H.v*Hq;
where P.v denotes the environmental data of the specified grid point P, A.v denotes the environmental data of grid point A, B.v denotes the environmental data of grid point B, C.v denotes the environmental data of grid point C, D.v denotes the environmental data of grid point D, E.v denotes the environmental data of grid point E, F.v denotes the environmental data of grid point F, G.v denotes the environmental data of grid point G, H.v denotes the environmental data of grid point H, Aq denotes the pollution weight of grid point A to the specified grid point P, Bq denotes the pollution weight of grid point B to the specified grid point P, Cq denotes the pollution weight of grid point C to the specified grid point P, Dq denotes the pollution weight of grid point D to the specified grid point P, Eq denotes the pollution weight of grid point E to the specified grid point P, Fq denotes the pollution weight of grid point F to the specified grid point P, Gq denotes the pollution weight of grid point G to the specified grid point P, hq represents the pollution weight of grid point H to the designated grid point P;
wherein the content of the first and second substances,
Aq=(1–Qx)*(1–Qy)*(1–Qz);
Bq=Qx*(1–Qy)*(1–Qz);
Cq=Qx*Qy*(1–Qz);
Dq=(1–Qx)*Qy*(1–Qz);
Eq=(1–Qx)*(1–Qy)*Qz;
Fq=Qx*(1–Qy)*Qz;
Gq=Qx*Qy*Qz;
Hq=(1–Qx)*Qy*Qz;
wherein the content of the first and second substances,
Qx=(x–A.x)/Lx;
Qy=(y–A.y)/Ly;
Qz=(z–A.z)/Lz;
wherein the content of the first and second substances,
Lx=G.x–A.x;
Ly=G.y–A.y;
Lz=G.z–A.z;
lx, Ly and Lz are respectively the length, width and height of the designated three-dimensional grid, A.x, A.y and A.z are respectively the X-axis coordinate of the point A, the Y-axis coordinate of the point A and the Z-axis coordinate of the point A, G.x, G.y and G.z are respectively the X-axis coordinate of the point G, the Y-axis coordinate of the point G and the Z-axis coordinate of the point G, and X, Y and Z are respectively the X, Y and Z coordinate values of the designated grid point P in the three-dimensional coordinates.
Optionally, the apparatus further comprises:
the sixth acquisition module is used for acquiring a first data display area parallel to the ground and a second data display area which is correspondingly vertical to the first data display area according to the longitude and latitude coordinates of the designated area;
the drawing module is used for correspondingly drawing the first data display area and the second data display area to a virtual earth;
a seventh obtaining module, configured to obtain, by the NAQPMS system, environment data of each grid in the first data presentation area and the second data presentation area;
and the second display module is used for displaying the environment data of each grid in the first data display area and the second data display area in the virtual earth.
In a third aspect of the present invention, there is provided an electronic device comprising: at least one processor; and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the at least one processor to cause the at least one processor to perform the method for three-dimensional visualization of environmental data as set forth in any one of the above first aspects.
In a fourth aspect of the present invention, a computer-readable storage medium is provided, on which computer instructions are stored, which when executed by a processor, implement the method for three-dimensional visualization of environmental data according to any one of the first aspect above.
The technical scheme of the embodiment of the invention has the following advantages:
the embodiment of the invention provides a three-dimensional visualization method and a three-dimensional visualization device for environmental data, wherein the method comprises the following steps: constructing a three-dimensional space coordinate system corresponding to the designated area; constructing a three-dimensional grid model according to a NAQPMS system based on a three-dimensional space coordinate system; wherein, the range of coordinate axes of the three-dimensional grid model and the grid size correspond to the NAQPMS system; acquiring environmental data of each grid in the three-dimensional grid model through a NAQPMS system; acquiring first environment data of a first designated grid at a first moment and second environment data of the first designated grid at a second moment; acquiring environmental data of a third moment according to the first environmental data and the second environmental data through the following formula; the third moment is a moment between the first moment and the second moment;
Figure BDA0002845693330000091
a represents the first time, B represents the second time, T represents the third time, Tn represents the environment data at the third time, An represents the first environment data, and Bn represents the second environment data; and displaying the environment data of each grid through the three-dimensional grid model. The method solves the problems that the air quality forecast data is usually displayed on a two-dimensional map, the pollutant concentration distribution cannot be reflected through three-dimensional visualization, and the air quality forecast data at any grid point and any moment cannot be acquired in the prior art, realizes the three-dimensional visualization of the air quality mode forecast data at any moment, and can more accurately and intuitively display the distribution of the environmental data in the geographic space.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flow chart of a method for three-dimensional visualization of environmental data according to an embodiment of the invention;
FIG. 2 is a schematic diagram of a three-dimensional data grid according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a three-dimensional data grid employing tri-linear interpolation in accordance with an embodiment of the present invention;
FIG. 4 is a display diagram of a base path in a virtual earth, according to an embodiment of the invention;
FIG. 5 is a path splitting build vertical show surface according to an embodiment of the present invention;
FIG. 6 is a schematic diagram illustrating a transition of ground layer data from two dimensions to three dimensions in accordance with an embodiment of the invention;
FIG. 7 is a schematic diagram of a multidimensional visualization of pollution situation based on a three-dimensional virtual earth according to an embodiment of the invention;
FIG. 8 is a schematic diagram of a spatial interpolation structure for elevation data according to an embodiment of the present invention;
FIG. 9 is a PM2.5 rendering illustration according to an embodiment of the invention;
fig. 10 is a block diagram of a three-dimensional visualization apparatus for environmental data according to an embodiment of the present invention;
fig. 11 is a schematic diagram of a hardware structure of an electronic device according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
In the description of the present application, it is to be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", and the like indicate orientations or positional relationships based on those shown in the drawings, and are used merely for convenience of description and for simplicity of description, and do not indicate or imply that the referenced device or element must have a particular orientation, be constructed in a particular orientation, and be operated, and thus should not be considered as limiting the present application. Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more features. In the description of the present application, "a plurality" means two or more unless specifically limited otherwise.
In this application, the word "exemplary" is used to mean "serving as an example, instance, or illustration. Any embodiment described herein as "exemplary" is not necessarily to be construed as preferred or advantageous over other embodiments. The following description is presented to enable any person skilled in the art to make and use the application. In the following description, details are set forth for the purpose of explanation. It will be apparent to one of ordinary skill in the art that the present application may be practiced without these specific details. In other instances, well-known structures and processes are not set forth in detail in order to avoid obscuring the description of the present application with unnecessary detail. Thus, the present application is not intended to be limited to the embodiments shown, but is to be accorded the widest scope consistent with the principles and features disclosed herein.
In accordance with an embodiment of the present invention, there is provided an embodiment of a method for three-dimensional visualization of environmental data, it being noted that the steps illustrated in the flowchart of the figure may be performed in a computer system such as a set of computer-executable instructions and that, although a logical order is illustrated in the flowchart, in some cases, the steps illustrated or described may be performed in an order different than that illustrated herein.
In addition, the technical features involved in the different embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
In this embodiment, a three-dimensional visualization method for environmental data is provided, which may be used in an environmental monitoring system, an environmental forecast early warning system, and the like, fig. 1 is a flowchart of the three-dimensional visualization method for environmental data according to an embodiment of the present invention, as shown in fig. 1, the flowchart includes the following steps:
step S101, a three-dimensional space coordinate system corresponding to the designated area is constructed. For example, the three-dimensional space coordinate system has an axis parallel to the longitudinal direction as an X-axis, an axis parallel to the latitudinal direction as a Y-axis, and an elevation axis as a Z-axis.
Step S102, constructing a three-dimensional grid model according to a Nested grid Air Quality Prediction mode System (NAQPMS for short) System based on a three-dimensional space coordinate System; wherein the range of coordinate axes and the mesh size of the three-dimensional mesh model correspond to the NAQPMS system. The NAQPMS mode system is widely applied to sulfur oxide transnational transportation, sand dust sand raising, transportation and settlement simulation, research on the influence of acid rain on the environment, ozone simulation and air quality simulation research on the scales of cities, regions and the like, and successfully realizes wide application in business. The NAQPMS is independently developed and developed by the atmospheric physics research institute of Chinese academy of sciences. The model system has gone through the development of recent 20 years and is developed by integrating a series of city and region scale air quality models which are independently developed. The mode can be used for researching the air pollution problem of regional scales, researching the occurrence mechanism and the change rule of the problems of air quality and the like of urban scales, and researching the mutual influence process among different scales. The mode is an important tool for researching the interaction among pollutant discharge amount, meteorological conditions, chemical conversion and dry-wet removal, and can provide scientific pollution discharge control strategies for an environmental decision part. The air quality prediction subsystem (NAQPM) is the core of the whole model system and mainly treats the physical and chemical processes of pollutant emission, advection transportation, diffusion, dry and wet sedimentation, gas phase, liquid phase and heterogeneous reaction and the like. The spatial structure of the device is a three-dimensional Euler conveying mode, and the vertical coordinate adopts a terrain following coordinate. The horizontal structure is a multi-nested grid, a unidirectional and bidirectional nesting technology is adopted, the resolution is 3-81 km, and the vertical unequally-spaced structure is divided into 20 layers.
And step S103, acquiring environment data of each grid in the three-dimensional grid model through the NAQPMS system. The environmental data may be, for example, pollutant concentration data, or any other data that may represent environmental information. The contaminant distribution data used in this alternative embodiment may be generated by the NAQPMS mode. Specifically, the pollutant forecast data generated by the NAQPMS mode is a two-dimensional grid structure with one layer parallel to the ground, and the data of each layer represents the concentration grid distribution of pollutants of different height layers.
Specifically, the pattern prediction data is stored on the Web server, and requests all layer heights of the pattern prediction area to be visualized as required, for example, data of 20 layers of pattern prediction PM2.5 is used, and the sequence number of the layer height is pollutant concentration data of 1-20. According to the returned mode forecast data and mode basic information (the geographic spatial range, the number of layers, the average height represented by the layers, the pollutant indexes and the like), the number of grids and the size of each grid are determined in the directions of longitude, latitude and height, the position of each grid point and the pollutant concentration data on the grid points are set, a three-dimensional data grid of the pollutant concentration is constructed, and as shown in fig. 2, the organized three-dimensional data grid is stored in a computer memory.
Step S104, acquiring first environment data of the first designated grid at a first time and second environment data of the first designated grid at a second time.
Step S105, acquiring environmental data of a third moment according to the first environmental data and the second environmental data through the following formula; the third moment is a moment between the first moment and the second moment;
Figure BDA0002845693330000131
a denotes a first time, B denotes a second time, T denotes a third time, Tn denotes environment data of the third time, An denotes first environment data, and Bn denotes second environment data.
Because the pollutant forecast data has three dimensions, large time and space span, and many kinds of data (the forecast data of each pollutant is a kind of data, such as PM2.5, PM10, dust and sand, etc.), the amount of finally generated data is very large, and if the air quality distribution is continuously and dynamically shown according to the normal data request mode and the display processing method, the effect of continuously and dynamically showing the air quality distribution is not achieved. For example, the geographic range of the three-dimensional data grid constructed in the embodiment of the present invention is: the longitude is between 26.03304 degrees and 178.0110065 degrees, the latitude is between-6.59913 degrees (south latitude) and 68.11035717 degrees, and the layer height is 1 to 20 layers; there are 750 grids in the longitudinal direction, one data grid every 0.20270270 degrees; there are 371 grids in the latitudinal direction, one data grid every 0.20270270 degrees. Therefore, the grid point total number of the density forecast data of PM2.5 at one time is calculated as follows: sum Xn Yn Zn; wherein Xn, Yn and Zn represent the data grid points in the longitude direction, the latitude direction and the altitude direction respectively. Therefore, at one point in time in this embodiment, there are 750 × 371 × 20 — 5565000 grid points for a single contamination indicator, each grid point storing contaminant concentration data. In this way, different contaminants are displayed, and the amount of data may be large at different successive times.
Simulating middle time data by using 2 time data in consideration of data request efficiency and continuous visualization effect of a Web browserData transmission quantity can be reduced, and the dynamic effect of the pollutants changing according to the trend can be improved. For example, the pollutant forecast concentration data at any time of the whole day can be obtained by interpolation in the time dimension by using the pollutant forecast concentration data at 0 point, 4 points, 8 points, 12 points, 16 points, 20 points and 6 times of each day. By this step, the environmental data of each grid at any time can be obtained. For example, the predicted PM2.5 concentration data at 8 and 12 points are known to be 36 μ g/m3And 68. mu.g/m3Then the concentration data for PM2.5 at 11 points can be modeled as:
Figure BDA0002845693330000141
Figure BDA0002845693330000142
therefore, the data transmission pressure can be reduced, resources can be saved for computer processing data, the processed data amount is small, and continuous and dynamic display of the trend change of the pollutant concentration in the geographic space along with the time is realized.
And step S106, displaying the environment data of each grid through the three-dimensional grid model.
Through the steps, the NAQPMS system is used for obtaining pollutant forecast data, the pollutant forecast data are displayed in each grid of the corresponding three-dimensional grid model, and the pollutant concentration of each grid at any moment can be obtained through a simple algorithm, so that the problems that in the prior art, pollutant concentration display is usually performed on a two-dimensional map, pollutant concentration distribution cannot be reflected through three-dimensional visualization, and pollutant concentrations of any grid point and any moment cannot be obtained are solved, the three-dimensional visualization of the pollutant concentration at any moment is realized, and the distribution of environment data in a geographic space can be displayed more accurately, efficiently and visually.
In order to acquire pollutant concentration data of any grid point, traversing coordinate points of a grid of a horizontal display surface and all grid points of a vertical display surface from a three-dimensional data grid, and setting a concentration value of the position of each grid point of the horizontal display surface and the vertical display surface by means of spatial interpolation (trilinear interpolation). In an alternative embodiment, as shown in FIG. 3, a specified three-dimensional grid and specified grid points P within the specified three-dimensional grid are obtained; wherein, the four grid points of one face of the designated three-dimensional grid are A, B, C, D, and the four grid points of the corresponding face of the one face are E, F, G, H;
acquiring environmental data of a grid point A, a grid point B, a grid point C, a grid point D, a grid point E, a grid point F, a grid point G and a grid point H through the NAQPMS system;
the environmental data of the specified grid point P is calculated by the following formula:
P.v=A.v*Aq+B.v*Bq+C.v*Cq+D.v*Dq+E.v*Eq+F.v*Fq+G.v*Gq+H.v*Hq;
where P.v denotes environment data of the specified grid point P, A.v denotes environment data of the grid point A, B.v denotes environment data of the grid point B, C.v denotes environment data of the grid point C, D.v denotes environment data of the grid point D, E.v denotes environment data of the grid point E, F.v denotes environment data of the grid point F, G.v denotes environment data of the grid point G, H.v denotes environment data of the grid point H, Aq denotes a contamination weight of the specified grid point P by the grid point A, Bq denotes a contamination weight of the specified grid point P by the grid point B, Cq denotes a contamination weight of the specified grid point P by the grid point C, Dq denotes a contamination weight of the specified grid point P by the grid point D, Eq denotes a contamination weight of the specified grid point P by the grid point E, Fq denotes a contamination weight of the specified grid point P by the grid point F, Gq denotes a contamination weight of the specified grid point P by the grid point G, hq represents the pollution weight of grid point H to the designated grid point P;
wherein the content of the first and second substances,
Aq=(1–Qx)*(1–Qy)*(1–Qz);
Bq=Qx*(1–Qy)*(1–Qz);
Cq=Qx*Qy*(1–Qz);
Dq=(1–Qx)*Qy*(1–Qz);
Eq=(1–Qx)*(1–Qy)*Qz;
Fq=Qx*(1–Qy)*Qz;
Gq=Qx*Qy*Qz;
Hq=(1–Qx)*Qy*Qz;
the weight of the mesh point pair P point in each direction is calculated as follows:
qx ═ x-A.x)/Lx; the x value is larger than the weight of the grid point at one side of the X value of the P point in the x direction, and the corresponding grid point at the other side is 1-Qx;
qy ═ y-A.y)/Ly; the y value is larger than the weight of the grid point at one side of the y value of the P point in the y direction, and the corresponding grid point at the other side is 1-Qy;
qz ═ z-A.z/Lz; the z value is larger than the weight of the grid point on one side of the z value of the P point in the z direction, and the corresponding grid point on the other side is 1-Qz;
wherein the content of the first and second substances,
Lx=G.x–A.x;
Ly=G.y–A.y;
Lz=G.z–A.z;
lx, Ly and Lz are respectively the length, width and height of the designated three-dimensional grid, A.x, A.y and A.z are respectively the X-axis coordinate of the point A, the Y-axis coordinate of the point A and the Z-axis coordinate of the point A, G.x, G.y and G.z are respectively the X-axis coordinate of the point G, the Y-axis coordinate of the point G and the Z-axis coordinate of the point G, and X, Y and Z are respectively the X, Y and Z coordinate values of the designated grid point P in the three-dimensional coordinates. Through this optional embodiment, can acquire the pollutant concentration data of arbitrary grid point in the region, promoted the richness and the degree of accuracy of data.
In an optional embodiment, a first data display area parallel to the ground and a second data display area perpendicular to the first data display area are obtained according to longitude and latitude coordinates of a designated area, the first data display area and the second data display area are correspondingly drawn to a virtual earth, environmental data of each grid in the first data display area and the second data display area are obtained through a NAQPMS system, and the environmental data of each grid in the first data display area and the second data display area are displayed in the virtual earth.
The Virtual globe (Virtual globe) is a three-dimensional software model representing the globe or another world, and can provide users with functions of freely moving the environment and changing the viewing angle and position. In contrast to a traditional globe, the virtual earth can provide a user with different views of the surface of the earth for different people, such as geographic features, artificial features (e.g., roads, buildings), or abstract data similar to population quantities.
In recent years, as the performance of a Graphics Processing Unit (GPU) is improved remarkably, the 3D Processing capability of a general computer is increased with the increase of water, so that more and more virtual reality technologies are applied to people's real life, such as VR, 3D printing, three-dimensional maps, and the like. As the virtual reality technology is gradually mature, virtual Earth such as WorldWind, Google Earth, and ceium is also rapidly developed and applied to various industries. Compared with a common two-dimensional map, the virtual earth simulates the earth environment where people are located more truly, and is applied to the industries of aerospace, environmental protection, weather and the like.
As shown in fig. 4, any input curve (any Line with many intermediate points, which may be a boundary Line of a heavily polluted area, a boundary Line of a province, etc.) parallel to the spherical surface and having n intermediate points is used as a basic path for displaying pollutants in the vertical direction, and a curve Line may be expressed as: line ═ point0, point1, …, point ], point (x, y, z) is a point on a spherical surface, which is a three-dimensional coordinate point in a cisium Cartesian space rectangular coordinate system (Cartesian 3). As shown in fig. 5, a path from beijing-chongqing-guangzhou is input, and the curve Line can be expressed as: line ═ a, B, C ], where a stands for beijing, B stands for chongqing, C stands for guangzhou. The cartesian coordinates of all points in the curve Line are converted into longitude and latitude coordinates, then the curve Line is split into a certain number of Line segments, for example, the Line segments can be split into 100 segments, and the grid of the vertical display surface is constructed according to dynamic adjustment of a data range and a drawing range. According to the first layer mode data grid, wherein the first layer mode data is the data of the ground layer, the pollutant forecast data of each place is generally expressed by the mode data of the first layer, the grid structure of the horizontal display surface is constructed, and the data is set for the grid structure. The grid data of the horizontal display surface is finally drawn to the surface of the sphere of the virtual earth, and the shape of the area range is changed from a rectangle in the two-dimensional map to the surface of the sphere, so that the grid data of the horizontal display surface is distributed in a sector shape, as shown in fig. 6.
By the method for visualizing the environmental data based on the three-dimensional virtual earth, disclosed by the embodiment of the invention, pollution situations shown in the virtual earth are shown in fig. 7.
The above step S103 involves acquiring environment data of each mesh in the three-dimensional mesh model through the NAQPMS system, and in an optional embodiment, the environment data of the second specified mesh is acquired through the following formula:
Mvalue=Tvalue*QT+Bvalue*QB;
QT=(dataHgt-BottomHgt)/(TopHgt–BottomHgt);
QB=(TopHgt-dataHgt)/(TopHgt–BottomHgt);
wherein Mvalue represents environment data of the second specified grid, Tvalue represents environment data of the first grid, QT represents concentration contribution factor of the first grid, Bvalue represents environment data of the second grid, QB represents concentration contribution factor of the second grid, dataHgt represents height of the second specified grid, TopHgt represents height of the first grid, and BottomHgt represents height of the second grid; wherein the layer height of the second designated grid is between the layer heights of the first and second grids, and the layer height of the first grid is higher than the layer height of the second grid. In a specific optional embodiment, the elevation data grids of the forecast area are copied as the data display layer (the data display layer can be drawn by setting density values and colors), and according to the conditions input by the user (for example, the minimum value is 1, the maximum value is 100, the height values are increased at intervals of 0.01 and represent coefficients of height increase of the data display layer), the elevation values of all grid points are integrally increased by corresponding elevation values, so that the data display layer is raised. And (3) making a straight line parallel to the Z axis through the longitude and latitude coordinates of the data display grid points, and determining the elevation value of the intersection point of the line and each layer of elevation data. And calculating a spatial linear interpolation factor through the elevation value, and calculating a pollutant forecast concentration value of the data display grid point. As shown in fig. 8, the M point is a grid point in the data presentation layer, which is located at an altitude dataHgt of 1000M, a mode forecast data T point higher than the M point, a altitude TopHgt of 1200M, and a mode forecast number lower than the M pointAccording to the point B, the altitude is BottomHgt which is 500 m. It is assumed that the concentration Tvalue of PM2.5 at point T is 10. mu.g/m3The concentration Bvalue of PM2.5 at point B is 62. mu.g/m3. Firstly, calculating concentration contribution factors of the T point and the B point to M according to the height difference, wherein the farther the height distance is, the smaller the contribution is:
QB=(TopHgt-dataHgt)/(TopHgt–BottomHgt);
QT=(dataHgt-BottomHgt)/(TopHgt–BottomHgt);
therefore, the pollutant concentration value Mvalue at point M can be calculated by the above-mentioned factors,
Mvalue=Tvalue*QT+Bvalue*QB;
in this example, the calculation result is that Mvalue 10 × 0.714+62 × 0.286 × 24.9(μ g/m)3) (ii) a By the above method, the density value of the point can be set to the grid points of all the data presentation layers.
The step S106 is related to displaying the environment data of each grid through the three-dimensional grid model, obtaining a color index corresponding to the environment data of each grid, and obtaining a color corresponding to the environment data of each grid according to the color index, wherein the color index corresponding to the environment data of each grid is obtained through the following formula:
cIndex=((V–Nmin)/(Nmax–Nmin))*CLength;
wherein, cndex represents a color index, V represents environment data corresponding to each grid, Nmin represents a minimum value of the environment data, Nmax represents a maximum value of the environment data, and CLength represents the number of colors of the rendering legend of the environmental pollution. The color index and the color of the environment data of each grid are preset, for example, the number of 1-10 is set for 10 grids, 10 colors are put in turn, and the 6 th color is obtained each time the color is corresponding to the calculated number, for example, the number is 6.
Specifically, for example, the concentration of the pollutant (PM2.5) at any one grid point P is N, and 0 in this optional embodiment<=N<=350μg/m3By setting the rendering color of the concentration of the pollutant and rendering the concentration of the pollutant in the optional embodimentThe degree distribution is shown in fig. 9.
In this embodiment, a three-dimensional environment data visualization apparatus is further provided, and the apparatus is used to implement the foregoing embodiments and preferred embodiments, which have already been described and are not described again. As used hereinafter, the term "module" is a combination of software and/or hardware that can implement a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated.
The present embodiment provides an environmental data three-dimensional visualization apparatus, as shown in fig. 10, including:
a first constructing module 101, configured to construct a three-dimensional space coordinate system corresponding to the designated area;
a second constructing module 102, configured to construct a three-dimensional mesh model according to the NAQPMS system based on the three-dimensional spatial coordinate system; wherein, the range of coordinate axes and the grid size of the three-dimensional grid model correspond to the NAQPMS system;
a first obtaining module 103, configured to obtain environment data of each grid in the three-dimensional grid model through the NAQPMS system;
a second obtaining module 104, configured to obtain first environment data of the first specified grid at a first time and second environment data of the first specified grid at a second time;
a third obtaining module 105, configured to obtain environment data at a third time according to the first environment data and the second environment data through the following formula; wherein the third moment is a moment between the first moment and the second moment;
Figure BDA0002845693330000211
a represents the first time, B represents the second time, T represents the third time, Tn represents the environment data at the third time, An represents the first environment data, and Bn represents the second environment data;
a first display module 106, configured to display the environment data of each grid through the three-dimensional grid model.
Optionally, the apparatus further comprises:
the fourth acquisition module is used for acquiring the specified three-dimensional grid and the specified grid point P in the specified three-dimensional grid; wherein, the four grid points of one face of the designated three-dimensional grid are A, B, C, D, and the four grid points of the corresponding face of the one face are E, F, G, H;
a fifth obtaining module, configured to obtain, by using the NAQPMS system, environment data of the grid point a, the grid point B, the grid point C, the grid point D, the grid point E, the grid point F, the grid point G, and the grid point H;
a calculating module, configured to calculate the environment data of the specified grid point P by the following formula:
P.v=A.v*Aq+B.v*Bq+C.v*Cq+D.v*Dq+E.v*Eq+F.v*Fq+G.v*Gq+H.v*Hq;
where P.v denotes environment data of the specified grid point P, A.v denotes environment data of the grid point A, B.v denotes environment data of the grid point B, C.v denotes environment data of the grid point C, D.v denotes environment data of the grid point D, E.v denotes environment data of the grid point E, F.v denotes environment data of the grid point F, G.v denotes environment data of the grid point G, H.v denotes environment data of the grid point H, Aq denotes a contamination weight of the specified grid point P by the grid point A, Bq denotes a contamination weight of the specified grid point P by the grid point B, Cq denotes a contamination weight of the specified grid point P by the grid point C, Dq denotes a contamination weight of the specified grid point P by the grid point D, Eq denotes a contamination weight of the specified grid point P by the grid point E, Fq denotes a contamination weight of the specified grid point P by the grid point F, Gq denotes a contamination weight of the specified grid point P by the grid point G, hq represents the pollution weight of grid point H to the designated grid point P;
wherein the content of the first and second substances,
Aq=(1–Qx)*(1–Qy)*(1–Qz);
Bq=Qx*(1–Qy)*(1–Qz);
Cq=Qx*Qy*(1–Qz);
Dq=(1–Qx)*Qy*(1–Qz);
Eq=(1–Qx)*(1–Qy)*Qz;
Fq=Qx*(1–Qy)*Qz;
Gq=Qx*Qy*Qz;
Hq=(1–Qx)*Qy*Qz;
wherein the content of the first and second substances,
Qx=(x–A.x)/Lx;
Qy=(y–A.y)/Ly;
Qz=(z–A.z)/Lz;
wherein the content of the first and second substances,
Lx=G.x–A.x;
Ly=G.y–A.y;
Lz=G.z–A.z;
lx, Ly and Lz are respectively the length, width and height of the designated three-dimensional grid, A.x, A.y and A.z are respectively the X-axis coordinate of the point A, the Y-axis coordinate of the point A and the Z-axis coordinate of the point A, G.x, G.y and G.z are respectively the X-axis coordinate of the point G, the Y-axis coordinate of the point G and the Z-axis coordinate of the point G, and X, Y and Z are respectively the X, Y and Z coordinate values of the designated grid point P in the three-dimensional coordinates.
Optionally, the apparatus further comprises:
the sixth acquisition module is used for acquiring a first data display area parallel to the ground and a second data display area which is correspondingly vertical to the first data display area according to the longitude and latitude coordinates of the designated area;
the drawing module is used for drawing the first data display area and the second data display area to the virtual earth correspondingly;
a seventh obtaining module, configured to obtain, by the NAQPMS system, environment data of each grid in the first data presentation area and the second data presentation area;
and the second display module is used for displaying the environment data of each grid in the first data display area and the second data display area in the virtual earth.
The three-dimensional visualization device of the environment data in this embodiment is presented in the form of a functional unit, where the unit refers to an ASIC circuit, a processor and a memory executing one or more software or fixed programs, and/or other devices that can provide the above-mentioned functions.
Further functional descriptions of the modules are the same as those of the corresponding embodiments, and are not repeated herein.
An embodiment of the present invention further provides an electronic device, which includes the three-dimensional visualization apparatus for environment data shown in fig. 10.
Referring to fig. 11, fig. 11 is a schematic structural diagram of an electronic device according to an alternative embodiment of the present invention, and as shown in fig. 11, the terminal may include: at least one processor 1101, such as a CPU (Central Processing Unit), at least one communication interface 1103, memory 1104, and at least one communication bus 1102. Wherein a communication bus 1102 is used to enable connective communication between these components. The communication interface 1103 may include a Display screen (Display) and a Keyboard (Keyboard), and the optional communication interface 1103 may also include a standard wired interface and a standard wireless interface. The Memory 1104 may be a Random Access Memory (RAM) or a non-volatile Memory (non-volatile Memory), such as at least one disk Memory. The memory 1104 may optionally be at least one memory device located remotely from the processor 1101. Wherein the processor 1101 may be combined with the apparatus described in fig. 10, the memory 1104 stores an application program therein, and the processor 1101 calls the program code stored in the memory 1104 for executing any one of the above-mentioned three-dimensional visualization methods of the environment data.
The communication bus 1102 may be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus. The communication bus 1102 may be divided into an address bus, a data bus, a control bus, and the like. For ease of illustration, only one thick line is shown in FIG. 11, but this is not intended to represent only one bus or type of bus.
The memory 1104 may include a volatile memory (RAM), such as a random-access memory (RAM); the memory may also include a non-volatile memory (english: non-volatile memory), such as a flash memory (english: flash memory), a hard disk (english: hard disk drive, abbreviated: HDD) or a solid-state drive (english: SSD); the memory 1104 may also comprise a combination of memories of the kind described above.
The processor 1101 may be a Central Processing Unit (CPU), a Network Processor (NP), or a combination of a CPU and an NP.
The processor 1101 may further include a hardware chip. The hardware chip may be an application-specific integrated circuit (ASIC), a Programmable Logic Device (PLD), or a combination thereof. The PLD may be a Complex Programmable Logic Device (CPLD), a field-programmable gate array (FPGA), a General Array Logic (GAL), or any combination thereof.
Optionally, memory 1104 is also used to store program instructions. The processor 1101 may call program instructions to implement the method for three-dimensional visualization of environmental data as shown in the embodiment of fig. 1 of the present application.
The embodiment of the invention also provides a non-transitory computer storage medium, wherein the computer storage medium stores computer executable instructions, and the computer executable instructions can execute the three-dimensional visualization method of the environmental data in any method embodiment. The storage medium may be a magnetic Disk, an optical Disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a Flash Memory (Flash Memory), a Hard Disk (Hard Disk Drive, abbreviated as HDD) or a Solid State Drive (SSD), etc.; the storage medium may also comprise a combination of memories of the kind described above.
Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the invention, and such modifications and variations fall within the scope defined by the appended claims.

Claims (10)

1. A three-dimensional visualization method for environmental data is characterized by comprising the following steps:
constructing a three-dimensional space coordinate system corresponding to the designated area;
constructing a three-dimensional grid model according to a NAQPMS system based on the three-dimensional space coordinate system; wherein the range of coordinate axes and the mesh size of the three-dimensional mesh model correspond to the NAQPMS system;
acquiring environmental data of each grid in the three-dimensional grid model through the NAQPMS system;
acquiring first environment data of a first designated grid at a first moment and second environment data of the first designated grid at a second moment;
acquiring environmental data of a third moment according to the first environmental data and the second environmental data through the following formula; the third moment is a moment between the first moment and the second moment;
Figure FDA0002845693320000011
a represents the first time, B represents the second time, T represents the third time, Tn represents the environment data at the third time, An represents the first environment data, and Bn represents the second environment data;
and displaying the environment data of each grid through the three-dimensional grid model.
2. The method for three-dimensional visualization of environmental data according to claim 1, further comprising:
acquiring a designated three-dimensional grid and a designated grid point P in the designated three-dimensional grid; wherein, the four grid points of one face of the specified three-dimensional grid are A, B, C, D, and the four grid points of the corresponding face of the one face are E, F, G, H;
acquiring environment data of a grid point A, a grid point B, a grid point C, a grid point D, a grid point E, a grid point F, a grid point G and a grid point H through the NAQPMS system;
the environment data of the specified grid point P is calculated by the following formula:
P.v=A.v*Aq+B.v*Bq+C.v*Cq+D.v*Dq+E.v*Eq+F.v*Fq+G.v*Gq+H.v*Hq;
where P.v denotes the environmental data of the specified grid point P, A.v denotes the environmental data of grid point A, B.v denotes the environmental data of grid point B, C.v denotes the environmental data of grid point C, D.v denotes the environmental data of grid point D, E.v denotes the environmental data of grid point E, F.v denotes the environmental data of grid point F, G.v denotes the environmental data of grid point G, H.v denotes the environmental data of grid point H, Aq denotes the pollution weight of grid point A to the specified grid point P, Bq denotes the pollution weight of grid point B to the specified grid point P, Cq denotes the pollution weight of grid point C to the specified grid point P, Dq denotes the pollution weight of grid point D to the specified grid point P, Eq denotes the pollution weight of grid point E to the specified grid point P, Fq denotes the pollution weight of grid point F to the specified grid point P, Gq denotes the pollution weight of grid point G to the specified grid point P, hq represents the pollution weight of grid point H to the designated grid point P;
wherein the content of the first and second substances,
Aq=(1–Qx)*(1–Qy)*(1–Qz);
Bq=Qx*(1–Qy)*(1–Qz);
Cq=Qx*Qy*(1–Qz);
Dq=(1–Qx)*Qy*(1–Qz);
Eq=(1–Qx)*(1–Qy)*Qz;
Fq=Qx*(1–Qy)*Qz;
Gq=Qx*Qy*Qz;
Hq=(1–Qx)*Qy*Qz;
wherein the content of the first and second substances,
Qx=(x–A.x)/Lx;
Qy=(y–A.y)/Ly;
Qz=(z–A.z)/Lz;
wherein the content of the first and second substances,
Lx=G.x–A.x;
Ly=G.y–A.y;
Lz=G.z–A.z;
lx, Ly and Lz are respectively the length, width and height of the designated three-dimensional grid, A.x, A.y and A.z are respectively the X-axis coordinate of the point A, the Y-axis coordinate of the point A and the Z-axis coordinate of the point A, G.x, G.y and G.z are respectively the X-axis coordinate of the point G, the Y-axis coordinate of the point G and the Z-axis coordinate of the point G, and X, Y and Z are respectively the X, Y and Z coordinate values of the designated grid point P in the three-dimensional coordinates.
3. The method for three-dimensional visualization of environmental data according to claim 1, further comprising:
acquiring a first data display area parallel to the ground and a second data display area correspondingly perpendicular to the first data display area according to the longitude and latitude coordinates of the designated area;
correspondingly drawing the first data display area and the second data display area to a virtual earth;
acquiring environment data of each grid in the first data presentation area and the second data presentation area through the NAQPMS system;
and displaying the environment data of each grid in the first data display area and the second data display area in the virtual earth.
4. The method for three-dimensional visualization of environmental data according to claim 1, further comprising:
obtaining environmental data for the second specified grid by:
Mvalue=Tvalue*QT+Bvalue*QB;
QT=(dataHgt-BottomHgt)/(TopHgt–BottomHgt);
QB=(TopHgt-dataHgt)/(TopHgt–BottomHgt);
wherein Mvalue represents environment data of the second specified grid, Tvalue represents environment data of a first grid, QT represents concentration contribution factor of the first grid, Bvalue represents environment data of a second grid, QB represents concentration contribution factor of the second grid, dataHgt represents height of the second specified grid, TopHgt represents height of the first grid, and BottomHgt represents height of the second grid; wherein the layer height of the second designated mesh is between the layer heights of the first mesh and the second mesh, and the layer height of the first mesh is higher than the layer height of the second mesh.
5. The method for three-dimensional visualization of environment data according to any of claims 1 to 4, wherein the displaying of environment data of each mesh by the three-dimensional mesh model comprises:
acquiring color indexes corresponding to the environment data of each grid;
acquiring colors corresponding to the environment data of each grid according to the color indexes;
wherein the color index corresponding to the environment data of each grid is obtained by the following formula:
cIndex=((V–Nmin)/(Nmax–Nmin))*CLength;
wherein, cndex represents a color index, V represents environment data corresponding to each grid, Nmin represents a minimum value of the environment data, Nmax represents a maximum value of the environment data, and CLength represents the number of colors of the rendering legend of the environmental pollution.
6. An apparatus for three-dimensional visualization of environmental data, comprising:
the first construction module is used for constructing a three-dimensional space coordinate system corresponding to the designated area;
the second construction module is used for constructing a three-dimensional grid model based on the three-dimensional space coordinate system according to a NAQPMS system; wherein the range of coordinate axes and the mesh size of the three-dimensional mesh model correspond to the NAQPMS system;
a first obtaining module, configured to obtain, by the NAQPMS system, environment data of each grid in the three-dimensional grid model;
the second acquisition module is used for acquiring first environmental data of the first designated grid at a first moment and second environmental data of the first designated grid at a second moment;
the third obtaining module is used for obtaining the environmental data of a third moment according to the first environmental data and the second environmental data through the following formula; the third moment is a moment between the first moment and the second moment;
Figure FDA0002845693320000051
a represents the first time, B represents the second time, T represents the third time, Tn represents the environment data at the third time, An represents the first environment data, and Bn represents the second environment data;
and the first display module is used for displaying the environment data of each grid through the three-dimensional grid model.
7. The apparatus for three-dimensional visualization of environmental data according to claim 6, wherein said apparatus further comprises:
the fourth acquisition module is used for acquiring a specified three-dimensional grid and a specified grid point P in the specified three-dimensional grid; wherein, the four grid points of one face of the specified three-dimensional grid are A, B, C, D, and the four grid points of the corresponding face of the one face are E, F, G, H;
a fifth obtaining module, configured to obtain, by the NAQPMS system, environment data of the grid point a, the grid point B, the grid point C, the grid point D, the grid point E, the grid point F, the grid point G, and the grid point H;
a calculating module, configured to calculate the environment data of the specified grid point P by using the following formula:
P.v=A.v*Aq+B.v*Bq+C.v*Cq+D.v*Dq+E.v*Eq+F.v*Fq+G.v*Gq+H.v*Hq;
where P.v denotes the environmental data of the specified grid point P, A.v denotes the environmental data of grid point A, B.v denotes the environmental data of grid point B, C.v denotes the environmental data of grid point C, D.v denotes the environmental data of grid point D, E.v denotes the environmental data of grid point E, F.v denotes the environmental data of grid point F, G.v denotes the environmental data of grid point G, H.v denotes the environmental data of grid point H, Aq denotes the pollution weight of grid point A to the specified grid point P, Bq denotes the pollution weight of grid point B to the specified grid point P, Cq denotes the pollution weight of grid point C to the specified grid point P, Dq denotes the pollution weight of grid point D to the specified grid point P, Eq denotes the pollution weight of grid point E to the specified grid point P, Fq denotes the pollution weight of grid point F to the specified grid point P, Gq denotes the pollution weight of grid point G to the specified grid point P, hq represents the pollution weight of grid point H to the designated grid point P;
wherein the content of the first and second substances,
Aq=(1–Qx)*(1–Qy)*(1–Qz);
Bq=Qx*(1–Qy)*(1–Qz);
Cq=Qx*Qy*(1–Qz);
Dq=(1–Qx)*Qy*(1–Qz);
Eq=(1–Qx)*(1–Qy)*Qz;
Fq=Qx*(1–Qy)*Qz;
Gq=Qx*Qy*Qz;
Hq=(1–Qx)*Qy*Qz;
wherein the content of the first and second substances,
Qx=(x–A.x)/Lx;
Qy=(y–A.y)/Ly;
Qz=(z–A.z)/Lz;
wherein the content of the first and second substances,
Lx=G.x–A.x;
Ly=G.y–A.y;
Lz=G.z–A.z;
lx, Ly and Lz are respectively the length, width and height of the designated three-dimensional grid, A.x, A.y and A.z are respectively the X-axis coordinate of the point A, the Y-axis coordinate of the point A and the Z-axis coordinate of the point A, G.x, G.y and G.z are respectively the X-axis coordinate of the point G, the Y-axis coordinate of the point G and the Z-axis coordinate of the point G, and X, Y and Z are respectively the X, Y and Z coordinate values of the designated grid point P in the three-dimensional coordinates.
8. The apparatus for three-dimensional visualization of environmental data according to claim 6, wherein said apparatus further comprises:
the sixth acquisition module is used for acquiring a first data display area parallel to the ground and a second data display area which is correspondingly vertical to the first data display area according to the longitude and latitude coordinates of the designated area;
the drawing module is used for correspondingly drawing the first data display area and the second data display area to a virtual earth;
a seventh obtaining module, configured to obtain, by the NAQPMS system, environment data of each grid in the first data presentation area and the second data presentation area;
and the second display module is used for displaying the environment data of each grid in the first data display area and the second data display area in the virtual earth.
9. An electronic device, comprising: at least one processor; and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the at least one processor to cause the at least one processor to perform the method for three-dimensional visualization of environmental data as set forth in any one of claims 1-5.
10. A computer-readable storage medium having stored thereon computer instructions, which when executed by a processor, implement the method for three-dimensional visualization of environmental data according to any of the claims 1-5.
CN202011508725.3A 2020-12-18 2020-12-18 Three-dimensional visualization method and device for environmental data Pending CN112598801A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011508725.3A CN112598801A (en) 2020-12-18 2020-12-18 Three-dimensional visualization method and device for environmental data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011508725.3A CN112598801A (en) 2020-12-18 2020-12-18 Three-dimensional visualization method and device for environmental data

Publications (1)

Publication Number Publication Date
CN112598801A true CN112598801A (en) 2021-04-02

Family

ID=75199488

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011508725.3A Pending CN112598801A (en) 2020-12-18 2020-12-18 Three-dimensional visualization method and device for environmental data

Country Status (1)

Country Link
CN (1) CN112598801A (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109948840A (en) * 2019-03-08 2019-06-28 宁波市气象台 A kind of Urban Air Pollution Methods
CN112037290A (en) * 2020-08-28 2020-12-04 中科三清科技有限公司 Data processing method, device and equipment

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109948840A (en) * 2019-03-08 2019-06-28 宁波市气象台 A kind of Urban Air Pollution Methods
CN112037290A (en) * 2020-08-28 2020-12-04 中科三清科技有限公司 Data processing method, device and equipment

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
朱帅: "虚拟地球技术在三维空气质量数据可视化中的应用尝试", 《2017中国环境科学学会科学与技术年会论文集(第一卷)》 *

Similar Documents

Publication Publication Date Title
KR102360660B1 (en) Map data processing method, computer device and storage medium
CN112559665A (en) Air quality mode forecast data visualization method and device
KR101994317B1 (en) Navigation device, method of determining a height coordinate and method of generating a database
CN102117500A (en) Three-dimensional modeling method and system
CN107918957B (en) Three-dimensional building model simplification method capable of keeping structure and texture characteristics
CN111090716A (en) Vector tile data processing method, device, equipment and storage medium
CN102117494A (en) Method and device for delaminating and blocking data in three-dimensional modeling
EP2589933B1 (en) Navigation device, method of predicting a visibility of a triangular face in an electronic map view
CN112581615A (en) Environment data visualization method and device based on three-dimensional virtual earth
CN112598800A (en) Air quality mode forecast data acquisition method and device based on terrain
CN114138265B (en) Visualization method based on digital twinning
CN114239271A (en) Atmospheric pollution early warning method, device and equipment
CN113157988B (en) Method and device for representing geographic information by OFD format file
CN112598801A (en) Three-dimensional visualization method and device for environmental data
CN113066176A (en) Map data processing method, device, equipment and storage medium
CN114840896B (en) Waterlogging and ponding simulation analysis method based on urban road BIM
CN115187709A (en) Geographic model processing method and device, electronic equipment and readable storage medium
CN113495935A (en) Terrain slicing method and system based on irregular triangularization mesh generation
CN111383334B (en) System and method for rendering objects
CN113066179A (en) Map data processing method, device, equipment and storage medium
JP2003228599A (en) Cae device for supporting analytical work using numerical map data information
CN113066150A (en) Map data processing method, device, equipment and storage medium
CN112446959B (en) Oblique photography model optimization method
Liu et al. Fusing multiscale charts into 3D ENC systems based on underwater topography and remote sensing image
CN116363327B (en) Voxel map generation method and system

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210402