1.一種基于多光譜影像的湖泊水儲量遙感估算方法，其特征在于，包括：

S1、根據(jù)湖泊水體在多光譜遙感影像不同波段上的光譜響應(yīng)特性提取反映湖泊水面分布特征的水體指數(shù)；

S2、根據(jù)提取的水體指數(shù)獲取湖泊水體邊界，并將實測的湖泊水位信息賦給湖泊水體邊界，根據(jù)不同時期的湖泊水體邊界生成湖泊等水位線數(shù)據(jù)集；

S3、根據(jù)湖泊等水位線數(shù)據(jù)集模擬湖泊水下地形；

S4、根據(jù)模擬的湖泊水下地形和實測的水位數(shù)據(jù)計算不同時期的湖泊水儲量。

2.根據(jù)權(quán)利要求1所述的基于多光譜影像的湖泊水儲量遙感估算方法，其特征在于，所述步驟S1進(jìn)一步包括：

收集和遴選遙感影像；

將收集和遴選的遙感影像進(jìn)行輻射糾正和幾何糾正；

根據(jù)糾正后的遙感影像提取反映湖泊水面分布特征的水體指數(shù)。

3.根據(jù)權(quán)利要求1所述的基于多光譜影像的湖泊水儲量遙感估算方法，其特征在于，步驟S1中所述根據(jù)湖泊水體在多光譜遙感影像不同波段上的光譜響應(yīng)特性提取反映湖泊水面分布特征的水體指數(shù)包括：

根據(jù)湖泊水體在Landsat?MSS或HJ-1A/B影像不同波段上的光譜響應(yīng)特性提取歸一化差異水體指數(shù)；

或者，

根據(jù)湖泊水體在Landsat?TM/ETM+影像不同波段上的光譜響應(yīng)特性提取增強(qiáng)型歸一化差異水體指數(shù)。

4.根據(jù)權(quán)利要求3所述的基于多光譜影像的湖泊水儲量遙感估算方法，其特征在于，所述歸一化差異水體指數(shù)NDWI采用下式計算：

NDWI=ρGreen-ρNIRρGreen+ρNIR]]>

其中，ρ_Green為Landsat?MSS或HJ-1A/B影像上的綠光波段，ρ_NIR為Landsat?MSS或HJ-1A/B影像上的近紅外波段。

5.根據(jù)權(quán)利要求3所述的基于多光譜影像的湖泊水儲量遙感估算方法，其特征在于，所述增強(qiáng)型歸一化差異水體指數(shù)MNDWI采用下式計算：

MNDWI=ρGreen-ρSWIRρGreen+ρSWIR]]>

其中，ρ_Green為Landsat?TM/ETM+影像上的綠光波段，ρ_SWIR為Landsat?TM/ETM+影像上的短波紅外波段。

6.根據(jù)權(quán)利要求1所述的基于多光譜影像的湖泊水儲量遙感估算方法，其特征在于，步驟S2中根據(jù)提取的水體指數(shù)獲取湖泊水體邊界包括：

根據(jù)提取的水體指數(shù)采用改進(jìn)的最大類間方差閾值法提取湖泊水面柵格數(shù)據(jù)；

通過ArcMap9.3軟件中柵格轉(zhuǎn)矢量工具將提取的湖泊水面柵格數(shù)據(jù)轉(zhuǎn)換成湖泊水體邊界矢量數(shù)據(jù)。