1.一种无需大气校正的近海浒苔绿潮遥感监测方法,其特征在于,包括如下步骤:s1、收集时空匹配的卫星遥感影像X和Landsat8‑OLI影像对,并分别计算得到所述卫星影像X的大气层顶反射率ρx,i和所述Landsat8‑OLI影像的大气层顶反射率ρOLI,i,再计算Landsat8‑OLI的湿度分量TCWOLI;
s2、基于上述步骤s1中收集的时空匹配的卫星遥感影像X和Landsat8‑OLI影像对,采集若干种目标地物的像元样本,得到所述卫星影像X的大气层顶反射率ρx,i和所述Landsat8‑OLI的湿度分量TCWOLI组成的时空匹配样本数据集;
s3、基于上述步骤s2中收集的时空匹配样本数据集,采用多元线性拟合方法和施密特正交方法推导出面向所述卫星影像X的大气层顶反射率ρx,i的类缨帽变换矩阵c=[c1,i c2,i c3,i],其中,i=1,2,3,4;
s4、构建绿度指数算法;
s5、将上述步骤s4中的绿度指数算法应用于长时间序列卫星遥感影像,计算并得到每一副卫星影像的绿度遥感数据,进而得到卫星影像的浒苔遥感监测产品,最终基于长时序的浒苔遥感监测产品集分析并得到浒苔绿潮爆发过程的分布信息和漂移路径的动态信息。
2.根据权利要求1所述的一种无需大气校正的近海浒苔绿潮遥感监测方法,其特征在于,所述步骤s1中所述卫星影像X的大气层顶反射率ρx,i和所述Landsat8‑OLI影像的大气层顶反射率ρOLI,i的计算公式分别如下:其中,DNx,i表示所述卫星影像X在第i波段处的亮度值,Gaini、Biasi和Esuni分别表示第i波段的定标增益值、定标偏移值和大气顶部太阳辐照度,d为日地天文单位距离,θs为太阳天顶角;
ρOLI,i=(Mi·DNOLI,i+Ai)/cos(θs);
其中,DNOLI,i表示所述Landsat8‑OLI影像在第i波段处的亮度值,Mi和Ai分别为第i波段的反射率乘法调整系数和加法调整系数。
3.根据权利要求1所述的一种无需大气校正的近海浒苔绿潮遥感监测方法,其特征在于,所述步骤s1中的Landsat8‑OLI的湿度分量TCWOLI的计算公式如下:TCWOLI=0.13ρOLI,2+
0.06ρOLI,3+0.03ρOLI,4‑0.61ρOLI,5‑0.55ρOLI,6‑0.55ρOLI,7,其中,ρOLI,i为所述Landsat8‑OLI影像的第i波段大气层顶反射率。
4.根据权利要求3所述的一种无需大气校正的近海浒苔绿潮遥感监测方法,其特征在于,所述步骤s3中的类缨帽变换矩阵c的推导公式具体如下:首先,采用多元线性拟合方法计算湿度分量系数c3,i,计算公式:其中,wi为拟合系数,ρx,i表示所述卫星遥感影像X第i波段的大气层顶反射率;
接着,采用施密特正交方法计算亮度分量系数c1,i,计算公式:其中, 表示裸地的大气层顶反射率光谱平均值,D13和bi为中间变量;
最后,同样采用施密特正交方法计算绿度分量系数c2,i,计算公式:其中, 表示浒苔像元的大气层顶反射率光谱平均值,D12、D32、和gi为中间变量。
5.根据权利要求4所述的一种无需大气校正的近海浒苔绿潮遥感监测方法,其特征在于,所述步骤s4中的绿度指数算法计算公式如下:TCG=c2,1·ρx,1+c2,2·ρx,2+c2,3·ρx,3+c2,4·ρx,4;
其中,TCG为绿度分量。