CN106372599A - Method and system for extracting silt arrester in water and soil retaining period - Google Patents

Method and system for extracting silt arrester in water and soil retaining period Download PDF

Info

Publication number
CN106372599A
CN106372599A CN201610787008.6A CN201610787008A CN106372599A CN 106372599 A CN106372599 A CN 106372599A CN 201610787008 A CN201610787008 A CN 201610787008A CN 106372599 A CN106372599 A CN 106372599A
Authority
CN
China
Prior art keywords
water
water body
region
silt arrester
dam
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.)
Granted
Application number
CN201610787008.6A
Other languages
Chinese (zh)
Other versions
CN106372599B (en
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.)
Beijing Datum Science & Technology Development Co Ltd
Soil And Water Conservation Monitoring Center Of Ministry Of Water Resources
Original Assignee
Beijing Datum Science & Technology Development Co Ltd
Soil And Water Conservation Monitoring Center Of Ministry Of Water Resources
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 Beijing Datum Science & Technology Development Co Ltd, Soil And Water Conservation Monitoring Center Of Ministry Of Water Resources filed Critical Beijing Datum Science & Technology Development Co Ltd
Priority to CN201610787008.6A priority Critical patent/CN106372599B/en
Publication of CN106372599A publication Critical patent/CN106372599A/en
Application granted granted Critical
Publication of CN106372599B publication Critical patent/CN106372599B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/13Satellite images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/25Determination of region of interest [ROI] or a volume of interest [VOI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/44Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/194Terrestrial scenes using hyperspectral data, i.e. more or other wavelengths than RGB

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Astronomy & Astrophysics (AREA)
  • Remote Sensing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)

Abstract

The invention provides a method and a system for extracting a silt arrester in a water and soil retaining period. The method and the system are for extracting the silt arrester in the water and soil retaining period. The method specifically comprises steps that remote sensing images are processed, and water body areas are extracted according to water body spectrum characteristics, space characteristics and geometric characteristics; a dam body position scope is determined according to water body and dam body space relationships; overlapping analysis on the water body area and the dam body position scope is carried out to acquire the information of the slit arrester. After the information of the slit arrester is acquired, the quantity and distribution conditions of slit arresters in a pre-determined preliminary exploration area can be acquired, statistics and field investigation of construction data archives can be carried out not in a manual mode, manpower and material resources are greatly saved, and statistics efficiency is effectively improved.

Description

A kind of extracting method of retaining water and soil phase silt arrester and system
Technical field
The present invention relates to territory investigation and water conservancy generally investigate technical field, more particularly to a kind of retaining water and soil phase silt arrester Extracting method and system.
Background technology
Silt arrester refers to the dam work building in the raceway groove at different levels of erosion area, built for the purpose of blocking mud alluvial plain Thing, it blocks the ground that mud becomes silted up on dam ground.In the raceway groove of basin, silt arrester or produce dam for the dam that alluvial plain produces.Loess is high Original area's soil layer is deep, and loess blazons, and has the spies such as homogeneous, short texture, water penetration are strong, easy disintegrating, dewatering and coagulating are fast Point, is good dam material, can gather materials on the spot.Block mud alluvial plain using silt arrester, for raising groove channel erosion datum level, anti- Control soil erosion, flood detention, block mud, alluvial plain (dam ground), reduce steady into yellow mud sand, the local production and living conditions of improvement, construction high yield The basic farmland that produces, aspect has highly important meaning to promote the local masses to shake off poverty and set out on the road to prosperity etc., is comprehensive treatment of small watersheds One important measures.
Wherein, silt arrester according to run time divided, be divided into retaining water and soil phase silt arrester and production run phase alluvial plain Dam two class.Retaining water and soil phase silt arrester refers to the silt arrester in flood retention water-retention period, the footpath that upstream is mainly collected by silt arrester at this moment Stream and silt are blocked in raceway groove, have the large area water surface in raceway groove.Production run phase silt arrester is primarily referred to as completing to block mud alluvial plain Silt arrester, the silt arrester in this period intercepts, through long-term, the silt that upstream sweeps away, and groove channel erosion datum level is raised, and forms dam Ground, no longer plays flood retention Water storage.
In order that government department can formulate water-and-soil conservation measures according to the concrete condition of silt arrester it is necessary to silt arrester Quantity and distribution counted, traditionally, silt arrester quantity and distribution situation are entered to building profile using manual type The mode of row statistics and on-site inspection obtains, and this mode needs to expend substantial amounts of manpower and materials, causes statistical efficiency poor.
Content of the invention
In view of this, the invention provides a kind of extracting method of retaining water and soil phase silt arrester and system, artificial to solve The poor problem of the statistical method efficiency of on-site inspection.
In order to solve the above problems, the invention discloses a kind of extracting method of retaining water and soil phase silt arrester, specifically include Step:
Remote sensing image is processed, according to water spectral feature, space characteristics and Extraction of Geometrical Features water outlet body region;
According to the spatial relationship of water body and dam body, dam body position range is confined;
Described water body region and described dam body position range are overlapped analyzing, obtain silt arrester information.
Optionally, described according to water spectral feature, space characteristics and Extraction of Geometrical Features water outlet body region, comprising:
Described water body region is extracted according to normalization water body index, digital elevation model and coverage of water.
Optionally, the described spatial relationship according to water body region and dam body is confined to dam body position range, comprising:
Region of interest is created based on closing operation of mathematical morphology principle;
Water-outlet body upstream termination is extracted in described region of interest according to preset upstream termination extraction algorithm;
Water-outlet body downstream end is extracted in described region of interest according to preset downstream end extraction algorithm;
Described dam body position is confined out according to the spatial relationship of described water body upstream termination, described water body downstream end and dam body Put scope.
Optionally, further comprise the steps of:
According to the data correction instruction of user input, described silt arrester information is optimized perfect.
Optionally, further comprise the steps of:
According to described silt arrester information, the quantity of the silt arrester in predetermined census area and distribution are counted.
Accordingly, in order to ensure the enforcement of said method, present invention also offers the carrying of a kind of retaining water and soil phase silt arrester Take system, specifically include:
Clean water withdraw module, for processing to remote sensing image, special according to water spectral feature, space characteristics and geometry Levy extraction water-outlet body region;
Dam body confines module, for being confined to dam body position range according to the spatial relationship of water body and dam body;
Information extraction modules, for being overlapped analyzing described water body region and described dam body position range, are become silted up Ground dam information.
Optionally, described Clean water withdraw module includes:
Water body area extracting unit, for extracting institute according to normalization water body index, digital elevation model and coverage of water State water body region.
Optionally, described dam body is confined module and is included:
Region of interest creating unit, for creating region of interest based on closing operation of mathematical morphology principle;
Water body upstream termination extraction unit, for extracting in described region of interest according to preset upstream termination extraction algorithm Water-outlet body upstream termination;
Water body downstream end extraction unit, for extracting in described region of interest according to preset downstream end extraction algorithm Water-outlet body downstream end;
Dam body position range confines unit, for the sky according to described water body upstream termination, water body downstream end and dam body Between relation confine out described dam body position range.
Optionally, also include:
Optimize module, for the data correction instruction according to user input, described silt arrester information has been optimized Kind.
Optionally, also include:
Information Statistics module, for and dividing to the quantity of the silt arrester in predetermined census area according to described silt arrester information Cloth is counted.
From technique scheme as can be seen that embodiments providing a kind of extraction side of retaining water and soil phase silt arrester Method and system, the method and system are used for the silt arrester being in the retaining water and soil phase is extracted, and specially remote sensing image are entered Row is processed, according to water spectral feature, space characteristics and Extraction of Geometrical Features water outlet body region;Space according to water body and dam body Relation pair dam body position range is confined;Water body region and dam body position range are overlapped analyzing, obtain silt arrester letter Breath.After obtaining silt arrester information by simple statistical can obtain in predetermined census area the quantity of silt arrester and point Cloth situation, and without manual type, construction profile is counted and on-site inspection, thus saving substantial amounts of manpower thing Power, can effectively improve statistical efficiency.
Brief description
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing Have technology description in required use accompanying drawing be briefly described it should be apparent that, drawings in the following description be only this Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, acceptable Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is a kind of flow chart of steps of the extracting method embodiment of retaining water and soil phase silt arrester of the present invention;
Fig. 2 is the flow chart of steps of the extracting method embodiment of another kind of retaining water and soil phase silt arrester of the present invention;
Fig. 2 a is the water body extracted region design sketch of the present invention;
Fig. 2 b is the design sketch of the region of interest in water body region of the present invention;
Fig. 2 c is the design sketch before the water body upstream termination extraction of the present invention;
Fig. 2 d is the design sketch after the water body upstream termination extraction of the present invention;
Fig. 2 e is the design sketch before the water body downstream end extraction of the present invention;
Fig. 2 f is the design sketch after the water body downstream end extraction of the present invention;
Fig. 2 g is the design sketch of the dam body position range of the present invention;
Fig. 2 h is the design sketch of the silt arrester information of the present invention;
Fig. 2 i is the actual extracting design sketch of the present invention;
Fig. 3 is a kind of structured flowchart of the extraction system embodiment of retaining water and soil phase silt arrester of the present invention;
Fig. 4 is the structured flowchart of the extraction system embodiment of another kind of retaining water and soil phase silt arrester of the present invention.
Specific embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete Site preparation description is it is clear that described embodiment is only a part of embodiment of the present invention, rather than whole embodiments.It is based on Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of not making creative work Embodiment, broadly falls into the scope of protection of the invention.
The silt arrester of production run phase due to Sediment Siltation, is connected dam body and dam, is soil property, base on remote sensing image Originally cannot see under the conditions of visually.And the silt arrester of retaining water and soil phase, dam body is connected with water body, and water body has on remote sensing image There is obvious feature.The application is exactly using the principle of Clean water withdraw, dam body to be extracted, and then extracts the retaining water and soil phase Complete silt arrester.Specific embodiment is as described below.
Embodiment one
Fig. 1 is a kind of flow chart of steps of the extracting method embodiment of retaining water and soil phase silt arrester of the present invention.
It is used for the silt arrester being in the retaining water and soil phase is carried with reference to the extracting method that shown in Fig. 1, the present embodiment provides Take, specifically include following steps:
S101: according to water spectral feature, space characteristics and Extraction of Geometrical Features water outlet body region.
By remote sensing images analysis, to be distributed in mountain area and the hills area of northern Loess Plateau more than silt arrester, Build in trench bottom;The water body of silt arrester and general water body have similar spectral signature, but the water surface is less, many strips Shape, and water body downstream is connected with dam body.Therefore, according to the water spectral feature of water body, space characteristics and geometric properties to pre- The remote sensing image determining census area is processed, and therefrom proposes water body region.
Specific to being using normalization water body index ndwi (normalized difference water in the present embodiment Index), digital elevation model dem and coverage of water extract this water body region.
Normalization water body index ndwi is the normalized ratio index based on green light band and near infrared band.This formula For:
Ndwi=(green-nir)/(green+nir)
In formula, green green light band;Nir near infrared band.
S102: dam body position range is confined according to the spatial relationship in water body region and dam body.
Water body region comprises water body upstream termination and water body downstream end, and wherein water body downstream end refers to the dam of silt arrester , near the part of this dam body, water body upstream termination is corresponding with water body downstream end for water body more than body, and away from dam body according to This principle is confined to the dam body position range residing for dam body.
S103: water body region and dam body position range are overlapped analyzing, obtain silt arrester information.
By to the water body region extracted above and dam body position range overlay analysis, obtaining final silt arrester information.
From technique scheme as can be seen that present embodiments providing a kind of extracting method of retaining water and soil phase silt arrester, The method is used for the silt arrester being in the retaining water and soil phase is extracted, and specially remote sensing image is processed, according to water body Spectral signature, space characteristics and Extraction of Geometrical Features water outlet body region;Spatial relationship according to water body and dam body is to dam body position Scope is confined;Water body region and dam body position range are overlapped analyzing, obtain silt arrester information.Obtaining silt arrester After information, the quantity of silt arrester and distribution situation in predetermined census area can be obtained by simple statistical, and without people Work mode is counted and on-site inspection to building profile, thus saving substantial amounts of manpower and materials, can effectively improve Statistical efficiency.
Embodiment two
Fig. 2 is the flow chart of steps of the extracting method embodiment of another kind of retaining water and soil phase silt arrester of the present invention.
It is used for the silt arrester being in the retaining water and soil phase is carried with reference to the extracting method that shown in Fig. 2, the present embodiment provides Take, specifically include following steps:
S201: according to water spectral feature, space characteristics and Extraction of Geometrical Features water outlet body region.
By remote sensing images analysis, to be distributed in mountain area and the hills area of northern Loess Plateau more than silt arrester, Build in trench bottom;The water body of silt arrester and general water body have similar spectral signature, but the water surface is less, many strips Shape, and water body downstream is connected with dam body.Therefore, according to the water spectral feature of water body, space characteristics and geometric properties to pre- The remote sensing image determining census area is processed, and therefrom proposes water body region.
Specific to being using normalization water body index ndwi (normalized difference water in the present embodiment Index), digital elevation model dem and coverage of water extract this water body region.
Normalization water body index ndwi is the normalized ratio index based on green light band and near infrared band.This formula For:
Ndwi=(green-nir)/(green+nir)
In formula, green green light band;Nir near infrared band.
The water body extracted region effect being obtained by said method is as shown in Figure 2 a.
S202: dam body position range position range is confined according to the spatial relationship in water body region and dam body.
Water body region comprises water body upstream termination and water body downstream end, and wherein water body downstream end refers to the dam of silt arrester , near the part of this dam body, water body upstream termination is corresponding with water body downstream end for water body more than body, and away from dam body according to This principle is confined to the dam body position range residing for dam body..Process of specifically confining is:
Step 1: create upstream region of interest and downstream region of interest.
In the achievement of Clean water withdraw, based on closing operation of mathematical morphology principle, first expand and shrink again, extract a pixel wide Region of interest, respectively as the extraction basis of water body upstream termination and water body downstream end.Constructed region of interest reference picture 2b Shown.
Step 2: water body upstream termination is extracted from upstream region of interest according to default upstream termination extraction algorithm.
Water body typically becomes ribbon to be distributed, and usual upstream breadth is less than downstream, and water body centroid position is partial to downstream position, Therefore in region of interest pixel to water body barycenter distance be not less than water body length half.Meanwhile, in order to avoid in computed range When, other impacts that closely water body is likely to result in upstream Clean water withdraw, introduce inheritance here, by water body figure spot and upper Trip region of interest is as the parent object of water body.Then, in the region of interest of upstream, the centroid distance of pixel to water body parent object is not less than accordingly The half of water body length.
Here the upstream termination extraction algorithm upstream building is:
Upstream=[distance to superobject] -1/2 × [length (water)]
In formula, upstream is water body upstream termination;Distance to superobject is apart from water body parent object matter The distance of the heart;Length (water) is the length of corresponding water body figure spot.Extraction water body upstream termination can be reached by this formula Purpose, wherein water body upstream termination extract before effect reference picture 2c shown in, the effect after extraction is as shown in Figure 2 d.
Correlation distance (distance to) principle therein be select one kind to classify (parasang is for atural object and distance Pixel), with the center of atural object figure spot as initial point, calculate classified atural object and the unfiled atural object figure spot central point in all directions away from From figure spot in selected scope for the distance can be divided into the classification of formulation.As follows apart from the computing formula of d:
d = σ f ( v f ( s ) - v f ( o ) σ f ) 2
In formula:
D is the distance between sample figure spot s and target figure spot o;
vfS () is the eigenvalue (such as figure spot central point pixel coordinate value) of sample figure spot f feature;
vfO () is the eigenvalue (such as figure spot central point pixel coordinate value) of target figure spot f feature;
σfStandard deviation for f eigenvalue.
The computing formula of length length therein is by the length splitting the approximate acquisition of primitive object bounding box and wide ratio It is calculated, formula is as follows:
l = # p v · γ v
In formula:
The atural object figure spot object that v is formed through over-segmentation;
#pvThe pixel total number of atural object figure spot v;
γvImage object v length and wide ratio;
The object length that l application v length-width ratio calculates.
Step 3: water body downstream end is extracted from region of interest according to default downstream end extraction algorithm.
Distance according to water body upstream termination to downstream end is approximate with this water body length, builds water body downstream end and extracts Algorithm downstream, determines the position of water body downstream end.Wherein, the effect such as Fig. 2 e institute before water body downstream end extracts Show, the effect after extraction is as shown in figure 2f.
Downstream=[distance to (upstream)]/[length (water)]
In formula, downstream is downstream end water body;Distance to (upstream) is apart from upstream termination water body The distance of figure spot;Length is the length of water body figure spot.
Step 4: the spatial relationship being connected with dam body according to water body water body downstream end, using morphological dilations principle, Water body upstream termination creates relief area, confines the dam body position range of silt arrester.The effect confined is as shown in Figure 2 g.
S203: water body region and dam body position range are overlapped analyzing, obtain silt arrester information.
By to the water body region extracted above and dam body position range overlay analysis, obtaining final silt arrester information, This silt arrester information is as shown in fig. 2h.
S204: silt arrester information is optimized perfect.
In order to obtain more accurate silt arrester information, it is possible to use manually the silt arrester information tentatively obtaining is looked into See, and according to visual and experience input data revision directive, this data correction instructs for above-mentioned preliminary silt arrester information It is optimized perfect, thus obtaining the higher silt arrester information of precision.
S205: the quantity and distribution of silt arrester is counted.
In the case of necessary, can according to above-mentioned silt arrester information to the quantity of the silt arrester in predetermined census area and Distribution is counted such that it is able to provide useful decision information for government department.
From technique scheme as can be seen that present embodiments providing a kind of extracting method of retaining water and soil phase silt arrester, The method is used for the silt arrester being in the retaining water and soil phase is extracted, and specially remote sensing image is processed, according to water body Spectral signature, space characteristics and Extraction of Geometrical Features water outlet body region;Spatial relationship according to water body and dam body is to dam body position Scope is confined;Water body region and dam body position range are overlapped analyzing, obtain silt arrester information.Obtaining silt arrester After information, the quantity of silt arrester and distribution situation in predetermined census area can be obtained by simple statistical, and without people Work mode is counted and on-site inspection to building profile, thus saving substantial amounts of manpower and materials, can effectively improve Statistical efficiency.
Wherein, send out when a scape gf-1pms remote sensing image of 2014 Nian8Yue Fen Shaanxi Province Yulin City Hengshan County being tested Existing, technical scheme that the present embodiment provides and extracts effect to the extraction accuracy rate of retaining water and soil phase silt arrester up to more than 90% Rate is higher, up to 21min/1000km2.Shown in concrete extraction effect reference picture 2i.
It should be noted that for embodiment of the method, in order to be briefly described, therefore it is all expressed as a series of action group Close, but those skilled in the art should know, the embodiment of the present invention is not limited by described sequence of movement, because according to According to the embodiment of the present invention, some steps can be carried out using other orders or simultaneously.Secondly, those skilled in the art also should Know, embodiment described in this description belongs to preferred embodiment, the involved action not necessarily present invention is implemented Example is necessary.
Embodiment three
Fig. 3 is a kind of structured flowchart of the extraction system embodiment of retaining water and soil phase silt arrester of the present invention.
It is used for the silt arrester being in the retaining water and soil phase is carried with reference to the extraction system that shown in Fig. 3, the present embodiment provides Take, specifically include Clean water withdraw module 10, dam body confines module 20 and information extraction modules 30.
Clean water withdraw module 10 is used for according to water spectral feature, space characteristics and Extraction of Geometrical Features water outlet body region.
By remote sensing images analysis, to be distributed in mountain area and the hills area of northern Loess Plateau more than silt arrester, Build in trench bottom;The water body of silt arrester and general water body have similar spectral signature, but the water surface is less, many strips Shape, and water body downstream is connected with dam body.Therefore, according to the water spectral feature of water body, space characteristics and geometric properties to pre- The remote sensing image determining census area is processed, and therefrom proposes water body region.
Specific to the present embodiment, this Clean water withdraw module 10 includes water body area extracting unit 11, and this unit is used for adopting Normalization water body index ndwi (normalized difference water index), digital elevation model dem and water body face Long-pending this water body region of extraction.
Normalization water body index ndwi is the normalized ratio index based on green light band and near infrared band.This formula For:
Ndwi=(green-nir)/(green+nir)
In formula, green green light band;Nir near infrared band.
Dam body confines unit 20 for confining to dam body position range according to the spatial relationship in water body region and dam body.
Water body region comprises water body upstream termination and water body downstream end, and wherein water body downstream end refers to the dam of silt arrester , near the part of this dam body, water body upstream termination is corresponding with water body downstream end for water body more than body, and away from dam body according to This principle is confined to the dam body position range residing for dam body.
Information extraction modules 30 are used for being overlapped analyzing by water body region and dam body position range, obtain silt arrester letter Breath.
This module is by the water body region extracted above and dam body position range overlay analysis, obtaining final silt arrester Information.
From technique scheme as can be seen that present embodiments providing a kind of extraction system of retaining water and soil phase silt arrester, This system is used for the silt arrester being in the retaining water and soil phase is extracted, and specially remote sensing image is processed, according to water body Spectral signature, space characteristics and Extraction of Geometrical Features water outlet body region;Spatial relationship according to water body and dam body is to dam body position Scope is confined;Water body region and dam body position range are overlapped analyzing, obtain silt arrester information.Obtaining silt arrester After information, the quantity of silt arrester and distribution situation in predetermined census area can be obtained by simple statistical, and without people Work mode is counted and on-site inspection to building profile, thus saving substantial amounts of manpower and materials, can effectively improve Statistical efficiency.
Example IV
Fig. 4 is the structured flowchart of the extraction system embodiment of another kind of retaining water and soil phase silt arrester of the present invention.
It is used for the silt arrester being in the retaining water and soil phase is carried with reference to the extraction system that shown in Fig. 4, the present embodiment provides Take, specifically include Clean water withdraw module 10, dam body confines module 20, information extraction modules 30, optimize module 40 and information Statistical module 50.
Clean water withdraw module 10 is used for according to water spectral feature, space characteristics and Extraction of Geometrical Features water outlet body region.
By remote sensing images analysis, to be distributed in mountain area and the hills area of northern Loess Plateau more than silt arrester, Build in trench bottom;The water body of silt arrester and general water body have similar spectral signature, but the water surface is less, many strips Shape, and water body downstream is connected with dam body.Therefore, according to the water spectral feature of water body, space characteristics and geometric properties to pre- The remote sensing image determining census area is processed, and therefrom proposes water body region.
Specific in the present embodiment, this module includes water body area extracting unit 11, and water body area extracting unit 11 is used for Using normalization water body index ndwi (normalized difference water index), digital elevation model dem and water Bulk area extracts this water body region.
Normalization water body index ndwi is the normalized ratio index based on green light band and near infrared band.This formula For:
Ndwi=(green-nir)/(green+nir)
In formula, green green light band;Nir near infrared band.
Dam body confines module 20 for confining to dam body position range according to the spatial relationship in water body region and dam body.
Water body region comprises water body upstream termination and water body downstream end, and wherein water body downstream end refers to the dam of silt arrester , near the part of this dam body, water body upstream termination is corresponding with water body downstream end for water body more than body, and away from dam body according to This principle is confined to the dam body position range residing for dam body.Dam body is confined module 20 and is specifically included region of interest creating unit 21st, end extraction unit 22, water body downstream end extraction unit 23 and dam body position range is had to confine unit 24 on water body.
Region of interest creating unit 21 is used for creating region of interest.
In the achievement of Clean water withdraw, based on closing operation of mathematical morphology principle, first expand and shrink again, extract a pixel wide Region of interest, respectively as the extraction basis of water body upstream termination and water body downstream end.
Water body upstream termination extraction unit 22 is used for extracting water body according to default upstream termination extraction algorithm from region of interest Upstream termination.
Water body typically becomes ribbon to be distributed, and usual upstream breadth is less than downstream, and water body centroid position is partial to downstream position, Therefore in region of interest pixel to water body barycenter distance be not less than water body length half.Meanwhile, in order to avoid in computed range When, other impacts that closely water body is likely to result in upstream Clean water withdraw, introduce inheritance here, by water body figure spot and upper Trip region of interest is as the parent object of water body.Then, in region of interest, the centroid distance of pixel to water body parent object is not less than corresponding water body The half of length.
Here the upstream termination extraction algorithm upstream building is:
Upstream=[distance to superobject] -1/2 × [length (water)]
In formula, upstream is water body upstream termination;Distance to superobject is apart from water body parent object matter The distance of the heart;Length (water) is the length of corresponding water body figure spot.Extraction water body upstream termination can be reached by this formula Purpose.
Correlation distance (distance to) principle therein be select one kind to classify (parasang is for atural object and distance Pixel), with the center of atural object figure spot as initial point, calculate classified atural object and the unfiled atural object figure spot central point in all directions away from From figure spot in selected scope for the distance can be divided into the classification of formulation.As follows apart from the computing formula of d:
d = σ f ( v f ( s ) - v f ( o ) σ f ) 2
In formula:
D is the distance between sample figure spot s and target figure spot o;
vfS () is the eigenvalue (such as figure spot central point pixel coordinate value) of sample figure spot f feature;
vfO () is the eigenvalue (such as figure spot central point pixel coordinate value) of target figure spot f feature;
σfStandard deviation for f eigenvalue.
The computing formula of length length therein is by the length splitting the approximate acquisition of primitive object bounding box and wide ratio It is calculated, formula is as follows:
l = # p v · γ v
In formula:
The atural object figure spot object that v is formed through over-segmentation;
#pvThe pixel total number of atural object figure spot v;
γvImage object v length and wide ratio;
The object length that l application v length-width ratio calculates.
Water body downstream end extraction unit 23 is used for extracting water body according to default downstream end extraction algorithm from region of interest Downstream end.
Distance according to water body upstream termination to downstream end is approximate with this water body length, builds water body downstream end and extracts Algorithm downstream, determines the position of water body downstream end.
Downstream=[distance to (upstream)]/[length (water)]
In formula, downstream is downstream end water body;Distance to (upstream) is apart from upstream termination water body The distance of figure spot;Length is the length of water body figure spot.
Dam body position range confines unit 24 for the sky being connected with dam body according to water body upstream termination, water body downstream end Between relation, using morphological dilations principle, create relief area in water body upstream termination, confine the dam body position range of silt arrester.
Information extraction modules 30 are used for
By to the water body region extracted above and dam body position range overlay analysis, obtaining final silt arrester information.
Optimize module 40 perfect for being optimized to the silt arrester information that information extraction modules 30 obtain.
In order to obtain more accurate silt arrester information, it is possible to use manually the silt arrester information tentatively obtaining is looked into See, and according to visual and experience input data revision directive, this data correction instructs for above-mentioned preliminary silt arrester information It is optimized perfect, thus obtaining the higher silt arrester information of precision.
Information Statistics module 50 is used for the quantity of silt arrester and distribution are counted.
In the case of necessary, can according to above-mentioned silt arrester information to the quantity of the silt arrester in predetermined census area and Distribution is counted such that it is able to provide useful decision information for government department.
From technique scheme as can be seen that present embodiments providing a kind of extraction system of retaining water and soil phase silt arrester, This system is used for the silt arrester being in the retaining water and soil phase is extracted, and specially remote sensing image is processed, according to water body Spectral signature, space characteristics and Extraction of Geometrical Features water outlet body region;Spatial relationship according to water body and dam body is to dam body position Scope is confined;Water body region and dam body position range are overlapped analyzing, obtain silt arrester information.Obtaining silt arrester After information, the quantity of silt arrester and distribution situation in predetermined census area can be obtained by simple statistical, and without people Work mode is counted and on-site inspection to building profile, thus saving substantial amounts of manpower and materials, can effectively improve Statistical efficiency.
For device embodiment, due to itself and embodiment of the method basic simlarity, so description is fairly simple, related Part illustrates referring to the part of embodiment of the method.
Each embodiment in this specification is all described by the way of going forward one by one, what each embodiment stressed be with The difference of other embodiment, between each embodiment identical similar partly mutually referring to.
Those skilled in the art are it should be appreciated that the embodiment of the embodiment of the present invention can be provided as method, device or calculate Machine program product.Therefore, the embodiment of the present invention can using complete hardware embodiment, complete software embodiment or combine software and The form of the embodiment of hardware aspect.And, the embodiment of the present invention can using one or more wherein include computer can With in the computer-usable storage medium (including but not limited to disk memory, cd-rom, optical memory etc.) of program code The form of the computer program implemented.
The embodiment of the present invention is with reference to method according to embodiments of the present invention, terminal unit (system) and computer program The flow chart of product and/or block diagram are describing.It should be understood that can be by computer program instructions flowchart and/or block diagram In each flow process and/or the flow process in square frame and flow chart and/or block diagram and/or square frame combination.These can be provided Computer program instructions set to general purpose computer, special-purpose computer, Embedded Processor or other programmable data processing terminals Standby processor is to produce a machine so that holding by the processor of computer or other programmable data processing terminal equipments The instruction of row produces for realizing in one flow process of flow chart or multiple flow process and/or one square frame of block diagram or multiple square frame The device of the function of specifying.
These computer program instructions may be alternatively stored in and can guide computer or other programmable data processing terminal equipments So that the instruction being stored in this computer-readable memory produces bag in the computer-readable memory working in a specific way Include the manufacture of command device, this command device is realized in one flow process of flow chart or multiple flow process and/or one side of block diagram The function of specifying in frame or multiple square frame.
These computer program instructions also can be loaded in computer or other programmable data processing terminal equipments so that Series of operation steps is executed on computer or other programmable terminal equipments to produce computer implemented process, thus On computer or other programmable terminal equipments, the instruction of execution is provided for realizing in one flow process of flow chart or multiple flow process And/or the step of the function of specifying in one square frame of block diagram or multiple square frame.
Although having been described for the preferred embodiment of the embodiment of the present invention, those skilled in the art once know base This creative concept, then can make other change and modification to these embodiments.So, claims are intended to be construed to Including preferred embodiment and fall into being had altered and changing of range of embodiment of the invention.
Last in addition it is also necessary to explanation, herein, such as first and second or the like relational terms be used merely to by One entity or operation are made a distinction with another entity or operation, and not necessarily require or imply these entities or operation Between there is any this actual relation or order.And, term " inclusion ", "comprising" or its any other variant meaning Covering comprising of nonexcludability, so that include a series of process of key elements, method, article or terminal unit not only wrapping Include those key elements, but also include other key elements being not expressly set out, or also include for this process, method, article Or the intrinsic key element of terminal unit.In the absence of more restrictions, by wanting that sentence "including a ..." limits It is not excluded that also there is other identical element in process, method, article or the terminal unit including described key element in element.
Above technical scheme provided by the present invention is described in detail, specific case used herein is to this Bright principle and embodiment are set forth, the explanation of above example be only intended to help understand the method for the present invention and its Core concept;Simultaneously for one of ordinary skill in the art, according to the thought of the present invention, in specific embodiment and application All will change in scope, in sum, this specification content should not be construed as limitation of the present invention.

Claims (10)

1. a kind of extracting method of retaining water and soil phase silt arrester is it is characterised in that specifically include step:
Remote sensing image is processed, according to water spectral feature, space characteristics and Extraction of Geometrical Features water outlet body region;
According to the spatial relationship of water body and dam body, dam body position range is confined;
Described water body region and described dam body position range are overlapped analyzing, obtain silt arrester information.
2. extracting method as claimed in claim 1 it is characterised in that described according to water spectral feature, space characteristics and several What feature extraction water outlet body region, comprising:
Described water body region is extracted according to normalization water body index, digital elevation model and coverage of water.
3. extracting method as claimed in claim 1 is it is characterised in that the described spatial relationship pair according to water body region and dam body Dam body position range is confined, comprising:
Region of interest is created based on closing operation of mathematical morphology principle;
Water-outlet body upstream termination is extracted in described region of interest according to preset upstream termination extraction algorithm;
Water-outlet body downstream end is extracted in described region of interest according to preset downstream end extraction algorithm;
Described dam body position model is confined out according to the spatial relationship of described water body upstream termination, described water body downstream end and dam body Enclose.
4. the extracting method as described in any one of claims 1 to 3 is it is characterised in that further comprise the steps of:
According to the data correction instruction of user input, described silt arrester information is optimized perfect.
5. extracting method as claimed in claim 4 is it is characterised in that further comprise the steps of:
According to described silt arrester information, the quantity of the silt arrester in predetermined census area and distribution are counted.
6. a kind of extraction system of retaining water and soil phase silt arrester is it is characterised in that specifically include:
Clean water withdraw module, for processing to remote sensing image, carries according to water spectral feature, space characteristics and geometric properties Take out water body region;
Dam body confines module, for being confined to dam body position range according to the spatial relationship of water body and dam body;
Information extraction modules, for being overlapped analyzing described water body region and described dam body position range, obtain silt arrester Information.
7. extraction system as claimed in claim 6 is it is characterised in that described Clean water withdraw module includes:
Water body area extracting unit, for extracting described water according to normalization water body index, digital elevation model and coverage of water Body region.
8. extraction system as claimed in claim 6 is it is characterised in that described dam body confines module includes:
Region of interest creating unit, for creating region of interest based on closing operation of mathematical morphology principle;
Water body upstream termination extraction unit, for extracting water outlet according to preset upstream termination extraction algorithm in described region of interest Body upstream termination;
Water body downstream end extraction unit, for extracting water outlet according to preset downstream end extraction algorithm in described region of interest Body downstream end;
Dam body position range confines unit, closes for the space according to described water body upstream termination, water body downstream end and dam body System confines out described dam body position range.
9. the extraction system as described in any one of claim 6~8 is it is characterised in that also include:
Optimize module, perfect for being optimized to described silt arrester information according to the data correction instruction of user input.
10. extraction system as claimed in claim 9 is it is characterised in that also include:
Information Statistics module, for according to described silt arrester information to the quantity of the silt arrester in predetermined census area and be distributed into Row statistics.
CN201610787008.6A 2016-08-30 2016-08-30 A kind of extracting method and system of retaining water and soil phase silt arrester Active CN106372599B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610787008.6A CN106372599B (en) 2016-08-30 2016-08-30 A kind of extracting method and system of retaining water and soil phase silt arrester

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610787008.6A CN106372599B (en) 2016-08-30 2016-08-30 A kind of extracting method and system of retaining water and soil phase silt arrester

Publications (2)

Publication Number Publication Date
CN106372599A true CN106372599A (en) 2017-02-01
CN106372599B CN106372599B (en) 2019-08-30

Family

ID=57899491

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610787008.6A Active CN106372599B (en) 2016-08-30 2016-08-30 A kind of extracting method and system of retaining water and soil phase silt arrester

Country Status (1)

Country Link
CN (1) CN106372599B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109597868A (en) * 2018-11-29 2019-04-09 中国水利水电科学研究院 A kind of a wide range of silt arrester location determining method
CN114529814A (en) * 2022-01-26 2022-05-24 南京信息工程大学 Loess plateau silt dam extraction method based on multi-source data
CN117809172A (en) * 2023-11-20 2024-04-02 北京师范大学 Dam identification method and device based on remote sensing picture

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103034863A (en) * 2012-12-24 2013-04-10 重庆市勘测院 Remote-sensing image road acquisition method combined with kernel Fisher and multi-scale extraction
CN103440655A (en) * 2013-08-27 2013-12-11 西北工业大学 Crossing bridge and offshore ship joint detection method in onboard remote sensing image
CN104912041A (en) * 2015-06-17 2015-09-16 中国电建集团贵阳勘测设计研究院有限公司 Fish passing system for high dam of canyon river reach
CN105654488A (en) * 2015-12-31 2016-06-08 西北农林科技大学 Method for estimating silt interception amount of silt dam based on GIS

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103034863A (en) * 2012-12-24 2013-04-10 重庆市勘测院 Remote-sensing image road acquisition method combined with kernel Fisher and multi-scale extraction
CN103440655A (en) * 2013-08-27 2013-12-11 西北工业大学 Crossing bridge and offshore ship joint detection method in onboard remote sensing image
CN104912041A (en) * 2015-06-17 2015-09-16 中国电建集团贵阳勘测设计研究院有限公司 Fish passing system for high dam of canyon river reach
CN105654488A (en) * 2015-12-31 2016-06-08 西北农林科技大学 Method for estimating silt interception amount of silt dam based on GIS

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109597868A (en) * 2018-11-29 2019-04-09 中国水利水电科学研究院 A kind of a wide range of silt arrester location determining method
CN109597868B (en) * 2018-11-29 2022-05-10 中国水利水电科学研究院 Large-range check dam position determination method
CN114529814A (en) * 2022-01-26 2022-05-24 南京信息工程大学 Loess plateau silt dam extraction method based on multi-source data
CN117809172A (en) * 2023-11-20 2024-04-02 北京师范大学 Dam identification method and device based on remote sensing picture

Also Published As

Publication number Publication date
CN106372599B (en) 2019-08-30

Similar Documents

Publication Publication Date Title
Hayakawa et al. GIS analysis of fluvial knickzone distribution in Japanese mountain watersheds
CN103940974A (en) Spatio-temporal dynamic analysis method of soil erosion in meso-scale watershed based on GIS
Barrett et al. Lakes and rivers in the landscape: a quasi-variational inequality approach
Zhang et al. Effects of initial step height on the headcut erosion of bank gullies: A case study using a 3D photo-reconstruction method in the dry-hot valley region of southwest China
CN106372599A (en) Method and system for extracting silt arrester in water and soil retaining period
Ghunowa et al. Stream power index for networks (SPIN) toolbox for decision support in urbanizing watersheds
Muluneh et al. Morphometric analysis of Didessa River catchment in blue nile basin, western Ethiopia
Pisaturo et al. A procedure for human safety assessment during hydropeaking events
Couturier et al. Improvement of surface flow network prediction for the modeling of erosion processes in agricultural landscapes
CN104750936A (en) GIS (Geographic Information System)-based urban hydrological run-off calculating method and system
CN109190593B (en) Mountain road side slope stability preliminary discrimination method along line based on concave-convex classification
CN106447510A (en) Terrace field information extracting method and terrace field information extracting system
Elagca Application of Arc-GIS, HEC-GeoHMS and HEC-HMS in a holistic sense for estimation of rainfall-runoff process: case study over Ballikaya Basin.
Alexakis et al. Assessing soil erosion rate in a catchment area in Cyprus using remote sensing and GIS
Malekani et al. Application of GIS in modeling Zilberchai Basin runoff
Overton et al. Flood modelling and vegetation mapping in large river systems
Rigon et al. Application of the new Morphological Quality Index in the Cordevole river (BL, Italy)
Senadeera Morphmetric analysis of Randenigala reservoir catchment using GIS
Wu et al. Spatial variation of channel head curvature in small mountainous watersheds
Šamanović et al. Influence of pit removal algorithms on surface runoff simulation
Liliang et al. Application of radar-measured rain data in hydrological processes modeling during the intensified observation period of HUBEX
CN109724570A (en) The drop amount of underground drop, drop width, on bank water layer thickness calculation method
Abdelsalheen et al. Evaluation of Depression Storage Using Grid-Based GIS Model
Moretto et al. Integration of colour bathymetry, LiDAR and dGPS surveys for assessing fluvial changes after flood events in the Tagliamento River (Italy)
Kang Geospatial analysis of specific degradation in South Korea

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant