CN103363962B - 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, particularly relate to a kind of Lake Water reserves remote sensing estimation method based on multispectral image.

Background technology

Lake (comprising natural lake and the pool, artificial 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 affecting global sea change.Lake evolution and environmental change is closely related: on the one hand, and the water yield change that seiche causes and sedimentary environment change and can indicate climate change sensitively; On the other hand, the increase and decrease of Lake Water area can change land surface condition, thus has an impact to climate change.But for a long time, by investigation and observational data, Simulation and analysis method restriction, people accurately cannot obtain water resources in lake reserves describe its dynamic change characterization on a large scale, had a strong impact on people land table water reserve change is circulated with the whole world or regional water aggravate, the understanding of the problem in science such as the relation of global sea-level changes.In water resources management application aspect, this present situation also limit the reply timeliness of people to the drought and waterlogging event taken place frequently in global range in recent years.

Therefore, need to adopt the water reserve of effective method to lake accurately to estimate.Existing evaluation method is: to survey the water level-volume curve method built based on underwater topography data.The method mainly carries out the calculating of volume according to actual measurement underwater topography data and the corresponding mathematical model of level measuring data acquisition, thus obtain water level-volume curve equation, after knowing lake real time water level, just can learn according to water level-volume curve the volumetric values that this water level is corresponding.

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

Conventional storage capacity computing method have the method for section, level line volumetric method, square mesh method and Triangular meshes method, and by the impact of lake form, size and complexity, inappropriate mathematical model can bring the larger error of calculation.This problem can utilize current high-performance calculation resource, solves by increasing model differential magnitude.For the problem of gaging station measurement result spatial representative difference, then can solve by optimizing observation site location and increasing research station point quantity.

And for underwater topography data, there is early stage l:5000 or 1:10000 topomap in general area, by collecting scanning vector, gathers altitude figures and generating digital terrain model (DTM).The DTM data duration that the method generates is short, but the general comparatively difficulty of the collection of these basic datas, the storage capacity data precision obtained is subject to the impact of the precision of topomap own, and cannot reflect the impact that drawing later stage upper water sand and mankind's activity are transformed underwater topography.The way addressed this problem is comprehensive utilization 3S technology, i.e. GPS, GIS, RS, regularly carry out lake bathymetric surveying, as fully utilized the underwater topography drawing of GPS and laser radar data, underwater topography based on GPS field study data is simulated, and the digital elevation model (DEM) directly utilizing laser radar to obtain simulates lake underwater topography etc.Although these methods can partly solve underwater topography data acquisition and replacement problem, by the restriction of measuring duration and cost, cannot apply on a large scale.

Summary of the invention

(1) technical matters that will solve

The object of the invention is to propose a kind of Lake Water estimation method of reserve, make the method can on a large scale in apply.

(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 multi-spectrum remote sensing image different-waveband;

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

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

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

Wherein, before step S1, first carry out the collection of remotely-sensed data and select; And the remote sensing image collected and select is carried out radiation correcting and geometric correction; Then the water body index of reflection lake distribution characteristics is extracted according to the remote sensing image after correction.

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

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

Or,

Enhancement mode normalization difference water body index is extracted 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, ρ _greenfor the green light band on Landsat MSS or HJ-1A/B image, ρ _nIRfor the near-infrared band on Landsat MSS or 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, ρ _greenfor the green light band on Landsat TM/ETM+ image, ρ _sWIRfor the short infrared wave band on Landsat TM/ETM+ image.

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

Wherein, according to extract water body index C tonal range 0,1 ..., d-1}, is divided into C by C ₁and C ₂two classes, the gray threshold of segmentation is t, i.e. the water boundary threshold value in lake, and t meets:

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

Wherein, S ²for C ₁and C ₂inter-class variance, S ₁ ²for C ₁variance within clusters, S ₂ ²for C ₂variance within clusters;

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 ₁for C ₁middle pixel number accounts for the ratio of total pixel number in water body index C, P ₂for C ₂middle pixel number accounts for the ratio of total pixel number in water body index C, and i is pixel gray-scale value, and A is the pixel average gray of C, A ₁for C ₁pixel average gray, A ₂for C ₂pixel average gray, p _ifor the gray-scale value pixel number that is i accounts for the ratio of total pixel number in water body index C.

Wherein, TIN simulation lake underwater topography is adopted according to lake water table contour data set, namely the reticulate texture that lake underwater topography is made up 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 underwater topography of simulation and the waterlevel data of actual measurement, water body in lake is simplified to the set of a series of triangular prism, the volume of whole water body in lake be then 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-th triangular prism water body upper surface, h _i, h _i+1and h _i+2for the water depth value of triangular prism three incline positions.

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

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

(3) beneficial effect

Multispectral satellite remote sensing date is utilized to monitor lake underwater topography in embodiments of the invention, because the satellite remote sensing date used substantially can Free Acquisition, and the spatial coverage of remotely-sensed data is wide, renewal speed is fast, therefore, compare traditional fieldwork method, to have monitoring cost low for the method that the embodiment of the present invention proposes, upgrade convenient and can on a large scale in carry out the advantage applied.

In addition, the Lake Water reserves remote sensing estimation method under Lake Water based on terrain remote sensing monitoring result, effectively can eliminate the underwater topography change of different times lake itself to the impact of estimation result, it is more reliable that its precision compares classic method.

Accompanying drawing explanation

In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.

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

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

Embodiment

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

In order to grasp lake water reserve and situation of change thereof on a large scale better, the embodiment of the present invention proposes a kind of Lake Water reserves remote sensing estimation method based on multispectral image, and see 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 multi-spectrum remote sensing image different-waveband;

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

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

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

In the above embodiment of the present invention, multispectral satellite remote sensing date is utilized to monitor lake underwater topography, because the satellite remote sensing date used substantially can Free Acquisition, and the spatial coverage of remotely-sensed data is wide, renewal speed is fast, therefore, compare traditional fieldwork method, to have monitoring cost low for the method that the embodiment of the present invention proposes, upgrade convenient and can on a large scale in carry out the advantage applied.

In the above embodiment of the present invention, because different satellite can provide the multi-spectrum remote sensing image of different times, different quality, therefore, method needs to collect required remote sensing image and select before using; In addition, in order to eliminate the impact on satellite imagery such as cloud and mist covering, attitude of satellite change and topographic relief, needed to carry out radiation correcting and geometric correction to remote sensing image before use remote sensing image, then the water body index that can reflect lake distribution characteristics is extracted according to the remote sensing image after correction, see Fig. 2.

The water body index used 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 formula, ρ _greenfor the green light band on Landsat MSS or HJ-1A/B image, ρ _nIRfor the near-infrared band on Landsat MSS or HJ-1A/B image.

And MNDWI is applicable to Landsat TM/ETM+ image, 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 formula, ρ _greenfor the green light band on Landsat TM/ETM+ image, ρ _sWIRfor the short infrared wave band on Landsat TM/ETM+ image.

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

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

According to extract water body index C tonal range 0,1 ..., d-1}, is divided into C by C ₁and C ₂two classes, the gray threshold of segmentation is t, i.e. the water boundary threshold value in lake, and t meets:

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

Wherein, S ²for C ₁and C ₂inter-class variance, S ₁ ²for C ₁variance within clusters, S ₂ ²for C ₂variance within clusters;

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 ₁for C ₁middle pixel number accounts for the ratio of total pixel number in water body index C, P ₂for C ₂middle pixel number accounts for the ratio of total pixel number in water body index C, and i is pixel gray-scale value, and A is the pixel average gray of C, A ₁for C ₁pixel average gray, A ₂for C ₂pixel average gray, p _ifor the gray-scale value pixel number that is i accounts for the ratio of total pixel number in 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 lake underwater topography more complicated, the method simulation lake underwater topography of TIN (TIN) can be adopted, namely the reticulate texture that lake underwater topography is made up 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, using the input data of the waterlevel data collection such as the lake that obtains as software, after operating software, obtain the lake underwater topography using Triangulated irregular network model simulation.

Under the Lake Water obtained based on terrain data, water body in lake can be simplified to the set of a series of triangular prism by the waterlevel data in conjunction with actual measurement, the volume of whole water body in lake be then 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-th triangular prism water body upper surface, h _i, h _i+1and h _i+2for the water depth value of triangular prism three incline positions.

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

When lake underwater topography is relatively simple, square grid can be adopted to simulate lake underwater topography.

Therefore, when needing to estimate the water reserve in certain lake, can carry out by the following method:

1, the multispectral satellite remote-sensing image of the different time sections covering lake region is obtained, as Landsat MSS/TM/ETM+ and HJ-1A/B;

2, radiant correction and geometry correction are carried out to the remote sensing image obtained;

3, utilize the image data after correcting to calculate the water body index of lake country different times, comprise NDWI and MNDWI, (1) and (2) formula of employing calculates;

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

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

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

In the above embodiment of the present invention, for different regions, the satellite remote-sensing image data of employing may be different, as remote sensing images such as CBERS, ASTER, MODIS, SPOT; And the using method of surveying lake level data is also different, in embodiments of the invention, use the average of lake country measured water level, and for the region that water-level observation data are enriched, same water boundary can adopt different water level value.

Therefore, the invention has the beneficial effects as follows:

Multispectral satellite remote sensing date is utilized to monitor lake underwater topography in embodiments of the invention, because the satellite remote sensing date used substantially can Free Acquisition, and remotely-sensed data spatial coverage is wide, renewal speed is fast, therefore, compare traditional fieldwork method, to have monitoring cost low for the method that the embodiment of the present invention proposes, upgrade convenient and can on a large scale in carry out the advantage applied.

In addition, the Lake Water reserves remote sensing estimation method under Lake Water based on terrain remote sensing monitoring result, effectively can eliminate the underwater topography change of different times lake itself to the impact of estimation result, it is more reliable that its precision compares classic method.

Embodiments of the invention provide in order to 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 in order to principle of the present invention and practical application are better described, and those of ordinary skill in the art can understand the present invention thus design the various embodiments with various amendment being suitable for special-purpose.

Claims (8)

1., based on a Lake Water reserves remote sensing estimation method for multispectral image, it 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 multi-spectrum remote sensing image different-waveband;

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

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

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

Obtain water body in lake border according to the water body index extracted in step S2 to comprise:

The maximum between-cluster variance threshold method improved is adopted to extract lake raster data according to the water body index extracted;

Turn tool vector by grid in ArcMap 9.3 software and convert the lake raster data of extraction to water body in lake boundary vector data;

Wherein, the described water body index according to extracting adopts the maximum between-cluster variance threshold method extraction lake raster data improved to comprise:

According to extract water body index C tonal range 0,1 ..., d-1}, is divided into C by C ₁and C ₂two classes, the gray threshold of segmentation is t, i.e. the water boundary threshold value in lake, and t meets:

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

S ²for C ₁and C ₂inter-class variance, S ₁ ²for C ₁variance within clusters, S ₂ ²for C ₂variance within clusters;

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 ₁for C ₁middle pixel number accounts for the ratio of total pixel number in water body index C, P ₂for C ₂middle pixel number accounts for the ratio of total pixel number in water body index C, and i is pixel gray-scale value, and A is the pixel average gray of C, A ₁for C ₁pixel average gray, A ₂for C ₂pixel average gray, p _ifor the gray-scale value pixel number that is i accounts for the ratio of total pixel number in water body index C.

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

Collect and select remote sensing image;

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

The water body index reflecting lake distribution characteristics is extracted 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 according to the spectral response characteristic of water body in lake on multi-spectrum remote sensing image different-waveband described in step S1 and reflect that the water body index of lake distribution characteristics comprises:

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

Or,

Enhancement mode normalization difference water body index is extracted 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, ρ _greenfor the green light band on Landsat MSS or HJ-1A/B image, ρ _nIRfor the near-infrared band on Landsat MSS or 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, ρ _greenfor the green light band on Landsat TM/ETM+ image, ρ _sWIRfor the short infrared wave band on Landsat TM/ETM+ image.

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

TIN simulation lake underwater topography is adopted according to lake water table contour data set, namely the reticulate texture that lake underwater topography is made up 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;

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

According to the underwater topography of simulation and the waterlevel data of actual measurement, water body in lake is simplified to the set of a series of triangular prism, the volume of whole water body in lake is then 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-th triangular prism water body upper surface, h _i, h _i+1and h _i+2for the water depth value of triangular prism three incline positions.

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

According to the TIN creation module simulation lake underwater topography of the waterlevel data centralized procurements such as lake in the three dimensional analysis instrument of ArcMap 9.3 software.

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