Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7520649B2 - Weather information processing device, weather information processing system, weather information processing method, and weather information processing program - Google Patents
[go: Go Back, main page]

JP7520649B2 - Weather information processing device, weather information processing system, weather information processing method, and weather information processing program - Google Patents

Weather information processing device, weather information processing system, weather information processing method, and weather information processing program Download PDF

Info

Publication number
JP7520649B2
JP7520649B2 JP2020147794A JP2020147794A JP7520649B2 JP 7520649 B2 JP7520649 B2 JP 7520649B2 JP 2020147794 A JP2020147794 A JP 2020147794A JP 2020147794 A JP2020147794 A JP 2020147794A JP 7520649 B2 JP7520649 B2 JP 7520649B2
Authority
JP
Japan
Prior art keywords
cloud
image
celestial
weather information
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.)
Active
Application number
JP2020147794A
Other languages
Japanese (ja)
Other versions
JP2022042380A (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.)
Central Research Institute of Electric Power Industry
Sky Perfect Jsat Corp
Original Assignee
Central Research Institute of Electric Power Industry
Sky Perfect Jsat Corp
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 Central Research Institute of Electric Power Industry, Sky Perfect Jsat Corp filed Critical Central Research Institute of Electric Power Industry
Priority to JP2020147794A priority Critical patent/JP7520649B2/en
Publication of JP2022042380A publication Critical patent/JP2022042380A/en
Application granted granted Critical
Publication of JP7520649B2 publication Critical patent/JP7520649B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Processing (AREA)

Description

本発明は、気象情報処理装置、気象情報処理システム、気象情報処理方法及び気象情報処理プログラムに関する。 The present invention relates to a weather information processing device, a weather information processing system, a weather information processing method, and a weather information processing program.

気象観測は、静止気象衛星、レーダー及び各地の観測所等を用いて行われる。そして、気象観測により得られた各種の気象データは、一般に公開される等して、事業者や個人を問わず広く利用される。 Weather observations are conducted using geostationary meteorological satellites, radar, and observation stations around the country. Various types of meteorological data obtained through weather observations are made publicly available and are widely used by businesses and individuals alike.

近年、気象データは様々な分野での活用が進んでいる。例えば、得られた気象データを用いて将来の気象を予測することで、天気予報、災害予報、河川管理、登山者向けの情報提供及び農業生産管理等を行うことが可能である。また、気象データにより雲の動き等を予測することで、太陽光発電の出力予測等を行うことも可能となる。さらには、気象データから得られる太陽光発電の出力予測等を用いて、発電所からの電力供給量の調整等といった運用管理に用いることも可能である。また、雲の状況を把握することで、気象データは航空機の運航管理等にも用いられる。 In recent years, meteorological data has been increasingly used in a variety of fields. For example, by using the obtained meteorological data to predict future weather, it is possible to carry out weather forecasts, disaster forecasts, river management, providing information to hikers, and agricultural production management. In addition, by predicting cloud movements using meteorological data, it is possible to predict the output of solar power generation. Furthermore, it is possible to use the output predictions of solar power generation obtained from meteorological data for operational management such as adjusting the amount of power supplied from power plants. In addition, by understanding the cloud conditions, meteorological data can also be used for aircraft operation management, etc.

特に、気象データにおける雲の情報は、防災や太陽光発電等において重要な役割を担う情報である。気象庁は、気象衛星による観測データや数値予報データを用いて求めた雲頂高度、雲型及び雲の品質情報等の雲情報の提供を行っている。また、気象衛星画像を用いて気象学的な解釈に基づいて雲画像に含まれる各種の気象情報を抽出する技術がある。この技術では、例えば、可視(VIS:Visible)画像を使用することで、雲の輝度から雲を識別したり雲の厚みを求めたりできる。また、赤外(IR:InfraRed)画像を用いることで、雲頂温度及び雲頂高度を求めることができ、且つ、雲の識別も行うことができる(特許文献1)。 Cloud information in meteorological data in particular plays an important role in disaster prevention, solar power generation, and the like. The Japan Meteorological Agency provides cloud information such as cloud top altitude, cloud type, and cloud quality information obtained using observation data and numerical forecast data from meteorological satellites. There is also technology that uses meteorological satellite images to extract various weather information contained in cloud images based on meteorological interpretation. With this technology, for example, visible (VIS) images can be used to identify clouds and determine cloud thickness from their brightness. Infrared (IR) images can also be used to determine cloud top temperature and cloud top altitude, and clouds can also be identified (Patent Document 1).

また、雲の厚さに関しては、雲底高度を測定する雲底高度計(シーロメータ)を用いて雲量判定を行なう雲量判定装置が提案されている。雲底高度計は、鉛直上方に向けてレーザー光を発射し、発射された該レーザー光の雲による反射光を受光素子により受光し、発射から受光までに要した時間に基づいて雲の最下部の高度である雲底高度を求める装置である。 As for cloud thickness, a cloud amount determination device has been proposed that uses a cloud base altimeter (ceilometer) to determine cloud amount. A cloud base altimeter emits a laser beam vertically upward, receives the reflected light of the emitted laser beam from clouds using a light receiving element, and determines the cloud base height, which is the height of the bottom of the cloud, based on the time required from emission to reception.

このような、気象データを取り扱う技術として、気象衛星画像の可視画像データ(VIS)から取得した反射情報及び赤外画像データから取得した輝度温度情報を用いて雲の種類を判定し、日射量を推定する従来技術がある。また、2台のカメラで撮影した天頂を中央とする広角画像から雲の経度緯度を求めて雲の移動速度や風速を求める従来技術がある(特許文献2)。また、衛星画像データにおける雲やエアロゾルの影響を表すパラメータを調整して、調整後の衛星画像データから日射量を推定する従来技術がある(特許文献3)。他にも、日射強度から雲種を判定し、雲種から日射強度を補正して日射強度の推定値を算出する従来技術がある(特許文献4)。 As a technology for handling such meteorological data, there is a conventional technology that uses reflectance information obtained from visible image data (VIS) of meteorological satellite images and brightness temperature information obtained from infrared image data to determine the type of cloud and estimate the amount of solar radiation. There is also a conventional technology that obtains the longitude and latitude of clouds from a wide-angle image centered on the zenith captured by two cameras, and determines the cloud movement speed and wind speed (Patent Document 2). There is also a conventional technology that adjusts parameters that represent the effects of clouds and aerosols in satellite image data, and estimates the amount of solar radiation from the adjusted satellite image data (Patent Document 3). There is also a conventional technology that determines the cloud type from solar radiation intensity, corrects the solar radiation intensity from the cloud type, and calculates an estimate of the solar radiation intensity (Patent Document 4).

特開2017-90442号公報JP 2017-90442 A 特開2019-60754号公報JP 2019-60754 A 特開平11-211560号公報Japanese Patent Application Publication No. 11-211560 特開2015-137903号公報JP 2015-137903 A

しかしながら、衛星画像データを用いて雲の様相(雲の規模、雲底高度、雲頂高度、雲の厚み)を求める場合、例えば輝度から雲の厚みを予測することになり、雲の高度等による影響を受けるため正確な厚みを求めることは困難である。また、雲底高度計を用いることで雲底高度測定を行なうことはできるが、雲底高度計は高価であり、且つ、鉛直上の雲底の計測に制限されるため設置性が悪い。そのため、設置場所によっては設置することが困難であり、また、多数の装置を設置した場合には多額のコストがかかることになり、広い範囲で雲の厚みを測定することが困難である。したがって、雲の様相を正確に同定することは困難である。 However, when using satellite image data to determine cloud appearance (cloud scale, cloud base altitude, cloud top altitude, cloud thickness), the cloud thickness is predicted from brightness, for example, and it is difficult to determine the exact thickness because it is affected by the cloud altitude, etc. Also, while it is possible to measure cloud base altitude using a cloud base altimeter, cloud base altimeters are expensive and difficult to install because they are limited to measuring the cloud base vertically. As a result, they can be difficult to install in some locations, and installing multiple devices incurs a large cost, making it difficult to measure cloud thickness over a wide area. Therefore, it is difficult to accurately identify the cloud appearance.

また、気象衛星画像から雲の種類を判定する従来技術では、雲の正確な厚みを求めることは困難である。これは、天頂を中心とする広角画像から雲の移動速度や風速を求める従来技術や、雲やエアロゾルの影響を表すパラメータを調整して日射量を推定する従来技術や、日射強度から雲種を判定する従来技術であっても同様である。したがって、雲の様相を正確に同定することは困難であり、気象データを用いるシステムの性能を向上させることは難しい。 Furthermore, with conventional technology that determines the type of cloud from meteorological satellite images, it is difficult to determine the exact thickness of clouds. This is true for conventional technology that determines the cloud movement speed and wind speed from wide-angle images centered on the zenith, conventional technology that estimates solar radiation by adjusting parameters that represent the effects of clouds and aerosols, and conventional technology that determines cloud types from solar radiation intensity. Therefore, it is difficult to accurately identify the state of clouds, and it is difficult to improve the performance of systems that use meteorological data.

開示の技術は、上記に鑑みてなされたものであって、気象データを活用するシステムの性能を向上させる気象情報処理装置、気象情報処理システム、気象情報処理方法及び気象情報処理プログラムを提供することを目的とする。 The disclosed technology has been developed in consideration of the above, and aims to provide a weather information processing device, a weather information processing system, a weather information processing method, and a weather information processing program that improve the performance of a system that utilizes weather data.

本願の開示する気象情報処理装置、気象情報処理システム、気象情報処理方法及び気象情報処理プログラムの一つの態様において、画像照合部は、天頂方向を中心とした平面の天球画像を撮影するカメラを有する地上センサ及び気象衛星で撮影された衛星画像を配信する気象データ配信システムから受信した、前記天球画像及び前記衛星画像を基に、前記天球画像に映った雲分布と前記衛星画像に映った雲分布との同定を行って画像照合を実行する。範囲取得部は、前記画像照合の結果を基に、前記衛星画像の水平方向の距離情報から前記天球画像における水平方向の距離情報を取得する。雲底高度算出部は、前記天球画像における前記水平方向の距離情報及び前記天球画像における雲の天頂角を基に前記雲の雲底高度を算出する。 In one aspect of the weather information processing device, weather information processing system, weather information processing method, and weather information processing program disclosed in the present application, the image matching unit performs image matching by identifying the cloud distribution captured in the celestial image and the cloud distribution captured in the satellite image based on the celestial image and the satellite image received from a ground sensor having a camera that captures a planar celestial image centered in the zenith direction and a weather data distribution system that distributes satellite images captured by a meteorological satellite. The range acquisition unit acquires horizontal distance information in the celestial image from horizontal distance information in the satellite image based on the result of the image matching. The cloud base altitude calculation unit calculates the cloud base altitude of the cloud based on the horizontal distance information in the celestial image and the zenith angle of the cloud in the celestial image.

1つの側面では、本発明は、気象データを活用するシステムの性能を向上させることができる。 In one aspect, the present invention can improve the performance of systems that utilize weather data.

図1は、実施例1に係る気象情報提供システムのブロック図である。FIG. 1 is a block diagram of a weather information providing system according to a first embodiment. 図2は、天球カメラ及び気象情報処理装置の詳細を表すブロック図である。FIG. 2 is a block diagram showing the details of the celestial camera and the weather information processing device. 図3は、天球画像の平面画像への変換の一例を示す図である。FIG. 3 is a diagram showing an example of converting a celestial spherical image into a planar image. 図4は、天球カメラにおける画像取得の概略を説明するための図である。FIG. 4 is a diagram for explaining an outline of image acquisition by the celestial camera. 図5は、気象衛星画像の一例を表す図である。FIG. 5 is a diagram showing an example of a meteorological satellite image. 図6は、雲種の判定テーブルの一例を示す図である。FIG. 6 is a diagram showing an example of a cloud type determination table. 図7は、画像照合解析処理の一例を示す図である。FIG. 7 is a diagram showing an example of the image matching and analysis process. 図8は、雲のマッチング処理を説明するための図である。FIG. 8 is a diagram for explaining the cloud matching process. 図9は、天頂角と方位角とを説明するための図である。FIG. 9 is a diagram for explaining the zenith angle and the azimuth angle. 図10は、雲底高度の算出の一例を示す図である。FIG. 10 is a diagram showing an example of calculation of the cloud base altitude. 図11は、雲の厚さと雲頂高度及び雲底高度との関係を示す図である。FIG. 11 is a diagram showing the relationship between cloud thickness and the cloud top and cloud base altitudes. 図12は、雲の様相を含む気象情報提供のフローチャートである。FIG. 12 is a flow chart for providing weather information including cloud appearance. 図13は、気象情報処理装置による雲の様相の識別処理のフローチャートである。FIG. 13 is a flowchart showing the process of identifying cloud appearances by the weather information processing device. 図14は、高度の違う雲の平面画像における位置関係を説明するための図である。FIG. 14 is a diagram for explaining the positional relationship in a planar image of clouds at different altitudes. 図15は、実施例1の変形例に係る気象情報処理装置による高さの異なる雲のそれぞれの雲底高度の算出を説明するための図である。FIG. 15 is a diagram for explaining the calculation of the cloud base altitudes of clouds with different heights by the weather information processing device according to the modified example of the first embodiment. 図16は、実施例2に係る気象情報処理装置のブロック図である。FIG. 16 is a block diagram of a weather information processing device according to the second embodiment. 図17は、移動方向及び移動速度の算出の一例を説明するための図である。FIG. 17 is a diagram for explaining an example of calculation of the movement direction and the movement speed. 図18は、実施例2に係る気象情報処理装置による雲及び雲域の移動方向及び移動速度の算出処理のフローチャートである。FIG. 18 is a flowchart of a calculation process of the moving direction and moving speed of clouds and cloud areas by the weather information processing device according to the second embodiment. 図19は、気象情報処理装置のハードウェア構成図である。FIG. 19 is a hardware configuration diagram of the weather information processing device.

以下に、本願の開示する気象情報処理装置、気象情報処理システム、気象情報処理方法及び気象情報処理プログラムの実施例を図面に基づいて詳細に説明する。なお、以下の実施例により本願の開示する気象情報処理装置、気象情報処理システム、気象情報処理方法及び気象情報処理プログラムが限定されるものではない。 Below, examples of the weather information processing device, weather information processing system, weather information processing method, and weather information processing program disclosed in the present application will be described in detail with reference to the drawings. Note that the weather information processing device, weather information processing system, weather information processing method, and weather information processing program disclosed in the present application are not limited to the following examples.

図1は、実施例1に係る気象情報提供システムのブロック図である。また、図2は、天球カメラ及び気象情報処理装置の詳細を表すブロック図である。気象情報提供システム100は、図1に示すように、気象情報処理装置1、地上センサ2、気象情報データベース3及び情報活用基盤4を有する。 Figure 1 is a block diagram of a weather information providing system according to the first embodiment. Also, Figure 2 is a block diagram showing the details of the celestial camera and the weather information processing device. As shown in Figure 1, the weather information providing system 100 has a weather information processing device 1, a ground sensor 2, a weather information database 3, and an information utilization platform 4.

地上センサ2は、気象情報処理装置1に接続される。地上センサ2は、図2に示すように、天球カメラ21及び気象データ取得部22を有する。また、地上センサ2は、自装置の設置場所の緯度及び経度と時刻とを測定する測定機能を有する。 The ground sensor 2 is connected to the weather information processing device 1. As shown in FIG. 2, the ground sensor 2 has a celestial camera 21 and a weather data acquisition unit 22. The ground sensor 2 also has a measurement function for measuring the latitude and longitude of the location where the ground sensor 2 is installed and the time.

地上センサ2には、風が通り抜ける個所が設けられる。気象データ取得部22は、地上センサ2を通り抜ける風を用いて、地上の気温、気圧及び相対湿度を含む地上気象データを測定する。 The ground sensor 2 is provided with a portion through which wind passes. The weather data acquisition unit 22 uses the wind passing through the ground sensor 2 to measure ground weather data including ground temperature, air pressure, and relative humidity.

地上センサ2は、天球カメラ21により撮影された天球画像201及び気象データ取得部22による測定結果を気象情報処理装置1へ送信する。この際、地上センサ2は、自装置の設置場所の緯度及び経度、並びに、撮影時刻も気象情報処理装置1へ送信する。 The ground sensor 2 transmits the celestial image 201 captured by the celestial camera 21 and the measurement results by the meteorological data acquisition unit 22 to the weather information processing device 1. At this time, the ground sensor 2 also transmits the latitude and longitude of the location where the ground sensor 2 is installed, and the time of capture to the weather information processing device 1.

天球カメラ21は、予め設定された時間毎に、可視光及び赤外線で天頂方向を中心とした半球面画像である図3に示す天球画像201を撮影する。図3は、天球画像の平面画像への変換の一例を示す図である。ここで、天球画像201について具体的に説明する。 The celestial camera 21 captures the celestial image 201 shown in FIG. 3, which is a hemispherical image centered in the zenith direction, using visible light and infrared light at preset time intervals. FIG. 3 is a diagram showing an example of converting the celestial image into a planar image. Here, the celestial image 201 will be described in detail.

天球カメラ21は、図4に示すイメージセンサ23及び魚眼レンズ24を有する。図4は、天球カメラにおける画像取得の概略を説明するための図である。魚眼レンズ24は、半球の形をしており、天頂を中心にした円周の全周にわたって例えば画角180度の光を取り込む。魚眼レンズ24に取り込まれた光は屈折される入射ルートを経由してイメージセンサ23により受光される。これにより、イメージセンサ23は、天頂を中心にした円周の全周にわたって画角180度の画像を形成する。 The celestial camera 21 has an image sensor 23 and a fisheye lens 24 as shown in FIG. 4. FIG. 4 is a diagram for explaining an outline of image acquisition in the celestial camera. The fisheye lens 24 is hemispherical in shape and captures light with a field angle of, for example, 180 degrees all around the circumference centered on the zenith. The light captured by the fisheye lens 24 is received by the image sensor 23 via an incident route where it is refracted. As a result, the image sensor 23 forms an image with a field angle of 180 degrees all around the circumference centered on the zenith.

ここで、天球画像201は、地上から撮影された画像であるため、図3に示すように東西が反転された画像となる。天球画像201には、雲底部分が映される。 Here, the celestial image 201 is an image taken from the ground, so the image is inverted east-west as shown in FIG. 3. The celestial image 201 shows the cloud base.

図1に戻って説明を続ける。気象データ配信システム5は、ひまわり8号等の気象衛星6により撮影された気象衛星画像を取得する。そして、気象データ配信システム5は、取得した気象衛星画像を蓄積する。気象衛星画像には、可視画像、赤外画像及び近赤外画像が含まれる。 Returning to Figure 1, we will continue the explanation. The weather data distribution system 5 acquires meteorological satellite images taken by meteorological satellites 6 such as Himawari-8. The meteorological data distribution system 5 then stores the acquired meteorological satellite images. The meteorological satellite images include visible images, infrared images, and near-infrared images.

気象情報処理装置1は、地上センサ2から取得した天球画像201及び地上気象データ、並びに、気象データ配信システム5から取得した衛星画像を用いて、天球画像201に映った雲と衛星画像に映った雲との同定を行う。そして、気象情報処理装置1は、同定し天球画像から雲の様相の情報を取得する。その後、気象情報処理装置1は、取得した雲の様相の情報を気象情報データベース3に格納する。気象情報処理装置1は、地上センサ2、気象情報データベース3及び気象データ配信システム5が接続される。以下に、図2を参照して気象情報処理装置1の詳細について説明する。 The weather information processing device 1 uses the celestial image 201 and ground weather data acquired from the ground sensor 2, and the satellite image acquired from the weather data distribution system 5 to identify clouds captured in the celestial image 201 and clouds captured in the satellite image. The weather information processing device 1 then acquires information on the appearance of clouds from the identified celestial image. The weather information processing device 1 then stores the acquired information on the appearance of clouds in the weather information database 3. The weather information processing device 1 is connected to the ground sensor 2, the weather information database 3, and the weather data distribution system 5. The weather information processing device 1 is described in detail below with reference to FIG. 2.

図2に示すように、本実施例に係る気象情報処理装置1は、画像解析部101、変換部102、雲頂高度算出部103、画像照合部104、範囲取得部105、雲厚算出部106、雲底高度算出部107及び様相識別部108を有する。 As shown in FIG. 2, the weather information processing device 1 according to this embodiment has an image analysis unit 101, a conversion unit 102, a cloud top height calculation unit 103, an image matching unit 104, a range acquisition unit 105, a cloud thickness calculation unit 106, a cloud base height calculation unit 107, and an appearance identification unit 108.

画像解析部101は、天球カメラ21で撮影された天球画像201を地上センサ2から受信する。この時、画像解析部101は、地上センサ2の設置場所の緯度及び経度、撮影時刻及び地上気象データも地上センサ2から受信する。また、画像解析部101は、天球画像203の撮影時刻の通知を地上センサ2から受ける。そして、画像解析部101は、通知された撮影時刻の気象衛星画像を気象データ配信システム5から取得する。画像解析部101が取得する気象衛星画像には、図5に示すような可視画像301及び赤外画像302が含まれる。図5は、気象衛星画像の一例を表す図である。 The image analysis unit 101 receives the celestial image 201 captured by the celestial camera 21 from the ground sensor 2. At this time, the image analysis unit 101 also receives the latitude and longitude of the installation location of the ground sensor 2, the capture time, and ground meteorological data from the ground sensor 2. The image analysis unit 101 also receives notification of the capture time of the celestial image 203 from the ground sensor 2. The image analysis unit 101 then acquires the meteorological satellite image of the notified capture time from the meteorological data distribution system 5. The meteorological satellite image acquired by the image analysis unit 101 includes a visible image 301 and an infrared image 302 as shown in FIG. 5. FIG. 5 is a diagram showing an example of a meteorological satellite image.

ここで、図5を参照して気象衛星画像について説明する。例えば、気象衛星6として気象衛星ひまわり8号を用いる場合について説明する。その場合の気象衛星6が可視画像301を撮影する場合の衛星直下点における空中分解能は、0.5~1.0kmである。また、気象衛星6が赤外画像302を撮影する場合の空中分解能は、1~2kmである。気象衛星6が近赤外画像を撮影する場合の空中分解能も同様に2kmである。また、気象衛星6は、可視画像301を生成する観測バンドを3個、近赤外の画像を生成するバンドを3個、赤外画像302を生成する観測バンドを10個有する。そして、気象衛星6は、2.5分の観測時間間隔で、可視画像301及び赤外画像302を取得する。 Now, referring to FIG. 5, we will explain meteorological satellite images. For example, we will explain the case where the meteorological satellite Himawari-8 is used as the meteorological satellite 6. In this case, the aerial resolution at the satellite's nadir when the meteorological satellite 6 captures a visible image 301 is 0.5 to 1.0 km. When the meteorological satellite 6 captures an infrared image 302, the aerial resolution is 1 to 2 km. When the meteorological satellite 6 captures a near-infrared image, the aerial resolution is also 2 km. The meteorological satellite 6 has three observation bands that generate the visible image 301, three bands that generate near-infrared images, and ten observation bands that generate the infrared image 302. The meteorological satellite 6 acquires the visible image 301 and the infrared image 302 at an observation time interval of 2.5 minutes.

気象衛星画像には、雲頂部分が映される。そして、気象衛星画像では、画像の任意の点に対して、緯度及び経度で表される地上における位置が得られる。したがって、気象衛星画像を用いることで、気象衛星画像内の雲の水平方向の距離を把握することが可能である。 Weather satellite images show cloud tops. For any point on the image, the satellite image provides a location on the ground expressed as latitude and longitude. Therefore, by using the satellite image, it is possible to determine the horizontal distance of clouds in the satellite image.

可視画像301は、太陽光の反射であるアルベドを観測する画像である。可視画像301において雲がある部分は白く映り、雲がない部分は黒く映る。可視画像301では、それぞれの雲の高さの違いは表現されない。また、可視画像301では、雲の厚さの違いは表現されない。 Visible image 301 is an image that observes the albedo, which is the reflection of sunlight. In visible image 301, areas with clouds appear white, and areas without clouds appear black. Visible image 301 does not show differences in the height of each cloud. Also, visible image 301 does not show differences in the thickness of the clouds.

一方、赤外画像302は、雲の放射輝度温度を表す画像である。ここで、雲の最も高い位置は雲頂と呼ばれる。赤外画像302において雲頂が高い雲は白く映り、雲頂が低い部分は黒く映る。すなわち、赤外画像302を用いることで、雲の高さを把握することができる。 On the other hand, the infrared image 302 is an image that represents the radiance temperature of the cloud. Here, the highest point of the cloud is called the cloud top. In the infrared image 302, clouds with high tops appear white, and parts with low tops appear black. In other words, the height of the clouds can be grasped by using the infrared image 302.

また、図5の紙面に向かって右端の輝度温度差画像303は、異なる観測バンドを用いて撮影された赤外画像302の差分の画像である。輝度温度差画像303は、雲が厚い場合は白く映り、雲が薄い場合は黒く映る。すなわち、赤外画像302を用いて輝度温度差画像303を求めることで、雲のおおよその厚さを把握することができる。 The brightness temperature difference image 303 on the far right side of the page in Figure 5 is an image of the difference between infrared images 302 taken using different observation bands. The brightness temperature difference image 303 appears white when the clouds are thick, and appears black when the clouds are thin. In other words, by obtaining the brightness temperature difference image 303 using the infrared image 302, the approximate thickness of the clouds can be grasped.

画像解析部101は、気象衛星画像である可視画像301及び赤外画像302を用いて画像解析を行い、各領域の雲のありなしの判定、個々の雲の識別を行う。また、画像解析部101は、時間毎の衛星画像を蓄積して個別の雲の時系列変化の把握も行う。 The image analysis unit 101 performs image analysis using meteorological satellite images, such as visible images 301 and infrared images 302, to determine whether there are clouds in each area and to identify individual clouds. The image analysis unit 101 also accumulates satellite images taken at each hour to grasp the time series changes of individual clouds.

また、画像解析部101は、天球画像201を用いて画像解析を行い、各領域の雲のありなしの判定、個々の雲の識別及び雲の種類の特定を行う。また、画像解析部101は、時間毎の天球画像201を蓄積して個別の雲の時系列変化の把握も行う。さらに、画像解析部101は、識別した天球画像201における個々の雲の端部の数点を用いて、個々の雲の地上センサ2からの方位角及び天頂角を算出する。ここでの天頂角とは、天頂からの角度を指す。 The image analysis unit 101 also performs image analysis using the celestial image 201 to determine whether or not there are clouds in each area, identify each cloud, and specify the type of cloud. The image analysis unit 101 also accumulates the celestial image 201 for each hour to grasp the time series changes of each cloud. Furthermore, the image analysis unit 101 calculates the azimuth angle and zenith angle of each cloud from the ground sensor 2 using several points on the edge of each cloud in the identified celestial image 201. The zenith angle here refers to the angle from the zenith.

さらに、画像解析部101は、天球画像又は気象衛星画像のいずれかもしくは双方を用いて画像解析を行い、天球画像及び気象衛星画像に映った雲の種類を特定する。画像解析部101による雲の種類の判定方法は様々な方法が存在する。以下に雲の種類判定方法のいくつかの例を説明する。 Furthermore, the image analysis unit 101 performs image analysis using either or both of the celestial image and the meteorological satellite image to identify the type of clouds shown in the celestial image and the meteorological satellite image. There are various methods for the image analysis unit 101 to determine the type of cloud. Some examples of the cloud type determination methods are described below.

例えば、画像解析部101は、図6に示す雲種の判定テーブル304を有する。図6は、雲種の判定テーブルの一例を示す図である。判定テーブル304に示すように、雲の種類は、輝度温度と予め決められたバンドを用いた場合の輝度温度差との組み合わせから識別可能である。図6の輝度温度差におけるB13は10.4μmの波長帯(バンド)の赤外画像を表す、B15は12.3μmの波長帯の赤外画像を表す。すなわち、輝度温度差B13-B15は、10.4μmの波長帯(バンド)の赤外画像と12.3μmの波長帯の赤外画像との差から得られる輝度温度差である。そこで、画像解析部101は、赤外画像302から輝度温度及び輝度温度差を赤外画像302から求める。そして、画像解析部101は、求めた輝度温度及び輝度温度差から雲の種類を特定する。 For example, the image analysis unit 101 has a cloud type determination table 304 shown in FIG. 6. FIG. 6 is a diagram showing an example of the cloud type determination table. As shown in the determination table 304, the type of cloud can be identified from a combination of the brightness temperature and the brightness temperature difference when a predetermined band is used. B13 in the brightness temperature difference in FIG. 6 represents an infrared image in the 10.4 μm wavelength band, and B15 represents an infrared image in the 12.3 μm wavelength band. That is, the brightness temperature difference B13-B15 is the brightness temperature difference obtained from the difference between the infrared image in the 10.4 μm wavelength band and the infrared image in the 12.3 μm wavelength band. Then, the image analysis unit 101 obtains the brightness temperature and the brightness temperature difference from the infrared image 302. Then, the image analysis unit 101 identifies the type of cloud from the obtained brightness temperature and brightness temperature difference.

他にも、画像解析部101は、既存の天球画像や衛星画像およびそれぞれの画像における雲の種類及び状態の情報を用いて、深層学習を行い天球画像や衛星画像と雲の種類及び状態との関係を予め学習する。そして、画像解析部101は、画像解析実行時に取得した天球画像や衛星画像を入力として学習結果を用いて、その天球画像や衛星画像に映された雲の種類を特定する。この方法では、画像解析部101は、天球画像又は衛星画像のいずれか一方を用いてもよいし双方を用いてもよい。 Additionally, the image analysis unit 101 performs deep learning using existing celestial images and satellite images and information on the types and conditions of clouds in each image to learn in advance the relationship between the celestial images and satellite images and the types and conditions of clouds. The image analysis unit 101 then uses the learning results with the celestial images and satellite images acquired when performing image analysis as input to identify the types of clouds shown in the celestial images and satellite images. In this method, the image analysis unit 101 may use either the celestial images or the satellite images, or both.

以下では、画像解析部101が可視画像301及び赤外画像302の画像解析により求めた各雲に紐づけされた雲の種類の判定結果を含む画像情報を、「衛星画像解析情報」と呼ぶ。また、画像解析部101が天球画像201の画像解析により求めた情報を、「天球画像解析情報」と呼ぶ。本実施例では、衛星画像解析情報と天球画像解析情報ともに雲の種類の判定結果が含まれるとして説明する。 In the following, the image information including the cloud type determination results associated with each cloud obtained by the image analysis unit 101 through image analysis of the visible image 301 and the infrared image 302 is referred to as "satellite image analysis information". Also, the information obtained by the image analysis unit 101 through image analysis of the celestial image 201 is referred to as "celestial image analysis information". In this embodiment, both the satellite image analysis information and the celestial image analysis information are described as including the cloud type determination results.

画像解析部101は、天球画像201を変換部102へ出力する。また、画像解析部101は、可視画像301及び赤外画像302を含む気象衛星画像、衛星画像解析情報、天球画像解析情報、並びに、地上気象データを画像照合部104へ出力する。このとき、画像解析部101は、地上センサ2の緯度経度及び平面画像202の元となった天球画像203の撮影時刻も画像照合部104へ出力する。さらに、画像解析部101は、地上気象データ及び赤外画像302を雲頂高度算出部103へ出力する。 The image analysis unit 101 outputs the celestial image 201 to the conversion unit 102. The image analysis unit 101 also outputs meteorological satellite images including a visible image 301 and an infrared image 302, satellite image analysis information, celestial image analysis information, and ground meteorological data to the image matching unit 104. At this time, the image analysis unit 101 also outputs the latitude and longitude of the ground sensor 2 and the shooting time of the celestial image 203 that was the basis of the planar image 202 to the image matching unit 104. Furthermore, the image analysis unit 101 outputs the ground meteorological data and the infrared image 302 to the cloud top height calculation unit 103.

変換部102は、天球画像201の入力を画像解析部101から受ける。魚眼レンズ24は、図4に示すように画角が180度に近づく、すなわち天頂から離れるにしたがい入射ルートの屈折率が大きくなる。このため、魚眼レンズ24を用いて撮影された天球画像201は、天頂から離れるにしたがい歪みが大きくなる。これに対して、気象衛星6により撮影される気象衛星画像は平面として撮影された画像である。そのため、天球画像と気象衛星画像とを比較照合するには、天球画像201を平面として撮影された画像に変換することが好ましい。そこで、変換部102は、図4に示すように、天球画像201を平面展開して平面画像202に変換する。 The conversion unit 102 receives an input of a celestial image 201 from the image analysis unit 101. As shown in FIG. 4, the angle of view of the fisheye lens 24 approaches 180 degrees, i.e., the further away from the zenith, the greater the refractive index of the incident route. For this reason, the celestial image 201 captured using the fisheye lens 24 becomes more distorted the further away from the zenith. In contrast, the meteorological satellite image captured by the meteorological satellite 6 is an image captured as a plane. For this reason, in order to compare and match the celestial image with the meteorological satellite image, it is preferable to convert the celestial image 201 into an image captured as a plane. Therefore, the conversion unit 102 converts the celestial image 201 into a plane image 202 by expanding it into a plane, as shown in FIG. 4.

ここで、変換部102による平面展開の方法の一例について説明する。変換部102は、天球カメラ21により調整用の平面を撮影することで得られた複数の調整用画像を予め蓄積する。そして、変換部102は、複数の調整用画像を画集データとしてキャリブレーションを行い、求めた曲値を基に魚眼レンズ24で撮影した画像を平面展開して、天球画像201を平面展開するためのパラメータの最適値を求める。そして、変換部102は、天球画像201を平面展開する場合、天球画像201の端部付近の歪みが大きい箇所を切り捨てる。これにより、変換部102が生成する平面画像202に含まれる範囲は、天球画像201よりも狭まる。例えば、変換部102は、天頂を0度として天頂を中心とする天頂角が45度の範囲の平面画像202を生成する。例えば、天頂角が45度の場合、上空で雲が存在する高度が1kmから10kmとすれば、天頂角内に納まる円形画像は、地上の水平距離で半径0.5kmから5.0kmの雲の画像となる。 Here, an example of a method of planar development by the conversion unit 102 will be described. The conversion unit 102 accumulates in advance a plurality of images for adjustment obtained by photographing a plane for adjustment with the celestial camera 21. The conversion unit 102 then performs calibration using the plurality of images for adjustment as art book data, and performs planar development of the image photographed by the fisheye lens 24 based on the obtained curvature value, thereby obtaining optimal values of parameters for planar development of the celestial image 201. When the conversion unit 102 performs planar development of the celestial image 201, it cuts off the portions near the ends of the celestial image 201 that have large distortion. As a result, the range included in the planar image 202 generated by the conversion unit 102 is narrower than that of the celestial image 201. For example, the conversion unit 102 generates a planar image 202 with a zenith angle of 45 degrees centered on the zenith, with the zenith at 0 degrees. For example, if the zenith angle is 45 degrees and the altitude at which clouds exist in the sky is between 1 km and 10 km, then a circular image that falls within the zenith angle will be an image of clouds with a radius of 0.5 km to 5.0 km in horizontal distance on the ground.

ただし、ここでは一例として天頂角が45度である平面画像202を用いて説明するが、パラメータをさらに増やして高次元での計算による学習を行うことで、天頂角を広げることができる。例えば、変換部102は、平面画像202の天頂角を60度にすることも可能である。 However, here, as an example, the planar image 202 with a zenith angle of 45 degrees is used for explanation, but the zenith angle can be widened by further increasing the parameters and performing learning through high-dimensional calculations. For example, the conversion unit 102 can set the zenith angle of the planar image 202 to 60 degrees.

そして、変換部102は、学習で得られたパラメータを用いて天球カメラ21により撮影された天球画像201を平面展開して平面画像202を生成する。図4に示すように、変換部102は、天球画像201における端部付近の画像を除外し、45度の天頂角の平面画像202を生成する。ここで、天球画像201を変換して生成された平面画像202も、天球画像201と同様に東西が反転された画像となる。変換部102は、生成した平面画像202を画像照合部104へ出力する。 Then, the conversion unit 102 uses the parameters obtained by learning to expand the celestial sphere image 201 captured by the celestial sphere camera 21 onto a plane to generate a planar image 202. As shown in FIG. 4, the conversion unit 102 removes the image near the edge of the celestial sphere image 201 to generate a planar image 202 with a zenith angle of 45 degrees. Here, the planar image 202 generated by converting the celestial sphere image 201 is also an image with east-west inversion, just like the celestial sphere image 201. The conversion unit 102 outputs the generated planar image 202 to the image matching unit 104.

画像照合部104は、可視画像301及び赤外画像302を含む気象衛星画像、衛星画像解析情報、天球画像解析情報、並びに、地上気象データの入力を画像解析部101から受ける。また、画像照合部104は、地上センサ2の緯度経度及び平面画像202の元となった天球画像203の撮影時刻の入力を画像解析部101から受ける。また、画像照合部104は、平面画像202の入力を変換部102から受ける。 The image matching unit 104 receives inputs of meteorological satellite images including visible images 301 and infrared images 302, satellite image analysis information, celestial image analysis information, and ground meteorological data from the image analysis unit 101. The image matching unit 104 also receives inputs of the latitude and longitude of the ground sensor 2 and the capture time of the celestial image 203 that is the basis of the planar image 202 from the image analysis unit 101. The image matching unit 104 also receives inputs of the planar image 202 from the conversion unit 102.

ここで、平面画像202は、天頂角に応じた範囲の画像であることから、雲の高さによって映る領域が異なる。すなわち、高度が高い雲であれば広い範囲の雲が平面画像202に納まり、高度が低い雲であれば狭い範囲の雲が平面画像202に納まることになる。そのため、雲の高さをある程度推定することで、平面画像202に映った雲の範囲が把握でき、可視画像301に映った雲と平面画像202に映った雲とを照合することが可能となる。 Here, since the planar image 202 is an image of a range according to the zenith angle, the area captured varies depending on the height of the clouds. That is, if the clouds are at a high altitude, a wide area of the clouds will be captured in the planar image 202, and if the clouds are at a low altitude, a narrow area of the clouds will be captured in the planar image 202. Therefore, by estimating the cloud height to a certain extent, it is possible to grasp the range of the clouds captured in the planar image 202, and it becomes possible to compare the clouds captured in the visible image 301 with the clouds captured in the planar image 202.

そこで、画像照合部104は、衛星画像解析情報及び天球画像解析情報に含まれる可視画像301及び天球画像201に映った雲の種類から雲の高度を推定する。次に、平面画像202は東西が反転された画像であるので、画像照合部104は、平面画像202の東西を反転して可視画像301と平面画像202との方位を合わせる。次に、画像照合部104は、地上センサ2の緯度及び経度の情報から、図7に示すように、平面画像202が可視画像301のどの領域に位置するかを特定する。図7は、画像照合解析処理の一例を示す図である。次に、画像照合部104は、推定した雲の高さを用いて、平面画像202に映された雲が含まれる可視画像301の領域における同程度の範囲で形状等が類似する雲を平面画像202に映った雲として抽出する。 The image matching unit 104 estimates the cloud altitude from the type of clouds shown in the visible image 301 and the celestial image 201 included in the satellite image analysis information and the celestial image analysis information. Next, since the planar image 202 is an image in which east and west are reversed, the image matching unit 104 reverses the east and west of the planar image 202 to align the orientations of the visible image 301 and the planar image 202. Next, the image matching unit 104 specifies in which area of the visible image 301 the planar image 202 is located, as shown in FIG. 7, from the latitude and longitude information of the ground sensor 2. FIG. 7 is a diagram showing an example of an image matching analysis process. Next, the image matching unit 104 uses the estimated cloud altitude to extract clouds with similar shapes, etc., within the same range in the area of the visible image 301 that includes the clouds shown in the planar image 202, as clouds shown in the planar image 202.

その後、図8に示すように、画像照合部104は、平面画像202に対して拡大、回転及び平行移動を施して、可視画像301に映った雲と平面画像202に映った雲とをマッチングする。図8は、雲のマッチング処理を説明するための図である。図8における点線で表される各枠が拡大及び平行移動が加えられた平面画像202を表す。図8では、説明の都合上、回転は除いて図示したが、実際のマッチングでは平面画像202に対して回転が加えられてもよい。画像照合部104は、平面画像202における雲の距離が推定されていることから、このマッチングを実行することが可能となる。 Then, as shown in FIG. 8, the image matching unit 104 enlarges, rotates, and translates the planar image 202 to match the clouds captured in the visible image 301 with the clouds captured in the planar image 202. FIG. 8 is a diagram for explaining the cloud matching process. Each frame represented by a dotted line in FIG. 8 represents the planar image 202 to which enlargement and translation have been applied. For convenience of explanation, FIG. 8 shows the planar image 202 without rotation, but in actual matching, the planar image 202 may be rotated. The image matching unit 104 is able to perform this matching because the distance of the clouds in the planar image 202 has been estimated.

このように、画像照合部104は、平面画像202に映った雲と可視画像301に映った雲とを一致させて映った雲を同定する画像照合解析処理を行う。その後、画像照合部104は、画像照合解析処理を行った平面画像202及び可視画像301を範囲取得部105へ出力する。また、画像照合部104は、衛星画像解析情報、天球画像解析情報及び地上気象データを範囲取得部105へ出力する。また、画像照合部104は、画像照合解析処理を施した可視画像301と平面画像202との対応関係から、赤外画像302における平面画像202に対応する領域を特定する。そして、画像照合部104は、赤外画像302における平面画像202に対応する領域の情報を雲頂高度算出部103へ出力する。さらに、画像照合部104は、衛星画像解析情報及び天球画像解析情報を雲厚算出部106へ出力する。 In this way, the image matching unit 104 performs an image matching analysis process to identify the clouds captured in the planar image 202 by matching the clouds captured in the visible image 301. After that, the image matching unit 104 outputs the planar image 202 and the visible image 301 that have been subjected to the image matching analysis process to the range acquisition unit 105. The image matching unit 104 also outputs satellite image analysis information, celestial image analysis information, and ground meteorological data to the range acquisition unit 105. The image matching unit 104 also identifies the area in the infrared image 302 that corresponds to the planar image 202 from the correspondence between the visible image 301 and the planar image 202 that have been subjected to the image matching analysis process. The image matching unit 104 then outputs information on the area in the infrared image 302 that corresponds to the planar image 202 to the cloud top height calculation unit 103. Furthermore, the image matching unit 104 outputs the satellite image analysis information and the celestial image analysis information to the cloud thickness calculation unit 106.

範囲取得部105は、画像照合解析処理が施された平面画像202及び可視画像301の入力を画像照合部104から受ける。また、範囲取得部105は、衛星画像解析情報、天球画像解析情報及び地上気象データの入力を画像照合部104から受ける。範囲取得部105は、気象衛星画像である可視画像301を用いることで水平方向の正確な距離を算出することができる。そこで、範囲取得部105は、図8に示すように、可視画像301を用いて画像照合解析処理が施された平面画像202に映った雲の水平方向の距離を求める。衛星画像の雲と天球の平面画像の雲との同定を行うことで、平面画像に距離情報が加わる。同じ種類の雲であれば雲底の高度は同じと仮定することができるため、範囲取得部105は、雲の規模(大きさ)を測定するのに必要な任意の2地点間の距離の情報を求めることができる。図8における距離LHが、平面画像202に映った雲の水平方向の距離にあたる。 The range acquisition unit 105 receives inputs of the planar image 202 and the visible image 301 that have been subjected to image matching analysis processing from the image matching unit 104. The range acquisition unit 105 also receives inputs of satellite image analysis information, celestial image analysis information, and ground meteorological data from the image matching unit 104. The range acquisition unit 105 can calculate accurate horizontal distances by using the visible image 301, which is a meteorological satellite image. Therefore, as shown in FIG. 8, the range acquisition unit 105 uses the visible image 301 to obtain the horizontal distance of clouds reflected in the planar image 202 that has been subjected to image matching analysis processing. By identifying clouds in the satellite image and clouds in the planar image of the celestial sphere, distance information is added to the planar image. Since it can be assumed that the altitude of the cloud base is the same for the same type of cloud, the range acquisition unit 105 can obtain information on the distance between any two points required to measure the scale (size) of the cloud. The distance LH in FIG. 8 corresponds to the horizontal distance of the cloud reflected in the planar image 202.

その後、範囲取得部105は、平面画像202に映った雲の水平方向の距離の情報とともに、衛星画像解析情報、天球画像解析情報及び地上気象データを雲底高度算出部107へ出力する。 Then, the range acquisition unit 105 outputs the satellite image analysis information, the celestial image analysis information, and the ground meteorological data, together with information on the horizontal distance of the clouds shown in the planar image 202, to the cloud base height calculation unit 107.

雲底高度算出部107は、天球画像解析情報から、平面画像202に映った雲の端部の地上センサ2の天頂からの方位角及び天頂角を取得する。そして、雲底高度算出部107は、平面画像202に映った雲の水平方向の距離、並びに、平面画像202に映った雲の端部の地上センサ2の天頂からの方位角及び天頂角を用いて雲の最下部の高度である雲底高度を算出する。 The cloud base altitude calculation unit 107 obtains the azimuth angle and zenith angle of the edge of the cloud captured in the planar image 202 from the zenith of the ground sensor 2 from the zenith, from the celestial image analysis information. The cloud base altitude calculation unit 107 then calculates the cloud base altitude, which is the altitude of the bottom of the cloud, using the horizontal distance of the cloud captured in the planar image 202, and the azimuth angle and zenith angle of the edge of the cloud captured in the planar image 202 from the zenith of the ground sensor 2.

雲底高度算出部107は、方位角により雲の天頂角を求める。図9は、天頂角と方位角とを説明するための図である。αは方位角であり、方向と真南とのなす角であり、西周りに計測した角度である。すなわち真南が方位角0度の方位である。そして、雲底高度算出部107は、図9に示すように、雲の位置を表す3次元座標が(Rsinθcosα,Rsinθsinα,Rcosθ)と表されることから天頂角θを算出する。 The cloud base altitude calculation unit 107 determines the zenith angle of the cloud from the azimuth angle. Figure 9 is a diagram for explaining the zenith angle and azimuth angle. α is the azimuth angle, which is the angle between a direction and due south, and is an angle measured around the west. In other words, due south is the direction with an azimuth angle of 0 degrees. Then, as shown in Figure 9, the cloud base altitude calculation unit 107 calculates the zenith angle θ from the three-dimensional coordinates representing the position of the cloud, which are expressed as (R sin θ cos α, R sin θ sin α, R cos θ).

図10は、雲底高度の算出の一例を示す図である。図10は、地上センサ2の真上の雲の雲底高度を算出する場合を示す。ここでは、天頂角がθ度である場合で説明する。この場合、雲底高度算出部107は、雲底高度を次の数式(1)で算出する。Hは、雲底高度を表す。 Figure 10 is a diagram showing an example of calculating the cloud base height. Figure 10 shows a case where the cloud base height of the cloud directly above the ground sensor 2 is calculated. Here, a case where the zenith angle is θ degrees is explained. In this case, the cloud base height calculation unit 107 calculates the cloud base height using the following formula (1). H represents the cloud base height.

Figure 0007520649000001
Figure 0007520649000001

例えば、平面画像202における水平方向の雲の範囲が18kmであり、天頂角が45度である場合、雲底高度算出部107は、雲底高度を約9kmと算出する。 For example, if the horizontal cloud range in the planar image 202 is 18 km and the zenith angle is 45 degrees, the cloud base altitude calculation unit 107 calculates the cloud base altitude to be approximately 9 km.

さらに、雲底高度算出部107は、地上センサ2により取得された地上気象データに含まれる地上における気温、気圧及び相対湿度を用いて、持ち上げ凝結高度を求める。そして、雲底高度算出部107は、求めた持ち上げ凝結高度から雲底高度の推定値を算出する。ここで、雲底高度の推定値とは、求めた持ち上げ凝結高度から求まるおおよその雲底高度である。そして、雲底高度算出部107は、雲底高度の推定値を用いて、可視画像301から求めた雲底高度の検証及び補正を行う。その後、雲底高度算出部107は、検証や補正を行った雲底高度を雲厚算出部106へ出力する。 Furthermore, the cloud base height calculation unit 107 calculates the lifting condensation height using the air temperature, air pressure, and relative humidity on the ground contained in the ground weather data acquired by the ground sensor 2. The cloud base height calculation unit 107 then calculates an estimate of the cloud base height from the calculated lifting condensation height. Here, the estimate of the cloud base height is an approximate cloud base height calculated from the calculated lifting condensation height. The cloud base height calculation unit 107 then uses the estimate of the cloud base height to verify and correct the cloud base height calculated from the visible image 301. The cloud base height calculation unit 107 then outputs the verified and corrected cloud base height to the cloud thickness calculation unit 106.

図2に戻って説明を続ける。雲頂高度算出部103は、赤外画像302の入力を画像解析部101から受ける。また、雲頂高度算出部103は、赤外画像302における平面画像202に対応する領域の情報の入力を画像照合部104から受ける。次に、雲頂高度算出部103は、平面画像202に対応する赤外画像302の領域に含まれる雲の輝度から雲頂の温度を算出する。また、雲頂高度算出部103は、気象データ取得部22から地上の温度を取得する。 Returning to FIG. 2, the explanation will be continued. The cloud top height calculation unit 103 receives an input of the infrared image 302 from the image analysis unit 101. The cloud top height calculation unit 103 also receives an input of information on the area in the infrared image 302 that corresponds to the planar image 202 from the image matching unit 104. Next, the cloud top height calculation unit 103 calculates the temperature of the cloud top from the brightness of the clouds included in the area of the infrared image 302 that corresponds to the planar image 202. The cloud top height calculation unit 103 also acquires the ground temperature from the meteorological data acquisition unit 22.

次に、雲頂高度算出部103は、地上の温度と雲頂の温度との温度差を算出する。そして、雲頂高度算出部103は、算出した温度差に対して気温減率を用いて雲頂高度を算出する。例えば、標準大気の気温減率は6.5K/kmである。その後、雲頂高度算出部103は、算出した雲頂高度を雲厚算出部106へ出力する。 Next, the cloud top height calculation unit 103 calculates the temperature difference between the temperature on the ground and the temperature at the cloud top. The cloud top height calculation unit 103 then calculates the cloud top height using the temperature lapse rate for the calculated temperature difference. For example, the temperature lapse rate of standard atmosphere is 6.5 K/km. The cloud top height calculation unit 103 then outputs the calculated cloud top height to the cloud thickness calculation unit 106.

例えば、地上センサ2で計測された地上の温度が20℃であり、赤外画像302から求められた雲頂の温度が220Kの場合で説明する。この場合、雲頂高度算出部103は、220-(273.15+20)=-73.15K/-6.5K=11.25から、雲頂高度を約11kmと算出する。 For example, let us assume that the ground temperature measured by the ground sensor 2 is 20°C, and the cloud top temperature obtained from the infrared image 302 is 220K. In this case, the cloud top height calculation unit 103 calculates the cloud top height to be approximately 11km from 220-(273.15+20)=-73.15K/-6.5K=11.25.

雲厚算出部106は、雲底高度の入力を雲底高度算出部107から受ける。また、雲厚算出部106は、雲頂高度の入力を雲頂高度算出部103から受ける。図11は、雲の厚さと雲頂高度及び雲底高度との関係を示す図である。図11において、距離Hは雲底高度を表し、距離H’は雲頂高度を表し、距離LTは雲の厚さを表し、距離LWは雲の規模にあたる雲の水平距離を表す。図11に示すように、距離LT=距離H’-距離Hである。そこで、雲厚算出部106は、雲頂高度から雲底高度を減算して、雲の厚さを算出する。例えば、雲頂高度が11kmであり、雲底高度が9kmの場合、雲厚算出部106は、雲の厚さを2kmと算出する。その後、雲厚算出部106は、求めた雲の厚さを様相識別部108へ出力する。 The cloud thickness calculation unit 106 receives an input of the cloud base height from the cloud base height calculation unit 107. The cloud thickness calculation unit 106 also receives an input of the cloud top height from the cloud top height calculation unit 103. FIG. 11 is a diagram showing the relationship between the cloud thickness and the cloud top height and the cloud base height. In FIG. 11, distance H represents the cloud base height, distance H' represents the cloud top height, distance LT represents the cloud thickness, and distance LW represents the horizontal distance of the cloud, which corresponds to the scale of the cloud. As shown in FIG. 11, distance LT = distance H' - distance H. Therefore, the cloud thickness calculation unit 106 subtracts the cloud base height from the cloud top height to calculate the cloud thickness. For example, if the cloud top height is 11 km and the cloud base height is 9 km, the cloud thickness calculation unit 106 calculates the cloud thickness to be 2 km. The cloud thickness calculation unit 106 then outputs the calculated cloud thickness to the aspect identification unit 108.

ここで、雲厚算出部106は、衛星画像解析情報及び天球画像解析情報から取得できる雲の厚さの情報を用いて、求めた雲の厚さの検証及び補正を行ってもよい。例えば、雲厚算出部106は、衛星画像解析情報及び天球画像解析情報を画像照合部104から取得する。次に、雲厚算出部106は、衛星画像解析情報及び天球画像解析情報に含まれる雲の種類の情報から雲の厚さの推定値を求める。そして、雲厚算出部106は、雲の厚さの推定値を用いて、雲底高度及び雲頂高度から求めた雲の厚さの検証や補正を行う。 Here, the cloud thickness calculation unit 106 may verify and correct the calculated cloud thickness using cloud thickness information that can be obtained from the satellite image analysis information and celestial image analysis information. For example, the cloud thickness calculation unit 106 obtains the satellite image analysis information and celestial image analysis information from the image matching unit 104. Next, the cloud thickness calculation unit 106 obtains an estimate of the cloud thickness from the cloud type information contained in the satellite image analysis information and celestial image analysis information. Then, the cloud thickness calculation unit 106 uses the estimate of the cloud thickness to verify and correct the cloud thickness calculated from the cloud base altitude and cloud top altitude.

様相識別部108は、雲の厚さの情報の入力を雲厚算出部106から受ける。そして、様相識別部108は、雲の厚さの情報を基に、衛星画像等を用いて、雲の規模LWといった雲の様相を識別する。その後、様相識別部108は、地上センサ2の天球カメラ21により撮影された雲の様相を気象情報データベース3へ送信して格納させる。 The appearance identification unit 108 receives input of cloud thickness information from the cloud thickness calculation unit 106. Then, based on the cloud thickness information, the appearance identification unit 108 identifies the cloud appearance, such as the cloud scale LW, using satellite images, etc. Then, the appearance identification unit 108 transmits the cloud appearance photographed by the celestial camera 21 of the ground sensor 2 to the meteorological information database 3 for storage.

図1に戻って説明を続ける。気象情報データベース3は、各種気象情報を保持管理するデータベースである。気象情報データベース3は、気象情報処理装置1の様相識別部108により識別された地上センサ2の天球カメラ21により撮影された雲の様相を取得して格納する。また、気象情報データベース3は、雲底高度、雲頂高度、雲の厚さを格納する。他にも、気象情報データベース3は、気象庁高解像度降水ナウキャスト等から得られるレーダー雨量、気象庁LIDEN(Lightning Detection Network system)から得られる雷情報、及び、気象庁の竜巻等のトップデータベースから得られる突風情報等の各種気象情報を格納する。 Returning to FIG. 1, we will continue the explanation. The weather information database 3 is a database that holds and manages various weather information. The weather information database 3 acquires and stores the cloud appearance photographed by the celestial camera 21 of the ground sensor 2 that is identified by the appearance identification unit 108 of the weather information processing device 1. The weather information database 3 also stores cloud base altitude, cloud top altitude, and cloud thickness. In addition, the weather information database 3 stores various weather information such as radar rainfall obtained from the Japan Meteorological Agency's high-resolution precipitation nowcast, lightning information obtained from the Japan Meteorological Agency's LIDEN (Lightning Detection Network system), and wind gust information obtained from the Japan Meteorological Agency's top database of tornadoes, etc.

気象情報データベース3に、雲の種類、雲底高度、雲頂高度、雲の厚さ及び雲の規模を蓄積していき、それらの情報を用いて機械学習を行うことで、いずれかの情報が得られない場合でも、他の情報から得られない情報を補うことが可能となる。この情報の補填は、情報活用基盤4が実行してもよいし、他のコンピュータにより実行されてもよい。 By storing cloud type, cloud base altitude, cloud top altitude, cloud thickness, and cloud scale in the meteorological information database 3 and using this information to perform machine learning, even if one piece of information cannot be obtained, it is possible to supplement the missing information from other information. This information supplementation may be performed by the information utilization platform 4 or by another computer.

情報活用基盤4は、気象情報データベース3に格納された雲の様相を含む各種気象情報を取得して、各種気象情報活用システム7で利用し易い形に加工する。例えば、情報活用基盤4は、雲の種類、規模、形状と、レーダー雨量、雷情報及び突風情報等の情報と関連付ける。これにより、情報活用基盤4は、様々な施設や設備の運用及び管理活用することができるデータに各種気象情報を加工することができる。情報活用基盤4は、気象情報を活用する各種気象情報活用システム7へ加工したデータを提供する。 The information utilization platform 4 acquires various weather information, including cloud appearances, stored in the weather information database 3, and processes the information into a form that is easy to use in the various weather information utilization systems 7. For example, the information utilization platform 4 associates the type, size, and shape of clouds with information such as radar rainfall, lightning information, and wind gust information. This allows the information utilization platform 4 to process the various weather information into data that can be used to operate and manage various facilities and equipment. The information utilization platform 4 provides the processed data to the various weather information utilization systems 7 that utilize the weather information.

各種気象情報活用システム7は、例えば、太陽光発電システム71、発電所運用管理システム72、河川管理システム73、航空運航管理システム74、山岳情報提供システム75、気象予報システム76及び気象情報視覚化システム77等を含む。ただし、図1に示したシステムは一例であり、気象情報を活用することが可能なシステムであれば、他のシステムであってもよい。 The various weather information utilization systems 7 include, for example, a solar power generation system 71, a power plant operation management system 72, a river management system 73, an air traffic control system 74, a mountain information provision system 75, a weather forecast system 76, and a weather information visualization system 77. However, the system shown in FIG. 1 is only an example, and other systems may be used as long as they are capable of utilizing weather information.

例えば、太陽光発電システム71であれば、情報活用基盤4から提供された雲の様相を利用することで、より正確な太陽光発電による出力予測等を行うことが可能となる。また、発電所運用管理システム72であれば、情報活用基盤4から提供された雲の様相を利用することで、電力供給量の調整をより適切に行うことが可能となる。また、河川管理システム73であれば、情報活用基盤4から提供された雲の様相を利用することで、より正確な雨の予測ができ、川の氾濫等を正確に予測して適切な河川管理を行うことが可能となる。また、航空運航管理システム74であれば、情報活用基盤4から提供された雲の様相を利用することで、航路上の雲の状態等を正確に把握することができ、航空運航の管理をより適切に行うことができる。また、山岳情報提供システム75や気象予報システム76であれば、情報活用基盤4から提供された雲の様相を利用することで、より正確に将来の気象予測を行うことができ、適切な情報提供を行うことができる。また、気象情報視覚化システム77であれば、情報活用基盤4から提供された雲の様相を利用することで、より正確な雲の画像を提供することが可能となる。 For example, the photovoltaic power generation system 71 can use the cloud aspects provided by the information utilization platform 4 to more accurately predict the output of photovoltaic power generation. The power plant operation management system 72 can use the cloud aspects provided by the information utilization platform 4 to more appropriately adjust the amount of power supply. The river management system 73 can use the cloud aspects provided by the information utilization platform 4 to more accurately predict rain, and can accurately predict river flooding and perform appropriate river management. The air traffic management system 74 can use the cloud aspects provided by the information utilization platform 4 to accurately grasp the cloud conditions on the route, and can more appropriately manage air traffic. The mountain information provision system 75 and the weather forecast system 76 can use the cloud aspects provided by the information utilization platform 4 to more accurately predict future weather and provide appropriate information. The weather information visualization system 77 can use the cloud aspects provided by the information utilization platform 4 to more accurately provide cloud images.

次に、図12を参照して、雲の様相を含む気象情報提供の流れについて説明する。図12は、雲の様相を含む気象情報提供のフローチャートである。 Next, the flow of providing weather information including cloud appearance will be described with reference to FIG. 12. FIG. 12 is a flowchart of providing weather information including cloud appearance.

地上センサ2は、天球カメラ21を用いて天球画像201を撮影する。また、地上センサ2は、気象データ取得部22を用いて地上の気温、気圧及び相対湿度を含む地上気象データを取得する(ステップS1)。そして、地上センサ2は、天球画像201及び地上気象データを気象情報処理装置1へ送信する。 The ground sensor 2 captures a celestial image 201 using the celestial camera 21. The ground sensor 2 also acquires ground meteorological data including ground temperature, air pressure, and relative humidity using the meteorological data acquisition unit 22 (step S1). The ground sensor 2 then transmits the celestial image 201 and the ground meteorological data to the weather information processing device 1.

気象情報処理装置1は、天球画像201及び地上気象データを地上センサ2から受信する。また、気象情報処理装置1は、可視画像301及び赤外画像302を含む気象衛星画像を気象データ配信システム5から取得する。そして、気象情報処理装置1は、天球画像201、地上気象データ、並びに、可視画像301及び赤外画像302を含む気象衛星画像を用いて天球画像201に映った雲の厚さを求めてその雲の様相を識別する(ステップS2)。 The weather information processing device 1 receives a celestial image 201 and ground meteorological data from the ground sensor 2. The weather information processing device 1 also acquires meteorological satellite images including a visible image 301 and an infrared image 302 from the meteorological data distribution system 5. The weather information processing device 1 then uses the celestial image 201, the ground meteorological data, and the meteorological satellite images including the visible image 301 and the infrared image 302 to determine the thickness of clouds shown in the celestial image 201 and identify the state of the clouds (step S2).

気象情報処理装置1は、識別した雲の様相の情報を気象情報データベース3に格納する(ステップS3)。 The weather information processing device 1 stores the identified cloud appearance information in the weather information database 3 (step S3).

情報活用基盤4は、気象情報データベース3に格納された雲の様相の情報を含む各種気象情報を、各種気象情報活用システム7が利用し易いように加工する(ステップS4)。 The information utilization platform 4 processes the various weather information, including information on cloud appearance stored in the weather information database 3, so that it can be easily used by the various weather information utilization systems 7 (step S4).

その後、情報活用基盤4は、加工した気象情報を各種気象情報活用システム7へ提供する(ステップS5)。 Then, the information utilization platform 4 provides the processed weather information to various weather information utilization systems 7 (step S5).

各種気象情報活用システム7は、加工された気象情報を情報活用基盤4から取得する。そして、各種気象情報活用システム7は、天気予報や太陽光発電による出力予測等に加工された気象情報を活用する(ステップS6)。 The various weather information utilization system 7 acquires the processed weather information from the information utilization base 4. The various weather information utilization system 7 then utilizes the processed weather information for weather forecasts, solar power generation output predictions, etc. (step S6).

次に、図13を参照して、気象情報処理装置1による雲の様相の識別処理の流れについて説明する。図13は、気象情報処理装置による雲の様相の識別処理のフローチャートである。 Next, the flow of the cloud appearance identification process by the weather information processing device 1 will be described with reference to FIG. 13. FIG. 13 is a flowchart of the cloud appearance identification process by the weather information processing device.

画像解析部101は、天球画像201及び気象データを地上センサ2から取得する(ステップS101)。 The image analysis unit 101 acquires a celestial image 201 and meteorological data from the ground sensor 2 (step S101).

また、画像解析部101は、気象衛星6により撮影された可視画像301及び赤外画像302を含む気象衛星画像を気象データ配信システム5から取得する(ステップS102)。 The image analysis unit 101 also acquires meteorological satellite images, including a visible image 301 and an infrared image 302, taken by a meteorological satellite 6 from the meteorological data distribution system 5 (step S102).

次に、画像解析部101は、天球画像201の画像解析を行い、各領域の雲のありなしの判定、個々の雲の識別、雲の種類の特定、個別の雲の時系列変化、並びに、個々の雲の地上センサ2からの方位角及び天頂角を含む天球画像解析情報を取得する。また、画像解析部101は、気象衛星画像の画像解析を行い、各領域の雲のありなしの判定、個々の雲の識別、雲の種類の特定及び個別の雲の時系列変化を含む衛星画像解析情報を取得する(ステップS103)。その後、画像解析部101は、天球画像201を変換部102へ出力する。また、画像解析部101は、可視画像301及び赤外画像302を含む気象衛星画像、衛星画像解析情報、天球画像解析情報、並びに、地上気象データを画像照合部104へ出力する。このとき、画像解析部101は、地上センサ2の緯度経度及び平面画像202の元となった天球画像203の撮影時刻も画像照合部104へ出力する。さらに、画像解析部101は、地上気象データ及び赤外画像302を雲頂高度算出部103へ出力する。 Next, the image analysis unit 101 performs image analysis of the celestial image 201, and obtains celestial image analysis information including a determination of the presence or absence of clouds in each area, identification of each cloud, identification of the type of cloud, time series changes of each cloud, and the azimuth angle and zenith angle of each cloud from the ground sensor 2. The image analysis unit 101 also performs image analysis of the meteorological satellite image, and obtains satellite image analysis information including a determination of the presence or absence of clouds in each area, identification of each cloud, identification of the type of cloud, and time series changes of each cloud (step S103). The image analysis unit 101 then outputs the celestial image 201 to the conversion unit 102. The image analysis unit 101 also outputs the meteorological satellite image including the visible image 301 and the infrared image 302, the satellite image analysis information, the celestial image analysis information, and the ground meteorological data to the image matching unit 104. At this time, the image analysis unit 101 also outputs the latitude and longitude of the ground sensor 2 and the capture time of the celestial image 203 that was the basis of the planar image 202 to the image matching unit 104. Furthermore, the image analysis unit 101 outputs the ground meteorological data and the infrared image 302 to the cloud top height calculation unit 103.

変換部102は、天球画像201の入力を画像解析部101から受ける。そして、変換部102は、天球画像201を平面展開して平面画像202に変換する(ステップS104)。その後、変換部102は、生成した平面画像202を画像照合部104へ出力する。 The conversion unit 102 receives an input of the celestial sphere image 201 from the image analysis unit 101. The conversion unit 102 then expands the celestial sphere image 201 onto a plane to convert it into a planar image 202 (step S104). The conversion unit 102 then outputs the generated planar image 202 to the image matching unit 104.

画像照合部104は、可視画像301及び赤外画像302を含む気象衛星画像、衛星画像解析情報、天球画像解析情報、地上気象データ、地上センサ2の緯度経度及び平面画像202の元となった天球画像201の撮影時刻の入力を画像解析部101から受ける。そして、画像照合部104は、衛星画像解析情報及び天球画像解析情報から雲の種類を取得する。次に、画像照合部104は、取得した雲の種類から雲の高さを推定する。次に、推定した雲の高さを用いて、推定した雲の高さから平面画像202に納まった雲の水平方向の距離の推定値を求める。さらに、画像照合部104は、平面画像202の東西を反転する。次に、画像照合部104は、地上センサ2の緯度及び経度の情報から、平面画像202が可視画像301のどの領域に位置するかを特定し、可視画像301における平面画像202に映った雲を抽出する。次に、画像照合部104は、平面画像202に対して平行移動、回転及び拡大縮小等を行い、平面画像202と可視画像301とのマッチングを行う。画像照合部104は、このような画像照合解析処理を行うことで、可視画像301に映された雲の中から平面画像202に映された雲を同定する(ステップS105)。その後、画像照合部104は、画像照合解析処理を施した平面画像202及び可視画像301を範囲取得部105へ出力する。 The image matching unit 104 receives inputs of meteorological satellite images including the visible image 301 and the infrared image 302, satellite image analysis information, celestial image analysis information, ground meteorological data, the latitude and longitude of the ground sensor 2, and the shooting time of the celestial image 201 that was the basis of the planar image 202 from the image analysis unit 101. Then, the image matching unit 104 acquires the type of cloud from the satellite image analysis information and the celestial image analysis information. Next, the image matching unit 104 estimates the cloud height from the acquired cloud type. Next, using the estimated cloud height, an estimate of the horizontal distance of the cloud contained in the planar image 202 is obtained from the estimated cloud height. Furthermore, the image matching unit 104 inverts the east-west of the planar image 202. Next, the image matching unit 104 identifies in which area of the visible image 301 the planar image 202 is located from the latitude and longitude information of the ground sensor 2, and extracts the clouds reflected in the planar image 202 from the visible image 301. Next, the image matching unit 104 performs parallel translation, rotation, zooming, etc. on the planar image 202, and matches the planar image 202 with the visible image 301. By performing such image matching analysis processing, the image matching unit 104 identifies clouds captured in the planar image 202 from among clouds captured in the visible image 301 (step S105). After that, the image matching unit 104 outputs the planar image 202 and the visible image 301 that have been subjected to the image matching analysis processing to the range acquisition unit 105.

範囲取得部105は、画像照合解析処理が施された平面画像202及び可視画像301を画像照合部104から取得する。次に、範囲取得部105は、可視画像301を用いて、平面画像202に映された雲の水平方向の距離を求める(ステップS106)。その後、範囲取得部105は、平面画像202に映った雲の水平方向の距離を雲底高度算出部107へ出力する。 The range acquisition unit 105 acquires the planar image 202 and the visible image 301 that have been subjected to the image matching analysis process from the image matching unit 104. Next, the range acquisition unit 105 uses the visible image 301 to determine the horizontal distance of the clouds reflected in the planar image 202 (step S106). After that, the range acquisition unit 105 outputs the horizontal distance of the clouds reflected in the planar image 202 to the cloud base altitude calculation unit 107.

雲底高度算出部107は、平面画像202に映った雲の水平方向の距離の入力を範囲取得部105から受ける。そして、雲底高度算出部107は、平面画像202に映された雲の端部の地上センサ2に対する方位角及び天頂角を用いて、雲の水平方向の距離から平面画像202に映った雲の雲底高度を算出する(ステップS107)。この時、雲底高度算出部107は、地上の温度、気圧及び相対温度から持ち上げ凝結高度を求めて、雲底高度の推定値を算出する。そして、雲底高度算出部107は、雲底高度の推定値を用いて算出した雲底高度の検証及び補正を行う。その後、雲底高度算出部107は、算出した雲底高度を雲厚算出部106へ出力する。 The cloud base height calculation unit 107 receives input of the horizontal distance of the cloud captured in the planar image 202 from the range acquisition unit 105. The cloud base height calculation unit 107 then calculates the cloud base height of the cloud captured in the planar image 202 from the horizontal distance of the cloud using the azimuth angle and zenith angle of the edge of the cloud captured in the planar image 202 relative to the ground sensor 2 (step S107). At this time, the cloud base height calculation unit 107 calculates an estimate of the cloud base height by determining the lift condensation height from the ground temperature, air pressure, and relative temperature. The cloud base height calculation unit 107 then verifies and corrects the calculated cloud base height using the estimated cloud base height. The cloud base height calculation unit 107 then outputs the calculated cloud base height to the cloud thickness calculation unit 106.

雲頂高度算出部103は、赤外画像302及び地上気象データの入力を画像解析部101から受ける。次に、雲頂高度算出部103は、地上気象データ及び赤外画像302から地上の温度及び雲頂温度を取得する(ステップS108)。 The cloud top height calculation unit 103 receives the infrared image 302 and ground meteorological data from the image analysis unit 101. Next, the cloud top height calculation unit 103 obtains the ground temperature and cloud top temperature from the ground meteorological data and the infrared image 302 (step S108).

次に、雲頂高度算出部103は、地上の温度と雲頂温度との温度差に気温減率を用いて雲頂高度を算出する(ステップS109)。その後、雲頂高度算出部103は、算出した雲頂高度を雲厚算出部106へ出力する。 Next, the cloud top height calculation unit 103 calculates the cloud top height using the temperature lapse rate for the temperature difference between the ground temperature and the cloud top temperature (step S109). After that, the cloud top height calculation unit 103 outputs the calculated cloud top height to the cloud thickness calculation unit 106.

雲厚算出部106は、雲底高度の入力を雲底高度算出部107から受ける。また、雲厚算出部106は、雲頂高度の入力を雲頂高度算出部103から受ける。そして、雲厚算出部106は、雲頂高度から雲底高度を減算して、雲の厚さを算出する(ステップS110)。この時、雲厚算出部106は、衛星画像解析情報及び天球画像解析情報に含まれる雲の種類の情報を用いて雲の厚さの推定値を求めて、算出した雲の厚さの検証及び補正を行ってもよい。その後、雲厚算出部106は、算出した雲の厚さを様相識別部108へ出力する。 The cloud thickness calculation unit 106 receives an input of the cloud base altitude from the cloud base altitude calculation unit 107. The cloud thickness calculation unit 106 also receives an input of the cloud top altitude from the cloud top altitude calculation unit 103. The cloud thickness calculation unit 106 then subtracts the cloud base altitude from the cloud top altitude to calculate the cloud thickness (step S110). At this time, the cloud thickness calculation unit 106 may obtain an estimate of the cloud thickness using information on the cloud type contained in the satellite image analysis information and the celestial image analysis information, and verify and correct the calculated cloud thickness. The cloud thickness calculation unit 106 then outputs the calculated cloud thickness to the appearance identification unit 108.

様相識別部108は、雲の厚さの情報を雲厚算出部106から受ける。そして、様相識別部108は、雲の厚さを基に、可視画像301、赤外画像302及び平面画像202を用いて雲の種類、規模及び形状等の平面画像202に映った雲の様相を識別する(ステップS111)その後、様相識別部108は、識別した雲の様相の情報を気象情報データベース3に格納する。 The appearance identification unit 108 receives information on the cloud thickness from the cloud thickness calculation unit 106. Then, based on the cloud thickness, the appearance identification unit 108 identifies the cloud appearance reflected in the planar image 202, such as the type, size, and shape of the cloud, using the visible image 301, the infrared image 302, and the planar image 202 (step S111). After that, the appearance identification unit 108 stores the information on the identified cloud appearance in the weather information database 3.

以上に説明したように、本実施例に係る気象情報処理装置は、地上から撮影された天球画像を平面画像に変換し、平面画像に映る雲を可視画像に映る雲の中から同定する。さらに、気象情報処理装置は、同定した結果を用いて平面画像に映る雲の水平方向の距離を求め、水平方向の距離を基に雲底高度を求める。また、気象情報処理装置は、赤外画像及び地上の温度から雲頂高度を求める。そして、気象情報処理装置は、求めた雲底高度及び雲頂高度から雲の厚さを算出し、算出した雲の厚さを用いて雲の様相を識別する。 As described above, the weather information processing device according to this embodiment converts a celestial image taken from the ground into a planar image, and identifies clouds that appear in the planar image from among the clouds that appear in the visible image. Furthermore, the weather information processing device uses the identification results to determine the horizontal distance of the clouds that appear in the planar image, and determines the cloud base height based on the horizontal distance. The weather information processing device also determines the cloud top height from the infrared image and ground temperature. The weather information processing device then calculates the cloud thickness from the determined cloud base height and cloud top height, and identifies the cloud appearance using the calculated cloud thickness.

これにより、気象情報処理装置は、正確に雲の様相を識別することができ、気象データを活用するシステムの性能を向上させることが可能となる。また、雲の種類、雲底高度、雲頂高度及び雲の厚さを気象情報データベースに蓄積していき、それらの情報を用いて機械学習を行うことで、いずれかの情報が得られない場合でも、他の情報から得られない情報を補うことが可能となる。したがって、情報が十分に得られない場合でも、適切な気象情報を各種気象情報活用システムに提供することが可能となる。 This allows the weather information processing device to accurately identify the state of clouds, improving the performance of systems that utilize weather data. In addition, by accumulating cloud type, cloud base altitude, cloud top altitude, and cloud thickness in a weather information database and using this information to perform machine learning, even if one piece of information is not available, it is possible to compensate for the missing information from the other information. Therefore, even if sufficient information is not available, it is possible to provide appropriate weather information to various weather information utilization systems.

(変形例)
次に、実施例1の変形例について説明する。図14は、高度の違う雲の平面画像における位置関係を説明するための図である。例えば、図14に示すように、低い高度の雲C1と高い高度の雲C2とが存在する場合について説明する。平面画像202では、高度の差は表現されずに、高度の異なる雲C1及びC2はいずれも同一平面に表される。そのため、雲C1及びC2が同じ高度に存在するようにも考えられ、その場合、雲C2は、見かけの位置が雲C1と同じ高度となり、雲C21の位置に存在するように見える。この時、雲C2を雲C1と同じ高度と考えた場合、雲C2の天頂からの距離は、距離L1と算出される。しかし、実際には、雲C2からの天頂からの距離は、距離L2であり、算出された距離L1とは大きく異なる。例えば、雲C1の高度が1000mであり、雲C2の高度が10000mであり、天頂と雲C2との角度φが30度である場合、距離L2は5800mであるのに対して、距離L1は580mとなる。このように、雲の構造が複雑な場合に、平面画像202に映った雲の高度を考慮せずに同一高度として算出すると、雲の水平方向の距離に誤差が発生してしまう。
(Modification)
Next, a modified example of the first embodiment will be described. FIG. 14 is a diagram for explaining the positional relationship in a planar image of clouds at different altitudes. For example, as shown in FIG. 14, a case where a cloud C1 at a low altitude and a cloud C2 at a high altitude exist will be described. In the planar image 202, the difference in altitude is not expressed, and the clouds C1 and C2 at different altitudes are both displayed on the same plane. Therefore, it is also considered that the clouds C1 and C2 exist at the same altitude, and in that case, the apparent position of the cloud C2 is the same altitude as the cloud C1, and it appears to exist at the position of the cloud C21. At this time, if the cloud C2 is considered to be at the same altitude as the cloud C1, the distance from the zenith of the cloud C2 is calculated as the distance L1. However, in reality, the distance from the zenith of the cloud C2 is the distance L2, which is significantly different from the calculated distance L1. For example, if the altitude of cloud C1 is 1000 m, the altitude of cloud C2 is 10000 m, and the angle φ between the zenith and cloud C2 is 30 degrees, then distance L2 is 5800 m, while distance L1 is 580 m. In this way, when the cloud structure is complex, if the altitudes of the clouds shown in the planar image 202 are not taken into consideration and the clouds are calculated as being at the same altitude, an error will occur in the horizontal distance of the clouds.

そこで、本変形例に係る気象情報処理装置1は、雲の構成が複雑な場合に、高さの異なる雲を分類して、それぞれについて水平方向の距離を求める。以下に、変形例に係る気象情報処理装置1による雲の水平方向の距離の求め方について説明する。図15は、実施例1の変形例に係る気象情報処理装置による高さの異なる雲のそれぞれの雲底高度の算出を説明するための図である。 Therefore, when the cloud structure is complex, the weather information processing device 1 according to this modification classifies clouds of different heights and calculates the horizontal distance for each. Below, we will explain how the weather information processing device 1 according to the modification calculates the horizontal distance of clouds. Figure 15 is a diagram for explaining the calculation of the cloud base altitude of each cloud of different heights by the weather information processing device according to the modification of Example 1.

範囲取得部105は、衛星画像解析情報及び天球画像解析情報から、識別された個々の雲の種類の情報を取得する。そして、範囲取得部105は、雲の種類に応じて平面画像202に映った雲を分類することで、雲の高度に応じた分類を行う。ここで、範囲取得部105は、等しい高度に存在する雲については同じ分類としてもよい。そして、範囲取得部105は、分類毎に、天頂からの雲の端部の方位角及び天頂角を用いて雲の水平距離を求める。そして、範囲取得部105は、雲の高度に応じた分類毎の雲の水平距離を雲底高度算出部107へ出力する。 The range acquisition unit 105 acquires information on each identified cloud type from the satellite image analysis information and the celestial image analysis information. The range acquisition unit 105 then classifies the clouds shown in the planar image 202 according to the cloud type, thereby classifying the clouds according to their altitude. Here, the range acquisition unit 105 may assign the same classification to clouds that exist at the same altitude. The range acquisition unit 105 then determines the horizontal distance of the cloud for each classification using the azimuth angle and zenith angle of the cloud edge from the zenith. The range acquisition unit 105 then outputs the horizontal distance of the cloud for each classification according to the cloud altitude to the cloud base altitude calculation unit 107.

例えば、範囲取得部105は、図15に示す雲C1が属する雲のグループと、雲C2が属する雲のグループに雲を分類する。そして、範囲取得部105は、雲C1が属するグループの雲の水平方向の距離を距離LH1と求める。また、範囲取得部105は、雲C2が属するグループの雲の水平方向の距離を距離LH2と求める。このように、範囲取得部105は、異なる高さの雲C1と雲C2とのそれぞれについて水平方向の距離を求める。 For example, the range acquisition unit 105 classifies clouds into a cloud group to which cloud C1 belongs and a cloud group to which cloud C2 belongs, as shown in FIG. 15. The range acquisition unit 105 then determines the horizontal distance of the clouds in the group to which cloud C1 belongs as distance LH1. The range acquisition unit 105 also determines the horizontal distance of the clouds in the group to which cloud C2 belongs as distance LH2. In this way, the range acquisition unit 105 determines the horizontal distance for each of clouds C1 and C2, which are at different heights.

雲底高度算出部107は、雲の高度に応じた分類毎の雲の水平距離を範囲取得部105から取得する。そして、雲底高度算出部107は、雲の高度に応じた分類毎に雲底高度を算出する。例えば、図15に示すように、雲底高度算出部107は、雲C1の雲底高度を距離H1と求める。また、雲底高度算出部107は、雲C2の雲底高度を距離H2と求める。その後、雲底高度算出部107は、雲の高度に応じた分類毎の雲底高度を雲厚算出部106へ出力する。 The cloud base height calculation unit 107 acquires the horizontal distance of the cloud for each classification according to the cloud height from the range acquisition unit 105. Then, the cloud base height calculation unit 107 calculates the cloud base height for each classification according to the cloud height. For example, as shown in FIG. 15, the cloud base height calculation unit 107 calculates the cloud base height of cloud C1 as distance H1. The cloud base height calculation unit 107 also calculates the cloud base height of cloud C2 as distance H2. Then, the cloud base height calculation unit 107 outputs the cloud base height for each classification according to the cloud height to the cloud thickness calculation unit 106.

その後、雲厚算出部106は、雲の高度に応じた分類毎の雲底高度の入力を雲底高度算出部107から受ける。そして、雲厚算出部106は、雲の高度に応じた分類毎の雲底高度を用いて、雲の高度に応じた分類毎に雲の厚さを算出する。 Then, the cloud thickness calculation unit 106 receives an input of the cloud base altitude for each classification according to the cloud altitude from the cloud base altitude calculation unit 107. Then, the cloud thickness calculation unit 106 calculates the cloud thickness for each classification according to the cloud altitude using the cloud base altitude for each classification according to the cloud altitude.

ここで、本実施例では、異なる高度の雲のそれぞれについての様相を識別する場合で説明したが、特定の雲について注目する場合には、気象情報処理装置1は、特定の雲についての情報を取得すればよい。例えば、範囲取得部105は、雲の高度に応じて雲を分類したうえで、注目する特定の雲以外の雲は除外し、特定の雲が含まれる分類に制限して水平方向の算出を行ってもよい。 In this embodiment, the case of identifying the appearance of each cloud at different altitudes has been described, but when focusing on a specific cloud, the weather information processing device 1 can acquire information on the specific cloud. For example, the range acquisition unit 105 may classify clouds according to their altitude, exclude clouds other than the specific cloud of interest, and perform horizontal calculations limited to classifications that include the specific cloud.

以上に説明したように、本変形例に係る気象情報処理装置は、雲の種類を用いて雲を分類することで高度に応じた雲の分類を行う。そして、気象情報処理装置は、雲の高度に応じた分類毎に雲の水平方向の距離を求めて雲の厚さを算出する。これにより、平面画像に映った高度の異なる雲の水平方向の距離の誤差を低減することができる。すなわち、雲の構成が複雑な場合にも、各雲の様相を適切に識別することができ、気象データを活用するシステムの性能を向上させることが可能となる。 As described above, the weather information processing device according to this modified example classifies clouds according to altitude by classifying them using the cloud type. The weather information processing device then calculates the cloud thickness by determining the horizontal distance of the cloud for each classification according to the cloud altitude. This makes it possible to reduce errors in the horizontal distance of clouds at different altitudes captured in a planar image. In other words, even when the cloud structure is complex, it is possible to properly identify the appearance of each cloud, thereby improving the performance of a system that utilizes weather data.

図16は、実施例2に係る気象情報処理装置のブロック図である。本実施例に係る気象情報処理装置1は、雲の様相の識別に加えて、雲及び雲域の移動方向及び移動速度を求める。そこで、以下では、気象情報処理装置1による雲及び雲域の移動方向及び移動速度の算出について主に説明する。本実施例に係る気象情報処理装置1は、実施例1の各部に加えて、移動速度算出部109を有する。以下の説明では、実施例1と同様の各部の動作については説明を省略する。 Figure 16 is a block diagram of a weather information processing device according to Example 2. The weather information processing device 1 according to this example determines the direction and speed of movement of clouds and cloud areas in addition to identifying the state of clouds. Therefore, the following mainly describes the calculation of the direction and speed of movement of clouds and cloud areas by the weather information processing device 1. The weather information processing device 1 according to this example has a speed calculation unit 109 in addition to the components of Example 1. In the following description, the operation of the components similar to those of Example 1 will be omitted.

画像解析部101は、衛星画像解析情報及び天球画像解析情報を移動速度算出部109へ出力する。また、画像照合部104は、画像照合解析処理を施した平面画像202及び可視画像301を移動速度算出部109へ出力する。 The image analysis unit 101 outputs the satellite image analysis information and the celestial image analysis information to the moving speed calculation unit 109. In addition, the image matching unit 104 outputs the planar image 202 and the visible image 301 that have been subjected to the image matching analysis process to the moving speed calculation unit 109.

移動速度算出部109は、衛星画像解析情報及び天球画像解析情報の入力を画像解析部101から受ける。また、移動速度算出部109は、画像照合解析処理を施した平面画像202及び可視画像301の入力を画像照合部104から受ける。 The moving speed calculation unit 109 receives inputs of satellite image analysis information and celestial image analysis information from the image analysis unit 101. The moving speed calculation unit 109 also receives inputs of the planar image 202 and the visible image 301 that have been subjected to image matching analysis processing from the image matching unit 104.

そして、移動速度算出部109は、衛星画像解析情報、天球画像解析情報、並びに、画像照合解析処理を施した平面画像202及び可視画像301を用いて、個別の雲や雲域の時系列変化をトレースする。そして、移動速度算出部109は、トレースした時系列変化を用いて、個別の雲や雲域の移動方向及び移動速度を算出する。 Then, the moving speed calculation unit 109 uses the satellite image analysis information, the celestial image analysis information, and the planar image 202 and the visible image 301 that have been subjected to image matching analysis processing to trace the time series changes of individual clouds and cloud areas.Then, the moving speed calculation unit 109 uses the traced time series changes to calculate the moving direction and moving speed of individual clouds and cloud areas.

図17は、移動方向及び移動速度の算出の一例を説明するための図である。例えば、移動速度算出部109は、12:00に撮影された可視画像401において同定された位置411の雲を抽出する。画像解析部101は、2.5分~10分程度で次の衛星画像を取得する。そこで、移動速度算出部109は、12:10に撮影された可視画像402において同定された位置412の雲を抽出する。そして、移動速度算出部109は、雲の種類や形状等から位置411に存在した雲が位置412まで移動した判定する。そこで、移動速度算出部109は、10分間で位置411が位置412まで雲が移動したことから、矢印P1で表される雲の移動速度及び移動方向を求める。 Figure 17 is a diagram for explaining an example of calculation of the moving direction and moving speed. For example, the moving speed calculation unit 109 extracts the cloud at the position 411 identified in the visible image 401 captured at 12:00. The image analysis unit 101 acquires the next satellite image in about 2.5 to 10 minutes. Then, the moving speed calculation unit 109 extracts the cloud at the position 412 identified in the visible image 402 captured at 12:10. Then, the moving speed calculation unit 109 determines that the cloud that existed at the position 411 has moved to the position 412 based on the type and shape of the cloud. Then, the moving speed calculation unit 109 calculates the moving speed and moving direction of the cloud represented by the arrow P1 based on the fact that the cloud has moved from the position 411 to the position 412 in 10 minutes.

次に、移動速度算出部109は、現在までの個別の雲や雲域の移動方向及び移動速度の履歴から、個別の雲や雲域の将来の移動速度及び移動位置を算出して、個別の雲や雲域の将来の位置を予測する。例えば、図17に示すように、移動速度算出部109は、矢印P1で表される移動方向及び移動速度を含むそれまでの雲の移動方向及び移動速度の履歴から、雲のその後の移動方向及び移動速度を算出する。そして、移動速度算出部109は、その後の雲の移動位置である位置413~415を予測する。 Next, the movement speed calculation unit 109 calculates the future movement speed and position of each individual cloud or cloud area from the history of the movement direction and movement speed of each individual cloud or cloud area up to the present, and predicts the future position of each individual cloud or cloud area. For example, as shown in FIG. 17, the movement speed calculation unit 109 calculates the future movement direction and movement speed of the cloud from the history of the movement direction and movement speed of the cloud up to that point, including the movement direction and movement speed represented by the arrow P1. Then, the movement speed calculation unit 109 predicts positions 413 to 415, which are the future movement positions of the cloud.

その後、移動速度算出部109は、個別の雲や雲域の移動方向及び移動速度の履歴を気象情報データベース3へ送信して格納させる。また、移動速度算出部109は、個別の雲や雲域の移動速度及び移動方向の予測値、並びに、個別の雲や雲域の予測した位置の情報を気象情報データベース3へ送信して格納させる。 Then, the movement speed calculation unit 109 transmits the history of the movement direction and movement speed of each cloud or cloud area to the weather information database 3 for storage. In addition, the movement speed calculation unit 109 transmits the predicted values of the movement speed and movement direction of each cloud or cloud area, as well as information on the predicted position of each cloud or cloud area, to the weather information database 3 for storage.

情報活用基盤4は、個別の雲や雲域の移動方向及び移動速度、並びに、将来の位置等を用いて、データ加工を行う。そして、情報活用基盤4は、個別の雲や雲域の移動方向及び移動速度、並びに、将来の位置等を用いて加工した気象情報を各種気象情報活用システム7へ提供する。 The information utilization platform 4 processes data using the direction and speed of movement of individual clouds and cloud areas, as well as future positions, etc. Then, the information utilization platform 4 provides the weather information processed using the direction and speed of movement of individual clouds and cloud areas, as well as future positions, etc. to the various weather information utilization systems 7.

次に、図18を参照して、気象情報処理装置1による雲及び雲域の移動方向及び移動速度の算出処理の流れについて説明する。図18は、実施例2に係る気象情報処理装置による雲及び雲域の移動方向及び移動速度の算出処理のフローチャートである。 Next, the flow of the calculation process of the moving direction and moving speed of clouds and cloud areas by the weather information processing device 1 will be described with reference to FIG. 18. FIG. 18 is a flowchart of the calculation process of the moving direction and moving speed of clouds and cloud areas by the weather information processing device according to the second embodiment.

画像解析部101は、天球画像201及び気象データを地上センサ2から取得する(ステップS201)。 The image analysis unit 101 acquires a celestial image 201 and meteorological data from the ground sensor 2 (step S201).

また、画像解析部101は、気象衛星6により撮影された可視画像301及び赤外画像302を含む気象衛星画像を気象データ配信システム5から取得する(ステップS202)。 The image analysis unit 101 also acquires meteorological satellite images, including a visible image 301 and an infrared image 302, taken by a meteorological satellite 6 from the meteorological data distribution system 5 (step S202).

次に、画像解析部101は、天球画像201の画像解析を行い、各領域の雲のありなしの判定、個々の雲の識別、雲の種類の特定、個別の雲の時系列変化、並びに、個々の雲の地上センサ2からの方位角及び天頂角を含む天球画像解析情報を取得する。また、画像解析部101は、気象衛星画像の画像解析を行い、各領域の雲のありなしの判定、個々の雲の識別、雲の種類の特定及び個別の雲の時系列変化を含む衛星画像解析情報を取得する(ステップS203)。その後、画像解析部101は、天球画像201を変換部102へ出力する。また、画像解析部101は、可視画像301及び赤外画像302を含む気象衛星画像、衛星画像解析情報、天球画像解析情報、並びに、地上気象データを画像照合部104へ出力する。このとき、画像解析部101は、地上センサ2の緯度経度及び平面画像202の元となった天球画像203の撮影時刻も画像照合部104へ出力する。さらに、画像解析部101は、衛星画像解析情報及び天球画像解析情報を移動速度算出部109へ出力する。 Next, the image analysis unit 101 performs image analysis of the celestial image 201, and obtains celestial image analysis information including a determination of the presence or absence of clouds in each area, identification of each cloud, identification of the type of cloud, time series changes of each cloud, and the azimuth angle and zenith angle of each cloud from the ground sensor 2. The image analysis unit 101 also performs image analysis of the meteorological satellite image, and obtains satellite image analysis information including a determination of the presence or absence of clouds in each area, identification of each cloud, identification of the type of cloud, and time series changes of each cloud (step S203). The image analysis unit 101 then outputs the celestial image 201 to the conversion unit 102. The image analysis unit 101 also outputs the meteorological satellite image including the visible image 301 and the infrared image 302, the satellite image analysis information, the celestial image analysis information, and the ground meteorological data to the image matching unit 104. At this time, the image analysis unit 101 also outputs the latitude and longitude of the ground sensor 2 and the capture time of the celestial image 203 that was the basis of the planar image 202 to the image matching unit 104. Furthermore, the image analysis unit 101 outputs the satellite image analysis information and the celestial image analysis information to the movement speed calculation unit 109.

変換部102は、天球画像201の入力を画像解析部101から受ける。そして、変換部102は、天球画像201を平面展開して平面画像202に変換する(ステップS204)。その後、変換部102は、生成した平面画像202を画像照合部104へ出力する。 The conversion unit 102 receives the input of the celestial sphere image 201 from the image analysis unit 101. The conversion unit 102 then expands the celestial sphere image 201 onto a plane to convert it into a planar image 202 (step S204). The conversion unit 102 then outputs the generated planar image 202 to the image matching unit 104.

画像照合部104は、可視画像301及び赤外画像302を含む気象衛星画像、衛星画像解析情報、天球画像解析情報、地上気象データ、地上センサ2の緯度経度及び平面画像202の元となった天球画像203の撮影時刻の入力を画像解析部101から受ける。そして、画像照合部104は、画像照合解析処理を行うことで、可視画像301に映された雲の中から平面画像202に映された雲を同定する(ステップS205)。その後、画像照合部104は、画像照合解析処理を施した平面画像202及び可視画像301を移動速度算出部109へ出力する。 The image matching unit 104 receives inputs from the image analysis unit 101 of meteorological satellite images including the visible image 301 and the infrared image 302, satellite image analysis information, celestial image analysis information, ground meteorological data, the latitude and longitude of the ground sensor 2, and the capture time of the celestial image 203 that is the basis of the planar image 202. The image matching unit 104 then performs image matching analysis processing to identify clouds captured in the planar image 202 from among the clouds captured in the visible image 301 (step S205). The image matching unit 104 then outputs the planar image 202 and the visible image 301 that have been subjected to image matching analysis processing to the movement speed calculation unit 109.

移動速度算出部109は、衛星画像解析情報及び天球画像解析情報の入力を画像解析部101から受ける。また、移動速度算出部109は、画像照合解析処理を施した平面画像202及び可視画像301の入力を画像照合部104から受ける。そして、移動速度算出部109は、衛星画像解析情報、天球画像解析情報、並びに、画像照合解析処理を施した平面画像202及び可視画像301を用いて、個別の雲や雲域の時系列変化をトレースして移動方向及び移動速度の履歴を求める(ステップS206)。 The moving speed calculation unit 109 receives input of satellite image analysis information and celestial image analysis information from the image analysis unit 101. The moving speed calculation unit 109 also receives input of the planar image 202 and visible image 301 that have been subjected to image matching and analysis processing from the image matching unit 104. The moving speed calculation unit 109 then uses the satellite image analysis information, celestial image analysis information, and the planar image 202 and visible image 301 that have been subjected to image matching and analysis processing to trace the time-series changes of individual clouds and cloud areas to determine the history of the moving direction and moving speed (step S206).

次に、移動速度算出部109は、移動方向及び移動速度の履歴から、個別の雲及び雲域の将来の移動を予測する(ステップS207)。その後、移動速度算出部109は、予測した個別の雲及び雲域の将来の移動の情報を気象情報データベース3に格納させる。 Next, the movement speed calculation unit 109 predicts the future movement of individual clouds and cloud areas from the history of the movement direction and movement speed (step S207). After that, the movement speed calculation unit 109 stores information on the predicted future movement of individual clouds and cloud areas in the weather information database 3.

以上に説明したように、本実施例に係る気象情報処理装置は、衛星画像解析情報、天球画像解析情報、並びに、雲の同定を行った平面画像及び可視画像を用いて、個別の雲及び雲域の移動方向及び移動速度の履歴を求める。また、気象情報処理装置は、個別の雲及び雲域の移動方向及び移動速度の履歴から、個別の雲及び雲域の将来の移動を予測する。これらの情報を加えて加工を行った気象情報を各種気象情報活用システムへ提供することでより、詳細で利用価値の高い情報を提供することができ、気象データを活用するシステムの性能を向上させることが可能となる。 As described above, the weather information processing device according to this embodiment uses satellite image analysis information, celestial image analysis information, and planar and visible images with cloud identification to determine the history of the movement direction and speed of individual clouds and cloud areas. In addition, the weather information processing device predicts the future movement of individual clouds and cloud areas from the history of the movement direction and speed of individual clouds and cloud areas. By adding this information and providing processed weather information to various weather information utilization systems, it is possible to provide more detailed and useful information, thereby improving the performance of systems that utilize weather data.

(ハードウェア構成)
図19は、気象情報処理装置のハードウェア構成図である。気象情報処理装置1は、CPU(Central Processing Unit)91、メモリ92、ハードディスク93及びネットワークインタフェース94を有する。CPU91は、メモリ92、ハードディスク93及びネットワークインタフェース94とバスで接続される。
(Hardware configuration)
19 is a hardware configuration diagram of the weather information processing device. The weather information processing device 1 has a CPU (Central Processing Unit) 91, a memory 92, a hard disk 93, and a network interface 94. The CPU 91 is connected to the memory 92, the hard disk 93, and the network interface 94 via a bus.

ネットワークインタフェース94は、例えば、地上センサ2、気象情報データベース3及び気象データ配信システム5と通信を行うためのインタフェースである。 The network interface 94 is an interface for communicating with, for example, the ground sensor 2, the weather information database 3, and the weather data distribution system 5.

ハードディスク93は、図2及び16に例示した画像解析部101、変換部102、雲頂高度算出部103、画像照合部104、範囲取得部105、雲厚算出部106、雲底高度算出部107、様相識別部108及び移動速度算出部109の機能を実現するためのプログラムを含む各種プログラムを格納する。 The hard disk 93 stores various programs including programs for realizing the functions of the image analysis unit 101, conversion unit 102, cloud top height calculation unit 103, image matching unit 104, range acquisition unit 105, cloud thickness calculation unit 106, cloud base height calculation unit 107, appearance identification unit 108, and movement speed calculation unit 109 illustrated in Figures 2 and 16.

CPU91は、図2及び16に例示した、各種プログラムをハードディスク93から読み出し、メモリ92上に展開して実行する。これにより、CPU91及びメモリ92は、図2及び16に例示した、画像解析部101、変換部102、雲頂高度算出部103、画像照合部104、範囲取得部105、雲厚算出部106、雲底高度算出部107、様相識別部108及び移動速度算出部109の機能を実現する。 The CPU 91 reads out the various programs exemplified in Figs. 2 and 16 from the hard disk 93, and deploys and executes them on the memory 92. As a result, the CPU 91 and memory 92 realize the functions of the image analysis unit 101, conversion unit 102, cloud top height calculation unit 103, image matching unit 104, range acquisition unit 105, cloud thickness calculation unit 106, cloud base height calculation unit 107, appearance identification unit 108, and movement speed calculation unit 109 exemplified in Figs. 2 and 16.

1 気象情報処理装置
2 地上センサ
3 気象情報データベース
4 情報活用基盤
5 気象データ配信システム
6 気象衛星
7 各種気象情報活用システム
101 画像解析部
102 変換部
103 雲頂高度算出部
104 画像照合部
105 範囲取得部
106 雲厚算出部
107 雲底高度算出部
108 様相識別部
109 移動速度算出部
21 天球カメラ
22 気象データ取得部
REFERENCE SIGNS LIST 1 Weather information processing device 2 Ground sensor 3 Weather information database 4 Information utilization platform 5 Weather data distribution system 6 Meteorological satellite 7 Various weather information utilization systems 101 Image analysis unit 102 Conversion unit 103 Cloud top height calculation unit 104 Image matching unit 105 Range acquisition unit 106 Cloud thickness calculation unit 107 Cloud base height calculation unit 108 Aspect identification unit 109 Travel speed calculation unit 21 Celestial camera 22 Weather data acquisition unit

Claims (13)

天頂方向を中心とした平面の天球画像を撮影するカメラを有する地上センサ及び気象衛星で撮影された衛星画像を配信する気象データ配信システムから受信した、前記天球画像及び前記衛星画像を基に、前記天球画像に映った雲分布と前記衛星画像に映った雲分布との同定を行って画像照合を実行する画像照合部と、
前記画像照合の結果を基に、前記衛星画像の水平方向の距離情報から前記天球画像における水平方向の距離情報を取得する範囲取得部と、
前記天球画像における前記水平方向の距離情報及び前記天球画像における雲の天頂角を基に前記雲の雲底高度を算出する雲底高度算出部と
を備えたことを特徴とする気象情報処理装置。
an image matching unit that performs image matching by identifying the cloud distribution captured in the celestial image and the cloud distribution captured in the satellite image based on the celestial image and the satellite image received from a ground sensor having a camera that captures a planar celestial image centered in the zenith direction and a meteorological data distribution system that distributes satellite images captured by a meteorological satellite;
a range acquisition unit that acquires horizontal distance information in the celestial image from horizontal distance information in the satellite image based on a result of the image matching;
and a cloud base altitude calculation unit that calculates a cloud base altitude of the cloud based on the horizontal distance information in the celestial image and the zenith angle of the cloud in the celestial image.
前記天球画像又は前記衛星画像のいずれかもしくは双方を用いて画像解析を行い、積雲、積乱雲及び巻雲を含む雲の種類を特定する画像解析部と、
前記天球画像を平面画像に変換する変換部と
前記衛星画像及び地上気温を基に前記雲の雲頂高度を算出する雲頂高度算出部と、
前記雲底高度と前記雲頂高度との差から前記雲の厚さを算出する雲厚算出部と、
前記雲の厚さを基に雲の規模を含む雲の様相を識別する様相識別部とをさらに備え、
画像照合部は、前記画像解析部により特定された前記雲の種類を基に前記画像照合を実行する
ことを特徴とする請求項1に記載の気象情報処理装置。
an image analysis unit that performs image analysis using either or both of the celestial image and the satellite image and identifies cloud types including cumulus, cumulonimbus, and cirrus;
a conversion unit for converting the celestial image into a planar image; and a cloud-top height calculation unit for calculating the cloud-top height based on the satellite image and the ground temperature.
A cloud thickness calculation unit that calculates the thickness of the cloud from the difference between the cloud base altitude and the cloud top altitude;
A cloud aspect identification unit that identifies a cloud aspect including a cloud scale based on the cloud thickness,
The weather information processing apparatus according to claim 1 , wherein the image matching unit executes the image matching based on the type of cloud identified by the image analysis unit.
前記画像解析部により特定された雲の種類、前記画像照合部による前記画像照合の結果及び前記衛星画像を基に、雲の移動情報を取得する移動情報取得部をさらに備えたことを特徴とする請求項2に記載の気象情報処理装置。 The weather information processing device according to claim 2, further comprising a movement information acquisition unit that acquires cloud movement information based on the type of cloud identified by the image analysis unit, the result of the image matching by the image matching unit, and the satellite image. 前記雲厚算出部は、前記画像解析部により特定された前記雲の種類から前記雲の厚さの推定を行い、推定結果を用いて算出した前記雲の厚さの検証及び補正を行うことを特徴とする請求項2に記載の気象情報処理装置。 The weather information processing device according to claim 2, characterized in that the cloud thickness calculation unit estimates the cloud thickness from the type of cloud identified by the image analysis unit, and verifies and corrects the calculated cloud thickness using the estimation result. 前記範囲取得部は、前記画像解析部により特定された前記雲の種類を基に、個々の前記雲の種類に応じた高度に前記雲を分類し、分類毎に前記雲の水平方方向の距離を求めることを特徴とする請求項2~4のいずれか一つに記載の気象情報処理装置。 The weather information processing device according to any one of claims 2 to 4, characterized in that the range acquisition unit classifies the clouds into altitudes corresponding to the individual cloud types based on the cloud types identified by the image analysis unit, and determines the horizontal distance of the clouds for each classification. 前記地上センサは、地上の温度、気圧及び相対湿度を含む地上気象データを取得する気象データ取得部をさらに備え、
前記雲底高度算出部は、前記気象データ取得部により取得された前記地上気象データから持ち上げ凝結高度を求め、求めた前記持ち上げ凝結高度を基に前記雲底高度の推定を行い、推定結果を用いて前記雲底高度の検証及び補正を行うことを特徴とする請求項1~5のいずれか一つに記載の気象情報処理装置。
The ground sensor further includes a meteorological data acquisition unit that acquires ground meteorological data including a ground temperature, an air pressure, and a relative humidity,
The weather information processing device according to any one of claims 1 to 5, characterized in that the cloud base altitude calculation unit calculates a lifting condensation altitude from the ground meteorological data acquired by the meteorological data acquisition unit, estimates the cloud base altitude based on the calculated lifting condensation altitude, and verifies and corrects the cloud base altitude using the estimated result.
地上センサ、気象データ配信システム及び気象情報処理装置を有する気象情報処理システムであって、
前記地上センサは、
天頂方向を中心とした平面の天球画像を撮影するカメラを備え、
前記気象データ配信システムは、気象衛星で撮影された衛星画像を配信し、
前記気象情報処理装置は、
前記地上センサ及び前記気象データ配信システムから受信した、前記天球画像及び前記衛星画像を基に、前記天球画像に映った雲分布と前記衛星画像に映った雲分布との同定を行って画像照合を実行する画像照合部と、
前記画像照合の結果を基に、前記衛星画像の水平方向の距離情報から前記天球画像における水平方向の距離情報を取得する範囲取得部と、
前記天球画像における前記水平方向の距離情報及び前記天球画像における雲の天頂角を基に前記雲の雲底高度を算出する雲底高度算出部と
を備えた
ことを特徴とする気象情報処理システム。
A weather information processing system having a ground sensor, a weather data distribution system, and a weather information processing device,
The ground sensor includes:
A camera is provided for capturing a flat celestial image centered on the zenith direction,
The meteorological data distribution system distributes satellite images taken by a meteorological satellite,
The weather information processing device includes:
an image matching unit that performs image matching by identifying a cloud distribution shown in the celestial image and a cloud distribution shown in the satellite image based on the celestial image and the satellite image received from the ground sensor and the weather data distribution system;
a range acquisition unit that acquires horizontal distance information in the celestial image from horizontal distance information in the satellite image based on a result of the image matching;
a cloud base altitude calculation unit that calculates a cloud base altitude of the cloud based on the horizontal distance information in the celestial image and the zenith angle of the cloud in the celestial image.
前記気象情報処理装置は、
前記天球画像又は前記衛星画像のいずれかもしくは双方を用いて画像解析を行い、積雲、積乱雲及び巻雲を含む雲の種類を特定する画像解析部と、
前記天球画像を平面画像に変換する変換部と
前記衛星画像及び地上気温を基に前記雲の雲頂高度を算出する雲頂高度算出部と、
前記雲底高度と前記雲頂高度との差から前記雲の厚さを算出する雲厚算出部と、
前記雲の厚さを基に雲の規模を含む雲の様相を識別する様相識別部とをさらに備え、
画像照合部は、前記画像解析部により特定された前記雲の種類を基に前記画像照合を実行する
ことを特徴とする請求項7に記載の気象情報処理システム。
The weather information processing device includes:
an image analysis unit that performs image analysis using either or both of the celestial image and the satellite image and identifies cloud types including cumulus, cumulonimbus, and cirrus;
a conversion unit for converting the celestial image into a planar image; and a cloud-top height calculation unit for calculating the cloud-top height based on the satellite image and the ground temperature.
A cloud thickness calculation unit that calculates the thickness of the cloud from the difference between the cloud base altitude and the cloud top altitude;
A cloud aspect identification unit that identifies a cloud aspect including a cloud scale based on the cloud thickness,
The weather information processing system according to claim 7 , wherein the image matching unit executes the image matching based on the type of cloud identified by the image analysis unit.
前記地上センサは、地上の温度、気圧及び相対湿度を含む地上気象データを取得する気象データ取得部をさらに備え、
前記雲頂高度算出部は、前記気象データ取得部により取得された前記地上の温度を基に前記雲の雲頂高度を算出し、
前記雲底高度算出部は、前記気象データ取得部により取得された前記地上気象データから持ち上げ凝結高度を求め、求めた前記持ち上げ凝結高度を基に前記雲底高度の推定を行い、推定結果を用いて前記雲底高度の検証及び補正を行うことを特徴とする請求項8に記載の気象情報処理システム。
The ground sensor further includes a meteorological data acquisition unit that acquires ground meteorological data including a ground temperature, an air pressure, and a relative humidity,
the cloud top altitude calculation unit calculates a cloud top altitude of the cloud based on the ground temperature acquired by the meteorological data acquisition unit;
The weather information processing system according to claim 8, characterized in that the cloud base altitude calculation unit calculates a lifting condensation altitude from the ground meteorological data acquired by the meteorological data acquisition unit, estimates the cloud base altitude based on the calculated lifting condensation altitude, and verifies and corrects the cloud base altitude using the estimated result.
予め取得した前記雲の種類、前記雲底高度、前記雲頂高度、前記雲の厚さ及び前記雲の規模の情報を用いて機械学習を行い、前記機械学習の結果を基に、前記雲の種類、前記雲底高度、前記雲頂高度、前記雲の厚さ及び前記雲の規模の情報に含まれる情報のうちのいくつかの情報を用いて残りの情報を推定する気象情報生活基盤をさらに備えたことを特徴とする請求項8又は9に記載の気象情報処理システム。 The weather information processing system according to claim 8 or 9, further comprising a weather information infrastructure that performs machine learning using information on the cloud type, the cloud base altitude, the cloud top altitude, the cloud thickness, and the cloud scale that has been previously acquired, and estimates remaining information using some of the information contained in the information on the cloud type, the cloud base altitude, the cloud top altitude, the cloud thickness, and the cloud scale based on the results of the machine learning. 前記気象情報生活基盤は、前記雲の様相を他の気象情報と関連付けて加工して提供することを特徴とする請求項10に記載の気象情報処理システム。 The weather information processing system according to claim 10, characterized in that the weather information life infrastructure provides the cloud appearance by associating it with other weather information and processing it. 天頂方向を中心とした平面の天球画像を撮影するカメラを有する地上センサ及び気象衛星で撮影された衛星画像を配信する気象データ配信システムから受信した、前記天球画像及び前記衛星画像を基に、前記天球画像に映った雲分布と前記衛星画像に映った雲分布との同定を行って画像照合を実行し、
前記画像照合の結果を基に、前記衛星画像の水平方向の距離情報から前記天球画像における水平方向の距離情報を取得し、
前記天球画像における前記水平方向の距離情報及び前記天球画像における雲の天頂角を基に前記雲の雲底高度を算出する
ことを特徴とする気象情報処理方法。
Based on the celestial image and the satellite image received from a ground sensor having a camera for capturing a planar celestial image centered in the zenith direction and a meteorological data distribution system for distributing satellite images captured by a meteorological satellite, the cloud distribution captured in the celestial image is identified with the cloud distribution captured in the satellite image, and an image matching is performed;
acquiring horizontal distance information in the celestial image from horizontal distance information in the satellite image based on a result of the image matching;
13. A weather information processing method comprising: calculating a cloud base altitude of the cloud based on the horizontal distance information in the celestial image and a zenith angle of the cloud in the celestial image.
天頂方向を中心とした平面の天球画像を撮影するカメラを有する地上センサ及び気象衛星で撮影された衛星画像を配信する気象データ配信システムから受信した、前記天球画像及び前記衛星画像を基に、前記天球画像に映った雲分布と前記衛星画像に映った雲分布との同定を行って画像照合を実行し、
前記画像照合の結果を基に、前記衛星画像の水平方向の距離情報から前記天球画像における水平方向の距離情報を取得し、
前記天球画像における前記水平方向の距離情報及び前記天球画像における雲の天頂角を基に前記雲の雲底高度を算出する
処理をコンピュータに実行させることを特徴とする気象情報処理プログラム。
Based on the celestial image and the satellite image received from a ground sensor having a camera for capturing a planar celestial image centered in the zenith direction and a meteorological data distribution system for distributing satellite images captured by a meteorological satellite, the cloud distribution captured in the celestial image is identified with the cloud distribution captured in the satellite image, and an image matching is performed;
acquiring horizontal distance information in the celestial image from horizontal distance information in the satellite image based on a result of the image matching;
Calculating the cloud base altitude of the cloud based on the horizontal distance information in the celestial image and the zenith angle of the cloud in the celestial image.
JP2020147794A 2020-09-02 2020-09-02 Weather information processing device, weather information processing system, weather information processing method, and weather information processing program Active JP7520649B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020147794A JP7520649B2 (en) 2020-09-02 2020-09-02 Weather information processing device, weather information processing system, weather information processing method, and weather information processing program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020147794A JP7520649B2 (en) 2020-09-02 2020-09-02 Weather information processing device, weather information processing system, weather information processing method, and weather information processing program

Publications (2)

Publication Number Publication Date
JP2022042380A JP2022042380A (en) 2022-03-14
JP7520649B2 true JP7520649B2 (en) 2024-07-23

Family

ID=80629462

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020147794A Active JP7520649B2 (en) 2020-09-02 2020-09-02 Weather information processing device, weather information processing system, weather information processing method, and weather information processing program

Country Status (1)

Country Link
JP (1) JP7520649B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2025203540A1 (en) * 2024-03-29 2025-10-02 三菱電機株式会社 Cloud cover forecasting device, power generation amount forecasting device, cloud cover forecasting method, and power generation amount forecasting method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007184354A (en) 2006-01-05 2007-07-19 Mitsubishi Electric Corp Solar power system
JP2013529051A (en) 2010-05-07 2013-07-11 アドバンスド エナージィ インダストリーズ,インコーポレイテッド Photovoltaic power generation prediction system and method
JP2014048130A (en) 2012-08-30 2014-03-17 Toshiba Corp Image processing device, method, and program
JP2015059821A (en) 2013-09-18 2015-03-30 株式会社東芝 Amount-of-insolation predicting device and method
JP2019015517A (en) 2017-07-03 2019-01-31 株式会社東芝 Solar radiation intensity estimation device, solar radiation intensity estimation system, and solar radiation intensity estimation method
WO2020080373A1 (en) 2018-10-15 2020-04-23 国立研究開発法人宇宙航空研究開発機構 Optical ground station operational management system, optical operation planning device, and optical ground station operational management method and program

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6035710A (en) * 1999-05-28 2000-03-14 Lockheed Martin Missiles & Space Co. Cloud base height and weather characterization, visualization and prediction based on satellite meteorological observation

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007184354A (en) 2006-01-05 2007-07-19 Mitsubishi Electric Corp Solar power system
JP2013529051A (en) 2010-05-07 2013-07-11 アドバンスド エナージィ インダストリーズ,インコーポレイテッド Photovoltaic power generation prediction system and method
JP2014048130A (en) 2012-08-30 2014-03-17 Toshiba Corp Image processing device, method, and program
JP2015059821A (en) 2013-09-18 2015-03-30 株式会社東芝 Amount-of-insolation predicting device and method
JP2019015517A (en) 2017-07-03 2019-01-31 株式会社東芝 Solar radiation intensity estimation device, solar radiation intensity estimation system, and solar radiation intensity estimation method
WO2020080373A1 (en) 2018-10-15 2020-04-23 国立研究開発法人宇宙航空研究開発機構 Optical ground station operational management system, optical operation planning device, and optical ground station operational management method and program

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Takashi Y. Nakajima et al.,Theoretical basis of the algorithms and early phase results of the GCOM-C(Shikisai) SGLI cloud products,Progress in Earth and Planetary Science,2019年07月22日,Vol.6,No.52,pp.1-25,doi.org/10.1186/s40645-019-0295-9

Also Published As

Publication number Publication date
JP2022042380A (en) 2022-03-14

Similar Documents

Publication Publication Date Title
Rodríguez-Benítez et al. Assessment of new solar radiation nowcasting methods based on sky-camera and satellite imagery
Chow et al. Intra-hour forecasting with a total sky imager at the UC San Diego solar energy testbed
US10989839B1 (en) Ground-based sky imaging and irradiance prediction system
Menzel et al. Introducing GOES-I: The first of a new generation of geostationary operational environmental satellites
JP6680788B2 (en) Detecting apparatus and method for identifying and monitoring clouds in the observation region of the sky
Nguyen et al. Stereographic methods for cloud base height determination using two sky imagers
US20170031056A1 (en) Solar Energy Forecasting
JP6599935B2 (en) Solar radiation intensity estimation device, solar radiation intensity estimation system, and solar radiation intensity estimation method
Dev et al. Estimation of solar irradiance using ground-based whole sky imagers
Meier et al. Atmospheric correction of thermal-infrared imagery of the 3-D urban environment acquired in oblique viewing geometry
JP6925619B2 (en) How to measure cloud altitude and wind speed using optical images
Seiz et al. Cloud mapping from the ground: Use of photogrammetric methods
CN113076865A (en) Method and system for inverting irradiance based on sky photographing image and satellite cloud image
Liandrat et al. Cloud cover forecast from a ground-based all sky infrared thermal camera
US9734396B2 (en) Method of refraction angle mapping of earth'S limb
Seiz et al. Cloud mapping with ground‐based photogrammetric cameras
Sánchez-Segura et al. Solar irradiance components estimation based on a low-cost sky-imager
JP7520649B2 (en) Weather information processing device, weather information processing system, weather information processing method, and weather information processing program
Öktem et al. Stereophotogrammetry of oceanic clouds
Liu et al. High-spatial-resolution nighttime light dataset acquisition based on volunteered passenger aircraft remote sensing
CN109945969B (en) Method and device for determining earth radiation balance based on meteorological satellite observation
CN115933007A (en) Weather monitoring method and device based on multi-temporal and spatial-scale satellite meteorology
Blum Nowcasting of solar irradiance and photovoltaic production using a network of all-sky imagers
CN118411399A (en) A cloud height estimation method and system based on ground-based cloud image data
Terrén-Serrano et al. Geospatial perspective reprojections for ground-based sky imaging systems

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230705

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20240214

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240220

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240416

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: 20240618

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240710

R150 Certificate of patent or registration of utility model

Ref document number: 7520649

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150