REMOTE SENSING-BASED EXTRACTION METHOD FOR TYPE OF RIVER CHANNEL

Abstract

The present disclosure discloses a remote sensing-based extraction method for a type of a river channel, including: obtaining multi-source remote sensing data of a target area including a river channel; preprocessing the multi-source remote sensing data, and obtaining corresponding reflectance data; according to the reflectance data, analyzing a water index and a vegetation index of the target area; and according to the water index and the vegetation index, constructing first preset conditions for determining whether there is water in the river channel and second preset conditions for determining whether the river channel is a non-dry river channel to determine the type of the river channel. Types of river channels can be divided into four types: non-dry river channels, seasonal dry river channels, temporary water channels, and dry river channels.

Claims

1. A remote sensing-based extraction method for a type of a river channel, comprising: obtaining multi-source remote sensing data of a target area comprising a river channel; preprocessing the multi-source remote sensing data, and obtaining corresponding reflectance data; according to the reflectance data, calculating a water index and a vegetation index of the target area; and according to the water index and the vegetation index, determining the type of the river channel.

2. The remote sensing-based extraction method for a type of a river channel according to claim 1, wherein the preprocessing the multi-source remote sensing data comprises: performing geometric, radiometric, and atmospheric correction processing on the multi-source remote sensing data.

3. The remote sensing-based extraction method for a type of a river channel according to claim 2, further comprising: if it is determined that there is water in the river channel, vectorizing raster data corresponding to the remote sensing data of the river channel to form a series of water patches; and calculating a length of the connected patches as a length of water existing in the river water within a certain time range.

4. The remote sensing-based extraction method for a type of a river channel according to claim 1, wherein the calculating a water index and a vegetation index of the target area comprises: if the multi-source remote sensing data is a remote sensing image with more than four bands, obtaining a normalized difference vegetation index (NDVI) and an modified normalized difference water index (MNDWI); and if the multi-source remote sensing data is a remote sensing image with four bands, obtaining an NDVI and a normalized difference water index (NDWI)-NDVI.

5. The remote sensing-based extraction method for a type of a river channel according to claim 4, further comprising: if it is determined that there is water in the river channel, vectorizing raster data corresponding to the remote sensing data of the river channel to form a series of water patches; and calculating a length of the connected patches as a length of water existing in the river water within a certain time range.

6. The remote sensing-based extraction method for a type of a river channel according to claim 4, wherein the MNDWI is expressed as: $MNDWI=\frac{R_g-R_{swir}}{R_{swir}+R_g}$ wherein in the formula (1), MNDWI is the modified normalized difference water index; and R.sub.swir is reflectance in a short-wave infrared (SWIR) band, and R.sub.9 is reflectance in a red light band; and the NDWI-NDVI is expressed as: $NDWI-NDVI=\frac{R_g-R_{nir}}{R_g+R_{nir}}-NDVI$ wherein in the formula (2), NDVI is the normalized difference vegetation index, NDWI is the normalized difference water index, and R.sub.nir is reflectance in a near-infrared (NIR) band.

7. The remote sensing-based extraction method for a type of a river channel according to claim 6, further comprising: if it is determined that there is water in the river channel, vectorizing raster data corresponding to the remote sensing data of the river channel to form a series of water patches; and calculating a length of the connected patches as a length of water existing in the river water within a certain time range.

8. The remote sensing-based extraction method for a type of a river channel according to claim 6, wherein the NDVI is expressed as: $NDVI=\frac{R_{nir}-R_r}{R_{nir}+R_r}$ wherein in the formula (3), NDVI is the normalized difference vegetation index, and R.sub.nir is the reflectance in the NIR band; and R.sub.r is reflectance in a red band.

9. The remote sensing-based extraction method for a type of a river channel according to claim 8, further comprising: if it is determined that there is water in the river channel, vectorizing raster data corresponding to the remote sensing data of the river channel to form a series of water patches; and calculating a length of the connected patches as a length of water existing in the river water within a certain time range.

10. The remote sensing-based extraction method for a type of a river channel according to claim 8, wherein the determining the type of the river channel comprises: according to the water index and the vegetation index, constructing first preset conditions for determining whether there is water in the river channel and second preset conditions for determining whether the river channel is a non-dry river channel to determine the type of the river channel.

11. The remote sensing-based extraction method for a type of a river channel according to claim 10, further comprising: if it is determined that there is water in the river channel, vectorizing raster data corresponding to the remote sensing data of the river channel to form a series of water patches; and calculating a length of the connected patches as a length of water existing in the river water within a certain time range.

12. The remote sensing-based extraction method for a type of a river channel according to claim 10, wherein the first preset conditions for determining whether there is water in the river channel comprise: 1) during identification of a pixel of a non-ice and snow-covered river channel, if MNDWI.sub.t0 ≥ a.sub.1 is satisfied, or MNDWI.sub.ti < a.sub.1 , NDVI.sub.ti ≥ a.sub.2 , MNDWI.sub.ti-tj ≥ a.sub.1 , and NDVI.sub.ti-tj < a.sub.2 are satisfied, determining that the river channel is a river channel with water, wherein MNDWI.sub.t0 , MNDWI.sub.ti , MNDWI.sub.ti-tj , NDVI.sub.ti , and NDVI.sub.ti-tj are corresponding MNDWIs and NDVIs at times t0, ti, and ti-tj respectively; and a.sub.1 and a.sub.2 are a first preset threshold and a second preset threshold respectively; and 2) during identification of a pixel of an ice and snow-covered river channel, constructing an ice and snow coverage model index and an ice and snow reflection intensity model $IS=\frac{R_g-R_{nir}}{R_{nir}+R_g}$ .sub.g index IT = R.sub.g.sup.2 + R.sub.r.sup.2 + R.sub.nir.sup.2, and if IS ≥ a.sub.3 and IT ≥ a.sub.4, determining that the ice and snow-covered river channel is a river channel with water, wherein a.sub.3 and a.sub.4 are a third preset threshold and a fourth preset threshold respectively.

13. The remote sensing-based extraction method for a type of a river channel according to claim 12, further comprising: if it is determined that there is water in the river channel, vectorizing raster data corresponding to the remote sensing data of the river channel to form a series of water patches; and calculating a length of the connected patches as a length of water existing in the river water within a certain time range.

14. The remote sensing-based extraction method for a type of a river channel according to claim 12, wherein the second preset conditions for determining whether the river channel is a non-dry river channel comprise: 3) establishing a non-dry river channel identification index for a non-ice and snow-covered season time series: $NISI=\sqrt{{\displaystyle \underset{i=1}{\overset{n_1}{.Math.}}\frac{{\left(ND{W}_i-\overline{NDW}\right)}^2}{n_1}}},\text{if}\overline{NDW}\ge {\text{a}}_5$ is satisfied, determining that the river channel of the non-ice and snow-covered season time series is a non-dry river channel, wherein in the formula, .sup.NISI is a non-dry pixel index; NDW is a non-negative water value, NDW = MNDWI + b.sub.1, wherein b.sub.1 is a preset threshold; NDW is an average value in a statistical time range; a.sub.5 is a fifth preset threshold; n.sub.1 is a number of statistical image scenes of the non-ice and snow-covered season time series; and i is an index value of the number of the statistical image scenes of the non-ice and snow-covered season time series; and 4) establishing a non-dry river channel identification model for an ice and snow-covered season time series: $ISI={\displaystyle \underset{i=1}{\overset{n_2}{.Math.}}\frac{{\left(I{S}_i-\overline{IS}\right)}^2}{n_2},}$ if $\overline{IS}\ge a_6$ is satisfied, and $ITI={\displaystyle \underset{i=1}{\overset{n_2}{.Math.}}\frac{{\left(I{T}_i-\overline{IT}\right)}^2}{n_2},}$ if $\overline{IT}\ge a_7$ is satisfied, determining that the river channel of the ice and snow-covered season time series is a non-dry river channel, wherein in the formula, ISIand ITIare an ice and snow coverage index and an ice and snow reflection intensity index respectively; IS and IT are an average ice and snow coverage model index and an average ice and snow reflection intensity model index in a statistical time range respectively; a.sub.6 is a sixth preset threshold, and a.sub.7 is a seventh preset threshold; and n.sub.2 is a number of statistical image scenes of the ice and snow-covered season time series.

15. The remote sensing-based extraction method for a type of a river channel according to claim 14, further comprising: if it is determined that there is water in the river channel, vectorizing raster data corresponding to the remote sensing data of the river channel to form a series of water patches; and calculating a length of the connected patches as a length of water existing in the river water within a certain time range.

16. The remote sensing-based extraction method for a type of a river channel according to claim 14, wherein the process of determining the type of the river channel comprises: A, if 3) and 4) in the second preset conditions for determining whether the river channel is a non-dry river channel are both satisfied, determining that the river channel is a non-dry river channel; B, if only one of 3) and 4) in the second preset conditions for determining whether the river channel is a non-dry river channel is satisfied, determining that the river channel is a seasonal dry river channel; or if 3) and 4) in the second preset conditions for determining whether the river channel is a non-dry river channel are not satisfied, but in a t1 to t2 time range of a precipitation season, the non-dry river channel identification index for the non-ice and snow-covered season time series is as follows: $NISI=\sqrt{{\displaystyle \underset{i=1}{\overset{n_3}{.Math.}}\frac{{\left(ND{W}_i-\overline{NDW}\right)}^2}{n_3},\text{if}\overline{NDW}\ge a_8}}$ is satisfied, determining that the river channel is a seasonal dry river channel of the non-ice and snow-covered season time series, wherein NISI is the non-dry river channel identification index; and n.sub.3 is a number of statistical pixels of the non-ice and snow-covered season time series in the t1 to t2 time range of the precipitation season; and C, during remote sensing monitoring of river water in a time series, if 1) in the first preset conditions for determining whether there is water in the river channel is satisfied within a certain time range, and a total time interval for monitoring water t is less than Δt1, and is not within the t1 to t2 time range of the precipitation season, determining that the river channel is a temporary water channel, wherein Δt1 is a preset time threshold.

17. The remote sensing-based extraction method for a type of a river channel according to claim 16, further comprising: if it is determined that there is water in the river channel, vectorizing raster data corresponding to the remote sensing data of the river channel to form a series of water patches; and calculating a length of the connected patches as a length of water existing in the river water within a certain time range.

18. The remote sensing-based extraction method for a type of a river channel according to claim 1, further comprising: if it is determined that there is water in the river channel, vectorizing raster data corresponding to the remote sensing data of the river channel to form a series of water patches; and calculating a length of the connected patches as a length of water existing in the river water within a certain time range.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0045] FIG. 1 is a flow chart of a remote sensing-based extraction method for a type of a river channel provided by an embodiment of the present disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0046] In order to make the technical means, creative features, objectives to be achieved, and efficacy achieved by the disclosure easy to understand, the present disclosure will be further described below in combination with implementations.

[0047] In the description of the present disclosure, it should be noted that orientation or position relationships indicated by terms such as “upper”, “lower”, “inner”, “outer”, “front end”, “rear end”, “both ends”, “one end” and “the other end” are described based on the accompanying drawings. These terms are merely intended to facilitate and simplify the description of the present disclosure, rather than to indicate or imply that the mentioned device or components must have a specific orientation or must be constructed and operated in a specific orientation. Therefore, these terms should not be understood as a limitation to the present disclosure. Moreover, terms such as “first” and “second” are merely intended for the purpose of description, and should not be construed as indicating or implying relative importance.

[0048] In the description of the present disclosure, it should be noted that, unless otherwise clearly specified and limited, meanings of terms “installed”, “disposed with”, and “connected with” should be understood in a board sense. For example, the “connection” may be a fixed connection, a removable connection, or an integral connection; may be a mechanical connection or an electrical connection; may be a direct connection or an indirect connection through an intermediate medium; or may be intercommunication between two components. Those of ordinary skill in the art may understand specific meanings of the foregoing terms in the present disclosure based on a specific situation.

[0049] With reference to FIG. 1, a remote sensing-based extraction method for a type of a river channel provided by the present disclosure includes the following steps.

[0050] S10, Multi-source remote sensing data of a target area including a river channel is obtained.

[0051] S20, The multi-source remote sensing data is preprocessed, and corresponding reflectance data is obtained.

[0052] S30, According to the reflectance data, a water index and a vegetation index of the target area are calculated.

[0053] S40, According to the water index and the vegetation index, the type of the river channel is determined.

[0054] The present disclosure aims to provide a seasonal dry river channel for monitoring and management of the ecological environment of river channels on the basis of fully considering change characteristics of an underlying surface of a seasonal dry river channel and a river channel in a wet season and a fluctuating period and spectral characterization of multi-source remote sensing data, and provides a remote sensing prediction method for predicting a type of a river channel on this basis, which has high identification accuracy of the type of the river channel.

[0055] The specific solution of the present disclosure will be described in detail below.

[0056] 1. Data preparation:

[0057] In steps S10 to S20, multi-source remote sensing data is collected, which may also include water system data, land use, and digital elevation model (DEM) data. Geometric, radiometric, and atmospheric correction processing is performed on the remote sensing data to obtain reflectance data, such as reflectance in a SWIR band, reflectance in a red light band, and reflectance in a NIR band.

[0058] 2. Water index and vegetation index:

[0059] In step S30, there are two cases.

[0060] First, for a remote sensing image with more than four bands:

[0061] an MNDWI model or an NDWI-NDVI model can be selected.

[0062] The MNDWI model is expressed as:

[00009]$\text{MNDWI}\text{=}\frac{{\text{R}}_{\text{g}}-{\text{R}}_{\text{swir}}}{{\text{R}}_{\text{swir}}+{\text{R}}_{\text{g}}}$

[0063] In the formula (1), MNDWI is the modified normalized difference water index. R.sub.swir is reflectance in a SWIR band, and R.sub.g is reflectance in a red light band.

[0064] The NDWI-NDVI model is expressed as:

[00010]$NDWI-NDVI=\frac{R_g-R_{nir}}{R_g+R_{nir}}-NDVI$

[0065] In the formula (2), NDVI is the normalized difference vegetation index, NDWI is the normalized difference water index, and R.sub.nir is reflectance in a NIR band.

[0066] Second, for a remote sensing image with four bands:

[0067] an NDVI model is selected.

[00011]$NDVI=\frac{R_{nir}-R_r}{R_{nir}+R_r}$

[0068] In the formula (3), NDVI is the normalized difference vegetation index, and R.sub.nir is the reflectance in the NIR band. R.sub.r is reflectance in a red band.

[0069] 3. Determination of whether there is water in the river channel or water system:

[0070] Two first preset conditions for determining whether there is water in the river channel are constructed according to the water index and the vegetation index in step S40.

[0071] (1) During identification of a pixel of a non-ice and snow-covered river channel:

[0072] If MNDWI.sub.t0 ≥ a.sub.1 is satisfied, it is determined that the river channel is a river channel with water.

[0073] Or considering the impact of vegetation, if MNDWI.sub.ti < a.sub.1, NDVI.sub.ti ≥ a.sub.2, MNDWI.sub.ti-tj ≥ a.sub.1, and NDVI.sub.ti-tj < a.sub.2 are satisfied, it is determined that the river channel is a river channel with water when covered by vegetation.

[0074] MNDWI.sub.t0, MNDWI.sub.ti, MNDWI.sub.ti-tj, NDVI.sub.ti, and NDVI.sub.ti-tj are corresponding MNDWIs and NDVIs at times t0, ti, and ti-tj respectively. t0 is a time starting point of monitoring, ti and ti-tj are both any time points after the time starting point t0 of monitoring, and ti≠ti-tj.

[0075] a.sub.1 and a.sub.2 are a first preset threshold and a second preset threshold respectively. The values of a.sub.1 and a.sub.2 are related to the type of the remote sensing image and the numerical index of the remote sensing image respectively. The subsequent a.sub.3-a.sub.8 are related to the type of the remote sensing image and the numerical index of the remote sensing image. The determination of these thresholds is related to specific remote sensing images, which is not limited in the embodiment of the present disclosure.

[0076] A river channel that satisfies the above conditions is considered as a river channel with water. In order to avoid the occurrence of non-negative values, the following is established:

[00012]$NDW=MNDWI+b_{1\cdot }$

[0077] In the formula, NDW is a non-negative water value, and b.sub.1 is a preset threshold. The preset threshold is related to the type of remote sensing and communication resolution, and the appropriate value can be selected according to the specific situation.

[0078] (2) During identification of a pixel of an ice and snow-covered river channel:

[0079] An ice and snow coverage model index

[00013]$\text{IS}\text{=}\frac{{\text{R}}_{\text{g}}-{\text{R}}_{\text{nir}}}{{\text{R}}_{\text{nir}}\text{+}{\text{R}}_{\text{g}}}$

and an ice and snow reflection intensity model index IT = R.sub.g.sup.2 + R.sub.r.sup.2 + R.sub.nir.sup.2 are constructed. If IS ≥ a.sub.3 and IT ≥ a.sub.4 are both satisfied, it is determined that the pixel is covered by ice and snow, that is, it is determined that the ice and snow-covered river channel is a river channel with water. a.sub.3 and a.sub.4 are a third preset threshold and a fourth preset threshold respectively.

[0080] 4. Determination of identification of pixel of non-dry river water:

[0081] Two second preset conditions for determining whether the river channel is a non-dry river channel are constructed according to the water index and the vegetation index in step S40.

[0082] (1) A non-dry river channel identification index for a non-ice and snow-covered season time series is established: is satisfied, [0084] it is determined that the pixel is a non-dry pixel, that is, it is determined that the river channel of the non-ice and snow-covered season time series is a non-dry river channel.

[0085] In the formula, NISI is a non-dry pixel index. NDW is an average value in a statistical time range. a.sub.5 is a fifth preset threshold. n.sub.1 is a number of statistical image scenes of the non-ice and snow-covered season time series. i is an index value of the number of the statistical image scenes of the non-ice and snow-covered season time series.

[0086] (2) A non-dry river channel identification model for an ice and snow-covered season time series is established: is satisfied, and is satisfied, [0089] it is determined that the pixel is a non-dry pixel in the ice and snow-covered season, that is, it is determined that the river channel of the ice and snow-covered season time series is a non-dry river channel.

[0090] In the formula, ISI and ITI are an ice and snow coverage index and an ice and snow reflection intensity index respectively. IS and IT are an average ice and snow coverage model index and an average ice and snow reflection intensity model index in a statistical time range respectively. a.sub.6 is a sixth preset threshold, and a.sub.7 is a seventh preset threshold. n.sub.2 is a number of statistical image scenes of the ice and snow-covered season time series.

[0091] If the above 4(1) and 4(2) are both satisfied, it is a non-dry river channel.

[0092] 5. Identification of pixel of seasonal river water:

[0093] A process of determining the type of the river channel in step S40 includes the following sub-steps.

[0094] A, If only one of 4(1) and 4(2) is satisfied, it is a seasonal water pixel, that is, it is determined that the channel is a non-dry river channel.

[0095] B, If 4(2) and 4(1) are not satisfied, but in a t1 to t2 time range of a precipitation season, the non-dry river channel identification index for the non-ice and snow-covered season time series is as follows: is satisfied, [0097] it is determined that the pixel is a non-dry pixel, that is, it is determined that the river channel is a seasonal dry river channel of the non-ice and snow-covered season time series. n.sub.3 is a number of statistical pixels of the non-ice and snow-covered season time series in the t1 to t2 time range of the precipitation season.

[0098] 6. Identification of pixel of temporary river water:

[0099] During remote sensing monitoring of river water in a time series, if 3(1) in the first preset conditions for determining whether there is water in the river channel is satisfied within a certain time range, and a total time interval for monitoring water t is less than Δt1, and is not within the t1 to t2 time range of the precipitation season, it is determined that the pixel is a temporary river water pixel, that is, it is determined that the river channel is a temporary water channel. Δt1 is a preset time threshold. The preset time threshold is related to the frequency of acquisition and monitoring of remote sensing images, for example, it may be 14 days. It is not within the t1 to t2 range of the precipitation season. For example, in Beijing, it is not between July and August.

[0100] 7. Identification of pixel of dry river channel

[0101] If 3(1) or 3(2) in the first preset conditions for determining whether there is water in the river channel do not satisfy the conditions within the monitoring time range, it is determined that the pixel is a dry river pixel.

[0102] 8. Identification of length of river water

[0103] Further, the method includes:

[0104] S50, If it is determined that there is water in the river channel, raster data corresponding to the remote sensing data of the river channel is vectorized to form a series of water patches. A length of the connected patches is calculated as a length of water existing in the river water within a certain time range.

[0105] If it is determined that there is water in the river channel, the above obtained raster data of the river channel is vectorized to form a series of water patches. Through the method of spatial analysis, if two adjacent patches exist:

[00018]$\Delta l=\sqrt{{\left(p_{x1}-p_{x2}\right)}^2+\left(p_{y1}-p_{y2}\right)}\left({}^2\right)\text{and}\Delta l\le {\text{a}}_{8_{,}}$

it is considered that the two patches are connected.

[0106] Δl is a judgment example, a.sub.8 is an eighth preset threshold, (p.sub.x1, p.sub.y1) is the coordinate value of the first patch, and (p.sub.x2, p.sub.y2) is the coordinate value of the second patch adjacent to the first patch. The length of the connected patches is calculated as the length of water existing in the river water within a certain time range.

[0107] The remote sensing-based extraction method for a type of a river channel provided by the present disclosure includes: obtaining multi-source remote sensing data of a target area including a river channel; preprocessing the multi-source remote sensing data, and obtaining corresponding reflectance data; according to the reflectance data, analyzing a water index and a vegetation index of the target area; and according to the water index and the vegetation index, determining the type of the river channel. The present disclosure has advantages that the present disclosure aims to provide a seasonal dry river channel for monitoring and management of the ecological environment of river channels on the basis of fully considering change characteristics of an underlying surface of a seasonal dry river channel and a river channel in a wet season and a fluctuating period and spectral characterization of multi-source remote sensing data, and provides a remote sensing prediction method for predicting a type of a river channel on this basis, which has high identification accuracy.

[0108] The above shows and describes the basic principles, main features, and advantages of the present disclosure. It should be understood by those skilled in the art that the present disclosure is not limited by the above embodiments, and the above embodiments and the description only illustrate the principle of the present disclosure. Various changes and modifications may be made to the present disclosure without departing from the spirit and scope of the present disclosure, and such changes and modifications all fall within the claimed scope of the present disclosure. The claimed protection scope of the present disclosure is defined by the appended claims and equivalents thereof.