CN115129401B - Method for realizing non-offset and non-deformation superposition of multi-source data based on Cesum - Google Patents

Method for realizing non-offset and non-deformation superposition of multi-source data based on Cesum Download PDF

Info

Publication number
CN115129401B
CN115129401B CN202211037107.4A CN202211037107A CN115129401B CN 115129401 B CN115129401 B CN 115129401B CN 202211037107 A CN202211037107 A CN 202211037107A CN 115129401 B CN115129401 B CN 115129401B
Authority
CN
China
Prior art keywords
data
coordinates
gaussian
coordinate system
cartesian
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202211037107.4A
Other languages
Chinese (zh)
Other versions
CN115129401A (en
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.)
Geospace Information Technology Co ltd
Original Assignee
Geospace Information 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 Geospace Information Technology Co ltd filed Critical Geospace Information Technology Co ltd
Priority to CN202211037107.4A priority Critical patent/CN115129401B/en
Publication of CN115129401A publication Critical patent/CN115129401A/en
Application granted granted Critical
Publication of CN115129401B publication Critical patent/CN115129401B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00Arrangements for program control, e.g. control units
    • G06F9/06Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/44Arrangements for executing specific programs
    • G06F9/448Execution paradigms, e.g. implementations of programming paradigms
    • G06F9/4482Procedural
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/27Replication, distribution or synchronisation of data between databases or within a distributed database system; Distributed database system architectures therefor
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00Arrangements for program control, e.g. control units
    • G06F9/06Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/44Arrangements for executing specific programs
    • G06F9/445Program loading or initiating
    • G06F9/44521Dynamic linking or loading; Link editing at or after load time, e.g. Java class loading
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00Arrangements for program control, e.g. control units
    • G06F9/06Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/44Arrangements for executing specific programs
    • G06F9/448Execution paradigms, e.g. implementations of programming paradigms
    • G06F9/4488Object-oriented
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics

Abstract

The invention is suitable for the field of geographic information technology Cesium, and provides a method for realizing non-offset and non-deformation superposition of multi-source data based on Cesium, which comprises the following steps: configuring a reference data range of a Gaussian coordinate system in the ceium scene mode; establishing a conversion interface between a Gaussian coordinate system and a Cartesian coordinate system; setting a pyramid loading scheme; for received multi-source data to be loaded, if the multi-source data is model data or vector data, preprocessing is carried out; and under a Gaussian coordinate system, selecting and calling a conversion interface or a pyramid loading scheme according to the type of the multi-source data to load the data. According to the invention, through establishing the conversion interface and setting the pyramid loading scheme, the corresponding scheme is adopted for loading according to different types of multi-source data, and finally the effects of no offset and no deformation of data superposition are achieved, and the requirements of actual engineering are met.

Description

Method for realizing non-offset and non-deformation superposition of multi-source data based on Cesum
Technical Field
The invention belongs to the technical field of geographic information, and particularly relates to a method for realizing non-offset and non-deformation superposition of multi-source data based on Cesum.
Background
WebGL (Web Graphics Library) is an open-source free and cross-platform three-dimensional Graphics API (Application Programming Interface), and is mainly used for rendering three-dimensional data at a browser end. The center rendering engine is a three-dimensional earth and map visualization library implemented based on WebGL, and is designed to display three-dimensional GIS (Geographic Information System) data, such as oblique photography data, vector data, raster data, and terrain data, and so on.
The Cesium supports 3D, 2.5D and 2D rendering modes, the default coordinate system is a WGS84 geodetic coordinate system, the adopted projection mode is longitude and latitude direct projection, and the projection has the characteristics of extremely large deformation, and the deformation is larger when the latitude is higher. Although Cesium has a 2.5D planar three-dimensional mode, it is limited by its WGS84 fundamental coordinate system and cannot load three-dimensional GIS data in a gaussian coordinate system. At present, a gaussian coordinate system is still widely applied to engineering projects due to the characteristic that gaussian projection is not deformed, so that loading of three-dimensional data under the coordinate system is a very urgent requirement, and a set of three-dimensional GIS rendering engine capable of loading the data of the gaussian coordinate system needs to be realized. Because of the characteristics of wide use and open source of Cesium, the function of realizing three-dimensional GIS data loading in a Gaussian coordinate system based on Cesium is selected.
Disclosure of Invention
In view of the above problems, an object of the present invention is to provide a method for implementing offset-free and deformation-free superposition of multi-source data based on cesum, which aims to solve the technical problem that the existing cesum rendering engine cannot load gaussian coordinate system data.
The invention adopts the following technical scheme:
the method for realizing the offset-free and deformation-free superposition of multi-source data based on the Cesum comprises the following steps:
s1, configuring a reference data range of a Gaussian coordinate system in a ceium scene mode;
s2, establishing a conversion interface between a Gaussian coordinate system and a Cartesian coordinate system;
s3, setting a pyramid loading scheme;
s4, for the received multi-source data to be loaded, if the multi-source data is model data or vector data, preprocessing is carried out;
and S5, selecting and calling a conversion interface or pyramid loading scheme according to the type of the multi-source data to load the data in a Gaussian coordinate system.
Further, in step S1, the reference data range represents a rectangular range, and one of the largest ranges in the multi-source data is taken, and includes four values, which are a minimum X value xmin, a minimum Y value ymin, a maximum X value xmax, and a maximum Y value ymax.
Further, in step S2, the conversion interface includes a forward calculation interface for converting the gaussian coordinate system into the cartesian coordinate system, and a backward calculation interface for converting the cartesian coordinate system into the gaussian coordinate system; the calculation interface processing process comprises the following steps: firstly, converting Gaussian coordinates into geographic coordinates by a longitude and latitude direct projection mode, and then converting a geographic coordinate system into Cartesian coordinates; the inverse calculation interface processing process comprises the following steps: firstly, converting Cartesian coordinates into geographic coordinates, and then converting the geographic coordinates into Gaussian coordinates in a longitude and latitude direct projection mode;
the forward interface processing process specifically comprises the following steps:
201. setting a translation reference point [0, 0], and carrying out projection calculation on Gaussian coordinates of multi-source data:
x = r * x0
y = r * y0
wherein r is the long radius of the earth, X0 and Y0 are the gaussian X-coordinate, Y-coordinate of the multi-source data;
202. and (3) calculating the coordinate translation amount:
Tx = (e.xmin+e.xmax) /2–x
Ty = (e.xmin+e.xmax) /2–y
wherein e is a reference data range;
203. calculating the translated coordinates:
Ax = Gx–Tx
Ay = Gy–Ty
204. and (3) back projection calculation:
Cx = Ax / r
Cy = Ay / r
205. calculating to obtain Cartesian coordinates:
cartesian coordinates c are calculated through the interface cartesian.
The reverse calculation interface processing process specifically comprises the following steps:
211. and (3) converting geographic coordinates:
tosradians (c) converts cartesian coordinates to geographic coordinates carto using the interface cartesian provided by ceium;
212. the geographic coordinates carto are subjected to projection calculations:
x’ = carto.x * r
y’ = carto.y * r
213. calculating Gaussian coordinates:
Gx = x’ + Tx
Gy = y’ + Ty 。
further, the specific process of step S3 is as follows:
s31, calculating the offset by combining the reference data range with the map center point:
311. carrying out projection calculation on the central point of the map:
x" = r * x1
y" = r * y1
wherein X1 and Y1 are the X coordinate and the Y coordinate of the central point;
312. and (3) calculating the coordinate translation amount:
Tx" = (e.xmin+e.xmax) /2–x"
Ty" = (e.ymin+e.ymax) /2–y"
wherein e is a reference data range;
s32, carrying out offset superposition operation on the incoming data range to obtain a new data range:
e1".xmin = e1.xmin – Tx"
e1".xmax = e1.xmax – Tx"
e1".ymin = e1.ymin – Ty"
e1".ymax = e1.ymax – Ty"
e1 is the incoming data range and e1 "is the new data range.
Further, in step S4, the model data is preprocessed as follows: converting the model data modeling vertex coordinates into Cartesian coordinates by adopting a forward calculation interface; the vector data is preprocessed as follows: a data range is set for the vector data.
Further, in step S5, if the multi-source data is raster data, obtaining a correct map tile data range according to the pyramid loading scheme, and initiating a request for loading according to a raster service address; if the multi-source data is terrain data, obtaining a correct terrain data range according to the pyramid loading scheme, and rendering the terrain according to a terrain service address; if the multi-source data are vector data, calling a forward calculation interface to convert the coordinates of the vector data into Cartesian coordinates according to the set data range; and if the multi-source data is model data, extracting translation values from a matrix in a requested model description file, combining the translation values into a Cartesian coordinate, and converting the Cartesian coordinate into a Gaussian coordinate through a back calculation interface to obtain a real coordinate of the model in a Gaussian coordinate system.
The beneficial effects of the invention are: according to the method, the conversion interface between the Gaussian coordinate system and the Cartesian coordinate system is established, the pyramid loading scheme is set, the scene is moved to the place with the map center coordinate of 0, the characteristic that no deformation exists at the equator can be fully utilized to achieve the effect of restoring the authenticity of data, finally, the corresponding scheme is adopted for loading according to different types of multi-source data, the effect that data are overlapped without deviation and deformation is finally achieved, and the requirements of actual engineering are met.
Drawings
Fig. 1 is a flowchart of a method for implementing offset-free and deformation-free superposition of multi-source data based on cesum according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The Cesium is limited by a specific coordinate system, the data of a Gaussian coordinate system cannot be loaded, and the Cesium adopts a longitude and latitude direct projection mode. In order to load three-dimensional data under a Gaussian coordinate system, the invention provides a method for superposing multisource three-dimensional GIS data under the Gaussian coordinate without offset and deformation based on Cesium.
In order to explain the technical means of the present invention, the following description will be given by way of specific examples.
Fig. 1 shows a flow of a method for implementing offset-free and deformation-free superposition of multi-source data based on cesum according to an embodiment of the present invention, and only the relevant parts of the embodiment of the present invention are shown for convenience of description.
The method for realizing the offset-free and deformation-free superposition of multi-source data based on the Cesum comprises the following steps:
step S1, a reference data range of a Gaussian coordinate system is configured in the cecum scene mode.
Firstly, a scene mode is newly added through a scene object in the step, and during specific operation, a localMode attribute is mounted through the scene object Viewer, and the user sets the attribute to indicate that the local mode is started, so that preparation is made for loading Gaussian coordinate system data on a service. After the localhode attribute is set and the local mode is started, the loading function of multi-source data can be realized by the embodiment.
In addition, a global reference data range is required to be configured in the step, and the reference data range is used as a basis for the subsequent translation calculation. In a specific operation, a localExtent attribute may be mounted under the scene object Viewer, which represents a rectangular range including four values: xmin, yminY, xmax, ymax, i.e. minimum X value, minimum Y value, maximum X value, maximum Y value. The value of the attribute is taken from the largest range of multi-source data, namely raster data, vector data, terrain data and model data, for example, the largest range of terrain data, and then the range thereof is taken as the reference data range.
And S2, establishing a conversion interface between the Gaussian coordinate system and the Cartesian coordinate system.
The coordinate frame of the Cesum is a space rectangular coordinate system, belongs to one of Cartesian coordinates, and therefore the input coordinate of the Cesum is the Cartesian coordinate, and the input of the Gaussian coordinate needs to be supported on the premise of loading the Gaussian coordinate coefficient data, so that the Gaussian coordinate system needs to be converted into the Cartesian coordinate system.
The Gaussian coordinate system is converted into the Cartesian coordinate system, and the Cartesian coordinate system is converted into the Cartesian coordinate system.
The interface Cartesian and fromGauss has three parameters, x0, y0 and e, the expression is Cartesian and fromGauss (x 0, y0, e), x0 and y0 are input multi-source data Gaussian coordinates, and e is the reference data range.
The calculation interface processing process is as follows:
201. setting a translation reference point [0, 0], and performing projection calculation on Gaussian coordinates of multi-source data:
x = r * x0
y = r * y0
wherein r is the long radius of the earth, and X0 and Y0 are the gaussian X-coordinate, Y-coordinate of the multi-source data;
202. and (3) calculating the coordinate translation amount:
Tx = (e.xmin+e.xmax) /2–x
Ty = (e.xmin+e.xmax) /2–y
wherein e is a reference data range;
203. calculating the translated coordinates:
Ax = Gx–Tx
Ay = Gy–Ty
204. calculating back projection;
Cx = Ax / r
Cy = Ay / r
205. calculating to obtain Cartesian coordinates:
cartesian coordinates are calculated via the interface cartesian. From radians is an interface for converting radian coordinates into cartesian coordinates, and a ceium system is provided with an interface.
The conversion of Cartesian coordinates into Gaussian coordinates belongs to inverse calculation, an interface Cartesian-ToGauss is added in Cartesian classes of Cartesian classes correspondingly during operation, namely, an inverse calculation interface is obtained, the steps of the inverse calculation interface are opposite to those of the forward calculation interface, namely, the Cartesian coordinates are converted into geographic coordinates firstly, and then the geographic coordinates are converted into the Gaussian coordinates through a longitude and latitude direct projection mode.
The expression for the back-calculation interface is cartesian.
The back calculation interface processing process is as follows:
211. converting the geographic coordinates:
tosradians (c) converts cartesian coordinates to geographic coordinates carto using the interface cartesian provided by ceium;
212. the geographic coordinates carto are subjected to projection calculations:
x' = carto.x * r
y' = carto.y * r
213. calculating Gaussian coordinates:
Gx = x’ + Tx
Gy = y’ + Ty 。
and S3, setting a pyramid loading scheme.
The loading of the grid service and the terrain service depends on a pyramid scheme, such as a common Web mercator projection pyramid scheme, and the map service based on WGS84 coordinates can be loaded.
The pyramid loading scheme set in the step is used for loading the grid data and the terrain data under the Gaussian coordinate system. During specific operation, a pyramid scheme class CustomTiling scheme is firstly stated, a constructor of the pyramid scheme class needs to be transmitted into a scene object Viewer, resolution arrays consisting of each level of resolution of a map slice, and a data range dataExtent of a map service is published. And then realizing projection project and back projection project interfaces, wherein the specific realization is a forward calculation interface and a backward calculation interface which adopt Gaussian coordinates and Cartesian coordinates. And finally recalculating the data range.
The specific operation of recalculating the data range is: firstly, moving a scene to a map central point, namely a place with coordinates of 0, namely calculating the offset by combining a reference data range and the map central point; and carrying out the superposition operation of the offset on the incoming data range to obtain a new data range. In this step, the offset is calculated by using the localExtent attribute in the scene object Viewer in combination with the map center point (i.e., the point with the coordinate of 0), and then the offset is superimposed on the data range dataextend transmitted by the pyramid loading scheme, so that a new data range is obtained. Moving the scene to the place with the map center coordinate of 0 can fully utilize the characteristic of no deformation at the equator to achieve the purpose of restoring the authenticity of the data.
Specifically, the method comprises the following steps:
s31, firstly, calculating the offset by combining the reference data range and the map center point.
The map central point position is also the position of the translated scene, the specific meaning of the map central point position is the intersection point of the meridian and the equator, the map near the point is not deformed, and the purpose of no deformation of the Gaussian scene loading is achieved by utilizing the characteristic. One of the more important concepts is offset, which means the moving distance required for moving the gaussian coordinate scene, and the specific calculation process is as follows:
311. performing projection calculation on the map center:
x" = r * x1
y" = r * y1
wherein X1 and Y1 are the X coordinate and the Y coordinate of the central point;
312. calculating the coordinate translation amount:
Tx" = (e.xmin+e.xmax) /2–x"
Ty" = (e.ymin+e.ymax) /2–y"
where e is the reference data range.
And S32, carrying out the superposition operation of the offset on the incoming data range to obtain a new data range.
When loading terrain data and raster data, a user-defined pyramid scheme object needs to be created, and when the object is created, a data range of specific data needs to be input for calculating the translated data range of the object. The method is to subtract the input data range and the translation amount. The matrix is calculated as follows:
e1".xmin = e1.xmin – Tx"
e1".xmax = e1.xmax – Tx"
e1".ymin = e1.ymin – Ty"
e1".ymax = e1.ymax – Ty"
e1 is the incoming data range, e1 "is the new data range.
And S4, for the received multi-source data to be loaded, if the multi-source data is model data or vector data, preprocessing is carried out.
In the step, only model data and vector data are preprocessed to adapt to a Gaussian scene, and raster data and terrain data are normally issued without special processing.
For model data, two modeling modes are provided for generating vertexes during modeling, wherein one mode is that a modeling central point is located at an original point, and the vertexes of the model data need to be modeled and assigned by using coordinates under Gaussian coordinates; another is to model that the center point is a true gaussian coordinate system coordinate and then the vertex is set to a relative coordinate with respect to that coordinate. The model data modeled by adopting the two modes needs to be processed on the vertex coordinates of the model, namely the vertex coordinates are converted into Cartesian coordinates by adopting a Gaussian coordinate-to-Cartesian coordinate conversion interface so as to recalculate the conversion matrix of the model bounding box and the parent-child nodes.
For vector data, an attribute needs to be added to the vector data, that is, an attribute localExtent is added, and the value of the attribute is the data range of the vector data. The method aims to judge whether data are Gaussian data or not during operation, and if so, the Gaussian coordinates are converted into Cartesian coordinates by means of the set localExtent attribute so as to obtain a correct rendering result.
And S5, selecting and calling a conversion interface or pyramid loading scheme according to the type of the multi-source data to load the data in a Gaussian coordinate system.
The step needs to adopt corresponding schemes for different types of multi-source data. The multi-source data includes raster data, vector data, terrain data, model data, and the like.
In the embodiment, during specific operation, a custom pyramid scheme, customTilligeryProvider class is set for loading raster data, and a constructor of the custom pyramid scheme needs to be transmitted into a custom pyramid scheme, customTillingScheme and a raster service address; setting CustomTinrainProvider class for loading terrain data, wherein the input parameters are custom pyramid scheme CustomTiling Scheme and terrain service address; setting a CustomGeoJsonDataSourceProvider class for loading vector data and needing to transmit a service address of the vector data; the CustomTiesetProvider class is set for loading of model data.
And if the multi-source data is raster data, the Cesium rendering engine obtains a correct map tile data range according to the pyramid loading scheme and initiates a request for loading according to the raster service address.
And if the multi-source data is terrain data, the Cesium rendering engine obtains a correct terrain data range according to the pyramid loading scheme and renders the terrain according to the terrain service address.
If the multi-source data are vector data, a pyramid scheme is not required to be transmitted, and a normal computation interface is called to convert the coordinates of the vector data into Cartesian coordinates according to the set data range. During specific operation, a gauss crsfunction interface is set, which is an interface for converting the coordinates of the existing vector data into cartesian coordinates, and the gauss crsfunction interface transmits the localeExtent attribute set during the previous preprocessing, and then calls a forward calculation interface for converting the coordinates of the vector data into the cartesian coordinates for the use of the procession rendering engine.
If the multi-source data is model data, a relocation interface is set to be used for calculating the real position of the model and positioning without transmitting the model data into a pyramid scheme, a translation value is extracted from a matrix in a requested model description file through specific operation, then the translation value is combined into a Cartesian coordinate, the Cartesian coordinate is converted into a Gaussian coordinate through a back calculation interface, and the real coordinate of the model in a Gaussian coordinate system is obtained. Json, which has an attribute node named transform, and represents the conversion matrix of the model, which is a 4 × 4 fourth-order matrix and is the integration of the model rotation and translation scaling, wherein the first three numbers of the fourth column in the matrix are translation values.
Therefore, the loading operation is executed by the four types of multi-source data through the corresponding definition classes in the step.
In conclusion, because the Cesium can only load three-dimensional data under the WGS84 coordinate system at present, and the 2.5D plane three-dimensional mode of the Cesium is greatly deformed due to the problem of the projection mode, the invention provides the method capable of loading the multi-source three-dimensional data of the Gaussian coordinate system, so that the positions of the data such as grid data, vector data, terrain data and model data after superposition can be ensured to be accurate and have no offset, the deformation of the data after loading can be ensured, and the actual engineering requirements can be met.
The above description is intended to be illustrative of the preferred embodiment of the present invention and should not be taken as limiting the invention, but rather, the intention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention.

Claims (3)

1. A method for realizing non-migration and non-deformation superposition multi-source data based on Cesium is characterized by comprising the following steps:
s1, configuring a reference data range of a Gaussian coordinate system in a ceium scene mode;
s2, establishing a conversion interface between a Gaussian coordinate system and a Cartesian coordinate system;
s3, setting a pyramid loading scheme;
s4, for the received multi-source data to be loaded, if the multi-source data is model data or vector data, preprocessing is carried out;
s5, selecting and calling a conversion interface or a pyramid loading scheme according to the type of multi-source data to load data under a Gaussian coordinate system;
in the step S1, the reference data range represents a rectangular range, and the largest range in the multi-source data is taken, which includes four values, namely, a minimum X value xmin, a minimum Y value ymin, a maximum X value xmax, and a maximum Y value ymax;
in the step S2, the conversion interface comprises a forward calculation interface for converting the Gaussian coordinate system into the Cartesian coordinate system and a backward calculation interface for converting the Cartesian coordinate system into the Gaussian coordinate system; the calculation interface processing process comprises the following steps: firstly, converting Gaussian coordinates into geographic coordinates by a longitude and latitude direct projection mode, and then converting a geographic coordinate system into Cartesian coordinates; the reverse calculation interface processing process comprises the following steps: firstly, converting Cartesian coordinates into geographic coordinates, and then converting the geographic coordinates into Gaussian coordinates in a longitude and latitude direct projection mode;
the forward interface processing process specifically comprises the following steps:
201. setting a translation reference point [0, 0], and performing projection calculation on Gaussian coordinates of multi-source data:
x = r * x0
y = r * y0
wherein r is the long radius of the earth, and X0 and Y0 are the gaussian X-coordinate, Y-coordinate of the multi-source data;
202. and (3) calculating the coordinate translation amount:
Tx = (e.xmin+e.xmax) /2–x
Ty = (e.ymin+e.ymax) /2–y
wherein e is a reference data range;
203. calculating the translated coordinates:
Ax = Gx–Tx
Ay = Gy–Ty
204. back projection calculation:
Cx = Ax / r
Cy = Ay / r
205. calculating to obtain Cartesian coordinates:
cartesian coordinates c are calculated through the interface cartesian.
The reverse calculation interface processing process is as follows:
211. and (3) converting geographic coordinates:
tosradians (c) converts cartesian coordinates to geographic coordinates carto using the interface cartesian provided by ceium;
212. the projection calculation is performed on the geographic coordinates carto:
x’ = carto.x * r
y’ = carto.y * r
213. calculating Gaussian coordinates:
Gx = x’+ Tx
Gy = y’+ Ty
the specific process of the step S3 is as follows:
s31, calculating the offset by combining the reference data range with the map center point:
carrying out projection calculation on the central point of the map:
x" = r * x1
y" = r * y1
wherein X1 and Y1 are the X coordinate and the Y coordinate of the central point;
and (3) calculating the coordinate translation amount:
Tx" = (e.xmin+e.xmax) /2–x"
Ty" = (e.ymin+e.ymax) /2–y"
wherein e is a reference data range;
s32, carrying out the superposition operation of the offset on the incoming data range to obtain a new data range:
e1".xmin = e1.xmin – Tx"
e1".xmax = e1.xmax – Tx"
e1".ymin = e1.ymin – Ty"
e1".ymax = e1.ymax – Ty"
e1 is the incoming reference data range and e1 "is the new data range.
2. The method for realizing shift-free and deformation-free superposition of multi-source data based on Cesium as claimed in claim 1, wherein in step S4, model data are preprocessed as follows: converting the model data modeling vertex coordinates into Cartesian coordinates by adopting a forward calculation interface; the vector data is preprocessed as follows: a reference data range is set for the vector data.
3. The method for realizing non-migration and non-deformation superposition of multi-source data based on the Cesum as claimed in claim 2, wherein in the step S5, if the multi-source data is raster data, a correct map tile data range is obtained according to a pyramid loading scheme, and a request is initiated according to a raster service address for loading; if the multi-source data is terrain data, obtaining a correct terrain data range according to the pyramid loading scheme, and rendering the terrain according to a terrain service address; if the multi-source data are vector data, calling a forward calculation interface to convert the coordinates of the vector data into Cartesian coordinates according to the set data range; and if the multi-source data is model data, extracting translation values from the matrix in the requested model description file, combining the translation values into a Cartesian coordinate, and converting the Cartesian coordinate into a Gaussian coordinate through a back calculation interface to obtain the real coordinate of the model in the Gaussian coordinate system.
CN202211037107.4A 2022-08-29 2022-08-29 Method for realizing non-offset and non-deformation superposition of multi-source data based on Cesum Active CN115129401B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211037107.4A CN115129401B (en) 2022-08-29 2022-08-29 Method for realizing non-offset and non-deformation superposition of multi-source data based on Cesum

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211037107.4A CN115129401B (en) 2022-08-29 2022-08-29 Method for realizing non-offset and non-deformation superposition of multi-source data based on Cesum

Publications (2)

Publication Number Publication Date
CN115129401A CN115129401A (en) 2022-09-30
CN115129401B true CN115129401B (en) 2022-11-15

Family

ID=83387238

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211037107.4A Active CN115129401B (en) 2022-08-29 2022-08-29 Method for realizing non-offset and non-deformation superposition of multi-source data based on Cesum

Country Status (1)

Country Link
CN (1) CN115129401B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111402390A (en) * 2020-02-20 2020-07-10 平安科技(深圳)有限公司 Model rendering method, device, equipment and storage medium
CN112256897A (en) * 2020-11-04 2021-01-22 重庆市地理信息和遥感应用中心 Vector tile loading method in three-dimensional scene
CN113192193A (en) * 2021-04-23 2021-07-30 安徽省皖北煤电集团有限责任公司 High-voltage transmission line corridor three-dimensional reconstruction method based on Cesium three-dimensional earth frame
CN113378367A (en) * 2021-06-03 2021-09-10 中国人民解放军陆军防化学院 Chemical simulation environment construction method based on complex terrain

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101532834A (en) * 2008-12-12 2009-09-16 北京理工大学 Method for switching coordinates of WGS84 ellipsoid and Clarke 80 ellipsoid
CN107003122A (en) * 2014-12-09 2017-08-01 巴斯夫欧洲公司 Fluorescence detector
US10559079B2 (en) * 2017-05-19 2020-02-11 Uih America, Inc. System and method for image reconstruction
CN109612445B (en) * 2018-12-17 2021-04-30 中国水利水电第十四工程局有限公司 High-precision terrain establishing method under WebGIS platform based on unmanned aerial vehicle
CN110851558B (en) * 2020-01-14 2020-05-19 武大吉奥信息技术有限公司 Universal loading method for map of coordinate system in earth
CN112489210A (en) * 2020-11-12 2021-03-12 广东省国土资源技术中心(广东省基础地理信息中心) Method for constructing autonomous controllable three-dimensional natural resource map system
CN112184798B (en) * 2020-11-27 2021-05-04 武大吉奥信息技术有限公司 WebGL three-dimensional-based surface area calculation method and device
CN114756937B (en) * 2022-04-21 2024-03-12 四川蓉电科技发展有限公司 Visualization system and method based on UE4 engine and Cesium framework

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111402390A (en) * 2020-02-20 2020-07-10 平安科技(深圳)有限公司 Model rendering method, device, equipment and storage medium
CN112256897A (en) * 2020-11-04 2021-01-22 重庆市地理信息和遥感应用中心 Vector tile loading method in three-dimensional scene
CN113192193A (en) * 2021-04-23 2021-07-30 安徽省皖北煤电集团有限责任公司 High-voltage transmission line corridor three-dimensional reconstruction method based on Cesium three-dimensional earth frame
CN113378367A (en) * 2021-06-03 2021-09-10 中国人民解放军陆军防化学院 Chemical simulation environment construction method based on complex terrain

Also Published As

Publication number Publication date
CN115129401A (en) 2022-09-30

Similar Documents

Publication Publication Date Title
CN108648269B (en) Method and system for singulating three-dimensional building models
US20230053462A1 (en) Image rendering method and apparatus, device, medium, and computer program product
WO2021164233A1 (en) Model rendering method, apparatus, and device, and computer readable storage medium
CN108230435B (en) Graphics processing using cube map textures
CN106990961B (en) Method for establishing WebGL graphic rendering engine
CN112102489A (en) Navigation interface display method and device, computing equipment and storage medium
CN115546371A (en) Point cloud optimization method and system, electronic device and storage medium
CN114419241A (en) Three-dimensional model construction method and device and storage medium
CN115129401B (en) Method for realizing non-offset and non-deformation superposition of multi-source data based on Cesum
CN104517313A (en) AO (ambient occlusion) method based on screen space
CN115035233B (en) Global multi-resolution space grid spherical ray projection visualization method and device
CN115619986B (en) Scene roaming method, device, equipment and medium
CN111243089B (en) Method, device and system for drawing three-dimensional island model based on island two-dimensional data
CN111861873B (en) Method and device for generating simulation image
CN109658495B (en) Rendering method and device for ambient light shielding effect and electronic equipment
CN112308778B (en) Method and terminal for assisting panoramic camera splicing by utilizing spatial three-dimensional information
US6188409B1 (en) 3D graphics device
CN106384372A (en) View synthesis method and device
CN113256771A (en) Graphics system and graphics processing method therefor
CN115774896B (en) Data simulation method, device, equipment and storage medium
Zhou et al. Object detection and spatial location method for monocular camera based on 3D virtual geographical scene
CN113674435A (en) Image processing method, electronic map display method and device and electronic equipment
CN109191556B (en) Method for extracting rasterized digital elevation model from LOD paging surface texture model
KR102077494B1 (en) Method and apparatus for providing around view
CN116992640A (en) Two-stage submarine topography scouring dynamic simulation method based on three.js

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
GR01 Patent grant
GR01 Patent grant