发明名称 基于高分辨率影像的区域自适应耕地提取方法
摘要 本发明属于遥感应用技术领域,具涉及一种基于高分辨率影像的区域自适应耕地提取方法,包括:对数据进行预处理;对高分辨率影像进行多尺度分割,确定耕地的最优分割尺度,得到耕地分类的最优分割层;在最优分割层中构建多个耕地样本集,形成每个耕地样本集的属性特征空间;进行最优影像特征空间选择;根据最优影像特征空间构建耕地提取的模糊分类器;基于耕地样本集和模糊分类器进行耕地提取。本发明通过耕地样本集与模糊分类器相结合的方式实现了耕地对象的自动识别。通过选择最优分割尺度和不同耕地样本集的属性特征空间,提高了高分辨率影像的分割对象的均质性和模糊分类器的提取精度;通过耕地样本集的构建,降低了数据源对时相和数量的限制要求。
申请公布号 CN105608474A 申请公布日期 2016.05.25
申请号 CN201511034394.3 申请日期 2015.12.30
申请人 二十一世纪空间技术应用股份有限公司 发明人 文强;丁媛;李丽;纪中奎;周会珍;沙漠泉;周淑芳;张强;任昊冬
分类号 G06K9/62(2006.01)I;G06T7/00(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 国防专利服务中心 11043 代理人 文伟能
主权项 一种基于高分辨率影像的区域自适应耕地提取方法,其特征在于,所述提取方法包括:S1,对数据进行预处理;S2,对高分辨率影像进行多尺度分割,确定耕地的最优分割尺度,得到耕地分类的最优分割层;S3,在最优分割层中构建多个耕地样本集,形成每个耕地样本集的属性特征空间;S4,对每一个耕地样本集,进行最优影像特征空间选择;S5,对每一个耕地样本集,根据最优影像特征空间构建耕地提取的模糊分类器;S6,基于耕地样本集和模糊分类器进行耕地提取:S61,根据耕地样本集的属性特征空间,逐一对步骤S2最优分割后的对象进行属性特征计算;S62,分别对S61得到的属性特征计算结果与步骤S3构建的耕地样本集的属性特征空间范围值进行匹配,如果对象的属性特征值在某一耕地样本集的属性特征范围内,则记为1,否则记为0,得到匹配矩阵D<sub>ij</sub><maths num="0001"><math><![CDATA[<mrow><msub><mi>D</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><msup><mfenced open = '[' close = ']'><mtable><mtr><mtd><mrow><msub><mi>l</mi><mn>1</mn></msub><mo>,</mo><msub><mi>l</mi><mn>2</mn></msub><mo>,</mo><msub><mi>l</mi><mn>3</mn></msub><mo>,</mo><msub><mi>l</mi><mn>4</mn></msub><mo>,</mo><mo>...</mo><msub><mi>l</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>p</mi><mn>1</mn></msub><mo>,</mo><msub><mi>p</mi><mn>2</mn></msub><mo>,</mo><msub><mi>p</mi><mn>3</mn></msub><mo>,</mo><msub><mi>p</mi><mn>4</mn></msub><mo>,</mo><mo>...</mo><msub><mi>p</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>,</mo><msub><mi>n</mi><mn>2</mn></msub><mo>,</mo><msub><mi>n</mi><mn>3</mn></msub><mo>,</mo><msub><mi>n</mi><mn>4</mn></msub><mo>,</mo><mo>...</mo><msub><mi>n</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>b</mi><mn>1</mn></msub><mo>,</mo><msub><mi>b</mi><mn>2</mn></msub><mo>,</mo><msub><mi>b</mi><mn>3</mn></msub><mo>,</mo><msub><mi>b</mi><mn>4</mn></msub><mo>,</mo><mo>...</mo><msub><mi>b</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mn>...................</mn></mtd></mtr><mtr><mtd><mrow><msub><mi>w</mi><mn>1</mn></msub><mo>,</mo><msub><mi>w</mi><mn>2</mn></msub><mo>,</mo><msub><mi>w</mi><mn>3</mn></msub><mo>,</mo><msub><mi>w</mi><mn>4</mn></msub><mo>,</mo><mo>...</mo><msub><mi>w</mi><mi>i</mi></msub></mrow></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>]]></math><img file="FSA0000125645100000011.GIF" wi="594" he="326" /></maths>其中,l,p,n,b,…,w分别为耕地样本集的属性特征,i为耕地样本集的个数,j为属性特征的个数,T表示矩阵的转置;S63,根据定义的耕地样本集中每个属性特征所占权重值与匹配矩阵D<sub>ij</sub>,分别计算对象与每一个耕地样本集的相似度,得到相似度矩阵S<sub>i</sub>,<maths num="0002"><math><![CDATA[<mrow><msub><mi>S</mi><mi>i</mi></msub><mo>=</mo><msub><mi>D</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&times;</mo><msub><mi>&delta;</mi><mi>j</mi></msub><mo>=</mo><msup><mfenced open = '[' close = ']'><mtable><mtr><mtd><mrow><msub><mi>l</mi><mn>1</mn></msub><mo>,</mo><msub><mi>l</mi><mn>2</mn></msub><mo>,</mo><msub><mi>l</mi><mn>3</mn></msub><mo>,</mo><msub><mi>l</mi><mn>4</mn></msub><mo>,</mo><mo>...</mo><msub><mi>l</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>p</mi><mn>1</mn></msub><mo>,</mo><msub><mi>p</mi><mn>2</mn></msub><mo>,</mo><msub><mi>p</mi><mn>3</mn></msub><mo>,</mo><msub><mi>p</mi><mn>4</mn></msub><mo>,</mo><mo>...</mo><msub><mi>p</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>,</mo><msub><mi>n</mi><mn>2</mn></msub><mo>,</mo><msub><mi>n</mi><mn>3</mn></msub><mo>,</mo><msub><mi>n</mi><mn>4</mn></msub><mo>,</mo><mo>...</mo><msub><mi>n</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>b</mi><mn>1</mn></msub><mo>,</mo><msub><mi>b</mi><mn>2</mn></msub><mo>,</mo><msub><mi>b</mi><mn>3</mn></msub><mo>,</mo><msub><mi>b</mi><mn>4</mn></msub><mo>,</mo><mo>...</mo><msub><mi>b</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>...</mo><mo>...</mo><mo>...</mo><mo>...</mo><mo>...</mo><mo>...</mo><mo>&CenterDot;</mo></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>w</mi><mn>1</mn></msub><mo>,</mo><msub><mi>w</mi><mn>2</mn></msub><mo>,</mo><msub><mi>w</mi><mn>3</mn></msub><mo>,</mo><msub><mi>w</mi><mn>4</mn></msub><mo>,</mo><mo>...</mo><msub><mi>w</mi><mi>i</mi></msub></mrow></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>&times;</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>&delta;</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>&delta;</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>&delta;</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><msub><mi>&delta;</mi><mn>4</mn></msub></mtd></mtr><mtr><mtd><mo>...</mo></mtd></mtr><mtr><mtd><msub><mi>&delta;</mi><mi>j</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>S</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>S</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>S</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><msub><mi>S</mi><mn>4</mn></msub></mtd></mtr><mtr><mtd><mo>...</mo></mtd></mtr><mtr><mtd><msub><mi>S</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FSA0000125645100000012.GIF" wi="1104" he="335" /></maths>其中,δ<sub>i</sub>为样本集属性特征权重矩阵,δ<sub>i</sub>根据区域情况进行自定义;S64,设定相似度矩阵S<sub>i</sub>的阈值判断条件,若对象与某个耕地样本集的相似度满足阈值判断条件,则对象与某个耕地样本集相匹配;S65,根据S64的匹配结果,对该对象执行相匹配的耕地样本集对应的模糊分类器,得到分类类别函数E(i)和类别隶属度ω(i),<img file="FSA0000125645100000021.GIF" wi="977" he="136" />S66,采用加权投票法进行耕地加权计算,计算公式如下:<maths num="0003"><math><![CDATA[<mrow><mi>E</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><mi>E</mi><mo>(</mo><mi>i</mi><mo>)</mo><mo>&times;</mo><mi>&omega;</mi><mo>(</mo><mi>i</mi><mo>)</mo><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000125645100000022.GIF" wi="514" he="163" /></maths>其中,E(x)为耕地加权判断函数,n为与对象x相匹配的模糊分类器总数,i为模糊分类器的个数,i=1,…,n;S67,对对象是否为耕地进行判断,判断条件如下:<img file="FSA0000125645100000023.GIF" wi="995" he="270" />
地址 100096 北京市海淀区西三旗建材城东路26号