Copernicus Browser数据来源,欧氏空间局需要注册登录,有点麻烦可以搜索教程多试几次。
实验数据:
2023年4月11日云南省玉溪市发生火灾,本次实验提取该处火灾面积。
第一步:
1、加载进两个文件夹下的“MTD_MSIL2A.xml”文件。之后修改RGB波段为843(近红、红、绿),通过目视解译就可以看出两个受灾地区,本次实验仅作教学,所以只裁取一个火灾区域做讲解。
2、裁剪研究区域,菜单栏有绘制感兴趣区域的快捷键,点击快捷键进行绘制火灾区域。
3、绘制好以后,File→save as→ENVI ,之后出现Data Selection窗口,源文件包括11号和21号,所以裁剪区域同样要裁剪两次不同日期的底图,波段选择4,mask选择上一步绘制的感兴趣区域。
左图21号,右图11号
第二步:
1、工具栏选择change Detection→Image Change Workflow(图像变化检测工作流),出现Save ROI选择否,之后对话框中选择前后时相图,直接next,下一步也直接next。
2、工具提供两种方法image difference(图像差值法) 和image transform(图像变换法),选择差值法。该方法提供三种方法去检测,每种方法出图也都不一样,可以都试试看效果,如果不懂可以选择左下角的蓝色问号,官方网站会介绍每一种方法,很详细的。本实验选择第三个波谱角。.
3、直接变化检测工作流工具能从结果中提取三种变化信息:Increase and Decrease:增加 (蓝色) 和减少(红色)变化信息;Increase Only:增加 (蓝色)变化信息;Decrease Only:减少(红色)变化信息。实验选择减少(红色)变化信息
自动设置阈值有四种可选方法:
Otsu's: 基于直方图形状的方法。
Tsai's: 基于力矩的方法。
Kapur's: 基于信息熵的方法。
Kittler's: 基于直方图形状的方法。
4、
下一步输出,Additional Export勾选Export Change Class Statistics(输出统计文本文件),Finish完成输出。
每一种方法提取出来的都不一样,如果对结果不满意可以叠加不同方法所提取的结果,选择最适合的方式。