传感器不变大气校正(siac)
SIAC的Python项目详细描述
传感器不变大气校正(SIAC)
冯寅
加州大学洛杉矶分校地理系
ucfafyi@ucl.ac.uk
该大气校正方法利用modis mcd43 brdf产品对地表进行粗分辨率模拟。建立了基于modis-psf的modis与哨兵2/landsat 8比例尺差异模型。以ecmwf-cams预报作为大气状态预报的先验,结合6s模式求解大气参数。我们没有地形校正,使用均匀曲面而不考虑brdf效应。
所需数据:
- MCD43:Sentinel 2/Landsat 8号卫星探测日期前16天和后16天
- ecmwf-cams近实时预测:时间步长为3小时,开始时间为00:00:00,从2015年4月1日开始的数据在ucl服务器中镜像到:http://www2.geog.ucl.ac.uk/~ucfafyi/cams/
- 全局dem:根据astgtm2 dem构建的全局dem vrt文件,陆地上的大部分dem都镜像到ucl服务器:http://www2.geog.ucl.ac.uk/~ucfafyi/eles/
- 模拟器:大气路径反射、总透射和单次散射反照率模拟器,以及使用6s.v2训练的哨兵2号、陆地卫星8号和modis模拟器,可在:http://www2.geog.ucl.ac.uk/~ucfafyi/emus/
安装:
- 直接从Github
pip install https://github.com/MarcYin/SIAC/archive/master.zip
- 使用pypi
pip install SIAC
- 使用水蟒
conda install -c f0xy -c conda-forge siac
为了节省安装gdal的时间:
conda install -c conda-forge gdal>2.1
哨兵2号和陆地卫星8号的典型用法:
fromSIACimportSIAC_S2SIAC_S2('/directory/where/you/store/S2/data/')# this can be either from AWS or Senitinel offical package
fromSIACimportSIAC_L8SIAC_L8('/directory/where/you/store/L8/data/')
为了更详细地说明使用情况,陆地卫星5号的校正示例如下所示here
示例和地图:
apage显示一些校正示例。
amap用于比较toa和boa。
引文:
Yin,F.,Lewis,P.E.,Gomez Dans,J.,&Wu,Q.(2019年2月21日)。一种传感器不变大气校正方法:在哨兵-2/msi和陆地卫星8/oli上的应用。https://doi.org/10.31223/osf.io/ps957
许可证
GNU通用公共许可v3