||
Landsat8 OLI的逐像元成像几何可以由官方提供的python包l8_angles来实现。
操作前的几个要点:
1,使用了Pillow包来把生成的数据存储为tif文件,所以请确保提前安装(conda install Pillow);2,l8_angles包不支持Python3.7,所以请使用3.6版。
综上,建议在anaconda中自己设置一个专门的虚拟环境。
详细步骤如下:
1,下载(https://github.com/DHI-GRAS/l8-angles);
2,安装
conda config --add channels DHI-GRAS
conda install py-l8angles
3,运行之,如:
python
import l8angles
from PIL import Image
data = l8angles.calculate_angles('E:/l8anglesmaster/test_ANG.txt', angle_type='SOLAR',subsample=1, bands=[5])
tt=data['sun_az']
img_t=tt[0]
im=Image.fromarray(img_t)
im.save('sun_az.tif')
注,以上得到的tif 没有坐标,可以将原始Landsat影像坐标赋给上述结果。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-24 13:37
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社