CN102411872A - Method for generating line element buffer zone on earth elliposoid based on GIS (geographic information system) component - Google Patents
Method for generating line element buffer zone on earth elliposoid based on GIS (geographic information system) component Download PDFInfo
- Publication number
- CN102411872A CN102411872A CN201110232771XA CN201110232771A CN102411872A CN 102411872 A CN102411872 A CN 102411872A CN 201110232771X A CN201110232771X A CN 201110232771XA CN 201110232771 A CN201110232771 A CN 201110232771A CN 102411872 A CN102411872 A CN 102411872A
- Authority
- CN
- China
- Prior art keywords
- earth
- buffer zone
- line feature
- line
- ellipsoid face
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Processing Or Creating Images (AREA)
Abstract
The invention discloses a method for generating a line element buffer zone on an earth elliposoid based on a GIS (geographic information system) component. The method comprises the following steps of: a, determining vertexes of line elements on an encrypted earth elliposoid; b, calculating out geodetic circles; c, calculating out geodetic rectangles; and d, adopting a CIS component to perform union set generation on the geodetic circles and rectangles established in the step b and step c, so as to obtain a polygon, namely a line element buffer zone on the earth elliposoid. By the method in the invention, the process for generating a line element buffer zone on the earth elliposoid is greatly simplified, and a convenient and efficient method for generating a line element buffer zone on the earth elliposoid on any map projection plane can be provided for users having GIS secondary development capability, an effective technical means for generating high-precision unilateral and bilateral ocean boundaries can be provided for ocean demarcation, and the method has an important significance on protecting national ocean rights and interests according to laws.
Description
Technical field
The present invention relates to GIS-Geographic Information System (GIS) and description field, ocean, relate in particular to a kind of based on line feature buffer zone generation method on the earth ellipsoid face of GIS assembly.
Background technology
At present, the method for on the map projection plane, setting up point, line, surface and compound key element (hereinafter to be referred as each key element) buffer zone is very ripe, and becomes the most basic a kind of spacial analytical method in various GIS-Geographic Information System (GIS) software already; And be that the user provides the application corresponding service with the form of GIS desktop software or GIS assembly (control); But, thereby can't satisfy the request for utilization in field such as ocean description because each key element buffer zone that this method is set up all can receive the influence of distortions of a map projection inevitably on precision; To this problem; Relevant scholar utilizes ellipsoid Geodetic Calculation method, has eliminated the influence of distortions of a map projection, and (Peng recognizes bright to have set up each buffer zone constructing technology based on earth's spheroid; Wang Jiayao. based on the buffer zone constructing technology research of earth's spheroid; The mapping journal, vol31,2002 (3): 270-273) (document 1); Peng recognizes bright, Chen Yi, Liu Guohui; Deng. MapInfo system line (face) destination buffer makes up the problem of model existence and improves one's methods, Wuhan University's journal information science version, and vol 32; 2007 (8): 719-722) (document 2), through being not difficult, the analysis of above-mentioned document finds that their summarys have provided the basic ideas of on the earth ellipsoid face, setting up each key element buffer zone; And about how utilizing the GIS assembly to set up the buffer zone of each key element on the earth ellipsoid face; Then also belong to technological gap, owing to put the special case that key element can be regarded line feature or face key element as, and face key element and compound key element can be formed by the point, line, surface factor combination; Therefore, only need research to get final product based on line feature buffer zone generation method on the earth ellipsoid face of GIS assembly.
Summary of the invention
The purpose of this invention is to provide a kind of based on line feature buffer zone generation method on the earth ellipsoid face of GIS assembly.
The present invention for realizing the technical scheme that above-mentioned purpose adopted is: a kind of based on line feature buffer zone generation method on the earth ellipsoid face of GIS assembly, may further comprise the steps:
A, determine each summit of line feature on the earth ellipsoid face after encrypting through the geodesic line interpolation model;
B, to calculate with line feature two-end-point on the earth ellipsoid face and each summit after encrypting through ellipsoid the earth computing formula be the center of circle, buffering apart from being the geodetic circle of radius;
C, to calculate with the line segment between each adjacent vertex after encrypting of line feature on the earth ellipsoid face through ellipsoid the earth computing formula be public base, its both sides with buffering apart from being high the earth rectangle;
D, with all geodetic circles set up among step b and the step c and the earth rectangle, utilize the CIS assembly to carry out union and generate to handle, the polygon that obtains is line feature buffer zone on the earth ellipsoid face.
On the earth ellipsoid face of encrypting without the geodesic interpolation model among the step a each apex coordinate of line feature be (b1, l1), (b2; L2) ..., (bm; Lm), then on the earth ellipsoid face after the geodesic interpolation model is encrypted each apex coordinate of line feature be (B1, L1); (B2, L2) ... (Bn; Ln), wherein, n >=m.
At first line feature two-end-point on the earth ellipsoid face and each summit are regarded as independently putting key element among the said step b, set up its geodetic circle on the earth ellipsoid face then respectively.
Among the said step c on the earth ellipsoid face line segment between line feature each adjacent vertex after encrypting be the segment geodesic line; By the positive solution of the earth problem respectively from each segment geodesic line two-end-point; The left side that to ask the position angle to become 90 ° and/or 270 ° of distances with it be the buffer zone radius and the latitude and longitude coordinates of parallel lines point position, right side, and pair warp and weft degree coordinate matches and is connected to form each required the earth rectangle.
The present invention is with the scale of measurement of geodesic line yardstick as the buffering distance; That is to say; Line feature is to the minor increment on its buffer zone border---and buffering is apart from calculating by the geodesic line length amount of carrying out; Correspondingly, be the center of circle, given buffering apart from being that closed region on the constructed earth ellipsoid face of radius is called geodetic circle with the end points of line feature on the earth ellipsoid face or summit; Line segment to form line feature on the earth ellipsoid face is public base, is called the earth rectangle with given buffering distance for the closed region on the constructed earth ellipsoid face of its both sides height.
The present invention is a kind of based on line feature buffer zone generation method on the earth ellipsoid face of GIS assembly; The various technical matterss of utilizing the GIS assembly to generate line feature buffer zone on the earth ellipsoid face have been solved; Promptly through ask on the earth ellipsoid face with line feature two-end-point and each summit be the center of circle, with buffering apart from being the geodetic circle of radius; And with the line segment between each adjacent vertex of line feature be public base, with buffering apart from method for the union of the high the earth rectangle in its both sides; Obtain the buffer zone of line feature on the earth ellipsoid face; The polygon union systematic function that the present invention directly utilizes general GIS assembly Arc Objects and SuperMap Objects to be provided; Simplified the process that generates line feature buffer zone on the earth ellipsoid face greatly; For the user with GIS secondary development ability provides a kind of ten minutes convenient, effective in the method that generates line feature buffer zone on the earth ellipsoid face on arbitrary map projection plane, utilize the present invention to demarcate a kind of effective technology means that generate the monolateral and bilateral ocean of high precision boundary are provided for the ocean, have great importance to safeguard the national marine rights and interests in accordance with the law.
Description of drawings
Fig. 1 is that the present invention is at the line feature synoptic diagram after the geodesic line interpolation model is encrypted on the cylindrical equidistant projection plane.
The line feature that Fig. 2 is the present invention after encrypting on the cylindrical equidistant projection plane fix a point pairing geodetic circle and the earth rectangle synoptic diagram.
Fig. 3 is the buffer zone synoptic diagram that the present invention utilizes SuperMap Objects assembly line feature on the earth ellipsoid face that generates on the cylindrical equidistant projection plane.
Fig. 4 is that the present invention utilizes line feature buffer zone synoptic diagram on the earth ellipsoid face that Arc Objects assembly generates.
Fig. 5 is that the present invention utilizes line feature buffer zone synoptic diagram on the earth ellipsoid face that SuperMap Objects assembly generates.
Embodiment
Like Fig. 1-shown in Figure 5, the present invention is based on line feature buffer zone generation method on the earth ellipsoid face of GIS assembly, its concrete steps are: a, determine each summit of line feature on the earth ellipsoid face after encrypting through the geodesic line interpolation model, suppose that each apex coordinate without line feature on the earth ellipsoid face of geodesic line interpolation model encryption is (b1; L1), (b2, l2) ... (then each apex coordinate of line feature is (B1 on the earth ellipsoid face after the geodesic line interpolation model is encrypted for bm, lm) (unit is degree, down together); L1), (B2, L2) ... (Bn, Ln), wherein, n >=m; B, to calculate with line feature two-end-point on the earth ellipsoid face and each summit after encrypting through ellipsoid the earth computing formula be the center of circle, buffering apart from being the geodetic circle of radius; This is exactly in fact at first line feature two-end-point on the earth ellipsoid face and each summit to be regarded as independently putting key element, sets up its geodetic circle on the earth ellipsoid face then respectively; C, to calculate with the line segment between each adjacent vertex after encrypting of line feature on the earth ellipsoid face through ellipsoid the earth computing formula be public base; Its both sides are high the earth rectangle with the buffering distance; Because the line segment on the earth ellipsoid face between line feature each adjacent vertex after encrypting is the segment geodesic line; Therefore can be by the positive solution of the earth problem respectively from each segment geodesic line two-end-point; The left side that to ask the position angle to become 90 ° and/or 270 ° of distances with it be the buffer zone radius and the latitude and longitude coordinates of parallel lines point position, right side, and it is matched be connected to form each required the earth rectangle; D, with all geodetic circles set up among step b and the step c and the earth rectangle, utilize the CIS assembly to carry out union and generate to handle, the polygon that obtains is line feature buffer zone on the earth ellipsoid face.
At present; Many powerful GIS assemblies have all been released both at home and abroad; Like the Arc Objects of U.S. ESRI company and China hypergraph Corporation's Super Map Objects etc., these assemblies all provide polygon union systematic function with different method of calling, under the VC++6.0 environment; The Arc Objects assembly of U.S. ESRI company provides a series of geometric set operation operators through the ItopologicalOperator interface; Component class such as GeometryBag, Multipoint, Point, Polygon, Polyline have all encapsulated this interface, and wherein method Union returns with reference to geometric figure and more geometric union according to the reference geometric figure and comparison geometric figure of appointment; Its prototype is: HRESULT Union (IGeometry * other, IGeometry * * resultGeom);
Rreturn value: the HRESULT type, show the execution result state, if be 0 to represent successfully, otherwise the expression failure, parameter: other:IGeometry * type variable, expression is geometric figure relatively; ResultGeom:IGeometry * * type variable; The union that expression is returned; Under the VC++6.0 environment; China's hypergraph Corporation's Super Map Objects assembly then provides the space topological operator through the CsoSpatialOperator interface, and it is encapsulated in the component class such as CsoGeoLine, CsoGeoLineM, CsoGeoRegion, CsoGeoText, and wherein method Union is used to calculate the union of two geometric objects; And generating a new geometric object, its prototype is: CsoGeometry CsoSpatialOperator.Union (CsoGeometry objOtherGeometry); Rreturn value: the CsoGeometry type, represent new geometric object, if for NULL represent the failure; Otherwise represent successfully parameter: objOtherGeometry:CsoGeometry type variable, expression participate in asking another geometric object also; Below utilize SuperMap Objects assembly in the main process that generates line feature buffer zone on the earth ellipsoid face on the cylindrical equidistant projection plane through 3 cut lines; Wherein, buffering is apart from being 1.5 nautical miles, and solid dot is an original vertices among Fig. 1; Its coordinate that with the degree is unit representes by (latitude, longitude) form is respectively (57.394402,120.489263) from left to right; (58.263160,121.687429), (57.387827,122.485861), (57.903221,123.667414); Hollow dots is for inserting summit (inserting a point in every 5000m) in encrypting, the coordinate of inserting the summit in the encryption between 1,2 original vertices is (57.430789,120.537996), (57.467157,120.586829), (57.503505,120.635756); (57.539833,120.684783), (57.576143,120.733906), (57.612434,120.783128), (57.648705,120.832448); (57.684957,120.881866), (57.721188,120.931383), (57.757402,120.980998), (57.793595,121.0307213); (57.829767,121.080528), (57.865921,121.130443), (57.902054,121.180457), (57.938168,121.230573); (57.974260,121.280789), (58.010335,121.331106), (58.046388,121.381524), (58.082421,121.432044); (58.118433,121.482666), (58.154426,121.53389), (58.190397,121.584217), (58.226350,121.635146); 2, the coordinate on slotting summit is (58.222907,121.725110) in the encryption between 3 original vertices, (58.182642,121.762708), (58.142367,121.800220), (58.102080,121.837647); (58.061783,121.874990), (58.021472,121.912250), (57.981153,121.949425), (57.940822,121.986517); (57.900478,122.023526), (57.860125,122.060451), (57.819762,122.097295); (57.779386,122.134056), (57.738999,122.170734), (57.698603,122.207333); (57.658194,122.243849), (57.617777,122.280284), (57.577347,122.316638); (57.536907,122.352912), (57.496456,122.389105), (57.455994,122.425218); The coordinate of inserting the summit in (57.415523,122.461252), the encryption between 3,4 original vertices is (57.416443,122.549945), (57.445029,122.614131); (57.473581,122.678417), (57.502100,122.742802), (57.530586,122.807288); (57.559038,122.871874), (57.587459,122.936563), (57.615846,123.001351); (57.644198,123.066240), (57.672519,123.131231), (57.700805,123.196323); (57.729057,123.261517), (57.757275,123.326812), (57.785461,123.392208); (57.813611,123.457706), (57.841729,123.523308), (57.869812,123.589012); The buffer zone that utilizes Union method line feature on the earth ellipsoid face that generates on the Mercator projection plane that Arc Objects assembly provides that is as shown in Figure 4; Wherein, The coordinate of 5 original vertices of line feature is respectively (80.985717 ,-47.870009) from left to right, (84.544953 ,-29.694633), (81.415986 ,-14.548422); (85.004481,6.584150), (82.625294,30.818086); The buffering distance is 20000 meters; The buffer zone that utilizes Union method line feature on the earth ellipsoid face that generates on the Mercator projection plane that SuperMap Objects assembly provides that is as shown in Figure 5, wherein, the coordinate of 5 original vertices of line feature is respectively (84.377886,61.759633) from left to right; (81.285814,72.650481); (84.117675,85.777197), (80.675369,106.332769), (85.419628,126.671967); The buffering distance is 20000 meters.
The polygon union systematic function that the present invention directly utilizes general GIS assembly Arc Objects and SuperMap Objects to be provided; Simplified the process that generates line feature buffer zone on the earth ellipsoid face greatly; For the user with GIS secondary development ability provides a kind of ten minutes convenient, effective in the method that generates line feature buffer zone on the earth ellipsoid face on arbitrary map projection plane; Utilize the present invention to demarcate a kind of effective technology means that generate the monolateral and bilateral ocean of high precision boundary are provided, have great importance safeguard the national marine rights and interests in accordance with the law for the ocean.
Claims (4)
1. one kind based on line feature buffer zone generation method on the earth ellipsoid face of GIS assembly, it is characterized in that, may further comprise the steps:
A, determine each summit of line feature on the earth ellipsoid face after encrypting through the geodesic line interpolation model;
B, to calculate with line feature two-end-point on the earth ellipsoid face and each summit after encrypting through ellipsoid the earth computing formula be the center of circle, buffering apart from being the geodetic circle of radius;
C, to calculate with the line segment between each adjacent vertex after encrypting of line feature on the earth ellipsoid face through ellipsoid the earth computing formula be public base, its both sides with buffering apart from being high the earth rectangle;
D, with all geodetic circles set up among step b and the step c and the earth rectangle, utilize the CIS assembly to carry out union and generate to handle, the polygon that obtains is line feature buffer zone on the earth ellipsoid face.
2. according to claim 1 a kind of based on line feature buffering area generation method on the earth ellipsoid face of GIS assembly, it is characterized in that: on the earth ellipsoid face of encrypting without the geodesic interpolation model among the step a each apex coordinate of line feature for (b1, l1); (b2, l2) ... (bm, lm), then each apex coordinate of line feature is (B1 on the earth ellipsoid face after the geodesic interpolation model is encrypted; L1), (B2, L2);, (Bn, Ln); Wherein, n >=m.
3. according to claim 1 a kind of based on line feature buffer zone generation method on the earth ellipsoid face of GIS assembly; It is characterized in that: at first line feature two-end-point on the earth ellipsoid face and each summit are regarded as independently putting key element among the said step b, set up its geodetic circle on the earth ellipsoid face then respectively.
4. according to claim 1 a kind of based on line feature buffer zone generation method on the earth ellipsoid face of GIS assembly; It is characterized in that: among the said step c on the earth ellipsoid face line segment between line feature each adjacent vertex after encrypting be the segment geodesic line; By the positive solution of the earth problem respectively from each segment geodesic line two-end-point; The left side that to ask the position angle to become 90 ° and/or 270 ° of distances with it be the buffer zone radius and the latitude and longitude coordinates of parallel lines point position, right side, and pair warp and weft degree coordinate matches and is connected to form each required the earth rectangle.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110232771XA CN102411872A (en) | 2011-08-15 | 2011-08-15 | Method for generating line element buffer zone on earth elliposoid based on GIS (geographic information system) component |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110232771XA CN102411872A (en) | 2011-08-15 | 2011-08-15 | Method for generating line element buffer zone on earth elliposoid based on GIS (geographic information system) component |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN102411872A true CN102411872A (en) | 2012-04-11 |
Family
ID=45913932
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110232771XA Pending CN102411872A (en) | 2011-08-15 | 2011-08-15 | Method for generating line element buffer zone on earth elliposoid based on GIS (geographic information system) component |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102411872A (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104008276A (en) * | 2014-05-08 | 2014-08-27 | 中国科学院南海海洋研究所 | Sea power indexes and method for calculating island reef marine rights and interests quantitatively |
| CN107729292A (en) * | 2017-11-16 | 2018-02-23 | 彭认灿 | The geodesic interpolating method of azimuth sagitta limit difference constraint based on Mercator projection |
| CN108427842A (en) * | 2018-03-13 | 2018-08-21 | 厦门理工学院 | Power pipe webmaster line core checking method based on interpolation dot buffer zone |
| CN111460685A (en) * | 2020-04-21 | 2020-07-28 | 中国地质科学院地质力学研究所 | Buffer area establishing method, analysis method and application thereof |
| CN113360837A (en) * | 2021-07-28 | 2021-09-07 | 国家海洋信息中心 | Sea area middle line demarcation method based on earth ellipsoid |
| CN115223371A (en) * | 2022-09-20 | 2022-10-21 | 深圳市城市交通规划设计研究中心股份有限公司 | Big data analysis system of electric bicycle and working method thereof |
-
2011
- 2011-08-15 CN CN201110232771XA patent/CN102411872A/en active Pending
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104008276A (en) * | 2014-05-08 | 2014-08-27 | 中国科学院南海海洋研究所 | Sea power indexes and method for calculating island reef marine rights and interests quantitatively |
| CN104008276B (en) * | 2014-05-08 | 2017-01-18 | 中国科学院南海海洋研究所 | Sea power indexes for calculating island reef marine rights and interests quantitatively and method for determining sea power indexes |
| CN107729292A (en) * | 2017-11-16 | 2018-02-23 | 彭认灿 | The geodesic interpolating method of azimuth sagitta limit difference constraint based on Mercator projection |
| CN107729292B (en) * | 2017-11-16 | 2020-11-03 | 中国人民解放军海军大连舰艇学院 | Large ground wire interpolation method based on azimuth arch height tolerance constraint of mercator projection |
| CN108427842A (en) * | 2018-03-13 | 2018-08-21 | 厦门理工学院 | Power pipe webmaster line core checking method based on interpolation dot buffer zone |
| CN108427842B (en) * | 2018-03-13 | 2022-02-22 | 厦门理工学院 | Power pipe network pipeline checking method based on interpolation point buffer area |
| CN111460685A (en) * | 2020-04-21 | 2020-07-28 | 中国地质科学院地质力学研究所 | Buffer area establishing method, analysis method and application thereof |
| CN111460685B (en) * | 2020-04-21 | 2020-12-01 | 中国地质科学院地质力学研究所 | Buffer area establishing method, analysis method and application thereof |
| CN113360837A (en) * | 2021-07-28 | 2021-09-07 | 国家海洋信息中心 | Sea area middle line demarcation method based on earth ellipsoid |
| CN113360837B (en) * | 2021-07-28 | 2023-03-10 | 国家海洋信息中心 | Sea area middle line demarcation method based on earth ellipsoid |
| CN115223371A (en) * | 2022-09-20 | 2022-10-21 | 深圳市城市交通规划设计研究中心股份有限公司 | Big data analysis system of electric bicycle and working method thereof |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105160197A (en) | Comprehensive geographical spatial data coordinate conversion method and system | |
| CN102411872A (en) | Method for generating line element buffer zone on earth elliposoid based on GIS (geographic information system) component | |
| CN109670005B (en) | Position matching method of BIM (building information modeling) model and three-dimensional geographic scene | |
| CN103196426A (en) | Building surveying method utilizing total station and three-dimensional laser scanner | |
| KR101659780B1 (en) | Method for Representing Global Spherical Space of 2 Dimentional Drawing data in GIS | |
| KR100686287B1 (en) | Distorting Modeling method for Transforming the Presize Position of Partial/Positional information | |
| KR102097416B1 (en) | An augmented reality representation method for managing underground pipeline data with vertical drop and the recording medium thereof | |
| Gümüş et al. | Evaluation of NRTK GNSS positioning methods for displacement detection by a newly designed displacement monitoring system | |
| CN104794314A (en) | Three-dimensional simulation model of power transmission line and construction method of three-dimensional simulation model | |
| Tong et al. | Modeling and expression of vector data in the hexagonal discrete global grid system | |
| CN102223709A (en) | Method for acquiring specific reference Gaussian plane rectangular coordinate based on continuous operational reference system (CORS) in real time | |
| CN101441075A (en) | Map projection method using remote sense map tool | |
| CN110363844B (en) | Coal mine tunnel three-dimensional modeling method and system | |
| KR101756858B1 (en) | Secure Trajectory Query Processing with Surrogate Vector Model | |
| Jaud et al. | Georeferencing within IFC: A Novel Approach for Infrastructure Objects | |
| CN113160403B (en) | Modeling method of high-precision highway information model | |
| Kim et al. | Efficient encoding and decoding extended geocodes for massive point cloud data | |
| CN111061824A (en) | Range judgment method, device and equipment based on improved quadtree | |
| CN111426329A (en) | Road generation method and device, computer equipment and storage medium | |
| CN111353203A (en) | Seismic acquisition arrangement optimization method and device and storage medium | |
| CN105976314A (en) | Reference elliptical projection plane coordinate system converting method based on identical central project longitude | |
| KR20200145257A (en) | Method and apparatus for generating path profile | |
| Karimi et al. | Grid-based geoprocessing for integrated Global Navigation Satellite System simulation | |
| SUN et al. | Detection of changes of the Earth's volume and geometry by using GPS and VLBI data | |
| Lu et al. | Reference ellipsoid and the geodetic coordinate system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C12 | Rejection of a patent application after its publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20120411 |