CN115760521A - Method for fusing population distribution data and three-dimensional model scene - Google Patents
Method for fusing population distribution data and three-dimensional model scene Download PDFInfo
- Publication number
- CN115760521A CN115760521A CN202211440604.9A CN202211440604A CN115760521A CN 115760521 A CN115760521 A CN 115760521A CN 202211440604 A CN202211440604 A CN 202211440604A CN 115760521 A CN115760521 A CN 115760521A
- Authority
- CN
- China
- Prior art keywords
- data
- signaling
- scene
- dimensional
- layer range
- 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
Links
Images
Landscapes
- Processing Or Creating Images (AREA)
Abstract
The invention relates to the technical field of population distribution data analysis, in particular to a method for fusing population distribution data and a three-dimensional model scene, which comprises the steps of building a three-dimensional city base, standardizing data, preprocessing data, rendering thermodynamic diagram of signaling data, rendering thermodynamic diagram of a three-dimensional scene building and rendering the signaling data in real time, wherein a vector geographic element file of a WGS84 coordinate system generates a model in modeling software, the obtained mobile phone signaling data is converted into a GeoJSON file of the WGS84 coordinate system through processing of a QG I S tool, and the standardization of the data is kept; according to the method for fusing the population distribution data with the three-dimensional scene through the program, the three-dimensional city scene base is built by accessing the real geographic basic element data, the population thermodynamic diagram is rendered by combining the mobile phone signaling data, and the three-dimensional city scene base is rendered on the three-dimensional building through the shader.
Description
Technical Field
The invention relates to a method, in particular to a method for fusing population distribution data and a three-dimensional model scene, and belongs to the technical field of population distribution data analysis.
Background
With the development of digital twin technology and the improvement of hardware performance, digital twin application is more and more extensive, and a method, a device and a storage medium for determining population distribution thermodynamic data are provided in the chinese patent of invention with the application number CN202110482704.7, high-resolution regional attribute data is used to determine high-resolution population weight data, and then low-resolution population distribution thermodynamic data is processed based on the high-resolution population weight data to obtain high-resolution population distribution thermodynamic data, so that the accuracy of population distribution thermodynamic data can be improved.
The invention provides a method for detecting the distribution density of dynamic population based on mobile phone positioning data in Chinese invention patent with application number CN201210020492.1, which fully depends on the existing mobile communication network resources, takes mass mobile phone signal data as input, and adopts a discrimination mechanism based on travel track continuous tracking and stay time in an area to realize the detection of the distribution density of the dynamic population. The invention can be used for automatically detecting the dynamic population distribution density information with high frequency.
However, the existing comparison files, including traditional population data distribution, are still based on a two-dimensional presentation mode, cannot meet the requirements of users in digital twin application scenes, and have the problems of flattening and non-intuition. Population distribution data is combined to a three-dimensional scene, city planning and traffic management field decision can be assisted more intuitively and in multiple visual angles, the existing three-dimensional scene fusion based on the population distribution data is mostly in a plane form, data cannot be conveniently and rapidly presented intuitively, and the fusion of three-dimensional buildings and the population distribution data can also be realized through manual modeling mapping, but the mode has the problems of low efficiency and difficulty in updating and maintaining.
The present invention has been made in view of this situation.
Disclosure of Invention
The invention aims to solve the problems and provides a method for fusing population distribution data and a three-dimensional model scene, which is characterized in that the population distribution data and the three-dimensional scene are fused through a program, a three-dimensional city scene base is built by accessing real geographic basic element data, a population thermodynamic diagram is rendered by combining mobile phone signaling data, and finally the population thermodynamic diagram is rendered on a three-dimensional building through a shader.
The invention realizes the aim through the following technical scheme, and a method for fusing population distribution data and a three-dimensional model scene comprises the following steps:
s1: building a three-dimensional city base: real coordinate data are accessed, and a vector geographic element file of a WGS84 coordinate system generates a model in modeling software, wherein the model comprises a road surface model, a water system model, a mountain model and a house surface model, so that mobile phone signaling data are accurately matched with a three-dimensional scene;
s2: data normalization: the acquired mobile phone signaling data is processed and converted into a GeoJSON file of a WGS84 coordinate system through a QGIS tool, and the standardization of the data is kept;
s3: data preprocessing: obtaining the signaling point location data and the vector file of the house roof obtained in the step S2 through a layer range extraction tool of a QGIS (QGIS), entering a step S4 if the layer range of the point location data contains the layer range of the house roof, and supplementing boundary points if the layer range does not belong to an inclusion relation so as to enable the layer range after the signaling point location is supplemented to contain the layer range of the house roof;
s4: thermodynamic rendering of signaling data: converting the standardized signaling data obtained in the step S3 into x and y values based on screen coordinates in a program based on a heatmap.
Firstly, through the signaling point position layer range obtained in s3, 4 vertexes are marked as a, B, C and D, in order to calculate the distance between the vertexes, the longitude and latitude data of the 4 vertexes needs to be converted into mercator plane projection coordinates, and a back projection conversion formula is as follows:
R＝D/(2*π*r*；360)
r = converted geographical radians; d = projection distance before conversion; r = mean radius of the earth; taking a constant 6371004;
secondly, the distances between the signaling layer range in the x direction and the y direction are obtained through the plane projection coordinates of the 4 vertexes, and the formula is as follows:
W＝B.x-A.x；H＝A.y-C.y
a.x = abscissa x of vertex a; a.y = ordinate y of vertex a; b.x = abscissa x of vertex B; c.y = ordinate y of vertex C; w = position distance in x direction of the signaling layer range; h = distance in y direction of signalling layer range
Thirdly, carrying out back projection conversion on the signaling data to obtain the x value and the y value of the plane coordinate of each point, and then carrying out normalization processing based on the screen coordinate, wherein the formula is as follows:
p.x＝(x/W)＊canvasWidth
p.y＝(1-(y/H)＊canvasHeight
p.x = the signaling data is based on the abscissa x of the screen coordinates; p.y = abscissa y of the signaling data based on screen coordinates; canvas width = width of the browser view window (in pixels); canvas height = the height of the browser view window (in pixels);
and fourthly, traversing the signaling data set, marking each point as a weight, and finally obtaining the screen thermodynamic diagram rendering based on the geographic position through heatmap.
S5: performing thermal rendering on the three-dimensional scene building, and replacing the material of the obtained model with the thermodynamic diagram obtained in the step S4 after the model is read and before the scene is rendered;
s6: and (3) rendering the signaling data in real time, and visually seeing the population work distribution condition and the population living condition based on the building in the three-dimensional scene by respectively using the mobile phone signaling data at the noon time period and the evening time period of the workday through the processes from S2 to S5.
WGS84 coordinate system: world geographic System 1984, a coordinate System established for use with the GPS global positioning System.
And (3) QGIS: the desktop GIS software is an extremely open-source desktop GIS software and can provide functions of data display, editing and analysis.
GeoJSON: a format for encoding various geographic data structures is a geographic spatial information data exchange format based on Javascript Object Notation (JSON).
Js is a thermodynamic diagram framework based on canvas open source, and the thermodynamic diagram can be conveniently realized by using the framework.
Mercator plane projection coordinates: mercator Projection (also called mcator Projection, orthoaxial equiangular cylindrical Projection) is an equiangular cylindrical map Projection method.
The invention has the technical effects and advantages that: according to the method for fusing the population distribution data with the three-dimensional scene through the program, the three-dimensional city scene base is built by accessing the real geographic basic element data, the population thermodynamic diagram is rendered by combining the mobile phone signaling data, and the three-dimensional city scene base is rendered on the three-dimensional building through the shader.
Drawings
FIG. 1 is a detailed flow chart of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
Referring to fig. 1, a method for fusing population distribution data with a three-dimensional model scene includes the following steps:
s1: building a three-dimensional city base: and accessing real coordinate data, and generating a model comprising a road surface model, a water system model, a mountain model and a house model from the vector geographic element file of the WGS84 coordinate system in modeling software, so that the mobile phone signaling data is accurately matched with the three-dimensional scene.
S2: data normalization: and (3) converting the acquired mobile phone signaling data into a GeoJSON file of a WGS84 coordinate system through QGIS tool processing, and keeping the standardization of the data.
S3: data preprocessing: and (3) obtaining the signaling point location data and the vector file of the house roof from the S2 through a layer range extraction tool of a QGIS (QGIS), entering S4 if the layer range of the point location data contains the layer range of the house roof, and supplementing boundary points if the layer range does not belong to the inclusion relation so as to enable the layer range after signaling point location is supplemented to contain the layer range of the house roof.
S4: thermodynamic rendering of signaling data: converting the standardized signaling data obtained in the step S3 into x and y values based on screen coordinates in a program based on a heatmap.
Firstly, through the signaling point position layer range obtained in s3, 4 vertexes are marked as a, B, C and D, in order to calculate the distance between the vertexes, the longitude and latitude data of the 4 vertexes needs to be converted into mercator plane projection coordinates, and a back projection conversion formula is as follows:
R＝D/(2*π*r*360)
r = converted geographical radians; d = projection distance before conversion; r = mean radius of the earth; take the constant 6371004.
Secondly, the distance between the signaling layer range in the x direction and the distance between the signaling layer range in the y direction are obtained through the plane projection coordinates of 4 vertexes, and the formula is as follows:
W＝B.x-A.xH＝A.y-C.y
a.x = abscissa x of vertex a; a.y = ordinate y of vertex a; b.x = abscissa x of vertex B; c.y = ordinate y of vertex C; w = position distance of the signaling layer range in the x direction; h = distance in y direction of the signaling layer range.
Thirdly, carrying out back projection conversion on the signaling data to obtain the x value and the y value of the plane coordinate of each point, and then carrying out normalization processing based on the screen coordinate, wherein the formula is as follows:
p.x＝(x/W)＊canvasWidth
p.y＝(1-(y/H)＊canvasHeight
p.x = abscissa x of the signaling data based on screen coordinates; p.y = abscissa y of the signaling data based on screen coordinates; canvas width = width of the browser view window (in pixels); canvas height = height of the viewer view window (in pixels).
And fourthly, traversing the signaling data set, marking each point as a weight, and finally obtaining the screen thermodynamic diagram rendering based on the geographic position through heatmap.
S5: and (4) performing thermal rendering on the three-dimensional scene building, and replacing the material of the obtained model with the thermodynamic diagram obtained in the step (S4) after the model is read and before the scene is rendered.
S6: and (4) rendering the signaling data in real time, and visually seeing the population work distribution condition and the population living condition based on the building in the three-dimensional scene through the processes from S2 to S5 by respectively using the mobile phone signaling data at the noon time and the evening time of the workday.
Compared with the prior art, the method fuses the population data and the building data, and provides a more intuitive decision basis for city planning and traffic management.
In the present invention, the WGS84 coordinate system: the World Geodetic System 1984 is a coordinate System established for use by the GPS global positioning System.
And (3) QGIS: the desktop GIS software is a very open-source desktop GIS software and can provide functions of displaying, editing and analyzing data.
GeoJSON: a format for coding various geographic data structures is a geographic space information data exchange format based on Javascript Object Notation (JSON).
Js is a thermodynamic diagram framework based on the canvas open source, and the thermodynamic diagram can be conveniently realized by using the framework.
Mercator plane projection coordinates: mercator Projection (also called mcator Projection, orthoaxial equiangular cylindrical Projection) is an equiangular cylindrical map Projection method.
It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned.
Furthermore, it should be understood that although the present description refers to embodiments, not every embodiment may contain only a single embodiment, and such description is for clarity only, and those skilled in the art should integrate the description, and the embodiments may be combined as appropriate to form other embodiments understood by those skilled in the art.
Claims (8)
1. A method for fusing population distribution data with a three-dimensional model scene is characterized by comprising the following steps,
s1: building a three-dimensional city base, and generating a model by using a vector geographic element file of a WGS84 coordinate system in modeling software;
s2: data standardization, namely processing and converting the acquired mobile phone signaling data into a GeoJSON file of a WGS84 coordinate system through a QGIS tool, and keeping the standardization of the data;
s3: preprocessing data, namely obtaining the signaling point location data and the vector file of the house roof obtained in the step S2 through a layer range extraction tool of a QGIS (QGIS), and obtaining respective range surfaces;
s4: rendering the thermodynamic diagrams of the signaling data, and converting the standardized signaling data obtained in the step s3 into x and y values of the signaling data in a program based on screen coordinates of heatmap.js based on the thermodynamic diagram rendering mode of heatmap.js;
s5: performing thermal rendering on the three-dimensional scene building, and replacing the material of the obtained model with the thermodynamic diagram obtained in the step S4 after the model is read and before the scene is rendered;
s6: and (3) rendering the signaling data in real time, and visually seeing the population work distribution condition and the population living condition based on the building in the three-dimensional scene by respectively using the mobile phone signaling data at the noon time period and the evening time period of the workday through the processes from S2 to S5.
2. The method as claimed in claim 1, wherein in step S1, the specific operation flow is accessing real coordinate data, and the vector geographic element file of WGS84 coordinate system is used to generate a model in modeling software, including road surface, water system surface, mountain body, and house surface models, so as to precisely match the mobile phone signaling data with the three-dimensional scene.
3. The method according to claim 1, wherein in step S3, if the layer range of the point location data includes the layer range of the roof, the step proceeds to step S4, and if the layer range does not belong to the inclusion relationship, the boundary point is supplemented, so that the layer range after the signaling point location is supplemented includes the layer range of the roof.
4. The method for fusing the population distribution data with the three-dimensional model scene as recited in claim 1, wherein in step S4, the process comprises the following four steps:
firstly, marking 4 vertexes as A, B, C and D through the signaling point position layer range obtained in s3, and in order to calculate the distance between the vertexes, converting longitude and latitude data of the 4 vertexes into a mercator plane projection coordinate;
secondly, obtaining the distance between the signaling layer range in the x direction and the distance between the signaling layer range in the y direction through the plane projection coordinates of the 4 vertexes;
carrying out back projection conversion on the signaling data to obtain a plane coordinate x value and a plane coordinate y value of each point, and then carrying out normalization processing based on screen coordinates;
and fourthly, traversing the signaling data set, marking each point as a weight, and finally obtaining the screen thermodynamic diagram rendering based on the geographic position through heatmap.
5. The method as claimed in claim 4, wherein in step one, the back projection transformation formula is as follows:
R＝D/(2*π*r*360)
r = converted geographical radians; d = projection distance before conversion; r = mean radius of the earth; take the constant 6371004.
6. The method for fusing the population distribution data with the three-dimensional model scene as claimed in claim 4, wherein in the second step, the specific formula is as follows:
W＝-B.x-A.x；H＝A.y-C.y
a.x = abscissa x of vertex a; a.y = ordinate y of vertex a; b.x = abscissa x of vertex B; c.y = ordinate y of vertex C; w = position distance of the signaling layer range in the x direction; h = distance in y direction of the signaling layer range.
7. The method for fusing the population distribution data with the three-dimensional model scene as claimed in claim 4, wherein in step three, the specific formula is as follows:
p.x＝(x/W)＊canvasWidth
p.y＝(1-(y/H）＊canvasHeight
p.x = abscissa x of the signaling data based on screen coordinates; p.y = abscissa y of the signaling data based on screen coordinates; canvas width = width of the browser view window (in pixels); canvas height = height (in pixels) of the browser view window.
8. The method of claim 1, wherein in step S2, the QGIS processing tool further uses other GIS desktop software including ArcMap or open source space geometry algorithm library including turf.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211440604.9A CN115760521A (en) | 2022-11-17 | 2022-11-17 | Method for fusing population distribution data and three-dimensional model scene |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211440604.9A CN115760521A (en) | 2022-11-17 | 2022-11-17 | Method for fusing population distribution data and three-dimensional model scene |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115760521A true CN115760521A (en) | 2023-03-07 |
Family
ID=85372667
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211440604.9A Pending CN115760521A (en) | 2022-11-17 | 2022-11-17 | Method for fusing population distribution data and three-dimensional model scene |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115760521A (en) |
-
2022
- 2022-11-17 CN CN202211440604.9A patent/CN115760521A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109493407A (en) | Realize the method, apparatus and computer equipment of laser point cloud denseization | |
| WO2021017675A1 (en) | Trajectory data processing method, device and storage medium | |
| Yan et al. | Integration of 3D objects and terrain for 3D modelling supporting the digital twin | |
| CN108921943B (en) | Road three-dimensional model modeling method based on lane-level high-precision map | |
| CN109801371B (en) | Network three-dimensional electronic map construction method based on Cesium | |
| KR101876114B1 (en) | Terminal, server, system for 3d modeling and 3d modeling method using the same | |
| CN110659385B (en) | Fusion method of multi-channel video and three-dimensional GIS scene | |
| TW202117674A (en) | Data processing method and apparatus, storage medium and computer device | |
| Khayyal et al. | Creation and spatial analysis of 3D city modeling based on GIS data | |
| CN115187647A (en) | Vector-based road three-dimensional live-action structured modeling method | |
| Rau et al. | A cost-effective strategy for multi-scale photo-realistic building modeling and web-based 3-D GIS applications in real estate | |
| Chen et al. | A virtual globe-based integration and visualization framework for aboveground and underground 3D spatial objects | |
| Kang et al. | The application of UAV oblique photogrammetry in smart tourism: a case study of Longji terraced scenic SPOT in Guangxi province | |
| CN117351162A (en) | Live-action three-dimensional planning design method based on digital twin | |
| CN113420109B (en) | Method for measuring permeability of street interface, computer and storage medium | |
| Yang et al. | Analysis of visualization technology of 3d spatial geographic information system | |
| CN116842122A (en) | Geographic information system for digital twin city platform | |
| CN115760521A (en) | Method for fusing population distribution data and three-dimensional model scene | |
| CN116091715A (en) | POI (point of interest) point aggregation method based on three-dimensional GIS engine | |
| CN114490907A (en) | Method and device for constructing famous city management database and storage medium | |
| Fu et al. | [Retracted] 3D City Online Visualization and Cluster Architecture for Digital City | |
| Fan et al. | Pano2Geo: An efficient and robust building height estimation model using street-view panoramas | |
| CN110866971B (en) | Method for producing true square orthographic image | |
| Chio et al. | The establishment of 3D LOD2 objectivization building models based on data fusion | |
| CN113627561B (en) | Data fusion 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 |