CN103363962A - Remote sensing evaluation method of lake water reserves based on multispectral images - Google Patents

Abstract

The invention provides a remote sensing evaluation method of lake water reserves based on multispectral images. The remote sensing evaluation method includes: according to lake water body spectral response characteristics in different bands of multispectral remote sensing images, extracting water body indexes which reflect lake water surface distribution characteristics; according to the extracted water body indexes, obtaining lake water body boundaries, giving actually-measured lake water level information to the lake water body boundaries, and generating lake equal water level line data sets according to different periods of the lake water body boundaries; according to the lake equal water level line data sets, simulating lake underwater topography; and calculating the lake water reserves in different periods according to the simulated lake underwater topography and actually-measured water level data. According to the remote sensing evaluation method, the lake underwater topography is monitored through utilization of multispectral satellite remote sensing data, the utilized satellite remote sensing data substantially can be obtained for free, and the remote sensing data has the advantages of wide space coverage areas and rapid refreshing speed, so that compared with conventional field actual-measurement methods, the remote sensing evaluation method provided by the embodiment of the invention has the advantages of low monitoring cost, convenient refreshing, and possibility of wide-range popularization and application.

Description

A kind of Lake Water reserves remote sensing estimation method based on multispectral image

Technical field

The present invention relates to satellite remote sensing and earth observation field, relate in particular to a kind of Lake Water reserves remote sensing estimation method based on multispectral image.

Background technology

Lake (comprising natural lake and the artificial pool, storehouse) is the main form of expression of surface water resources.Except glacier and permanent snow lid, Lake Water is second largest surface water resources type, and it is one of principal element that affects the global sea variation.Lake evolution and environmental change is closely related: on the one hand, the water yield that seiche causes changes and the sedimentary environment change can be indicated climate change sensitively; On the other hand, the increase and decrease of Lake Water area can change the underlying surface condition, thereby climate change is exerted an influence.Yet for a long time, investigated the restriction with observational data, Simulation and analysis method, people can't accurately obtain on a large scale the water resources in lake reserves and describe its dynamic change characterization, had a strong impact on people to land table water reserve variation and the whole world or regional water circulation aggravate, the understanding of the problem in science such as relation of global sea lifting.In the water resources management application facet, this present situation has also limited people to the reply timeliness of the drought and waterlogging event that takes place frequently in the global range in recent years.

Therefore, need to adopt effective method that the water reserve in lake is accurately estimated.Existing evaluation method is: the water level take actual measurement underwater topography data as fundamental construction-volume curve method.The method mainly is to carry out the calculating of volume according to actual measurement underwater topography data and the corresponding mathematical model of level measuring the data, thereby obtain water level-volume curve equation, after knowing the lake real time water level, just can learn the volumetric values that this water level is corresponding according to water level-volume curve.

The application of said method is subject to obtaining of storage capacity computational mathematics model, level measuring spatial representative and these three aspect data of underwater topography.

Storage capacity computing method commonly used have the method for section, level line volumetric method, square grid method and Triangular meshes method, are subjected to the impact of lake form, size and complexity, and inappropriate mathematical model can be brought the larger error of calculation.This problem can be utilized current high-performance calculation resource, solves by increasing model differential magnitude.For the poor problem of water-level observation station measurement result spatial representative, then can solve by optimizing the observation site location and increasing research station point quantity.

And for the underwater topography data, there are early stage l:5000 or 1:10000 topomap in general area, can by collecting sweep vector, gather altitude figures and generating digital relief block (DTM).The DTM data duration that the method generates is short, but the collection of these basic datas is generally comparatively difficult, the storage capacity data precision that obtains is subject to the impact of the precision of topomap own, and can't reflect the impact that drawing later stage upper water sand and mankind's activity are transformed underwater topography.The way that addresses this problem is comprehensive utilization 3S technology, be GPS, GIS, RS, regularly carry out the lake bathymetric surveying, as the underwater topography that fully utilizes GPS and laser radar data is charted, based on the underwater topography simulation of GPS field study data, and the digital elevation model (DEM) that directly utilizes laser radar to obtain is simulated lake underwater topography etc.Although these methods can partly solve underwater topography data acquisition and replacement problem, the restriction of being measured duration and cost can't be applied on a large scale.

Summary of the invention

The technical matters that (one) will solve

The objective of the invention is to propose a kind of Lake Water estimation method of reserve, so that the method is applied in can be on a large scale.

(2) technical scheme

In order to solve the problems of the technologies described above, the present invention proposes a kind of Lake Water reserves remote sensing estimation method based on multispectral image, the method comprises:

S1, extract the water body index of reflection lake distribution characteristics according to the spectral response characteristic of water body in lake on the multi-spectrum remote sensing image different-waveband;

S2, obtain the water body in lake border according to the water body index that extracts, and the lake level information of actual measurement is assigned to the water body in lake border, generate lake water table contour data set according to the water body in lake border of different times;

S3, according to lake water table contour data set simulation lake underwater topography;

S4, the Lake Water reserves that calculate different times according to lake underwater topography and the waterlevel data of actual measurement of simulation.

Wherein, before step S1, at first carry out the collection of remotely-sensed data and select; And the remote sensing image that will collect and select carries out radiation correcting and geometric correction; Then extract the water body index of reflection lake distribution characteristics according to the remote sensing image after correcting.

Wherein, the extraction of water body index can adopt two kinds of approach to realize:

Extract normalization difference water body index according to the spectral response characteristics of water body in lake on Landsat MSS or HJ-1A/B image different-waveband;

Perhaps,

Extract enhancement mode normalization difference water body index according to the spectral response characteristics of water body in lake on Landsat TM/ETM+ image different-waveband.

Wherein, normalization difference water body index NDWI adopts following formula to calculate:

$\mathrm{NDWI}=\frac{{\mathrm{\ρ}}_{\mathrm{Green}}-{\mathrm{\ρ}}_{\mathrm{NIR}}}{{\mathrm{\ρ}}_{\mathrm{Green}}+{\mathrm{\ρ}}_{\mathrm{NIR}}}$

Wherein, ρ _GreenBe the green light band on Landsat MSS or the HJ-1A/B image, ρ _NIRBe the near-infrared band on Landsat MSS or the HJ-1A/B image.

Wherein, enhancement mode normalization difference water body index MNDWI adopts following formula to calculate:

$\mathrm{MNDWI}=\frac{{\mathrm{\ρ}}_{\mathrm{Green}}-{\mathrm{\ρ}}_{\mathrm{SWIR}}}{{\mathrm{\ρ}}_{\mathrm{Green}}+{\mathrm{\ρ}}_{\mathrm{SWIR}}}$

Wherein, ρ _GreenBe the green light band on the Landsat TM/ETM+ image, ρ _SWIRBe the short infrared wave band on the Landsat TM/ETM+ image.

Preferably, adopt improved maximum between-cluster variance threshold method to extract the lake raster data according to the water body index that extracts; Turn tool vector by grid in the ArcMap9.3 software and convert the lake raster data that extracts to water body in lake boundary vector data.

Wherein, according to the tonal range of the water body index C that extracts 0,1 ..., d-1} is divided into C with C ₁And C ₂Two classes, the gray threshold of cutting apart are t, i.e. the water boundary threshold value in lake, and t satisfies:

$t=\arg \underset{0\≤t\≤d-1}{\max }\frac{S^2}{S_1^2+S_2^2}$

Wherein, S ²Be C ₁And C ₂Inter-class variance, S ₁ ²Be C ₁The class internal variance, S ₂ ²Be C ₂The class internal variance;

Wherein,

S ²=P ₁(A ₁-A) ²+P ₂(A ₂-A) ²

$S_1^2=\frac{{\mathrm{\Σ}}_{i=0}^t{(i-A_1)}^2{p}_i}{P_1}$

$S_2^2=\frac{{\mathrm{\Σ}}_{i=t+1}^{d-1}{(i-A_2)}^2{p}_i}{P_2}$

Wherein, P ₁Be C ₁Middle pixel number accounts for the ratio of total pixel number among the water body index C, P ₂Be C ₂Middle pixel number accounts for the ratio of total pixel number among the water body index C, and i is the pixel gray-scale value, and A is the pixel average gray of C, A ₁Be C ₁The pixel average gray, A ₂Be C ₂The pixel average gray, p _iFor gray-scale value is the ratio that the pixel number of i accounts for total pixel number among the water body index C.

Wherein, adopt TIN simulation lake underwater topography according to lake water table contour data set, the reticulate texture that to be the lake underwater topography be comprised of the sealene triangle of a series of non-overlapping copies represents, each leg-of-mutton each node comprises independently latitude and longitude coordinates, water depth value and the angle of gradient;

According to the set that underwater topography and the waterlevel data of actual measurement of simulation is simplified to water body in lake a series of triangular prisms, the volume of whole water body in lake then be each triangular prism volume with, computing formula is:

$V=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{S}_i\frac{h_i+h_{i+1}+h_{i+2}}{3}$

Wherein, n is the number of triangular prism, S _iBe the area of i triangular prism water body upper surface, h _i, h _I+1And h _I+2Water depth value for three incline positions of triangular prism.

Preferably, simulate the lake underwater topography according to waterlevel data centralized procurements such as lakes with the TIN creation module in the three dimensional analysis instrument of ArcMap9.3 software.

Preferably, the Lake Water reserves utilize area in the ArcMap9.3 software three dimensional analysis instrument-volume statistical module to calculate.

(3) beneficial effect

Utilize multispectral satellite remote sensing date that the lake underwater topography is monitored in the embodiments of the invention, because the satellite remote sensing date that uses basically can Free Acquisition, and the spatial coverage of remotely-sensed data is wide, renewal speed is fast, therefore, compare traditional fieldwork method, it is low that the method that the embodiment of the invention proposes has a monitoring cost, upgrade convenient and can be on a large scale in the advantage applied.

In addition, the terrain remote sensing monitoring result is the Lake Water reserves remote sensing estimation method on basis under Lake Water, can effectively eliminate different times lake underwater topography itself and change estimating result's impact, and it is more reliable that its precision is compared classic method.

Description of drawings

In order to be illustrated more clearly in the embodiment of the invention or technical scheme of the prior art, the below will do to introduce simply to the accompanying drawing of required use in embodiment or the description of the Prior Art, apparently, accompanying drawing in the following describes only is some embodiments of the present invention, for those of ordinary skills, under the prerequisite of not paying creative work, can also obtain according to these accompanying drawings other accompanying drawing.

Fig. 1 is based on the Lake Water estimation method of reserve concept map of multi-spectrum remote sensing image in the embodiments of the invention;

Fig. 2 is based on the process flow diagram of the Lake Water estimation method of reserve of multi-spectrum remote sensing image in the embodiments of the invention.

Embodiment

Below in conjunction with drawings and Examples embodiments of the present invention are described in further detail.Following examples are used for explanation the present invention, but can not be used for limiting the scope of the invention.

For lake water reserve and situation of change thereof in can grasping better on a large scale, the embodiment of the invention has proposed a kind of Lake Water reserves remote sensing estimation method based on multispectral image, and referring to Fig. 1, the method comprises:

S1, extract the water body index of reflection lake distribution characteristics according to the spectral response characteristic of water body in lake on the multi-spectrum remote sensing image different-waveband;

S2, obtain the water body in lake border according to the water body index that extracts, and the lake level information of actual measurement is assigned to the water body in lake border, generate lake water table contour data set according to the water body in lake border of different times;

S3, according to lake water table contour data set simulation lake underwater topography;

S4, the Lake Water reserves that calculate different times according to lake underwater topography and the waterlevel data of actual measurement of simulation.

In the above embodiment of the present invention, utilize multispectral satellite remote sensing date that the lake underwater topography is monitored, because the satellite remote sensing date that uses basically can Free Acquisition, and the spatial coverage of remotely-sensed data is wide, renewal speed is fast, therefore, compare traditional fieldwork method, it is low that the method that the embodiment of the invention proposes has a monitoring cost, upgrade convenient and can be on a large scale in the advantage applied.

In the above embodiment of the present invention, because different satellites can provide the multi-spectrum remote sensing image of different times, different quality, therefore, method needs needed remote sensing image is collected and selected before using; In addition, in order to eliminate cloud and mist covering, attitude of satellite variation and topographic relief etc. to the impact of satellite imagery, before using remote sensing image, need remote sensing image is carried out radiation correcting and geometric correction, then extract the water body index that can reflect the lake distribution characteristics according to the remote sensing image after correcting, referring to Fig. 2.

The water body index that uses in the above embodiment of the present invention comprises normalization difference water body index NDWI and enhancement mode normalization difference water body index MNDWI, and wherein NDWI is applicable to Landsat MSS and HJ-1A/B image, and computing formula is as follows:

$\mathrm{NDWI}=\frac{{\mathrm{\ρ}}_{\mathrm{Green}}-{\mathrm{\ρ}}_{\mathrm{NIR}}}{{\mathrm{\ρ}}_{\mathrm{Green}}+{\mathrm{\ρ}}_{\mathrm{NIR}}}---\left(1\right)$

In the formula, ρ _GreenBe the green light band on Landsat MSS or the HJ-1A/B image, ρ _NIRBe the near-infrared band on Landsat MSS or the HJ-1A/B image.

And MNDWI is applicable to Landsat TM/ETM+ image, and its computing formula is as follows:

$\mathrm{MNDWI}=\frac{{\mathrm{\ρ}}_{\mathrm{Green}}-{\mathrm{\ρ}}_{\mathrm{SWIR}}}{{\mathrm{\ρ}}_{\mathrm{Green}}+{\mathrm{\ρ}}_{\mathrm{SWIR}}}---\left(2\right)$

In the formula, ρ _GreenBe the green light band on the Landsat TM/ETM+ image, ρ _SWIRBe the short infrared wave band on the Landsat TM/ETM+ image.

In the above embodiment of the present invention, adopt improved maximum between-cluster variance threshold method to extract the lake raster data according to the water body index that extracts; Turn tool vector by grid in the ArcMap9.3 software and convert the lake raster data that extracts to water body in lake boundary vector data.

In the above embodiment of the present invention, the ultimate principle of improved maximum between-cluster variance threshold method is as follows:

According to the tonal range of the water body index C that extracts 0,1 ..., d-1} is divided into C with C ₁And C ₂Two classes, the gray threshold of cutting apart are t, i.e. the water boundary threshold value in lake, and t satisfies:

$t=\arg \underset{0\≤t\≤d-1}{\max }\frac{S^2}{S_1^2+S_2^2}---\left(3\right)$

Wherein, S ²Be C ₁And C ₂Inter-class variance, S ₁ ²Be C ₁The class internal variance, S ₂ ²Be C ₂The class internal variance;

Wherein,

S ²=P ₁(A ₁-A) ²+P ₂(A ₂-A) ² (4)

$S_1^2=\frac{{\mathrm{\Σ}}_{i=0}^t{(i-A_1)}^2{p}_i}{P_1}---\left(5\right)$

$S_2^2=\frac{{\mathrm{\Σ}}_{i=t+1}^{d-1}{(i-A_2)}^2{p}_i}{P_2}---\left(6\right)$

Wherein, P ₁Be C ₁Middle pixel number accounts for the ratio of total pixel number among the water body index C, P ₂Be C ₂Middle pixel number accounts for the ratio of total pixel number among the water body index C, and i is the pixel gray-scale value, and A is the pixel average gray of C, A ₁Be C ₁The pixel average gray, A ₂Be C ₂The pixel average gray, p _iFor gray-scale value is the ratio that the pixel number of i accounts for total pixel number among the water body index C.

In the above embodiment of the present invention, obtain the waterlevel data collection such as lake according to remote sensing image after, when the underwater topography more complicated of lake, can adopt the method simulation lake underwater topography of TIN (TIN), the reticulate texture that to be the lake underwater topography be comprised of the sealene triangle of a series of non-overlapping copies represents, each leg-of-mutton each node comprises independently latitude and longitude coordinates, water depth value and the angle of gradient; The simulation process of underwater topography can adopt the TIN creation module in the three dimensional analysis instrument of ArcMap9.3 software to simulate, with the waterlevel data collection such as lake that the obtain input data as software, behind the operating software, obtain using the lake underwater topography of TIN modeling.

Terrain data is the basis under the Lake Water that obtains, and water body in lake can be simplified to the set of a series of triangular prisms in conjunction with the waterlevel data of actual measurement, the volume of whole water body in lake then be each triangular prism volume with, computing formula is:

$V=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{S}_i\frac{h_i+h_{i+1}+h_{i+2}}{3}---\left(7\right)$

Wherein, n is the number of triangular prism, S _iBe the area of i triangular prism water body upper surface, h _i, h _I+1And h _I+2Water depth value for three incline positions of triangular prism.

Above-mentioned estimation to the Lake Water reserves can adopt area in the ArcMap9.3 software three dimensional analysis instrument-volume statistical module to calculate, and lake underwater topography and lake level data by input obtains behind the operating software, can obtain Lake Water reserves value.

When the lake underwater topography is relatively simple, can adopt square grid simulation lake underwater topography.

Therefore, in the time need to estimating the water reserve in certain lake, can carry out by the following method:

1, obtains the multispectral satellite remote-sensing image of the different time sections that covers lake region, such as Landsat MSS/TM/ETM+ and HJ-1A/B;

2, the remote sensing image that obtains is carried out radiant correction and geometry correction;

3, utilize the water body index of the image data calculating lake country different times after proofreading and correct, comprise NDWI and MNDWI, (1) and (2) formula of employing is calculated;

4, adopt improved maximum between-cluster variance threshold method to extract the lake raster dataset from the water body index data centralization, grid in the recycling ArcMap9.3 software turns tool vector and converts the lake raster dataset to water body in lake boundary vector data set, again the lake level data of the different times collected are assigned to the water boundary vector of corresponding time period, can obtain the lake country water table contour data set of different times;

5, according to the water table contour data set that obtains, adopt the TIN analogy method, utilize the TIN creation module simulation underwater topography in the ArcMap9.3 software three dimensional analysis instrument;

6, according to the underwater topography of simulation and the waterlevel data of actual measurement, utilize area in the ArcMap9.3 software three dimensional analysis instrument-volume statistical module to calculate the water reserve in lake.

In the above embodiment of the present invention, for different zones, the satellite remote-sensing image data of employing may be different, such as remote sensing images such as CBERS, ASTER, MODIS, SPOT; And the using method of actual measurement lake level data is also different, and what use in the embodiments of the invention is the average of lake country measured water level, and for the abundant zone of water-level observation data, can adopt different water level values on the same water boundary.

Therefore, the invention has the beneficial effects as follows:

Utilize multispectral satellite remote sensing date that the lake underwater topography is monitored in the embodiments of the invention, because the satellite remote sensing date that uses basically can Free Acquisition, and the remotely-sensed data spatial coverage is wide, renewal speed is fast, therefore, compare traditional fieldwork method, it is low that the method that the embodiment of the invention proposes has a monitoring cost, upgrade convenient and can be on a large scale in the advantage applied.

In addition, the terrain remote sensing monitoring result is the Lake Water reserves remote sensing estimation method on basis under the Lake Water, can effectively eliminate different times lake underwater topography itself and change estimating result's impact, and it is more reliable that its precision is compared classic method.

Embodiments of the invention provide for example with for the purpose of describing, and are not exhaustively or limit the invention to disclosed form.Many modifications and variations are apparent for the ordinary skill in the art.Selecting and describing embodiment is for better explanation principle of the present invention and practical application, thereby and those of ordinary skill in the art can understand the various embodiment with various modifications that the present invention's design is suitable for special-purpose.

Claims (10)

1. the Lake Water reserves remote sensing estimation method based on multispectral image is characterized in that, comprising:

S1, extract the water body index of reflection lake distribution characteristics according to the spectral response characteristic of water body in lake on the multi-spectrum remote sensing image different-waveband;

S2, obtain the water body in lake border according to the water body index that extracts, and the lake level information of actual measurement is assigned to the water body in lake border, generate lake water table contour data set according to the water body in lake border of different times;

S3, according to lake water table contour data set simulation lake underwater topography;

S4, the Lake Water reserves that calculate different times according to lake underwater topography and the waterlevel data of actual measurement of simulation.

2. the Lake Water reserves remote sensing estimation method based on multispectral image according to claim 1 is characterized in that, described step S1 further comprises:

Collect and select remote sensing image;

The remote sensing image of collecting and select is carried out radiation correcting and geometric correction;

Extract the water body index of reflection lake distribution characteristics according to the remote sensing image after correcting.

3. the Lake Water reserves remote sensing estimation method based on multispectral image according to claim 1, it is characterized in that, extract the water body index that reflects the lake distribution characteristics according to the spectral response characteristic of water body in lake on the multi-spectrum remote sensing image different-waveband described in the step S1 and comprise:

Extract normalization difference water body index according to the spectral response characteristic of water body in lake on Landsat MSS or HJ-1A/B image different-waveband;

Perhaps,

Extract enhancement mode normalization difference water body index according to the spectral response characteristic of water body in lake on Landsat TM/ETM+ image different-waveband.

4. the Lake Water reserves remote sensing estimation method based on multispectral image according to claim 3 is characterized in that, described normalization difference water body index NDWI adopts following formula to calculate:

$\mathrm{NDWI}=\frac{{\mathrm{\ρ}}_{\mathrm{Green}}-{\mathrm{\ρ}}_{\mathrm{NIR}}}{{\mathrm{\ρ}}_{\mathrm{Green}}+{\mathrm{\ρ}}_{\mathrm{NIR}}}$

Wherein, ρ _GreenBe the green light band on Landsat MSS or the HJ-1A/B image, ρ _NIRBe the near-infrared band on Landsat MSS or the HJ-1A/B image.

5. the Lake Water reserves remote sensing estimation method based on multispectral image according to claim 3 is characterized in that, described enhancement mode normalization difference water body index MNDWI adopts following formula to calculate:

$\mathrm{MNDWI}=\frac{{\mathrm{\ρ}}_{\mathrm{Green}}-{\mathrm{\ρ}}_{\mathrm{SWIR}}}{{\mathrm{\ρ}}_{\mathrm{Green}}+{\mathrm{\ρ}}_{\mathrm{SWIR}}}$

Wherein, ρ _GreenBe the green light band on the Landsat TM/ETM+ image, ρ _SWIRBe the short infrared wave band on the Landsat TM/ETM+ image.

6. the Lake Water reserves remote sensing estimation method based on multispectral image according to claim 1 is characterized in that, obtains the water body in lake border according to the water body index that extracts among the step S2 and comprises:

Adopt improved maximum between-cluster variance threshold method to extract the lake raster data according to the water body index that extracts;

Turn tool vector by grid in the ArcMap9.3 software and convert the lake raster data that extracts to water body in lake boundary vector data.

7. the Lake Water reserves remote sensing estimation method based on multispectral image according to claim 6 is characterized in that, described water body index according to extracting adopts improved maximum between-cluster variance threshold method extraction lake raster data to comprise:

According to the tonal range of the water body index C that extracts 0,1 ..., d-1} is divided into C with C ₁And C ₂Two classes, the gray threshold of cutting apart are t, i.e. the water boundary threshold value in lake, and t satisfies:

$t=\arg \underset{0\≤t\≤d-1}{\max }\frac{S^2}{S_1^2+S_2^2}$

Wherein, S ²Be C ₁And C ₂Inter-class variance, S ₁ ²Be C ₁The class internal variance, S ₂ ²Be C ₂The class internal variance;

Wherein,

S ²=P ₁(A ₁-A) ²+P ₂(A ₂-A) ²

$S_1^2=\frac{{\mathrm{\Σ}}_{i=0}^t{(i-A_1)}^2{p}_i}{P_1}$

$S_2^2=\frac{{\mathrm{\Σ}}_{i=t+1}^{d-1}{(i-A_2)}^2{p}_i}{P_2}$

Wherein, P ₁Be C ₁Middle pixel number accounts for the ratio of total pixel number among the water body index C, P ₂Be C ₂Middle pixel number accounts for the ratio of total pixel number among the water body index C, and i is the pixel gray-scale value, and A is the pixel average gray of C, A ₁Be C ₁The pixel average gray, A ₂Be C ₂The pixel average gray, p _iFor gray-scale value is the ratio that the pixel number of i accounts for total pixel number among the water body index C.

8. the described Lake Water reserves remote sensing estimation method based on multispectral image of any one is characterized in that according to claim 1～7, comprises according to lake water table contour data set simulation lake underwater topography described in the step S3:

Adopt TIN simulation lake underwater topography according to lake water table contour data set, the reticulate texture that to be the lake underwater topography be comprised of the sealene triangle of a series of non-overlapping copies represents, each leg-of-mutton each node comprises independently latitude and longitude coordinates, water depth value and the angle of gradient;

The Lake Water reserves that calculate different times according to lake underwater topography and the waterlevel data of actual measurement of simulation described in the step S4 comprise:

Water body in lake is simplified to the set of a series of triangular prisms according to the waterlevel data of the underwater topography of simulating and actual measurement, the volume of whole water body in lake then is each triangular prism volume sum, and computing formula is:

$V=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{S}_i\frac{h_i+h_{i+1}+h_{i+2}}{3}$

Wherein, n is the number of triangular prism, S _iBe the area of i triangular prism water body upper surface, h _i, h _I+1And h _I+2Water depth value for three incline positions of triangular prism.

9. the Lake Water reserves remote sensing estimation method based on multispectral image according to claim 8 is characterized in that, describedly comprises with TIN simulation lake underwater topography according to waterlevel data centralized procurements such as lakes:

Simulate the lake underwater topography according to waterlevel data centralized procurements such as lakes with the TIN creation module in the three dimensional analysis instrument of ArcMap9.3 software.

10. the Lake Water reserves remote sensing estimation method based on multispectral image according to claim 8 is characterized in that, described Lake Water reserves utilize area in the ArcMap9.3 software three dimensional analysis instrument-volume statistical module to calculate.

Images