CN111881333A - Method for determining gravity center of regional soil erosion space - Google Patents

Method for determining gravity center of regional soil erosion space Download PDF

Info

Publication number
CN111881333A
CN111881333A CN202010704812.XA CN202010704812A CN111881333A CN 111881333 A CN111881333 A CN 111881333A CN 202010704812 A CN202010704812 A CN 202010704812A CN 111881333 A CN111881333 A CN 111881333A
Authority
CN
China
Prior art keywords
erosion
soil erosion
space
intensity
soil
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
Application number
CN202010704812.XA
Other languages
Chinese (zh)
Inventor
王友胜
刘冰
马海宽
成晨
张晓明
王昭艳
张磊
解刚
殷小琳
赵阳
辛艳
张永娥
邢先双
董明明
杜鹏飞
焦剑
郭乾坤
任正龑
刘卉芳
孟琳
赵莹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN202010704812.XA priority Critical patent/CN111881333A/en
Publication of CN111881333A publication Critical patent/CN111881333A/en
Pending legal-status Critical Current

Links

Images

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/90—Details of database functions independent of the retrieved data types
    • G06F16/904—Browsing; Visualisation therefor
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

The invention discloses a method for determining the gravity center of a regional soil erosion space, which comprises the following steps: collecting soil erosion intensity grading space grid data in a measuring area and standardizing an erosion intensity code of the soil erosion intensity grading space grid data; normalizing the reference coordinate system thereof to a geographical coordinate system; converting soil erosion intensity hierarchical spatial grid data in a geographic coordinate system into an ASCII format; respectively calculating the gravity center coordinates of the regional soil erosion space according to a formula; and converting the calculated barycentric coordinates into a shp-format point file in ArcGIS, superposing the shp-format point file with the base map of the soil erosion intensity classification space grid data of different time periods, and performing mapping display to compare the barycentric of the soil erosion space of the region of different time periods. The invention can provide effective support for fast and quantitatively mastering the macroscopic level and the erosion distribution pattern. The spatial gravity center migration of regional erosion can be quantitatively explained, and scientific reference is provided for mastering the change of the soil erosion pattern.

Description

Method for determining gravity center of regional soil erosion space
Technical Field
The invention relates to the technical field of environmental protection, in particular to a method for determining the gravity center of a regional soil erosion space, and relates to the influence of different erosion strengths on the gravity center of the soil erosion space.
Background
The spatial center of gravity is a center of gravity point obtained by comprehensively considering the magnitude of the acting force in each direction, the moving direction of the center of gravity point is represented as moving towards the direction with the large acting force, and the moving direction reflects the direction of the changing pattern. Its spatial location is mainly influenced by two major factors, namely the geographical location and the attribute change of different directions.
The soil erosion intensity map comprehensively reflects the spatial distribution pattern of different erosion intensities. Because the influence of different erosion strengths on the surface environment of the earth is different, the main manifestations are that the soil loss amount of the high-level erosion strength is large, the nutrient loss is rapid, and the soil fertility is more obviously reduced. The soil erosion space heterogeneity is obvious and the difference of the erosion intensity among the areas is obvious under the influence of various aspects such as landform, soil type, surface vegetation, artificial farming, rainfall and the like.
In the prior art, only qualitative description exists for analyzing the gravity center of the regional soil erosion space, and quantitative representation is lacked. The method for calculating the gravity center of rainfall erosion force exists in the analysis of soil erosion influence factors, but the soil erosion strength and the erosion force belong to two different types of elements, the numerical value of the rainfall erosion force is a continuous numerical value, but for the soil erosion strength, the grading of the soil erosion strength relates to 6 grades of strength grades, and the influence of different strengths on erosion is different, so that the influence of different strengths needs to be reflected in the gravity center calculation, the original calculation formula is not applicable to the calculation of the gravity center of the regional soil erosion space, and the requirement of reflection of the gravity centers of different erosion strengths cannot be met.
Disclosure of Invention
The invention aims to provide a method for determining the gravity center of a regional soil erosion space by considering different erosion intensity weights.
The invention is realized by the following technical scheme:
a method for determining the gravity center of a regional soil erosion space,
the method comprises the following steps: collecting soil erosion intensity grading space grid data in a measuring area and standardizing an erosion intensity code of the soil erosion intensity grading space grid data: the respective soil erosion strength levels were: micro erosion, mild erosion, moderate erosion, strong erosion, extremely strong erosion and severe erosion, wherein erosion strength codes corresponding to all levels are respectively specified as 11, 12, 13, 14, 15 and 16;
step two: checking a reference coordinate system of the soil erosion intensity grading space grid data, and normalizing the reference coordinate system into a geographical coordinate system:
step three: converting soil erosion intensity hierarchical spatial grid data in a geographic coordinate system into an ASCII format;
step four: and (3) calculating the gravity center of the soil erosion space in the measurement area: respectively calculating the barycentric coordinates of the regional soil erosion space according to the formula (1) and the formula (2): including a longitude coordinate and a latitude coordinate,
Figure BDA0002594326490000021
wherein, the weight coefficients are respectively:
Figure BDA0002594326490000022
satisfy the requirement of
Figure BDA0002594326490000023
In the formula:
Figure BDA0002594326490000024
longitude coordinates representing the center of gravity of the erosion space in degrees (decimal); k is a radical ofiIs shown asThe weighting coefficients of the i grids; j represents the erosion intensity code of the ith grid, and j is equal to 11, 12, 13, 14, 15 and 16, which respectively represents micro erosion, light erosion, medium erosion, strong erosion and severe erosion; eiRepresents the longitude of the center of the ith grid cell in degrees (decimal); n represents the total number of grids;
Figure BDA0002594326490000025
in the formula:
Figure BDA0002594326490000026
latitude coordinate representing the center of gravity of the erosion space in deg. (decimal), NiLatitude coordinates representing the center of the ith grid cell in degrees (decimal); k is a radical ofiAnd n is as defined above.
Further scheme, still include step five: and converting the calculated barycentric coordinates into a shp-format point file in ArcGIS, superposing the shp-format point file with a base map of the soil erosion intensity classification space grid data, and performing drawing display.
Further, step five: and converting the calculated barycentric coordinates into a shp-format point file in ArcGIS, superposing the shp-format point file with the base map of the soil erosion intensity classification space grid data of different time periods, and performing mapping display to compare the barycentric of the soil erosion space of the region of different time periods.
Further, in the step one, the soil erosion intensity grading space grid data is in a TIFF format, and different soil erosion intensity grades are represented by different color blocks.
The invention has the beneficial effects that:
can provide effective support for fast and quantitatively mastering the macroscopic level and the erosion distribution pattern. On the basis of the first-stage soil erosion center, the change analysis of the multi-stage soil erosion space center can quantitatively explain the spatial center migration of regional erosion, and scientific reference is provided for mastering the change of the soil erosion pattern.
Drawings
FIG. 1 is an ASCII format file schematic;
FIG. 2 is a schematic view of the center of gravity of soil erosion in 2010 area (the bottom view is a 2010 soil erosion intensity distribution diagram, and asterisks indicate the position of the center of gravity);
FIG. 3 is a graph showing changes of the soil erosion space barycenter in 2010-2020 region (the bottom graph is a soil erosion intensity distribution graph in 2020, and asterisks indicate the barycenter positions).
Detailed Description
The invention relates to a method for determining the gravity center of a regional soil erosion space, which is implemented by combining the following specific operation steps:
the method comprises the following steps: collecting soil erosion intensity grading space grid data in a measuring area and standardizing intensity codes of the soil erosion intensity grading space grid data: the respective soil erosion strength levels were: micro erosion, mild erosion, moderate erosion, strong erosion, extremely strong erosion and severe erosion, wherein the strength codes corresponding to all levels are respectively specified as 11, 12, 13, 14, 15 and 16; for codes that are not normalized, they are normalized to six-level codes of 11-16.
Step two: and checking a reference Coordinate system of the soil erosion intensity grading space grid data, and if the Coordinate system is a projection Coordinate system or a Coordinate system is absent, normalizing the Coordinate system into a Geographic Coordinate system (Geographic Coordinate Systems).
Step three: converting soil erosion intensity space grid data in tiff format under a geographic coordinate system into ASCII format; the converted ASCII format is shown in figure 1.
Step four: and (3) calculating the gravity center of the soil erosion space in the measurement area: and (3) respectively calculating the spatial gravity center (longitude and latitude) of soil erosion according to the formula (1) and the formula (2) according to the coordinates of the starting point in the head file and the grid size of the grid file. The xllcorner and ylcerner in the header file represent the coordinates of the lower left corner of the raster file, and the cellsize represents the raster size (°).
Figure BDA0002594326490000041
Wherein, the weight coefficients are respectively:
Figure BDA0002594326490000042
satisfy the requirement of
Figure BDA0002594326490000043
In the formula:
Figure BDA0002594326490000044
longitude coordinates representing the center of gravity of the erosion space in degrees (decimal); k is a radical ofiA weight coefficient representing the ith grid; j represents the erosion intensity code of the ith grid, and j is equal to 11, 12, 13, 14, 15 and 16, which respectively represents micro erosion, light erosion, medium erosion, strong erosion and severe erosion; eiRepresents the longitude of the center of the ith grid cell in degrees (decimal); n represents the total number of grids.
Figure BDA0002594326490000045
In the formula:
Figure BDA0002594326490000046
latitude coordinate representing the center of gravity of the erosion space in deg. (decimal), NiLatitude coordinates representing the center of the ith grid cell in degrees (decimal); n, kiThe meaning is the same as above.
Fifthly, the calculated barycentric coordinates (longitude E and latitude N) are converted into a shp-format point file in ArcGIS, and the shp-format point file is overlaid with base maps of different soil erosion strengths for graphical display (fig. 2 and fig. 3).
The calculation result of the center of gravity in 2010: e113.5992 ° N41.6614 °
The calculation result of the center of gravity in 2020: e117.4354 ° N40.8273 °
As can be seen from the data of the gravity centers, the spatial gravity center of the soil erosion strength in 2010-2020 shifts to the southeast, 3.8362 degrees (117.4354-113.5992 degrees) to the east and 0.8341 degrees (41.6614-40.8273 degrees) to the south.
The method for reflecting the centers of gravity of different soil erosion spaces in the region is constructed, the quantitative analysis of the positions of the centers of gravity of the soil erosion spaces in the region can be realized, the migration change of the positions of the centers of gravity of the soil erosion spaces in different years can be further analyzed, the change direction of the soil erosion strength of the region in the space can be better, the space change of the soil erosion strength of the region can be better mastered, scientific support is provided for the arrangement of macro-level water and soil conservation measures and the establishment of regional water and soil conservation planning, and the soil erosion can be more accurately prevented and treated.

Claims (4)

1. A method for determining the gravity center of a regional soil erosion space is characterized by comprising the following steps: the method comprises the following steps:
the method comprises the following steps: collecting soil erosion intensity grading space grid data in a measuring area and standardizing an erosion intensity code of the soil erosion intensity grading space grid data: the respective soil erosion strength levels were: micro erosion, mild erosion, moderate erosion, strong erosion, extremely strong erosion and severe erosion, wherein erosion strength codes corresponding to all levels are respectively specified as 11, 12, 13, 14, 15 and 16;
step two: checking a reference coordinate system of the soil erosion intensity grading space grid data, and normalizing the reference coordinate system into a geographical coordinate system:
step three: converting soil erosion intensity hierarchical spatial grid data in a geographic coordinate system into an ASCII format;
step four: and (3) calculating the gravity center of the soil erosion space in the measurement area: respectively calculating the barycentric coordinates of the regional soil erosion space according to the formula (1) and the formula (2): including a longitude coordinate and a latitude coordinate,
Figure FDA0002594326480000011
wherein, the weight coefficients are respectively:
Figure FDA0002594326480000012
in the formula:
Figure FDA0002594326480000013
a longitudinal coordinate representing the center of gravity of the erosion space in °; k is a radical ofiA weight coefficient representing the ith grid; j represents the erosion intensity code of the ith grid, and j is equal to 11, 12, 13, 14, 15 and 16, which respectively represents micro erosion, light erosion, medium erosion, strong erosion and severe erosion; eiRepresents the longitude of the center of the ith grid cell in °; n represents the total number of grids;
Figure FDA0002594326480000014
in the formula:
Figure FDA0002594326480000015
latitude coordinate representing the center of gravity of the erosion space in degrees, NiA latitude coordinate representing the center of the ith grid cell in degrees; k is a radical ofiAnd n is as defined above.
2. The method for determining the gravity center of the regional soil erosion space according to claim 1, wherein the method comprises the following steps: further comprises the following steps: and converting the calculated barycentric coordinates into a shp-format point file in ArcGIS, superposing the shp-format point file with a base map of the soil erosion intensity classification space grid data, and performing drawing display.
3. The method for determining the gravity center of the regional soil erosion space according to claim 2, wherein the method comprises the following steps: step five: and converting the calculated barycentric coordinates into a shp-format point file in ArcGIS, superposing the shp-format point file with the base map of the soil erosion intensity classification space grid data of different time periods, and performing mapping display to compare the barycentric of the soil erosion space of the region of different time periods.
4. The method for determining the gravity center of the regional soil erosion space according to claim 1, wherein the method comprises the following steps: in the first step, the soil erosion intensity grading space grid data is in a TIFF format, and different soil erosion intensity grades are represented by different color blocks.
CN202010704812.XA 2020-07-21 2020-07-21 Method for determining gravity center of regional soil erosion space Pending CN111881333A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010704812.XA CN111881333A (en) 2020-07-21 2020-07-21 Method for determining gravity center of regional soil erosion space

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010704812.XA CN111881333A (en) 2020-07-21 2020-07-21 Method for determining gravity center of regional soil erosion space

Publications (1)

Publication Number Publication Date
CN111881333A true CN111881333A (en) 2020-11-03

Family

ID=73155777

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010704812.XA Pending CN111881333A (en) 2020-07-21 2020-07-21 Method for determining gravity center of regional soil erosion space

Country Status (1)

Country Link
CN (1) CN111881333A (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106599601A (en) * 2016-12-29 2017-04-26 中国科学院遥感与数字地球研究所 Remote sensing assessment method and system for ecosystem vulnerability
CN107886216A (en) * 2017-10-12 2018-04-06 天津大学 A kind of land carrying capacity analysis method based on Evaluation of Land Use Suitability
CN108229859A (en) * 2018-02-09 2018-06-29 中国环境科学研究院 A kind of method and system of the key area of determining bio-diversity conservation

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106599601A (en) * 2016-12-29 2017-04-26 中国科学院遥感与数字地球研究所 Remote sensing assessment method and system for ecosystem vulnerability
CN107886216A (en) * 2017-10-12 2018-04-06 天津大学 A kind of land carrying capacity analysis method based on Evaluation of Land Use Suitability
CN108229859A (en) * 2018-02-09 2018-06-29 中国环境科学研究院 A kind of method and system of the key area of determining bio-diversity conservation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王尧: "喀斯特地区土壤侵蚀模拟研究", 《中国博士学位论文全文数据库农业科技辑》 *

Similar Documents

Publication Publication Date Title
CN111666918B (en) Coastline change identification method based on multiple factors
Verdonck et al. Influence of neighbourhood information on ‘Local Climate Zone’mapping in heterogeneous cities
CN113283802B (en) Landslide risk assessment method for complicated and difficult mountain areas
Bunting et al. Modelling pollen dispersal and deposition using HUMPOL software, including simulating windroses and irregular lakes
JP4545219B1 (en) Analysis method of topographic change using topographic image and program thereof
Klouček et al. How does data accuracy influence the reliability of digital viewshed models? A case study with wind turbines
CN106780586B (en) A kind of solar energy potential evaluation method based on ground laser point cloud
CN101582215A (en) Multi-stage nine-grid locating method of spatial information
Tin‐Seong Integrating GIS and remote sensing techniques for urban land‐cover and land‐use analysis
Pili et al. A GIS tool for the calculation of solar irradiation on buildings at the urban scale, based on Italian standards
CN113743027B (en) Method and device for drawing wind resource map based on CFD technology
Sangawongse et al. Urban growth and land cover change in Chiang Mai and Taipei: results from the SLEUTH model
Khahro et al. Data preparation for GIS based land suitability modelling: A stepped approach
Sengupta et al. Assessment of population exposure and risk zones due to air pollution using the geographical information system
Williams Implementing LESA on a geographic information system: A case study
CN111881333A (en) Method for determining gravity center of regional soil erosion space
KR20050063615A (en) Method for providing surface roughness in geographic information system
CN115712978A (en) Site selection data processing method and system based on numerical simulation and population data
CN109376638A (en) A kind of literary ground rate calculation method based on remote sensing images and GIS-Geographic Information System
Al-Ahmadi et al. Validating spatial patterns of urban growth from a cellular automata model
Heyns Reduced target-resolution strategy for rapid multi-observer site location optimisation
CN107908858B (en) Coastal wetland waterfowl movement path simulation method based on GIS and RS
CN117523063B (en) Rendering method of high-performance ray tracing simulation based on three-dimensional space
CN115034435B (en) Target observation index forecasting method based on numerical model, storage medium and terminal
Cloete Optimizing visibility analysis for visual impact assessment

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20201103