劉 研，向 陽，李桂盛，雷 淦，聶 瑤，向 慧

（湖南農(nóng)業(yè)大學機電工程學院，湖南 長沙 410128）

GPS技術已深入國民經(jīng)濟的各個部門，在農(nóng)業(yè)中應用的日益廣泛[1]。相對傳統(tǒng)的人工測量面積方法，基于GPS數(shù)據(jù)計算面積可有效提高計算效率，降低作業(yè)成本[2]，因而應用GPS對農(nóng)機作業(yè)面積進行測量，成為獲取農(nóng)機作業(yè)面積和作業(yè)量的重要手段。

近年來，農(nóng)業(yè)科技工作者在農(nóng)田面積測量方面不斷研究，并且進行大量實驗，開發(fā)出一系列具有可操作性的系統(tǒng)及產(chǎn)品，為農(nóng)田面積測量提供了理論依據(jù)及技術支持。何勇、趙春江在《精細農(nóng)業(yè)》一書中，系統(tǒng)地論述了3S技術在農(nóng)業(yè)中的應用，詳述了農(nóng)田面積測繪的相關內(nèi)容[3]。夏友福應用CASIO fx-4500P編程器及手持GPS測量儀計算作業(yè)區(qū)域面積[4]，但此方法不適宜狹長地帶林地測量。熊中剛深入研究了插秧機智能作業(yè)遠程監(jiān)控系統(tǒng)，詳述了農(nóng)機作業(yè)軌跡擬合、作業(yè)區(qū)域及其邊界點識別、作業(yè)區(qū)域面積計算算法，但該方法處理數(shù)據(jù)量太大，在計算復雜不規(guī)則作業(yè)區(qū)域時誤差較大[5]。

文章基于高斯面積公式提出一種利用GPS數(shù)據(jù)進行農(nóng)機作業(yè)面積測量的算法，可有效實現(xiàn)農(nóng)機作業(yè)面積的測量，并具有算法簡單，易于實現(xiàn)，對硬件要求低的特點。

1 GPS數(shù)據(jù)獲取

GPS數(shù)據(jù)采集裝置主要由單片機系統(tǒng)、移動電源、GPS模塊、U盤讀寫模塊和U盤組成，系統(tǒng)結(jié)構如圖1所示。

圖1 GPS數(shù)據(jù)采集裝置系統(tǒng)結(jié)構

單片機選用Arduino MEGA2560系統(tǒng)，該系統(tǒng)核心處理器為ATmega2560，具有54路數(shù)字輸入/輸出接口和4路串行通信接口，片內(nèi)包含256KB Flash、8KB SRAM和4KB EEPROM，運行時鐘頻率為 16MHz。GPS模塊選用低成本、高性能的UBLOX M8N，該模塊具有22個跟蹤通道，支持GPS與北斗雙星系統(tǒng)，支持WAAS、EGNOS、MSAS和GAGAN廣域差分，典型定位精度為3m，最大更新速率為10Hz，其通訊協(xié)議支持標準的NMEA-0183協(xié)議。U盤讀寫模塊選取UV706，該模塊是一款無協(xié)議的U盤數(shù)據(jù)存儲模塊，整合了U盤規(guī)范和FAT32文件格式規(guī)范，采用ARM處理器，U盤接口類型為USB3.0，最大支持32GB的U盤。通過串口發(fā)送數(shù)據(jù)給模塊就能將數(shù)據(jù)存儲在U盤中的文件中，每次上電自動新建一個文件，文件名以D+7位數(shù)字自動命名，文件掉電不丟失，文件大小不超過10Mb。U盤選用金士頓32GU盤，移動電源額定容量為6000mAh，輸出電壓為5V。所有模塊均采用TTL串口通信，將UBLOX M8N連接至單片機的串口1，將U盤讀寫模塊連接至單片機的串口2，設置所有串口通信速率為115200bps。

采用UBLOX專用設置軟件將UBLOX M8N模塊的輸出頻率設置為1Hz，輸出內(nèi)容設置為僅輸出$GPRMC（Recommended Minimum Specific GPS/TRANSIT Data）數(shù)據(jù)，其數(shù)據(jù)格式為“$GPRMC，＜UTC 時間＞，＜定位狀態(tài)＞，＜緯度＞，＜緯度半球＞，＜經(jīng)度＞，＜經(jīng)度半球＞，＜地面速率＞，＜地面航向＞，＜UTC日期＞，＜磁偏角＞，＜磁偏角方向＞，＜模式指示＞*hh＜CR＞＜LF＞”。單片機系統(tǒng)接收到GPS模塊發(fā)送過來的$GPRMC數(shù)據(jù)后，對數(shù)據(jù)進行解析，提取數(shù)據(jù)中的經(jīng)緯度，通過U盤讀寫模塊將經(jīng)緯度數(shù)據(jù)保存在U盤中的TXT文件中。

采集GPS數(shù)據(jù)時，手持GPS數(shù)據(jù)采集裝置沿農(nóng)機作業(yè)區(qū)域邊界繞行一周，即可獲取農(nóng)機作業(yè)區(qū)域輪廓離散點的GPS位置數(shù)據(jù)。

2 作業(yè)面積計算算法研究

2.1 經(jīng)緯度轉(zhuǎn)換為平面直角坐標

由于采集到的GPS數(shù)據(jù)為經(jīng)緯度，無法直接用來計算面積，因此需要對經(jīng)緯度進行坐標轉(zhuǎn)換。文章采用高斯投影正算將經(jīng)緯度轉(zhuǎn)換到平面直角坐標系中進行相關計算[6]。

設中央子午線為 a（115.5°），緯度為 w，經(jīng)度為j，其單位均為度。將緯度轉(zhuǎn)換為平面坐標x，經(jīng)度轉(zhuǎn)換為平面坐標y，其單位均為米。則：

2.2 基于高斯面積公式的面積計算

應用GPS數(shù)據(jù)采集裝置沿作業(yè)區(qū)域邊界繞行一周，即可形成一個按時間順序排列的封閉多邊形。如圖2所示為一個6點序列形成的多邊形區(qū)域。

圖2 封閉多邊形

按傳統(tǒng)面積計算方法難以求解不規(guī)則多邊形的面積，文章選用高斯面積公式計算該多邊形面積，該算法的基本原理是在笛卡爾坐標中的平面中，對于多邊形的頂點，按順序交叉倍增相應的坐標，找到包含多邊形的區(qū)域，并從周圍的多邊形中減去它，以找到其中的多邊形區(qū)域。因為構成多邊形的坐標不斷交叉倍增，就像綁鞋帶一樣，因此該方法也被稱為鞋帶方法。其基本公式描述為：由點序 P1(x1,y1),P2(x2,y2),P3(x3,y3),……Pn(xn,yn)組成的多邊形，其面積S按以下公式計算：

相應的，圖2所示的多邊形面積就等于：

采用C語編制相應的面積計算函數(shù)為：

3 作業(yè)面積計量結(jié)果檢驗

為了驗證面積計量算法的精度，在湖南農(nóng)業(yè)大學東田徑場進行作業(yè)面積計量試驗。

試驗步驟如下：①在湖南農(nóng)業(yè)大學東田徑場的足球場以地面劃線為基準選取3個矩形，分別用皮尺測量各矩形的邊長并計算面積6；②手持GPS數(shù)據(jù)采集裝置分別繞行3個矩形邊界一周，得到邊界輪廓離散點序列的經(jīng)緯度并保存在U盤中；③將3個矩形邊界點的經(jīng)緯度轉(zhuǎn)化為平面坐標；④將平面坐標代入多邊形面積計算程序，得到矩形面積。

試驗結(jié)果如表1所示。其中，相對誤差=（測量面積-計算面積）/測量面積*100%。

表1 試驗結(jié)果

4 結(jié)論

提出一種基于高斯面積公式和GPS數(shù)據(jù)的農(nóng)機作業(yè)面積計量算法，可依據(jù)農(nóng)機作業(yè)區(qū)域邊界點序列計算出農(nóng)機作業(yè)區(qū)域面積，具有算法簡單，實現(xiàn)容易的優(yōu)點。

試驗數(shù)據(jù)表明，當農(nóng)機作業(yè)區(qū)域測量面積為2720㎡、3570㎡、7140㎡時，計量相對誤差分別為2.28%、1.82%、1.09%，基本滿足實際應用要求，且誤差隨著作業(yè)區(qū)域面積增大而減小。