JP2795454B2 - カラー情報による境界線の検出方法 - Google Patents
カラー情報による境界線の検出方法Info
- Publication number
- JP2795454B2 JP2795454B2 JP1065153A JP6515389A JP2795454B2 JP 2795454 B2 JP2795454 B2 JP 2795454B2 JP 1065153 A JP1065153 A JP 1065153A JP 6515389 A JP6515389 A JP 6515389A JP 2795454 B2 JP2795454 B2 JP 2795454B2
- Authority
- JP
- Japan
- Prior art keywords
- value
- threshold value
- color
- boundary line
- hough
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
- 238000001514 detection method Methods 0.000 title description 2
- 238000000034 method Methods 0.000 claims description 22
- 238000003384 imaging method Methods 0.000 claims description 16
- 238000005070 sampling Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 12
- 238000011156 evaluation Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 description 7
- 241000209094 Oryza Species 0.000 description 6
- 235000007164 Oryza sativa Nutrition 0.000 description 6
- 235000013339 cereals Nutrition 0.000 description 6
- 235000009566 rice Nutrition 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 4
- 239000003086 colorant Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 208000019300 CLIPPERS Diseases 0.000 description 1
- 244000068988 Glycine max Species 0.000 description 1
- 235000010469 Glycine max Nutrition 0.000 description 1
- 241000209140 Triticum Species 0.000 description 1
- 235000021307 Triticum Nutrition 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 208000021930 chronic lymphocytic inflammation with pontine perivascular enhancement responsive to steroids Diseases 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000003306 harvesting Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000009333 weeding Methods 0.000 description 1
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Image Processing (AREA)
- Closed-Circuit Television Systems (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Image Analysis (AREA)
- Guiding Agricultural Machines (AREA)
Description
【発明の詳細な説明】 〔産業上の利用分野〕 本発明は、圃場に植立された稲、麦、大豆等の当該列
に沿って走行機体を走行させて刈取り脱穀する等のいわ
ゆる汎用コンバインまたは普通形コンバイン等の収穫機
や、圃場を耕作する乗用型耕耘機等の操作機体を、未作
業地と既作業地との境界線に沿って走行させるために、
カラー撮像して得られたカラー画像データにて前記境界
線を検出する2値化処理方法に関するものである。
に沿って走行機体を走行させて刈取り脱穀する等のいわ
ゆる汎用コンバインまたは普通形コンバイン等の収穫機
や、圃場を耕作する乗用型耕耘機等の操作機体を、未作
業地と既作業地との境界線に沿って走行させるために、
カラー撮像して得られたカラー画像データにて前記境界
線を検出する2値化処理方法に関するものである。
圃場に植立した稲列に沿ってコンバインを前進走行さ
せ、当該稲列を複数条づつ刈取るように、当該コンバイ
ンを自動操舵走行させる際して、未刈取地である未作業
地と既刈取地である既作業地との境界線を検出する技術
として、特開昭62−122507号公報では、走行機体に搭載
した白黒ビデオカメラにて既刈取地と未刈取地との両者
に跨る圃場面を撮像し、この画像データをその平均明度
差に基づき2値化処理することにより、前記境界線を抽
出することを提案している。
せ、当該稲列を複数条づつ刈取るように、当該コンバイ
ンを自動操舵走行させる際して、未刈取地である未作業
地と既刈取地である既作業地との境界線を検出する技術
として、特開昭62−122507号公報では、走行機体に搭載
した白黒ビデオカメラにて既刈取地と未刈取地との両者
に跨る圃場面を撮像し、この画像データをその平均明度
差に基づき2値化処理することにより、前記境界線を抽
出することを提案している。
前記先行技術の2値化処理に際して、画像データにお
ける明度変化の微分値が一定以上の値をしきい値とする
処理を実行しているが、このようにしきい値を予め固定
的に設定すると、2値化の処理自体は簡単になるが、画
像データは白黒およびその中間階調(灰色の濃淡)であ
るため、情報量が少なく2値化処理に際して、境界線の
抽出等の特徴抽出に誤差が出るという問題があった。
ける明度変化の微分値が一定以上の値をしきい値とする
処理を実行しているが、このようにしきい値を予め固定
的に設定すると、2値化の処理自体は簡単になるが、画
像データは白黒およびその中間階調(灰色の濃淡)であ
るため、情報量が少なく2値化処理に際して、境界線の
抽出等の特徴抽出に誤差が出るという問題があった。
つまり、撮像したときの天候、稲の品種、稲の生育状
況が種々異なること等撮像条件が変動するから、未刈取
地と既刈取地との白黒および灰色階調に差異が表れ難い
という問題があった。
況が種々異なること等撮像条件が変動するから、未刈取
地と既刈取地との白黒および灰色階調に差異が表れ難い
という問題があった。
そこで、情報量の多いカラー画像によるデータにて2
値化するため、本出願人は先に特願昭63−263794号にお
いて、カラー撮像手段にて得られたカラー情報の赤色成
分、緑色成分、青色成分の3色成分のサンプルデータを
基準とする2値化条件式から2値化処理する技術を提案
した。
値化するため、本出願人は先に特願昭63−263794号にお
いて、カラー撮像手段にて得られたカラー情報の赤色成
分、緑色成分、青色成分の3色成分のサンプルデータを
基準とする2値化条件式から2値化処理する技術を提案
した。
本発明はこの技術を発展させ、より正確な2値化処理
を実行し、未作業地とそれ以外の地(既作業地等)との
境界線を迅速且つ正確に検出することができるようにす
ることを目的とするものである。
を実行し、未作業地とそれ以外の地(既作業地等)との
境界線を迅速且つ正確に検出することができるようにす
ることを目的とするものである。
そこで本発明は、カラー撮像手段にて未作業地とそれ
以外の地とその両者間の境界線を撮像して得られたカラ
ー画像データの3原色成分を、多階調のデジタル信号に
変換し、未作業地側のサンプリング領域における前記3
原色成分のデジタル信号からしきい値を設定して2値化
処理し、前記境界線を検出する方法において、前記2値
化したのち境界線を求めるためのハフ変換処理にて、そ
のハフ値のピークレベルが最大値または最大値の近傍と
なるように、前記しきい値を適宜変更する評価学習を実
行することを特徴とするものである。
以外の地とその両者間の境界線を撮像して得られたカラ
ー画像データの3原色成分を、多階調のデジタル信号に
変換し、未作業地側のサンプリング領域における前記3
原色成分のデジタル信号からしきい値を設定して2値化
処理し、前記境界線を検出する方法において、前記2値
化したのち境界線を求めるためのハフ変換処理にて、そ
のハフ値のピークレベルが最大値または最大値の近傍と
なるように、前記しきい値を適宜変更する評価学習を実
行することを特徴とするものである。
このように、本発明では、情報量の多いカラー画像デ
ータから3原色成分を多階調のデジタル信号に変換し、
これをもとにしきい値を設定するので、白黒の画像デー
タに比べて、より人間の認識に近いデータ量にて判断す
ることができる。
ータから3原色成分を多階調のデジタル信号に変換し、
これをもとにしきい値を設定するので、白黒の画像デー
タに比べて、より人間の認識に近いデータ量にて判断す
ることができる。
未作業地側のサンプリング領域におけるカラー画像デ
ータを使用して、しきい値を設定するので、サンプリン
グ時の状況や撮像の対象が変化することでカラー画像デ
ータに変動があっても、前記状況や対象に応じたしきい
値を設定でき、固定的なしきい値を設定するのに比べて
状況等の変化に即応して境界線の検出を正確に実行でき
る。
ータを使用して、しきい値を設定するので、サンプリン
グ時の状況や撮像の対象が変化することでカラー画像デ
ータに変動があっても、前記状況や対象に応じたしきい
値を設定でき、固定的なしきい値を設定するのに比べて
状況等の変化に即応して境界線の検出を正確に実行でき
る。
また、2値化したのち境界線を求めるためのハフ変換
処理にて、そのハフ値のピークレベルが最大値または最
大値の近傍となるように、前記しきい値を適宜変更する
評価学習を実行するのであるから、一旦設定したしきい
値のうちから、サンプリングの状況や撮像対象等に応じ
て境界線がより一層正確に計算できるような最適のしき
値を選択し、確定できるのである。
処理にて、そのハフ値のピークレベルが最大値または最
大値の近傍となるように、前記しきい値を適宜変更する
評価学習を実行するのであるから、一旦設定したしきい
値のうちから、サンプリングの状況や撮像対象等に応じ
て境界線がより一層正確に計算できるような最適のしき
値を選択し、確定できるのである。
次に本発明を汎用コンバインに適用した実施例につい
て説明すると、走行作業機であるコンバイン1は走行機
体2の下面に左右一対の履帯式の走行装置3,3を備え、
走行機体2の前進方向に向かって右側前部には座席付き
の操縦部4を配設し、その後方にはエンジン(図示せ
ず)及び穀粒蓄積用のタンク5を備え、走行機体2の左
側には、円周外面にダブルピッチのスクリュー板と適宜
個数の扱歯とを植設した前後長手の扱胴6a等を内蔵した
脱穀部6と、その下方に受け網8と、シーブ等による揺
動選別装置9と、唐箕フアン10の風による風選別装置と
を備える。
て説明すると、走行作業機であるコンバイン1は走行機
体2の下面に左右一対の履帯式の走行装置3,3を備え、
走行機体2の前進方向に向かって右側前部には座席付き
の操縦部4を配設し、その後方にはエンジン(図示せ
ず)及び穀粒蓄積用のタンク5を備え、走行機体2の左
側には、円周外面にダブルピッチのスクリュー板と適宜
個数の扱歯とを植設した前後長手の扱胴6a等を内蔵した
脱穀部6と、その下方に受け網8と、シーブ等による揺
動選別装置9と、唐箕フアン10の風による風選別装置と
を備える。
一方、走行機体2の前面には、前記脱穀部6の前部開
口に連通する角筒状のフイーダハウス11を昇降用油圧シ
リンダ12を介して走行機体に対して昇降自在に装着して
あり、該フイーダハウス11の内部には、左右長手の搬送
板を左右一対のチエンに適宜間隔で取付けたチエンコン
ベヤ13を配設し、このチエンコンベヤ13にて後述の刈取
部15からの刈取穀稈を前記脱穀部6に搬送するものであ
る。
口に連通する角筒状のフイーダハウス11を昇降用油圧シ
リンダ12を介して走行機体に対して昇降自在に装着して
あり、該フイーダハウス11の内部には、左右長手の搬送
板を左右一対のチエンに適宜間隔で取付けたチエンコン
ベヤ13を配設し、このチエンコンベヤ13にて後述の刈取
部15からの刈取穀稈を前記脱穀部6に搬送するものであ
る。
符号7aは、二番樋(図示せず)からの穀粒を脱穀部6
に還元するように搬送する二番還元筒、符号7bは前記タ
ンク5内の穀粒を外部に排出するための排出筒である。
に還元するように搬送する二番還元筒、符号7bは前記タ
ンク5内の穀粒を外部に排出するための排出筒である。
刈取部14は、走行機体2の全幅にわたって左右に延び
るバケット状のプラネットホーム17と、該プラットホー
ム17に走行機体2の全幅にわたって設けたタインバー18
付きのリール19とから成り、プラットホーム17の下面側
には同じく左右長手のバリカン状の刈刃20を有し、プラ
ットホーム17の底板上方には、矢印方向に回転する横長
の掻き込み用のオーガ22を備えてあり、伝動部を介して
前記リール19及びオーガ22を回転駆動し、且つ刈刃20を
駆動する。
るバケット状のプラネットホーム17と、該プラットホー
ム17に走行機体2の全幅にわたって設けたタインバー18
付きのリール19とから成り、プラットホーム17の下面側
には同じく左右長手のバリカン状の刈刃20を有し、プラ
ットホーム17の底板上方には、矢印方向に回転する横長
の掻き込み用のオーガ22を備えてあり、伝動部を介して
前記リール19及びオーガ22を回転駆動し、且つ刈刃20を
駆動する。
なお、符号23は刈取部15の左右両端から前向きにする
分草体であり、また、前記リール19は、穀稈の倒伏状態
等に応じて前後移動調節及び上下揺動調節自在に構成さ
れている。
分草体であり、また、前記リール19は、穀稈の倒伏状態
等に応じて前後移動調節及び上下揺動調節自在に構成さ
れている。
符号25は、第1図等に示すように、圃場に適宜間隔に
て縦横に植立した被検出物である稲の穀桿で、この縦横
に適宜間隔で植立した穀桿25に向かってコンバイン1を
前進させてこれらを刈取るものである。
て縦横に植立した被検出物である稲の穀桿で、この縦横
に適宜間隔で植立した穀桿25に向かってコンバイン1を
前進させてこれらを刈取るものである。
符号24は前記操縦部4のキヤビン前部上端等に取付く
カラービデオカメラ等の撮像手段である。
カラービデオカメラ等の撮像手段である。
また、第3図は、本発明の制御手段のブロック図を示
し、符号27は前記撮像手段24のアナログ入力信号をデジ
タル信号に変換するA/D変換器、符号28は画像メモリ、
符号29は画像処理プロセッサ、符号30は中央処理装置
(CPU)で、符号31は読み出し専用メモリ(ROM)、符号
32は読み書き可能メモリ(RAM)32である。
し、符号27は前記撮像手段24のアナログ入力信号をデジ
タル信号に変換するA/D変換器、符号28は画像メモリ、
符号29は画像処理プロセッサ、符号30は中央処理装置
(CPU)で、符号31は読み出し専用メモリ(ROM)、符号
32は読み書き可能メモリ(RAM)32である。
なお、中央制御装置30の指令信号に従って駆動される
油圧切換弁等の電磁ソレノイドで代表される操舵アクチ
ュエータ33、34にて走行機体2を左右に操舵することが
きる。
油圧切換弁等の電磁ソレノイドで代表される操舵アクチ
ュエータ33、34にて走行機体2を左右に操舵することが
きる。
前記カラービデオカメラ等のカラー撮像手段24は、レ
ンズ系を通して得られた画像(結像)を電気信号に置き
換えるCCD固体撮像素子等から成る二次元撮像デバイス
とからなり、該カラー撮像手段24における二次元撮像デ
バイスは、例えば、1画面につき64×64ドット(縦×
横)の画素を有し、その各画素ごとにRGB表色系〔赤色
(R),緑色(G),青色(B)の原色光とし、加色に
より白色が得られる〕の各色成分の信号に分解すること
ができ、このアナログの色信号をA/D変換器27にて64階
調のデジタル信号に変換する。つまりこのデジタル信号
は、1画素あたり各色成分ごとに64階調(レベル)に量
子化するものであり、ここでのデジタル信号への変換
は、前記一回の撮像ごとにて得られた一つの画面35全体
について実行される(第4図参照)。
ンズ系を通して得られた画像(結像)を電気信号に置き
換えるCCD固体撮像素子等から成る二次元撮像デバイス
とからなり、該カラー撮像手段24における二次元撮像デ
バイスは、例えば、1画面につき64×64ドット(縦×
横)の画素を有し、その各画素ごとにRGB表色系〔赤色
(R),緑色(G),青色(B)の原色光とし、加色に
より白色が得られる〕の各色成分の信号に分解すること
ができ、このアナログの色信号をA/D変換器27にて64階
調のデジタル信号に変換する。つまりこのデジタル信号
は、1画素あたり各色成分ごとに64階調(レベル)に量
子化するものであり、ここでのデジタル信号への変換
は、前記一回の撮像ごとにて得られた一つの画面35全体
について実行される(第4図参照)。
この1画素ごとのデジタル信号は画像メモリ28に格納
された後、画像処理プロセッサ29で所定と範囲のサンプ
リング領域Sについてのみ、以下のしきい値の演算及び
ハフ(hough)変換処理を実行する。
された後、画像処理プロセッサ29で所定と範囲のサンプ
リング領域Sについてのみ、以下のしきい値の演算及び
ハフ(hough)変換処理を実行する。
次に第6図で示すフローチャートに従って、本発明の
画像処理方法について説明する。
画像処理方法について説明する。
まず、走行作業機1を未作業地Mと既作業地Kと境界
線Ko近傍に配置し、スタートに続くステップS1で、i=
0等初期値設定し、ステップS2にてカラー撮像手段24
で、画面35の左右に未作業地Mと既作業地Kとが写りそ
の中程部分に縦方向に前記境界線Koがあるようにして撮
像し、ステップS3では、前記R,G,B3色のデジタル信号を
画像メモリ28に読み込む。
線Ko近傍に配置し、スタートに続くステップS1で、i=
0等初期値設定し、ステップS2にてカラー撮像手段24
で、画面35の左右に未作業地Mと既作業地Kとが写りそ
の中程部分に縦方向に前記境界線Koがあるようにして撮
像し、ステップS3では、前記R,G,B3色のデジタル信号を
画像メモリ28に読み込む。
ステップS4では、画像35のうち未作業地Mの領域から
サンプリング領域Sを特定する。
サンプリング領域Sを特定する。
サンリング領域は、撮像された画面35のうち特定の小
さい領域箇所を予め設定するが、例えば、カラー撮像手
段24で撮像された画面35のうち左略半分の領域に未刈取
地Mが撮像され、右半分には既刈取地Kやその他の背景
が撮像されていると仮定すると、前記左半分の画面を3
×3画素ないし4×4画素の小さい領域に分割し、その
うちの1つないし複数のサンプリング領域Sを特定する
(第4図参照)。
さい領域箇所を予め設定するが、例えば、カラー撮像手
段24で撮像された画面35のうち左略半分の領域に未刈取
地Mが撮像され、右半分には既刈取地Kやその他の背景
が撮像されていると仮定すると、前記左半分の画面を3
×3画素ないし4×4画素の小さい領域に分割し、その
うちの1つないし複数のサンプリング領域Sを特定する
(第4図参照)。
ステップS5にて前記サンプリング領域Sでの前記R,G,
B3色成分デジタル信号の平均値(算術平均値)Rav,Gav,
Bavを計算する。
B3色成分デジタル信号の平均値(算術平均値)Rav,Gav,
Bavを計算する。
前記サンプリング領域Sの大きさが3×3画素である
場合には、その画素の数は9つであり、その各画素ごと
に赤色(R),緑色(G),青色(B)の各色信号のサ
ンプリング出力値を求め、その平均値(算術平均値)を
各々Rav,Gav,Bavとする。
場合には、その画素の数は9つであり、その各画素ごと
に赤色(R),緑色(G),青色(B)の各色信号のサ
ンプリング出力値を求め、その平均値(算術平均値)を
各々Rav,Gav,Bavとする。
ステップS6では前記各サンプリング出力値から、各々
の色成分における標準偏差σr,σg,σbを演算する。
の色成分における標準偏差σr,σg,σbを演算する。
そして、以下のステップS7からステップS13までの処
理実行して、画像処理における2値化処理の際の最適な
しきい値を確定するものである。
理実行して、画像処理における2値化処理の際の最適な
しきい値を確定するものである。
即ち、前記各色のしきい値を設定するとき、前記サン
プリング領域Sのデータはある程度のばらつきがあるの
で、そのばらつきの程度は標準偏差を計算することによ
り得ることができる。
プリング領域Sのデータはある程度のばらつきがあるの
で、そのばらつきの程度は標準偏差を計算することによ
り得ることができる。
この場合、サンプリング数(データ)が充分大きいと
きには、一般に偶然誤差の効果として、通常の正規分布
(ガウス分布)が生じる。
きには、一般に偶然誤差の効果として、通常の正規分布
(ガウス分布)が生じる。
今、一般に測定値の平均値Xavとすると、その平均値
の信頼区間は規定された確立Pで不確かである限界(Xa
v−aσ〜Xav+aσ)を与える(aは係数、σは標準偏
差)。
の信頼区間は規定された確立Pで不確かである限界(Xa
v−aσ〜Xav+aσ)を与える(aは係数、σは標準偏
差)。
例えば、確率P=95.4%のときには、前記係数a=2
となり、すべてのデータあは、95.4%の確立で、Xav−
2σ〜Xav+2σの範囲内に入るものである。
となり、すべてのデータあは、95.4%の確立で、Xav−
2σ〜Xav+2σの範囲内に入るものである。
しかし、サンプリング数が少ないときには、データの
信頼度、ひいてはしきい値の信頼度(正確さ)を通常の
正規分布(ガウス分布)における確率からもとめること
ができない。
信頼度、ひいてはしきい値の信頼度(正確さ)を通常の
正規分布(ガウス分布)における確率からもとめること
ができない。
そこで、前記各サンプリング出力値から、各々の色成
分における標準偏差σr,σg,σbを演算し、各色のしき
い値Rth,Gth,Bthを演算するに際して、 とし、係数aiを適宜間隔で変化させて、その各しきい値
を求めて実際に2値化処理を実行し、その2値「0」,
「1」の値にて画面35上の未作業地Mの領域と既作業地
Kの領域とに区分し、その境界線Koに相当する画面35上
の多数の候補点Sが略直線状等(ある程度のばらつきが
ある)に分布しているとき(第5図参照)、この多数の
候補点Sから仮想境界直線K′を求める直線当てはめ処
理を行って前記境界線を求めるハフ(hough)変換を実
行し、そのうち最も良い結果が得られたしきい値を採用
するのである。
分における標準偏差σr,σg,σbを演算し、各色のしき
い値Rth,Gth,Bthを演算するに際して、 とし、係数aiを適宜間隔で変化させて、その各しきい値
を求めて実際に2値化処理を実行し、その2値「0」,
「1」の値にて画面35上の未作業地Mの領域と既作業地
Kの領域とに区分し、その境界線Koに相当する画面35上
の多数の候補点Sが略直線状等(ある程度のばらつきが
ある)に分布しているとき(第5図参照)、この多数の
候補点Sから仮想境界直線K′を求める直線当てはめ処
理を行って前記境界線を求めるハフ(hough)変換を実
行し、そのうち最も良い結果が得られたしきい値を採用
するのである。
このときの最も良い結果であると評価するのは、前記
係数aiを所定の一定間隔で変化させ、各係数aiのときに
取る各色のしきい値Rth,Gth,Bthの組で2値化を実行
し、その時のハフ(hough)変換で得られるハフ値(ρ
とθが取る値のヒストグラム)のピーク値Pkiとする
と、前記係数aiを適宜間隔で変化させた場合に、そのピ
ーク値Pkiが最大またはその近傍となるときとするもの
である。
係数aiを所定の一定間隔で変化させ、各係数aiのときに
取る各色のしきい値Rth,Gth,Bthの組で2値化を実行
し、その時のハフ(hough)変換で得られるハフ値(ρ
とθが取る値のヒストグラム)のピーク値Pkiとする
と、前記係数aiを適宜間隔で変化させた場合に、そのピ
ーク値Pkiが最大またはその近傍となるときとするもの
である。
ここで、ハフ変換とは、第5図に示すようなx−y座
標系における任意の候補点Sの座標(xi,yi)を通るす
べての直線群は極座標系では、ρ=xicosθ+yisinθで
表現できる(第7図〜第9図参照)。つまり、ひとつの
座標(xi,yi)を通る直線群をρ−θ空間に写像すれ
ば、ρ−θ座標系では一本の軌跡としての曲線を描くこ
とになる。
標系における任意の候補点Sの座標(xi,yi)を通るす
べての直線群は極座標系では、ρ=xicosθ+yisinθで
表現できる(第7図〜第9図参照)。つまり、ひとつの
座標(xi,yi)を通る直線群をρ−θ空間に写像すれ
ば、ρ−θ座標系では一本の軌跡としての曲線を描くこ
とになる。
従って、前記境界線Koに対応して、所定のしきい値で
得られたすべての候補点S(xi,yi)(第5図参照)に
対して、第7図〜第9図に示すようにハフ軌跡線HBの曲
線群が得られる。
得られたすべての候補点S(xi,yi)(第5図参照)に
対して、第7図〜第9図に示すようにハフ軌跡線HBの曲
線群が得られる。
なお、第7図〜第9図では、ρの範囲を+44〜−44ま
で、θを+π/2〜−π/2まで取る。
で、θを+π/2〜−π/2まで取る。
そして、このハフ軌跡線HGの群がρ−θ平面を最も多
く交差する箇所の座標(ρ1,θ1)を見出すと(第8図
〜第9図参照)、ρ1=xicosθ1+yisinθ1としての
一本のハフ軌跡線が特定され、これをx−y座標に逆写
像すると、第5図に示すように、前記境界線Koに対応す
る1本の仮想境界線K′が求められる。
く交差する箇所の座標(ρ1,θ1)を見出すと(第8図
〜第9図参照)、ρ1=xicosθ1+yisinθ1としての
一本のハフ軌跡線が特定され、これをx−y座標に逆写
像すると、第5図に示すように、前記境界線Koに対応す
る1本の仮想境界線K′が求められる。
この場合、前記(1)式における係数aiの大きさによ
り、ρ−θ平面でのハフ軌跡線HGの群の分布状態は各々
異なり、前記(1)式におけるしきい値が適切でないと
きには、x−y平面上での候補点Sがばらつくので、第
7図のように、ハフ軌跡線HGの群がρ−θ平面で交差す
る箇所があまり無く、且つ交差する箇所が集中すること
もない。しきい値の適切さが中程度であるときには、第
8図のように交差する箇所がある程度集中的となり、し
きい値が最適であるときには、x−y平面上での候補点
Sが略直線状(帯状)に並ぶので、第9図のようにハフ
軌跡線HGの群がρ−θ平面で最も多く交差する箇所の座
標(ρ1,θ1)が微小の範囲内におさまり、逆写像すれ
ば前記境界線Koに近い仮想境界線K′としての直線を得
ることができることになる。
り、ρ−θ平面でのハフ軌跡線HGの群の分布状態は各々
異なり、前記(1)式におけるしきい値が適切でないと
きには、x−y平面上での候補点Sがばらつくので、第
7図のように、ハフ軌跡線HGの群がρ−θ平面で交差す
る箇所があまり無く、且つ交差する箇所が集中すること
もない。しきい値の適切さが中程度であるときには、第
8図のように交差する箇所がある程度集中的となり、し
きい値が最適であるときには、x−y平面上での候補点
Sが略直線状(帯状)に並ぶので、第9図のようにハフ
軌跡線HGの群がρ−θ平面で最も多く交差する箇所の座
標(ρ1,θ1)が微小の範囲内におさまり、逆写像すれ
ば前記境界線Koに近い仮想境界線K′としての直線を得
ることができることになる。
なお、第10図は、係数aiをパラメータとして、ある係
数aiであるときに、ハフ軌跡線HGの群がρの値とθの値
に応じて取る分布量V(ヒストグラム)を三次元グラフ
で示したものであり、第11図は、係数aiを変化させた場
合に前記分布量V(ヒストグラム)のうちの最大分布量
V(ρ,θ)が増減変化する状態曲線Pを示し、ある係
数ai(i=1〜5)における各状態曲線P中の頂点をピ
ークレベルPki(Pk1〜Pk5)として示す。
数aiであるときに、ハフ軌跡線HGの群がρの値とθの値
に応じて取る分布量V(ヒストグラム)を三次元グラフ
で示したものであり、第11図は、係数aiを変化させた場
合に前記分布量V(ヒストグラム)のうちの最大分布量
V(ρ,θ)が増減変化する状態曲線Pを示し、ある係
数ai(i=1〜5)における各状態曲線P中の頂点をピ
ークレベルPki(Pk1〜Pk5)として示す。
そしてこのピークレベルPkiが最大のもの、または最
大値近傍となる前記係数aiを検索すれば、その係数aiに
て(1)式に従って計算されるしきい値(Rth,Gth,Bt
h)の組を最適しきい値として確定するように評価学習
させるのである。
大値近傍となる前記係数aiを検索すれば、その係数aiに
て(1)式に従って計算されるしきい値(Rth,Gth,Bt
h)の組を最適しきい値として確定するように評価学習
させるのである。
なお、ステップS7は前記係数aiを所定の間隔で変化さ
せるために実行するものであり、ステップS8は2値化処
理実行、ステップS9はハフ(hough)変換処理実行、ス
テップS10はピークレベルPkiの記憶、ステップS11は前
記係数aiを所定の間隔で変化させるごとにそれに対応す
るピークレベルPkiと前回と係数aiで得られたピークレ
ベルPkiとの大小関係を判断するピークレベルの評価を
実行するものである。
せるために実行するものであり、ステップS8は2値化処
理実行、ステップS9はハフ(hough)変換処理実行、ス
テップS10はピークレベルPkiの記憶、ステップS11は前
記係数aiを所定の間隔で変化させるごとにそれに対応す
るピークレベルPkiと前回と係数aiで得られたピークレ
ベルPkiとの大小関係を判断するピークレベルの評価を
実行するものである。
このステップS11におけるピークレベルの評価におい
て、前記係数aiを例えば0.25の間隔で順次増加させつ
つ、しきい値計算し、その結果から各々ピークレベルPk
iを求め、すべてのピーク値レベルPkiを一旦記憶させた
後、そのうちの最大値を検索するという手法では、メモ
リ容量が増大すると共に計算時間が長引き、結果を迅速
に出すことができないので、いわゆる山登り法を利用し
てピークレベルPkiが最大または最大値の近傍となるし
きい値を検索する。
て、前記係数aiを例えば0.25の間隔で順次増加させつ
つ、しきい値計算し、その結果から各々ピークレベルPk
iを求め、すべてのピーク値レベルPkiを一旦記憶させた
後、そのうちの最大値を検索するという手法では、メモ
リ容量が増大すると共に計算時間が長引き、結果を迅速
に出すことができないので、いわゆる山登り法を利用し
てピークレベルPkiが最大または最大値の近傍となるし
きい値を検索する。
ここで、山登り法とは、ある係数ai(i=1)を出発
点とし、その値におけるしきい値からピークレベルPk1
を求め、次に前記係数aiの値を一つずらせ(適宜の間隔
分だけ増大させてi=2とする)、そのときしきい値か
ら求めたピークレベルPk2と前記Pk1の大小を比較し、Pk
2>Pk1であるときには、大きい方のピークレベルPk2を
記憶させておき、さらに係数ai(i=3)を一つずらせ
て得たピークレベルPk3と前記記憶させたピークレベルP
k2と大小比較を実行し、大きい方を採用するというよう
に、順次高い値のものを採用することにより、ピークレ
ベルの最大値(山の頂点)付近を見出すものであり、今
回のピークレベルが前回のものより小さい値のときには
前回の値の箇所が山の頂点付近であったと判断するもの
である。この方法ではピークレベルの最大値そのものを
正確に出すことはできないが、最大の近傍であることを
迅速に検索することができる。
点とし、その値におけるしきい値からピークレベルPk1
を求め、次に前記係数aiの値を一つずらせ(適宜の間隔
分だけ増大させてi=2とする)、そのときしきい値か
ら求めたピークレベルPk2と前記Pk1の大小を比較し、Pk
2>Pk1であるときには、大きい方のピークレベルPk2を
記憶させておき、さらに係数ai(i=3)を一つずらせ
て得たピークレベルPk3と前記記憶させたピークレベルP
k2と大小比較を実行し、大きい方を採用するというよう
に、順次高い値のものを採用することにより、ピークレ
ベルの最大値(山の頂点)付近を見出すものであり、今
回のピークレベルが前回のものより小さい値のときには
前回の値の箇所が山の頂点付近であったと判断するもの
である。この方法ではピークレベルの最大値そのものを
正確に出すことはできないが、最大の近傍であることを
迅速に検索することができる。
なお、前記ρ−θ平面上でのハフ軌跡線HGの曲線群が
最も多く集まる箇所を検出するため、ρ−θ平面を縦横
(二次元)に細かく分割して検索すると、この分割数に
相当する大きな容量の二次元配列メモリが必要となる
し、ρ−θ平面全体についてピークレベルPkiの記憶さ
れると、不必要な箇所の境界線(例えば中割り作業時の
他方の境界線)も記憶させるという不都合が生じると共
に演算時間が長引くので、ρの値をminρ〜maxρの範囲
内に限定するようにその最大値と最小値とを適宜限定す
るのが好ましい。
最も多く集まる箇所を検出するため、ρ−θ平面を縦横
(二次元)に細かく分割して検索すると、この分割数に
相当する大きな容量の二次元配列メモリが必要となる
し、ρ−θ平面全体についてピークレベルPkiの記憶さ
れると、不必要な箇所の境界線(例えば中割り作業時の
他方の境界線)も記憶させるという不都合が生じると共
に演算時間が長引くので、ρの値をminρ〜maxρの範囲
内に限定するようにその最大値と最小値とを適宜限定す
るのが好ましい。
ステップS12にて前記ステップS11のピークレベルの評
価、つまりピークレベルが最大値又はその近傍となる条
件が満足するか否かを判別し、条件満足しない(no)と
きには、でステップS6の後に戻して、係数aiを一定間
隔変化させてしきい値計算、2値化処理、ハフ(houg
h)変換処理、ピークレベルPkiの記憶、ピークレベルの
評価を繰り返す。
価、つまりピークレベルが最大値又はその近傍となる条
件が満足するか否かを判別し、条件満足しない(no)と
きには、でステップS6の後に戻して、係数aiを一定間
隔変化させてしきい値計算、2値化処理、ハフ(houg
h)変換処理、ピークレベルPkiの記憶、ピークレベルの
評価を繰り返す。
前記ステップS12にて条件満足する(yes)ときには、
ステップS13にてしきい値を確定させる。
ステップS13にてしきい値を確定させる。
その後は、走行作業機1を前進させつつ、撮像手段24
で読み込んだ未作業地と既作業地との画像情報を、前記
確定したしきい地でもって2値化処理し、未作業地Mと
既作業地Kとの境界線Koを決定すれば良いのである。
で読み込んだ未作業地と既作業地との画像情報を、前記
確定したしきい地でもって2値化処理し、未作業地Mと
既作業地Kとの境界線Koを決定すれば良いのである。
なお、2値化画像に含まれている雑音成分、穴、欠け
を埋める整形処理のため、注目する画素の外周に隣接す
る画素の連結関係(8連結論理和)で調べ、拡散、侵
食、孤立点除去の手法で処理した。
を埋める整形処理のため、注目する画素の外周に隣接す
る画素の連結関係(8連結論理和)で調べ、拡散、侵
食、孤立点除去の手法で処理した。
この2値化処理は、x−y平面の画面47の全画素につ
いて実行しても良いが、不必要な対象の箇所を2値化す
ることで、直線化の誤差が大きくなるから、一定のウイ
ンドウ処理(x−y平面である枠内の領域のみの2値化
処理)を実行しても良いのである。
いて実行しても良いが、不必要な対象の箇所を2値化す
ることで、直線化の誤差が大きくなるから、一定のウイ
ンドウ処理(x−y平面である枠内の領域のみの2値化
処理)を実行しても良いのである。
このようにすれば、画面35の画素全体にわたって2値
化する場合に要するメモリよりメモリ量を削減できる
し、境界線検出までの処理時間を大幅に短縮できるので
ある。
化する場合に要するメモリよりメモリ量を削減できる
し、境界線検出までの処理時間を大幅に短縮できるので
ある。
なお、撮像して検出する境界線の対象は、未作業地と
既作業地との境界線ばかりでなく、未作業地と畦(コン
クリート製畦を含む)との境界線等もある。
既作業地との境界線ばかりでなく、未作業地と畦(コン
クリート製畦を含む)との境界線等もある。
図面は本発明の実施例を示し、第1図はコンバインの平
面図、第2図は側面図、第3図は制御装置のブロック
図、第4図は既作業地と未作業地と境界線を撮像した画
面の説明図、第5図はx−y平面における候補点を示し
た画面の説明図、第6図はフローチャート、第7図と第
8図と第9図は各々ρ−θ平面におけるハフ軌跡線HGの
説明図、第10図はハフ値(ρとθ)に対する分布量V
(ヒストグラム)を示す図、第11図は係数aiの変化に対
応する分布量V(ヒストグラム)のピークレベルを示す
図である。 1……作業機、3,3……走行装置、6……脱穀部、15…
…刈取部、24……撮像手段、30……中央制御装置、27…
…A/D変換器、28……画像メモリ、29……画像プロセッ
サ、35……画面。
面図、第2図は側面図、第3図は制御装置のブロック
図、第4図は既作業地と未作業地と境界線を撮像した画
面の説明図、第5図はx−y平面における候補点を示し
た画面の説明図、第6図はフローチャート、第7図と第
8図と第9図は各々ρ−θ平面におけるハフ軌跡線HGの
説明図、第10図はハフ値(ρとθ)に対する分布量V
(ヒストグラム)を示す図、第11図は係数aiの変化に対
応する分布量V(ヒストグラム)のピークレベルを示す
図である。 1……作業機、3,3……走行装置、6……脱穀部、15…
…刈取部、24……撮像手段、30……中央制御装置、27…
…A/D変換器、28……画像メモリ、29……画像プロセッ
サ、35……画面。
───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.6,DB名) G05D 1/02 A01B 69/00 303 G06T 1/00 G06T 9/20
Claims (1)
- 【請求項1】カラー撮像手段にて未作業地とそれ以外の
地とその両者間の境界線を撮像して得られたカラー画像
データの3原色成分を、多階調のデジタル信号に変換
し、未作業地側のサンプリング領域における前記3原色
成分のデジタル信号からしきい値を設定して2値化処理
し、前記境界線を検出する方法において、前記2値化し
たのち境界線を求めるためのハフ変換処理にて、そのハ
フ値のピークレベルが最大値または最大の近傍となるよ
うに、前記しきい値を適宜変更する評価学習を実行する
ことを特徴とするカラー情報による境界線の検出方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP1065153A JP2795454B2 (ja) | 1989-03-16 | 1989-03-16 | カラー情報による境界線の検出方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP1065153A JP2795454B2 (ja) | 1989-03-16 | 1989-03-16 | カラー情報による境界線の検出方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH02242405A JPH02242405A (ja) | 1990-09-26 |
JP2795454B2 true JP2795454B2 (ja) | 1998-09-10 |
Family
ID=13278652
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP1065153A Expired - Lifetime JP2795454B2 (ja) | 1989-03-16 | 1989-03-16 | カラー情報による境界線の検出方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2795454B2 (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06176159A (ja) * | 1992-12-10 | 1994-06-24 | Nec Corp | カラー画像対応型物体認識装置 |
DE19634402C1 (de) * | 1996-08-26 | 1997-07-24 | Koenig & Bauer Albert Ag | Kurvenscheibe und Verfahren zu deren Herstellung |
WO2017048491A1 (en) | 2015-09-16 | 2017-03-23 | Thermo Electron Scientific Instruments Llc | Image analysis system and method |
WO2020134236A1 (zh) * | 2018-12-29 | 2020-07-02 | 丰疆智能科技股份有限公司 | 收割机及其自动驾驶方法 |
-
1989
- 1989-03-16 JP JP1065153A patent/JP2795454B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JPH02242405A (ja) | 1990-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20180084708A1 (en) | Agricultural work machine for avoiding anomalies | |
CN110243372B (zh) | 基于机器视觉的智能农机导航***及方法 | |
EP1655620B1 (en) | Obstacle detection using stereo vision | |
CN105989601B (zh) | 基于机器视觉的农业agv玉米行间导航基准线提取方法 | |
US8185275B2 (en) | System for vehicular guidance with respect to harvested crop | |
US7580549B2 (en) | Method and system for vehicular guidance using a crop image | |
AU2006202812A1 (en) | Method and system for vehicular guidance using a crop image | |
CN110969654A (zh) | 基于收割机的玉米高通量表型测量的方法及装置、收割机 | |
CN112395984A (zh) | 一种用于无人驾驶农业机械的苗列导航线检测方法 | |
CN114092822B (zh) | 图像处理方法、移动控制方法以及移动控制*** | |
JP2795454B2 (ja) | カラー情報による境界線の検出方法 | |
CN114387343A (zh) | 一种成簇平菇的菌柄采摘位置的检测方法 | |
JP3505372B2 (ja) | 月、惑星軟着陸のための画像処理装置及び画像処理方法 | |
JPH04293404A (ja) | 自動走行作業機における画像処理方法 | |
CN115451965B (zh) | 基于双目视觉的插秧机插植***相对航向信息检测方法 | |
Kaur et al. | Automatic crop furrow detection for precision agriculture | |
CN113065501A (zh) | 苗线识别的方法、装置和农机 | |
JP2740268B2 (ja) | 走行作業機の操舵制御における画像情報の2値化処理装置 | |
JPH04293403A (ja) | 自動走行作業機における画像処理方法 | |
Mei et al. | Combing color index and region growing with simple non-iterative clustering for plant segmentation | |
US20230403964A1 (en) | Method for Estimating a Course of Plant Rows | |
JP2828520B2 (ja) | 作業機用の作物列検出装置 | |
JP2513374Y2 (ja) | 収穫機の自動操向装置 | |
JP2583583B2 (ja) | 農作業機における自動燥舵制御装置 | |
Wang et al. | Detection of navigation route in greenhouse environment with machine vision |