JP2017163934A - Methods of calculating damage classification of pine wilt disease and apparatus for calculating damage classification of pine wilt disease - Google Patents

Methods of calculating damage classification of pine wilt disease and apparatus for calculating damage classification of pine wilt disease Download PDF

Info

Publication number
JP2017163934A
JP2017163934A JP2016054551A JP2016054551A JP2017163934A JP 2017163934 A JP2017163934 A JP 2017163934A JP 2016054551 A JP2016054551 A JP 2016054551A JP 2016054551 A JP2016054551 A JP 2016054551A JP 2017163934 A JP2017163934 A JP 2017163934A
Authority
JP
Japan
Prior art keywords
tree
damage
image data
pine
forest
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.)
Granted
Application number
JP2016054551A
Other languages
Japanese (ja)
Other versions
JP6544582B2 (en
Inventor
加藤 正人
Masato Kato
正人 加藤
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.)
Shinshu University NUC
Original Assignee
Shinshu University NUC
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 Shinshu University NUC filed Critical Shinshu University NUC
Priority to JP2016054551A priority Critical patent/JP6544582B2/en
Publication of JP2017163934A publication Critical patent/JP2017163934A/en
Application granted granted Critical
Publication of JP6544582B2 publication Critical patent/JP6544582B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide methods of calculating damage classification of pine wilt disease and apparatus for calculating damage classification of pine wilt disease which make the production of materials showing the damage situation of pine wilt disease in a Japanese red pine forest achievable in a short time and with high accuracy.SOLUTION: A method of calculating damage classification of pine wilt disease in which the number of Japanese red pines are calculated by the damage classification using the image data of a Japanese red pine forest produced based on spectral reflection properties of tree species included in photographed image data comprises a Japanese red pine forest tree apex imaging data production step S31 of producing imaging data of Japanese red pine forest tree apex, a Japanese red pine forest tree crown imaging data production step S32 of producing imaging data of Japanese red pine forest tree crown, a tree crown specific damage classification image data production step S33 of producing imaging data of damage classification by tree crown, a tree apex specific damage classification imaging data production step S34, a damage classification specific number calculating step S40 of calculating the number of Japanese red pine trees by a damage classification, and a damage ratio classification figure production step S50 of producing a damage ratio classification figure.SELECTED DRAWING: Figure 2

Description

本発明は、松くい虫の被害区分算定方法及び松くい虫の被害区分算定装置に関する。   The present invention relates to a method for calculating a damage category for pine worms and a device for calculating a damage category for pine worms.

松くい虫の被害は我が国最大の森林病虫害被害である。全国の松くい虫被害量は昭和54年をピークに減少傾向にあるが、高緯度・高標高地域にある長野県や東北地方においては被害量が増加し、森林被害が深刻化している。被害のまん延により、森林所有者の林業意欲とマツタケ生産の減退、枯れ木・倒木による危険性、道路法面の防災機能の低下、地域の松林や社寺林及び公園や環境林などの公益的機能の低下など、その対策は極めて重要な課題である。   The damage of pine weevil is the largest damage caused by forest pests in Japan. The amount of pine wilt damage nationwide has been decreasing since its peak in 1979, but the amount of damage has increased in Nagano Prefecture and Tohoku region in high latitudes and high altitude areas, and forest damage has become serious. Due to the spread of damage, forest owners' willingness to forestry and matsutake production decline, danger due to dead trees and fallen trees, reduced road disaster prevention function, local pine forests and shrines and temples, public facilities such as parks and environmental forests Countermeasures such as decline are extremely important issues.

松くい虫の被害と樹葉の変化は、マツノザイセンチュウという線虫が付着したマツノマダラカミキリが春に羽化して、健全なマツ(アカマツとする。)の若枝を食べる時にマツノザイセンチュウが樹体内に入り、夏に増殖し、秋から冬に感染した松が衰弱し、葉が黄緑から赤褐色に変色し枯損する。翌年の夏以降、色が薄れて茶色になり、やがて落葉する。マツノマダラカミキリは産卵するために枯損木を利用し、4〜6月にカミキリが蛹になると、マツノザイセンチュウは蛹の体に付着し、羽化の際に共生して移動することで別の健全なアカマツへと他力運搬され、被害が拡大していく。   The damage caused by the pine weevil and the changes in the leaves of the pinewood nematode are as follows. In the summer, it grows in the summer, the infected pine is weakened from autumn to winter, and the leaves turn from yellowish green to reddish brown and wither. From the summer of the following year, the color fades and turns brown and eventually leaves. The pine wood beetle uses dead trees to lay eggs, and when the beetle becomes a cocoon in April-June, the pine wood nematode attaches to the body of the cocoon and moves symbiotically when it emerges. The damage is spread by being transported to other pine trees.

松くい虫対策については、林野庁から提案されている(例えば、非特許文献1参照。)。非特許文献1によれば、公益的機能の高い松林を「保全すべき松林」、その周辺に位置する松林を「周辺松林」として都道府県知事等が定め、これらのマツ林を対象として重点的かつ総合的に対策を実施することとしている。また、現地での対策として、被害の状況を踏まえ、保全すべき松林において、薬剤等による「予防対策」や被害木の伐倒燻蒸等の「駆除対策」等を実施するとともに周辺松林では被害区分転換を推進することとすることとしている。   A countermeasure against pine weevil has been proposed by the Forestry Agency (see, for example, Non-Patent Document 1). According to Non-Patent Document 1, prefectural governors, etc., define pine forests with high public interest functions as “pine forests to be preserved” and pine forests located in the vicinity as “peripheral pine forests”, and focus on these pine forests. And comprehensive measures will be taken. In addition, as a local measure, based on the damage situation, in the pine forests to be preserved, “preventive measures” such as drugs and “extermination measures” such as felling and fumigating the damaged trees, etc. It is supposed to promote conversion.

上記した対策は、地上調査で調査員が被害の様子を見て行う必要があるが、松くい虫による森林被害対策を進めている都道府県や市町村においては、予算的にも人的にも対応しきれず、保全すべきマツ林と周辺マツ林のみの限定的な防除処理(薬剤散布と伐倒燻蒸処理)に留まっており、被害の進行に対応できず、被害のまん延が進んでいる。被害の最前線で対応する市町村は職員が少人数であり、しかも、専門の技術職員が配置されていないことが多いため、広域の森林域や道から離れた被害森林の様子を把握することは困難である。   The above measures need to be taken by the investigator to see how the damage occurred during ground surveys, but in prefectures and municipalities where measures against forest damage caused by pine weevil are being promoted, both budgetary and personal It is not limited, and only limited pine forests to be preserved and surrounding pine forests are limited (chemical spraying and felling and fumigation treatment). The municipalities that deal with the forefront of damage have a small number of staff, and there are often no specialized technical staff, so it is not possible to grasp the state of damaged forests away from wide area forest areas and roads Have difficulty.

その一方で、市町村単位で松くい虫被害の詳細な被害分布をつかむこと、被害の拡大を防ぐために、虫くいの被害を受けている感染木及び枯損木の集団をとらえて的確に防除することが望まれている。しかし、予算などが限られた状況においては、住宅近接地や道路沿いの松林などの狭い範囲での対策に限定されており、市町村単位での広域を対象とした的確な防除と対策ができていないという問題がある。   On the other hand, in order to grasp the detailed distribution of pine worm damage by municipalities and to prevent the spread of the damage, it is necessary to capture and accurately control the group of infected and dead trees that have been damaged by wormwood. It is desired. However, in a situation where the budget is limited, it is limited to measures in a narrow area such as a residential neighborhood or a pine forest along the road, and accurate control and measures for a wide area on a municipal level basis have been achieved. There is no problem.

調査委員などによる調査に頼らずに、感染木や枯損木の集団をとらえて的確に防除するためには、調査対象となるアカマツ林の現況を情報処理技術によって把握して、当該アカマツ林における松くい虫の被害状況を視覚的に容易に理解できる資料の作成を行えるようにすることが必要となる。なお、「感染木」というのは、当年生被害木であって、葉が黄緑に変色し元気のないアカマツを指しており、また、「枯損木」というのは、前年度被害木であって、葉が茶色になって枯れたアカマツを指している。   In order to capture and accurately control infected and dead tree groups without relying on surveys by survey committee members, etc., the current status of the Pinus densiflora forest to be surveyed is ascertained using information processing technology, and pine trees in the Pinus densiflora forest. It is necessary to be able to create materials that can easily understand the damage situation of caterpillars visually. “Infected trees” refers to damaged trees of the current year, whose leaves turn yellowish green, and are not energetic, and “dead trees” are damaged trees of the previous year. This refers to red pine with brown leaves.

一方、調査対象となるアカマツ林の所定範囲の森林地域の現況を情報処理技術によって把握する技術として、森林地域を上空から撮影して得られた撮影画像データに基づいて取得した森林画像情報から、森林の現況を把握するための技術が提案されている(例えば、特許文献1参照。)。   On the other hand, as a technique for grasping the current state of the forest area of the predetermined range of red pine forest to be investigated by information processing technology, from the forest image information obtained based on the captured image data obtained by photographing the forest area from the sky, Techniques for grasping the current state of forests have been proposed (see, for example, Patent Document 1).

特許文献1は、アカマツ林に限られるものではなく調査対象となる森林を含む地域を上空から撮影して得られた撮影画像データと森林の地理的情報とに基づいて抽出された森林域に関する森林域画像データに基づいて樹頂点に関するデータを抽出して、樹木の本数を高い精度で算定するものである。   Patent Document 1 is not limited to red pine forests, but forests related to forest areas extracted based on photographed image data obtained by photographing regions including forests to be surveyed from above and geographical information of the forests. Data on tree vertices is extracted based on area image data, and the number of trees is calculated with high accuracy.

林野庁 松くい虫被害、松くい虫対策について、[平成28年2月6日検索]、インターネット<URL:<http://www.rinya.maff.go.jp/j/hogo/higai/matukui.html>>Forestry Agency About pine weevil damage and pine weevil countermeasures [Search February 6, 2016], Internet <URL: <http://www.rinya.maff.go.jp/j/hogo/higai/matukui. html >>

特許第4858793号公報Japanese Patent No. 4858793

上述の非特許文献1に記載の松くい虫対策では、調査員が調査対象となるアカマツ林を歩いて被害の様子を実際に見て回ることによって松くい虫の被害状況を把握して、調査員が被害状況を示す資料の作成を行う必要があり、人件費などの多くの費用を要するとともに調査に多くの時間を要するといった課題がある。また、特許文献1に記載の樹木本数算定方法では、情報処理技術によって、樹木の本数を高い精度で算定することは可能であるが、調査対象となるアカマツ林における松くい虫の被害状況を示す資料の作成はできない。   In the pine worm countermeasures described in Non-Patent Document 1, the investigator walks through the red pine forest to be surveyed and looks around the state of the damage to find out the damage situation of the pine worm. There is a problem that it is necessary for the staff to prepare a document indicating the damage situation, which requires a lot of expenses such as labor costs and a lot of time for the investigation. In addition, in the method for calculating the number of trees described in Patent Document 1, it is possible to calculate the number of trees with high accuracy by using information processing technology, but it shows the damage situation of pine weevil in the red pine forest to be investigated. You cannot create materials.

そこで本発明は、アカマツ林における松くい虫の被害状況を示す資料の作成を安価な費用で短時間にかつ高精度に実現可能とする松くい虫の被害区分算定方法及び松くい虫の被害区分算定装置を提供することを目的とする。   Therefore, the present invention provides a method for calculating the damage category of pine worms and the damage category of pine worms, which enables the creation of data indicating the damage status of pine worms in a red pine forest to be realized in a short time and with high accuracy at a low cost. The purpose is to provide a calculation device.

本発明の発明者は、実在するアカマツ林において虫くいの被害を受けている感染木、枯損木及び虫くいの被害を受けていない健全なアカマツ(健全木という。)の各樹冠から感染葉、枯損葉、健全葉を採取して、人工衛星や航空機センサと同じ観測波長を計測できる分光反射計で計測したところ、感染葉、枯損葉及び健全葉の分光反射特性に違いがあることから、被害区分(感染木、枯損木及び健全木)を抽出することができることを発見した。本発明の発明者は、この点に着目して、鋭意研究を行うことにより本発明を完成させた。本発明は、以下の要素からなる。   The inventor of the present invention has an infected leaf, a dead leaf from each crown of an infected tree, a dead tree, and a healthy red pine tree (a healthy tree) that is not damaged by the insect bite in a real red pine forest. We collected healthy leaves and measured them with a spectroreflectometer that can measure the same observation wavelength as an artificial satellite or an aircraft sensor. Because there are differences in the spectral reflection characteristics of infected leaves, dead leaves, and healthy leaves, damage classification ( It was discovered that infected trees, dead trees, and healthy trees can be extracted. The inventor of the present invention paid attention to this point and completed the present invention by conducting intensive studies. The present invention comprises the following elements.

[1]本発明の松くい虫の被害区分算定方法は、調査対象となるアカマツ林を含む地域を上空から撮影して得られた撮影画像データに含まれる樹種の分光反射特性に基づいて作成されたアカマツ林画像データを用いて、前記アカマツ林に含まれるアカマツの本数を被害区分別に算定する松くい虫の被害区分算定方法であって、前記アカマツ林画像データに基づいて、前記アカマツの各樹頂点を抽出し、前記各樹頂点に関するアカマツ林樹頂点画像データを作成するアカマツ林樹頂点画像データ作成ステップと、前記アカマツ林画像データに基づいて、前記アカマツの各樹冠を抽出し、前記各樹冠に関するアカマツ林樹冠画像データを作成するアカマツ林樹冠画像データ作成ステップと、前記アカマツ林樹冠画像データと、松くい虫の被害を受けているアカマツから採取された被害葉及び健全なアカマツから採取された健全葉の各分光反射特性とに基づいて、前記各樹冠の領域における感染木、枯損木及び健全木の分光反射特性を作成し、作成した感染木、枯損木及び健全木の分光反射特性に基づいて樹冠別の被害区分を特定して、樹冠別被害区分画像データを作成する樹冠別被害区分画像データ作成ステップと、前記樹冠別被害区分画像データと前記アカマツ林樹頂点画像データとに基づいて、樹頂点別に前記被害区分を特定して、樹頂点別被害区分画像データを作成する樹頂点別被害区分画像データ作成ステップと、前記樹頂点別被害区分画像データと地理的情報とに基づいて、前記アカマツ林の全領域又は任意の領域に含まれるアカマツの本数を前記被害区分ごとに算定する被害区分本数算定ステップと、前記被害区分本数算定結果に基づいて、当該アカマツ林の全領域又は任意の領域に含まれる感染木と枯損木とを合計した被害木本数を、当該アカマツ林の全領域又は任意の領域に含まれるアカマツの全体本数で除することで得られる被害率を計算して被害率区分図を作成する被害率区分図作成ステップと、を有すること特徴とする。   [1] The method for calculating the damage classification of pine worms according to the present invention is created based on spectral reflectance characteristics of tree species included in photographed image data obtained by photographing an area including a red pine forest to be investigated from the sky. A method for calculating a damage category of a pine pine beetle by calculating the number of Pinus densiflora included in the Pinus densiflora forest according to the damage category using the Pinus densiflora image data, and each tree of the Pinus densiflora based on the Pinus densiflora image data A pine forest tree vertex image data creating step for extracting vertices and creating a red pine forest tree vertex image data for each tree vertex, and extracting each pine tree crown based on the pine forest image data; The pine forest canopy image data creating step for creating the red pine forest canopy image data, and the pine forest canopy image data, Spectral reflection characteristics of infected trees, dead trees and healthy trees in the area of each crown based on the spectral reflection characteristics of damaged leaves collected from red pine and healthy leaves collected from healthy red pine, The damage classification image data creation step by crown for identifying the damage classification by crown based on the spectral reflection characteristics of the created infected tree, dead tree and healthy tree, and creating the damage classification image data by crown, and the damage by the crown Based on the classification image data and the red pine forest tree vertex image data, the damage classification image data for each tree vertex is specified by identifying the damage classification for each tree vertex and creating the damage classification image data for each tree vertex; Damage zone that calculates the number of red pine contained in the whole area of the pine forest or in any area based on the image data and geographical information by vertex Based on the branch number calculation step and the damage category number calculation result, the total number of damaged trees in the entire area of the red pine forest or the infected trees and dead trees included in any area is calculated as the total area of the red pine forest or A damage rate classification diagram creation step of creating a damage rate classification diagram by calculating a damage rate obtained by dividing by the total number of red pine trees included in an arbitrary area.

本発明においては、調査対象となる森林を含む地域を上空から撮影して得られた撮影画像データに含まれる樹種の分光反射特性に基づいて抽出されたアカマツ林に関するアカマツ林画像データを用いて、上記した各ステップを実施することによって、実在のアカマツ林における松くい虫の被害区分を算定して、被害率区分図の作成を可能としている。なお、当該被害率区分図は、松くい虫の被害状況を示す資料として使用することができる。このように、本発明によれば、調査対象となるアカマツ林の現況を情報処理技術によって把握することができるため、アカマツ林における松くい虫の被害状況を示す資料の作成を安価な費用で短時間にかつ高精度に実現できる。   In the present invention, using red pine forest image data relating to red pine forest extracted based on the spectral reflection characteristics of the tree species included in the photographed image data obtained by photographing the region including the forest to be investigated from the sky, By implementing each of the above steps, it is possible to calculate the damage category of pine worms in the real red pine forest and to create a damage rate classification diagram. Note that the damage rate classification chart can be used as a material indicating the damage situation of pine weevil. As described above, according to the present invention, the current state of the red pine forest to be investigated can be grasped by the information processing technology. Therefore, it is possible to shorten the creation of the material indicating the pine weevil damage state in the red pine forest at a low cost. Realized in time and with high accuracy.

また、本発明において作成される被害率区分図は、被害率の凡例を色分けすることもできるので、被害が生じている区域を一目で把握できる。当該被害率区分図は、松くい虫による森林被害対策を進めている都道府県や市町村が求めている出力図である。市町村単位で松くい虫被害の詳細な被害分布をつかむことにおいて有効であるとともに、被害の拡大を防ぐため最前線の感染木や枯損木の集団をとらえ、的確に防除する上での計画立案や森林管理に有効な機能を有している。   In addition, since the damage rate classification diagram created in the present invention can be color-coded in the legend of the damage rate, the area where the damage is occurring can be grasped at a glance. The damage rate classification chart is an output chart required by prefectures and municipalities that are promoting countermeasures against forest damage caused by pine needles. It is effective in grasping the detailed distribution of pine weevil damage at the municipal level, and in order to prevent the spread of damage, it is necessary to plan for the prevention and proper control of the infected and dead tree groups at the forefront. It has an effective function for forest management.

なお、この明細書において、「樹冠」というのは、森林を構成する各樹木を上空からほぼ垂直方向に見た場合に、各樹木の枝と葉で構成される部分を指しており、周囲に他の樹木がない場合には、一本の樹木の枝と葉が1つの樹冠として抽出される場合もあるが、複数の樹木が近接して存在する場合には、複数の樹木の枝と葉が重なり合うため、複数の樹木の枝と葉で構成される広い範囲が1つの樹冠として抽出される(後述する図8参照。)。また、「樹頂点」というのは、樹冠内で尖った先端部分であり、これは、各樹木の頂点の部分に対応する(後述する図8参照。)。   In this specification, the term “tree crown” refers to a portion composed of branches and leaves of each tree when the trees constituting the forest are viewed from above in a substantially vertical direction. If there are no other trees, the branches and leaves of one tree may be extracted as one crown, but if there are multiple trees in close proximity, the branches and leaves of multiple trees Therefore, a wide range composed of a plurality of tree branches and leaves is extracted as one crown (see FIG. 8 described later). Further, the “tree vertex” is a pointed tip portion in the tree crown, and this corresponds to the vertex portion of each tree (see FIG. 8 described later).

[2]本発明の松くい虫の被害区分算定方法においては、前記撮影画像データに含まれる前記樹種の分光反射特性に基づいて、前記撮影画像データからアカマツ林を抽出することによってアカマツ林画像データを作成するアカマツ林画像データ作成ステップをさらに有し、当該アカマツ林画像データ作成ステップを実施した後に、前記アカマツ林ごとに、前記アカマツ林樹頂点画像データ作成ステップと、前記アカマツ林樹冠画像データ作成ステップと、前記被害区分画像データ作成ステップと、前記樹頂点別被害区分画像データ作成ステップと、前記被害区分本数算定ステップと、前記被害率区分図作成ステップとを実施することが好ましい。   [2] In the method for calculating the damage category of pine beetle of the present invention, red pine forest image data is extracted by extracting red pine forest from the photographed image data based on the spectral reflection characteristics of the tree species included in the photographed image data. A red pine forest image data creation step, and after performing the red pine forest image data creation step, for each red pine forest, the red pine forest tree vertex image data creation step and the red pine forest crown image data creation It is preferable to perform a step, the damage category image data creation step, the tree vertex-specific damage category image data creation step, the damage category number calculation step, and the damage rate category diagram creation step.

これにより、アカマツ林画像データを作成することができ、作成された当該アカマツ林画像データを用いて、上記各ステップを実施することができる。なお、アカマツ林画像データ作成ステップは、調査対象となるアカマツ林を含む地域を上空から撮影して得られた撮影画像データに含まれる反射輝度値が統計的に近似のものをクラス分けして行うものであり、公知の最尤法による土地被覆分類を用いて自動的に抽出することができる。   Thereby, red pine forest image data can be created, and the above steps can be performed using the created red pine forest image data. The red pine forest image data creation step is performed by classifying the reflection luminance values included in the photographed image data obtained by photographing the region including the red pine forest to be investigated from the sky into classes. It can be automatically extracted using land cover classification by the known maximum likelihood method.

[3]本発明の松くい虫の被害区分算定方法においては、前記撮影画像データは、秋から冬に撮影された画像データを用いることが好ましい。   [3] In the pine weevil damage classification calculation method of the present invention, it is preferable to use image data taken from autumn to winter as the photographed image data.

このように、秋から冬に撮影された画像データを用いることで、調査対象となるアカマツ林に落葉樹が混在している場合に、当該落葉樹は落葉しているため、被害区分すなわち感染木、枯損木及び健全木の区分がし易くなり、被害区分が高い精度で可能となる。   In this way, by using image data taken from autumn to winter, when deciduous trees are mixed in the red pine forest to be investigated, the deciduous trees are deciduous. It becomes easy to classify trees and healthy trees, and damage classification is possible with high accuracy.

[4]本発明の松くい虫の被害区分算定方法においては、前記感染木、枯損木及び健全木のうちの前記感染木と前記枯損木との区分は、前記松くい虫の被害を受けているアカマツから採取された被害葉と健全なアカマツから採取された健全葉の各分光反射特性に基づいて行い、前記感染木は、可視光の赤色域において前記枯損木及び前記健全木に比べて反射率が高いことで、前記枯損木及び前記健全木と区分し、前記枯損木は、近赤外域において前記感染木及び前記健全木に比べて反射率が低いことで、前記感染木及び前記健全木と区分する、ことが好ましい。   [4] In the method for calculating the damage classification of pine worms according to the present invention, the classification of the infected tree and the dead tree among the infected tree, the dead tree and the healthy tree is affected by the pine worm. The infected tree is reflected in the red region of visible light compared to the dead tree and the healthy tree. It is distinguished from the dead tree and the healthy tree by a high rate, and the dead tree has a lower reflectance than the infected tree and the healthy tree in the near infrared region, and thus the infected tree and the healthy tree. It is preferable to categorize.

このように、感染木、枯損木及び健全木の区分を、被害葉及び健全葉の各分光反射特性に基づいて行うことにより、感染木と枯損木と健全木との区分を高精度に行うことができ、特に、感染木と枯損木との区分を高精度に行うことができる。なお、「被害葉」というのは、感染木から採取された感染葉と枯損木から採取された枯損葉の両方を指すものである。   In this way, by classifying infected trees, dead trees, and healthy trees based on the spectral reflection characteristics of damaged leaves and healthy leaves, it is possible to accurately classify infected trees, dead trees, and healthy trees. In particular, it is possible to accurately distinguish between infected trees and dead trees. The “damaged leaf” refers to both an infected leaf collected from an infected tree and a dead leaf collected from a dead tree.

[5]本発明の松くい虫の被害区分算定方法においては、前記アカマツ林樹頂点画像データ作成ステップは、輝度値が極大値を持つ点を自動抽出する局所最大値フィルターを用いてアカマツの樹頂点を抽出することが好ましい。   [5] In the method for calculating the damage category of a pine pine beetle according to the present invention, the step of creating the pine forest tree vertex image data uses a local maximum value filter that automatically extracts a point having a maximum luminance value. It is preferable to extract vertices.

これにより、アカマツの樹頂点を精度良く抽出できる。なお、アカマツの樹頂点は、太陽光を強く反射する樹冠の尖った部分であり、前記アカマツ林画像データにおいて、反射値が最大になることから、局所最大値フィルター処理で自動的に抽出することができる。また、このような処理を行うことにより、全森林域から様々な樹冠の樹頂点を精度良く抽出できる。   As a result, the tree vertices of Pinus densiflora can be extracted with high accuracy. Note that the pinus of the red pine is the pointed portion of the crown that strongly reflects sunlight, and since the reflection value is maximum in the red pine forest image data, it is automatically extracted by local maximum value filtering. Can do. In addition, by performing such processing, tree vertices of various crowns can be extracted from the entire forest area with high accuracy.

[6]本発明の松くい虫の被害区分算定方法においては、前記アカマツ林樹冠画像データ作成ステップは、樹冠の縁部が樹冠部分よりも暗いことを利用して樹冠を抽出することが好ましい。   [6] In the pine beetle damage classification calculation method according to the present invention, it is preferable that the step of creating the red pine forest canopy image data extracts the crown using the fact that the edge of the crown is darker than the crown portion.

これにより、樹冠(アカマツ樹冠)を精度良く抽出できる。なお、樹冠の抽出は、前記アカマツ林画像データにおいて、樹冠の縁部(樹冠縁ともいう)が樹冠部分よりも暗いことを利用して行うものであり、これは、ブァレイフォローイング・アルゴリズム(Va11ey fo11owing a1gorithm)と自動マスク処理で自動的に抽出することにより、アカマツ林画像データから樹冠を精度良く抽出できる。なお、「樹冠の縁部」及び「樹冠部分」については後述の図8において説明する。   Thereby, a canopy (red pine canopy) can be extracted with high accuracy. In addition, the extraction of the tree crown is performed by utilizing the fact that the edge of the tree crown (also referred to as the tree crown edge) is darker than the tree crown part in the red pine forest image data, and this is represented by the valley following algorithm (Va11ey fo11owing a1gorithm) and automatic mask processing, tree crowns can be extracted accurately from Japanese red pine forest image data. Note that “the edge of the crown” and “the crown portion” will be described later with reference to FIG.

[7]本発明の松くい虫の被害区分算定方法においては、前記樹冠別被害区分画像データ作成ステップは、前記被害区分を感染木、枯損木及び健全木の3種類として、前記アカマツの各樹冠から感染木、枯損木及び健全木が抽出された場合、前記各樹冠から抽出された感染木、枯損木及び健全木のそれぞれの占有面積を算出し、算出された感染木、枯損木及び健全木の各占有面積のうち、占有面積が大きい被害区分を当該樹冠の被害区分とみなすことが好ましい。   [7] In the method for calculating the damage category of a pine beetle according to the present invention, the damage category image data creation step for each crown includes the three types of the damage category as an infected tree, a dead tree, and a healthy tree. When infected trees, dead trees, and healthy trees are extracted from the above, the occupied areas of the infected trees, dead trees, and healthy trees extracted from the respective crowns are calculated, and the calculated infected trees, dead trees, and healthy trees are calculated. Of these occupying areas, it is preferable to regard a damage category having a large occupancy area as a damage category of the crown.

これにより、1つの樹冠を1つの被害区分として区分することができる。例えば、ある樹冠において、感染木、枯損木及び健全木の各占有面積を算出した結果、当該樹冠においては、感染木の占有面積が枯損木及び健全木の各占有面積よりも大きい場合には、当該樹冠の被害区分は「感染木」とみなすということである。このように、1つの樹冠が1つの被害区分として区分されることにより、広域かつ詳細に、市町村ごとに1本単位で感染木の早期把握と、守るべき松林と防除エリアの選定を行うことができる。   Thereby, one crown can be classified as one damage category. For example, in a certain crown, as a result of calculating each occupation area of the infected tree, dead tree and healthy tree, in the crown, if the occupation area of the infected tree is larger than each occupation area of the dead tree and healthy tree, The damage classification of the tree crown is to be regarded as an “infected tree”. In this way, one tree crown is classified as one damage category, so that it is possible to quickly grasp the infected tree and select the pine forest and the pest control area to be protected for each municipality in a wide area and in detail. it can.

[8]本発明の松くい虫の被害区分算定方法においては、前記樹頂点別被害区分画像データ作成ステップは、前記地理的情報に基づいて、前記アカマツ林樹頂点画像データを前記樹冠別被害区分画像データに重ねることによって、樹頂点ごとに、前記被害区分を特定することが好ましい。   [8] In the method for calculating a damage category of a pine beetle according to the present invention, the step of creating the damage category image data for each vertex of the tree uses the top image data of the red pine forest tree as the damage category for each crown based on the geographical information. It is preferable to specify the damage category for each tree vertex by superimposing the image data.

このように、各樹頂点の被害区分を特定することで、個々の樹木ごとに被害区分を特定することができる。また、調査対象地域を上空から撮影して得られた撮影画像に被害区分を色分けした個々の樹木位置の分布画像を作成できるため、森林被害の様子をビジュアルに把握できる。なお、被害区分は感染木、枯損木及び健全木を指している。   Thus, by specifying the damage classification of each tree vertex, it is possible to specify the damage classification for each individual tree. In addition, it is possible to create a distribution image of individual tree positions by color-coding the damage classification on the photographed image obtained by photographing the survey area from above, so it is possible to visually grasp the state of forest damage. The damage category refers to infected trees, dead trees, and healthy trees.

[9]本発明の松くい虫の被害区分算定方法においては、前記被害区分本数算定ステップは、前記樹頂点別被害区分画像データに基づいて、前記アカマツ林の全領域又は任意の領域において前記地理的情報により小班と定義される森林所有者の境界データから小班区画を抽出し、抽出された小班において前記被害区分ごとに樹木の本数を算定することが好ましい。   [9] In the method for calculating the damage category of a pine beetle according to the present invention, the step of calculating the number of damage categories is based on the damage category image data classified by tree vertices in the entire region or an arbitrary region of the red pine forest. It is preferable to extract a subcompartment section from the boundary data of the forest owner defined as a subcompartment based on the target information, and to calculate the number of trees for each of the damage categories in the extracted subcompartment.

これにより、各小班区画において、被害区分すなわち感染木、枯損木及び健全木ごとに樹木の本数を算定することができる。   As a result, the number of trees can be calculated for each subdivision, for each damage category, that is, for each infected tree, dead tree, and healthy tree.

[10]本発明の松くい虫の被害区分算定装置は、調査対象となるアカマツ林を含む地域を上空から撮影して得られた撮影画像データに含まれる樹種の分光反射特性に基づいて作成されたアカマツ林画像データを用いて、前記アカマツ林に含まれるアカマツの本数を被害区分別に算定する松くい虫の被害区分算定装置であって、前記アカマツ林画像データに基づいて、前記アカマツの各樹頂点を抽出し、前記各樹頂点に関するアカマツ林樹頂点画像データを作成するアカマツ林樹頂点画像データ作成部と、前記アカマツ林画像データに基づいて、前記アカマツの各樹冠を抽出し、前記各樹冠に関するアカマツ林樹冠画像データを作成するアカマツ林樹冠画像データ作成部と、前記アカマツ林樹冠画像データと、松くい虫の被害を受けているアカマツから採取された被害葉及び健全なアカマツから採取された健全葉の各分光反射特性とに基づいて、前記各樹冠の領域における感染木、枯損木及び健全木の分光反射特性を作成し、作成した感染木、枯損木及び健全木の分光反射特性に基づいて樹冠別の被害区分を特定して、樹冠別被害区分画像データを作成する樹冠別被害区分画像データ作成部と、前記樹冠別被害区分画像データと前記アカマツ林樹頂点画像データとに基づいて、樹頂点別に前記被害区分を特定して、樹頂点別被害区分画像データを作成する樹頂点別被害区分画像データ作成部と、前記樹頂点別被害区分画像データと地理的情報とに基づいて、前記アカマツ林の全領域又は任意の領域に含まれるアカマツの本数を前記被害区分ごとに算定する被害区分本数算定部と、前記被害区分本数算定結果に基づいて、当該アカマツ林の全領域又は任意の領域に含まれる感染木と枯損木とを合計した被害木本数を、当該アカマツ林の全領域又は任意の領域に含まれるアカマツの全体本数で除することで得られる被害率を計算して被害率区分図を作成する被害率区分図作成部と、
を有すること特徴とする。
[10] The pine beetle damage classification calculating apparatus according to the present invention is created based on spectral reflection characteristics of tree species included in photographed image data obtained by photographing an area including a red pine forest to be investigated from the sky. An apparatus for calculating the damage category of a pine weevil that uses the red pine forest image data to calculate the number of red pine trees contained in the red pine forest according to the damage category, wherein each tree of the red pine tree is based on the red pine forest image data. A pine forest tree vertex image data creation unit that extracts vertices and creates red pine forest tree vertex image data relating to each tree vertex, and extracts each of the red pine tree crowns based on the red pine forest image data, Japanese pine forest canopy image data creation unit for creating red pine forest canopy image data, and the red pine forest canopy image data Based on the spectral reflection characteristics of the damaged leaves collected from the tree and the healthy leaves collected from the healthy red pine, the spectral reflection characteristics of the infected tree, the dead tree and the healthy tree in the area of each crown are created and created. Damage classification image data creation unit for identifying damage classification by crown based on the spectral reflection characteristics of infected infected trees, dead trees and healthy trees, and creating damage classification image data by crown, and said damage classification by crown Based on the image data and the red pine forest tree vertex image data, the damage classification image data creating unit for each tree vertex identifying the damage category for each tree vertex and creating the damage category image data for each tree vertex, and the tree vertex A damage category number calculating unit for calculating the number of red pine trees included in the entire area of the red pine forest or an arbitrary area based on the different damage category image data and geographical information; Based on the result of calculating the number of damage categories, the total number of damaged trees in the entire pine forest or in any area of the red pine forest is calculated as the total number of damaged trees in the entire pine forest or in any area. A damage rate classification diagram creation unit that calculates the damage rate obtained by dividing by the total number of
It is characterized by having.

本発明によれば、上記[1]に記載の松くい虫の被害区分算定方法と同様の効果が得られる。なお、本発明の松くい虫の被害区分算定装置は、当該松くい虫の被害区分算定装置に含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   According to the present invention, the same effect as the method for calculating the damage category of pine worms described in [1] above can be obtained. Note that the pine worm damage category calculating apparatus of the present invention has the functions of each of the above-described components included in the pine worm damage category calculating device installed as a computer program. By giving this data, the function of each component is executed on the software of the computer.

[11]本発明の松くい虫の被害区分算定装置においては、前記樹種の分光反射特性に基づいて前記撮影画像データからアカマツ林を抽出することによってアカマツ林画像データを作成するアカマツ林画像データ作成部をさらに有し、当該アカマツ林画像データ作成部によってアカマツ林画像データの作成を実施した後に、前記アカマツ林ごとに、前記アカマツ林樹頂点画像データ作成部による前記アカマツ林樹頂点画像データの作成と、前記アカマツ林樹冠画像データ作成部による前記アカマツ林樹冠画像データの作成と、前記樹冠別被害区分画像データ作成部による前記樹冠別被害区分画像データの作成と、前記樹頂点別被害区分画像データ作成部による前記樹頂点別被害区分画像データ作成と、前記被害区分本数算定部による前記被害区分本数算定と、前記被害率区分図作成部による前記被害率区分図の作成とを実施することが好ましい。   [11] In the pine beetle damage classification calculating apparatus according to the present invention, red pine forest image data is created by extracting red pine forest from the photographed image data based on the spectral reflection characteristics of the tree species. And generating the red pine forest tree vertex image data by the red pine forest tree vertex image data creating unit for each of the red pine forests after the red pine forest image data is created by the red pine forest image data creation unit. Creation of the red pine forest canopy image data by the red pine forest canopy image data creation unit, creation of the damage classification image data by the crown by the crown specific damage classification image data creation unit, and damage classification image data by the tree apex Creation of the damage classification image data for each tree vertex by the creation section and the damage section by the damage classification number calculation section And the number calculation, it is preferred to carry out the creation of the damage rate levels view by the damage rate levels diagram creating unit.

これにより、上記[2]に記載の松くい虫の被害区分算定方法と同様の効果が得られる。なお、本発明の松くい虫の被害区分算定装置は、当該松くい虫の被害区分算定装置に含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   Thereby, the effect similar to the damage classification calculation method of the pine weevil described in the above [2] is obtained. Note that the pine worm damage category calculating apparatus of the present invention has the functions of each of the above-described components included in the pine worm damage category calculating device installed as a computer program. By giving this data, the function of each component is executed on the software of the computer.

[12]本発明の松くい虫の被害区分算定装置においては、前記撮影画像データは、秋から冬に撮影された画像データを用いることが好ましい。   [12] In the pine weevil damage category calculating apparatus according to the present invention, it is preferable that the photographed image data is image data photographed from autumn to winter.

これにより、上記[3]に記載の松くい虫の被害区分算定方法と同様の効果が得られる。なお、本発明の松くい虫の被害区分算定装置は、当該松くい虫の被害区分算定装置に含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   Thereby, the same effect as the method for calculating the damage category of pine worms described in [3] above can be obtained. Note that the pine worm damage category calculating apparatus of the present invention has the functions of each of the above-described components included in the pine worm damage category calculating device installed as a computer program. By giving this data, the function of each component is executed on the software of the computer.

[13]本発明の松くい虫の被害区分算定装置においては、前記感染木、枯損木及び健全木のうちの前記感染木と前記枯損木との区分は、前記松くい虫の被害を受けているアカマツから採取された被害葉と健全なアカマツから採取された健全葉の各分光反射特性に基づいて行い、前記感染木は、可視光の赤色域において前記枯損木及び前記健全木に比べて反射率が高いことで、前記枯損木及び前記健全木と区分し、前記枯損木は、近赤外域において前記感染木及び前記健全木に比べて反射率が低いことで、前記感染木及び前記健全木と区分する、ことが好ましい。   [13] In the pine weevil damage classification calculating apparatus of the present invention, the classification of the infected tree and the dead tree among the infected tree, the dead tree, and the healthy tree is affected by the pine worm. The infected tree is reflected in the red region of visible light compared to the dead tree and the healthy tree. It is distinguished from the dead tree and the healthy tree by a high rate, and the dead tree has a lower reflectance than the infected tree and the healthy tree in the near infrared region, and thus the infected tree and the healthy tree. It is preferable to categorize.

これにより、上記[4]に記載の松くい虫の被害区分算定方法と同様の効果が得られる。なお、本発明の松くい虫の被害区分算定装置は、当該松くい虫の被害区分算定装置に含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   Thereby, the same effect as the method for calculating the damage category of pine worms described in [4] above can be obtained. Note that the pine worm damage category calculating apparatus of the present invention has the functions of each of the above-described components included in the pine worm damage category calculating device installed as a computer program. By giving this data, the function of each component is executed on the software of the computer.

[14]本発明の松くい虫の被害区分算定装置においては、前記アカマツ林樹頂点画像データ作成部は、輝度値が極大値を持つ点を自動抽出する局所最大値フィルターを用いて樹頂点を抽出することが好ましい。   [14] In the pine beetle damage classification calculating apparatus according to the present invention, the red pine forest tree vertex image data creating unit uses a local maximum value filter that automatically extracts a point having a maximum luminance value to determine a tree vertex. It is preferable to extract.

これにより、上記[5]に記載の松くい虫の被害区分算定方法と同様の効果が得られる。なお、本発明の松くい虫の被害区分算定装置は、当該松くい虫の被害区分算定装置に含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   Thereby, the same effect as the method for calculating the damage category of pine worms described in [5] above can be obtained. Note that the pine worm damage category calculating apparatus of the present invention has the functions of each of the above-described components included in the pine worm damage category calculating device installed as a computer program. By giving this data, the function of each component is executed on the software of the computer.

[15]本発明の松くい虫の被害区分算定装置においては、前記アカマツ林樹冠画像データ作成部は、樹冠の縁部が樹冠部分よりも暗いことを利用して行うものであり、これは、ブァレイフォローイング・アルゴリズムと自動マスク処理で自動的に抽出することが好ましい。   [15] In the pine beetle damage classification calculating apparatus according to the present invention, the red pine forest canopy image data creation unit is performed by utilizing the fact that the edge of the canopy is darker than the canopy portion. It is preferable that the extraction is automatically performed by the valley following algorithm and the automatic mask process.

これにより、上記[6]に記載の松くい虫の被害区分算定方法と同様の効果が得られる。なお、本発明の松くい虫の被害区分算定装置は、当該松くい虫の被害区分算定装置に含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   Thereby, the same effect as the method for calculating the damage category of pine worms described in [6] above can be obtained. Note that the pine worm damage category calculating apparatus of the present invention has the functions of each of the above-described components included in the pine worm damage category calculating device installed as a computer program. By giving this data, the function of each component is executed on the software of the computer.

[16]本発明の松くい虫の被害区分算定装置においては、前記樹冠別被害区分画像データ作成ステップは、前記被害区分を感染木、枯損木及び健全木の3種類として、前記アカマツの各樹冠から感染木、枯損木及び健全木が抽出された場合、前記各樹冠から抽出された感染木、枯損木及び健全木のそれぞれの占有面積を算出し、算出された感染木、枯損木及び健全木の各占有面積のうち、占有面積が大きい被害区分を当該樹冠の被害区分とみなすことが好ましい。   [16] In the pine beetle damage classification calculating apparatus according to the present invention, the crown-specific damage classification image data creating step includes setting the damage classification as an infected tree, a dead tree and a healthy tree, and each crown of the red pine. When infected trees, dead trees, and healthy trees are extracted from the above, the occupied areas of the infected trees, dead trees, and healthy trees extracted from the respective crowns are calculated, and the calculated infected trees, dead trees, and healthy trees are calculated. Of these occupying areas, it is preferable to regard a damage category having a large occupancy area as a damage category of the crown.

これにより、上記[7]に記載の松くい虫の被害区分算定方法と同様の効果が得られる。なお、本発明の松くい虫の被害区分算定装置は、当該松くい虫の被害区分算定装置に含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   Thereby, the same effect as the method for calculating the damage category of pine worms described in [7] above can be obtained. Note that the pine worm damage category calculating apparatus of the present invention has the functions of each of the above-described components included in the pine worm damage category calculating device installed as a computer program. By giving this data, the function of each component is executed on the software of the computer.

[17]本発明の松くい虫の被害区分算定装置においては、前記樹頂点別被害区分画像データ作成ステップは、前記地理的情報に基づいて、前記アカマツ林樹頂点画像データを前記樹冠別被害区分画像データに重ねることによって、樹頂点ごとに、前記被害区分を特定することが好ましい。   [17] In the pine beetle damage category calculating apparatus according to the present invention, the tree-based damage category image data creating step includes converting the red pine forest tree vertex image data into the tree-based damage category based on the geographical information. It is preferable to specify the damage category for each tree vertex by superimposing the image data.

これにより、上記[8]に記載の松くい虫の被害区分算定方法と同様の効果が得られる。なお、本発明の松くい虫の被害区分算定装置は、当該松くい虫の被害区分算定装置に含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   Thereby, the same effect as the method for calculating the damage category of pine worms described in [8] above can be obtained. Note that the pine worm damage category calculating apparatus of the present invention has the functions of each of the above-described components included in the pine worm damage category calculating device installed as a computer program. By giving this data, the function of each component is executed on the software of the computer.

[18]本発明の松くい虫の被害区分算定装置においては、前記被害区分本数算定部は、前記樹頂点別被害区分画像データに基づいて、前記アカマツ林の全領域又は任意の領域において前記地理的情報により小班と定義される森林所有者の境界データから小班区画を抽出し、抽出された小班において前記被害区分ごとに樹木の本数を算定することが好ましい。   [18] In the pine beetle damage category calculating apparatus according to the present invention, the damage category number calculating unit is configured to perform the geography in the entire area of the red pine forest or an arbitrary area based on the damage category image data classified by tree vertices. It is preferable to extract a subcompartment section from the boundary data of the forest owner defined as a subcompartment based on the target information, and to calculate the number of trees for each of the damage categories in the extracted subcompartment.

これにより、上記[9]に記載の松くい虫の被害区分算定方法と同様の効果が得られる。なお、本発明の松くい虫の被害区分算定装置は、当該松くい虫の被害区分算定装置に含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   Thereby, the same effect as the method for calculating the damage category of pine worms described in [9] above can be obtained. Note that the pine worm damage category calculating apparatus of the present invention has the functions of each of the above-described components included in the pine worm damage category calculating device installed as a computer program. By giving this data, the function of each component is executed on the software of the computer.

実施形態に係る松くい虫の被害区分算定装置100の構成を示す図である。It is a figure which shows the structure of the damage classification calculation apparatus 100 of the pine weevil which concerns on embodiment. 実施形態に係る松くい虫の被害区分算定方法を説明するフローチャートである。It is a flowchart explaining the damage classification calculation method of the pine weevil which concerns on embodiment. 松くい虫被害林における感染木の集団と枯損木の集団を示す図である。It is a figure which shows the group of the infected tree and the group of dead trees in a pine stake insect damage forest. 分光反射計による健全葉、感染葉、枯損葉の測定を行う際に使用する機材及び測定対象葉の一例を示す図である。It is a figure which shows an example of the equipment used when measuring the healthy leaf, an infected leaf, and a dead leaf by a spectral reflectometer, and a measurement object leaf. 健全葉、感染葉、枯損葉の分光反射特性を示す図である。It is a figure which shows the spectral reflection characteristic of a healthy leaf, an infected leaf, and a dead leaf. ディスプレイ上に表示された撮影画像を示す図である。It is a figure which shows the picked-up image displayed on the display. ディスプレイ上に表示されたアカマツ林画像を示す図である。It is a figure which shows the pine forest image displayed on the display. 樹木の樹頂点と樹冠を説明する図である。It is a figure explaining the tree vertex and tree crown of a tree. ディスプレイ上に表示されたアカマツ林画像から抽出したアカマツ林樹頂点画像データを白点で示す図である。It is a figure which shows the pine forest tree vertex image data extracted from the pine forest image displayed on the display with a white point. ディスプレイ上に表示されたアカマツ林樹冠画像を示す図である。It is a figure which shows the pine forest canopy image displayed on the display. ディスプレイ上に表示された樹冠別被害区分画像を示す図である。It is a figure which shows the damage classification image classified by tree crown displayed on the display. ディスプレイ上に表示された樹頂点別被害区分画像を示す図である。It is a figure which shows the damage classification image classified by tree | vertice displayed on the display. ディスプレイ上に表示された松くい虫の被害区分算定画像を示す図である。It is a figure which shows the damage classification calculation image of the pine weevil displayed on the display. ディスプレイ上に表示された被害状況データベースの一例を示す図である。It is a figure which shows an example of the damage condition database displayed on the display. ディスプレイ上に表示された被害率区分図を示す図である。It is a figure which shows the damage rate division figure displayed on the display.

以下、本発明の実施形態を図面に基づき記述する。   Hereinafter, embodiments of the present invention will be described with reference to the drawings.

図1は、実施形態における松くい虫の被害区分算定装置100の構成を示す図である。実施形態に係る松くい虫の被害区分算定装置100は、図1に示すように、撮影画像データ入力部10と、アカマツ林画像データ作成部20と、アカマツ林樹頂点画像データ作成部31と、アカマツ林樹冠画像データ作成部32と、樹冠別被害区分画像データ作成部33、樹頂点別被害区分画像データ作成部34と、被害区分本数算定部40と、被害率区分図作成部50とを有している。なお、アカマツ林樹頂点画像データ作成部31と、アカマツ林樹冠画像データ作成部32と、樹冠別被害区分画像データ作成部33と、樹頂点別被害区分画像データ作成部34とによって、アカマツ林樹頂点別被害区分システム30が構成される。   FIG. 1 is a diagram illustrating a configuration of a pine worm damage category calculation apparatus 100 according to the embodiment. As shown in FIG. 1, a pine beetle damage category calculating apparatus 100 according to the embodiment includes a photographed image data input unit 10, a red pine forest image data creation unit 20, a red pine forest tree vertex image data creation unit 31, There is a red pine forest canopy image data creation unit 32, a damage classification image data creation unit 33 by crown, a damage classification image data creation unit 34 by tree apex, a damage number calculation unit 40, and a damage rate classification diagram creation unit 50 doing. The red pine forest tree vertex image data creation unit 31, the red pine forest crown image data creation unit 32, the damage classification image data creation unit 33 by crown, and the damage classification image data creation unit 34 by tree vertex provide a red pine forest tree. The damage classification system 30 according to vertex is comprised.

撮影画像データ入力部10は、調査対象となるアカマツ林を含む地域(調査対象地域という。)を上空から撮影して得られた撮影画像データを入力する機能を有している。   The photographed image data input unit 10 has a function of inputting photographed image data obtained by photographing a region including a red pine forest to be surveyed (referred to as a survey target region) from above.

アカマツ林画像データ作成部20は、撮影画像データに含まれる樹種の分光反射特性に基づいて、撮影画像データからアカマツ林を抽出することによってアカマツ林画像データを作成する機能を有している。   The red pine forest image data creating unit 20 has a function of creating red pine forest image data by extracting red pine forest from the photographed image data based on the spectral reflection characteristics of the tree species included in the photographed image data.

アカマツ林樹頂点画像データ作成部31は、アカマツ林画像データに基づいて、アカマツの各樹頂点を抽出し、各樹頂点に関するアカマツ林樹頂点画像データを作成する機能を有している。具体的には、アカマツ林樹頂点画像データ作成部31は、輝度値が極大値を持つ点を自動抽出する局所最大値フィルターを用いてアカマツの樹頂点を抽出する。   The red pine forest tree vertex image data creation unit 31 has a function of extracting each tree vertices of red pine based on the red pine forest image data and creating red pine forest tree vertex image data regarding each tree vertex. Specifically, the red pine forest tree vertex image data creation unit 31 extracts the red pine tree vertices using a local maximum value filter that automatically extracts a point having a maximum luminance value.

アカマツ林樹冠画像データ作成部32は、アカマツ林画像データに基づいて、アカマツの各樹冠を抽出し、各樹冠に関するアカマツ林樹冠画像データを作成する機能を有している。具体的には、アカマツ林樹冠画像データ作成部32は、樹冠の縁部(樹冠縁)が樹冠部分よりも暗いことを利用して樹冠を抽出する。   The red pine forest canopy image data creation unit 32 has a function of extracting each red pine canopy based on the red pine forest image data and creating red pine forest canopy image data relating to each canopy. Specifically, the red pine forest canopy image data creation unit 32 extracts the crown using the fact that the edge (crown edge) of the crown is darker than the crown portion.

樹冠別被害区分画像データ作成部33は、アカマツ林樹冠画像データと、松くい虫の被害を受けているアカマツから採取された被害葉及び健全なアカマツから採取された健全葉の各分光反射特性とに基づいて、各樹冠の領域における感染木、枯損木及び健全木の分光反射特性を作成し、作成した感染木、枯損木及び健全木の分光反射特性に基づいて樹冠別の被害区分(感染木、枯損木及び健全木)を特定して、樹冠別被害区分画像データを作成する機能を有している。   The damage classification image data creation unit 33 according to canopy includes a red pine forest canopy image data, spectral reflection characteristics of damaged leaves collected from red pine that have been damaged by pine needles, and healthy leaves collected from healthy red pine. Based on the above, the spectral reflection characteristics of infected trees, dead trees and healthy trees in the area of each canopy are created, and damage classification (infected trees) by crown based on the spectral reflection characteristics of the infected trees, dead trees and healthy trees created. , Withered trees and healthy trees), and has a function of creating damage classification image data by crown.

具体的には、樹冠別被害区分画像データ作成部33は、松くい虫の被害区分を感染木、枯損木及び健全木の3種類として、アカマツの各樹冠から感染木、枯損木及び健全木が抽出された場合、各樹冠から抽出された感染木、枯損木及び健全木のそれぞれの占有面積を算出し、算出された感染木、枯損木及び健全木の各占有面積のうち、占有面積が大きい被害区分を当該樹冠の被害区分とみなすことによって、樹冠別被害区分画像データを作成する。   Specifically, the damage classification image data creation unit 33 classified by tree crown uses three types of damage categories of pine worms: infected tree, dead tree, and healthy tree. If extracted, calculate the occupied area of each infected tree, dead tree, and healthy tree extracted from each crown, and among the calculated occupied area of the infected tree, dead tree, and healthy tree, the occupied area is large. By considering the damage classification as the damage classification of the crown, the damage classification image data for each crown is created.

樹頂点別被害区分画像データ作成部34は、樹冠別被害区分画像データとアカマツ林樹頂点画像データとに基づいて、樹頂点別に被害区分(感染木、枯損木及び健全木)を特定して、樹頂点別被害区分画像データを作成する機能を有している。具体的には、樹頂点別被害区分画像データ作成部34は、地理的情報に基づいて、アカマツ林樹頂点画像データを樹冠別被害区分画像デ−タに重ねることによって、樹頂点ごとに被害区分を特定して樹頂点別被害区分画像データを作成する。なお、地理情報というのは、この場合、都道府県などで整備されている森林の管理台帳などに登録されている森林域の全領域又は任意の領域を指す。当該地理情報は、図1に示す地理情報システム(GISという。)60に登録されている。   The damage classification image data creation unit 34 classified by tree vertices identifies the damage classification (infected tree, dead tree and healthy tree) for each tree vertex based on the damage classification image data by tree crown and the red pine forest tree vertex image data, It has a function to create tree-based damage classification image data. Specifically, the tree-based damage classification image data creation unit 34 superimposes the red pine forest tree vertex image data on the tree-specific damage classification image data based on the geographical information, thereby determining the damage classification for each tree vertex. The damage classification image data classified by tree vertices is created. In this case, the geographical information refers to the entire area of the forest area or an arbitrary area registered in the forest management ledger or the like maintained in the prefecture. The geographic information is registered in the geographic information system (referred to as GIS) 60 shown in FIG.

被害区分本数算定部40は、樹頂点別被害区分画像データと地理的情報とに基づいて、アカマツ林の全領域又は任意の領域に含まれるアカマツの本数を被害区分(感染木、枯損木及び健全木)ごとに算定して被害区分本数の算定を行う機能を有している。具体的には、被害区分本数算定部は、松くい虫の被害区分を感染木、枯損木及び健全木の3種類として、各被害区分(感染木、枯損木及び健全木)ごとに、アカマツ林の全領域又は任意の領域において、小班と定義される森林所有者の境界データに基づいて小班区画を抽出し、抽出された小班に関する被害区分ごとに樹木の本数を算定する。   The damage classification number calculation unit 40 determines the number of red pine trees contained in the whole or arbitrary area of the red pine forest based on the damage classification image data for each vertex and geographical information. It has a function to calculate the number of damage categories by calculating each tree. Specifically, the damage category number calculation section uses three types of pine worm damage categories as infected trees, dead trees and healthy trees, and pine forests for each damage category (infected trees, dead trees and healthy trees). In all or any area of the above, subdivisions are extracted based on the boundary data of forest owners defined as subdivisions, and the number of trees is calculated for each damage category related to the extracted subdivisions.

被害率区分図作成部50は、被害区分本数算定結果に基づいて、アカマツ林の全領域又は任意の領域に含まれる感染木と枯損木とを合計した被害木本数を、当該アカマツ林の全領域又は任意の領域に含まれるアカマツの全体本数で除することで得られる被害率を計算して被害率区分図を作成する機能を有している。   Based on the result of calculating the number of damage categories, the damage rate category map creation unit 50 calculates the total number of damaged trees in the entire area of the red pine forest by adding the total number of infected trees and dead trees contained in any area of the red pine forest. Or, it has a function to create a damage rate classification chart by calculating the damage rate obtained by dividing by the total number of red pine contained in an arbitrary area.

なお、被害率区分図を作成する際においては、被害区分本数算定部40で算定された被害区分(感染木、枯損木及び健全木)のそれぞれの本数、感染木と枯損木とを合計した被害木本数、全体本数(合計)及び被害率などをデータベース化し、当該データベース化した内容に基づいて、作表したり、森林調査簿データベース(この場合、都道府県などで整備されている森林の管理台帳などのデータベース:例えば図14参照。)又は地理情報システム(G1S)に登録したりする。そして、当該作表した内容や、森林調査簿データベース又はG1Sに登録した内容に基づいて被害率区分図を作成する。   In preparing the damage rate classification diagram, the total number of damage categories (infected trees, dead trees and healthy trees) calculated by the damage category number calculation unit 40, and the total number of infected trees and dead trees Create a database of the number of trees, total number (total), damage rate, etc., based on the contents of the database, create a forest survey database (in this case, forest management ledgers etc. maintained in prefectures, etc.) (For example, see FIG. 14) or the geographic information system (G1S). Then, a damage rate classification diagram is created based on the tabled contents and the contents registered in the forest survey book database or G1S.

また、図1に示した松くい虫の被害区分算定装置100は、当該松くい虫の被害区分算定装置100に含まれる上記各構成要素が有する機能が、コンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   The pine worm damage category calculating apparatus 100 shown in FIG. 1 has the functions of the components included in the pine worm damage category calculating apparatus 100 installed as a computer program. By giving predetermined data to each component, the function of each component is executed on the software of the computer.

図2は、実施形態に係る松くい虫の被害区分算定方法における松くい虫の被害区分算定処理を説明するフローチャートである。松くい虫の被害区分算定処理の手順は、図2に示すように、調査対象となるアカマツ林を含む地域を上空から撮影して得られた撮影画像データに対応する撮影画像(カラー画像)の表示を行う(ステップS10)。また、撮影画像データに含まれる樹種の分光反射特性に基づいて抽出されたアカマツ林に関するアカマツ林画像データを作成する(ステップS20)。そして、アカマツ林画像データに基づいて、アカマツの各樹頂点を抽出し、各樹頂点に関するアカマツ林樹頂点画像データを作成する(ステップS31)とともに、アカマツ林画像データに基づいて、前記アカマツの各樹冠を抽出し、各樹冠に関するアカマツ林樹冠画像データを作成する(ステップS32)。   FIG. 2 is a flowchart for explaining a pine worm damage category calculation process in the pine worm damage category calculation method according to the embodiment. As shown in FIG. 2, the pine weevil damage classification calculation procedure is as follows. The photographed image (color image) corresponding to the photographed image data obtained by photographing the area including the red pine forest to be investigated from the sky is shown in FIG. Display is performed (step S10). Also, red pine forest image data relating to red pine forest extracted based on the spectral reflection characteristics of the tree species included in the photographed image data is created (step S20). Then, based on the red pine forest image data, each of the red pine tree vertices is extracted, and the red pine forest tree vertex image data relating to the respective tree vertices is created (step S31). Tree crowns are extracted, and pine forest crown image data relating to each tree crown is created (step S32).

その後、アカマツ林樹冠画像データと、松くい虫の被害を受けているアカマツから採取された被害葉及び健全なアカマツから採取された健全葉の各分光反射特性とに基づいて、各樹冠の領域における感染木、枯損木及び健全木の分光反射特性を作成し、樹冠別に被害区分(感染木、枯損木及び健全木)を特定して樹冠別被害区分画像データを作成する(ステップS33)。   After that, based on the image data of the red pine forest canopy and the spectral reflection characteristics of the damaged leaves collected from the red pine and the healthy leaves collected from the healthy red pine, Spectral reflection characteristics of the infected tree, dead tree, and healthy tree are created, the damage classification (infected tree, dead tree, and healthy tree) is specified for each tree crown, and the damage classification image data for each tree crown is created (step S33).

そして、樹冠別被害区分画像データとアカマツ林樹頂点画像データとに基づいて、樹頂点別に被害区分(感染木、枯損木及び健全木)を特定して樹頂点別被害区分画像データを作成する(ステップS34)。その後、樹頂点別被害区分画像データと地理的情報とに基づいて、アカマツ林の全領域又は任意の領域に含まれるアカマツの本数を被害区分(感染木、枯損木及び健全木)ごとに算定して被害区分本数算定データを作成する(ステップS40)。   Then, based on the damage classification image data by tree crown and the red pine forest tree vertex image data, the damage classification (infected tree, dead tree and healthy tree) is identified by tree vertex and the damage classification image data by tree vertex is created ( Step S34). After that, based on the damage classification image data by tree vertex and geographical information, the number of red pine contained in the whole area of the pine forest or any area is calculated for each damage classification (infected tree, dead tree and healthy tree). The damage category number calculation data is created (step S40).

その後、被害区分本数算定データに基づいて、当該アカマツ林の全領域又は任意の領域に含まれる感染木と枯損木とを合計した被害木本数を、当該アカマツ林の全領域又は任意の領域に含まれるアカマツの全体本数で除することで得られる被害率を計算して被害率区分図を作成する(ステップS50)。   After that, based on the data for calculating the number of damage categories, the total number of damaged trees and infected trees in all areas of the pine forest or any area is included in all areas or any area of the pine forest. The damage rate obtained by dividing by the total number of red pine is calculated and a damage rate classification diagram is created (step S50).

図3は、松くい虫被害林における感染木の集団と枯損木の集団とを示す図である。図3において、領域1は感染木の集団を示しており、領域2は枯損木の集団を示している。なお、図3はモノクロ画像であるため、図3からは、感染木の集団(領域1)及び枯損木の集団(領域2)は色別として読み取ることは困難であるが、実際には、図3の元となるカラー画像上では、感染木の集団(領域1)は、くすんだ緑から茶色に変わりつつある色で示され、枯損木の集団(領域2)は全体が茶色で示されているため、ディスプレイ上において表示される撮影画像(カラー画像)上では、感染木の集団(領域1)及び枯損木の集団(領域2)を容易に読み取ることができる。   FIG. 3 is a diagram showing a group of infected trees and a group of dead trees in a pine scabies damaged forest. In FIG. 3, a region 1 indicates a group of infected trees, and a region 2 indicates a group of dead trees. Since FIG. 3 is a monochrome image, from FIG. 3, it is difficult to read the group of infected trees (region 1) and the group of dead trees (region 2) by color. On the original color image of 3, the infected tree population (region 1) is shown in a color that is changing from dull green to brown, and the dead tree population (region 2) is shown in brown overall Therefore, the group of infected trees (area 1) and the group of dead trees (area 2) can be easily read on the captured image (color image) displayed on the display.

図4は、分光反射計による健全葉、感染葉、枯損葉の測定を行う際に使用する機材及び測定対象葉の一例を示す図である。測定を行う際に使用する機材としては、図4に示すように、分光反射計3、標準白色版4を用いる。また、図4において、測定対象葉としては健全葉5が示されているが、感染葉及び枯損葉(これらの図示は省略する。)も測定対処葉として用いる。分光反射特性は太陽から放射される電磁波の可視域から近赤外域までの波長ごとの反射の強さを表すもので、標準白色版4の反射の強さを100として、その比を分光反射率とする。なお、健全葉とは健全木から採取された葉であり、感染葉は感染木から採取された葉であり、枯損葉とは採取された葉である。なお、図4はモノクロ画像であるため、図4からは、アカマツの健全葉が緑色であることを読み取ることは困難であるが、実際には、健全葉は濃い緑色となっている。   FIG. 4 is a diagram illustrating an example of equipment and measurement target leaves used when measuring healthy leaves, infected leaves, and dead leaves with a spectroreflectometer. As equipment used for the measurement, a spectral reflectometer 3 and a standard white plate 4 are used as shown in FIG. Further, in FIG. 4, the healthy leaf 5 is shown as the measurement target leaf, but the infected leaf and the dead leaf (these illustrations are omitted) are also used as the measurement countermeasure leaf. Spectral reflection characteristics represent the intensity of reflection of the electromagnetic wave radiated from the sun for each wavelength from the visible range to the near infrared range. The reflection intensity of the standard white plate 4 is 100, and the ratio is the spectral reflectance. And In addition, healthy leaves are leaves collected from healthy trees, infected leaves are leaves collected from infected trees, and dead leaves are collected leaves. Since FIG. 4 is a monochrome image, it is difficult to read from FIG. 4 that the healthy leaves of red pine are green, but in reality, the healthy leaves are dark green.

図5は、松くい虫の被害葉と健全葉の分光反射特性を示す図である。なお、松くい虫の被害葉というのは、感染葉と枯損葉を指している。松くい虫の被害葉と健全葉の分光反射特性は、実在する地域のアカマツの被害林分から感染木(当年生被害木:葉が黄緑に変色し元気のない松)、枯損木(前年度被害木:葉が茶色になって枯れた松)、健全木の各樹冠から感染葉、枯損葉、健全葉を採取して、図4に示すように人工衛星や航空機センサと同じ観測波長を計測できる分光反射計3で計測した結果、図5に示すような分光反射特性が得られた。   FIG. 5 is a diagram showing the spectral reflection characteristics of damaged leaves and healthy leaves of a pine weevil. Note that the damaged leaves of pine worms refer to infected leaves and dead leaves. Spectral reflection characteristics of damaged and healthy leaves of pine weevil are from infected trees of red pine trees in the area where they are actually infected (current year damaged trees: pine trees with discolored leaves that turn yellow-green), dead trees (previous year) Damaged trees: pine trees with brown leaves), infected leaves, dead leaves, and healthy leaves from each crown of healthy trees, and measured the same observation wavelength as satellites and aircraft sensors as shown in Fig. 4. As a result of measurement with the spectroreflectometer 3, the spectral reflection characteristics as shown in FIG. 5 were obtained.

図5において、太い実線で示す曲線L1は健全葉の分光反射特性を示し、破線で示す曲線L2は枯損葉の分光反射特性を示し、白抜きの実践で示す曲線L3は感染葉の分光反射特性を示している。分光反射特性は太陽から放射される電磁波の可視域から近赤外域までの波長ごとの反射の強さを表すもので、標準白色版を100として、その比を反射率とする。   In FIG. 5, a curved line L1 indicated by a thick solid line indicates the spectral reflection characteristic of healthy leaves, a curved line L2 indicated by a broken line indicates the spectral reflection characteristics of dead leaves, and a curved line L3 indicated by a white outline indicates the spectral reflection characteristics of infected leaves. Is shown. The spectral reflection characteristic represents the intensity of reflection of each wavelength of the electromagnetic wave radiated from the sun from the visible range to the near infrared range. The standard white plate is defined as 100, and the ratio is defined as the reflectance.

図5に示すように、3種の葉(感染葉、枯損葉及び健全葉)の分光反射特性に違いがあることから、感染葉、枯損葉及び健全葉の区分ができることを発見した。すなわち、図5に示すように、枯損葉の分光反射特性(曲線L2)によれば、枯損葉は近赤外域(770nm〜830nm)7で反射率が極めて低いことで、感染葉及び健全葉と容易に区分できる。一方、感染葉の分光反射特性(曲線L3)によれば、感染葉は可視光の赤色域(670nm〜710nm)6でクロロフィルと水分の減少により光合成能力が劣り電磁波の吸収が減少することにより、反射率が健全葉及枯損葉より高くなる。このため、感染葉は枯損葉及び健全葉と区分できる。   As shown in FIG. 5, since there is a difference in the spectral reflection characteristics of the three types of leaves (infected leaves, dead leaves, and healthy leaves), it was found that infected leaves, dead leaves, and healthy leaves can be classified. That is, as shown in FIG. 5, according to the spectral reflection characteristics (curve L2) of the dead leaves, the dead leaves have a very low reflectance in the near-infrared region (770 nm to 830 nm) 7. Can be easily segmented. On the other hand, according to the spectral reflection characteristics of the infected leaves (curve L3), the infected leaves are inferior in photosynthetic ability due to a decrease in chlorophyll and water in the visible red region (670 nm to 710 nm) 6, thereby reducing the absorption of electromagnetic waves. Reflectivity is higher than healthy and dead leaves. For this reason, infected leaves can be distinguished from dead leaves and healthy leaves.

この結果は、図5に示すように、通常の人工衛星の観測波長(可視域の青・緑・赤、近赤外域)における4チャンネル(ch)の近赤外域で枯損葉を感染葉及び健全葉に対して区分でき、最新型の8チャンネル人工衛星では、図5に示すように、可視域の赤色域においては3チャンネルと多波長であることから、感染葉を枯損葉及健全葉に対して、高精度に区分できることを示している。   As shown in FIG. 5, this result shows that dead leaves are infected and healthy in the near-infrared region of 4 channels (ch) at the observation wavelengths of ordinary satellites (blue, green, red, near-infrared region). As shown in Fig. 5, the latest 8-channel satellite can be divided into 3 channels and multiple wavelengths in the red region of the visible region. Thus, it can be classified with high accuracy.

このように、3種の葉(感染葉、枯損葉及び健全葉)の分光反射特性に違いによって、感染葉、枯損葉及び健全葉の区分を行うことにより、感染木、枯損木及び健全木の区分を高精度に行うことができる。特に、感染木と枯損木との区分についても高精度に行うことができる。   Thus, by classifying infected leaves, dead leaves, and healthy leaves depending on the spectral reflection characteristics of the three types of leaves (infected leaves, dead leaves, and healthy leaves), infected trees, dead trees, and healthy trees Classification can be performed with high accuracy. In particular, it is possible to accurately classify infected trees and dead trees.

図6は、ディスプレイ上に表示された撮影画像を示す図である。図6に示す撮影画像は、調査対象となるアカマツ林を含む地域を上空から撮影して得られた撮影画像データに対応する画像である。当該撮影画像において、全体的に黒色で示されている領域は森林域であり、白色で示されている領域は、集落及び田畑などである。また、森林域においても、アカマツ林の領域とその他の樹木の領域とが存在する。   FIG. 6 is a diagram illustrating a captured image displayed on the display. The photographed image shown in FIG. 6 is an image corresponding to photographed image data obtained by photographing an area including a red pine forest to be investigated from the sky. In the photographed image, the area shown in black as a whole is a forest area, and the area shown in white is a village, a field, and the like. Also in the forest area, there are areas of red pine forest and other trees.

なお、図6はモノクロ画像であるため、図6からは、アカマツ林の領域と、その他の樹木の領域とは色別として読み取ることは困難であるが、実際には、図6の元となるカラー画像上では、アカマツ林は紫色で示され、その他の樹木による森林域は緑色で示されているため、ディスプレイ上において表示される撮影画像(カラー画像)上では、アカマツ林とその他の樹木とを容易に区別することができる。   Since FIG. 6 is a monochrome image, it is difficult to read the red pine forest area and the other tree areas separately by color from FIG. 6, but in practice it is the origin of FIG. On the color image, the red pine forest is shown in purple, and the forest area of other trees is shown in green. Therefore, on the photographed image (color image) displayed on the display, the red pine forest and other trees Can be easily distinguished.

図7は、ディスプレイ上に表示されたアカマツ林画像を示す図である。図7に示すアカマツ林画像は、アカマツ林画像データ作成部20によって作成されたアカマツ林画像データに対応する画像である。当該アカマツ林画像は、図6に示す撮影画像においてアカマツ林の領域を抽出するために、アカマツ林の領域を特定の色(橙色とする。)で表したものである。なお、図7はモノクロ画像であるため、図7からは、アカマツ林を橙色として読み取ることは困難であるが、実際には、図7の元となるカラー画像上では、アカマツ林は橙色で示されているため、ディスプレイ上において表示される撮影画像(カラー画像)上では、アカマツ林を容易に読み取ることができる。   FIG. 7 shows a red pine forest image displayed on the display. The red pine forest image shown in FIG. 7 is an image corresponding to the red pine forest image data created by the red pine forest image data creation unit 20. The red pine forest image is obtained by expressing the red pine forest area in a specific color (in orange) in order to extract the red pine forest area in the photographed image shown in FIG. Since FIG. 7 is a monochrome image, it is difficult to read red pine forest as orange from FIG. 7, but actually, the red pine forest is shown in orange on the original color image of FIG. Therefore, the red pine forest can be easily read on the photographed image (color image) displayed on the display.

図8は、樹木の樹頂点と樹冠を説明する図である。図8(a)は森林を地上の横(ほぼ水平方向)から見た場合を模式的に示す図であり、図8(b)は森林を上空からほぼ垂直方向に見た場合を模式的に示す図である。図8において、樹頂点を「Pe」で表し、樹冠を「Ti」で表している。複数の樹木が近接して存在するような場合、上空からみると、樹冠が3つの樹木の各樹冠が重なり合うため、図8(b)に示すように、1つの樹冠(図8(b)における実線で示す)ように撮影される。なお、個々の樹冠Tiにおいて当該樹冠Tiの輪郭線TiLを「樹冠の縁部(樹冠縁)」とし、輪郭線TiLの内部(輪郭線TiLで囲まれている部分)を「樹冠部分」とする。また、樹冠の縁部及び樹冠の内部を総称して「樹冠の領域」という。   FIG. 8 is a diagram illustrating tree vertices and tree crowns. FIG. 8A is a diagram schematically illustrating a case where the forest is viewed from the side (almost horizontal direction) on the ground, and FIG. 8B is a diagram schematically illustrating a case where the forest is viewed from the sky almost vertically. FIG. In FIG. 8, the tree apex is represented by “Pe”, and the tree crown is represented by “Ti”. When there are a plurality of trees close to each other, as seen from above, the crowns of the three trees overlap each other, so that as shown in FIG. 8B, one crown (in FIG. 8B) Taken as shown by the solid line). In each crown Ti, the contour line TiL of the crown Ti is defined as “the edge of the crown (crown edge)”, and the inside of the contour line TiL (the portion surrounded by the contour line TiL) is defined as the “crown portion”. . In addition, the edge of the crown and the inside of the crown are collectively referred to as a “crown area”.

図9は、ディスプレイ上に表示されたアカマツ林画像から抽出したアカマツ林樹頂点画像を示す図である。図9に示すアカマツ林樹頂点画像は、アカマツ林樹頂点画像データ作成部31によって作成されたアカマツ林樹頂点画像データに対応する画像である。当該アカマツ林樹頂点画像は、図7において説明したような処理(アカマツ林を抽出するための処理)を施した画像から、太陽光を強く反射し、アカマツ林樹冠内で最大値の各樹頂点を局所最大値フィルター処理で自動的に抽出することによって得ることができる。このような処理を施すことによって、肉眼判読では困難な樹頂点を、精度良く抽出することができる。なお、図9は図7に示す画像の所定の領域を拡大して示すものである。   FIG. 9 is a diagram showing a red pine forest tree vertex image extracted from the red pine forest image displayed on the display. The red pine forest tree vertex image shown in FIG. 9 is an image corresponding to the red pine forest tree vertex image data created by the red pine forest tree vertex image data creation unit 31. The peak image of the red pine forest tree reflects the sunlight strongly from the image that has been processed as described in FIG. 7 (process for extracting the red pine forest), and each peak of the maximum value in the crown of the red pine forest Can be obtained automatically by local maximum filtering. By performing such processing, tree vertices difficult to read with the naked eye can be extracted with high accuracy. FIG. 9 is an enlarged view of a predetermined area of the image shown in FIG.

図9はモノクロ画像であるため、図9においては、アカマツ林は黒色で示され、各アカマツの樹頂点は白い点で示されている。このため、図9からは、アカマツ林の領域及びその他の樹木による領域を色別として読み取ることは困難であるが、実際には、図9の元となるカラー画像上では、アカマツ林は特定の色(図9においては図7と異なり紫色)で示され、その他の樹木による森林域は緑色で示されているため、ディスプレイ上において表示される撮影画像(カラー画像)上では、アカマツ林とその他の樹木とを容易に区別することができる。   Since FIG. 9 is a monochrome image, the red pine forest is shown in black in FIG. 9, and the tree vertices of each red pine are shown as white dots. For this reason, from FIG. 9, it is difficult to read the pine forest area and other tree areas as different colors, but in reality, the pine forest is a specific area on the original color image of FIG. In FIG. 9, it is shown in purple (unlike in FIG. 7), and the forest area by other trees is shown in green. Therefore, on the photographed image (color image) displayed on the display, red pine forest and others Can be easily distinguished from other trees.

図10は、ディスプレイ上に表示されたアカマツ林樹冠画像を示す図である。図10に示すアカマツ林樹冠画像は、アカマツ林樹冠画像データ作成部32によって作成されたアカマツ林樹冠画像データに対応する画像である。図10に示すアカマツ林樹冠画像は、図9に示すアカマツ林樹頂点画像と同様に、アカマツ林は黒色で示されており、白抜きで囲まれた小領域が各々の樹冠を表している。   FIG. 10 is a diagram showing a red pine forest canopy image displayed on the display. The red pine forest canopy image shown in FIG. 10 is an image corresponding to the red pine forest canopy image data created by the red pine forest canopy image data creation unit 32. The red pine forest canopy image shown in FIG. 10 is similar to the red pine forest tree apex image shown in FIG. 9. The red pine forest is shown in black, and the small areas surrounded by white represent the respective crowns.

なお、図10はモノクロ画像であるため、図10からは、アカマツ林の領域とその他の樹木による領域は色別として読み取ることは困難であるが、実際には、図10の元となるカラー画像上では、アカマツ林は図9と同様に、紫色で示され、その他の樹木による森林域は緑色で示されているため、ディスプレイ上において表示される撮影画像(カラー画像)上では、アカマツ林とその他の樹木とを容易に区別することができる。   Since FIG. 10 is a monochrome image, it is difficult to read the red pine forest region and the other tree region from each other according to color from FIG. 10, but actually, the color image that is the original of FIG. Above, the red pine forest is shown in purple as in FIG. 9, and the forest area of the other trees is shown in green. Therefore, on the photographed image (color image) displayed on the display, It can be easily distinguished from other trees.

ところで、樹冠の抽出は、樹冠の縁部が樹冠部分よりも暗いことを利用して行うものであり、このように樹冠を抽出する技術は、公知のブァレイフォローイング・アルゴリズムと自動マスク処理によって可能となるものであり、例えば、下記の公知文献1が知られている。公知文献1は、ブァレイフォローイング・アルゴリズムを用いて、樹冠の抽出を求めるものである。   By the way, the extraction of the canopy is performed by utilizing the fact that the edge of the canopy is darker than the canopy portion, and the technique for extracting the canopy in this way is possible by using a well-known bray following algorithm and automatic mask processing For example, the following publicly known document 1 is known. In the known document 1, tree crown extraction is obtained using a valley following algorithm.

公知文献1:Masato Katoh、Francois Gougeon、Don Leckie “Application of high-resolution airborne data using individual tree crowns in Japanese conifer plantations“、Journal of Forestry Research 14(1)、2009年、P、10〜19   Known Document 1: Masato Katoh, Francois Gougeon, Don Leckie “Application of high-resolution airborne data using individual tree crowns in Japanese conifer plantations”, Journal of Forestry Research 14 (1), 2009, P, 10-19

図11は、ディスプレイ上に表示された樹冠別被害区分画像を示す図である。図11に示す被害区分画像は、樹冠別被害区分画像データ作成部33によって作成された樹冠別被害区分画像データに対応する画像である。   FIG. 11 is a diagram showing a damage classification image by tree crown displayed on the display. The damage classification image shown in FIG. 11 is an image corresponding to the damage classification image data for each crown created by the damage classification image data creation unit 33 for each crown.

なお、図11に示す樹冠別被害区分画像においては、感染木の集団が存在する領域8と枯損木の集団が存在する領域9とがそれぞれ複数箇所ずつ存在するが、図11はモノクロ画像であるため、図11からは、感染木の集団が存在する領域8、枯損木の集団が存在する領域9を色別として読み取ることは困難である。しかし、実際には、図11の元となるカラー画像上では、感染木の集団が存在する領域8と枯損木の集団が存在する領域9とを色別として読み取ることは容易である。この場合、図11の元となるカラー画像上では、感染木の集団が存在する領域8は明るい茶色で示され、枯損木の集団が存在する領域9は濃い茶色で示されているため、ディスプレイ上において表示される撮影画像(カラー画像)上では、感染木の集団が存在する領域8と枯損木の集団が存在する領域9とを容易に読み取ることができる。   In the damage classification image by crown shown in FIG. 11, there are a plurality of regions 8 where the infected tree group exists and a region 9 where the dead tree group exists, but FIG. 11 is a monochrome image. Therefore, from FIG. 11, it is difficult to read the area 8 where the infected tree population exists and the area 9 where the dead tree population exist as different colors. However, in practice, in the original color image in FIG. 11, it is easy to read the region 8 where the infected tree population exists and the region 9 where the dead tree population exist as different colors. In this case, in the original color image of FIG. 11, the area 8 where the infected tree group exists is shown in light brown, and the area 9 where the dead tree group exists is shown in dark brown. On the photographed image (color image) displayed above, it is possible to easily read the region 8 where the infected tree group exists and the region 9 where the dead tree group exists.

図12は、ディスプレイ上に表示された樹頂点別被害区分画像を示す図である。図12に示す樹頂点別被害区分画像は、樹頂点別被害区分画像データ作成部34によって作成された樹頂点別被害区分画像データに対応する画像である。   FIG. 12 is a diagram showing a damage classification image classified by tree vertices displayed on the display. The damage classification image classified by tree vertices shown in FIG. 12 is an image corresponding to the damage classification image data classified by tree vertices created by the damage classification image data creation part 34 classified by tree vertices.

図12において、各樹木(アカマツ)の樹頂点は、個々のアカマツごとに位置情報のX、Y座標を有している。このため、GPS(Global Positioning System)によって感染木、枯損木及び健全木の現地確認ができる。   In FIG. 12, the tree vertices of each tree (Pinus densiflora) have X and Y coordinates of position information for each Pinus densiflora. For this reason, on-site confirmation of infected trees, dead trees and healthy trees can be performed by GPS (Global Positioning System).

なお、図12はモノクロ画像であるため、感染木は樹頂点が灰色で示され、枯損木は樹頂点が黒色で示され、健全木は樹頂点が白色で示されている。このため、図12からは感染木、枯損木及び健全木を明確に区別することは困難である。しかし、図12の元となるカラー画像上では、感染木、枯損木及び健全木を読み取ることは容易である。この場合、図12の元となるカラー画像上では、感染木は樹頂点が明るい茶色で示され、枯損木は樹頂点が赤色で示され、健全木は樹頂点が白色で示されているため、ディスプレイ上において表示される撮影画像(カラー画像)上では、感染木、枯損木及び健全木を容易に読み取ることができる。   Since FIG. 12 is a monochrome image, an infected tree has a tree vertex shown in gray, a dead tree has a tree vertex shown in black, and a healthy tree has a tree vertex shown in white. For this reason, it is difficult to clearly distinguish an infected tree, a dead tree, and a healthy tree from FIG. However, it is easy to read an infected tree, a dead tree, and a healthy tree on the original color image of FIG. In this case, in the original color image of FIG. 12, the infected tree is shown with a light brown tree vertex, the dead tree is shown with a red tree vertex, and a healthy tree is shown with a white tree vertex. On the captured image (color image) displayed on the display, the infected tree, the dead tree, and the healthy tree can be easily read.

図13は、ディスプレイ上に表示された松くい虫の被害区分本数算定画像を示す図である。なお、図13は図12に示されている森林域よりも広い範囲の森林域が示されている。被害区分本数算定部40は、図13に示す被害区分本数算定画像に基づいて被害区分本数を算定する。   FIG. 13 is a diagram showing an image of calculating the number of damage categories of pine worms displayed on the display. Note that FIG. 13 shows a forest area in a wider range than the forest area shown in FIG. The damage category number calculation unit 40 calculates the number of damage categories based on the damage category number calculation image shown in FIG.

被害区分本数算定部40は、上記したように、樹頂点別被害区分画像データと地理的情報とに基づいて、アカマツ林の全領域又は任意の領域に含まれるアカマツの本数を被害区分(感染木、枯損木及び健全木)ごとに算定する。図13においては、小班と定義される森林所有者の境界データに基づいて小班区画(図13において白抜きの線で囲まれた区画)を抽出し、抽出された小班区画に関する被害区分ごとに樹木の本数を算定する。   As described above, the damage category number calculation unit 40 determines the number of red pine trees contained in the whole area of the pine forest or an arbitrary area based on the damage category image data classified by tree vertices and geographical information. , Dead and healthy trees). In FIG. 13, a subcompartment section (a section surrounded by a white line in FIG. 13) is extracted based on the boundary data of the forest owner defined as a subcompartment, and a tree is formed for each damage category related to the extracted subcompartment section. Calculate the number of

なお、図13はモノクロ画像であるため、図13からは、感染木、枯損木及び健全木を色別として読み取ることは困難であるが、実際には、図13の元となるカラー画像上では、感染木、枯損木及び健全木を読み取ることは容易である。この場合、図13の元となるカラー画像上では、感染木は明るい茶色で示され、枯損木は赤色で示され、健全木は白色で示されているため、ディスプレイ上において表示される撮影画像(カラー画像)上では、感染木、枯損木及び健全木を容易に読み取ることができる。   Since FIG. 13 is a monochrome image, it is difficult to read an infected tree, a dead tree, and a healthy tree as different colors from FIG. 13, but actually, on the color image that is the original of FIG. It is easy to read infected trees, dead trees and healthy trees. In this case, in the original color image of FIG. 13, the infected tree is shown in light brown, the dead tree is shown in red, and the healthy tree is shown in white. On (color image), it is possible to easily read infected trees, dead trees, and healthy trees.

図14は、ディスプレイ上に表示された被害状況データベースの一例を示す図である。例えば、被害区分本数算定部40で算定された感染木、枯損木及び健全木のそれぞれの本数、感染木と枯損木とを合計した被害木本数、合計(全体本数)、被害率などを、各林班の各小班に対応付けて記録する。なお、被害率は、各小班における被感染木と枯損木の合計した被害木本数を全体本数(合計)で除することで得られる。   FIG. 14 is a diagram illustrating an example of a damage status database displayed on the display. For example, the number of infected trees, dead trees and healthy trees calculated by the damage classification number calculation unit 40, the total number of damaged trees and the total number of damaged trees, the total number of damaged trees, the damage rate, etc. Record in association with each small group in the forest group. The damage rate can be obtained by dividing the total number of damaged trees by the total number (total) of infected trees and dead trees in each subdivision.

そして、図14に示すような被害状況データベースの内容に基づいて、作表したり、森林調査簿データベース(この場合、都道府県などで整備されている森林の管理台帳などのデータベース)又は地理情報システム(G1S)に登録したりする。これにより、アカマツ林の全領域又は任意の領域である小班における被害区分別の樹木本数を把握できるため、被害状況の診断と報告に有効である。   Then, based on the contents of the damage situation database as shown in FIG. 14, it can be created, forest survey database (in this case, a database such as a forest management ledger maintained in prefectures) or a geographic information system ( G1S). This makes it possible to grasp the number of trees by damage category in the entire area of the pine forest or an arbitrary area, which is effective in diagnosing and reporting the damage situation.

図15は、ディスプレイ上に表示された被害率区分図を示す図である。当該被害率区分図は、松くい虫による森林被害対策を進めている都道府県や市町村で求められている出力図である。当該被害率区分図は、アカマツ林における松くい虫の被害状況を視覚的に容易に理解できる資料として用いることができる。なお、図15において、被害率を示す凡例は被害率を5〜10パーセントごとに色分けし、どこの森林域や小班で被害があるのか一目で把握でき、松くい虫被害の対策に生かすことができる。   FIG. 15 is a diagram showing a damage rate classification diagram displayed on the display. The damage rate classification diagram is an output diagram required in prefectures and municipalities that are promoting forest damage countermeasures due to pine weevil. The damage rate classification chart can be used as a material for easily understanding the state of damage of pine weevil in a red pine forest. In addition, in Fig. 15, the legend indicating the damage rate is color-coded every 5 to 10%, and it is possible to grasp at a glance which forest area and sub-compartment is damaged, and it can be used for measures against pine weevil damage. it can.

なお、図15はモノクロ画像であるため、図15からは、被害率を色別として読み取ることは困難であるが、実際には、図15の元となるカラー画像上では、被害率を色別として読み取ることは容易である。一例を挙げると、濃い緑色で示されている領域は被害率が0%、黄緑色で示されている領域は被害率が6%〜10%であり、濃い赤色で示されている範囲は被害率が31%〜62%である。   Since FIG. 15 is a monochrome image, it is difficult to read the damage rate for each color from FIG. 15, but actually, the damage rate is displayed for each color on the original color image of FIG. It is easy to read as. For example, the area shown in dark green has a damage rate of 0%, the area shown in yellow green has a damage rate of 6% to 10%, and the area shown in dark red is damaged. The rate is 31% to 62%.

このような被害率区分図をアカマツ林における松くい虫の被害状況を示す資料として用いることにより、森林管理者や技術者などが市町村単位で松くい虫被害の詳細な被害分布をつかむことができ、被害の拡大を防ぐため最前線の感染木や枯損木の集団をとらえ、的確に防除する上での計画立案や森林管理に有効となる。従って、従来のように、調査員が森林を歩いて被害の様子を実際に見て回ることによって、松くい虫の被害状況などを把握する必要がなくなるため、調査に必要な人件費を削減することができるとともに調査に必要な時間も短縮することができ、松くい虫の被害区分の算定を高精度にかつ安価な費用で短時間に実現可能となる。   By using this damage rate classification chart as a material to show the damage situation of pine weevil in Pinus densiflora forest, forest managers and engineers can grasp the detailed distribution of pine weevil damage by municipality. In order to prevent the spread of damage, it will be effective for planning and forest management to capture the group of infected and dead trees at the forefront and to control them accurately. Therefore, it is no longer necessary for the investigator to walk through the forest and actually look at the damage as in the past, so it is no longer necessary to grasp the damage situation of the pine weevil, thus reducing the labor cost required for the investigation. In addition, the time required for the survey can be shortened, and the damage classification of pine worms can be calculated with high accuracy and at a low cost in a short time.

1・・・感染木の集団が存在する領域、2・・・枯損木の集団が存在する領域、3・・・分校放射計、4・・・標準白色板、5・・・測定対象となる葉(健全葉)、10・・・撮影画像データ入力部、20・・・アカマツ林画像データ作成部、30・・・アカマツ林樹頂点別被害区分システム、31・・・アカマツ林樹頂点画像データ作成部、32・・・アカマツ林樹冠画像データ作成部、33・・・樹冠別被害区分画像データ作成部、34・・・樹頂点別被害区分画像データ作成部、40・・・被害区分本数算定部、50・・・被害率区分図作成部、60・・・地理情報システム(GIS)、Pe・・・樹頂点、Ti・・・樹冠   DESCRIPTION OF SYMBOLS 1 ... Area | region where the group of infected trees exists, 2 ... Area | region where the group of dead trees exist, 3 ... Branch radiometer, 4 ... Standard white board, 5 ... It becomes a measuring object Leaf (healthy leaf), 10 ... photographed image data input unit, 20 ... red pine forest image data creation unit, 30 ... damage classification system by red pine forest tree vertices, 31 ... red pine forest tree vertex image data Creation part, 32 ... Red pine forest canopy image data creation part, 33 ... Damage classification image data creation part by crown, 34 ... Damage classification image data creation part by tree apex, 40 ... Damage classification number calculation Part, 50 ... damage rate classification diagram creation part, 60 ... Geographic Information System (GIS), Pe ... tree vertex, Ti ... tree crown

Claims (18)

調査対象となるアカマツ林を含む地域を上空から撮影して得られた撮影画像データに含まれる樹種の分光反射特性に基づいて作成されたアカマツ林画像データを用いて、前記アカマツ林に含まれるアカマツの本数を被害区分別に算定する松くい虫の被害区分算定方法であって、
前記アカマツ林画像データに基づいて、前記アカマツの各樹頂点を抽出し、前記各樹頂点に関するアカマツ林樹頂点画像データを作成するアカマツ林樹頂点画像データ作成ステップと、
前記アカマツ林画像データに基づいて、前記アカマツの各樹冠を抽出し、前記各樹冠に関するアカマツ林樹冠画像データを作成するアカマツ林樹冠画像データ作成ステップと、
前記アカマツ林樹冠画像データと、松くい虫の被害を受けているアカマツから採取された被害葉及び健全なアカマツから採取された健全葉の各分光反射特性とに基づいて、前記各樹冠の領域における感染木、枯損木及び健全木の分光反射特性を作成し、作成した感染木、枯損木及び健全木の分光反射特性に基づいて樹冠別の被害区分を特定して、樹冠別被害区分画像データを作成する樹冠別被害区分画像データ作成ステップと、
前記樹冠別被害区分画像データと前記アカマツ林樹頂点画像データとに基づいて、樹頂点別に前記被害区分を特定して、樹頂点別被害区分画像データを作成する樹頂点別被害区分画像データ作成ステップと、
前記樹頂点別被害区分画像データと地理的情報とに基づいて、前記アカマツ林の全領域又は任意の領域に含まれるアカマツの本数を前記被害区分ごとに算定する被害区分本数算定ステップと、
前記被害区分本数算定結果に基づいて、当該アカマツ林の全領域又は任意の領域に含まれる感染木と枯損木とを合計した被害木本数を、当該アカマツ林の全領域又は任意の領域に含まれるアカマツの全体本数で除することで得られる被害率を計算して被害率区分図を作成する被害率区分図作成ステップと、
を有すること特徴とする松くい虫の被害区分算定方法。
Using the red pine forest image data created based on the spectral reflection characteristics of the tree species included in the photographed image data obtained by photographing the region including the red pine forest to be surveyed from above, the red pine forest included in the red pine forest is used. Is a method for calculating the damage category of pine beetles,
Based on the red pine forest image data, each of the red pine tree vertices is extracted, and a red pine forest tree vertex image data creating step for creating red pine forest tree vertex image data relating to each tree vertex;
Based on the red pine forest image data, each of the red pine canopies is extracted, and a red pine forest canopy image data creation step for creating a red pine forest canopy image data for each of the canopies,
Based on the image data of the red pine forest canopy and the spectral reflection characteristics of the damaged leaves collected from the red pine and the healthy leaves collected from the healthy red pine, Create spectral reflection characteristics of infected, dead, and healthy trees, identify damage classifications by tree crown based on the spectral reflection characteristics of the infected, dead, and healthy trees, and generate damage classification image data by crown. Step of creating damage classification image data by tree crown,
Based on the damage classification image data by tree crown and the pine forest tree vertex image data, the damage classification image data creation step by tree vertex which identifies the damage classification by tree vertex and creates damage classification image data by tree vertex When,
Based on the damage classification image data by tree vertex and geographical information, the number of damage classification number calculating step for calculating the number of red pine included in the whole area or arbitrary area of the red pine forest for each damage classification,
Based on the result of calculating the number of damage categories, the total number of damaged trees and infected trees in all areas or arbitrary areas of the red pine forest is included in all areas or arbitrary areas of the red pine forest. A damage rate classification diagram creation step of calculating a damage rate classification diagram by calculating the damage rate obtained by dividing by the total number of red pine,
A method for calculating a damage classification of a pine weevil characterized by having
請求項1に記載の松くい虫の被害区分算定方法において、
前記撮影画像データに含まれる前記樹種の分光反射特性に基づいて、前記撮影画像データからアカマツ林を抽出することによってアカマツ林画像データを作成するアカマツ林画像データ作成ステップをさらに有し、当該アカマツ林画像データ作成ステップを実施した後に、前記アカマツ林ごとに、前記アカマツ林樹頂点画像データ作成ステップと、前記アカマツ林樹冠画像データ作成ステップと、前記被害区分画像データ作成ステップと、前記樹頂点別被害区分画像データ作成ステップと、前記被害区分本数算定ステップと、前記被害率区分図作成ステップと、を実施すること特徴とする松くい虫の被害区分算定方法。
In the method for calculating the damage category of a pine weevil according to claim 1,
And further comprising a red pine forest image data creating step for creating red pine forest image data by extracting red pine forest from the photographed image data based on the spectral reflection characteristics of the tree species included in the photographed image data, After performing the image data creation step, for each of the red pine forests, the red pine forest tree vertex image data creation step, the red pine forest crown image data creation step, the damage classification image data creation step, and the damage by tree vertex A method for calculating a damage category of pine worms, comprising performing a step of creating a category image data, a step of calculating the number of damage categories, and a step of creating a damage rate category diagram.
請求項1又は2に記載の松くい虫の被害区分算定方法において、
前記撮影画像データは、秋から冬に撮影された画像データを用いることを特徴とする松くい虫の被害区分算定方法。
In the method for calculating the damage category of the pine beetle according to claim 1 or 2,
A method for calculating the damage classification of a pine weevil, wherein the photographed image data is image data photographed from autumn to winter.
請求項1〜3のいずれかに記載の松くい虫の被害区分算定方法において、
前記感染木、枯損木及び健全木のうちの前記感染木と前記枯損木との区分は、前記松くい虫の被害を受けているアカマツから採取された被害葉と健全なアカマツから採取された健全葉の各分光反射特性に基づいて行い、
前記感染木は、可視光の赤色域において前記枯損木及び前記健全木に比べて反射率が高いことで、前記枯損木及び前記健全木と区分し、
前記枯損木は、近赤外域において前記感染木及び前記健全木に比べて反射率が低いことで、前記感染木及び前記健全木と区分する、
ことを特徴とする松くい虫の被害区分算定方法。
In the method for calculating the damage classification of pine beetles according to any one of claims 1 to 3,
Among the infected tree, the dead tree and the healthy tree, the infected tree and the dead tree are classified into a damaged leaf collected from the red pine damaged by the pine weevil and a healthy tree collected from the healthy red pine. Based on the spectral reflection characteristics of the leaves,
The infected tree has a higher reflectance than the dead tree and the healthy tree in the visible red region, and is classified as the dead tree and the healthy tree,
The dead tree has a lower reflectance than the infected tree and the healthy tree in the near-infrared region, and is classified from the infected tree and the healthy tree.
This is a method for calculating the damage classification of pine weevil.
請求項1〜4のいずれかに記載の松くい虫の被害区分算定方法において、
前記アカマツ林樹頂点画像データ作成ステップは、輝度値が極大値を持つ点を自動抽出する局所最大値フィルターを用いてアカマツの樹頂点を抽出することを特徴とする松くい虫の被害区分算定方法。
In the method for calculating the damage classification of a pine beetle according to any one of claims 1 to 4,
The pine pine forest tree vertex image data creating step includes extracting a pine tree tree vertex using a local maximum filter that automatically extracts a point having a maximum luminance value. .
請求項1〜5のいずれかに記載の松くい虫の被害区分算定方法において、
前記アカマツ林樹冠画像データ作成ステップは、樹冠の縁部が樹冠部分よりも暗いことを利用して樹冠を抽出することを特徴とする松くい虫の被害区分算定方法。
In the method for calculating the damage classification of a pine beetle according to any one of claims 1 to 5,
A method for calculating a damage classification of a pine weevil, wherein the step of creating the image data of a red pine forest canopy extracts the crown using the fact that the edge of the crown is darker than the crown portion.
請求項1〜6のいずれかに記載の松くい虫の被害区分算定方法において、
前記樹冠別被害区分画像データ作成ステップは、前記被害区分を感染木、枯損木及び健全木の3種類として、前記アカマツの各樹冠から感染木、枯損木及び健全木が抽出された場合、前記各樹冠から抽出された感染木、枯損木及び健全木のそれぞれの占有面積を算出し、算出された感染木、枯損木及び健全木の各占有面積のうち、占有面積が大きい被害区分を当該樹冠の被害区分とみなすことを特徴とする松くい虫の被害区分算定方法。
In the method for calculating the damage classification of a pine beetle according to any one of claims 1 to 6,
The damage classification image data creation step according to the tree crown includes the damage classification as an infected tree, a dead tree, and a healthy tree, and when the infected tree, the dead tree, and the healthy tree are extracted from each crown of the red pine, Calculate the occupation area of each infected tree, dead tree, and healthy tree extracted from the tree canopy, and select the damage category with the largest occupation area of the calculated infected tree, dead tree, and healthy tree. A method for calculating the damage category of a pine beetle characterized by being regarded as a damage category.
請求項1〜7のいずれかに記載の松くい虫の被害区分算定方法において、
前記樹頂点別被害区分画像データ作成ステップは、前記地理的情報に基づいて、前記アカマツ林樹頂点画像データを前記樹冠別被害区分画像データに重ねることによって、樹頂点ごとに、前記被害区分を特定することを特徴とする松くい虫の被害区分算定方法。
In the method for calculating the damage classification of a pine beetle according to any one of claims 1 to 7,
The damage classification image data creation step for each tree vertex identifies the damage classification for each tree vertex based on the geographical information by superimposing the vertex image data for the red pine forest tree on the damage classification image data for each tree crown A method for calculating the damage classification of a pine weevil characterized by:
請求項1〜8のいずれかに記載の松くい虫の被害区分算定方法において、
前記被害区分本数算定ステップは、前記樹頂点別被害区分画像データに基づいて、前記アカマツ林の全領域又は任意の領域において前記地理的情報により小班と定義される森林所有者の境界データから小班区画を抽出し、抽出された小班において前記被害区分ごとに樹木の本数を算定することを特徴とする松くい虫の被害区分算定方法。
In the method for calculating the damage classification of a pine beetle according to any one of claims 1 to 8,
The step of calculating the number of damage categories is based on the submarine data from the forest owner's boundary data defined by the geographical information in the entire area of the red pine forest or an arbitrary area based on the damage classification image data classified by tree vertices. And calculating the number of trees for each of the damage categories in the extracted subcompartment.
調査対象となるアカマツ林を含む地域を上空から撮影して得られた撮影画像データに含まれる樹種の分光反射特性に基づいて作成されたアカマツ林画像データを用いて、前記アカマツ林に含まれるアカマツの本数を被害区分別に算定する松くい虫の被害区分算定装置であって、
前記アカマツ林画像データに基づいて、前記アカマツの各樹頂点を抽出し、前記各樹頂点に関するアカマツ林樹頂点画像データを作成するアカマツ林樹頂点画像データ作成部と、
前記アカマツ林画像データに基づいて、前記アカマツの各樹冠を抽出し、前記各樹冠に関するアカマツ林樹冠画像データを作成するアカマツ林樹冠画像データ作成部と、
前記アカマツ林樹冠画像データと、松くい虫の被害を受けているアカマツから採取された被害葉及び健全なアカマツから採取された健全葉の各分光反射特性とに基づいて、前記各樹冠の領域における感染木、枯損木及び健全木の分光反射特性を作成し、作成した感染木、枯損木及び健全木の分光反射特性に基づいて樹冠別の被害区分を特定して、樹冠別被害区分画像データを作成する樹冠別被害区分画像データ作成部と、
前記樹冠別被害区分画像データと前記アカマツ林樹頂点画像データとに基づいて、樹頂点別に前記被害区分を特定して、樹頂点別被害区分画像データを作成する樹頂点別被害区分画像データ作成部と、
前記樹頂点別被害区分画像データと地理的情報とに基づいて、前記アカマツ林の全領域又は任意の領域に含まれるアカマツの本数を前記被害区分ごとに算定する被害区分本数算定部と、
前記被害区分本数算定結果に基づいて、当該アカマツ林の全領域又は任意の領域に含まれる感染木と枯損木とを合計した被害木本数を、当該アカマツ林の全領域又は任意の領域に含まれるアカマツの全体本数で除することで得られる被害率を計算して被害率区分図を作成する被害率区分図作成部と、
を有すること特徴とする松くい虫の被害区分算定装置。
Using the red pine forest image data created based on the spectral reflection characteristics of the tree species included in the photographed image data obtained by photographing the region including the red pine forest to be surveyed from above, the red pine forest included in the red pine forest is used. Is a pine weevil damage category calculation device that calculates the number of
Based on the red pine forest image data, each of the red pine tree vertices is extracted, and a red pine forest tree vertex image data creation unit for creating red pine forest tree vertex image data relating to each tree vertex;
Based on the red pine forest image data, each of the red pine canopies is extracted, and a red pine forest canopy image data creation unit for creating red pine forest canopy image data relating to each of the canopies;
Based on the image data of the red pine forest canopy and the spectral reflection characteristics of the damaged leaves collected from the red pine and the healthy leaves collected from the healthy red pine, Create spectral reflection characteristics of infected, dead, and healthy trees, identify damage classifications by tree crown based on the spectral reflection characteristics of the infected, dead, and healthy trees, and generate damage classification image data by crown. The damage classification image data creation part by tree to create,
Based on the damage classification image data by tree crown and the vertex image data of the red pine forest tree, the damage classification image data generation unit by tree apex that identifies the damage classification by tree vertex and creates damage classification image data by tree vertex When,
Based on the damage classification image data and geographical information according to the tree vertices, a damage classification number calculating unit that calculates the number of red pine contained in the whole area of the red pine forest or an arbitrary area for each damage classification,
Based on the result of calculating the number of damage categories, the total number of damaged trees and infected trees in all areas or arbitrary areas of the red pine forest is included in all areas or arbitrary areas of the red pine forest. A damage rate classification diagram creation unit that calculates the damage rate obtained by dividing by the total number of red pine and creates a damage rate classification diagram;
A device for calculating the damage classification of pine worms, characterized by comprising:
請求項10に記載の松くい虫の被害区分算定装置において、
前記樹種の分光反射特性に基づいて前記撮影画像データからアカマツ林を抽出することによってアカマツ林画像データを作成するアカマツ林画像データ作成部をさらに有し、当該アカマツ林画像データ作成部によってアカマツ林画像データの作成を実施した後に、
前記アカマツ林ごとに、前記アカマツ林樹頂点画像データ作成部による前記アカマツ林樹頂点画像データの作成と、前記アカマツ林樹冠画像データ作成部による前記アカマツ林樹冠画像データの作成と、前記樹冠別被害区分画像データ作成部による前記樹冠別被害区分画像データの作成と、前記樹頂点別被害区分画像データ作成部による前記樹頂点別被害区分画像データ作成と、前記被害区分本数算定部による前記被害区分本数算定と、前記被害率区分図作成部による前記被害率区分図の作成とを実施することを特徴とする松くい虫の被害区分算定装置。
In the pine weevil damage classification calculating apparatus according to claim 10,
It further comprises a red pine forest image data creation unit that creates red pine forest image data by extracting red pine forest from the photographed image data based on the spectral reflection characteristics of the tree species, and the red pine forest image data creation unit further comprises a red pine forest image data creation unit. After creating the data,
For each of the red pine forests, the red pine forest tree vertex image data creation unit by the red pine forest tree vertex image data creation unit, the red pine forest crown image data creation unit by the red pine forest crown image data creation unit, and the damage by the crown Creation of the damage classification image data for each crown by the classification image data creation section, creation of the damage classification image data for each tree vertex by the damage classification image data creation section by the tree apex, and the number of damage classification by the damage classification number calculation section An apparatus for calculating the damage category of pine worms, characterized in that the calculation and the creation of the damage rate category map by the damage rate category diagram creation unit are performed.
請求項10又は11に記載の松くい虫の被害区分算定装置において、
前記撮影画像データは、秋から冬に撮影された画像データを用いることを特徴とする松くい虫の被害区分算定装置。
In the pine weevil damage classification calculation apparatus according to claim 10 or 11,
An apparatus for calculating the damage category of a pine weevil, wherein the photographed image data uses image data photographed from autumn to winter.
請求項10〜12のいずれかに記載の松くい虫の被害区分算定装置において、
前記感染木、枯損木及び健全木のうちの前記感染木と前記枯損木との区分は、前記松くい虫の被害を受けているアカマツから採取された被害葉と健全なアカマツから採取された健全葉の各分光反射特性に基づいて行い、
前記感染木は、可視光の赤色域において前記枯損木及び前記健全木に比べて反射率が高いことで、前記枯損木及び前記健全木と区分し、
前記枯損木は、近赤外域において前記感染木及び前記健全木に比べて反射率が低いことで、前記感染木及び前記健全木と区分する、
ことを特徴とする松くい虫の被害区分算定装置。
In the pine beetle damage classification calculation apparatus according to any one of claims 10 to 12,
Among the infected tree, the dead tree and the healthy tree, the infected tree and the dead tree are classified into a damaged leaf collected from the red pine damaged by the pine weevil and a healthy tree collected from the healthy red pine. Based on the spectral reflection characteristics of the leaves,
The infected tree has a higher reflectance than the dead tree and the healthy tree in the visible red region, and is classified as the dead tree and the healthy tree,
The dead tree has a lower reflectance than the infected tree and the healthy tree in the near-infrared region, and is classified from the infected tree and the healthy tree.
An apparatus for calculating the damage classification of pine weevil, which is characterized by that.
請求項10〜13のいずれかに記載の松くい虫の被害区分算定装置において、
前記アカマツ林樹頂点画像データ作成部は、輝度値が極大値を持つ点を自動抽出する局所最大値フィルターを用いて樹頂点を抽出することを特徴とする松くい虫の被害区分算定装置。
In the pine beetle damage classification calculation apparatus according to any one of claims 10 to 13,
The Pinus densiflora damage category calculation apparatus, wherein the Pinus densiflora vertices image data creation unit extracts tree vertices using a local maximum value filter that automatically extracts a point having a maximum luminance value.
請求項10〜14のいずれかに記載の松くい虫の被害区分算定装置において、
前記アカマツ林樹冠画像データ作成部は、樹冠の縁部が樹冠部分よりも暗いことを利用して樹冠を抽出することを特徴とする松くい虫の被害区分算定装置。
In the pine weevil damage classification calculation apparatus according to any one of claims 10 to 14,
The pine pine beetle damage classification calculating apparatus, wherein the red pine forest canopy image data creation unit extracts a tree crown using the fact that the edge of the tree crown is darker than the crown portion.
請求項10〜15のいずれかに記載の松くい虫の被害区分算定装置において、
前記樹冠別被害区分画像データ作成ステップは、前記被害区分を感染木、枯損木及び健全木の3種類として、前記アカマツの各樹冠から感染木、枯損木及び健全木が抽出された場合、前記各樹冠から抽出された感染木、枯損木及び健全木のそれぞれの占有面積を算出し、算出された感染木、枯損木及び健全木の各占有面積のうち、占有面積が大きい被害区分を当該樹冠の被害区分とみなすことを特徴とする松くい虫の被害区分算定装置。
In the pine beetle damage classification calculation apparatus according to any one of claims 10 to 15,
The damage classification image data creation step according to the tree crown includes the damage classification as an infected tree, a dead tree, and a healthy tree, and when the infected tree, the dead tree, and the healthy tree are extracted from each crown of the red pine, Calculate the occupation area of each infected tree, dead tree, and healthy tree extracted from the tree canopy, and select the damage category with the largest occupation area of the calculated infected tree, dead tree, and healthy tree. An apparatus for calculating the damage category of pine beetles characterized by being regarded as a damage category.
請求項10〜16のいずれかに記載の松くい虫の被害区分算定装置において、
前記樹頂点別被害区分画像データ作成ステップは、前記地理的情報に基づいて、前記アカマツ林樹頂点画像データを前記樹冠別被害区分画像データに重ねることによって、樹頂点ごとに、前記被害区分を特定することを特徴とする松くい虫の被害区分算定装置。
In the pine weevil damage classification calculation apparatus according to any one of claims 10 to 16,
The damage classification image data creation step for each tree vertex identifies the damage classification for each tree vertex based on the geographical information by superimposing the vertex image data for the red pine forest tree on the damage classification image data for each tree crown A device for calculating damage classification of pine weevil characterized by
請求項10〜17のいずれかに記載の松くい虫の被害区分算定装置において、
前記被害区分本数算定部は、前記樹頂点別被害区分画像データに基づいて、前記アカマツ林の全領域又は任意の領域において前記地理的情報により小班と定義される森林所有者の境界データから小班区画を抽出し、抽出された小班において前記被害区分ごとに樹木の本数を算定することを特徴とする松くい虫の被害区分算定装置。
In the pine beetle damage classification calculation apparatus according to any one of claims 10 to 17,
The damage classification number calculation unit is configured to calculate subdivisions from the boundary data of forest owners defined as subdivisions by the geographical information in the whole area of the pine forest or in any area based on the damage classification image data classified by tree vertices. And calculating the number of trees for each damage category in the extracted sub-compartment.
JP2016054551A 2016-03-17 2016-03-17 Calculation method of damage classification of pine pileworm and damage classification calculation device of pine pileworm Active JP6544582B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016054551A JP6544582B2 (en) 2016-03-17 2016-03-17 Calculation method of damage classification of pine pileworm and damage classification calculation device of pine pileworm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016054551A JP6544582B2 (en) 2016-03-17 2016-03-17 Calculation method of damage classification of pine pileworm and damage classification calculation device of pine pileworm

Publications (2)

Publication Number Publication Date
JP2017163934A true JP2017163934A (en) 2017-09-21
JP6544582B2 JP6544582B2 (en) 2019-07-17

Family

ID=59908548

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016054551A Active JP6544582B2 (en) 2016-03-17 2016-03-17 Calculation method of damage classification of pine pileworm and damage classification calculation device of pine pileworm

Country Status (1)

Country Link
JP (1) JP6544582B2 (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109360376A (en) * 2018-09-13 2019-02-19 航天信德智图(北京)科技有限公司 Infect the withered masson pine satellite monitoring spatial database of pine nematode
CN112577907A (en) * 2020-11-18 2021-03-30 上海市园林科学规划研究院 Urban green land tree crown loss rate calculation method
CN112753456A (en) * 2020-12-30 2021-05-07 山东农业大学 Accurate prevention and control method and system for pine wood nematode disease based on space-time law
WO2021187231A1 (en) * 2020-03-16 2021-09-23 株式会社電通国際情報サービス Affected tree diagnosis system, affected tree diagnosis method, and affected tree diagnosis program
CN114051853A (en) * 2020-07-31 2022-02-18 华南农业大学 Method for rapidly identifying high-sensitivity pine wood nematode disease masson pine in vitro spring shoots and application of method
KR20230130895A (en) * 2022-03-04 2023-09-12 충남대학교산학협력단 Method and apparatus for evaluating activity vitality of conifer based on infection symptom of wilt disease
CN116824200A (en) * 2022-11-24 2023-09-29 武汉很多鱼钓具有限公司 Forestry pest intelligent identification detection method based on target detection technology
CN116977865A (en) * 2023-09-22 2023-10-31 北京林业大学 Pine wood nematode disease risk prediction method based on single wood scale

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011103098A (en) * 2009-11-12 2011-05-26 Shinshu Univ Tree counting method and tree counting apparatus

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011103098A (en) * 2009-11-12 2011-05-26 Shinshu Univ Tree counting method and tree counting apparatus

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
平田 泰雅: "高分解能衛星データの森林モニタリングへの利用", 日本森林学会誌, vol. 2008 年 91 巻 2 号, JPN7019001104, 2008, JP, pages 136 - 146, ISSN: 0004028669 *
油井 竜太, 龍原 哲, 阿部 信行: "LandsatTMデータを使用したパターン展開法によるマツ枯れ被害林分の抽出", 日本森林学会誌, vol. 2001 年 83 巻 4 号, JPN7019001101, 2001, JP, pages 299 - 306, ISSN: 0004028668 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109360376A (en) * 2018-09-13 2019-02-19 航天信德智图(北京)科技有限公司 Infect the withered masson pine satellite monitoring spatial database of pine nematode
WO2021187231A1 (en) * 2020-03-16 2021-09-23 株式会社電通国際情報サービス Affected tree diagnosis system, affected tree diagnosis method, and affected tree diagnosis program
CN114051853A (en) * 2020-07-31 2022-02-18 华南农业大学 Method for rapidly identifying high-sensitivity pine wood nematode disease masson pine in vitro spring shoots and application of method
CN114051853B (en) * 2020-07-31 2023-01-10 华南农业大学 Method for rapidly identifying high-sensitivity pine wood nematode disease masson pine in vitro spring shoots and application of method
CN112577907A (en) * 2020-11-18 2021-03-30 上海市园林科学规划研究院 Urban green land tree crown loss rate calculation method
CN112753456A (en) * 2020-12-30 2021-05-07 山东农业大学 Accurate prevention and control method and system for pine wood nematode disease based on space-time law
CN112753456B (en) * 2020-12-30 2021-12-21 山东农业大学 Accurate prevention and control method and system for pine wood nematode disease based on space-time law
KR20230130895A (en) * 2022-03-04 2023-09-12 충남대학교산학협력단 Method and apparatus for evaluating activity vitality of conifer based on infection symptom of wilt disease
KR102646217B1 (en) 2022-03-04 2024-03-08 충남대학교산학협력단 Method and apparatus for evaluating activity vitality of conifer based on infection symptom of wilt disease
CN116824200A (en) * 2022-11-24 2023-09-29 武汉很多鱼钓具有限公司 Forestry pest intelligent identification detection method based on target detection technology
CN116977865A (en) * 2023-09-22 2023-10-31 北京林业大学 Pine wood nematode disease risk prediction method based on single wood scale
CN116977865B (en) * 2023-09-22 2024-02-06 北京林业大学 Pine wood nematode disease risk prediction method based on single wood scale

Also Published As

Publication number Publication date
JP6544582B2 (en) 2019-07-17

Similar Documents

Publication Publication Date Title
JP6544582B2 (en) Calculation method of damage classification of pine pileworm and damage classification calculation device of pine pileworm
Navarro et al. The application of Unmanned Aerial Vehicles (UAVs) to estimate above-ground biomass of mangrove ecosystems
Chong et al. A review of remote sensing applications for oil palm studies
Mucher et al. Land cover characterization and change detection for environmental monitoring of pan-Europe
Lucas et al. Classification of Australian forest communities using aerial photography, CASI and HyMap data
Rapinel et al. Combined use of LiDAR data and multispectral earth observation imagery for wetland habitat mapping
JP4900356B2 (en) Tree species classification method and tree species classification system, forest status information creation method and forest status information creation system, thinning target area selection method and thinning target area selection system
Kim et al. Distinguishing between live and dead standing tree biomass on the North Rim of Grand Canyon National Park, USA using small-footprint lidar data
JP4858793B2 (en) Tree number calculation method and tree number calculation device
Shouse et al. Identification of understory invasive exotic plants with remote sensing in urban forests
Han et al. Extraction of multilayer vegetation coverage using airborne LiDAR discrete points with intensity information in urban areas: A case study in Nanjing City, China
Zharikova et al. Methods of remote sensing in forest fund assessment problems
Gómez‐Sapiens et al. Improving the efficiency and accuracy of evaluating aridland riparian habitat restoration using unmanned aerial vehicles
Stoyanova et al. Potential of multispectral imaging technology for assessment coniferous forests bitten by a bark beetle in Central Bulgaria
Kopeć et al. How to effectively use long-term remotely sensed data to analyze the process of tree and shrub encroachment into open protected wetlands.
Fryskowska et al. Effective detection of sub-surface archeological features from laser scanning point clouds and imagery data
Jusoff Sustainable management of a matured oil palm plantation in UPM campus, Malaysia using airborne remote sensing
Stone et al. Integrating plantation health surveillance and wood resource inventory systems using remote sensing
Jaya et al. Spatial analysis of vegetation density classification in determining environmental impacts using uav imagery
Shi A Remote Sensing and GIS-based Wetland Analysis In Canaan Valley, West Virginia
Johenneken et al. Towards multi-class forest floor analysis
Burke Shelterbelt density dynamics and their driving forces in Grand Forks County, North Dakota, 1962 to 2014
Fadaei et al. Estimation of tree density with high-resolution imagery in the zarbin forest of North Iran (Cupressus sempervirence var. horzontalis)
Steinberg Assessment of hazard tree/snag detection using drone-based, multi-spectral sensors
Apan et al. Mapping olive varieties and within-field spatial variability using high resolution QuickBird imagery

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20181226

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20181228

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20181227

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190404

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20190405

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190409

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20190409

TRDD Decision of grant or rejection written
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20190409

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20190514

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190605

R150 Certificate of patent or registration of utility model

Ref document number: 6544582

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250