CN116402974A - Geonet-like geolevel refinement model construction method for geographic boundary constraint - Google Patents
Geonet-like geolevel refinement model construction method for geographic boundary constraint Download PDFInfo
- Publication number
- CN116402974A CN116402974A CN202310010593.9A CN202310010593A CN116402974A CN 116402974 A CN116402974 A CN 116402974A CN 202310010593 A CN202310010593 A CN 202310010593A CN 116402974 A CN116402974 A CN 116402974A
- Authority
- CN
- China
- Prior art keywords
- model
- grid
- geographic boundary
- geographic
- point
- 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
- 238000010276 construction Methods 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 19
- 238000005259 measurement Methods 0.000 claims abstract description 14
- 230000001788 irregular Effects 0.000 claims abstract description 13
- 238000011084 recovery Methods 0.000 claims abstract description 5
- 230000005484 gravity Effects 0.000 claims description 20
- 230000002159 abnormal effect Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 claims description 3
- 238000010845 search algorithm Methods 0.000 claims description 3
- 230000005856 abnormality Effects 0.000 claims 1
- 238000013178 mathematical model Methods 0.000 abstract description 2
- 238000012876 topography Methods 0.000 abstract description 2
- 230000007704 transition Effects 0.000 abstract description 2
- 229920006068 Minlon® Polymers 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- 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/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
Abstract
The invention discloses a construction method of a geodesic level surface refinement model constrained by geographic boundaries, which specifically comprises the following steps: step one, boundary establishment; step two, measuring control point data; thirdly, building a model; step four, cutting constraint; fifthly, generating a model; the invention relates to the technical field of reservoir measurement. The geonet-like ground level refinement model construction method of geographic boundary constraint has geographic space applicability, avoids infinite space expansion caused by adopting a single mathematical model fitting method, causes uncertain model errors, takes regional index variation functions as weighting factors, considers influences of reservoir management ranges, topography factors and the like, realizes smooth transition of an elevation reference model of an irregular reservoir management range, simulates reality conditions more truly and accurately, and simultaneously realizes removal and recovery functions easily by adopting space addition and subtraction operations of continuous grid models.
Description
Technical Field
The invention relates to the technical field of reservoir measurement, in particular to a grid-type geoid refinement model construction method for geographic boundary constraint.
Background
In a complex irregular reservoir area, which is limited by factors such as geographical environment, traffic condition and basic data, a traditional elevation control measurement method cannot be adopted to construct an elevation reference model with high precision in a reservoir management and protection range, at present, the area elevation reference model is mainly constructed by realizing a refined regional geoid, the main method comprises a gravity method, a geometric method and a combination method, the normal height and the geoid of GNSS points are mainly utilized, the whole regional geoid is fitted by means of a single mathematical function in combination with gravity data, and particularly in mountain areas and hilly areas, a continuous smooth grid-type geoid cannot be formed, and the elevation reference plane of a complex reservoir cannot be expressed more truly.
Disclosure of Invention
(one) solving the technical problems
Aiming at the defects of the prior art, the invention provides a grid-type geoid refinement model construction method of geographic boundary constraint, which solves the problem that the elevation datum plane of a complex reservoir cannot be expressed more truly.
(II) technical scheme
In order to achieve the above purpose, the present invention provides the following technical solutions: the construction method of the geonet-like geolevel refinement model with geographic boundary constraint specifically comprises the following steps:
step one, boundary establishment: basic geographic information data in the irregular water body area is collected and researched, the design flood level of the irregular water body is calculated by combining the topographic data and the hydrologic data, and the boundary line of the management protection range of the irregular water body is defined and used as the calculated geographic boundary range;
step two, measuring control point data: uniformly distributing GNSS control points within a geographic boundary range according to a 10km point distance, acquiring plane coordinates (B, L) of each point by adopting GNSS plane control measurement, and acquiring normal height H of each point by adopting a height control measurement mode;
step three, building a model: acquiring EGM2008 earth gravity field model data, forming a continuous geospatial grid surface which completely covers a geographic boundary range line area, and searching xi of each GNSS control point in the EGM2008 earth gravity field model data file by adopting a minimum distance proximity search algorithm according to the measurement data of each GNSS control point GM By means of the data normalization removal process,calculating the elevation anomaly residual value of each GNSS control point;
step four, clipping constraint: according to the calculation result of the step three, a continuous grid elevation abnormal residual model is realized by using a Kriging interpolation method of an exponential variation function, and then cutting constraint is carried out according to the management protection range boundary line defined in the step one as a geographic boundary line, so as to obtain the continuous grid model of the elevation abnormal residual;
step five, generating a model: and (3) carrying out superposition recovery operation on the continuous grid model of the elevation anomaly residual error obtained in the step four and the gravity field model grid model to generate an elevation anomaly continuous grid model of the region.
The invention is further provided with: and step three, the EGM2008 earth gravity field model data obtained in the step three comprises Gao Chengyi constant values with the spatial resolution of 5 'x 5',2.5 'x 2.5' and 1 'x 1', and are stored in a discrete point data file format.
Wherein the file format is as follows:
first row: description dec
Second row: maximum column number NumCols; maximum number of rows NumRows
Third row: minimum latitude MinLat; maximum latitude MaxLat;
fourth row: minimum longitude MinLon; maximum longitude MaxLon;
fifth line: minimum value MinG, maximum value MaxG;
sixth line: G1G 1
Seventh row: G2G 2
………。
The invention is further provided with: the construction method of the geospatial grid surface in the third step specifically comprises the following steps: and calculating longitude max (lon), longitude min (lon), latitude max (lat) and latitude min (lat) according to the geographic boundary range, setting a grid distance d by taking (min (lon), min (lat)) at the left side of a view provided by the gravity field model data as a starting point, and extending grid lines in the east and north directions to form a continuous geographic space grid surface which completely covers the geographic boundary range line area.
The invention is further provided with: the calculation formula of the elevation abnormal residual value of the control point in the third step is as follows:
ξ C =H-h-ξ GM 。
the invention is further provided with: the implementation method of the continuous grid elevation anomaly residual error model in the step four is as follows:
a. suppose that the variable ζ is studied in the geospatial grid plane C (x) Attribute value ζ at point xi e a (i=1, 2, … …, n) C (x i ) Then the attribute value xi at the point x0 e A to be calculated C (x 0 ) And (2) the Kriging interpolation result xi * C (x 0 ) Is the known sample point attribute value ζ C (x i ) A weighted sum of (i=1, 2, … …, n), i.e.:
wherein lambda is i Is a pending weight coefficient;
b. on the premise of no bias, the variance is minimum, and the coefficient lambda of the weight to be determined is obtained i Is set of equations:
the above-described matrix form of the system of equations is represented as follows:
c. according to two conditions of the second order stationarity assumption, the following relationship exists between the covariance function and the variation function:
C(i,j)=C(0)-γ(Dis(i,j))
wherein, gamma (Dis (i, j)) is a variation function corresponding to the zoned variable when the distance is Dis (i, j), C (0) is a covariance value corresponding to the zoned variable when the distance is 0, and C (i, j) is a covariance value corresponding to the zoned variable when the distance is Dis (i, j);
d. the variation function is expressed as:
e. each cell is assigned a different independent value ζ * C (x 0 ) Obtaining a grid model with continuous elevation anomaly residual values in a research area, and calculating the elevation anomaly residual values of each cell according to the following calculation formula:
ξ * C (x 0 )=λ 1 ξ C (x 1 )+λ 2 ξ C (x 2 )+λ 3 ξ C (x 3 )+λ 4 ξ C (x 4 );
f. establishing a geographic boundary constraint matrix according to the geographic boundary range line and the grid surface, when judging whether to cross the grid boundary line according to the geographic boundary line, assigning 1 when crossing, assigning 0 when not crossing,
(III) beneficial effects
The invention provides a grid-type geoid refinement model construction method for geographic boundary constraint. The beneficial effects are as follows:
(1) According to the geodetic level refinement model construction method of geographic boundary constraint, when the geodetic level model is constructed, complex reservoir topographic data and hydrologic data are fused, the constructed continuous geodetic level model has applicability to geographic space based on the defined reservoir management and protection range as boundary constraint conditions, infinite space expansion caused by a single mathematical model fitting method is avoided, model errors are uncertain, and meanwhile, removal and recovery functions are realized easily through adoption of space addition and subtraction operation of the continuous geodetic model.
(2) According to the geonet-type geolike level refinement model construction method of geographic boundary constraint, the geonet-type model can be formed by simulating and inverting the geographic discrete data such as gravity field coefficients, elevation anomaly residual errors and the like into a continuous geonet-type model by adopting a Kriging interpolation method, and the model is stored in a grid file, so that the model can be well coupled with the conventional GNSS observation equipment, and the normal height of the grid point can be conveniently obtained in real time.
(3) According to the geonet-type geoid refinement model construction method based on geographic boundary constraint, the regional index variation function is used as a weighting factor by the Kriging interpolation method based on the index variation function, influences of reservoir management range, topography factors and the like are considered, smooth transition of an elevation reference model of an irregular reservoir management range is realized, and reality conditions are simulated more truly and accurately.
Drawings
FIG. 1 is a schematic flow chart of the present invention;
FIG. 2 is a schematic view of the division of the space grid surface according to the present invention;
FIG. 3 is a schematic view of a continuous grid model of the present invention;
FIG. 4 is a GNSS point location distribution diagram of a Danjiangkou reservoir area in accordance with an embodiment of the present invention;
FIG. 5 is a chart showing the precision of the level surface of the Danjiangkou reservoir similar to the earth in the embodiment of the invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Referring to fig. 1-5, the embodiment of the present invention provides a technical solution: the construction method of the geonet-like geolevel refinement model with geographic boundary constraint specifically comprises the following steps:
step one, boundary establishment: basic geographic information data in the irregular water body area is collected and researched, the design flood level of the irregular water body is calculated by combining the topographic data and the hydrologic data, and the boundary line of the management protection range of the irregular water body is defined and used as the calculated geographic boundary range;
step two, measuring control point data: uniformly distributing GNSS control points within a geographic boundary range according to a 10km point distance, acquiring plane coordinates (B, L) of each point by adopting GNSS plane control measurement, and acquiring normal height H of each point by adopting a height control measurement mode;
step three, building a model: acquiring global ultra-high-order gravity field model data of EGM2008, wherein the model is global ultra-high-order gravity field model data which can provide 2159 times and 2190-order spherical harmonic coefficients by combining ground gravity, satellite height measurement and satellite gravity, can provide Gao Chengyi constant values with spatial resolution of 5 'x 5',2.5 'x 2.5',1 'x 1', and is stored in a discrete point data file format, wherein the file format is as follows:
first row: description dec
Second row: maximum column number NumCols; maximum number of rows NumRows
Third row: minimum latitude MinLat; maximum latitude MaxLat;
fourth row: minimum longitude MinLon; maximum longitude MaxLon;
fifth line: minimum value MinG, maximum value MaxG;
sixth line: G1G 1
Seventh row: G2G 2
………;
Forming a continuous geospatial grid surface which completely covers a geographical boundary range line area, specifically, as shown in fig. 2, calculating longitude max (lon), longitude min (lon), latitude max (lat) and latitude min (lat) according to the geographical boundary range, taking (min (lon), min (lat)) at the left side of a view provided by gravity field model data as a starting point, setting grid spacing d, extending grid lines to the east and north directions, and forming the continuous geospatial grid surface which completely covers the geographical boundary range line area;
according to the measurement data of each GNSS control point, wherein the measurement data comprises the plane coordinates (B, L) of the GNSS control point, the ground height H and the normal height H, and the minimum distance proximity search algorithm is adopted to search each GNSS control point in the EGM2008 earth gravity field model data fileIs xi of (2) GM Through data normalization removal processing, the elevation anomaly residual value of each GNSS control point is calculated, and a specific calculation formula of the elevation anomaly residual value is as follows:
ξ C =H-h-ξ GM ;
step four, clipping constraint: according to the calculation result of the step three, a continuous grid elevation anomaly residual model is realized by using a kriging interpolation method of an exponential variation function, and the realization method is as follows:
a. suppose that the variable ζ is studied in the geospatial grid plane C (x) Attribute value ζ at point xi e a (i=1, 2, … …, n) C (x i ) Then the attribute value xi at the point x0 e A to be calculated C (x 0 ) And (2) the Kriging interpolation result xi * C (x 0 ) Is the known sample point attribute value ζ C (x i ) A weighted sum of (i=1, 2, … …, n), i.e.:
wherein lambda is i Is a pending weight coefficient;
b. taking unbiased premise, the variance is minimum, and obtaining the undetermined weight coefficient lambda according to the step a i Is set of equations:
the above-described matrix form of the system of equations is represented as follows:
c. according to two conditions of the second order stationarity assumption, the following relationship exists between the covariance function and the variation function:
C(i,j)=C(0)-γ(Dis(i,j))
wherein, gamma (Dis (i, j)) is a variation function corresponding to the zoned variable when the distance is Dis (i, j), C (0) is a covariance value corresponding to the zoned variable when the distance is 0, and C (i, j) is a covariance value corresponding to the zoned variable when the distance is Dis (i, j);
d. the variation function in step c is expressed by the following formula:
e. in connection with the formula presented in step a, each cell is assigned a different independent value ζ * C (x 0 ) The mesh model with continuous elevation anomaly residual values in the research area is obtained, as shown in fig. 3 (each cell is endowed with different elevation anomaly residual values, the colored cells are distributed over the whole mesh), and the elevation anomaly residual values of each cell are calculated according to the following calculation formula:
ξ * C (x 0 )=λ 1 ξ C (x 1 )+λ 2 ξ C (x 2 )+λ 3 ξ C (x 3 )+λ 4 ξ C (x 4 );
f. establishing a geographic boundary constraint matrix according to the geographic boundary range line and the grid surface, when judging whether to cross the grid boundary line according to the geographic boundary line, assigning 1 when crossing, assigning 0 when not crossing,
then cutting and restraining according to the management protection range boundary line defined in the first step as a geographic boundary line to obtain a continuous grid model of the elevation abnormal residual error;
step five, generating a model: and (3) carrying out superposition recovery operation on the continuous grid model of the elevation anomaly residual error obtained in the step four and the gravity field model grid model to generate an elevation anomaly continuous grid model of the region.
Examples: selecting a reservoir area of a Danjiang mouth reservoir, wherein the reservoir has more branches and larger topographic relief, and the total area of the area is about 2000km 2 This regionAnd (3) carrying out domain setting and measuring 40C-level GNSS control points, carrying out three-level joint measurement on each control point, wherein the highest point is about 355m in height, the lowest point is Cheng Yaowei m in height, different fitting points are selected according to the distribution condition of the C-level GNSS control points, different interpolation fitting methods are adopted to construct a model, and the rest points are selected as checking points to carry out precision statistics, as shown in fig. 4 and fig. 5, the method shows better precision in constructing a DAN river mouth reservoir elevation reference model, and can greatly improve the precision and applicability of the irregular large reservoir geoid model.
Although embodiments of the present invention have been shown and described, it will be understood by those skilled in the art that various changes, modifications, substitutions and alterations can be made therein without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.
Claims (6)
1. The method for constructing the geonet-like geolevel refinement model of the geographic boundary constraint is characterized by comprising the following steps of: the method specifically comprises the following steps:
step one, boundary establishment: basic geographic information data in the irregular water body area is collected and researched, the design flood level of the irregular water body is calculated by combining the topographic data and the hydrologic data, and the boundary line of the management protection range of the irregular water body is defined and used as the calculated geographic boundary range;
step two, measuring control point data: uniformly distributing GNSS control points within a geographic boundary range according to a 10km point distance, acquiring plane coordinates (B, L) of each point by adopting GNSS plane control measurement, and acquiring normal height H of each point by adopting a height control measurement mode;
step three, building a model: acquiring EGM2008 earth gravity field model data, forming a continuous geospatial grid surface which completely covers a geographic boundary range line area, and searching xi of each GNSS control point in the EGM2008 earth gravity field model data file by adopting a minimum distance proximity search algorithm according to the measurement data of each GNSS control point GM Calculating the elevation abnormality of each GNSS control point through data normalization removal processingResidual value;
step four, clipping constraint: according to the calculation result of the step three, a continuous grid elevation abnormal residual model is realized by using a Kriging interpolation method of an exponential variation function, and then cutting constraint is carried out according to the management protection range boundary line defined in the step one as a geographic boundary line, so as to obtain the continuous grid model of the elevation abnormal residual;
step five, generating a model: and (3) carrying out superposition recovery operation on the continuous grid model of the elevation anomaly residual error obtained in the step four and the gravity field model grid model to generate an elevation anomaly continuous grid model of the region.
2. The geodetic refinement model construction method of geographic boundary constraints of claim 1, wherein: and step three, the EGM2008 earth gravity field model data obtained in the step three comprises Gao Chengyi constant values with the spatial resolution of 5 'x 5',2.5 'x 2.5' and 1 'x 1', and are stored in a discrete point data file format.
3. The geodetic refinement model construction method of geographic boundary constraints of claim 1, wherein: the construction method of the geospatial grid surface in the third step specifically comprises the following steps: and calculating longitude max (lon), longitude min (lon), latitude max (lat) and latitude min (lat) according to the geographic boundary range, setting a grid distance d by taking (min (lon), min (lat)) at the left side of a view provided by the gravity field model data as a starting point, and extending grid lines in the east and north directions to form a continuous geographic space grid surface which completely covers the geographic boundary range line area.
4. The geodetic refinement model construction method of geographic boundary constraints of claim 1, wherein: the calculation formula of the elevation abnormal residual value of the control point in the third step is as follows:
ξ C =H-h-ξ GM 。
5. the geodetic refinement model construction method of geographic boundary constraints of claim 1, wherein: the implementation method of the continuous grid elevation anomaly residual error model in the step four is as follows:
a. suppose that the variable ζ is studied in the geospatial grid plane C (x) Attribute value ζ at point xi e a (i=1, 2, … …, n) C (x i ) Then the attribute value xi at the point x0 e A to be calculated C (x 0 ) And (2) the Kriging interpolation result xi * C (x 0 ) Is the known sample point attribute value ζ C (x i ) A weighted sum of (i=1, 2, … …, n), i.e.:
wherein lambda is i Is a pending weight coefficient;
b. on the premise of no bias, the variance is minimum, and the coefficient lambda of the weight to be determined is obtained i Is set of equations:
the above-described matrix form of the system of equations is represented as follows:
c. according to two conditions of the second order stationarity assumption, the following relationship exists between the covariance function and the variation function:
C(i,j)=C(0)-γ(Dis(i,j))
wherein, gamma (Dis (i, j)) is a variation function corresponding to the zoned variable when the distance is Dis (i, j), C (0) is a covariance value corresponding to the zoned variable when the distance is 0, and C (i, j) is a covariance value corresponding to the zoned variable when the distance is Dis (i, j);
d. the variation function is expressed as:
e. each cell is assigned a different independent value ζ * C (x 0 ) Obtaining a grid model with continuous elevation abnormal residual values in a research area, and calculating the elevation abnormal residual values of each cell;
f. establishing a geographic boundary constraint matrix according to the geographic boundary range line and the grid surface, when judging whether to cross the grid boundary line according to the geographic boundary line, assigning 1 when crossing, assigning 0 when not crossing,
6. the geodetic refinement model construction method of geographic boundary constraints of claim 7, wherein: the calculation formula of the elevation anomaly residual value of each cell is as follows:
ξ * C (x 0 )=λ 1 ξ C (x 1 )+λ 2 ξ C (x 2 )+λ 3 ξ C (x 3 )+λ 4 ξ C (x 4 )。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310010593.9A CN116402974A (en) | 2023-01-05 | 2023-01-05 | Geonet-like geolevel refinement model construction method for geographic boundary constraint |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310010593.9A CN116402974A (en) | 2023-01-05 | 2023-01-05 | Geonet-like geolevel refinement model construction method for geographic boundary constraint |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116402974A true CN116402974A (en) | 2023-07-07 |
Family
ID=87008131
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310010593.9A Pending CN116402974A (en) | 2023-01-05 | 2023-01-05 | Geonet-like geolevel refinement model construction method for geographic boundary constraint |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116402974A (en) |
-
2023
- 2023-01-05 CN CN202310010593.9A patent/CN116402974A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kumler | An intensive comparison of triangulated irregular networks (TINs) and digital elevation models (DEMs) | |
Panhalkar et al. | Assessment of spatial interpolation techniques for river bathymetry generation of Panchganga River basin using geoinformatic techniques | |
CN110544305B (en) | Terrain steeple ridge line information fusion method constructed by DEM (digital elevation model) oriented to regular grids | |
CN110335355A (en) | A kind of high automatic calculating method of large-scale shallow water lake | |
CN112242003B (en) | City sub-catchment area division method considering land type and flow direction | |
CN117195603B (en) | Flood disaster deduction method, equipment and medium based on high-resolution remote sensing elements | |
Neupane et al. | Landslide susceptibility analysis using frequency ratio and weight of evidence approaches along the Lakhandehi Khola watershed in the Sarlahi District, southern Nepal | |
CN108416082A (en) | A kind of marine site pour point external disturbance gravitation horizontal component without unusual computational methods | |
CN116402974A (en) | Geonet-like geolevel refinement model construction method for geographic boundary constraint | |
Radoš et al. | Morfometrijske značajke šireg područja Duvanjskog polja, Bosna i Hercegovina | |
CN115018973A (en) | Low-altitude unmanned-machine point cloud modeling precision target-free evaluation method | |
Gerey et al. | Groundwater Single‐and Multiobjective Optimization Using Harris Hawks and Multiobjective Billiards‐Inspired Algorithm | |
CN115112090A (en) | Mountain area digital elevation datum establishing method based on few GNSS leveling points | |
CN112733215A (en) | Earthwork engineering quantity calculation method based on improved constraint irregular triangulation network construction method | |
Hou et al. | Effects of sampling and interpolation methods on accuracy of extracted watershed features | |
Abdurakhmonov et al. | Application of GIS technologies in the methods of geodetic work in the mapping of demographic processes | |
Yakubu et al. | Modelling Uncertainties in Differential Global Positioning System Dataset | |
Beaubouef et al. | Rough set based uncertainty management for spatial databases and geographical information systems | |
TWI786874B (en) | Method of the digital grid model and system thereof | |
Refadah | Digital Elevation Model and Drainage Basin Delineation | |
Prakasam et al. | Evaluation of geomorphic resources using GIS technology: a case study of selected villages in Ausgram Block, Burdwam District, West Bengal, India | |
Singh et al. | River network identification using Remote Sensing and GIS | |
Abdel-Bary Ebrahim et al. | Planning road networks in new cities using GIS: The case of New Sohag, Egypt | |
Lomsdalen | Viewscapes and cosmology in the prehistoric temples of Malta | |
CN107194820B (en) | Method for determining representative area and representative degree of crop water demand monitoring site |
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 |