(57)【要約】【課題】  航空写真や衛星画像から海岸線のポイントデータを抽出する岸線抽出装置に係り、地図や海図の作成時の負荷を軽減し、更に、さんご礁のように輝度が高い海域部分の影響を除去することを課題とする。【解決手段】  バンド分解処理Ｓ２０で、衛星写真等をＲ，Ｇ，Ｂ等のバンドに分解し、バンド間演算２値化処理Ｓ４０で、緑と赤のバンドの値を比較し、緑の画素値が大きい場合に水面領域の候補と判定し、赤の画素値が大きい場合に陸領域の候補と判定し、２値化画像を生成する。輪郭追跡Ｓ５０で、２値化画像中の閉じた領域の輪郭追跡を行い、冗長ポイント除去Ｓ６０で、画像の端の座標を排除し、岸線ポイント出力Ｓ７０で、緯度及び経度のポイントを出力する。【選択図】        図１
2005148906.tif 000002
【特許請求の範囲】
【請求項１】
以下の要素を有することを特徴とする岸線抽出装置
（１）高度上方から撮影した水面領域及び陸領域を含むカラー画像又はモノクロ画像、あるいは水面領域又は陸領域からなるカラー画像又はモノクロ画像を入力する画像入力部
（２）入力したカラー画像を２つ以上のバンドに分解するバンド分解部
（３）分解した２以上のバンドの画素値を演算し、演算結果を２値化して２値化画像を生成するバンド間演算２値化部
（４）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する輪郭追跡部。
【請求項２】
以下の要素を有することを特徴とする岸線抽出装置
（１）高度上方から撮影した水面領域及び陸領域を含む入力画像、あるいは水面領域又は陸領域からなる入力画像を入力する画像入力部
（２）入力画像の輝度について、所定の閾値で２値化して２値化画像を生成する２値化生成部
（３）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する輪郭追跡部。
【請求項３】
更に、抽出した座標のうち、画像の端の直線部分の座標を除外する冗長ポイント除去部を有することを特徴とする請求項１または２に記載の岸線抽出装置。
【請求項４】
更に、抽出した座標を緯度及び経度に変換し、変換した緯度及び経度を岸線ポイントとして出力する岸線ポイント出力部を有することを特徴とする請求項１から３までのいずれかに記載の岸線抽出装置。
【請求項５】
バンド分解部は、入力したカラー画像を、少なくとも緑のバンドと赤のバンドに分解し、
バンド間演算２値化部は、各画素について緑のバンドと赤のバンドの値を比較し、緑のバンドの画素値が大きい場合に、水面領域の候補と判定し、赤のバンドの画素値が大きい場合に、陸領域の候補と判定することを特徴とする請求項１記載の岸線抽出装置。
【請求項６】
以下の要素を有することを特徴とする岸線抽出方法
（１）高度上方から撮影した水面領域及び陸領域を含むカラー画像又はモノクロ画像、あるいは水面領域又は陸領域からなるカラー画像又はモノクロ画像を入力する工程
（２）入力したカラー画像を２つ以上のバンドに分解する工程
（３）分解した２以上のバンドの画素値を演算し、演算結果を２値化して２値化画像を生成する工程
（４）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する工程。
【請求項７】
以下の要素を有することを特徴とする岸線抽出方法
（１）高度上方から撮影した水面領域及び陸領域を含む入力画像、あるいは水面領域又は陸領域からなる入力画像を入力する工程
（２）入力画像の輝度について、所定の閾値で２値化して２値化画像を生成する工程
（３）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する工程。
【請求項８】
岸線抽出装置となるコンピュータに、以下の処理を実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体
（１）高度上方から撮影した水面領域及び陸領域を含むカラー画像又はモノクロ画像、あるいは水面領域又は陸領域からなるカラー画像又はモノクロ画像を入力する処理
（２）入力したカラー画像を２つ以上のバンドに分解する処理
（３）分解した２以上のバンドの画素値を演算し、演算結果を２値化して２値化画像を生成する処理
（４）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する処理。
【請求項９】
岸線抽出装置となるコンピュータに、以下の処理を実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体
（１）高度上方から撮影した水面領域及び陸領域を含む入力画像、あるいは水面領域又は陸領域からなる入力画像を入力する処理
（２）入力画像の輝度について、所定の閾値で２値化して２値化画像を生成する処理
（３）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する処理。
【請求項１０】
岸線抽出装置となるコンピュータに、以下の手順を実行させるためのプログラム
（１）高度上方から撮影した水面領域及び陸領域を含むカラー画像又はモノクロ画像、あるいは水面領域又は陸領域からなるカラー画像又はモノクロ画像を入力する手順
（２）入力したカラー画像を２つ以上のバンドに分解する手順
（３）分解した２以上のバンドの画素値を演算し、演算結果を２値化して２値化画像を生成する手順
（４）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する手順。
【請求項１１】
岸線抽出装置となるコンピュータに、以下の手順を実行させるためのプログラム
（１）高度上方から撮影した水面領域及び陸領域を含む入力画像、あるいは水面領域又は陸領域からなる入力画像を入力する手順
（２）入力画像の輝度について、所定の閾値で２値化して２値化画像を生成する手順
（３）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する手順。
【請求項１２】
岸線抽出装置は、輪郭座標を抽出した後、少なくとも面積及び周囲長に基づいて岸線判定を行うことを特徴とする請求項１または２に記載の岸線抽出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、航空写真や衛星画像から海岸線のポイントデータを抽出する岸線抽出装置に係り、地図や海図の作成時の負荷を軽減する技術に関する。
【背景技術】
【０００２】
従来、地図を作成するためには、現地に赴いて測地調査を行う必要あったが、近年では航空写真や衛星画像が用いられるようになっている。地図作成において、陸と海の境界である岸線部分を作成するには、航空写真や衛星画像から人の手によって海岸線を判断し、各岸線ポイントを指定しなければならなかった。
【０００３】
また、通常、海と陸では輝度レベルが異なるため、２値化することにより岸線部分を抽出できる。しかしながら２値化したのみではラスターデータのままであるため、岸線のポイントデータとしての利用はできない。
【０００４】
また、２値化の際、さんご礁が群生するような地域では海域部分の輝度が高くなることがあるため、海と陸が上手く分離できない場合がある。
【０００５】
図１６は、陸域と海域を含む入力画像を示す図である。陸域と海域の輝度は、陸域＞海域であるため、ある閾値で２値化することにより、陸域と海域を分離することができる。しかしながら、さんご礁が群生するような低緯度地域においては、海域の輝度が高いため、単純な２値化では陸域と海域を上手く分離できない場合がある。
【０００６】
図１７は、２値化画像を示す図である。図に示すように、輝度がＡ≒Ｂ＞＞Ｃであると２値化した場合、陸と海を上手く分離できない。
【特許文献１】特開平８−０１６７７７号公報
【特許文献２】特開平９−３３０４１８号公報
【特許文献３】特開２０００−１３２６６６号公報
【特許文献４】特開２００１−０６６９８４号公報
【特許文献５】特開２００２−２５１６１８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
本発明は、上記した従来技術の欠点を除くためになされたものであって、その目的とするところは、航空写真や衛星画像を用いて岸線ポイントを自動的に抽出することにより、上記の手間を軽減する装置を提供することである。
【０００８】
特に、２値化画像から輪郭を抽出する技法を用いることにより、岸線のポイントデータを得ることを課題とする。
【０００９】
更に、２値化処理時に所定のバンド間演算を行うことにより、さんご礁のように輝度が高い海域部分の影響を除去することを課題とする。
【課題を解決するための手段】
【００１０】
本発明に係る岸線抽出装置は、
以下の要素を有することを特徴とする
（１）高度上方から撮影した水面領域及び陸領域を含むカラー画像又はモノクロ画像、あるいは水面領域又は陸領域からなるカラー画像又はモノクロ画像を入力する画像入力部
（２）入力したカラー画像を２つ以上のバンドに分解するバンド分解部
（３）分解した２以上のバンドの画素値を演算し、演算結果を２値化して２値化画像を生成するバンド間演算２値化部
（４）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する輪郭追跡部。
【００１１】
本発明に係る岸線抽出装置は、
以下の要素を有することを特徴とする
（１）高度上方から撮影した水面領域及び陸領域を含む入力画像、あるいは水面領域又は陸領域からなる入力画像を入力する画像入力部
（２）入力画像の輝度について、所定の閾値で２値化して２値化画像を生成する２値化生成部
（３）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する輪郭追跡部。
【００１２】
更に、抽出した座標のうち、画像の端の直線部分の座標を除外する冗長ポイント除去部を有することを特徴とする。
【００１３】
更に、抽出した座標を緯度及び経度に変換し、変換した緯度及び経度を岸線ポイントとして出力する岸線ポイント出力部を有することを特徴とする。
【００１４】
バンド分解部は、入力したカラー画像を、少なくとも緑のバンドと赤のバンドに分解し、
バンド間演算２値化部は、各画素について緑のバンドと赤のバンドの値を比較し、緑のバンドの画素値が大きい場合に、水面領域の候補と判定し、赤のバンドの画素値が大きい場合に、陸領域の候補と判定することを特徴とする。
【００１５】
本発明に係る岸線抽出方法は、
以下の要素を有することを特徴とする
（１）高度上方から撮影した水面領域及び陸領域を含むカラー画像又はモノクロ画像、あるいは水面領域又は陸領域からなるカラー画像又はモノクロ画像を入力する工程
（２）入力したカラー画像を２つ以上のバンドに分解する工程
（３）分解した２以上のバンドの画素値を演算し、演算結果を２値化して２値化画像を生成する工程
（４）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する工程。
【００１６】
本発明に係る岸線抽出方法は、
以下の要素を有することを特徴とする
（１）高度上方から撮影した水面領域及び陸領域を含む入力画像、あるいは水面領域又は陸領域からなる入力画像を入力する工程
（２）入力画像の輝度について、所定の閾値で２値化して２値化画像を生成する工程
（３）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する工程。
【００１７】
本発明に係るプログラムを記録したコンピュータ読み取り可能な記録媒体は、
岸線抽出装置となるコンピュータに、以下の処理を実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体であることを特徴とする
（１）高度上方から撮影した水面領域及び陸領域を含むカラー画像又はモノクロ画像、あるいは水面領域又は陸領域からなるカラー画像又はモノクロ画像を入力する処理
（２）入力したカラー画像を２つ以上のバンドに分解する処理
（３）分解した２以上のバンドの画素値を演算し、演算結果を２値化して２値化画像を生成する処理
（４）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する処理。
【００１８】
本発明に係るプログラムを記録したコンピュータ読み取り可能な記録媒体は、
岸線抽出装置となるコンピュータに、以下の処理を実行させるためのプログラムを記録したコンピュータ読み取り可能な記録媒体であることを特徴とする
（１）高度上方から撮影した水面領域及び陸領域を含む入力画像、あるいは水面領域又は陸領域からなる入力画像を入力する処理
（２）入力画像の輝度について、所定の閾値で２値化して２値化画像を生成する処理
（３）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する処理。
【００１９】
本発明に係るプログラムは、
岸線抽出装置となるコンピュータに、以下の手順を実行させるためのプログラムであることを特徴とする
（１）高度上方から撮影した水面領域及び陸領域を含むカラー画像又はモノクロ画像、あるいは水面領域又は陸領域からなるカラー画像又はモノクロ画像を入力する手順
（２）入力したカラー画像を２つ以上のバンドに分解する手順
（３）分解した２以上のバンドの画素値を演算し、演算結果を２値化して２値化画像を生成する手順
（４）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する手順。
【００２０】
本発明に係るプログラムは、
岸線抽出装置となるコンピュータに、以下の手順を実行させるためのプログラムであることを特徴とする
（１）高度上方から撮影した水面領域及び陸領域を含む入力画像、あるいは水面領域又は陸領域からなる入力画像を入力する手順
（２）入力画像の輝度について、所定の閾値で２値化して２値化画像を生成する手順
（３）生成した２値化画像中の閉じた領域の輪郭追跡を行い、輪郭を構成する座標を抽出する手順。
【００２１】
本発明に係る岸線抽出装置は、
岸線抽出装置は、輪郭座標を抽出した後、少なくとも面積及び周囲長に基づいて岸線判定を行うことを特徴とする。
【発明の効果】
【００２２】
本発明においては、２値化画像から輪郭を抽出する技法を用いるので、岸線のポイントデータを自動的に得ることができる。更に、２値化処理時に所定のバンド間演算を行うので、さんご礁のように輝度が高い海域部分の影響を除去することができる。
【発明を実施するための最良の形態】
【００２３】
実施の形態１．
以下本発明を図面に示す実施例に基づいて説明する。図１は、岸線抽出装置の処理フローを示す図である。
【００２４】
まず、Ｓ１０で画像入力部により、モノクロあるいはカラー（Ｔｒｕｅ，Ｆａｌｓｅ）画像を読み込む。次に、Ｓ２０でバンド分解処理部により、カラー画像の場合、各バンド（Ｒ（Ｒｅｄ：赤），Ｇ（Ｇｒｅｅｎ：緑），Ｂ（Ｂｌｕｅ：青），ＮＩＲ（近赤外））に分解する。Ｓ３０で２値化部により、入力画像に対して所定のしきい値で２値化を行い、２値画像を生成する。更に、Ｓ４０でバンド間演算２値化処理部により、バンド間の演算を伴う２値化を行い、２値画像を生成する。Ｓ５０で輪郭追跡部により、２値画像中の閉じた領域（図形と呼称）の輪郭追跡を行い、輪郭を構成する座標を抽出する。そして、画像の端の直線部分を岸線座標として抽出した場合には、Ｓ６０で冗長ポイント除去部により、当該両端の座標を除去する。最後に、Ｓ７０で岸線ポイント出力部により、抽出座標を緯度・経度に変換後、岸線ポイントして出力する。Ｓ５０は、Ｓ３０とＳ４０の処理により得られる２値化画像のいずれかを用いる。
【００２５】
続いて、処理の詳細について説明する。まず、バンド間演算２値化（Ｓ４０）について詳述する。この処理では、さんご礁域のような２値化に影響する海域の輝度の高い部分を除去するために、バンド間演算を伴った２値化処理を行っている。
【００２６】
図２は、バンド間演算２値化の処理フローを示す図である。まず、Ｓ４０１で比画像作成（Ｇ画像／Ｒ画像）を行う。この処理では、Ｇ画像／Ｒ画像による比画像を作成する。つまり、各画素毎に、Ｇ値をＲ値で割った商を求め、その商の値を画素値とする演算画像（比画像と呼ぶ。）を生成する。
【００２７】
図３は、比較画像作成の概要を示す図である。入力画像中の対象物の特徴を以下に示す。
（１）陸域は、Ｇ（緑）バンドより、Ｒ（赤）バンドの輝度が高い
（２）さんご礁域は、Ｒ（赤）バンドより、Ｇ（緑）バンドの輝度が高い
従って、Ｇ／Ｒの演算により、陸の部分は輝度＜１となり、それ以外は輝度＞１となる画像を作成することができる。
【００２８】
この演算画像では上図のように海の部分に輝度＜１の部分があってもかまわない。この部分は、後の処理で除去されるからである。また、Ｒ画像の代わりに、ＮＩＲ（近赤外）画像を用いてもよい。つまり、Ｇ／ＮＩＲ＞１の領域とＧ／ＮＩＲ＜１の領域に分けた演算画像を生成する。
【００２９】
そして、Ｓ４０２で比画像の２値化を行う。この処理では、比画像を所定の閾値で２値化し、陸とそれ以外の２値画像を作成する。つまり、１を閾値として、Ｇ／Ｒ＞１の領域とＧ／Ｒ＜１の領域に分けた画像を生成する。図４は、比較画像の２値化の概要を示す図である。この例によらず、１以外の閾値を設定することも考えられる。
【００３０】
Ｓ４０３でＲ画像の２値化を行う。この処理では、Ｒ画像を所定のしきい値で２値化し、海とそれ以外の２値画像を作成する。この例によらず、Ｒ画像の代わりにＮＩＲ画像を用いてもよい。図５は、Ｒ画像の２値化の概要を示す図である。
【００３１】
次に、Ｓ４０４で反転画像作成を行う。この処理では、Ｓ４０３で生成した２値画像の反転画像を作成する。図６は、反転画像作成の概要を示す図である。この例では、ここでは、黒を輝度０、白を輝度１としている。但し、Ｓ４０３の２値化の際、輝度の低い方を輝度１で２値化した場合は、本処理は必要ない。
【００３２】
そして、Ｓ４０５で合成２値画像作成（ＯＲ演算）を行う。この処理では、Ｓ４０２で生成した画像とＳ４０４で生成した画像をＯＲ演算し、合成２値画像を作成する。図７は、合成２値画像作成の概要を示す図である。ここでは、黒を輝度０、白を輝度１としている。
【００３３】
Ｓ４０５の合成処理で除去できないノイズは、公知の方法（メディアンフィルタ等）を用いて除去する（Ｓ４０６）。
【００３４】
続いて、輪郭追跡（Ｓ５０）の処理について詳述する。図８は、輪郭追跡の処理フローを示す図である。
【００３５】
まず、Ｓ５０１で２値図形輪郭追跡を行う。図９は、２値画像の走査方式を示す図である。２値画像を図９のような走査方式で画像を走査し、図形の始点を見つけたら、１本の輪郭線を追跡する。その後、再び走査を続け、次の図形の始点を探す。この操作を繰り返し、画像中の全ての図形（内部図形も含む）の輪郭の座標を抽出し、テンポラリ配列に格納する。さらに周囲長あるいは面積も併せて算出する。なお具体的な追跡方法、周囲長、面積算出については、公知な技術であるため省略する。図１０は、輪郭追跡後の輪郭座標を示す図である。
【００３６】
Ｓ５０２のループ１の処理により、上述の図形の中から、周囲長や面積を基に岸線として妥当でない図形を除去する。従って、ループカウンタは図形数に対応する。具体的には、Ｓ５０２で周囲長（あるいは面積）の判定を行う。つまり、Ｓ５０１で算出した周囲長（あるいは面積）を設定値と比較する。そして、設定値よりも小さい場合には、Ｓ５０４の座標削除により、図形の座標が格納されているテンポラリ配列から、当該図形の全座標を削除する。
【００３７】
Ｓ５０５のループ２の処理により、１画素ずつ抽出された座標を間引くことにより、全体のデータ数を削減する。ループカウンタは、図形数に対応する。そのため、Ｓ５０６のカウンタ変数初期化で、カウンタ変数Ｋ＝１に設定する。
【００３８】
Ｓ５０７のループ３の処理により、各図形毎に間引いた座標を岸線配列１に格納する。ステップは、設定間引き数（＝１の場合、間引きなし）であり、ループカウンタは、図形の座標数に対応する。Ｓ５０８の座標格納で、岸線第一配列（Ｋ）のＸ座標にテンポラリ配列（Ｊ）のＸ座標を格納し、岸線第二配列（Ｋ）のＹ座標にテンポラリ配列（Ｊ）のＹ座標を格納する。そして、Ｓ５０９のカウンタインクリメントで、カウンタ変数をインクリメント（＋１）する。
【００３９】
続いて、冗長ポイント除去（Ｓ６０）について詳述する。前述の輪郭追跡部による処理（Ｓ５０）では、２値図形の輪郭を全て追跡するため、画像の端も抽出する。しかしながら、画像の端は、岸線を構成するものではないため、詳細な座標情報は必要とせず、それぞれのコーナー座標があれば十分である。そこで本装置では、冗長ポイント除去部により、画像の端の余分な抽出点を削除している。図１１と図１２は、冗長ポイント削除前の輪郭座標と冗長ポイント削除後の輪郭座標を示している。
【００４０】
図１３から図１５は、冗長ポイント除去の処理フローを示す図である。まず、Ｓ６０１のフラグ変数初期化で、冗長点除去用のフラグ変数を初期化する。ＦＬＧに０を、ＦＬＧ＿ＯＬＤに０を設定する。
【００４１】
Ｓ６０２のループ１の処理では、各図形の座標の中から、画像の端の部分を除去する。従って、ループカウンタは図形数に対応する。続いて、Ｓ６０３のカウンタ変数初期化で、カウンタ変数Ｋに１を設定する。
【００４２】
そして、Ｓ６０４のループ３で、岸線第一配列の座標を、画像の端を除いて岸線第二配列に格納する。ループカウンタは、図形の座標数に対応する。Ｓ６０５の座標判定で、Ｘ、Ｙの座標が、画像の端（上下左右）か否かを判定する。端である場合には、Ｓ６０６のフラグ設定で、座標に応じた値をフラグ変数に設定する。
Ｘ＝０（左端）→     ＦＬＧ＝１
Ｘ＝Ｘサイズ−１（右端）    →     ＦＬＧ＝２
Ｙ＝０（上端）→     ＦＬＧ＝３
Ｙ＝Ｙサイズ−１（下端）    →     ＦＬＧ＝４
Ｓ６０７のフラグ判定で、ＦＬＧとＦＬＧ＿ＯＬＤの値を比較する。ＦＬＧ＝ＦＬＧ＿ＯＬＤならば、画像の端部分の続きであり、ＦＬＧ≠ＦＬＧ＿ＯＬＤならば、画像の端部分の開始座標となる。
【００４３】
画像の端部分の開始座標である場合には、Ｓ６０８の座標格納で、岸線第二配列（Ｋ）のＸ座標に岸線第一配列（Ｊ）のＸ座標を設定し、岸線第二配列（Ｋ）のＹ座標に岸線第一配列（Ｊ）のＹ座標を設定する。そして、Ｓ６０９のカウンタインクリメントで、カウンタ変数をインクリメント（＋１）する。そして、６０１０のフラグ変数格納で、次回のフラグ判定のため、現在の値を格納する。つつまり、ＦＬＧ＿ＯＬＤにＦＬＧの値を設定する。
【００４４】
Ｓ６０５でＸ、Ｙの座標が画像の端（上下左右）でないと判定した場合には、図１４に示した処理を行う。まず、Ｓ６０１１のフラグ変数設定判定で、フラグ変数が設定されているか否か判定する。設定されている場合は、１つ前の座標が画像の端にあったことになるため、Ｊ−１の座標も格納する必要がある。従って、Ｓ６０１２の座標格納で、画像の端の最後の座標を格納する。つまり、岸線第二配列（Ｋ）のＸ座標に、岸線第一配列（Ｊ−１）のＸ座標を設定し、岸線第二配列（Ｋ）のＹ座標に、岸線第一配列（Ｊ−１）のＹ座標を設定する。そして、Ｓ６０１３のカウンタインクリメントで、カウンタ変数をインクリメント（＋１）し、Ｓ６０１４のフラグ変数をクリアする。ＦＬＧに０を設定し、ＦＬＧ＿ＯＬＤに０を設定する。
【００４５】
そして、Ｓ６０１５の座標格納で、岸線第二配列（Ｋ）のＸ座標に岸線第一配列（Ｊ）のＸ座標を設定し、岸線第二配列（Ｋ）のＹ座標に岸線第一配列（Ｊ）のＹ座標を設定し、Ｓ６０１６のカウンタインクリメントで、カウンタ変数をインクリメント（＋１）する。
【００４６】
岸線抽出装置は、コンピュータであり、各要素はプログラムにより処理を実行することができる。また、プログラムを記憶媒体に記憶させ、記憶媒体からコンピュータに読み取られるようにすることができる。
【００４７】
上述の例では、比画像作成の際に、Ｒ（赤）とＧ（緑）を比較したが、Ｒ（赤）とＢ（青）を比較する形態や、Ｂ（青）とＧ（緑）を比較する形態も有効である。また、Ｒ（赤）に代えてＮＩＲ（近赤外）を用いることも有効である。
【図面の簡単な説明】
【００４８】
【図１】岸線抽出装置の処理フローを示す図である。
【図２】バンド間演算２値化の処理フローを示す図である。
【図３】比較画像作成の概要を示す図である。
【図４】比較画像の２値化の概要を示す図である。
【図５】Ｒ画像の２値化の概要を示す図である。
【図６】反転画像作成の概要を示す図である。
【図７】合成２値画像作成の概要を示す図である。
【図８】輪郭追跡の処理フローを示す図である。
【図９】２値画像の走査方式を示す図である。
【図１０】輪郭追跡後の輪郭座標を示す図である。
【図１１】冗長ポイント削除前の輪郭座標を示す図である。
【図１２】冗長ポイント削除後の輪郭座標を示す図である。
【図１３】冗長ポイント除去の処理フロー（１／３）を示す図である。
【図１４】冗長ポイント除去の処理フロー（２／３）を示す図である。
【図１５】冗長ポイント除去の処理フロー（３／３）を示す図である。
【図１６】陸域と海域を含む入力画像を示す図である。
【図１７】２値化画像を示す図である。
【図１】
2005148906.tif 000003
【図２】
2005148906.tif 000004
【図３】
2005148906.tif 000005
【図４】
2005148906.tif 000006
【図５】
2005148906.tif 000007
【図６】
2005148906.tif 000008
【図７】
2005148906.tif 000009
【図８】
2005148906.tif 000010
【図９】
2005148906.tif 000011
【図１０】
2005148906.tif 000012
【図１１】
2005148906.tif 000013
【図１２】
2005148906.tif 000014
【図１３】
2005148906.tif 000015
【図１４】
2005148906.tif 000016
【図１５】
2005148906.tif 000017
【図１６】
2005148906.tif 000018
【図１７】
2005148906.tif 000019