JP6184901B2 - Satellite image data processing apparatus, satellite image data processing system, satellite image data processing method and program - Google Patents

Satellite image data processing apparatus, satellite image data processing system, satellite image data processing method and program Download PDF

Info

Publication number
JP6184901B2
JP6184901B2 JP2014098689A JP2014098689A JP6184901B2 JP 6184901 B2 JP6184901 B2 JP 6184901B2 JP 2014098689 A JP2014098689 A JP 2014098689A JP 2014098689 A JP2014098689 A JP 2014098689A JP 6184901 B2 JP6184901 B2 JP 6184901B2
Authority
JP
Japan
Prior art keywords
satellite image
image data
grid
cloud coverage
satellite
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 - Fee Related
Application number
JP2014098689A
Other languages
Japanese (ja)
Other versions
JP2015215252A (en
Inventor
力 草場
力 草場
太一 渥美
太一 渥美
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Solutions Ltd
Original Assignee
Hitachi Solutions Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hitachi Solutions Ltd filed Critical Hitachi Solutions Ltd
Priority to JP2014098689A priority Critical patent/JP6184901B2/en
Publication of JP2015215252A publication Critical patent/JP2015215252A/en
Application granted granted Critical
Publication of JP6184901B2 publication Critical patent/JP6184901B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、大量の衛星画像データを自動的に接合し、広範囲の接合済み衛星画像データを生成するシステムに関する。   The present invention relates to a system for automatically joining a large amount of satellite image data to generate a wide range of joined satellite image data.

光学の地球観測衛星は人の目に見える情報を宇宙空間から観測するため、雲の影響を必ず受ける。地球観測衛星の可視光センサーにより取得された衛星画像データには雲が写り込むことが多く、雲にかかってしまった領域の情報は衛星画像から確認することが出来ない。このような場合、通常は他の日時に撮影された代替の衛星画像を別途用意し、雲に覆われた領域をこの代替の衛星画像から抜き出したデータで置き換える処理が行われている(特許文献1、2)。   Optical earth observation satellites are affected by clouds in order to observe human-visible information from outer space. Clouds often appear in the satellite image data acquired by the visible light sensor of the Earth observation satellite, and information on the area covered by the clouds cannot be confirmed from the satellite image. In such a case, normally, an alternative satellite image captured at another date and time is separately prepared, and a process of replacing an area covered with a cloud with data extracted from the alternative satellite image is performed (Patent Document). 1, 2).

特開2001-143054号公報JP 2001-143054 特開2010-218434号公報JP 2010-218434 A

昨今、地球観測衛星で観測されるデータ量は膨大であり、大量の衛星画像データから広範囲の接合済み衛星画像データを作成することが求められている。しかしながら、上記特許文献1、2に記載の技術では、単一の衛星画像データの修正のみを対象としており、多くの衛星画像候補から最適な衛星画像を抽出して接合し、広範な衛星画像を生成することは考慮されていない。   Recently, the amount of data observed by earth observation satellites is enormous, and it is required to create a wide range of joined satellite image data from a large amount of satellite image data. However, the techniques described in Patent Documents 1 and 2 are intended only for correction of a single satellite image data, and an optimum satellite image is extracted from a large number of satellite image candidates and joined to obtain a wide range of satellite images. It is not considered to generate.

大量の衛星画像データから広範囲の接合済み衛星画像データを作成しようとする場合、目視と手動選択で行われることが多く、作業が繁雑になる上に見落としのリスクがある。   When trying to create a wide range of joined satellite image data from a large amount of satellite image data, it is often performed by visual inspection and manual selection, and the work becomes complicated and there is a risk of oversight.

一方、自動的な処理では、雲を避けられない可能性が高いという問題がある。
本発明は、雲の影響を抑制し衛星画像データの作成精度を向上することを目的とする。
On the other hand, there is a high possibility that clouds cannot be avoided by automatic processing.
An object of the present invention is to suppress the influence of clouds and improve the creation accuracy of satellite image data.

本発明の一観点によれば、撮影条件の異なる複数の衛星画像データを接合して前記複数の衛星画像データの領域を包含する衛星画像データを自動的に生成する衛星画像データ処理装置であって、接合を実施したい衛星画像データと、衛星画像データの接合処理条件である被雲率条件(及び撮影期間)と、衛星画像を複数の領域に区画するグリッドのサイズを示すグリッドパラメータと、を受け取り、前記グリッドパラメータに基づいて、前記衛星画像データの全てを覆うグリッドデータを生成するグリッドデータ生成部と、前記衛星画像データについて、前記グリッドデータにより画定されたグリッド内の被雲率を計算する被雲率計算処理部と、前記被雲率条件を満たすように前記衛星画像データを平面方向に並べ替え、重ね合わせ方向に被雲率が最小になるように前記衛星画像データを重畳することにより前記衛星画像を接合する衛星画像データ接合処理部と、を有することを特徴とする衛星画像データ処理装置が提供される。   According to one aspect of the present invention, there is provided a satellite image data processing device that automatically generates satellite image data including a region of the plurality of satellite image data by joining a plurality of satellite image data having different imaging conditions. Receiving satellite image data to be joined, cloud coverage condition (and photographing period) which is a joining process condition of the satellite image data, and grid parameters indicating the size of the grid dividing the satellite image into a plurality of regions. A grid data generating unit that generates grid data covering all of the satellite image data based on the grid parameters; and a target for calculating a cloud coverage in the grid defined by the grid data for the satellite image data. The satellite image data is rearranged in the plane direction so as to satisfy the cloud coverage ratio and the cloud coverage calculation processing unit, and the coverage in the overlay direction. Incidence and satellite image data combining processing section for joining the satellite images, satellite image data processing apparatus characterized by having a is provided by the superimposed said satellite image data so as to minimize.

衛星画像をより小さいグリッド単位で画像処理して被雲率が最小になるようにすることで、広範囲にわたって被雲率の小さな衛星画像を得ることができるため、陸地等の様子を知る上で雲の影響を低減することができる。   By processing satellite images in smaller grid units so that the cloud coverage is minimized, satellite images with a low cloud coverage can be obtained over a wide range. Can be reduced.

前記衛星画像データの画素値に基づいて海域と陸域とを分類する海域陸域分類処理部を有し、前記衛星画像データ接合処理部は、前記海域陸域分類処理部において分類された陸域の被雲率が最小となるように前記衛星画像データを重畳することにより前記衛星画像を接合することを特徴とする。
陸域における雲の影響を低減することができる。
A land area classification processing unit for classifying a sea area and a land area based on a pixel value of the satellite image data, and the satellite image data joining processing unit is a land area classified in the sea area land classification processing unit The satellite images are joined by superimposing the satellite image data so that the cloud coverage is minimized.
The influence of clouds on land can be reduced.

前記海域陸域分類処理部は、各ピクセル値でNDWIを計算し、衛星画像データにおいてNDWIの値のヒストグラムを作成して陸側と海側とのピークを分けるNDWIの値を海域と陸域とを分類するしきい値として用いることを特徴とする。   The sea area land classification processing unit calculates NDWI for each pixel value, creates a histogram of the NDWI values in the satellite image data, and determines the NDWI values that divide the peak between the land side and the sea side into the sea area and the land area. Is used as a threshold value for classifying.

前記被雲率計算処理部は、計算した被雲率が、(雲の影響が抑制されて陸地が見えるようになる)許容値まで下がったか否かを判定する被雲率判定部と、許容値まで被雲率が下っていないと判定された場合に、前記衛星画像データの接合処理条件の再設定を行う再設定部と、を有し、再設定された前記衛星画像データの接合処理条件を用いて再処理を行うことを特徴とする。
これにより、自動的に良好な画像を得ることができる。
The cloud coverage calculation processing unit is configured to determine whether or not the calculated cloud coverage is reduced to an allowable value (the influence of the cloud is suppressed and the land can be seen), and an allowable value. A re-setting unit that resets the joint processing conditions of the satellite image data when it is determined that the cloud coverage is not lowered to the above, and the joint processing conditions of the reset satellite image data are It is characterized by using and reprocessing.
Thereby, a good image can be obtained automatically.

前記被雲率計算処理部は、計算した被雲率が、(雲の影響が抑制されて陸地が見えるようになる)許容値まで下がったか否かを判定する被雲率判定部と、許容値まで被雲率が下っていないと判定された場合に、前記グリッドパラメータの再設定を行う再設定部と、を有し、再設定された前記グリッドパラメータを用いて再処理を行うことを特徴とする。
これにより、自動的に良好な画像を得ることができる。
The cloud coverage calculation processing unit is configured to determine whether or not the calculated cloud coverage is reduced to an allowable value (the influence of the cloud is suppressed and the land can be seen), and an allowable value. And a resetting unit that resets the grid parameter when it is determined that the cloud coverage is not lowered until the cloud coverage is reduced, and reprocessing is performed using the reset grid parameter. To do.
Thereby, a good image can be obtained automatically.

前記接合処理条件として、前記被雲率条件に代えて又は前記被雲率条件とともに撮影期間を用いることを特徴とする。   As the joining process condition, an imaging period is used instead of the cloud coverage ratio condition or together with the cloud coverage ratio condition.

また、本発明は、撮影条件の異なる複数の衛星画像データを接合して前記複数の衛星画像データの領域を包含する衛星画像データを自動的に生成する衛星画像データ処理装置であって、接合を実施したい衛星画像データと、衛星画像データの接合処理条件である被雲率条件(及び撮影期間)と、衛星画像を複数の領域に区画するグリッドのサイズを示すグリッドパラメータと、を受け取り、前記グリッドパラメータに基づいて、前記衛星画像データの全てを覆うグリッドデータを生成するグリッドデータ生成部と、前記衛星画像データについて、前記グリッドデータにより画定されたグリッド内の被雲率を計算する被雲率計算処理部と、前記被雲率条件を満たすように前記衛星画像データを平面方向に並べ替え、重ね合わせ方向に被雲率が最小になるように前記衛星画像データを重畳することにより前記衛星画像を接合する衛星画像データ接合処理部と、を有することを特徴とする衛星画像データ処理装置である。   The present invention also provides a satellite image data processing device that automatically generates satellite image data including a region of the plurality of satellite image data by combining a plurality of satellite image data having different imaging conditions. Receiving the satellite image data to be executed, the cloud coverage rate condition (and the imaging period) that is a joining process condition of the satellite image data, and the grid parameter indicating the size of the grid that divides the satellite image into a plurality of regions, A grid data generation unit that generates grid data that covers all of the satellite image data based on the parameters, and a cloud coverage calculation that calculates a cloud coverage in the grid defined by the grid data for the satellite image data The satellite image data is rearranged in the plane direction so as to satisfy the cloud coverage condition with the processing unit, and the cloud coverage is maximized in the overlay direction. A satellite image data combining processing section for joining the satellite images by the superimposing satellite image data so that a satellite image data processing apparatus characterized by having a.

前記衛星画像データ接合処理部は、重ね合わせ方向に被雲率が最小になるように前記衛星画像データを重畳することにより前記衛星画像を接合する処理の直前に衛星画像をグリッドの範囲に切り出すことを特徴とする。   The satellite image data joint processing unit cuts out the satellite image into a grid range immediately before the process of joining the satellite images by superimposing the satellite image data so that the cloud coverage is minimized in the overlay direction. It is characterized by.

前記衛星画像データの接合処理条件のうち撮影年月日の条件に合わない衛星画像データを処理対象から予め除外する除外処理部を有することを特徴とする。
これにより、無駄な処理を行わないで済む。
It has an exclusion processing unit for excluding satellite image data that does not meet the shooting date and time conditions among the joint processing conditions of the satellite image data in advance.
This eliminates unnecessary processing.

また、本発明は、撮影条件の異なる複数の衛星画像データを接合して前記複数の衛星画像データの領域を包含する衛星画像データを自動的に生成する衛星画像データ処理方法であって、接合を実施したい衛星画像データと、衛星画像データの接合処理条件である被雲率条件(及び撮影期間)と、衛星画像を複数の領域に区画するグリッドのサイズを示すグリッドパラメータと、を受け取り、前記グリッドパラメータに基づいて、前記衛星画像データの全てを覆うグリッドデータを生成するグリッドデータ生成ステップと、前記衛星画像データについて、前記グリッドデータにより画定されたグリッド内の被雲率を計算する被雲率計算処理ステップと、前記被雲率条件を満たすように前記衛星画像データを平面方向に並べ替え、重ね合わせ方向に被雲率が最小になるように前記衛星画像データを重畳することにより前記衛星画像を接合する衛星画像データ接合処理ステップと、を有することを特徴とする衛星画像データ処理方法である。   The present invention also provides a satellite image data processing method for automatically generating satellite image data including a region of the plurality of satellite image data by joining a plurality of satellite image data having different imaging conditions. Receiving the satellite image data to be executed, the cloud coverage rate condition (and the imaging period) that is a joining process condition of the satellite image data, and the grid parameter indicating the size of the grid that divides the satellite image into a plurality of regions, A grid data generation step for generating grid data covering all of the satellite image data based on the parameters, and a cloud coverage calculation for calculating a cloud coverage in the grid defined by the grid data for the satellite image data Processing step, and rearranging the satellite image data in a plane direction so as to satisfy the cloud coverage ratio, a superposition direction Cloud coverage is a satellite image data processing method characterized by having a satellite image data bonding process step of bonding the satellite image by superimposing the satellite image data so as to minimize.

本発明は、コンピュータに上記に記載の衛星画像データ処理方法を実行させるためのプログラムであっても良い。   The present invention may be a program for causing a computer to execute the satellite image data processing method described above.

また、本発明は、撮影条件の異なる複数の衛星画像データを接合して前記複数の衛星画像データの領域を包含する衛星画像データを自動的に生成する衛星画像データ処理装置であって、接合を実施したい衛星画像データと、衛星画像データの接合処理条件である被雲率条件(及び撮影期間)と、衛星画像を複数の領域に区画するグリッドのサイズを示すグリッドパラメータと、を受け取り、前記グリッドパラメータに基づいて、前記衛星画像データの全てを覆うグリッドデータを生成するグリッドデータ生成部と、前記衛星画像データについて、前記グリッドデータにより画定されたグリッド内の被雲率を計算する被雲率計算処理部と、前記被雲率条件を満たすように前記衛星画像データを平面方向に並べ替え、重ね合わせ方向に被雲率が最小になるように前記衛星画像データを重畳することにより前記衛星画像を接合する衛星画像データ接合処理部と、を有することを特徴とする衛星画像データ処理装置と;
前記衛星画像データの接合処理条件、前記グリッドパラメータを前記衛星画像データ処理装置に送信する端末装置と、を有することを特徴とする衛星画像データ処理システムである。
The present invention also provides a satellite image data processing device that automatically generates satellite image data including a region of the plurality of satellite image data by combining a plurality of satellite image data having different imaging conditions. Receiving the satellite image data to be executed, the cloud coverage rate condition (and the imaging period) that is a joining process condition of the satellite image data, and the grid parameter indicating the size of the grid that divides the satellite image into a plurality of regions, A grid data generation unit that generates grid data that covers all of the satellite image data based on the parameters, and a cloud coverage calculation that calculates a cloud coverage in the grid defined by the grid data for the satellite image data The satellite image data is rearranged in the plane direction so as to satisfy the cloud coverage condition with the processing unit, and the cloud coverage is maximized in the overlay direction. And satellite image data combining processing section for joining the satellite images by the superimposing satellite image data so as to, a satellite image data processing apparatus characterized by having;
A satellite image data processing system comprising: a joint processing condition for the satellite image data; and a terminal device that transmits the grid parameters to the satellite image data processing device.

本発明によれば、作業者の目視による衛星画像データの選択や並び替え作業が少なくなる。また、衛星画像データ選択の誤判断や最適画像の見落としのリスク、作業者の負担を軽減することが可能である。     According to the present invention, the selection and rearrangement work of satellite image data by the operator's visual observation is reduced. In addition, it is possible to reduce the risk of misjudgment in selecting satellite image data, the risk of overlooking the optimum image, and the burden on the operator.

本発明の第1の実施の形態による衛星画像データ処理システムの一構成例を示す機能ブロック図である。It is a functional block diagram which shows the example of 1 structure of the satellite image data processing system by the 1st Embodiment of this invention. 被雲率計算処理部の一構成例を示す機能ブロック図である。It is a functional block diagram which shows the example of 1 structure of a cloud coverage calculation process part. 本発明の第1及び第2の実施の形態による衛星画像データ処理の流れを示すフローチャート図である。It is a flowchart figure which shows the flow of the satellite image data processing by the 1st and 2nd embodiment of this invention. 本発明の第2の実施の形態によるパラメータの自動調整処理の例を示す図である。It is a figure which shows the example of the automatic adjustment process of the parameter by the 2nd Embodiment of this invention. グリッドデータ格納部に格納されているデータの一例を示す図である。It is a figure which shows an example of the data stored in the grid data storage part. 本発明の実施の形態による第1の画像処理の様子を示す図である。It is a figure which shows the mode of the 1st image processing by embodiment of this invention. 本発明の実施の形態による第2の画像処理の様子を示す図である。It is a figure which shows the mode of the 2nd image processing by embodiment of this invention. 本発明の実施の形態による第3の画像処理の様子を示す図である。It is a figure which shows the mode of the 3rd image processing by embodiment of this invention. 本発明の実施の形態による第4の画像処理の様子を示す図である。It is a figure which shows the mode of the 4th image processing by embodiment of this invention. 本発明の実施の形態による第5の画像処理の様子を示す図である。It is a figure which shows the mode of the 5th image processing by embodiment of this invention. 本発明の実施の形態による衛星画像毎に作成したヒストグラムの例を示す図である。It is a figure which shows the example of the histogram produced for every satellite image by embodiment of this invention. 本発明の実施の形態による第6の画像処理の様子を示す図である。It is a figure which shows the mode of the 6th image processing by embodiment of this invention. 本発明の第3の実施の形態による画像処理の様子を示す図である。It is a figure which shows the mode of the image processing by the 3rd Embodiment of this invention. 本発明の第3の実施の形態による衛星画像データ処理の流れを示すフローチャート図である。It is a flowchart figure which shows the flow of the satellite image data process by the 3rd Embodiment of this invention.

以下、本発明の実施の形態について、図面を参照しながら詳細に説明する。   Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.

本実施の形態においては、例えば、多数の衛星画像データを利用する。
図1Aは、本発明の第1の実施形態における衛星画像データ処理システムの一構成例を示す機能ブロック図である。
In the present embodiment, for example, a large number of satellite image data is used.
FIG. 1A is a functional block diagram illustrating a configuration example of a satellite image data processing system according to the first embodiment of the present invention.

図1Aに示すように、衛星画像データ処理システムAは、衛星画像データ処理装置1と、入力装置11と、出力装置である表示装置12と、を有する。これらは一体の装置であっても良い。衛星画像データ処理装置1は、サーバ装置、入力装置は端末装置(例えば携帯端末など)、出力装置も端末装置(例えば、携帯端末など)であっても良い。   As shown in FIG. 1A, the satellite image data processing system A includes a satellite image data processing device 1, an input device 11, and a display device 12 that is an output device. These may be integrated devices. The satellite image data processing device 1 may be a server device, the input device may be a terminal device (for example, a portable terminal), and the output device may be a terminal device (for example, a portable terminal).

衛星画像データ処理装置1は、衛星画像データを格納する衛星画像データ格納部2、衛星画像データを、例えば格子状(方眼状)の小さな領域に区画するグリッドデータを生成するグリッドデータ生成部3(グリッドデータ及びグリッド単位での画像データを格納するグリッドデータ格納部3aを有する)、衛星画像データをグリッドデータに基づいて切り出す切り出し処理部4(処理からある衛星画像データを除外する除外処理部4aを含む)、陸域、海域の分類を行う陸域、海域分類処理部5、衛星画像データ並べ替え処理部6、被雲率を計算する被雲率計算部7、衛星画像データ接合処理部8、接合処理済み衛星画像データ格納部9を有している。   The satellite image data processing device 1 includes a satellite image data storage unit 2 that stores satellite image data, and a grid data generation unit 3 that generates grid data that divides the satellite image data into, for example, small grid (grid-like) regions. A grid data storage unit 3a that stores grid data and image data in units of grids), a cutout processing unit 4 that cuts out satellite image data based on the grid data (an exclusion processing unit 4a that excludes certain satellite image data from the processing) Land area, land area for classifying sea areas, sea area classification processing unit 5, satellite image data rearrangement processing unit 6, cloud cover rate calculating unit 7 for calculating cloud cover rate, satellite image data joining process unit 8, It has a jointed satellite image data storage unit 9.

衛星画像データ格納部2は接合処理する衛星画像データを格納するが、接合処理対象の衛星画像データは2つ以上必要である。グリッドデータ生成部3は衛星画像データの並べ替えや接合処理を実施する領域を表すグリッドを生成する部分であり、システムの利用者が入力装置11からグリッドの大きさを指定することができる。   The satellite image data storage unit 2 stores satellite image data to be joined, but two or more satellite image data to be joined are required. The grid data generation unit 3 is a part that generates a grid representing a region where the satellite image data is rearranged and joined, and the system user can specify the size of the grid from the input device 11.

衛星画像データ切り出し処理部4はグリッドデータの領域で衛星画像データの切り出しを実施する部分である。除外処理部4aは、衛星画像データ切り出し処理部4への入力の前に、入力装置11から入力された撮影年月日などの条件に合わない衛星画像データを処理対象から除外する。   The satellite image data cutout processing unit 4 is a part that cuts out the satellite image data in the grid data area. The exclusion processing unit 4a excludes satellite image data that does not meet the conditions such as the shooting date input from the input device 11 from the processing target before input to the satellite image data cutout processing unit 4.

陸域、海域分類処理部5は、衛星画像データにおける陸域と海域とを分類(区別)する処理をする部分である。例えば、衛星画像の各画素(ピクセル)においてNDWI(正規化水指数)を計算し、しきい値をもって陸域と海域を分類する。   The land area / sea area classification processing unit 5 is a part that performs processing for classifying (distinguishing) land areas and sea areas in the satellite image data. For example, NDWI (normalized water index) is calculated for each pixel (pixel) of the satellite image, and the land area and the sea area are classified based on threshold values.

衛星画像データ並び替え処理部6は、衛星画像データをグリッド領域ごとに、優先順位パラメータに従って並び替え処理をする部分である。被雲率計算処理部7は、陸域、海域分類処理部5で取得した陸域、海域情報を考慮して被雲率を計算し、ユーザーが入力した被雲率を満たす衛星画像データの組み合わせを検索する。   The satellite image data rearrangement processing unit 6 is a part for rearranging the satellite image data according to the priority order parameter for each grid area. The cloud coverage calculation processing unit 7 calculates the cloud coverage in consideration of the land and sea area information acquired by the land and sea classification processing unit 5, and a combination of satellite image data satisfying the cloud coverage input by the user. Search for.

衛星画像データ接合処理部8は、衛星画像データをグリッド領域ごとに接合処理する部分である。接合処理済み衛星画像データ格納部9は、接合処理された衛星画像データを格納する部分である。衛星画像データ処理装置1は、表示装置12と接続可能とされており、衛星画像データを表示することで視覚的に認識することができる。   The satellite image data joining processing unit 8 is a part that joins satellite image data for each grid area. The spliced satellite image data storage unit 9 is a part for storing spliced satellite image data. The satellite image data processing device 1 can be connected to the display device 12 and can be visually recognized by displaying the satellite image data.

図2は、本発明の第1の実施の形態による衛星画像データ処理システムAにおいて、衛星画像データ格納部2に保持されている衛星画像データを自動的に接合処理するための処理の流れの一例を示すフローチャート図である。衛星画像データの自動な画像処理過程は大きく分けて3つであり、衛星画像データをグリッドデータに基づいてグリッド領域に切り出す第1の処理、衛星画像データを優先順位パラメータに従って、グリッド領域単位で検索(選別)する第2の処理、グリッド領域ごとに検索された衛星画像データを接合してより広域の広域衛星画像データを生成する第3の処理を含む。   FIG. 2 shows an example of a processing flow for automatically joining the satellite image data held in the satellite image data storage unit 2 in the satellite image data processing system A according to the first embodiment of the present invention. FIG. There are roughly three automatic image processing processes for satellite image data. The first process is to extract satellite image data into a grid area based on grid data. Search for satellite image data in units of grid areas according to priority parameters. A second process of (selecting) includes a third process of generating satellite image data of a wider area by joining the satellite image data searched for each grid area.

図4は、グリッドデータ格納部3aに格納されているデータの一例を示す図である。図4に示すように、グリッドデータ格納部3aには、グリッド番号と、対象画像管理テーブルNoと、後述する被雲率計算処理(S11でYesの場合)のフラグと、グリッド内の被雲率の計算結果と、が格納されている。このデータにより、どの画像データを利用するか、被雲率計算処理が済んでいるかどうか、被雲率計算処理によるグリッド内の被雲率の計算結果が得られる。これにより、後述する衛星画像の接合処理において、最小の被雲率のグリッドを簡単に得ることができるため、接合するグリッドを迅速に取得することができる。   FIG. 4 is a diagram illustrating an example of data stored in the grid data storage unit 3a. As shown in FIG. 4, the grid data storage unit 3a includes a grid number, a target image management table No, a flag for cloud coverage calculation processing (in the case of Yes in S11) described later, and a cloud coverage in the grid. And the calculation result is stored. With this data, it is possible to obtain which image data is used, whether the cloud coverage calculation processing has been completed, and the calculation result of the cloud coverage in the grid by the cloud coverage calculation processing. Thereby, in the satellite image joining process described later, a grid with the minimum cloud coverage can be easily obtained, so that the grid to be joined can be quickly acquired.

図2に示すように、処理が開始され、ステップS1では、入力装置11から、接合処理を行う衛星画像データ候補、接合処理の際の並び替えの指標となる撮影年月日、被雲率、グリッドデータなどのパラメータを入力する。例えば、「ある年月日以降もしくは以前の衛星画像データを利用する」、「被雲率0.01以下を接合対象とする」、「グリッドのサイズを10m×10mとする」(衛星画像のサイズの1/10などのサイズとする。)といったような条件が入力装置11から指定される。   As shown in FIG. 2, the process is started, and in step S <b> 1, from the input device 11, the candidate satellite image data to be joined, the shooting date and cloud coverage, which are indices for rearrangement during the joining process, Enter parameters such as grid data. For example, “use satellite image data after or before a certain date”, “use cloud coverage of 0.01 or less”, “use grid size of 10 m × 10 m” (satellite image size The input device 11 designates a condition such as 1/10 of the size.

図5において衛星画像データ21a、21b等で示すように、衛星画像は完全に同じ領域が撮影されることは少なく、撮影領域と撮影日とが異なる画像が幾層にも重なっているのが一般的である。衛星画像データ21a、21b等は、その領域(緯度経度に基づく領域の情報)と、撮影年月日などの情報を保持している。従って、緯度経度(位置)がある程度が近い領域の衛星画像データ群を集めてくることで、図5に示すような、衛星画像データ群を全て包含するような所望の広範囲の領域の画像データを得ることができる。   As shown by satellite image data 21a, 21b, etc. in FIG. 5, it is rare that the same area is captured in the satellite image, and images with different shooting areas and shooting dates overlap in many layers. Is. The satellite image data 21a, 21b and the like hold information such as the region (region information based on latitude and longitude) and the shooting date. Therefore, by collecting satellite image data groups in areas where latitude and longitude (position) are close to a certain extent, image data in a desired wide range area including all the satellite image data groups as shown in FIG. Can be obtained.

ここで、雲23aが領域を覆っている割合を示す被雲率のパラメータは各グリッド領域に雲量が含まれる割合である。地表の様子を分析する場合において、海域の雲は陸域の雲に比べて影響が少ないことから、陸域、海域ごとの被雲率を設定することが可能となっており、それぞれの域において雲量の許容度を変更することが出来るようになっている。   Here, the cloud coverage parameter indicating the ratio of the cloud 23a covering the area is the ratio of the cloud amount included in each grid area. When analyzing the state of the surface of the earth, the cloud in the sea area has less influence than the cloud in the land area, so it is possible to set the cloud coverage for each land area and sea area. The tolerance of the cloud cover can be changed.

グリッドデータのパラメータは、例えば、生成するグリッドデータの各グリッド領域の一辺の長さ等で示される。   The parameter of the grid data is indicated by, for example, the length of one side of each grid area of the generated grid data.

これらのパラメータをもとに、衛星画像データの接合の順番を自動的に決定し、接合処理を行なうことができる。   Based on these parameters, it is possible to automatically determine the joining order of the satellite image data and perform the joining process.

ステップS2では、入力されたグリッドデータのパラメータをもとに、グリッドデータ生成部3が、図6に示すように、可能な限り全ての衛星画像候補21a、…を覆うようにグリッドデータ25を生成する。グリッドデータ25はグリッドデータ格納部3aに各像される。   In step S2, based on the input grid data parameters, the grid data generation unit 3 generates grid data 25 so as to cover all the satellite image candidates 21a as much as possible as shown in FIG. To do. The grid data 25 is imaged in the grid data storage unit 3a.

ステップS3では、除外処理部4aが、衛星画像データ候補の中からユーザーにより入力された撮影年月日の条件に合わない衛星画像データを接合処理から除外する。   In step S3, the exclusion processing unit 4a excludes, from the joining process, satellite image data that does not meet the photographing date input by the user from the satellite image data candidates.

ステップS4では、衛星画像データ切り出し処理部4が、ステップS2で生成されたグリッドデータ25を使って除外されなかった全ての衛星画像データを、図7に示すように、各グリッド領域25aで切り出す。   In step S4, the satellite image data cutout processing unit 4 cuts out all satellite image data not excluded using the grid data 25 generated in step S2 in each grid area 25a as shown in FIG.

ステップS5では、図8に示すように、陸域、海域分類処理部5が、衛星画像データについて、陸域と海域の分類を行うために各ピクセル値で正規化植生指数NDWI(Normalized Difference Water Index)を計算する。NDWIにより、海域と陸域とを識別することが可能である。その他の指数を用いても良い。   In step S5, as shown in FIG. 8, the land area / sea area classification processing unit 5 uses the normalized vegetation index NDWI (Normalized Difference Water Index) for each pixel value to classify the land area and the sea area for the satellite image data. ). With NDWI, it is possible to distinguish between sea area and land area. Other indices may be used.

図8の左図は、陸域33と海域35と、その境界にまたがる雲37とが写っている衛星画像データ31の様子を示す図である。ここで、式(1)に示すように、NDWIを衛星画像データ31毎に計算し、その計算結果をピクセル毎に再配置した結果が右図に示される。   The left figure of FIG. 8 is a diagram showing a state of the satellite image data 31 in which the land area 33, the sea area 35, and the clouds 37 straddling the boundary are shown. Here, as shown in Expression (1), the NDWI is calculated for each satellite image data 31, and the result of rearranging the calculation results for each pixel is shown in the right figure.

Figure 0006184901
Figure 0006184901

次いで、ステップS6では、陸域、海域分類処理部5が、図9に示すように、各衛星画像データ31においてNDWIの値のヒストグラムを作成する。図10は、衛星画像毎に作成したヒストグラムの例を示す図である。図10に示すように、横軸にNDWIを、縦軸に度数(画素数)をとると、陸側のヒストグラムは、NDWIが小さい方にピークをもつ度数分布となり、海側のヒストグラムは、NDWIが大きい方にピークをもつ度数分布となる。そして、2つのヒストグラム間の谷の部分を陸域と海域とを分けるしきい値とする。尚、しきい値は衛星画像データ31ごとに算出することが好ましい。   Next, in step S6, the land area / sea area classification processing unit 5 creates a histogram of NDWI values in each satellite image data 31, as shown in FIG. FIG. 10 is a diagram illustrating an example of a histogram created for each satellite image. As shown in FIG. 10, when NDWI is taken on the horizontal axis and frequency (number of pixels) is taken on the vertical axis, the land-side histogram has a frequency distribution with a peak at the smaller NDWI, and the sea-side histogram is NDWI. The frequency distribution has a peak in the larger one. A valley portion between the two histograms is set as a threshold value for dividing the land area and the sea area. The threshold value is preferably calculated for each satellite image data 31.

次いで、ステップS7では、ステップS6で求めたしきい値を用いて、図9の右図のように陸域(黒)33bと海域(白)35bの分類を行う。   Next, in step S7, the land area (black) 33b and the sea area (white) 35b are classified as shown in the right figure of FIG. 9 using the threshold values obtained in step S6.

図1Bは、被雲率を計算する処理を行う被雲率計算処理部7の一構成例を示す機能ブロック図である。被雲率計算処理部7は、被雲率計算部7−1と、被雲率判定部7−2と、再設定部7−3と、計算判定部7−4と、を有する。   FIG. 1B is a functional block diagram illustrating a configuration example of a cloud coverage calculation processing unit 7 that performs a process of calculating the cloud coverage. The cloud coverage calculation processing unit 7 includes a cloud coverage calculation unit 7-1, a cloud coverage determination unit 7-2, a resetting unit 7-3, and a calculation determination unit 7-4.

ステップS8で、衛星画像データ並び替え処理部6は、衛星画像データをグリッド領域ごとに、優先順位パラメータに従って並び替え処理をする。すなわち、衛星画像データの並べ替え処理の候補を全数検索する。   In step S8, the satellite image data rearrangement processing unit 6 rearranges the satellite image data for each grid area according to the priority order parameter. That is, all the candidates for the rearrangement processing of the satellite image data are searched.

ステップS9において、被雲率計算処理部7が、グリッドデータ毎の被雲率を計算する。より詳細には、被雲率計算部7−1が、ステップS1で指定した条件でグリッド内の被雲率を算出する。算出結果の例が図4に示される。ここで、条件とは、例えば、優先順位パラメータとしては、撮影年月日が2013年1月1日以降、地物毎の被雲率が、海側が0.05(5%)、陸側が0.01(1%)などである。   In step S9, the cloud coverage calculation processing unit 7 calculates the cloud coverage for each grid data. More specifically, the cloud coverage calculation unit 7-1 calculates the cloud coverage in the grid under the conditions specified in step S1. An example of the calculation result is shown in FIG. Here, the conditions are, for example, as priority order parameters, the shooting date is January 1, 2013 or later, and the cloud coverage for each feature is 0.05 (5%) on the sea side and 0 on the land side. .01 (1%).

ステップS10において、被雲率判定部7−2が、雲の影響が抑制されて陸地などが見えるようになる許容値まで被雲率(ステップS1における設定値など)が下がったか否かを判定する。その結果、許容値以下まで被雲率が下がったと判定された場合には(Yes)、ステップS11に進み、許容値まで下っていないと判定された場合には(No)、ステップS16において、再設定部7−3がパラメータの再設定等を行い、再設定された条件で、ステップS2に戻って処理を行う。この再設定処理については、第2の実施の形態で詳細に説明する。   In step S10, the cloud coverage determination unit 7-2 determines whether or not the cloud coverage (such as the set value in step S1) has decreased to an allowable value at which the influence of the clouds is suppressed and the land and the like can be seen. . As a result, when it is determined that the cloud coverage rate has decreased to the allowable value or less (Yes), the process proceeds to step S11. When it is determined that the cloud coverage has not decreased to the allowable value (No), in step S16, the process is repeated. The setting unit 7-3 resets parameters and the like, and returns to step S2 to perform processing under the reset conditions. This resetting process will be described in detail in the second embodiment.

次いで、ステップS11において、計算判定部7−4が、全重畳パターン、全グリッドで衛星画像データの並べ替えを検索し被雲率を計算し終わったか否かを判定し、Noであれば、ステップS9に戻り処理を継続し、Yesであれば、ステップS12に進む。   Next, in step S11, the calculation determination unit 7-4 determines whether or not the calculation of the cloud coverage is completed by searching for the rearrangement of the satellite image data in all the superimposition patterns and all the grids. Returning to S9 and continuing the process, if Yes, the process proceeds to Step S12.

被雲率の計算は、図11に示すように、衛星画像を重ね合わせて(例:L1−L6)、グリッド内の被雲率を前通りの衛星画像の並べ替えパターンで計算する(41)。並べ替えると、符号43のようになる。   As shown in FIG. 11, the cloud coverage is calculated by superimposing the satellite images (for example, L1-L6), and calculating the cloud coverage in the grid by the rearrangement pattern of the satellite images (41). . When rearranged, reference numeral 43 is obtained.

次いで、ステップS12において、衛星画像の接合処理が開始され、衛星画像データ接合処理部8がステップS13において、図12に示すように、グリッド領域45は位置が既知であるため、グリッド領域45ごとに被雲率が最も少ない衛星画像データの組み合わせで接合処理を実施することで、最適な接合処理の結果47が得られる。ステップS14では、全グリッド領域で接合処理が完了したか判定し、Noであれば、ステップS13に戻り、YesであればステップS14で、処理を終了する。接合済みの衛星画像データは、位置を含めて互いに関連付けなどをした状態で、例えば、接合処理済み衛星画像データ格納部9などに格納される。接合済みの衛星画像データは、被雲率の少ない衛星画像データであり、これにより、広域に雲の影響が少ない良好な衛星画像を得ることができる。接合済み衛星画像は、ステップS15において、表示装置12に表示させることができる。   Next, in step S12, the satellite image joining process is started. In step S13, the satellite image data joining processing unit 8 has the position of the grid area 45 known as shown in FIG. By performing the joining process using a combination of satellite image data with the least cloud coverage, an optimum joining process result 47 can be obtained. In step S14, it is determined whether the joining process has been completed in all grid areas. If No, the process returns to step S13, and if Yes, the process ends in step S14. The joined satellite image data is stored in, for example, the joined satellite image data storage unit 9 in a state of being associated with each other including the position. The joined satellite image data is satellite image data with a low cloud coverage rate, whereby a good satellite image with little influence of clouds can be obtained over a wide area. The joined satellite image can be displayed on the display device 12 in step S15.

本実施の形態によれば作業者の目視による衛星画像データの選択や並び替え作業がなくなり、衛星画像データ選択の誤判断や最適画像の見落としのリスク、作業者の負担を軽減することができる。   According to the present embodiment, there is no need to select or rearrange satellite image data by the operator's visual observation, and it is possible to reduce the risk of erroneous determination of satellite image data selection, oversight of the optimum image, and the burden on the operator.

尚、海域、陸域に依存する異なる被雲率に基づく計算は、任意であり、海域、陸域に依存しない被雲率に基づいて、グリッド毎に処理を行うようにしても良い。   The calculation based on the different cloud coverage depending on the sea area and the land area is arbitrary, and the processing may be performed for each grid based on the cloud coverage not depending on the sea area and the land area.

以上のように本実施の形態では、大量の衛星画像データの中から撮影年月日や被雲率を優先度としたユーザーにとって最適な広域接合処理済み衛星画像データを自動的に生成することで、ユーザーは生成したい衛星画像データの撮影年月日や被雲率の情報を入力すれば、システムが自動的に衛星画像データの候補から最適な画像を選択し、接合処理を実行するようになっている。   As described above, according to the present embodiment, it is possible to automatically generate satellite image data that has been subjected to wide-area joining processing, which is optimal for a user who prioritizes shooting date and cloud coverage from a large amount of satellite image data. When the user inputs the shooting date and cloud coverage information of the satellite image data that he / she wants to generate, the system automatically selects the optimum image from the candidate satellite image data and executes the joining process. ing.

次に、本発明の第2の実施の形態について説明する。
本実施の形態では、図2のステップS10において、被雲率判定部7−2が、雲の影響が抑制されて陸地などが見えるようになる許容値まで被雲率(ステップS1における設定値など)が下がったか否かを判定した結果、許容値まで被雲率が下っていないと判定された場合には(No)、ステップS16において、再設定部7−3がパラメータの再設定等を行い、ステップS16において再設定された条件で、ステップS2に戻って処理を行う。
Next, a second embodiment of the present invention will be described.
In the present embodiment, in step S10 in FIG. 2, the cloud coverage determination unit 7-2 causes the cloud coverage to reach an allowable value at which the influence of clouds is suppressed and the land and the like can be seen (the set value in step S1 and the like). As a result of determining whether or not the cloud coverage rate has decreased to the allowable value (No), the resetting unit 7-3 performs parameter resetting in step S16. Under the conditions reset in step S16, the process returns to step S2 for processing.

再設定処理に関しては、ユーザーの入力処理による第1の再設定処理と、自動的に再設定値を求める第2の再設定処理がある。   Regarding the resetting process, there are a first resetting process by a user input process and a second resetting process for automatically obtaining a reset value.

第1の再設定処理においては、例えば、ステップS1における優先順位パラメータやグリッドパラメータの設定値において、許容値まで被雲率が下らないことから、例えば、ステップS1における撮影年月日を2013年1月1日まで遡る設定であった場合に、より以前まで遡れば被雲率の条件をクリアできる可能性が高いため、2012年1月1日まで遡る設定とするなどの対応が可能である。   In the first resetting process, for example, since the cloud coverage does not fall to the allowable value in the setting values of the priority order parameter and the grid parameter in step S1, for example, the shooting date in step S1 is set to January 2013. If the setting is traced back to one day, it is highly possible that the condition of the cloud coverage will be cleared if traced back to a previous time. Therefore, it is possible to cope with the setting by going back to January 1, 2012.

図3は、自動的に再設定値を求める第2の再設定処理の例であって、パラメータの再設定値を自動的に求める一例を示す図である。図3(a)は、被雲率のグリッドサイズ依存性を示す図であり、実際の計算を繰り返してプロットした値に基づいて、所望の被雲率が得られるグリッドサイズを求めてその値に再設定する。この場合には、グリッドサイズをより小さい方向にすると被雲率の条件をクリアできる可能性が高いため、クリアできるグリッドサイズを求めてその値を自動的に再設定して処理を行うようにする。   FIG. 3 is a diagram showing an example of the second resetting process for automatically obtaining the reset value, and showing an example of automatically obtaining the parameter reset value. FIG. 3A is a diagram showing the dependence of the cloud coverage on the grid size. Based on the values plotted by repeating the actual calculation, the grid size for obtaining the desired cloud coverage is obtained and set to that value. Reset it. In this case, if the grid size is made smaller, it is highly possible that the cloud coverage condition can be cleared. Therefore, a grid size that can be cleared is obtained, and the value is automatically reset and processed. .

図3(b)は、被雲率の撮影年月日依存性を示す図であり、実際の計算を繰り返してプロットした値に基づいて、所望の被雲率が得られる撮影年月日までの遡及した年月日の値を自動的に求めて再設定して処理を行うようにする。   FIG. 3B is a diagram showing the shooting date dependency of the cloud coverage, and based on the values plotted by repeating actual calculation up to the shooting date when a desired cloud coverage is obtained. The retroactive date is automatically determined and reset to be processed.

以上に説明したように、本実施の形態によれば、所望の被雲率が得られるまで自動的に処理を継続して良好な衛星画像を得ることができるという利点がある。   As described above, according to the present embodiment, there is an advantage that a good satellite image can be obtained by automatically continuing processing until a desired cloud coverage is obtained.

次に、本発明の第3の実施の形態について説明する。
図13は、本実施の形態による衛星画像データ処理の流れを示すフローチャート図である。図13は、図2に対応する図であり、主として図2との相違点について説明する。本実施の形態では、グリッドデータを生成した段階では(ステップS22)、衛星画像をグリッドの範囲で切り出さず、ステップS31において、グリッド内で、陸域と海域の被雲率の条件を満たす重畳パターンを抽出した後に、図12に示すように、ステップS32において、グリッドデータを使って該当する衛星画像をグリッドの範囲で切り出し、ステップS33で、陸域の被雲率が最小のパターンで接合処理を行う。
Next, a third embodiment of the present invention will be described.
FIG. 13 is a flowchart showing the flow of satellite image data processing according to the present embodiment. FIG. 13 is a diagram corresponding to FIG. 2, and differences from FIG. 2 will be mainly described. In the present embodiment, at the stage where the grid data is generated (step S22), the satellite image is not cut out in the range of the grid, and in step S31, the superimposed pattern that satisfies the conditions of the cloud coverage of the land area and the sea area in the grid. As shown in FIG. 12, in step S32, the corresponding satellite image is cut out in the range of the grid in step S32, and in step S33, the joining process is performed with the pattern having the smallest cloud coverage in the land area. Do.

すなわち、第3の実施の形態では、陸域と海域の被雲率の条件を満たすグリッドのみを切り出して陸域の被雲率が最小のパターンでグリッドの接合処理を行う点で、予めグリッドに切り出しておき、陸域と海域の被雲率の条件を満たす重畳パターンを抽出した後に、陸域の被雲率が最小のパターンで衛星画像の接合処理を行う第1の実施の形態とは処理が異なる。第3の実施の形態では、図4のようなグリッド毎の衛星画像の管理は不要である。   That is, in the third embodiment, only grids that satisfy the conditions of the cloud coverage of the land area and the sea area are cut out, and grid joining processing is performed in a pattern with the minimum cloud coverage of the land area. What is the first embodiment that cuts out and extracts a superposition pattern that satisfies the conditions of the cloud coverage in the land and the sea, and then performs a satellite image joining process with a pattern with the minimum cloud coverage in the land? Is different. In the third embodiment, management of satellite images for each grid as shown in FIG. 4 is unnecessary.

尚、第1の実施の形態、第3の実施の形態とも、グリッドの範囲での切り出し処理は、実際の情報処理においては、グリッドの範囲内におけるデータの管理と処理とを意味し、グリッド毎の切り出しグリッド毎の接合を実際に行う必ずしも必要ない。   In both the first embodiment and the third embodiment, the cut-out processing in the grid range means the management and processing of data in the grid range in actual information processing. It is not always necessary to actually perform the joining for each cutting grid.

図1の入力装置11は、利用者が所持する携帯端末(PCやスマートフォンなど)などでも良い。これにより、ユーザーは、任意の位置から各グリッドの大きさを指定することができる。また、処理後の衛星画像も、携帯端末において得ることができることは言うまでもない。   The input device 11 in FIG. 1 may be a mobile terminal (such as a PC or a smartphone) possessed by a user. Thereby, the user can designate the size of each grid from an arbitrary position. Needless to say, the processed satellite image can also be obtained by the mobile terminal.

上記の実施の形態において、添付図面に図示されている構成等については、これらに限定されるものではなく、本発明の効果を発揮する範囲内で適宜変更することが可能である。その他、本発明の目的の範囲を逸脱しない限りにおいて適宜変更して実施することが可能である。ここでは、雲の影響について説明したが、その他の画像への悪影響(下が見えにくい状況)を及ぼす要因を排除するものではない。   In the above-described embodiment, the configuration and the like illustrated in the accompanying drawings are not limited to these, and can be appropriately changed within a range in which the effect of the present invention is exhibited. In addition, various modifications can be made without departing from the scope of the object of the present invention. Although the influence of clouds has been described here, it does not exclude other factors that have an adverse effect on images (a situation where it is difficult to see the bottom).

また、本発明の各構成要素は、任意に取捨選択することができ、取捨選択した構成を具備する発明も本発明に含まれるものである。   Each component of the present invention can be arbitrarily selected, and an invention having a selected configuration is also included in the present invention.

本発明は、衛星画像処理装置に利用可能である。   The present invention is applicable to a satellite image processing apparatus.

A…衛星画像データ処理システム、1…衛星画像データ処理装置、2…衛星画像データ格納部、3…グリッドデータ生成部、4…衛星画像データ切り出し処理部、5…陸域、海域分類処理部、6…衛星画像データ並べ替え処理部、7…衛星画像データ接合処理部、8…接合処理済み衛星画像データ格納部、11…入力装置、12…表示装置(結果表示部)。   A ... Satellite image data processing system, 1 ... Satellite image data processing device, 2 ... Satellite image data storage unit, 3 ... Grid data generation unit, 4 ... Satellite image data cut-out processing unit, 5 ... Land area, sea area classification processing unit, 6 ... Satellite image data rearrangement processing unit, 7 ... Satellite image data joining processing unit, 8 ... Joined satellite image data storage unit, 11 ... Input device, 12 ... Display device (result display unit).

Claims (12)

撮影条件の異なる複数の衛星画像データを接合して前記複数の衛星画像データの領域を包含する衛星画像データを自動的に生成する衛星画像データ処理装置であって、
接合を実施したい衛星画像データと、衛星画像データの接合処理条件である被雲率条件と、衛星画像を複数の領域に区画するグリッドのサイズを示すグリッドパラメータと、を受け取り、前記グリッドパラメータに基づいて、前記衛星画像データの全てを覆うグリッドデータを生成するグリッドデータ生成部と、
前記衛星画像データについて、前記グリッドデータにより画定されたグリッド内の被雲率を計算する被雲率計算処理部と、
前記被雲率条件を満たすように前記衛星画像データを平面方向に並べ替え、重ね合わせ方向に被雲率が最小になるように前記衛星画像データを重畳することにより前記衛星画像を接合する衛星画像データ接合処理部と、
前記衛星画像データの画素値に基づいて海域と陸域とを分類する海域陸域分類処理部と、を有し、
前記衛星画像データ接合処理部は、
前記海域陸域分類処理部において分類された陸域の被雲率が最小となるように前記衛星画像データを重畳することにより前記衛星画像を接合することを特徴とする衛星画像データ処理装置。
A satellite image data processing device that automatically generates satellite image data including a region of the plurality of satellite image data by joining a plurality of satellite image data with different imaging conditions,
Satellite image data that is to be joined, cloud coverage conditions that are joining processing conditions of the satellite image data, and grid parameters that indicate the size of the grid that divides the satellite image into a plurality of regions are received and based on the grid parameters A grid data generating unit for generating grid data covering all of the satellite image data;
With respect to the satellite image data, a cloud coverage calculation processing unit that calculates a cloud coverage in a grid defined by the grid data;
The satellite image data is arranged in the plane direction so as to satisfy the cloud coverage condition, and the satellite image is joined by superimposing the satellite image data so that the cloud coverage is minimized in the overlapping direction. A data junction processing unit;
A sea area land classification processing unit that classifies a sea area and a land area based on a pixel value of the satellite image data,
The satellite image data joining processor is
A satellite image data processing apparatus for joining the satellite images by superimposing the satellite image data so that the cloud coverage of the land area classified by the sea area land area classification processing unit is minimized .
前記海域陸域分類処理部は、
各ピクセル値でNDWIを計算し、衛星画像データにおいてNDWIの値のヒストグラムを作成して陸側と海側とのピークを分けるNDWIの値を海域と陸域とを分類するしきい値として用いることを特徴とする請求項に記載の衛星画像データ処理装置。
The sea area land classification processing unit,
Calculate the NDWI for each pixel value, create a histogram of the NDWI values in the satellite image data, and use the NDWI values that divide the peak between the land side and the sea side as threshold values for classifying the sea area and the land area. The satellite image data processing apparatus according to claim 1 .
撮影条件の異なる複数の衛星画像データを接合して前記複数の衛星画像データの領域を包含する衛星画像データを自動的に生成する衛星画像データ処理装置であって、
接合を実施したい衛星画像データと、衛星画像データの接合処理条件である被雲率条件と、衛星画像を複数の領域に区画するグリッドのサイズを示すグリッドパラメータと、を受け取り、前記グリッドパラメータに基づいて、前記衛星画像データの全てを覆うグリッドデータを生成するグリッドデータ生成部と、
前記衛星画像データについて、前記グリッドデータにより画定されたグリッド内の被雲率を計算する被雲率計算処理部と、
前記被雲率条件を満たすように前記衛星画像データを平面方向に並べ替え、重ね合わせ方向に被雲率が最小になるように前記衛星画像データを重畳することにより前記衛星画像を接合する衛星画像データ接合処理部と、を有し、
前記被雲率計算処理部は、
計算した被雲率が、許容値まで下がったか否かを判定する被雲率判定部と、
許容値まで被雲率が下っていないと判定された場合に、前記衛星画像データの接合処理条件の再設定を行う再設定部と、を有し、
再設定された前記衛星画像データの接合処理条件を用いて再処理を行うことを特徴とする衛星画像データ処理装置。
A satellite image data processing device that automatically generates satellite image data including a region of the plurality of satellite image data by joining a plurality of satellite image data with different imaging conditions,
Satellite image data that is to be joined, cloud coverage conditions that are joining processing conditions of the satellite image data, and grid parameters that indicate the size of the grid that divides the satellite image into a plurality of regions are received and based on the grid parameters A grid data generating unit for generating grid data covering all of the satellite image data;
With respect to the satellite image data, a cloud coverage calculation processing unit that calculates a cloud coverage in a grid defined by the grid data;
The satellite image data is arranged in the plane direction so as to satisfy the cloud coverage condition, and the satellite image is joined by superimposing the satellite image data so that the cloud coverage is minimized in the overlapping direction. A data joining processing unit,
The cloud coverage calculation processing unit,
A cloud coverage determination unit that determines whether or not the calculated cloud coverage is reduced to an allowable value;
When it is determined that the cloud coverage is not lowered to an allowable value, a resetting unit that resets the joint processing conditions of the satellite image data, and
You characterized by reprocessing performed using a bonding process conditions of the satellite image data that has been re-set satellite image data processing apparatus.
前記被雲率計算処理部は、The cloud coverage calculation processing unit,
計算した被雲率が、許容値まで下がったか否かを判定する被雲率判定部と、A cloud coverage determination unit that determines whether or not the calculated cloud coverage is reduced to an allowable value;
許容値まで被雲率が下っていないと判定された場合に、前記衛星画像データの接合処理条件の再設定を行う再設定部と、を有し、When it is determined that the cloud coverage is not lowered to an allowable value, a resetting unit that resets the joint processing conditions of the satellite image data, and
再設定された前記衛星画像データの接合処理条件を用いて再処理を行うことを特徴とする請求項2に記載の衛星画像データ処理装置。3. The satellite image data processing apparatus according to claim 2, wherein reprocessing is performed using the joint processing conditions of the reset satellite image data.
撮影条件の異なる複数の衛星画像データを接合して前記複数の衛星画像データの領域を包含する衛星画像データを自動的に生成する衛星画像データ処理装置であって、
接合を実施したい衛星画像データと、衛星画像データの接合処理条件である被雲率条件と、衛星画像を複数の領域に区画するグリッドのサイズを示すグリッドパラメータと、を受け取り、前記グリッドパラメータに基づいて、前記衛星画像データの全てを覆うグリッドデータを生成するグリッドデータ生成部と、
前記衛星画像データについて、前記グリッドデータにより画定されたグリッド内の被雲率を計算する被雲率計算処理部と、
前記被雲率条件を満たすように前記衛星画像データを平面方向に並べ替え、重ね合わせ方向に被雲率が最小になるように前記衛星画像データを重畳することにより前記衛星画像を接合する衛星画像データ接合処理部と、を有し、
前記被雲率計算処理部は、
計算した被雲率が、許容値まで下がったか否かを判定する被雲率判定部と、
許容値まで被雲率が下っていないと判定された場合に、前記グリッドパラメータの再設定を行う再設定部と、を有し、
再設定された前記グリッドパラメータを用いて再処理を行うことを特徴とする衛星画像データ処理装置。
A satellite image data processing device that automatically generates satellite image data including a region of the plurality of satellite image data by joining a plurality of satellite image data with different imaging conditions,
Satellite image data that is to be joined, cloud coverage conditions that are joining processing conditions of the satellite image data, and grid parameters that indicate the size of the grid that divides the satellite image into a plurality of regions are received and based on the grid parameters A grid data generating unit for generating grid data covering all of the satellite image data;
With respect to the satellite image data, a cloud coverage calculation processing unit that calculates a cloud coverage in a grid defined by the grid data;
The satellite image data is arranged in the plane direction so as to satisfy the cloud coverage condition, and the satellite image is joined by superimposing the satellite image data so that the cloud coverage is minimized in the overlapping direction. A data joining processing unit,
The cloud coverage calculation processing unit,
A cloud coverage determination unit that determines whether or not the calculated cloud coverage is reduced to an allowable value;
A resetting unit that resets the grid parameters when it is determined that the cloud coverage is not reduced to an allowable value,
You characterized by reprocessing performed using the grid parameters is reconfigured satellite image data processing apparatus.
前記被雲率計算処理部は、The cloud coverage calculation processing unit,
計算した被雲率が、許容値まで下がったか否かを判定する被雲率判定部と、A cloud coverage determination unit that determines whether or not the calculated cloud coverage is reduced to an allowable value;
許容値まで被雲率が下っていないと判定された場合に、前記グリッドパラメータの再設定を行う再設定部と、を有し、A resetting unit that resets the grid parameters when it is determined that the cloud coverage is not reduced to an allowable value,
再設定された前記グリッドパラメータを用いて再処理を行うことを特徴とする請求項2又は3に記載の衛星画像データ処理装置。4. The satellite image data processing apparatus according to claim 2, wherein reprocessing is performed using the reset grid parameter.
前記接合処理条件として、前記被雲率条件に代えて又は前記被雲率条件とともに撮影期間を用いることを特徴とする請求項1からまでのいずれか1項に記載の衛星画像データ処理装置。 The satellite image data processing apparatus according to any one of claims 1 to 6 , wherein an imaging period is used as the joining processing condition instead of the cloud coverage ratio condition or together with the cloud coverage ratio condition. 前記衛星画像データ接合処理部は、重ね合わせ方向に被雲率が最小になる条件を満たす重畳パターンを抽出した後に、グリッドデータを使って該当する衛星画像をグリッドの範囲に切り出し、前記衛星画像を接合することを特徴とする請求項1から3までのいずれか1項に記載の衛星画像データ処理装置。 The satellite image data joint processing unit superimposes after extracting satisfy superposition pattern cloud coverage is minimized in the direction to cut out a satellite image corresponding with the grid data within the range of the grid, the satellite The satellite image data processing apparatus according to any one of claims 1 to 3 , wherein the images are joined . 前記衛星画像データの接合処理条件のうち撮影年月日の条件に合わない衛星画像データを処理対象から予め除外する除外処理部を有することを特徴とする請求項1から8までのいずれか1項に記載の衛星画像データ処理装置。   9. The apparatus according to claim 1, further comprising: an exclusion processing unit configured to exclude in advance satellite image data that does not meet the shooting date and time conditions from among the joint processing conditions of the satellite image data. The satellite image data processing apparatus described in 1. 撮影条件の異なる複数の衛星画像データを接合して前記複数の衛星画像データの領域を包含する衛星画像データを自動的に生成する衛星画像データ処理方法であって、
接合を実施したい衛星画像データと、衛星画像データの接合処理条件である被雲率条件と、衛星画像を複数の領域に区画するグリッドのサイズを示すグリッドパラメータと、を受け取り、前記グリッドパラメータに基づいて、前記衛星画像データの全てを覆うグリッドデータを生成するグリッドデータ生成ステップと、
前記衛星画像データについて、前記グリッドデータにより画定されたグリッド内の被雲率を計算する被雲率計算処理ステップと、
前記被雲率条件を満たすように前記衛星画像データを平面方向に並べ替え、重ね合わせ方向に被雲率が最小になるように前記衛星画像データを重畳することにより前記衛星画像を接合する衛星画像データ接合処理ステップと、
前記衛星画像データの画素値に基づいて海域と陸域とを分類する海域陸域分類処理ステップと、を有し、
前記衛星画像データ接合処理ステップは、
前記海域陸域分類処理ステップにおいて分類された陸域の被雲率が最小となるように前記衛星画像データを重畳することにより前記衛星画像を接合することを特徴とする衛星画像データ処理方法。
A satellite image data processing method for automatically generating satellite image data including a region of the plurality of satellite image data by joining a plurality of satellite image data having different imaging conditions,
Satellite image data that is to be joined, cloud coverage conditions that are joining processing conditions of the satellite image data, and grid parameters that indicate the size of the grid that divides the satellite image into a plurality of regions are received and based on the grid parameters Grid data generating step for generating grid data covering all of the satellite image data;
For the satellite image data, a cloud coverage calculation processing step for calculating a cloud coverage in a grid defined by the grid data;
The satellite image data is arranged in the plane direction so as to satisfy the cloud coverage condition, and the satellite image is joined by superimposing the satellite image data so that the cloud coverage is minimized in the overlapping direction. A data joining process step;
A marine land classification process step for classifying a marine area and a land area based on a pixel value of the satellite image data,
The satellite image data joining process step includes:
A satellite image data processing method characterized by joining the satellite images by superimposing the satellite image data so that the cloud coverage of the land area classified in the sea area land area classification processing step is minimized .
コンピュータに請求項10に記載の衛星画像データ処理方法を実行させるためのプログラム。   A program for causing a computer to execute the satellite image data processing method according to claim 10. 撮影条件の異なる複数の衛星画像データを接合して前記複数の衛星画像データの領域を包含する衛星画像データを自動的に生成する衛星画像データ処理装置であって、
接合を実施したい衛星画像データと、衛星画像データの接合処理条件である被雲率条件と、衛星画像を複数の領域に区画するグリッドのサイズを示すグリッドパラメータと、を受け取り、前記グリッドパラメータに基づいて、前記衛星画像データの全てを覆うグリッドデータを生成するグリッドデータ生成部と、
前記衛星画像データについて、前記グリッドデータにより画定されたグリッド内の被雲率を計算する被雲率計算処理部と、
前記被雲率条件を満たすように前記衛星画像データを平面方向に並べ替え、重ね合わせ方向に被雲率が最小になるように前記衛星画像データを重畳することにより前記衛星画像を接合する衛星画像データ接合処理部と、
前記衛星画像データの画素値に基づいて海域と陸域とを分類する海域陸域分類処理部を有し、
前記衛星画像データ接合処理部は、
前記海域陸域分類処理部において分類された陸域の被雲率が最小となるように前記衛星画像データを重畳することにより前記衛星画像を接合することを特徴とする衛星画像データ処理装置と;
前記衛星画像データの接合処理条件、前記グリッドパラメータを前記衛星画像データ処理装置に送信する端末装置と、
を有することを特徴とする衛星画像データ処理システム。
A satellite image data processing device that automatically generates satellite image data including a region of the plurality of satellite image data by joining a plurality of satellite image data with different imaging conditions,
Satellite image data that is to be joined, cloud coverage conditions that are joining processing conditions of the satellite image data, and grid parameters that indicate the size of the grid that divides the satellite image into a plurality of regions are received and based on the grid parameters A grid data generating unit for generating grid data covering all of the satellite image data;
With respect to the satellite image data, a cloud coverage calculation processing unit that calculates a cloud coverage in a grid defined by the grid data;
The satellite image data is arranged in the plane direction so as to satisfy the cloud coverage condition, and the satellite image is joined by superimposing the satellite image data so that the cloud coverage is minimized in the overlapping direction. A data junction processing unit;
A sea area land classification processing unit that classifies the sea area and the land area based on the pixel value of the satellite image data,
The satellite image data joining processor is
A satellite image data processing apparatus that joins the satellite images by superimposing the satellite image data so that the cloud coverage of the land areas classified by the sea area land area classification processing unit is minimized ;
A joint processing condition for the satellite image data, and a terminal device for transmitting the grid parameter to the satellite image data processing device;
A satellite image data processing system.
JP2014098689A 2014-05-12 2014-05-12 Satellite image data processing apparatus, satellite image data processing system, satellite image data processing method and program Expired - Fee Related JP6184901B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2014098689A JP6184901B2 (en) 2014-05-12 2014-05-12 Satellite image data processing apparatus, satellite image data processing system, satellite image data processing method and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2014098689A JP6184901B2 (en) 2014-05-12 2014-05-12 Satellite image data processing apparatus, satellite image data processing system, satellite image data processing method and program

Publications (2)

Publication Number Publication Date
JP2015215252A JP2015215252A (en) 2015-12-03
JP6184901B2 true JP6184901B2 (en) 2017-08-23

Family

ID=54752265

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014098689A Expired - Fee Related JP6184901B2 (en) 2014-05-12 2014-05-12 Satellite image data processing apparatus, satellite image data processing system, satellite image data processing method and program

Country Status (1)

Country Link
JP (1) JP6184901B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI726396B (en) * 2019-08-23 2021-05-01 經緯航太科技股份有限公司 Environmental inspection system and method

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2019028560A (en) * 2017-07-26 2019-02-21 エスゼット ディージェイアイ テクノロジー カンパニー リミテッドSz Dji Technology Co.,Ltd Mobile platform, image composition method, program and recording medium
KR102632862B1 (en) * 2023-09-12 2024-02-05 주식회사 퍼픽셀 System and method for dynamic periodic composition of satellite images

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001143054A (en) * 1999-11-16 2001-05-25 Hitachi Ltd Satellite image-processing method
JP3961432B2 (en) * 2003-02-04 2007-08-22 日立ソフトウエアエンジニアリング株式会社 Image processing program
JP4145259B2 (en) * 2004-02-26 2008-09-03 富士通株式会社 Program, satellite adjustment device, control method
JP4788447B2 (en) * 2006-04-05 2011-10-05 株式会社日立製作所 Image processing program and image processing apparatus
JP4365887B1 (en) * 2009-03-18 2009-11-18 株式会社パスコ Method and apparatus for generating ground surface image data
WO2013111552A1 (en) * 2012-01-27 2013-08-01 パナソニック株式会社 Image processor, imaging device, and image processing method

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI726396B (en) * 2019-08-23 2021-05-01 經緯航太科技股份有限公司 Environmental inspection system and method

Also Published As

Publication number Publication date
JP2015215252A (en) 2015-12-03

Similar Documents

Publication Publication Date Title
US11144758B2 (en) System and method for object detection and classification in aerial imagery
JP6510247B2 (en) Survey data processing apparatus, survey data processing method and program
US20200110847A1 (en) System and Method for Construction Estimation Using Aerial Images
US10032267B2 (en) Automating the assessment of damage to infrastructure assets
Perko et al. Advanced DTM generation from very high resolution satellite stereo images
Castillo et al. SF3M software: 3-D photo-reconstruction for non-expert users and its application to a gully network
JP5535025B2 (en) Outdoor feature detection system, program for outdoor feature detection system, and recording medium for program for outdoor feature detection system
KR101510206B1 (en) Urban Change Detection Method Using the Aerial Hyper Spectral images for Digital Map modify Drawing
JP5542530B2 (en) Sampling position determination device
KR20070007790A (en) Forming a single image from overlapping images
JP6184901B2 (en) Satellite image data processing apparatus, satellite image data processing system, satellite image data processing method and program
US11543080B2 (en) Automated pipeline construction modelling
JP2010128565A (en) Building change detecting device, building change detecting method, and program
JP2007188117A (en) Cave-in area extraction method, device and program
JP2017162179A5 (en)
JPH10269347A (en) Method for eliminating shade shadow element on geographic image, geographic image processor and recording medium
Tahar et al. Capabilities of UAV-Based Watershed Segmentation Method for Estimating Tree Crown: A Case Study of Oil Palm Tree
KR101510203B1 (en) Land and City Monitoring Method Using Hyper Spectral Images
JP2012243226A (en) Ground height data generation device, ground height data generation method, and program
KR20170108552A (en) Information system for analysis of waterfront structure damage
JP2020098442A (en) Tree width setting device and program
Ali-Sisto et al. Comparison of 3D point clouds from aerial stereo images and lidar for forest change detection
JP2020144686A (en) Model creating system, and model creating method
Castillo et al. The SF3M approach to 3-D photo-reconstruction for non-expert users: application to a gully network
Zeng et al. A stereo image matching method to improve the DSM accuracy inside building boundaries

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160531

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170303

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170321

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170413

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20170711

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170726

R150 Certificate of patent or registration of utility model

Ref document number: 6184901

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees