发明名称 一种自动生成烟区检测结果的方法
摘要 本发明公开了一种自动生成烟区检测结果的方法,首先利用经预处理后的高分1号卫星宽幅数据,生成用于后续烟区与非烟区检测的7种指数;基于上述的7种指数,当某像元的值满足预设的条件时,将其判定为云或雪或过火像元,并去除该像元;当某像元的值同时满足特定条件,且为非云、雪和非过火像元的值时,将其判定为烟像元,并予以保留;将多个烟像元组成的烟区的像素按栅格转矢量的方式转换为矢量;将矢量化后的烟区与合成影像叠加显示,并输出该烟区中每个烟像元的中心经纬度、所在省市和地物类型信息。该方法能利用卫星数据通过烟的识别来发现林火,从而解决森林火灾预防扑救过程中急需林火卫星早期发现的问题。
申请公布号 CN105096511A 申请公布日期 2015.11.25
申请号 CN201510468951.6 申请日期 2015.08.04
申请人 中国林业科学研究院资源信息研究所 发明人 覃先林;祖笑锋;尹凌宇;孙桂芬
分类号 G08B17/10(2006.01)I 主分类号 G08B17/10(2006.01)I
代理机构 代理人
主权项 一种自动生成烟区检测结果的方法,其特征在于,利用经预处理后的高分1号卫星宽幅数据,生成用于后续烟区与非烟区检测的7种指数;其中,所述7种指数分别为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>SI</mi><mo>=</mo><mfrac><mrow><msub><mi>&rho;</mi><mi>Nir</mi></msub><mo>_</mo><msub><mi>&rho;</mi><mi>blue</mi></msub></mrow><mrow><msub><mi>&rho;</mi><mi>Nir</mi></msub><mo>+</mo><msub><mi>&rho;</mi><mi>blue</mi></msub></mrow></mfrac><mo>,</mo><mi>NDWI</mi><mo>=</mo><mfrac><mrow><msub><mi>&rho;</mi><mi>Green</mi></msub><mo>-</mo><msub><mi>&rho;</mi><mi>R</mi></msub></mrow><mrow><msub><mi>&rho;</mi><mi>Green</mi></msub><mo>+</mo><msub><mi>&rho;</mi><mi>R</mi></msub></mrow></mfrac><mo>,</mo><mi>BAI</mi><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mrow><mo>(</mo><mn>0.1</mn><mo>-</mo><msub><mi>&rho;</mi><mi>R</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mn>0.06</mn><mo>-</mo><msub><mi>&rho;</mi><mi>Nir</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>,</mo><mi>R</mi><mfrac><msub><mi>&rho;</mi><mi>blue</mi></msub><msub><mi>&rho;</mi><mi>R</mi></msub></mfrac><mo>,</mo></mrow>]]></math><img file="FSA0000119732450000011.GIF" wi="1569" he="138" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>GEMI</mi><mo>=</mo><mi>&eta;</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>0.25</mn><mi>&eta;</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><msub><mi>&rho;</mi><mi>R</mi></msub><mo>-</mo><mn>0.125</mn></mrow><mrow><mn>1</mn><mo>-</mo><msub><mi>&rho;</mi><mi>R</mi></msub></mrow></mfrac><mo>,</mo><mi>NDVI</mi><mo>=</mo><mfrac><mrow><msub><mi>&rho;</mi><mi>Nir</mi></msub><mo>_</mo><msub><mi>&rho;</mi><mi>R</mi></msub></mrow><mrow><msub><mi>&rho;</mi><mi>Nir</mi></msub><mo>+</mo><msub><mi>&rho;</mi><mi>R</mi></msub></mrow></mfrac><mo>,</mo><mi>SVI</mi><mo>=</mo><mfrac><mrow><msub><mi>&rho;</mi><mi>Nir</mi></msub><mo>-</mo><msub><mi>&rho;</mi><mi>R</mi></msub></mrow><mrow><msub><mi>&rho;</mi><mi>Nir</mi></msub><mo>+</mo><msub><mi>&rho;</mi><mi>R</mi></msub></mrow></mfrac><mo>&times;</mo><msub><mi>&rho;</mi><mi>Nir</mi></msub><mo>,</mo></mrow>]]></math><img file="FSA0000119732450000012.GIF" wi="1563" he="129" /></maths>上式中,ρ<sub>R</sub>为高分1号卫星宽幅数据第3波段的反射率,ρ<sub>Nir</sub>为高分1号卫星宽幅数据第4波段的反射率,ρ<sub>Green</sub>为高分1号卫星宽幅数据第2波段的反射率,ρ<sub>blue</sub>为高分1号卫星宽幅数据第1波段的反射率;SI为归一化烟指数;NDWI为归一化水体指数;BAI为燃烧指数;GEMI为全球环境植被指数;NDVI为归一化植被指数;SVI为土壤调节植被指数;R为比值指数;且<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&eta;</mi><mo>=</mo><mfrac><mrow><mn>2</mn><mrow><mo>(</mo><msup><msub><mi>&rho;</mi><mi>Nir</mi></msub><mn>2</mn></msup><mo>-</mo><msup><msub><mi>&rho;</mi><mi>R</mi></msub><mn>2</mn></msup><mo>)</mo></mrow><mo>+</mo><mn>1.5</mn><msub><mi>&rho;</mi><mi>Nir</mi></msub><mo>+</mo><mn>0.5</mn><msub><mi>&rho;</mi><mi>R</mi></msub></mrow><mrow><msub><mi>&rho;</mi><mi>Nir</mi></msub><mo>+</mo><msub><mi>&rho;</mi><mi>R</mi></msub><mo>+</mo><mn>0.5</mn></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FSA0000119732450000013.GIF" wi="753" he="161" /></maths>基于上述的7种指数,当某像元的值满足预设的条件时,将其判定为云或雪或过火像元,并去除该像元;当某像元的值同时满足0.17<ρ<sub>blue</sub><0.30,‑0.1<SI<0.02,0.2<GEMI<0.4,‑1.0<NDVI<0.2和1.0<R<1.3的条件,且为非云、雪和非过火像元的值时,将其判定为烟像元,并予以保留;将多个烟像元组成的烟区的像素按栅格转矢量的方式转换为矢量;将矢量化后的烟区与合成影像叠加显示,并输出该烟区中每个烟像元的中心经纬度、所在省市和地物类型信息。
地址 100091 北京市海淀区东小府2号