CN114936261B - Plane coordinate conversion method for power grid homologous maintenance tool and GIS graph migration - Google Patents
Plane coordinate conversion method for power grid homologous maintenance tool and GIS graph migration Download PDFInfo
- Publication number
- CN114936261B CN114936261B CN202210707929.2A CN202210707929A CN114936261B CN 114936261 B CN114936261 B CN 114936261B CN 202210707929 A CN202210707929 A CN 202210707929A CN 114936261 B CN114936261 B CN 114936261B
- Authority
- CN
- China
- Prior art keywords
- coordinates
- gis
- coordinate
- cgsc2000
- plane
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000012423 maintenance Methods 0.000 title claims abstract description 70
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000013508 migration Methods 0.000 title claims abstract description 18
- 230000005012 migration Effects 0.000 title claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 230000009466 transformation Effects 0.000 claims abstract description 25
- 230000009191 jumping Effects 0.000 claims description 10
- 230000001131 transforming effect Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 238000013519 translation Methods 0.000 claims description 3
- 230000002457 bidirectional effect Effects 0.000 abstract description 6
- 238000012937 correction Methods 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 4
- 230000002411 adverse Effects 0.000 abstract description 3
- 238000004422 calculation algorithm Methods 0.000 description 20
- 238000010586 diagram Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000013507 mapping Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Operations Research (AREA)
- Remote Sensing (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Graphics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Instructional Devices (AREA)
Abstract
The invention discloses a plane coordinate conversion method for power grid homologous maintenance tool and GIS graph migration, belonging to the technical field of calculation, calculation or counting. For GIS graphic data to be migrated, according to the equipment type described by the GIS graphic data, carrying out projection transformation on GIS coordinates of equipment in a power distribution network station through a grid to obtain local plane coordinates, revising the double-precision value of the local plane coordinates after projection transformation, eliminating adverse effects of projection transformation deformation, ensuring coordinate precision reserved in an original WGS-84 system after newly generated graphic data of a homologous maintenance tool is written back to the original GIS system, and realizing smooth migration of bidirectional graphic data between the grid GIS system and the homologous maintenance tool on the distribution network side without being influenced by grid alignment correction operation.
Description
Technical Field
The invention relates to an electric power automation technology, in particular to a plane coordinate conversion method for power grid homologous maintenance tools and GIS graph migration, and belongs to the technical field of calculation, calculation or counting.
Background
In the construction of a power grid resource middle station, the homologous maintenance tool needs to migrate stock data in the existing system first, and switch on line after the track is operated for a period of time, and in order not to influence the existing production system, incremental data in the homologous maintenance tool also needs to be written back to the original system during the track-parallel operation, and the bidirectional data synchronous migration relates to a basic GIS graph part. The national grid geographic information system (Geographic Information System, GIS) platform adopts unified WGS-84 geodetic coordinates for the power transmission network and distribution network service of a plurality of provinces and cities of the whole country, the homologous maintenance tool adopts urban independent plane coordinates when facing the urban distribution network service, the latest urban independent plane coordinates adopt CGCS2000 geodetic coordinates, and the CGCS2000 geodetic coordinates and the GIS geodetic coordinates have different reference ellipsoid parameters, so that the geodetic coordinates are converted into plane coordinates through a projection algorithm when migrating the GIS graph to the homologous maintenance tool, and projection deformation inevitably exists after the projection transformation operation. The deformation analysis of the commonly used gaussian-kriging projections is as follows: (1) The central meridian has no deformation, and the condition of unchanged length ratio after projection is satisfied; (2) The length ratio of any other point is larger than 1 except the length ratio of the central meridian is 1; (3) On the same weft, the farther from the central warp, the larger the deformation is, and the maximum value is positioned at the edge of the projection belt; (4) On the same meridian, the lower the latitude is, the larger the deformation is, and the maximum value is positioned on the equator; (5) The equiangular projection has no angular deformation, and the area ratio is the square of the length ratio; (6) The length-ratio isomorphous lines are parallel to the central axis meridian.
The geographical edge layout of the distribution system drawn by the power grid GIS platform comprises transformer substations, distribution stations, transmission lines, distribution lines, on-pole equipment and in-station equipment, wherein the power station equipment is not abstracted into simple points, the power station equipment is drawn into rectangular frames near site sites, internal electric wiring of the power station equipment is unfolded in the rectangular frames, and the equipment on the outside of the power station is distributed along the geographical direction. Non-uniform distortion and deformation exist after direct projection conversion of the in-station line drawing, so that the flat and vertical tidy style of the original picture surface is destroyed, and the projection deformation is unacceptable. Therefore, the coordinate grid alignment algorithm of the equipment in the power station needs to be expanded on the basis of the basic geographic coordinate conversion algorithm and the flow, and the electric wiring diagram in the power station is ensured not to be deformed after the homologous maintenance tool is led into the power station equipment for drawing. In addition, the newly generated graphic data of the homologous maintenance tool may need to be written back to the original GIS system, and how to ensure the accuracy of the write-back coordinates is also a problem that needs to be considered in the bidirectional smooth migration of the graphic between the homologous maintenance tool of the power distribution network and the GIS of the power distribution network.
In summary, the invention aims to provide a plane coordinate conversion method for graph migration of a power grid homologous maintenance tool and a GIS so as to realize the purpose of bi-directional seamless graph data migration between the power grid homologous maintenance tool and the GIS.
Disclosure of Invention
The invention aims to solve the problems in the background, and provides a plane coordinate conversion method for migrating a power grid homologous maintenance tool and a GIS graph, which expands a coordinate grid alignment algorithm of equipment in a station on the basis of a traditional geographic coordinate conversion algorithm and a traditional process, so that an electric wiring diagram in the station is not deformed after the homologous maintenance tool is led in, and simultaneously, by revising a plane coordinate double-precision value and adding necessary grid revision parameters, the original coordinate precision can be ensured to be restored when graphic data are written back to an original power grid GIS system, the adverse effect caused by projection deformation is eliminated, the graphic data is migrated in a two-way seamless manner between the homologous maintenance tool and the power grid GIS graph system, and the technical problem of deformation after the wiring diagram in the station of the traditional GIS system is migrated to the power grid homologous maintenance tool is solved.
The invention adopts the following technical scheme for realizing the purposes of the invention:
The coordinate conversion method for the power grid homologous maintenance tool and GIS graph migration comprises two parts: the GIS graph is transferred to the coordinate conversion method of the power grid homologous maintenance tool, and the graph generated by the power grid homologous maintenance tool is written back to the coordinate conversion method of the GIS.
The coordinate conversion method for migrating the GIS graph to the power grid homologous maintenance tool comprises the following 7 steps:
Step 1, judging the equipment type of GIS graphic description to be migrated: if the equipment to be migrated of the GIS graphic description is off-site equipment of the power distribution website, skipping to the step 2; if the equipment to be migrated of the GIS graphic description is power distribution station equipment, skipping to the step 3; if the equipment to be migrated of the GIS graphic description is equipment in the power distribution station, skipping to the step 4; if the equipment to be migrated of the GIS graphic description is bridging equipment, skipping to the step 5;
step 2, performing geographic coordinate conversion on the GIS graphic coordinates, converting the GIS graphic coordinates into CGSC2000 plane coordinates, and jumping to the step 7;
Step 3, acquiring angular point GIS graphic coordinates describing a station frame of the power distribution station equipment, acquiring GIS grid parameters, acquiring the size of a grid unit after projection transformation according to CGSC2000 plane coordinates corresponding to the angular point GIS graphic coordinates, and jumping to step 6;
step 4, projectively transforming the GIS graphic coordinates to the CGSC2000 plane coordinate system in a mode of aligning the GIS grids and the CGSC2000 plane coordinate system grids, and jumping to the step 6;
Step 5, if the GIS graphic coordinates fall into the station frame, jumping to the step 4; otherwise, jumping to the step 2;
step 6, revising a coordinate double-precision value format, and adding grid parameters to the CGSC2000 plane coordinate after the projection transformation in the step 4 to obtain a power distribution network homologous maintenance graphic coordinate;
And 7, outputting the homologous maintenance graphic coordinates of the power distribution network.
The specific method for converting the geographic coordinates of the GIS graphic coordinates in the step 2 is as follows:
Step 2-1, converting the GIS graphic coordinates (B D,LD,HD) into WGS-84 ellipsoid space rectangular coordinates (X D,YD,ZD), ND、/>A radius of curvature of a circle of mortise for a WGS84 ellipsoid, a first eccentricity;
Step 2-2, converting WGS-84 ellipsoid space rectangular coordinate (X D,YD,ZD) into CGSC2000 ellipsoid space rectangular coordinate (X C,YC,ZC) by using a Boolean seven-parameter conversion formula, T x,Ty,Tz is a 3-axis translation parameter, epsilon x,εy,εz is a 3-axis rotation parameter, and m is a scale parameter;
step 2-3, converting the space rectangular coordinate (X C,YC,ZC) of the CGSC2000 ellipsoid into the geodetic coordinate (B C,LC,HC) of the CGSC2000 ellipsoid, And (3) calculating B in an iterative mode until |B i+1-Bi | < epsilon, wherein the calculation formula of an initial value B 0 of B is as follows: /(I)Wherein, a C、NC,/>A long half shaft and a mortise unitary circle curvature radius which are CGSC2000 ellipsoids, and a first eccentricity;
Step 2-4, carrying out Gaussian-Kelvin projection forward calculation on the CGSC2000 ellipsoidal geodetic coordinate (B C,LC,HC) to obtain a CGSC2000 plane coordinate (x, y),
Wherein, t=tan B C,L=l C-L0,L0 is the central meridian, L 0 is the central meridian longitude, and X is the central meridian arc length.
In the step 3, GIS graphic coordinates (B 0,L0)、(B1,L1) describing the lower left corner and the upper right corner of a station frame of the power distribution station equipment are read, and the GIS transverse grid number g n,gn＝(L1-L0)/g is obtained; calculating the size of grid cells after projection transformation according to planar coordinates (x 0,y0) and (x 1,y1) after projection transformation of GIS graphic coordinates describing the lower left corner and the upper right corner of a station frame of power distribution station equipment
In step 4, the method for projectively transforming the GIS graphic coordinates to the CGSC2000 plane coordinate system by aligning the GIS grid and the CGSC2000 plane coordinate system grid comprises the following steps:
Calculating the coordinates (x g,yg) of the equipment GIS in the power distribution station relative to the reference point, i.e. relative to the lower left corner of the station frame,
Translating the GIS coordinates (B, L) of the equipment in the power distribution station relative to the coordinates (x g,yg) of the reference point according to the size of the grid unit after projection transformation to obtain CGSC2000 plane coordinates (x, y) after projection transformation,
In step 6, revising the coordinate double-precision value format, and adding grid parameters to the CGSC2000 plane coordinate after the projection transformation in step 4 comprises the following specific methods:
For the CGSC2000 plane coordinate (x, y) obtained in the step 4, intercepting the integer part of the plane coordinate and the value of three bits after the decimal point as an original coordinate value, and taking the data of 3 bits after the decimal point as a millimeter-level deviation value;
Comparing the CGSC2000 plane coordinate obtained in the step 4 with the plane coordinate obtained by directly converting the GIS graphic coordinate: when the two are free from deviation, the original coordinate value obtained by intercepting the CGSC2000 plane coordinate obtained in the step 4 is the final homologous maintenance graphic coordinate of the power distribution network; when the two have deviation, a millimeter-level deviation value needs to be added to the original coordinate value, and the CGSC2000 plane coordinate after the millimeter-level deviation value is added is the final homologous maintenance graphic coordinate of the power distribution network.
A coordinate conversion method for writing back a graph generated by a power distribution network homologous maintenance tool to a GIS comprises the steps of A to C.
Step A, newly generating plane coordinates of a graph for a power distribution network maintenance tool, and taking a integer value according to millimeter precision;
Step B, calculating a millimeter-level deviation value, correcting the plane coordinates of the newly generated graph of the power distribution network maintenance tool according to the millimeter-level deviation value, accumulating the millimeter-level deviation value for the plane coordinates of the newly generated graph of the power distribution network maintenance tool when the absolute value of the millimeter-level deviation value is larger than 1, and not processing the plane coordinates of the newly generated graph of the power distribution network maintenance tool when the absolute value of the millimeter-level deviation value is not larger than 1;
and C, carrying out coordinate conversion on the plane coordinates of the newly generated graph of the corrected power distribution network maintenance tool, carrying out Gaussian-Kelvin back calculation on the plane coordinates of the newly generated graph of the power distribution network maintenance tool to obtain CGSC2000 ellipsoidal geodetic coordinates, converting the CGSC2000 ellipsoidal geodetic coordinates into CGSC2000 ellipsoidal space rectangular coordinates, converting the CGSC2000 ellipsoidal space rectangular coordinates into WGS-84 ellipsoidal space rectangular coordinates by adopting a Boolean seven-parameter formula, and converting the WGS-84 ellipsoidal space rectangular coordinates into WGS-84 geodetic coordinates which are written back to a GIS system, wherein the whole conversion process is the inverse process of geographic coordinate conversion on the GIS graph coordinates.
In the step C, the specific method for carrying out Gaussian-Kelvin back calculation on the plane coordinates of the newly generated graph of the power distribution network maintenance tool comprises the following steps:
step C-1, iteratively solving the bottom point latitude B f, and calculating an initial value of the bottom point latitude B f: Up to/>
And C-2, back calculating GIS coordinates (B, L) corresponding to the plane coordinates (x, y) of the newly generated graph of the power distribution network maintenance tool according to the following formula:
L＝l+L0，
wherein, t f＝tan Bf is a total number of the components, aD、/>The long half shaft, the first eccentricity and the second eccentricity of the WGS-84 ellipsoid.
The invention adopts the technical scheme and has the following beneficial effects: aiming at the characteristics of the geographical edge layout of the power grid, the invention expands the coordinate grid alignment algorithm of the equipment in the station on the basis of the traditional geographic coordinate conversion algorithm and flow, so that the electric wiring diagram in the station is not deformed after the homologous maintenance tool is imported, and meanwhile, by revising the plane coordinate double-precision value and adding necessary grid revision parameters, the original coordinate precision can be ensured to be restored when the graphic data is written back to the original power grid GIS system, the adverse effect caused by projection deformation is eliminated, and the graphic data is transferred between the homologous maintenance tool and the power grid GIS graphic system in a bidirectional seamless manner.
Drawings
Fig. 1 is a flowchart of a planar coordinate conversion method for migrating a power grid homologous maintenance tool and a GIS graph.
Fig. 2 is a flowchart for calculating seven parameters of the boolean sa model.
FIG. 3 is a geographic coordinate transformation algorithm and flow.
Fig. 4 is a flowchart of a grid GIS WGS-84 longitude and latitude coordinate conversion homologous maintenance tool CGSC2000 plane coordinate.
Fig. 5 is a flowchart of the homologous maintenance tool CGSC2000 plane coordinate conversion grid GIS system WGS-84 latitude and longitude coordinates.
Detailed Description
The technical scheme of the invention is described in detail below with reference to the accompanying drawings.
The application provides a plane coordinate conversion method for power grid homologous maintenance tool and GIS graph migration, which comprises the following four steps as shown in figure 1: firstly, calculating seven parameters of a Boolean sha model through coincident points provided by a mapping department, and solving a reference plane conversion formula between a WGS-84 ellipsoid and a CSCG2000 ellipsoid; secondly, synthesizing a Gaussian-Kelvin projection forward and backward calculation formula, a geodetic coordinate and space rectangular coordinate interconversion formula and a reference plane conversion formula, and determining an algorithm and a flow of geographic coordinate conversion; thirdly, determining a flow of a grid GIS system WGS-84 longitude and latitude coordinate conversion homologous maintenance tool CGSC2000 plane coordinate: the grid alignment algorithm of the equipment coordinates in the station uses the lower left corner of the frame of the station before projection as a reference point coordinate, calculates the coordinate relative to the reference point, divides the size of the grid cells of the station to obtain grid relative coordinates, and calculates new plane coordinates based on the grid relative coordinates, the coordinates of the projected reference point and the mapped grid size; respectively converting and calculating coordinates of inter-station and inter-station cross-border equipment; revising the plane coordinate double-precision value and adding grid revision parameters; and fourthly, determining a flow of converting the grid GIS system WGS-84 longitude and latitude coordinates by the homologous maintenance tool CGSC2000 plane coordinates, and realizing the reverse restoration calculation of the plane coordinates to obtain the geodetic coordinates according to the double-precision revised value of the plane coordinates.
Step one, determining an ellipsoidal reference surface conversion parameter
According to three or more coincident point simultaneous observation equations, the seven parameters of the Boolean sha model are calculated through a least square method, the calculation flow is shown in figure 2, and the related parameters and formulas are as follows:
1) Ellipsoid parameters
2) Boolean sha model
The boolean sevof parameter model is a three-dimensional coordinate transformation model in a space rectangular coordinate system, and has a strict transformation relationship between two space rectangular coordinate systems, and comprises 7 transformation parameters, namely 3 translation parameters T x,Ty,Tz, 3 rotation parameters epsilon x,εy,εz and 1 scale parameter m. Let (X D,YD,ZD) denote the source coordinates and (X C,YC,ZC) denote the destination coordinates, then the directly available coordinate transformation equation is as follows:
XC＝XD+Tx-ZD*εy+YD*εz+XD*m
YC＝YD+Ty+ZD*εx-XD*εz+YD*m
ZC＝ZD+Tz-YD*εx+XD*εy+ZD*m
the Boolean seven-parameter conversion formula can be obtained by arranging according to a linear equation:
When the WGS-84 longitude and latitude coordinates are source coordinates, the CGSC2000 plane coordinates are target coordinates; when the CGSC2000 plane coordinate is a source coordinate, the WGS-84 longitude and latitude coordinate is a target coordinate, and the Boolean seven-parameter conversion formula can be used for bidirectional conversion between the space rectangular coordinate under the WGS-84 longitude and latitude coordinate system and the space rectangular coordinate under the CGSC2000 plane coordinate system. The space rectangular coordinates under the WGS-84 longitude and latitude coordinate system and the CGSC2000 plane coordinate system can be directly obtained or obtained by conversion according to the WGS-84 longitude and latitude coordinate (B, L) of the coincident point and the coordinates (x, y) of the CGSC2000 plane coordinate system.
3) Observation equation
Assuming that the WGS-84 ellipsoid space rectangular coordinate system and the CGSC2000 ellipsoid space rectangular coordinate system have n coincident points, and 7 parameters correspond to 7 necessary observation numbers, at least 3 points are needed to form a sufficient condition to complete equation solution, and an error equation matrix form can be listed as follows:
the simplified form is: l is the observed value, V is the correction number,/> Is an independent parameter (7 parameters in the present invention), and B is an observation matrix.
4) Adjustment calculation
The least squares solution can be used to:
the coordinates of each point are generally regarded as the independent observation value with the same precision, and then the weight matrix P can take an identity matrix to obtain a simplified formula:
Step two: determining geographic coordinate transformation algorithm and process
The bidirectional conversion flow between the WGS-84 longitude and latitude coordinates and the local plane coordinates of the power grid homologous maintenance tool is shown in fig. 3, wherein the local plane coordinates refer to the local plane coordinates of the power grid homologous maintenance tool, and are also plane coordinates of the CGSC2000 ellipsoidal geodetic coordinate system projected by Gauss-gram Lv Gezheng. The Boolean seven parameter formula can realize the conversion between the space rectangular coordinates of different ellipsoid references, and the process of ' WGS84 space rectangular coordinates ' is ' way, CGSC2000 space rectangular coordinates ' is ' way, and the mutual conversion between the space rectangular coordinates of different ellipsoid references is realized, so that the bridging effect is realized; the Gaussian-Gauss projection forward and backward calculation formula can realize the mutual conversion between ' the geodetic coordinates of the reference of the designated ellipsoid ' and the local plane coordinates '.
The correlation formula is as follows:
1) The geodetic coordinates BLH are converted into spatial rectangular coordinates XYZ (BLH→XYZ)
X＝(N+H)cos B cos L
Y＝(N+H)cos B sin L
Z＝[N(1-e2)+H]cos B sin L
Wherein N is the radius of curvature of the unitary mortise ring,E is the first eccentricity of the ellipsoid, a is the major half axis of the ellipsoid, B is the minor half axis of the ellipsoid, B is the latitude, L is the longitude, H is the geodetic altitude, and (X, Y, Z) is the space rectangular coordinate.
2) Conversion of space rectangular coordinate XYZ into geodetic coordinate BLH (XYZ→BLH)
The calculation formula for calculating the initial value B 0 of B in an iterative mode is as follows: The iteration is repeated until |b i+1-Bi | < epsilon, which is typically 2.78 x 10-9 degrees, i.e. 0.00001 seconds.
3) Gauss-Gauss projection forward calculation formula
Wherein: l 0 represents the central meridian longitude, t=tan B, η=e' 2cos2B,l＝L-L0, X represents the central meridian arc length, the calculation formula of which is:
A 1、A2、A3、A4 is a constant, and the calculation formula is as follows:
4) Gauss-Gauss projection back calculation formula
L＝l+L0
B f is the base point latitude, and when x=x, the meridian length corresponds to the latitude. Wherein the angle sign f is a function value related to the latitude of the bottom point, t f＝tan Bf,ηf＝e′2cos2Bf,
The bottom point latitude B f is solved by adopting an iteration method, and the method comprises the following steps:
(1) Calculating an initial value:
(2) Each iteration:
(3) Repeating the iteration until Epsilon is generally 2.78 x 10-9 degrees, i.e. 0.00001 seconds
Step three, determining a grid GIS system WGS-84 longitude and latitude coordinate conversion homologous maintenance tool CGSC2000 plane coordinate flow
As shown in fig. 4, the flow of the grid GIS WGS-84 longitude and latitude coordinate conversion homologous maintenance tool CGSC2000 plane coordinate includes the following 3 steps:
1) Grid alignment algorithm for determining equipment coordinates in a station
Firstly, selecting the lower left corner of a rectangular station frame of power station equipment as a reference point coordinate, calculating a coordinate relative to the reference point, dividing the coordinate by the size of a grid cell of the station equipment to obtain a grid relative coordinate, and then calculating a new plane coordinate based on the grid relative coordinate, the coordinate projected by the reference point and the mapped grid size.
Let the lower left corner coordinate of the station frame be (B 0,L0), the upper right corner coordinate be (B 1,L1), the grid cell size be g in the GIS system, then the transverse grid number be: g n＝(L1-L0)/g;
Let the coordinates after projection transformation of the lower left corner and the upper right corner of the station frame be (x 0,y0) and (x 1,y1), respectively, then the size of the grid cell after mapping is:
If the coordinates of the equipment points in the station are (B, L) and the lower left corner is taken as a reference point, the relative coordinates of the grids are as follows: The coordinates after projection transformation are: /(I)
2) Coordinate conversion calculation for in-station and out-of-station bridging equipment
The equipment outside the station directly converts according to a geographic coordinate conversion algorithm; the power station equipment reserves the size parameters of a rectangular station frame and a grid unit; calculating projected coordinates of the equipment coordinates in the station according to a grid alignment algorithm; the internal and external bridging equipment of the power station only needs to calculate the endpoint coordinates falling into the power station according to a grid alignment algorithm, and other point coordinates are directly converted according to geographic coordinates.
The rectangular station frame calculation algorithm of the power station equipment is as follows:
The lower left corner coordinate of the station frame is set as (B 0,L0), the upper right corner coordinate is set as (B 1,L1), the lower left corner plane coordinate (x 0,y0) can be obtained by direct conversion according to the geographic coordinates, and the rest points are calculated according to a grid alignment algorithm as endpoints falling into the power station.
3) Revising plane coordinate double-precision value
And the necessary grid correction parameters are added by revising the format of the plane coordinate double-precision numerical value, so that the original precision of the coordinate is ensured to be recovered during the write-back.
The 64-bit double-precision numerical value consists of a 1-bit sign bit, an 11-bit exponent part and a 52-bit mantissa part from high to low, and 15 or 16 decimal significant digits are at most used. The range of the urban east-west or north-south is not more than 450km, if the geometric center of the city is taken as the origin and the unit is taken as the meter, the plane coordinate value is shown as +/-123456.1234567890, 9-10 valid digits are arranged after the decimal point, and if the digits are accurate to the millimeter (3 digits after the decimal point), only 9 valid digits are occupied, and redundant 6-7 valid digits are needed.
Considering actual mapping achievements and power application requirements in China, the plane coordinate value is enough to be reserved to millimeter precision. Through calculation, the difference between the plane coordinate obtained by the grid alignment algorithm in the GIS system and the plane coordinate obtained after direct conversion is less than 10 meters, namely the deviation value is less than or equal to 9999 millimeters, so that the plane coordinate value can be divided into 2 parts to be expressed as XXX, XXX.XXXYYYYYY, for example: 318445.313009223, wherein the first 9 bits 318445.313 are the original coordinate values and the last 6 bits 009223 are used to represent the coordinate deviation values of 9233 mm that need to be corrected.
The decimal point is firstly intercepted to millimeter unit, such as: 318445.313, -318445.313;
When the projected transformed coordinates obtained by the grid alignment algorithm have no deviation from plane coordinates converted from direct coordinates, the coordinates are directly intercepted without any processing;
when the projected transformed coordinates obtained by the grid alignment algorithm have deviation values with plane coordinates converted from direct coordinates, namely, when the coordinates in the station need to be corrected due to grid alignment, a deviation part is added on the basis of intercepting the coordinates, such as:
318445.313+9223÷1000,000,000＝318445.313009220，
318445.313-9223÷1000,000,000＝318445.312990777，
-318445.313+9223÷1000,000,000＝-318445.312990777，
-318445.313-9223÷1000,000,000＝-318445.313009220。
After the deviation value is added, the numerical value is slightly changed, and the millimeter precision is unchanged.
Step four, determining a process flow of determining the longitude and latitude coordinates of a WGS-84 of a grid GIS system of the CGSC2000 plane coordinate conversion of the homologous maintenance tool
Because whether the unknown plane coordinate is revised or not, firstly, a integer value, such as-318445.312990777 meters, is obtained after millimeter rounding, and corresponding meter unit coordinate value-318445.313 is equivalent to the intercepting coordinate in the step three; then, calculating the deviation part (10 9 times of the obtained coordinates are amplified after the obtained coordinates are removed), and [ (318445.312990777- (-318445.313) ] multiplied by 1000,000,000 =9223, namely, the last 6 digits in the 9 digits after the decimal point of the original coordinates, if the absolute value of the deviation part is larger than 1, indicating that an additional correction value is needed, -318445.313+9223 ≡1000= -318436.090, otherwise, indicating that no additional correction value is needed, obtaining the obtained coordinates, and obtaining the needed coordinates without processing. The reverse recovery calculation flow is shown in fig. 5.
The foregoing description of the preferred embodiment of the invention is merely illustrative of the invention and is not intended to be limiting. Many changes, modifications and equivalents will occur to those skilled in the art and it is intended to cover all such changes, modifications and equivalents as fall within the true spirit and scope of the invention as defined by the appended claims.
Claims (7)
1. The planar coordinate conversion method for the power grid homologous maintenance tool and GIS graph migration is characterized by comprising a coordinate conversion process for obtaining CGSC2000 planar coordinates by carrying out coordinate conversion on GIS graph data to be migrated, wherein the coordinate conversion process for obtaining the CGSC2000 planar coordinates by carrying out coordinate conversion on the GIS graph data to be migrated specifically comprises the following steps:
Step 1, judging the equipment type of GIS graphic description to be migrated: if the equipment to be migrated of the GIS graphic description is off-site equipment of the power distribution website, skipping to the step 2; if the equipment to be migrated of the GIS graphic description is power distribution station equipment, skipping to the step 3; if the equipment to be migrated of the GIS graphic description is equipment in the power distribution station, skipping to the step 4; if the equipment to be migrated of the GIS graphic description is bridging equipment, skipping to the step 5;
step 2, performing geographic coordinate conversion on the GIS graphic coordinates, converting the GIS graphic coordinates into CGSC2000 plane coordinates, and jumping to the step 7;
Step 3, acquiring angular point GIS graphic coordinates describing a station frame of the power distribution station equipment, acquiring GIS grid parameters, acquiring the size of a grid unit after projection transformation according to CGSC2000 plane coordinates corresponding to the angular point GIS graphic coordinates, and jumping to step 6;
step 4, projectively transforming the GIS graphic coordinates to the CGSC2000 plane coordinate system in a mode of aligning the GIS grids and the CGSC2000 plane coordinate system grids, and jumping to the step 6;
Step 5, if the GIS graphic coordinates fall into the station frame, jumping to the step 4; otherwise, jumping to the step 2;
step 6, revising a coordinate double-precision value format, and adding grid parameters to the CGSC2000 plane coordinate after the projection transformation in the step 4 to obtain a power distribution network homologous maintenance graphic coordinate;
And 7, outputting the homologous maintenance graphic coordinates of the power distribution network.
2. The method for converting plane coordinates of power grid homologous maintenance tool and GIS graphic migration according to claim 1, wherein the specific method for converting geographic coordinates of GIS graphic coordinates in step 2 is as follows:
Step 2-1, converting the GIS graphic coordinates (B D,LD,HD) into WGS-84 ellipsoid space rectangular coordinates (X D,YD,ZD), ND、/>A radius of curvature of a circle of mortise for a WGS84 ellipsoid, a first eccentricity;
Step 2-2, converting WGS-84 ellipsoid space rectangular coordinate (X D,YD,ZD) into CGSC2000 ellipsoid space rectangular coordinate (X C,YC,ZC) by using a Boolean seven-parameter conversion formula, T x,Ty,Tz is a 3-axis translation parameter, epsilon x,εy,εz is a 3-axis rotation parameter, and m is a scale parameter;
step 2-3, converting the space rectangular coordinate (X C,YC,ZC) of the CGSC2000 ellipsoid into the geodetic coordinate (B C,LC,HC) of the CGSC2000 ellipsoid, And (3) calculating B in an iterative mode until |B i+1-Bi | < epsilon, wherein the calculation formula of an initial value B 0 of B is as follows: /(I)Wherein, a C、NC,/>A long half shaft and a mortise unitary circle curvature radius which are CGSC2000 ellipsoids, and a first eccentricity;
Step 2-4, carrying out Gaussian-Kelvin projection forward calculation on the CGSC2000 ellipsoidal geodetic coordinate (B C,LC,HC) to obtain a CGSC2000 plane coordinate (x, y),
Wherein, t= tanB C,L=l C-L0,L0 is the central meridian, L 0 is the central meridian longitude, and X is the central meridian arc length.
3. The method for converting planar coordinates of power grid homology maintenance tool and GIS graphic migration according to claim 1, wherein the GIS graphic coordinates of the corner points of the station frame describing the power distribution station equipment obtained in the step 3 are GIS graphic coordinates of the lower left corner of the station frame (B 0,L0) and GIS graphic coordinates of the upper right corner of the station frame (B 1,L1), the obtained GIS grid parameters are GIS transverse grid numbers g n,gn＝(L1-L0)/g, and the size of the grid unit after projective transformation is the size of the grid unit after projective transformationFor/>
4. The method for transforming planar coordinates of power grid homologous maintenance tool and GIS graphic migration according to claim 3, wherein in step 4, the specific method for projectively transforming GIS graphic coordinates into CGSC2000 planar coordinate system by aligning GIS grid and CGSC2000 planar coordinate system grid is as follows:
Calculating the coordinates (x g,yg) of GIS coordinates (B, L) of equipment in the power distribution station relative to a reference point by taking the lower left corner of the station frame as the reference point,
Translating the GIS coordinates (B, L) of the equipment in the power distribution station relative to the coordinates (x g,yg) of the reference point according to the size of the grid unit after projection transformation to obtain CGSC2000 plane coordinates (x, y) after projection transformation,
5. The method for converting plane coordinates of power grid homology maintenance tool and GIS graphic migration according to claim 4, wherein in step 6, revising a coordinate double-precision value format, and adding grid parameters to the CGSC2000 plane coordinates after the projection transformation in step 4 comprises the following specific steps:
For the CGSC2000 plane coordinate (x, y) obtained in the step 4, intercepting the integer part of the plane coordinate and the value of three bits after the decimal point as an original coordinate value, and taking the data of 3 bits after the decimal point as a millimeter-level deviation value;
when the CGSC2000 plane coordinate obtained in the step 4 has no deviation with the plane coordinate obtained by directly converting the GIS graph coordinate, the original coordinate value obtained by the CGSC2000 plane coordinate obtained in the step 4 is taken as the final power distribution network homologous maintenance graph coordinate; when the CGSC2000 plane coordinate obtained in the step 4 is deviated from the plane coordinate obtained by directly converting the GIS graphic coordinate, a millimeter-level deviation value is added to the original coordinate value, and the CGSC2000 plane coordinate after the millimeter-level deviation value is added is the final power distribution network homologous maintenance graphic coordinate.
6. The planar coordinate conversion method for graph migration between a power grid homologous maintenance tool and a GIS according to claim 1, wherein the method further comprises a coordinate conversion process of writing back the graph data newly generated by the power grid homologous maintenance tool to the GIS, and the coordinate conversion process of writing back the graph data newly generated by the power grid homologous maintenance tool to the GIS specifically comprises the following steps:
step A, newly generating plane coordinates of a graph for a power distribution network maintenance tool, and taking a integer value according to millimeter precision;
Step B, calculating a millimeter-level deviation value, correcting the plane coordinates of the newly generated graph of the power distribution network maintenance tool according to the millimeter-level deviation value, accumulating the millimeter-level deviation value for the plane coordinates of the newly generated graph of the power distribution network maintenance tool when the absolute value of the millimeter-level deviation value is larger than 1, and not processing the plane coordinates of the newly generated graph of the power distribution network maintenance tool when the absolute value of the millimeter-level deviation value is not larger than 1;
And C, carrying out coordinate conversion on the plane coordinates of the newly generated graph of the corrected power distribution network maintenance tool, carrying out Gaussian-Kelvin back calculation on the plane coordinates of the newly generated graph of the power distribution network maintenance tool to obtain CGSC2000 ellipsoidal geodetic coordinates, converting the CGSC2000 ellipsoidal geodetic coordinates into CGSC2000 ellipsoidal space rectangular coordinates, converting the CGSC2000 ellipsoidal space rectangular coordinates into WGS-84 ellipsoidal space rectangular coordinates by adopting a Boolean seven-parameter formula, and converting the WGS-84 ellipsoidal space rectangular coordinates into WGS-84 geodetic coordinates which are written back to a GIS system.
7. The method for converting plane coordinates of a graph migration between a power grid homologous maintenance tool and a GIS according to claim 6, wherein in step C, the specific method for performing gaussian-k-g inverse calculation on the plane coordinates of the graph newly generated by the power grid maintenance tool is as follows:
step C-1, iteratively solving the bottom point latitude B f, and calculating an initial value of the bottom point latitude B f: Up to/>
And C-2, back calculating GIS coordinates (B, L) corresponding to the plane coordinates (x, y) of the newly generated graph of the power distribution network maintenance tool according to the following formula:
L＝l+L0，
Wherein, t f＝tanBf is a total number of the components, aD、/>The long half shaft, the first eccentricity and the second eccentricity of the WGS-84 ellipsoid.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210707929.2A CN114936261B (en) | 2022-06-21 | 2022-06-21 | Plane coordinate conversion method for power grid homologous maintenance tool and GIS graph migration |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210707929.2A CN114936261B (en) | 2022-06-21 | 2022-06-21 | Plane coordinate conversion method for power grid homologous maintenance tool and GIS graph migration |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114936261A CN114936261A (en) | 2022-08-23 |
| CN114936261B true CN114936261B (en) | 2024-04-19 |
Family
ID=82868775
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210707929.2A Active CN114936261B (en) | 2022-06-21 | 2022-06-21 | Plane coordinate conversion method for power grid homologous maintenance tool and GIS graph migration |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114936261B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114001682B (en) * | 2021-11-19 | 2023-08-01 | 天津博迈科海洋工程有限公司 | Flatness judging method for mounting surface of heavy door frame of electric room module |
| CN117251518B (en) * | 2023-09-28 | 2024-06-21 | 广东省城乡规划设计研究院科技集团股份有限公司 | Urban safety planning design method for synchronous operation of GIS platform and CAD platform |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20060100157A (en) * | 2005-03-16 | 2006-09-20 | 윤홍식 | Distorting modeling method for transforming the presize position of partial/positional information |
| CN103441988A (en) * | 2013-08-02 | 2013-12-11 | 广东电网公司电力科学研究院 | Data migration method crossing GIS platforms |
| CN112052908A (en) * | 2020-09-17 | 2020-12-08 | 公安部交通管理科学研究所 | Traffic accident location clustering method and system |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9883337B2 (en) * | 2015-04-24 | 2018-01-30 | Mijix, Inc. | Location based services for RFID and sensor networks |
-
2022
- 2022-06-21 CN CN202210707929.2A patent/CN114936261B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20060100157A (en) * | 2005-03-16 | 2006-09-20 | 윤홍식 | Distorting modeling method for transforming the presize position of partial/positional information |
| CN103441988A (en) * | 2013-08-02 | 2013-12-11 | 广东电网公司电力科学研究院 | Data migration method crossing GIS platforms |
| CN112052908A (en) * | 2020-09-17 | 2020-12-08 | 公安部交通管理科学研究所 | Traffic accident location clustering method and system |
Non-Patent Citations (1)
| Title |
|---|
| 昆仑山口西Ms8.1地震的地壳变形特征;乔学军, 王琪, 杜瑞林, 陈宗时, 游新兆, 谭凯;大地测量与地球动力学;20021130(第04期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114936261A (en) | 2022-08-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114936261B (en) | Plane coordinate conversion method for power grid homologous maintenance tool and GIS graph migration | |
| CN103389080B (en) | A kind of based on geographic information application Urban Independent Coordinate System parameter acquiring method | |
| CN105160197A (en) | Comprehensive geographical spatial data coordinate conversion method and system | |
| CN105158760B (en) | Method for inverting underground fluid volume change and three dimension surface deformation using InSAR | |
| Markič et al. | Enabling geodetic coordinate reference systems in building information modeling for infrastructure | |
| CN111696156A (en) | Control point-free remote sensing image coordinate conversion method | |
| CN111307125A (en) | Inclined-axis cylindrical projection method based on GNSS and ground ranging combined adjustment | |
| CN110989886B (en) | Three-dimensional space grid selection method and device based on space map | |
| Markič et al. | Requirements for geo-locating transnational infrastructure BIM models | |
| CN115270522B (en) | Method and device for simulating and tracking target equipment based on WGS84 coordinates | |
| Habib et al. | A low-cost spatial tool for transforming feature positions of CAD-based topographic mapping | |
| CN102324096B (en) | Method and device for splicing characteristic point correction partitions | |
| Jaud et al. | Georeferencing within IFC: A Novel Approach for Infrastructure Objects | |
| Yakubu et al. | Ramification of datum and ellipsoidal parameters on post processed differential global positioning system (DGPS) data–A case study | |
| CN114091165B (en) | Method for fusing grid data of BIM application engineering coordinate system in GIS scene | |
| Habib et al. | Developing low-cost automated tool for integrating maps with GNSS satellite positioning data | |
| Kadir et al. | Geocentric datum GDM2000 for Malaysia: Implementation and implications | |
| Klu | Determination of a geoid model for Ghana using the Stokes-Helmert method | |
| CN114964169B (en) | Remote sensing image adjustment method for image space object space cooperative correction | |
| CN110542410A (en) | Geodetic survey result conversion method based on grid model | |
| Portnov | Research to justify the perspective use of individual coordinate systems for contiguous objects modelling with the case study of Moscow saints Petersburg railway | |
| CN113327012B (en) | Urban public transport index calculation method based on RGB color space and Monte Carlo method | |
| CN114490902B (en) | Multi-dimensional space self-adaptive subdivision and coding method and system for two-dimensional geographic entity | |
| Ngo | Research on terrestrial laser scanning technology in terrain and non-terrain fields in Vietnam | |
| CN112434393B (en) | Marketing water meter positioning method using hundred-degree map API and ArcGIS |
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 |