CN101320488B - Global ocean triangular net construction method - Google Patents
Global ocean triangular net construction method Download PDFInfo
- Publication number
- CN101320488B CN101320488B CN2008101382169A CN200810138216A CN101320488B CN 101320488 B CN101320488 B CN 101320488B CN 2008101382169 A CN2008101382169 A CN 2008101382169A CN 200810138216 A CN200810138216 A CN 200810138216A CN 101320488 B CN101320488 B CN 101320488B
- Authority
- CN
- China
- Prior art keywords
- triangular
- mesh unit
- triangular mesh
- numbering
- apex
- 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.)
- Expired - Fee Related
Links
Images
Abstract
The present invention relates to a global ocean triangular gridding constructing method in the marine science and ocean engineering numerical simulation field. Firstly, the global ocean is zoned continuously, and each subarea is processed for geographical projection transformation respectively; then the prior plane triangular gridding manufacturing method is used for constructing triangular gridding for each subarea respectively to obtain triangular griddings which are dissevered mutually to form a unit; finally, triangle vertexes and triangles in each triangular gridding unit are numbered uniformly to finish the combination of the triangular gridding unit, and a global ocean triangular gridding is constructed. The present invention overcomes the difficulty for constructing nonplanar triangular gridding, achieves the construction of the global ocean triangular gridding, avoids open boundary condition successfully, realizes the aim of manufacturing high-resolution triangular gridding at concerned sea area and sets a solid foundation for exact ocean numerical simulation and ocean forecast.
Description
Technical field
The present invention relates to the construction method of the global ocean triangular net in Marine Sciences and the oceanographic engineering numerical simulation field.
Background technology
China has more than 300 ten thousand square kilometres blue territory, and area vast, aboundresources in ocean is badly in need of exploitation and management at present.Along with Global climate change, the especially variation of marine climate, Oceanic disasters such as coastal generation storm tide of China are frequent further in recent years, and Oceanic disasters bring about great losses all for every year China people's life and property.In addition, along with the fast development of China's economy and trade, marine environmental pollution incidents such as ship spill accident and red tide disaster frequently take place, and have brought very big pressure for China's fish production and the marine eco-environment.Therefore, research Global climate change and coastal ocean disaster are forecast climate change and coastal ocean disaster, forecast significant to the distribution and the motion path of marine pollutant.Realize the accurate marine forecasting in marine site on a large scale, using a computer, to carry out the ocean numerical simulation be effective method.Yet present high precision, high-resolution ocean numerical model is regional ocean model, and this model depends critically upon out choosing of boundary condition, again can't objective, accurate acquisition but open boundary condition, can only artificially adjust, applicability is very poor, subjective, the value of forecasting is bad.On the other hand, rectangular node is all adopted in global ocean simulation at present, for reducing calculated amount, satisfies the requirement of computing power, and grid resolution is thicker, can't carry out careful simulation to being concerned about the marine site.Therefore, construct and have " being concerned about the marine site " local refinement, and " other zone " grid is sparse, need not out the global ocean triangular net of boundary condition, approach necessitates.For accurate marine forecasting, this is a crucial difficult problem of being badly in need of solution, but existing triangular mesh constructing technology all is only limited to the plane and makes up, can't finish have connectedness, the structure of the global ocean triangular net of sphere feature and become the bottleneck of accurate marine forecasting.
Summary of the invention
The objective of the invention is to remedy the deficiency of existing triangular mesh method for making, technology---the construction method of global ocean triangular net that makes up global ocean triangular net is provided, overcome and have difficulty connective and that the sphere triangular mesh makes up, open under the prerequisite on border the purpose that adopts the high resolving power triangular mesh to carry out accurate numerical simulation and marine forecasting to " care marine site " so that be implemented in not have.
The present invention is the improvement to existing plane triangle grid method for making, be sphere at global surface, the characteristics that the ocean is interconnected, carry out the Geographical projections conversion respectively with the continuous subregion of global ocean and to each subregion, utilize existing plane triangle grid method for making that above-mentioned each subregion is made grid respectively then, the triangular form grid component units that is isolated mutually, by merging " boundary line altogether ", at last with triangular apex in each triangular mesh unit and the unified numbering of triangle, after finishing the combination of triangular mesh unit, promptly be built into global ocean triangular net.
Concrete grammar of the present invention or step are as follows:
(1) with the continuous subregion of global ocean and carry out the corresponding geographical projective transformation of each subregion: the continuous subregion of global ocean is at least two non-overlapping marine sites, should there be the identical boundary line of geographic coordinate junction, adjacent marine site, this line is hereinafter referred to as " boundary line altogether ", each marine site basis is carried out projective transformation in the geographic position separately, is about to terrestrial coordinate (latitude and longitude coordinates) and is converted into local projection coordinate.
(2) utilize existing plane triangle grid method for making in above-mentioned each continuous subregion, to make triangular mesh respectively, form a plurality of " triangular mesh unit " that comprise triangular apex numbering and triangle numbering, and must satisfy when making the triangular mesh unit: on common boundary line, adjacent marine site, the identical triangular apex in geographic position is set.
(3) each triangular mesh unit intermediate cam shape apex coordinate of making in the step (2) is carried out and the middle corresponding projection inverse transformation of step (1), be about to local projection coordinate and be converted into unified terrestrial coordinate (latitude and longitude coordinates), so that carry out the combination of triangular mesh unit.
(4) utilize again in the step (2), adjacent marine site is total to the same geographic location triangular apex that is provided with on the boundary line, by being total to bordering compounding marine site triangular mesh unit successively, boundary line, and carry out triangular apex and leg-of-mutton renumbeing, so that each triangular mesh unit unification in same triangular mesh numbering system, is realized the structure of global ocean triangular net.
The present invention has overcome the difficulty that the on-plane surface triangular mesh makes up, finished the structure of global ocean triangular net, successfully avoided introducing and opened boundary condition, realized making the purpose of high resolving power triangular mesh, laid solid foundation for carrying out accurate ocean numerical simulation and marine forecasting in " care marine site ".
Description of drawings
Fig. 1, three continuous subregion example schematic of global ocean.Wherein a is triangular mesh unit, marine site, the arctic, and b is triangular mesh unit, marine site, the South Pole, and c is triangular mesh unit, middle latitude marine site; Digital 1-5 among Fig. 1, indicated boundary line is different common boundary line, the common boundary line of same numbers indication is expression identical common boundary line, geographic position.
Fig. 2, the global ocean triangular net that is built into are from the observation figure (example) of three different visual angles.
Embodiment
The present invention has provided the final detailed step that generates from the continuous subregion of global ocean to global ocean triangular net, describe and understanding for convenient simultaneously, as an embodiment, global ocean is divided into three continuous subregions carries out global triangular mesh and make up as follows:
(1) with the continuous subregion of global ocean and carry out the corresponding geographical projective transformation of each subregion-according to research or needs such as business, the continuous subregion of planning global ocean, for example by continuous subregion shown in Figure 1, global ocean is divided into three continuous marine sites, be respectively marine site, the arctic (shown in Fig. 1-a), marine site, the South Pole (shown in Fig. 1-b) and middle latitude marine site (shown in Fig. 1-c), each marine site can be selected suitable geographic coordinate projection pattern as required, adopt the polar region tomographic projection as marine site, the arctic and marine site, the South Pole, conic projection etc. is adopted in the middle latitude marine site.Zero lap marine site, adjacent marine site, junction only have same boundary line (boundary line altogether), represent the common boundary line of junction, adjacent marine site as the identical line of label among Fig. 1.Utilize above-mentioned continuous subregion and adopt the projecting method that is fit to particular geographic location, solved that global ocean has connective and sphere stereoscopic feature and the difficulty that can't in one plane finish the global ocean triangular net laying at different subregions.
(2) utilize existing plane triangle grid method for making, in each subregion, make triangular mesh respectively, be formed for making up the triangular mesh unit of global ocean triangular net, comprise triangular apex numbering and triangle numbering etc. in the triangular mesh unit.Fig. 1-a for example, Fig. 1-b shown in Fig. 1-c, is respectively the making result of the triangular mesh unit of marine site, the arctic, three subregions in marine site, the South Pole and middle latitude marine site.The requirement of making the triangular mesh unit is: the needs that 1, satisfy local refinement grids such as research and business; 2, satisfying the identical condition in geographic position that adjacent marine site is total to triangular apex on the boundary line, more than is the necessary condition of combined triangular shape grid cell in the following step.
(3) the triangular mesh unit intermediate cam shape apex coordinate in above-mentioned (2) is carried out and the middle corresponding projection inverse transformation of step (1), be about to local projection coordinate and be converted into terrestrial coordinate (latitude and longitude coordinates).
(4) each the triangular mesh unit that forms in the combination (3) at last is about to each triangular mesh unit and makes up successively, forms complete global triangular mesh gradually.Detailed process is earlier the triangular mesh unit to be merged boundary line altogether, then intermediate cam shape summit, triangular mesh unit and triangle are renumberd, finish the unification of global ocean triangular net triangle summit numbering and triangle numbering, finally construct global ocean triangular net.
Merge the triangular mesh unit altogether the method in boundary line be: with some triangular mesh unit (being called unit A) altogether the triangular apex numbering correspondence on the boundary line replace with another adjacent triangular mesh unit (being called unit B) has same geographic location altogether on the boundary line with it triangular apex numbering.
A kind of method that intermediate cam shape summit, triangular mesh unit renumbers is: according to the order of A intermediate cam shape summit, unit numbering, add 1 from unit B intermediate cam shape summit numbering sum, with among the unit A except that unit A and unit B altogether on the boundary line triangular apex the triangular apex number and increase by 1 successively, like this A intermediate cam shape summit, unit is renumberd, after triangular apex renumberd and finishes, the triangular apex sum that the maximum of the triangular apex among the unit A is numbered among the unit A added that the triangular apex sum in the unit B deducts triangular apex sum on both common boundary lines.
A kind of method of triangle numbering is: according to the order of unit A intermediate cam shape numbering, triangle numbering sum from unit B adds 1 beginning, unit A intermediate cam shape numbering is increased by 1 successively, like this unit A intermediate cam shape is renumberd, after triangle renumberd and finishes, the leg-of-mutton maximum numbering among the unit A should be triangle sum among the unit A and adds triangle sum in the unit B.
So just, finished the combination of triangular mesh unit A and triangular mesh unit B, form a new triangular mesh unit, according to as above process other triangular mesh unit are integrated with above-mentioned new triangular mesh unit successively below, just be built into global triangular mesh at last.
After for example the triangular mesh unit shown in Fig. 1-b and Fig. 1-c being merged, again Fig. 1-a is integrated with wherein, just can finish the structure of global ocean triangular net.If total m triangular apex in Fig. 1-b triangular mesh unit, p triangle, total n triangular apex in Fig. 1-c network of triangle lattice file, q triangle, all total k triangular apex on both common boundary lines 4.
No. 4 altogether merging in boundary line: with in Fig. 2-b triangular mesh unit No. 4 altogether the numbering of all triangular apex on the boundary line replace with in Fig. 1-c triangular mesh unit the corresponding identical triangular apex numbering in geographic position on the boundary line altogether No. 4.
Triangular apex and leg-of-mutton renumbeing: again the triangular apex on the non-No. 4 common boundary lines in Fig. 1-b triangular mesh unit is numbered according to existing sequencing, from the n+1 open numbering, number n+m-k number, triangle is according to existing sequencing, from the p+1 open numbering, number p+q.So just, finish the merging of triangular mesh unit shown in Fig. 1-b and Fig. 1-c, formed new triangular mesh unit.
Again with 1,2,3, No. 5 boundary line merging altogether, triangular apex and triangle renumber according to above same method, and Fig. 1-a integrates with above triangular mesh unit, has just finished the structure of global ocean triangular net, as shown in Figure 2.
Claims (4)
1. the construction method of global ocean triangular net, concrete grammar is as follows:
(1) with the continuous subregion of global ocean and carry out the corresponding geographical projective transformation of each subregion: the continuous subregion of global ocean is at least two non-overlapping marine sites, should there be the identical boundary line of geographic coordinate junction, adjacent marine site, this line is hereinafter referred to as " boundary line altogether ", each marine site basis is carried out projective transformation in the geographic position separately, and being about to greatly, latitude and longitude coordinates is converted into local projection coordinate;
(2) utilize existing plane triangle grid method for making in above-mentioned each continuous subregion, to make triangular mesh respectively, form a plurality of " triangular mesh unit " that comprise triangular apex numbering and triangle numbering, the condition that must satisfy when making the triangular mesh unit is: on the boundary line the identical triangular apex in geographic position is set altogether in adjacent marine site;
(3) each triangular mesh unit intermediate cam shape apex coordinate of making in the step (2) is carried out and the middle corresponding projection inverse transformation of step (1), be about to local projection coordinate and be converted into the earth latitude and longitude coordinates, so that carry out the combination of triangular mesh unit;
(4) at last each the triangular mesh unit that forms in (3) is made up successively, and form complete global triangular mesh gradually: earlier the triangular mesh unit is merged boundary line altogether, then intermediate cam shape summit, triangular mesh unit and triangle are renumberd, and finish the structure of global ocean triangular net triangle gridding.
2. the construction method of global ocean triangular net as claimed in claim 1, it is characterized in that above-mentioned merging triangular mesh unit altogether the method in boundary line be: the triangular apex numbering on the boundary line is corresponding altogether replaces with another adjacent triangular mesh unit B has same geographic location altogether on the boundary line with it triangular apex numbering with a triangular mesh unit A.
3. the construction method of global ocean triangular net as claimed in claim 1, it is characterized in that above-mentioned triangular apex method for numbering serial is as follows: according to the order of A intermediate cam shape summit, triangular mesh unit numbering, add 1 from triangular mesh unit B intermediate cam shape summit numbering sum, with among the triangular mesh unit A except that triangular mesh unit A and triangular mesh unit B altogether on the boundary line triangular apex the triangular apex number and increase by 1 successively, diabolo grid cell A intermediate cam shape summit renumbers like this, and triangular apex renumbers the triangular apex sum that the maximum of the triangular apex among the Vee formation shape grid cell A that finishes is numbered among the triangular mesh unit A and adds that the triangular apex sum in the triangular mesh unit B deducts triangular apex sum on both common boundary lines.
4. the construction method of global ocean triangular net as claimed in claim 1, it is characterized in that above-mentioned triangle method for numbering serial is as follows: according to the order of triangular mesh unit A intermediate cam shape numbering, triangle numbering sum from the triangular mesh unit B adds 1 beginning, triangular mesh unit A intermediate cam shape numbering is increased by 1 successively, diabolo grid cell A intermediate cam shape renumbers like this, after triangle renumberd and finishes, the leg-of-mutton maximum numbering among the triangular mesh unit A should be triangle sum among the triangular mesh unit A and adds triangle sum in the triangular mesh unit B.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2008101382169A CN101320488B (en) | 2008-06-27 | 2008-06-27 | Global ocean triangular net construction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2008101382169A CN101320488B (en) | 2008-06-27 | 2008-06-27 | Global ocean triangular net construction method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101320488A CN101320488A (en) | 2008-12-10 |
| CN101320488B true CN101320488B (en) | 2011-07-13 |
Family
ID=40180519
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2008101382169A Expired - Fee Related CN101320488B (en) | 2008-06-27 | 2008-06-27 | Global ocean triangular net construction method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101320488B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101826317B (en) * | 2010-01-26 | 2012-01-18 | 北京水晶石数字科技股份有限公司 | Sea-sky mergence real-time simulation method |
| CN105761309B (en) * | 2014-12-16 | 2019-02-05 | 高德软件有限公司 | A kind of method and device handling digital terrain model seam |
| CN106156472A (en) * | 2015-04-20 | 2016-11-23 | 中国科学院遥感与数字地球研究所 | A kind of remote sensing satellite covering power over the ground analyzes method |
| CN107886468A (en) * | 2016-09-29 | 2018-04-06 | 阿里巴巴集团控股有限公司 | Mapping method, reconstruction, processing method and the corresponding intrument and equipment of panoramic video |
| US10922783B2 (en) * | 2018-03-02 | 2021-02-16 | Mediatek Inc. | Cube-based projection method that applies different mapping functions to different square projection faces, different axes, and/or different locations of axis |
| CN109754449B (en) * | 2018-11-22 | 2020-04-03 | 清华大学 | Triangularization determination method for two-dimensional grid graph |
| CN111753451A (en) * | 2020-06-23 | 2020-10-09 | 中国水利水电科学研究院 | Unstructured grid splitting and merging method suitable for water conservancy relevant numerical simulation |
-
2008
- 2008-06-27 CN CN2008101382169A patent/CN101320488B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN101320488A (en) | 2008-12-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101320488B (en) | Global ocean triangular net construction method | |
| CN103236086B (en) | One takes the contextual multiple dimensioned DEM modeling method of the earth's surface hydrology into account | |
| Xiao et al. | Saltwater intrusion into groundwater systems in the Mekong Delta and links to global change | |
| CN103218517A (en) | GIS (Geographic Information System)-based region-meshed spatial population density computing method | |
| CN116069882B (en) | Airspace grid diagram generating method | |
| CN104573234B (en) | A kind of basin storage capacity density space method based on ArcGIS | |
| CN108053458A (en) | A kind of watershed generalization figure based on GIS-Geographic Information System makes and methods of exhibiting | |
| CN105354882A (en) | Method for constructing big data architecture based three-dimensional panoramic display platform for large-spatial-range electricity transmission | |
| CN111221925A (en) | Method and device for monitoring and networking wind-waterlogging disaster of power distribution network | |
| CN114881466A (en) | Multi-source data-based population space partition fitting method | |
| Duan et al. | The spatial form of traditional villages in Fuzhou area of Jiangxi Province determined via GIS methods | |
| Liu et al. | Construction and Analysis of Ecological Security Pattern of Qingdao Based on MSPA and MCR Models. | |
| Zhang et al. | Hydrologic information extraction based on arc hydro tool and DEM | |
| Xu et al. | Urban flooding resilience evaluation with coupled rainfall and flooding models: a small area in Kunming City, China as an example | |
| Lv et al. | Spatial patterns of China's major cities and their evolution mechanisms during the past decades of reform and opening up | |
| CN107122779A (en) | A kind of space basis geographical data acquisition system | |
| CN108009562B (en) | Method for identifying characteristic spatial variability of hydrology and water resource | |
| Chu et al. | Remote sensing-based life cycle analysis of land reclamation processes: case study on Tianjin Binhai New area | |
| Huang et al. | Ecological red line zoning of the Tibet autonomous region based on ecosystem services and ecological sensitivity | |
| CN102810157B (en) | TIN (Triangulated Irregular Network) data-based river network extraction method | |
| Zhang et al. | Spatial patterns and controlling factors of the evolution process of karst depressions in Guizhou province, China | |
| Raj et al. | Analysis of landuse landcover changes of Kazhakuttam block based on GIS | |
| Garau et al. | Investigating the socio-spatial logic of historic urban areas through space syntax. A comparative analysis of the Roman towns Cosa, Nora, Timgad and Thuburbo Majus | |
| Sulaiman et al. | Integration of geographic information system (GIS) and hydraulic modelling to simulate floodplain inundation level for Bandar Segamat | |
| Li | Study of urbanization trend prediction in Nanchang based on CA-Markov model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110713 Termination date: 20200627 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |