摘 "要： 以農(nóng)業(yè)用地與人工建筑的分割為研究目的，提出一種可在衛(wèi)星或航拍遙感圖像中自動(dòng)分割農(nóng)用田地與人工建筑的圖像處理算法。該算法使用圖像區(qū)域合成、濾波、灰度統(tǒng)計(jì)、背景分割等處理技術(shù)，可實(shí)現(xiàn)大面積農(nóng)用田地中人工建筑及道路的分割，并給出面積的估算。在此給出了解決思路，并給出試驗(yàn)效果。該技術(shù)對(duì)提高農(nóng)地統(tǒng)計(jì)自動(dòng)化水平、農(nóng)業(yè)用地侵占普查效率有一定工程意義。

關(guān)鍵詞： 耕地; 人工建筑; 閾值; 圖像分割

中圖分類(lèi)號(hào)： TN964?34 " " " " " " " " " 文獻(xiàn)標(biāo)識(shí)碼： A " " " " " " " " " " " " " "文章編號(hào)： 1004?373X（2015）06?0105?03

Research on automatic segmentation of agricultural land and artificial buildings in remote sensing image

WANG Min?yan1， ZHAO Kun2

（1. Commission for Discipline Inspection and Supervision Bureau of Jining， Jining 272045， China;

2. The 27th Research Institute of China Electronics Technology Group Corporation， Zhengzhou 450047， China）

Abstract： An image processing algorithm which can automatically divide up agricultural fields and artificial buildings in remote sensing images from satellite or aerial photo is proposed. With the technologies of image area synthesis， filtering， histogram statistics and background segmentation， the artificial buildings and roads in large?scale agricultural lands can be segmented by the algorithm， and the area calculation can be realized. Some solutions and test results are given in this paper. It is significant for this algorithm to improve the automation statistic level of agricultural land.

Keywords： agricultural field; artificial building; threshold; image segmentation

0 "引 "言

作為人口大國(guó)，農(nóng)林業(yè)對(duì)我國(guó)家安全和社會(huì)穩(wěn)定有重要意義，隨著經(jīng)濟(jì)的高速發(fā)展，鄉(xiāng)鎮(zhèn)及城市建設(shè)對(duì)土地需求日趨增長(zhǎng)，合理有效地使用土地資源，監(jiān)控土地實(shí)際使用情況是國(guó)家研究的重要課題。隨著現(xiàn)代航天技術(shù)、航空技術(shù)、探測(cè)器技術(shù)高速發(fā)展，國(guó)家投入大量衛(wèi)星、飛機(jī)及無(wú)人機(jī)進(jìn)行耕地普查，新技術(shù)的廣泛應(yīng)用給農(nóng)業(yè)用地統(tǒng)計(jì)、普查、災(zāi)害預(yù)警帶來(lái)革命性發(fā)展，提供了豐富的數(shù)據(jù)資源。但社會(huì)各界對(duì)農(nóng)業(yè)統(tǒng)計(jì)信息的需求越來(lái)越大，對(duì)統(tǒng)計(jì)質(zhì)量要求越來(lái)越高，按照傳統(tǒng)的方法，僅僅依靠手工作業(yè)是遠(yuǎn)遠(yuǎn)不能適應(yīng)實(shí)際海量數(shù)據(jù)篩選需要的。本文提出一種可在衛(wèi)星或航拍遙感圖像中自動(dòng)分割農(nóng)用田地與人工建筑的圖像處理算法，該算法使用圖像區(qū)域合成、濾波、灰度統(tǒng)計(jì)、背景分割等處理技術(shù)，可實(shí)現(xiàn)農(nóng)用田地中疑似建筑及道路進(jìn)行分割，并給出面積的估算。為后續(xù)土地和人工建筑確認(rèn)標(biāo)注提供預(yù)先提示。

1 "遙感圖像預(yù)處理

遙感圖像在形成、傳輸、記錄過(guò)程中，成像系統(tǒng)、傳輸介質(zhì)和記錄設(shè)備都會(huì)引起圖像質(zhì)量的下降[1?2]，典型表現(xiàn)為圖像模糊、失真、有噪聲等。圖像預(yù)處理的目的就是盡可能恢復(fù)退化圖像的本來(lái)面目。引起圖像退化的因素很多，圖像預(yù)處理將根據(jù)退化的原因，建立相應(yīng)的數(shù)學(xué)模型，并沿著圖像降質(zhì)的逆過(guò)程恢復(fù)原圖像。

農(nóng)用土地遙感圖像所關(guān)注的信息大部分集中在低頻和中頻段，為了有效抑制圖像的噪聲并盡可能多的保留有用的信息，采用均值濾波、中值濾波、低通濾波、匹配濾波等多種濾波器相結(jié)合，綜合線形濾波器和非線形濾波器的優(yōu)點(diǎn)，研究具有實(shí)時(shí)特性的濾波網(wǎng)絡(luò)，最大程度濾除噪聲，提高圖像的信噪比，減少圖像的失真[3?4]。

原始圖像中各點(diǎn)的灰度值f（x，y）用該點(diǎn)鄰域的灰度平均值代替，以抑制具有突變性質(zhì)的噪聲點(diǎn)。設(shè)g（x，y）為處理后的灰度值，則該方法可表示為：

[g（x，y）=1MN（m，n）∈S（x，y）f（m，n）] " " （1）

式中：[S（x，y）]為[（x，y）]的鄰域;m，n為[S（x，y）]中的像素?cái)?shù)。

圖像的預(yù)處理在空域中進(jìn)行，也可在頻域中進(jìn)行，將一維信號(hào)低通濾波器在二維圖像中應(yīng)用，圖像經(jīng)過(guò)二維傅里葉變換后，噪聲頻譜位于空間頻率較高的區(qū)域，而圖像中關(guān)鍵信息處于空間頻率較低的區(qū)域，從而實(shí)現(xiàn)對(duì)高頻分量的預(yù)制[5]。

轉(zhuǎn)移函數(shù)為：

[H（x，y）=11+[D（x，y）D0]2n] （2）

式中：[D0]為截止頻率;[n]為階數(shù)。

2 "遙感圖像分割與處理

嚴(yán)格意義上的圖像分割是指每一個(gè)物體都能和它的背景區(qū)分開(kāi)來(lái)而成為一個(gè)獨(dú)立完整的區(qū)域，人的視覺(jué)就包括這個(gè)過(guò)程。但是，在實(shí)際應(yīng)用中這是很困難的，因?yàn)閳D像分割沒(méi)有統(tǒng)一的判別標(biāo)準(zhǔn)，也沒(méi)有通用的分割方法，一般來(lái)說(shuō)，只要能把感興趣的部分從背景中分割出來(lái)就可以接受[6]。對(duì)遙感圖像進(jìn)行分割就是把圖像中不同于農(nóng)業(yè)用地的人工建筑或道路等具有特殊意義的區(qū)域分割開(kāi)來(lái)，這些區(qū)域本身并不相關(guān)，但每一個(gè)區(qū)域都不同于農(nóng)業(yè)用地[7?8]。

2.1 "大面積耕地直方圖分析

由于大面積農(nóng)用耕地的主背景的分布相對(duì)均勻，從圖像的直方圖分析明顯存在單峰（背景），如果有道路或人工建筑，則直方圖中將存在背景峰之外的多個(gè)。

lt;E：＼王芳＼現(xiàn)代電子技術(shù)201506＼現(xiàn)代電子技術(shù)15年38卷第6期＼Image＼05t1.tifgt;

圖1 大面積農(nóng)用田地圖及直方圖

2.2 "基于大面積耕地背景的分割閾值提取

[H（x）]為圖像直方圖分布，則圖像的平均灰度為[9?10]：

[hx=x·hxhx] （3）

圖像的灰度分布誤差為：

[σ=h（x）-h（x）?x2h（x）] （4）

將[h（x）]-ks作為圖像分割的第一個(gè)閾值，[h（x）]+ks作為圖像分割的第二個(gè)閾值。將灰度在[[h（x）]-ks，[h（x）]+ks]內(nèi)的像素判為背景，其他灰度圖像像素判為目標(biāo)，則就可以將目標(biāo)從背景中“分割”出來(lái)。如圖2，圖3為原始圖像及經(jīng)過(guò)預(yù)處理后的圖像，應(yīng)用式（4）求解出[σ]及分割門(mén)限，經(jīng)二值化后分割圖像如圖4所示。

lt;E：＼王芳＼現(xiàn)代電子技術(shù)201506＼現(xiàn)代電子技術(shù)15年38卷第6期＼Image＼05t2.tifgt;

圖2 原始圖像

lt;E：＼王芳＼現(xiàn)代電子技術(shù)201506＼現(xiàn)代電子技術(shù)15年38卷第6期＼Image＼05t3.tifgt;

圖3 圖像預(yù)處理

lt;E：＼王芳＼現(xiàn)代電子技術(shù)201506＼現(xiàn)代電子技術(shù)15年38卷第6期＼Image＼05t4.tifgt;

圖4 圖像分割

根據(jù)分割出的圖像可對(duì)目標(biāo)像素面積的統(tǒng)計(jì)，采用投影方法，可獲得圖像不同列中非農(nóng)用地的像素?cái)?shù)，如圖5所示，通過(guò)累加計(jì)算可獲得設(shè)定區(qū)域非農(nóng)用地占用比例，從而獲得實(shí)際農(nóng)用地面積。

2.3 "應(yīng)用

通過(guò)對(duì)遙感圖像處理，可對(duì)耕地或建筑面積進(jìn)行估算、統(tǒng)計(jì)，對(duì)疑似建筑物、道路進(jìn)行標(biāo)注，提高農(nóng)業(yè)用地統(tǒng)計(jì)自動(dòng)化水平，并提高土地管理部門(mén)、農(nóng)業(yè)部門(mén)、農(nóng)業(yè)統(tǒng)計(jì)部門(mén)及其他有關(guān)部門(mén)在技術(shù)手段和行政管理方面的優(yōu)勢(shì)，采用先進(jìn)的遙感技術(shù)、全面統(tǒng)計(jì)和抽樣調(diào)查等各種方法，對(duì)耕地及其他農(nóng)業(yè)用地的數(shù)量進(jìn)行全面的監(jiān)測(cè)和管理，及時(shí)對(duì)不利于耕地和其他農(nóng)業(yè)用地保護(hù)及使用的情況進(jìn)行處理和糾正。

lt;E：＼王芳＼現(xiàn)代電子技術(shù)201506＼現(xiàn)代電子技術(shù)15年38卷第6期＼Image＼05t5.tifgt;

圖5 目標(biāo)面積統(tǒng)計(jì)

3 "結(jié) "語(yǔ)

本文根據(jù)大面積農(nóng)用地統(tǒng)計(jì)需求，結(jié)合其圖像特性，提出一種大面積耕地背景的耕地及人工建筑分割方法，可以有效提高耕地統(tǒng)計(jì)自動(dòng)化水平，并輔助監(jiān)測(cè)農(nóng)用耕地占用情況。

參考文獻(xiàn)

[1] 章孝燦，黃智才，趙元洪.遙感數(shù)字圖像處理[M].杭州：浙江大學(xué)出版社，1997.

[2] 何紅艷，楊居奎，齊文雯.大氣對(duì)遙感衛(wèi)星圖像品質(zhì)的影響分析[J].航天返回與遙感，2011（2）：42?47.

[3] 汪哲慎，李翠華.基于目標(biāo)特征的圖像預(yù)處理技術(shù)[J].廈門(mén)大學(xué)學(xué)報(bào)：自然科學(xué)版，2004（4）：473?478.

[4] 沃焱，韓國(guó)強(qiáng).基于自適應(yīng)預(yù)處理的圖像分割方法[J].電子與信息學(xué)報(bào)，2007（1）：87?91.

[5] 楊枝開(kāi)，王開(kāi).數(shù)字圖像獲取處理及實(shí)踐應(yīng)用[M].北京：人民郵電出版社，2003.

[6] 章毓晉.圖像分割[M].北京：科學(xué)出版社，2001.

[7] 殷瑞娟， 施潤(rùn)和.一種高分辨率遙感影像的最優(yōu)分割尺度自動(dòng)選取方法[J].地球信息科學(xué)學(xué)報(bào)，2013（6）：904?910.

[8] 張彥.遙感圖像中建筑物分割方法研究[D].天津：河北工業(yè)大學(xué)，2012.

[9] 韓思奇，王蕾.圖像分割的閾值法綜述[J].系統(tǒng)工程與電子技術(shù)，2002（6）：91?94.

[10] 杜曉晨，劉建平.改進(jìn)的模糊閾值圖像分割方法[J].光電工程，2005（10）：51?53.