JP2010086276A - Method and system for tree classification, method and system for generating current status information of forest, and method and system for selecting target area for tree thinning - Google Patents

Method and system for tree classification, method and system for generating current status information of forest, and method and system for selecting target area for tree thinning Download PDF

Info

Publication number
JP2010086276A
JP2010086276A JP2008254520A JP2008254520A JP2010086276A JP 2010086276 A JP2010086276 A JP 2010086276A JP 2008254520 A JP2008254520 A JP 2008254520A JP 2008254520 A JP2008254520 A JP 2008254520A JP 2010086276 A JP2010086276 A JP 2010086276A
Authority
JP
Japan
Prior art keywords
image data
forest
tree
area
canopy
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
JP2008254520A
Other languages
Japanese (ja)
Other versions
JP4900356B2 (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 JP2008254520A priority Critical patent/JP4900356B2/en
Publication of JP2010086276A publication Critical patent/JP2010086276A/en
Application granted granted Critical
Publication of JP4900356B2 publication Critical patent/JP4900356B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

<P>PROBLEM TO BE SOLVED: To properly select target area for tree thinning by generating current status information by classifying trees properly and using the same. <P>SOLUTION: The system for selecting target area for tree thinning includes: a unit 310 for generating sub-compartment image data of current status of forest, the unit 310 generating sub-compartment image data of current status of the forest based on sub-compartment image data and the image data of the current status information of the forest, wherein the sub-compartment image data for dividing a forest into sub-compartments is generated based on image data of the current status information of the forest relating to the current status information of the forest; a unit 320 for generating image data of sub-compartment attribute information, the unit 320 extracting sub-compartment attribute information representing features of each sub-compartment and generating image data of sub-compartment attribute information relating to each sub-compartment attribute information; and a unit 330 for generating image data of target area for tree thinning, the unit 330 selecting target area for tree thinning based on the image data of the sub-compartment attribute information and geographical advantage of the forest and generating an image data of the target area for tree thinning representing the selected target area for tree thinning. <P>COPYRIGHT: (C)2010,JPO&INPIT

Description

本発明は、森林の樹種を分類する樹種分類方法及び樹種分類システム、森林現況情報の作成方法及び森林現況情報の作成システム、間伐対象区域の選定方法及び間伐対象区域の選定システムに関する。   The present invention relates to a tree species classification method and a tree species classification system for classifying forest tree species, a forest present state information creation method and a forest present state information creation system, a thinning target area selection method, and a thinning target area selection system.

森林の現況を正確に捉えることによって間伐すべき区域(間伐対象区域という)を適切に選定することは森林管理おいて極めて重要な課題である。一般に、森林は災害や気象害、大型獣による食害、不要な樹木の侵入や天然更新などの影響を受け、多くの森林は、植栽した樹種がそのまま育つことはなく、年月の経過とともにその状況は大きく変化する。   Proper selection of areas to be thinned (called thinning target areas) by accurately grasping the current state of forests is an extremely important issue in forest management. In general, forests are affected by disasters, meteorological damage, feeding damage by large animals, invasion of unnecessary trees, natural regeneration, etc., and many forests do not grow as planted trees as they are. The situation changes greatly.

しかし、森林の管理台帳である森林基本図や森林調査簿は植栽時の情報のままであり、現状との乖離が大きな問題である。さらに、近年、人工林では災害の起きやすい不健全な森林が増加している。このような森林は間伐対象区域として人の手によって間伐を行う必要がある。   However, the basic forest map and the forest survey book, which are the forest management ledgers, remain the information at the time of planting, and the deviation from the current situation is a big problem. Furthermore, in recent years, unhealthy forests that are prone to disasters are increasing in planted forests. Such forests need to be thinned by hand as a thinning target area.

しかしながら、間伐対象区域の選定は、従来は現場調査と技術者の経験に依存しており、多大な労力と時間を要する。また、森林は一般的には地形が急峻で面積が広く、アクセスが悪いことから、労働作業の関係から、人の目に見える範囲や林道の近くの森林が間伐地として選定されることが多く、間伐対象区域が適切に選定されているとは言い難い。   However, selection of areas to be thinned depends on field surveys and experience of engineers in the past, and requires a lot of labor and time. In addition, forests generally have steep terrain, large areas, and poor access. Therefore, forests that are visible to humans or near forest roads are often selected as thinning areas because of labor. It is hard to say that the areas to be thinned are properly selected.

これは、森林の現況を正確で捉え、間伐すべき箇所を客観的に判別することは極めて困難であるからである。すなわち、森林の現況を正確に捉えるには、森林の現況を表す情報(森林現況情報という)を適切に作成する必要がある。ここでいう森林現況情報というのは、例えば、当該森林の樹冠に関する情報、樹種に関する情報、樹冠の面積、全樹種の本数や樹種別の本数などが含まれる。   This is because it is extremely difficult to accurately grasp the current state of the forest and objectively determine the location to be thinned. That is, in order to accurately grasp the current state of the forest, it is necessary to appropriately create information representing the current state of the forest (referred to as forest present state information). The forest status information here includes, for example, information on the crown of the forest, information on the tree species, the area of the canopy, the number of all tree species, the number of tree types, and the like.

森林現況情報を取得して森林を管理する技術は従来から提案されている(例えば、特許文献1及び特許文献2参照。)。   Techniques for acquiring forest status information and managing forests have been proposed (see, for example, Patent Document 1 and Patent Document 2).

特許文献1に開示された技術(第1従来技術という。)は、調査対象となる森林を上空から撮影した撮影画像(カラー画像)を白黒画像に変換し、該白黒画像を画像処理して得られた樹冠の形状と元のカラー画像の分光反射(スペクトル)分析を行うことによって得られる樹木の色彩・輝度などとから、森林を構成する樹種など森林の植生の調査・評価を行うというものである。   The technique disclosed in Patent Document 1 (referred to as the first conventional technique) is obtained by converting a captured image (color image) obtained by photographing a forest to be investigated from the sky into a black and white image, and performing image processing on the black and white image. The tree vegetation, such as the tree species that make up the forest, is investigated and evaluated based on the shape of the tree crown and the color and brightness of the tree obtained by spectral reflection (spectrum) analysis of the original color image. is there.

また、特許文献2に開示された技術(第2従来技術という。)は、調査対象となる森林地域を上空から撮影した撮影画像を画像処理して、森林の樹冠を評価する技術であって、この第2従来技術には、樹冠の形状を求めたのち、当該樹冠の形状のテクスチャ特徴量を求め、求められたテクスチャ特徴量と既知の樹種のテクスチャ特徴量とを比較して当該樹冠の樹種を特定する樹種分類技術が開示されている。   In addition, the technique disclosed in Patent Document 2 (referred to as second conventional technique) is a technique for evaluating a forest canopy by performing image processing on a captured image obtained by photographing a forest area to be investigated from the sky, In the second prior art, after obtaining the shape of the crown, the texture feature amount of the shape of the crown is obtained, and the obtained texture feature amount is compared with the texture feature amount of a known tree species to determine the tree species of the crown. A tree species classification technique for identifying a tree is disclosed.

特開2001−357380号公報JP 2001-357380 A 特開2006−285310号公報JP 2006-285310 A

第1従来技術は、上記したように、カラー画像のスペクトル分析を行うことによって得られる樹木の色彩及び輝度などから、森林の樹種などの調査など森林の植生の調査・評価を行うようにしている。確かに、樹種によって光の波長に対する輝度の違いが生じるので、樹種の特定が可能となると考えられるが、第1従来技術においては、日照の状況(例えば、日向部と日陰部など)が考慮されていないので、高精度な分類はできす、誤分類が生じやすいという問題がある。   As described above, according to the first prior art, forest vegetation is investigated and evaluated, such as a survey of forest tree species, from the color and brightness of trees obtained by performing spectrum analysis of color images. . Certainly, the difference in brightness with respect to the wavelength of light occurs depending on the tree species, so it is considered that the tree species can be specified, but in the first conventional technique, the situation of sunlight (for example, the sun and shade areas) is considered. Therefore, there is a problem that high-precision classification is possible and misclassification is likely to occur.

また、第2従来技術は、樹種の特定を行う際、樹冠の形状のテクスチャ特徴量を求め、求められたテクスチャ特徴量と既知の樹種のテクスチャ特徴量とを比較して当該樹冠の樹種を特定するようにしている。このため、第2従来技術においては、既知の樹種のテクスチャ特徴量を個々の樹種ごとに予め作成しておく必要があるとともに、予め用意された既知の樹種のテクスチャ特徴量と求められたテクスチャ特徴量とを比較するといった演算量の多い処理を行う必要がある。また、樹種によってはテクスチャ特徴量が類似する場合もあり、第1従来技術と同様、誤分類が生じやすいという問題がある。   In addition, when the second conventional technology specifies a tree species, the texture feature amount of the shape of the crown is obtained, and the obtained texture feature amount is compared with the texture feature amount of a known tree species to identify the tree species of the tree crown. Like to do. For this reason, in the second prior art, it is necessary to create a texture feature amount of a known tree species in advance for each tree species, and a texture feature amount of a known tree species prepared in advance and a texture feature obtained. It is necessary to perform processing with a large amount of calculation such as comparing the amount. Also, depending on the tree type, the texture feature amount may be similar, and there is a problem that misclassification is likely to occur as in the first conventional technique.

そこで本発明は、樹木を樹種別に適切に分類可能とする樹種分類方法及び樹種分類システムを提供するとともに森林の現況情報を適切に取得可能な森林現況情報の作成方法及び森林現況情報の作成システムを提供することを目的とする。また、森林現況情報を用いて間伐対象区域を適切に選定可能とする間伐対象区域の選定方法及び間伐対象区域の選定システムを提供することを目的とする。   Accordingly, the present invention provides a tree species classification method and a tree species classification system that can appropriately classify trees according to tree types, and a forest status information creation method and a forest status information creation system that can appropriately acquire forest status information. The purpose is to provide. Another object of the present invention is to provide a method for selecting a thinning target area and a selection system for the thinning target area so that the thinning target area can be appropriately selected using the current forest information.

(1)本発明の樹種分類方法は、調査対象となる森林を含む地域を上空から撮影して得られた撮影画像データと前記森林の地理的情報とに基づいて抽出された森林域に関する森林域画像データを用いて前記森林域に含まれる各樹木を種類別に分類する樹種分類方法であって、前記森林域画像データに基づいて前記森林域に含まれる各樹冠を抽出し、前記各樹冠に関する樹冠画像データを作成するステップと、前記樹冠画像データに基づいて前記各樹冠の日向部を含む情報を樹冠情報として抽出し、前記樹冠情報に関する樹冠情報画像データを作成するステップと、前記樹冠情報画像データに基づいて前記各樹冠の日向部とされた領域における分光反射特性を作成し、前記分光反射特性に基づいて当該樹冠の樹種を特定して樹種分類画像データを作成するステップとを有すること特徴とする。   (1) The tree species classification method of the present invention includes a forest area related to a forest area extracted based on photographed image data obtained by photographing an area including a forest to be investigated from the sky and geographical information of the forest. A tree classification method for classifying each tree included in the forest area by using image data, wherein each tree crown included in the forest area is extracted based on the forest area image data, A step of creating image data, a step of extracting information including a sunny part of each canopy based on the canopy image data as canopy information, creating canopy information image data relating to the canopy information, and the canopy information image data A spectral reflection characteristic in a region that is the sunny part of each canopy based on the tree, and identifying a tree species of the crown based on the spectral reflection characteristic to identify tree species classification image data And wherein further comprising the step of creating.

本発明の樹種分類方法によれば、各樹冠の日向部における分光反射特性を用いて樹木を樹種別に分類することにより、誤分類が生じにくく高い精度で分類することができる。   According to the tree species classification method of the present invention, by classifying trees into tree types using the spectral reflection characteristics at the sunny part of each crown, misclassification is unlikely to occur and classification can be performed with high accuracy.

(2)本発明の樹種分類方法においては、前記各樹冠の日向部は、前記樹冠画像データにおける輝度値に基づいて抽出され、前記輝度値が当該樹冠において所定値以上である箇所を前記日向部とすることが好ましい。   (2) In the tree species classification method of the present invention, the sunshine portion of each crown is extracted based on the luminance value in the crown image data, and the location where the luminance value is a predetermined value or more in the tree crown is the sunshine portion. It is preferable that

これにより、各樹冠における日向部を適切に抽出することができる。なお、前記樹冠画像データにおける輝度値により、日向部だけではなく、日陰部及び樹頂部をも抽出することが好ましい。これら日向部、日陰部及び樹頂部の抽出は、輝度値が当該樹冠において所定値以上(例えば、当該樹冠における輝度値の平均値以上)である箇所を日向部とし、輝度値が当該樹冠において所定値未満(例えば、当該樹冠における輝度値の平均値未満)である箇所を日陰部とし、また、輝度値が当該樹冠における極大値である箇所を樹頂部とすることによって可能となる。   Thereby, the sunny part in each crown can be extracted appropriately. In addition, it is preferable to extract not only the sunny part but also the shade part and the tree top part by the luminance value in the tree crown image data. In the extraction of the sunny part, the shade part, and the top part of the tree, a part having a luminance value equal to or greater than a predetermined value in the crown (for example, an average value of luminance values in the crown) is defined as the sunny part, and the luminance value is predetermined in the crown. This is possible by setting a portion that is less than the value (for example, less than the average value of luminance values in the tree crown) as a shaded portion, and a point where the luminance value is a maximum value in the tree crown as a tree top portion.

(3)本発明の樹種分類方法においては、前記各樹冠から複数の樹種が抽出された場合、前記各樹冠から抽出された複数の樹種のそれぞれの占有面積を各樹種ごとに算出し、占有面積の多い樹種を当該樹冠の樹種とみなすことが好ましい。   (3) In the tree species classification method of the present invention, when a plurality of tree species are extracted from each canopy, the occupied area of each of the plurality of tree species extracted from each canopy is calculated for each tree species, and the occupied area It is preferable to regard the tree species having a large number as the tree species of the crown.

これにより、1つの樹冠を1つ樹種として分類することができる。1つの樹冠が1つ樹種として分類されることにより、森林の植生の調査・評価などを容易かつ適切に行うことができる。   Thereby, one crown can be classified as one tree species. By classifying one tree crown as one tree species, it is possible to easily and appropriately conduct surveys and evaluations of forest vegetation.

(4)本発明の樹種分類システムは、調査対象となる森林を含む地域を上空から撮影して得られた撮影画像データと前記森林の地理的情報とに基づいて抽出された森林域に関する森林域画像データを用いて前記森林域に含まれる各樹木の種類を分類する樹種分類システムであって、前記森林域画像データに基づいて前記森林域に含まれる各樹冠を抽出し、前記各樹冠に関する樹冠画像データを作成する樹冠画像データ作成部と、前記樹冠画像データに基づいて前記各樹冠の日向部を含む情報を樹冠情報として抽出し、前記樹冠情報に関する樹冠情報画像データを作成する樹冠情報画像データ作成部と、前記樹冠情報画像データに基づいて前記各樹冠の日向部とされた領域における分光反射特性を作成し、前記分光反射特性に基づいて当該樹冠の樹種を特定して樹種分類画像データを作成する樹種分類画像データ作成部とを有すること特徴とする。   (4) The tree species classification system of the present invention is a forest area relating to a forest area extracted based on photographed image data obtained by photographing an area including a forest to be investigated from the sky and geographical information of the forest. A tree type classification system for classifying the types of trees included in the forest area using image data, wherein each tree crown included in the forest area is extracted based on the forest area image data, Crown image data creation unit for creating image data, and information including the sunny part of each crown based on the crown image data is extracted as crown information, and crown information image data for creating the crown information image data related to the crown information is created. Creating a spectral reflection characteristic in a region that is a sunny part of each canopy based on the creation information and the canopy information image data, and the crown based on the spectral reflection characteristic The species identified and characterized by having a tree species classification image data creation unit that creates a tree species classification image data.

本発明の樹種分類システムにおいても前記樹種分類方法と同様の効果を得ることができる。なお、本発明の樹種分類システムは、当該樹種分類システムに含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。また、本発明の樹種分類システムにおいても、前記(2)及び(3)に記載の特徴を有することが好ましい。   In the tree species classification system of the present invention, the same effect as the tree species classification method can be obtained. The tree species classification system of the present invention is configured such that the functions of each component included in the tree species classification system are installed as a computer program, and each component is provided with predetermined data. The functions of the elements are executed on the computer software. Also, the tree species classification system of the present invention preferably has the characteristics described in (2) and (3) above.

(5)本発明の森林現況情報の作成方法は、調査対象となる森林を含む地域を上空から撮影して得られた撮影画像データと前記森林の地理的情報とに基づいて抽出された森林域に関する森林域画像データを用いて前記森林の現況に関する森林現況情報を作成するための森林現況情報の作成方法であって、前記森林域画像データに基づいて前記森林域に含まれる各樹冠を抽出し、前記各樹冠に関する樹冠画像データを作成するステップと、前記樹冠画像データに基づいて前記各樹冠の日向部を含む情報を樹冠情報として抽出し、前記樹冠情報に関する樹冠情報画像データを作成するステップと、前記樹冠情報画像データに基づいて前記各樹冠の日向部とされた領域における分光反射特性を作成し、前記分光反射特性に基づいて当該樹冠の樹種を特定して樹種分類画像データを作成するステップと、前記樹種分類画像データに基づいて単位面積当たりの全樹種樹冠面積占有度を算出し、前記全樹種樹冠面積占有度に基づく全樹種樹冠占有画像データを作成するステップと、前記樹種分類画像データに基づいて単位面積当たりの樹種別樹冠面積占有度を算出し、前記樹種別樹冠面積占有度に基づく樹種別樹冠占有画像データを作成するステップと、前記樹種分類画像データに基づいて算出された全樹種の樹冠本数から単位面積当たりの全樹種樹冠本数分布画像データを作成するステップと、前記樹種分類画像データに基づいて算出された樹種別の樹冠本数から単位面積当たりの樹種別樹冠本数分布画像データを作成するステップと、前記全樹種樹冠占有画像データ、前記樹種別樹冠占有画像データ、前記全樹種樹冠本数分布画像データ、前記樹種別樹冠本数分布画像データ及び前記樹種分類画像データに基づいて森林の現況情報に関する森林現況情報画像データを作成するステップとを有することを特徴とする。   (5) The method for creating forest status information according to the present invention includes a forest area extracted based on photographed image data obtained by photographing an area including a forest to be surveyed from above and geographical information of the forest. A method for creating forest status information for creating forest status information relating to the current status of the forest using forest area image data relating to, wherein each crown included in the forest area is extracted based on the forest area image data A step of creating canopy image data relating to each canopy; a step of extracting information including the sunny part of each canopy based on the canopy image data as canopy information and creating canopy information image data relating to the canopy information; , Creating a spectral reflection characteristic in a region of the crown of each crown based on the crown information image data, and determining a tree species of the crown based on the spectral reflection characteristic A step of creating tree species classification image data, calculating a total tree species crown area occupancy per unit area based on the tree species classification image data, and calculating all tree species crown occupancy image data based on the total tree species crown area occupancy Calculating a tree type crown area occupancy per unit area based on the tree species classification image data, creating a tree type crown occupancy image data based on the tree type crown area occupancy, and From the number of crowns of all tree species calculated based on the tree species classification image data, the step of creating the total tree species crown number distribution image data per unit area, and the number of tree types of crowns calculated based on the tree species classification image data A step of creating tree-type tree crown distribution image data per unit area, the tree-type tree-occupied image data, the tree-type tree-occupied image Creating forest present state information image data relating to the present state information of forest based on the data, the total tree species crown number distribution image data, the tree type canopy number distribution image data and the tree species classification image data. .

このように、本発明の森林現況情報の作成方法は、各樹冠の日向部における分光反射特性を用いて樹種分類を行うことによって得られた高精度な樹種分類画像データを用いて、森林現況情報画像データを作成するようにしている。これにより、本発明の森林現況情報の作成方法によれば、森林の現況を適切に表した森林現況情報画像データを得ることができる。   As described above, the present forest status information creation method uses the high-accuracy tree species classification image data obtained by classifying the tree species using the spectral reflection characteristics in the sunny portion of each crown. Image data is created. Thus, according to the method for creating forest status information of the present invention, it is possible to obtain forest status information image data that appropriately represents the forest status.

(6)本発明の森林現況情報の作成システムは、調査対象となる森林を含む地域を上空から撮影して得られた撮影画像データと前記森林の地理的情報とに基づいて抽出された森林域に関する森林域画像データを用いて前記森林の現況に関する森林現況情報を作成するための森林現況情報の作成システムであって、前記森林域画像データに基づいて前記森林域に含まれる各樹冠を抽出し、前記各樹冠に関する樹冠画像データを作成する樹冠画像データ作成部と、前記樹冠画像データに基づいて前記各樹冠の日向部を含む情報を樹冠情報として抽出し、前記樹冠情報に関する樹冠情報画像データを作成する樹冠情報画像データ作成部と、前記樹冠情報画像データに基づいて前記各樹冠の日向部とされた領域における分光反射特性を作成し、前記分光反射特性に基づいて当該樹冠の樹種を特定して樹種分類画像データを作成する樹種分類画像データ作成部と、前記樹種分類画像データに基づいて全樹種の樹冠本数を算出する機能と、単位面積当たりの全樹種樹冠面積占有度を算出し、前記全樹種樹冠面積占有度に基づく全樹種樹冠占有画像データを作成する機能とを有する全樹種樹冠占有画像データ作成部と、前記樹種分類画像データに基づいて樹種別の樹冠本数を算出する機能と、単位面積当たりの樹種別樹冠面積占有度を算出し、前記樹種別樹冠面積占有度に基づく樹種別樹冠占有画像データを作成する機能とを有する樹種別樹冠占有画像データ作成部と、前記全樹種の樹冠本数から単位面積当たりの全樹種樹冠本数分布画像データを作成する全樹種樹冠本数分布画像データ作成部と、前記樹種別の樹冠本数から単位面積当たりの樹種別樹冠本数分布画像データを作成する樹種別樹冠本数分布画像データ作成部と、前記全樹種樹冠占有画像データ、前記樹種別樹冠占有画像データ、前記全樹種樹冠本数分布画像データ、前記樹種別樹冠本数分布画像データ及び前記樹種分類画像データに基づいて森林の現況情報に関する森林現況情報画像データを作成する森林現況情報画像データ作成部とを有することを特徴とする。   (6) The forest status information creation system according to the present invention is a forest area extracted based on photographed image data obtained by photographing an area including a forest to be investigated from the sky and geographical information of the forest. A forest present state information creation system for creating forest present state information relating to the present state of the forest using the forest region image data relating to, wherein each tree crown included in the forest region is extracted based on the forest region image data. A crown image data creation unit for creating crown image data relating to each crown, and extracting information including the sunny part of each crown based on the crown image data as crown information, and crown information image data relating to the crown information A crown information image data creation unit to be created, and a spectral reflection characteristic in a region that is a sunlit portion of each crown based on the crown information image data, A tree species classification image data creation unit that identifies tree species of the tree crown based on reflection characteristics and creates tree species classification image data, a function that calculates the number of tree crowns of all tree species based on the tree species classification image data, and per unit area A total tree species crown occupancy image data creation unit having a function of calculating all tree species crown area occupancy and creating all tree species crown occupancy image data based on the all tree species crown area occupancy, and based on the tree species classification image data A tree type having a function of calculating the number of crowns of a tree type and a function of calculating a tree type crown area occupancy per unit area and creating tree type crown occupancy image data based on the tree type crown area occupancy A crown occupancy image data creation unit; and a total tree species crown number distribution image data creation unit for creating all tree species crown number distribution image data per unit area from the number of crowns of all tree species; Tree-type crown number distribution image data creation unit for creating tree-type crown number distribution image data per unit area from the number of tree-type crowns, all-tree-type crown-occupied image data, said tree-type canopy-occupied image data, and all A forest present state information image data creation unit that creates forest present state information image data related to the present state of forest based on the tree species crown number distribution image data, the tree type crown number distribution image data, and the tree species classification image data; And

本発明の森林現況情報の作成システムにおいても前記森林現況情報の作成方法と同様の効果を得ることができる。なお、本発明の森林現況情報の作成システムは、当該森林現況情報の作成システムに含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   Also in the forest present status information creation system of the present invention, the same effect as the forest present status information creation method can be obtained. In addition, in the forest present status information creation system according to the present invention, the functions of each of the components included in the forest present status information creation system are installed as a computer program, and give predetermined data to the respective components. Thus, the functions of each component are executed on the software of the computer.

(7)本発明の間伐対象区域の選定方法は、調査対象となる森林の現況に関する森林現況情報画像データに基づいて間伐すべき間伐対象区域を選定する間伐対象区域の選定方法であって、前記森林現況情報画像データに基づいて前記森林を小班区画に分別するための小班区画画像データを作成して、前記小班区画画像データと前記森林現況情報画像データとに基づいて森林現況小班区画画像データを作成するステップと、前記各小班区画の特徴を現す小班区画属性情報を抽出し、前記各小班区画属性情報に関する小班区画属性情報画像データを作成するステップと、前記小班区画属性情報画像データと当該森林の地利条件とに基づいて間伐対象区域を選定し、選定された間伐対象区域を表す間伐対象区域画像データを作成するステップとを有することを特徴とする。   (7) The method for selecting a thinning target area of the present invention is a method for selecting a thinning target area for selecting a thinning target area to be thinned based on the forest present state information image data relating to the current state of the forest to be investigated, Based on the forest present state information image data, a subcompartment image data for dividing the forest into subcompartments is created, and based on the subcompartment image data and the forest present state information image data, the forest present state subcompartment image data is obtained. A step of creating, extracting a piece segment attribute information representing the characteristics of each group segment, creating segment segment attribute information image data relating to each segment segment attribute information, the segment segment attribute information image data and the forest Selecting a thinning target area based on the land condition and creating a thinning target area image data representing the selected thinning target area. And wherein the Rukoto.

このように、本発明の間伐対象区域の選定方法は、森林現況情報画像データに基づいて分別された各小班区画における各樹冠の特徴を現す小班区画属性情報を抽出し、抽出された各小班区画属性情報と当該森林の地理条件とに基づいて間伐対象区域を選定するようにしている。これにより、従来においては、個人の労力と熟練に頼っていた間伐対象区域の選定を、本発明の間伐対象区域の選定方法によれば、調査対象地域の中から客観的かつ適切に選定することできる。   As described above, the thinning target area selection method of the present invention extracts the subcompartment attribute information indicating the characteristics of each crown in each subcompartment sorted based on the forest present state information image data, and extracts each subcompartment extracted. The thinning target area is selected based on the attribute information and the geographical conditions of the forest. As a result, according to the method of selecting the thinning target area of the present invention, the selection of the thinning target area, which has conventionally relied on the labor and skill of the individual, should be selected objectively and appropriately from the survey target area. it can.

なお、間伐対象区域が選定されたとしても、実際に間伐すべき区域は、選定された間伐対象区域全体であるとは限らず、例えば、都道府県などで整備されている森林の管理台帳などに登録されている森林境界に基づいて、選定された間伐対象区域全体のうちの所定の区域内のみを実際の間伐対象区域とするということもできる。   Even if an area to be thinned is selected, the area to be thinned is not necessarily the entire area to be thinned. For example, in the forest management ledger established in prefectures, etc. Based on the registered forest boundary, it can also be said that only a predetermined area of the selected thinning target area is the actual thinning target area.

(8)本発明の間伐対象区域の選定方法においては、前記森林現況情報画像データは、前記(5)に記載の森林現況情報の作成方法によって作成された森林現況情報画像データであることが好ましい。   (8) In the method for selecting a thinning target area of the present invention, the forest present state information image data is preferably forest present state information image data created by the forest present state information creating method described in (5). .

このように、森林現況情報画像データとして、前記(5)に記載の森林現況情報の作成方法によって作成された森林現況情報画像データを用いることにより、間伐対象区域の選定を、より適切に行うことができる。すなわち、前記(5)に記載の森林現況情報の作成方法によって作成された森林現況情報画像データは、各樹冠の日向部における分光反射特性を用いて樹種分類を行うことによって得られた高精度な樹種分類画像データを用いて作成された森林現況情報画像データであり、森林の現況を適切に表した森林現況情報画像データである。このような森林現況情報画像データを用いて間伐対象区域の選定を行うことにより、間伐対象区域の選定をより適切に行うことができる。   As described above, by using the forest present state information image data created by the forest present state information creation method described in (5) above as the forest present state information image data, more appropriate selection of the thinning target area is performed. Can do. That is, the forest present state information image data created by the method for creating forest present state information described in (5) above is a high-accuracy obtained by classifying tree species using the spectral reflection characteristics in the sunny part of each crown. This is forest status information image data created using tree species classification image data, and is forest status information image data appropriately representing the current status of the forest. By selecting the area to be thinned using such forest status information image data, the area to be thinned can be selected more appropriately.

(9)本発明の間伐対象区域の選定方法においては、前記小班区画属性情報は、当該小班区画における樹種別の樹冠本数、単位面積あたりの樹種別の樹冠本数を含むことが好ましい   (9) In the method for selecting a thinning target area of the present invention, it is preferable that the subcompartment section attribute information includes the number of tree types of crowns in the subcompartment section and the number of tree types of crowns per unit area.

これにより、各小班区画の特徴をより適切に表すことができる。このような小班区画属性情報と当該森林の地利条件とに基づいて間伐対象区域を選定することにより、現時点における実際の森林の状況(森林現況として例えば、樹冠面積が大きい、占有面積割合が高い、立木本数が多いなど)が考慮され、かつ、間伐作業のし易さ(例えば、林道から近い、標高があまり高くない、傾斜が緩いなど)が考慮された間伐対象区域を選定することができる。   Thereby, the characteristic of each subdivision section can be expressed more appropriately. By selecting the area to be thinned based on such sub-compartment attribute information and the land condition of the forest, the actual forest situation at the present time (for example, a large canopy area, a high occupied area ratio, It is possible to select a thinning target area in consideration of the ease of thinning work (for example, close to a forest road, not so high in altitude, loose slope, etc.).

(10)本発明の間伐対象区域の選定方法においては、前記森林の地利条件は、当該小班区画の林道からの距離、標高、傾斜を含むことが好ましい。   (10) In the method for selecting a thinning target area of the present invention, it is preferable that the land condition of the forest includes a distance from the forest road of the subcompartment section, an altitude, and a slope.

これにより、間伐作業のし易さ(例えば、林道から近い、標高があまり高くない、傾斜が緩いなど)を適切に評価することができる。   Thereby, it is possible to appropriately evaluate the ease of thinning work (for example, close to the forest road, the altitude is not so high, the slope is gentle, etc.).

(11)本発明の間伐対象区域の選定システムは、調査対象となる森林の現況に関する森林現況情報画像データに基づいて間伐すべき間伐対象区域を選定する間伐対象区域の選定システムであって、前記森林現況情報画像データに基づいて前記森林を小班区画に分別するための小班区画画像データを作成して、前記小班区画画像データと前記森林現況情報画像データとに基づいて森林現況小班区画画像データを作成する森林現況小班区画画像データ作成部と、前記各小班区画の特徴を現す小班区画属性情報を抽出し、前記各小班区画属性情報に関する小班区画属性情報画像データを作成する小班区画属性情報画像データ作成部と、前記小班区画属性情報画像データと当該森林の地利条件とに基づいて間伐対象区域を選定し、選定された間伐対象区域を表す間伐対象区域画像データを作成する間伐対象区域画像データを作成部とを有することを特徴とする。   (11) A thinning target area selection system according to the present invention is a thinning target area selection system for selecting a thinning target area to be thinned based on forest present state information image data relating to a current state of a forest to be investigated, Based on the forest present state information image data, a subcompartment image data for dividing the forest into subcompartments is created, and based on the subcompartment image data and the forest present state information image data, the forest present state subcompartment image data is obtained. Forest current subcompartment image data creation section to be created, and subcompartment attribute information image data for extracting subcompartment attribute information representing the characteristics of each subcompartment and creating subcompartment attribute information image data related to each subcompartment attribute information A thinning target area is selected based on the creation unit, the subcompartment section attribute information image data, and the land condition of the forest, and the selected thinning pair And having a preparation unit thinning target area image data to create a thinning target area image data representing the area.

本発明の間伐対象区域の選定システムにおいても前記間伐対象区域の選定方法と同様の効果を得ることができる。なお、本発明の間伐対象区域の選定システムは、当該間伐対象区域の選定システムに含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。また、本発明の間伐対象区域の選定システムにおいても、前記(9)及び(10)に記載の特徴を有することが好ましい。   The thinning target area selection system of the present invention can achieve the same effects as the thinning target area selection method. In the thinning target area selection system of the present invention, the functions of each component included in the thinning target area selection system are installed as a computer program, and predetermined data is given to each component. Thus, the functions of each component are executed on the software of the computer. Moreover, also in the selection system of the thinning object area of this invention, it is preferable to have the characteristics as described in said (9) and (10).

以下、本発明の実施形態について説明する。   Hereinafter, embodiments of the present invention will be described.

[実施形態1]
図1は、実施形態1に係る樹種分類システムの構成を示す図である。実施形態1に係る樹種分類システム100は、図1に示すように、撮影画像データ入力部110と、森林域画像データ作成部120と、樹冠画像データ作成部130と、樹冠情報画像データ作成部140と、樹種分類画像データ作成部150とを有している。
[Embodiment 1]
FIG. 1 is a diagram illustrating a configuration of a tree species classification system according to the first embodiment. As shown in FIG. 1, the tree classification system 100 according to the first embodiment includes a captured image data input unit 110, a forest area image data creation unit 120, a canopy image data creation unit 130, and a canopy information image data creation unit 140. And a tree species classification image data creation unit 150.

撮影画像データ入力部110は、調査対象となる森林を含む地域(調査対象地域という。)を上空から撮影して得られた撮影画像データを入力する機能を有している。
森林域画像データ作成部120は、撮影画像データと地理情報システム(GISという。)500に格納されている地理的情報(この場合、都道府県などで整備されている森林の管理台帳などに登録されている森林域の境界データ)とに基づいて森林域を抽出し、抽出された森林域に関する森林域画像データを作成する機能を有している。
The photographed image data input unit 110 has a function of inputting photographed image data obtained by photographing a region including a forest to be surveyed (referred to as a survey target region) from above.
The forest area image data creation unit 120 is registered in the captured image data and the geographic information stored in the geographic information system (GIS) 500 (in this case, the forest management ledger or the like maintained in the prefecture or the like). The forest area is extracted based on the boundary data of the forest area, and the forest area image data relating to the extracted forest area is created.

樹冠画像データ作成部130は、森林域画像データに基づいて森林域に含まれる各樹冠を抽出し、抽出された各樹冠に関する樹冠画像データを作成する機能を有している。
樹冠情報画像データ作成部140は、樹冠画像データに基づいて各樹冠の日向部、日陰部及び樹頂部などの情報を樹冠情報として抽出し、抽出された樹冠情報に関する樹冠情報画像データを作成する機能を有している。
樹種分類画像データ作成部150は、樹冠情報画像データに基づいて各樹冠の日向部とされた領域における分光反射特性を作成し、作成された分光反射特性に基づいて当該樹冠の樹種を特定して樹種分類画像データを作成する機能を有している。
The canopy image data creation unit 130 has a function of extracting each canopy included in the forest area based on the forest area image data and creating canopy image data regarding each extracted canopy.
The canopy information image data creation unit 140 is a function for extracting information on the crown, shade, and top of each canopy as canopy information based on the canopy image data, and creating canopy information image data regarding the extracted canopy information. have.
The tree species classification image data creation unit 150 creates a spectral reflection characteristic in the area defined as the sunny part of each crown based on the canopy information image data, and identifies a tree type of the crown based on the created spectral reflection characteristic It has a function to create tree species classification image data.

なお、この明細書において、「樹冠」とは森林を構成する各樹木を上部から見たときに個々の樹木の枝と葉で構成された部分をいう。また、上記撮影画像データは、高解像度の撮影機器によって撮影された高解像度画像データであって、撮影画像データの1画素が数10cm×数10cm〜100cm×100cm程度の面積に対応しているものとする。例えば、1つの樹冠の面積が100平方メートルであって、撮影画像データの1画素が仮に100cm×100cmの面積に対応するとすれば、1つの樹冠に対する撮影画像データの画素数は100画素となる。   In this specification, “tree canopy” refers to a portion composed of branches and leaves of individual trees when each tree constituting the forest is viewed from above. The captured image data is high-resolution image data captured by a high-resolution imaging device, and one pixel of the captured image data corresponds to an area of about several tens cm × several tens cm to 100 cm × 100 cm. And For example, if the area of one tree crown is 100 square meters and one pixel of the photographed image data corresponds to an area of 100 cm × 100 cm, the number of pixels of the photographed image data for one tree crown is 100 pixels.

また、図1に示した樹種分類システムは、当該樹種分類システムに含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   Further, in the tree species classification system shown in FIG. 1, the functions of each component included in the tree species classification system are installed as a computer program, and by giving predetermined data to each component, The functions of these components are executed on the software of the computer.

図2は、実施形態1に係る樹種分類システムにおける樹種分類処理を説明するフローチャートである。樹種分類処理は、図2に示すように、まずは、調査対象地域Aを撮影して得られた撮影画像データから得られる調査対象地域Aに地理情報システム500に記憶されている森林域の境界データを重ね合わせた森林域境界画像データを作成する(ステップS110)。続いて、当該森林域境界画像データに基づいて森林域と非森林域とを区別するための森林域画像データを作成する(ステップS120)。   FIG. 2 is a flowchart for explaining tree species classification processing in the tree species classification system according to the first embodiment. As shown in FIG. 2, the tree species classification process is first performed by forest area boundary data stored in the geographic information system 500 in the survey target area A obtained from the photographed image data obtained by photographing the survey target area A. Forest area boundary image data is created by superimposing (step S110). Subsequently, forest area image data for distinguishing between a forest area and a non-forest area is created based on the forest area boundary image data (step S120).

ステップS110,S120の処理は、森林域画像データ作成部120によって行われ、森林域画像データ作成部120によって森林域境界画像データが作成されることにより、当該森林域境界画像データに対応する森林域境界画像(図3参照。)をディスプレイ上に表示することができ、また、森林域画像データが作成されることにより、当該森林域画像データに対応する森林域画像(図4参照。)をディスプレイ上に表示させることができる。   The processes in steps S110 and S120 are performed by the forest area image data creation unit 120, and the forest area boundary image data is created by the forest area image data creation unit 120, so that the forest area corresponding to the forest area boundary image data is generated. The boundary image (see FIG. 3) can be displayed on the display, and the forest area image data (see FIG. 4) corresponding to the forest area image data is displayed by creating the forest area image data. Can be displayed above.

図3は、ディスプレイ上に表示された森林域境界画像を示す図である。なお、図3は調査対象地域Aを撮影して得られた撮影画像であり、当該調査対象地域Aに森林域の境界を重ね合わせた状態を示している。図3において、白線L1が各々の森林域の境界線であり、白線L1で囲まれた各領域(符号1,2,3,・・・が付された領域)が各々の森林域である。なお、図3においては、調査対象地域A内に存在する森林域のうちの一部の森林域のみに森林域の境界線(白線L1)が重ね合わされた例が示されている。なお、図3はモノクロ画像であるため、図3からは、森林域の境界を識別しにくいが、実際には図3の元となるカラー画像上では森林域の境界を容易に読み取ることができる。   FIG. 3 is a diagram showing a forest area boundary image displayed on the display. FIG. 3 is a photographed image obtained by photographing the survey target area A, and shows a state where the boundary of the forest area is superimposed on the survey target area A. In FIG. 3, the white line L1 is a boundary line of each forest area, and each area (areas denoted by reference numerals 1, 2, 3,...) Surrounded by the white line L1 is each forest area. Note that FIG. 3 shows an example in which a forest area boundary line (white line L1) is superimposed only on a part of the forest area in the survey area A. Since FIG. 3 is a monochrome image, it is difficult to identify the boundary of the forest area from FIG. 3, but in practice, the boundary of the forest area can be easily read on the original color image of FIG. .

図4は、ディスプレイ上に表示された森林域画像を示す図である。なお、図4において、調査対象地域A内において灰色で示されている領域が森林域Bであり、当該調査対象地域A内において白抜きで示されている領域が非森林域Cである。なお、図4はモノクロ画像であるため、図4からは、森林域Bと非森林域Cとが色別として読み取ることは困難であるが、実際には図4の元となるカラー画像上では森林域Bと非森林域Cとを色別で読み取ることができる。この場合、図4の元となるカラー画像上では森林域Bが赤色及び紫色で示され、非森林域Cが緑色で示されているため、ディスプレイ上において表示されるカラー画像上では抽出された森林域を容易に読み取ることができる。   FIG. 4 is a diagram showing a forest area image displayed on the display. In FIG. 4, the area shown in gray in the survey target area A is the forest area B, and the area shown in white in the survey target area A is the non-forest area C. Since FIG. 4 is a monochrome image, it is difficult to read the forest area B and the non-forest area C as different colors from FIG. 4, but in actuality, on the color image that is the original of FIG. The forest area B and the non-forest area C can be read by color. In this case, since the forest area B is shown in red and purple on the original color image of FIG. 4 and the non-forest area C is shown in green, it was extracted on the color image displayed on the display. The forest area can be easily read.

上記のような森林域の抽出は公知の技術であり、例えば、下記の公知文献1及び公知文献2が知られている。   Extraction of the forest area as described above is a known technique. For example, the following known document 1 and known document 2 are known.

公知文献1:加藤正人、“衛星リモートセンシングと地図の重ね合わせによる森林現況の把握”、日本林学会北海道支部論文集、1991年2月、第39号(別冊)、p,117〜118   Known Document 1: Masato Kato, “Understanding of Forest Current Status by Satellite Remote Sensing and Map Overlay”, The Japanese Forestry Society Hokkaido Chapter, 1991, No. 39 (separate volume), p, 117-118

公知文献2:加藤正人、“衛星リモートセンンシンングとGISの組み合わせ技術によるトドマツ小班データベースへの適用”、日本林学会誌別冊、1993年3月、Vol.75(No.2)、p,154〜158   Known Document 2: Masato Kato, “Application to Todomatsu Subcombination Database by Combination Technology of Satellite Remote Sensing and GIS”, Journal of the Japanese Forestry Society, March 1993, Vol. 75 (No. 2), p, 154 to 158

上記ステップS120によって森林域画像データが作成されると、森林域画像データに基づいて森林域に含まれる各樹冠を抽出し、抽出された各樹冠に関する樹冠画像データを作成する(ステップS130)。続いて、樹冠画像データに基づいて各樹冠ごとの日向部、日陰部及び樹頂部などの情報を樹冠情報として抽出し、樹冠情報画像データを作成する(ステップS140)。   When the forest area image data is created in step S120, each canopy included in the forest area is extracted based on the forest area image data, and canopy image data relating to each extracted canopy is created (step S130). Subsequently, information such as the sunshine part, the shade part, and the tree top part for each canopy is extracted as canopy information based on the canopy image data to create canopy information image data (step S140).

樹冠は太陽光の影響を受け、樹冠の日向部、日陰部及び樹頂部においてそれぞれ色調が異なる。したがって、樹冠画像データの輝度値によって樹冠の日向部、日陰部及び樹頂部に区分することができる。ここでは、輝度値が平均値以上である箇所を日向部とし、樹冠画像データにおける輝度値が平均値未満である箇所を日陰部とし、樹冠画像データにおける輝度値が極大値である箇所を樹頂部として各樹冠を細区分した樹冠情報画像データを作成する。   The canopy is affected by sunlight, and the shade of the sun is different from that of the sun, shade and top of the canopy. Accordingly, it is possible to classify the tree into the sun, shade and top of the tree according to the brightness value of the canopy image data. Here, a portion where the luminance value is equal to or higher than the average value is a sunburst portion, a portion where the luminance value in the crown image data is less than the average value is a shade portion, and a portion where the luminance value in the crown image data is a maximum value is the top portion As a result, crown information image data is created by subdividing each crown.

ステップS130の処理は樹冠画像データ作成部130によって行われ、樹冠画像データ作成部130によって樹冠画像データが作成されることにより、当該樹冠画像データに対応する樹冠画像(図5参照。)をディスプレイ上に表示させることができる。また、ステップS140の処理は樹冠情報画像データ作成部140によって行われ、樹冠情報画像データ作成部140によって樹冠情報画像データが作成されることにより、当該樹冠情報画像データに対応する樹冠情報画像(図6参照。)をディスプレイ上に表示させることができる。   The processing of step S130 is performed by the canopy image data creation unit 130, and the canopy image data is created by the canopy image data creation unit 130, whereby a canopy image (see FIG. 5) corresponding to the canopy image data is displayed on the display. Can be displayed. Further, the processing of step S140 is performed by the canopy information image data creation unit 140, and the canopy information image data is created by the canopy information image data creation unit 140, so that the canopy information image corresponding to the canopy information image data (FIG. 6) can be displayed on the display.

図5は、ディスプレイ上に表示された樹冠画像を示す図である。図5において、全体的に白で表示されている領域が図4に示す森林域Bに対応しており、白で表示されている領域をよく見ると、多数の白い粒が寄せ集まっているかのように見える。図5における1つ1つの白い粒が各樹冠を表している。図5ではモノクロ画像で示されているが、実際にはカラー画像としてディスプレイ上に表示される。図5によれば、樹冠域と非樹冠域とを区分することができる。なお、非樹冠域は上空から見た場合、上層樹冠の樹木がない箇所であり、陰部や高木の樹木がない下層植生に覆われている箇所(ギャップ)である。各樹冠域の画素数を積算することで、樹冠面積を求めることができる。   FIG. 5 is a diagram showing a tree crown image displayed on the display. In FIG. 5, the area displayed in white as a whole corresponds to the forest area B shown in FIG. 4, and if the area displayed in white is viewed closely, a large number of white grains are gathered together. looks like. Each white grain in FIG. 5 represents each tree crown. Although shown as a monochrome image in FIG. 5, it is actually displayed on the display as a color image. According to FIG. 5, the canopy area and the non-crown area can be divided. In addition, when viewed from above, the non-crown area is a place where there is no tree of the upper tree canopy, and is a place (gap) covered with lower vegetation where there is no genital area or tree of the high tree. By integrating the number of pixels in each canopy area, the canopy area can be obtained.

ところで、樹冠の抽出は、樹冠の縁部(樹冠縁という)が樹冠部分よりも暗いことを利用して行うものであり、これは、公知のブァレイフォローイング・アルゴリズム(Valley following algorithm)と自動マスク処理によって可能となる。また、各樹冠が抽出されることにより、抽出された各樹冠に対応する画素を積算することによって、各樹冠の面積を求めることができる。   By the way, the extraction of the canopy is performed by utilizing the fact that the edge of the canopy (called the canopy edge) is darker than the canopy, and this is based on the known valley following algorithm and automatic masking. It becomes possible by processing. Further, by extracting each canopy, the area of each canopy can be obtained by integrating the pixels corresponding to each extracted canopy.

図6は、ディスプレイ上に表示された樹冠情報画像を部分的に拡大して示す図である。樹冠情報画像は、図6に示すように、各樹冠において日向部、日陰部及び樹頂部によって色調が異なる。すなわち、図6において、ある1つの樹冠T1を例にとって説明すると、当該樹冠T1においては、日向部Su、日陰部Sh及び樹頂部Peが存在する。   FIG. 6 is a diagram showing a partially enlarged tree crown information image displayed on the display. As shown in FIG. 6, the tree crown information image has a different color tone depending on the sun, the shade, and the top of each tree. In other words, in FIG. 6, a certain crown T1 will be described as an example. In the crown T1, there are a sunshine portion Su, a shade portion Sh, and a treetop portion Pe.

なお、図6はモノクロ画像であるため、図6からは、日向部Su、日陰部Sh及び樹頂部Peを色別として読み取ることは困難であるが、実際には図6の元となるカラー画像上では日向部Su、日陰部Sh及び樹頂部Peを読み取ることは容易である。この場合、図6の元となるカラー画像上では日向部Suが桃色で示され、日陰部Shが黄色で示され、樹頂部Peが青色で示されているため、ディスプレイ上において表示されるカラー画像上では、日向部Su、日陰部Sh及び樹頂部Peを容易に読み取ることができる。   Since FIG. 6 is a monochrome image, it is difficult to read the sunshine portion Su, the shaded portion Sh, and the treetop portion Pe by color from FIG. 6, but actually, the color image that is the original of FIG. In the above, it is easy to read the sun portion Su, the shade portion Sh and the tree top portion Pe. In this case, in the original color image of FIG. 6, the sun portion Su is shown in pink, the shade portion Sh is shown in yellow, and the tree top portion Pe is shown in blue. Therefore, the color displayed on the display On the image, the sun portion Su, the shade portion Sh and the tree top portion Pe can be easily read.

次に、樹種分類画像データ作成部150により、樹冠の日向部の分光反射特性(スペクトル特性)を用いて各樹冠ごとに樹種を抽出し、抽出された樹種に基づいて樹種分類画像データを作成する(ステップS150)。   Next, the tree species classification image data creation unit 150 extracts the tree species for each canopy using the spectral reflection characteristics (spectral characteristics) of the sunny part of the crown, and creates tree species classification image data based on the extracted tree species. (Step S150).

なお、上空から調査対象地域Aを撮影して得られた撮影画像データの分光反射特性に関しては、以下の公知文献3に記載されている。ただし、公知文献3は、樹冠全体の輝度値を用いて分光反射特性を求めるものである。   Note that the spectral reflection characteristics of the captured image data obtained by imaging the survey area A from the sky are described in the following known document 3. However, in the known document 3, the spectral reflection characteristic is obtained using the luminance value of the entire tree crown.

公知文献3:Masato Katoh、“Classifying tree species in a northern mixed forest using high-resolution IKONOS data”、Jour. Forestry Research9、2004年、p,7〜14   Known Document 3: Masato Katoh, “Classifying tree species in a northern mixed forest using high-resolution IKONOS data”, Jour. Forestry Research 9, 2004, p. 7-14

図7は、樹冠の日向部の分光反射特性の一例を示す図である。図7において、横軸は光の波長を表し、縦軸は輝度値を表している。ここでは、ヒノキCo、アカマツPd、カラマツLk、スギCj、メタセコイアMg、広葉樹Blの6種類の樹木をサンプルとし、光の波長域として青(バンド1)、緑(バンド2)、赤(バンド3)、近赤外(バンド4)における輝度値を折れ線グラフで示している。   FIG. 7 is a diagram illustrating an example of the spectral reflection characteristics of the sunlight portion of the tree crown. In FIG. 7, the horizontal axis represents the wavelength of light, and the vertical axis represents the luminance value. Here, six kinds of trees of hinoki Co, red pine Pd, larch Lk, cedar Cj, metasequoia Mg, and broad-leaved tree Bl are used as samples, and light (waveband 1), green (band 2), and red (band 3) are used. ), Luminance values in the near infrared (band 4) are shown by line graphs.

図7において、輝度値は赤色(バンド3)で全ての樹種において低くなり、近赤外域(バンド4)で高くなる。また、樹種による輝度値の差は、近赤外域(バンド4)で、より大きくなる。このような特性は樹種の分光反射特性(スペクトル特性)と呼ばれており、この分光反射特性により樹種を特定することができる。   In FIG. 7, the luminance value is red (band 3) and decreases in all tree species, and increases in the near-infrared region (band 4). In addition, the difference in the luminance value depending on the tree species becomes larger in the near infrared region (band 4). Such a characteristic is called a spectral reflection characteristic (spectral characteristic) of the tree species, and the tree species can be specified by this spectral reflection characteristic.

このように、日向部における分光反射特性を用いることにより、誤分類が生じにくく高い精度で樹種を特定することができる。ちなみに、樹冠全体の輝度値を用いて分光反射特性を求めると、輝度値の標準偏差が大きく、誤分類が生じやすいことが知られている。これは、太陽光の日陰部分である樹冠縁は高さが低く、隣接する樹木の枝葉や陰の影響を受けやすいので、純粋な樹種の輝度値を反映していないことが起因するものと考えられる。   In this way, by using the spectral reflection characteristics in the sunny area, it is possible to specify the tree species with high accuracy that is unlikely to cause misclassification. Incidentally, it is known that when the spectral reflection characteristic is obtained using the luminance value of the entire tree crown, the standard deviation of the luminance value is large and misclassification is likely to occur. This is thought to be due to the fact that the crown edge, which is the shaded part of sunlight, is low in height and is susceptible to the branches and shades of adjacent trees, so it does not reflect the brightness value of pure tree species. It is done.

なお、ステップS150の処理は、樹種分類画像データ作成部150によって行われ、樹種分類画像データ作成部150によって樹種分類画像データが作成されることにより、当該樹種分類画像データに対応する樹種分類画像(図8参照)をディスプレイ上に表示させることができる。   The processing of step S150 is performed by the tree species classification image data creation unit 150, and the tree species classification image data is created by the tree species classification image data creation unit 150, so that the tree species classification image (corresponding to the tree species classification image data ( 8) can be displayed on the display.

図8は、ディスプレイ上に表示された森林域全体の樹種分類画像を示す図である。なお、図8はモノクロ画像であるため、図8からは、各樹種を色別として読み取ることは困難であるが、実際には図8の元となるカラー画像上ではヒノキCo、アカマツPd、カラマツLk、スギCj、メタセコイアMg、広葉樹Blが色別で示されているので、これら各樹種を容易に読み取ることができる。この場合、図8の元となるカラー画像上では、ヒノキCoが赤色、アカマツPdが橙色、カラマツLkが黄緑色、スギCjが濃緑色、メタセコイアMgが黄色、広葉樹Blが灰色で示されている。   FIG. 8 is a diagram showing a tree species classification image of the entire forest area displayed on the display. Since FIG. 8 is a monochrome image, it is difficult to read each tree species by color from FIG. 8, but actually, on the color image that is the original of FIG. 8, cypress Co, red pine Pd, and larch. Since Lk, cedar Cj, metasequoia Mg, and broad-leaved tree Bl are shown by color, each of these tree species can be easily read. In this case, cypress Co is red, red pine Pd is orange, larch Lk is yellow green, cedar Cj is dark green, metasequoia Mg is yellow, and broadleaf tree Bl is gray on the original color image of FIG. .

なお、実施形態1に係る樹種分類画像データ作成部150においては、実際には、画素(前述したように、1つの樹冠について100画素程度存在する。)ごとに樹種分類処理を実行し、各樹冠において、各画素ごとの樹種分類を行い(図9参照。)、分類された樹種の面積を計算し、面積の多い樹種を当該樹冠における樹種とする(図10参照。)。   In addition, in the tree classification image data creation unit 150 according to the first embodiment, actually, tree classification processing is executed for each pixel (as described above, there are about 100 pixels for one tree crown), and each tree crown is executed. Then, the tree species classification for each pixel is performed (see FIG. 9), the area of the classified tree species is calculated, and the tree species having a large area is set as the tree species in the crown (see FIG. 10).

図9は、各樹冠の画素ごとに樹種分類を行った樹種分類結果を示す図である。なお、図9の樹種分類は、各樹冠において日向部を基準として樹種分類を行った結果である。各樹冠における画素ごとに樹種分類を行った樹種分類結果は、図9に示すように、各樹冠内において、複数の樹種が混在した状態となっている。例えば、図9において、ある1つの樹冠T2を例にとって説明すると、当該樹冠T2においては、当該樹冠T2に対応する画素のうちの多くの画素がカラマツ(図9における白色の領域)として分類されているが、一部の画素においてはスギ(図9における灰色の領域)として分類されている。このように、1つの樹冠内で複数の樹種が混在して分類されるのは、樹冠表面が枝の伸長差により一様でないため、小さな日陰部や葉量の違いがあること、林冠縁では隣接する樹冠同士において、枝や葉が入り混じる場合があることなどが要因であると考えられる。   FIG. 9 is a diagram illustrating a tree species classification result obtained by performing tree species classification for each pixel of each crown. Note that the tree species classification in FIG. 9 is the result of tree species classification with respect to the sunny area at each crown. As shown in FIG. 9, the tree species classification result obtained by classifying the tree species for each pixel in each canopy is in a state where a plurality of tree species are mixed in each canopy. For example, in FIG. 9, a certain crown T2 will be described as an example. In the crown T2, many of the pixels corresponding to the crown T2 are classified as larch (white area in FIG. 9). However, some of the pixels are classified as cedar (gray area in FIG. 9). In this way, multiple tree species are classified and mixed within one canopy because the surface of the canopy is not uniform due to the difference in branch elongation, there are small shades and differences in the amount of leaves, It is considered that there are cases where branches and leaves are mixed between adjacent tree crowns.

なお、図9はモノクロ画像であるため、図9からは、各樹種を色別として読み取ることは困難であるが、実際には図9の元となるカラー画像上ではカラマツLk、スギCjを容易に区別して読み取ることができる。この場合、図9の元となるカラー画像上では、カラマツLkが黄緑色、スギCjが濃緑色で示されている。   Since FIG. 9 is a monochrome image, it is difficult to read each tree species by color from FIG. 9, but in reality, larch Lk and cedar Cj are easily displayed on the original color image of FIG. Can be read separately. In this case, the larch Lk is shown in yellow green and the cedar Cj is shown in dark green on the original color image of FIG.

図10は、図9に示した画素ごとの樹種分類結果に基づいて1つの樹冠を1つ樹種として分類した樹種分類結果を示す図である。1つの樹冠を1つ樹種として分類する処理は、各樹冠において、まずは、図9に示すように、各画素ごとの樹種分類を行い、分類された樹種の面積を計算し、面積の多い樹種を当該樹冠における樹種とする。このような処理を行うことにより、各樹冠はそれぞれ1つの樹種として分類される。   FIG. 10 is a diagram illustrating a tree species classification result obtained by classifying one crown as one tree species based on the tree species classification result for each pixel illustrated in FIG. 9. The process of classifying one crown as one tree type is as follows. First, in each crown, as shown in FIG. 9, the tree type is classified for each pixel, the area of the classified tree type is calculated, and the tree type having a large area is selected. The tree species in the crown. By performing such processing, each tree crown is classified as one tree species.

なお、図10はモノクロ画像であるため、図10からは、各樹種を色別として読み取ることは困難であるが、実際には図10の元となるカラー画像上では樹種を容易に区別して読み取ることができる。この場合、図10の元となるカラー画像上では、ヒノキCoが濃緑色、カラマツLkが黄色、アカマツPdが赤色、スギCjが緑色で示されている。このように、面積の多い樹種を当該樹冠における樹種とするという処理を行うには、より正確な分光反射特性を用いて樹種分類を行う必要があり、これを実現するために、実施形態1に係る樹種分類システムにおいては、各樹冠の日向部とされた領域における分光反射特性によって樹種分類を行っている。   Since FIG. 10 is a monochrome image, it is difficult to read each tree type for each color from FIG. 10, but actually, the tree type is easily distinguished and read on the original color image of FIG. be able to. In this case, cypress Co is shown in dark green, larch Lk in yellow, red pine Pd in red, and cedar Cj in green on the original color image of FIG. As described above, in order to perform a process of setting a tree species having a large area as a tree species in the crown, it is necessary to classify the tree species using more accurate spectral reflection characteristics. In such a tree species classification system, tree species classification is performed based on spectral reflection characteristics in a region defined as a sunny part of each crown.

以上説明したように実施形態1に係る樹種分類システムによれば、各樹冠の日向部における分光反射特性を用いて樹種分類を行うことにより、誤分類が生じにくく高い精度で樹種を分類することができる。樹種が高い精度で分類されることによって、森林の植生の調査及び評価を行う場合、調査結果及び評価結果の信頼性を高めることができる。   As described above, according to the tree species classification system according to the first embodiment, by performing tree species classification using the spectral reflection characteristics in the sunlight portion of each crown, it is possible to classify the tree species with high accuracy that is unlikely to cause misclassification. it can. By classifying the tree species with high accuracy, the reliability of the survey results and the evaluation results can be improved when surveying and evaluating the forest vegetation.

[実施形態2]
実施形態2は、実施形態1に係る樹種分類システムにおいて得られた樹種分類画像データを用いて森林の現況を示す森林現況情報画像データを作成するものである。
[Embodiment 2]
In the second embodiment, forest status information image data indicating the current status of the forest is created using the tree species classification image data obtained in the tree species classification system according to the first embodiment.

図11は、実施形態2に係る森林現況情報の作成システムの構成を示す図である。実施形態2に係る森林現況情報の作成システム200は、実施形態1に係る樹種分類システム100に加えて、全樹種樹冠占有画像データ作成部210と、樹種別樹冠占有画像データ作成部220と、全樹種樹冠本数分布画像データ作成部230と、樹種別樹冠本数分布画像データ作成部240と、森林現況情報画像データ作成部250とを有している。   FIG. 11 is a diagram showing a configuration of a system for creating forest status information according to the second embodiment. In addition to the tree species classification system 100 according to the first embodiment, the forest present situation information creation system 200 according to the second embodiment includes an all-tree species crown-occupied image data creation unit 210, a tree-type crown-occupied image data creation unit 220, It has a tree species crown number distribution image data creation unit 230, a tree type crown number distribution image data creation unit 240, and a forest present state information image data creation unit 250.

全樹種樹冠占有画像データ作成部210は、樹種分類システム100において作成された樹種分類画像データに基づいて全樹種の樹冠本数を算出する機能と、単位面積当たりの全樹種の樹冠面積占有度(全樹種樹冠面積占有度という)を算出して、算出された単位面積当たりの全樹種樹冠面積占有度に基づく全樹種樹冠占有画像データを作成する機能とを有している。   The all-tree species crown occupancy image data creation unit 210 has a function of calculating the number of crowns of all tree species based on the tree species classification image data created by the tree species classification system 100, and the crown area occupancy degree of all tree species per unit area (total A tree species canopy occupancy image data based on the calculated total tree species crown area occupancy per unit area.

樹種別樹冠占有画像データ作成部220は、樹種分類システム100において作成された樹種分類画像データに基づいて樹種別の樹冠本数を算出する機能と、単位面積当たりの樹種別の面積占有度(樹種別樹冠面積占有度という)を算出して、算出された単位面積あたりの樹種別樹冠面積占有度に基づく樹種別樹冠占有画像データを作成する機能とを有している。   The tree-type crown occupancy image data creation unit 220 has a function for calculating the number of tree-type crowns based on the tree-type classification image data created by the tree-type classification system 100, and the area occupancy of the tree type per unit area (tree type And a function of creating tree-type crown occupancy image data based on the calculated tree-type crown area occupancy per unit area.

全樹種樹冠本数分布画像データ作成部230は、全樹種の樹冠本数に基づいて単位面積当たりの全樹種樹冠本数分布画像データを作成する機能を有している。
樹種別樹冠本数分布画像データ作成部240は、樹種別の樹冠本数に基づいて単位面積当たりの樹種別樹冠本数分布画像データを作成する機能を有している。
The all tree species crown number distribution image data creation unit 230 has a function of creating all tree species crown number distribution image data per unit area based on the number of crowns of all tree species.
The tree-type crown number distribution image data creation unit 240 has a function of creating tree-type crown number distribution image data per unit area based on the number of tree-type crowns.

森林現況情報画像データ作成部250は、樹種分類画像データ、全樹種樹冠占有画像データ、樹種別樹冠占有画像データ、全樹種樹冠本数分布画像データ及び樹種別樹冠本数分布画像データに基づいて森林現況情報画像データを作成する機能を有している。   The forest present state information image data creation unit 250 is configured to provide forest present state information based on tree species classification image data, all tree species canopy occupation image data, tree type canopy occupation image data, all tree species canopy number distribution image data, and tree type canopy number distribution image data. It has a function to create image data.

なお、実施形態2に係る森林現況情報の作成システムは、当該森林現況情報の作成システムに含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   In addition, in the forest present status information creation system according to the second embodiment, the functions of each component included in the forest present status information creation system are installed as a computer program, and predetermined data is stored in each component. By giving the above, the function of each component is executed on the software of the computer.

図12は、実施形態2に係る森林現況情報の作成システムにおける森林現況情報作成処理を説明するフローチャートである。森林現況情報画像データ作成処理は、図12に示すように、樹種分類システム100において作成された樹種分類画像データに基づいて全樹種の樹冠本数を算出する(ステップS210)。また、樹種分類画像データに基づいて単位面積当たりの全樹種樹冠面積占有度を算出し(ステップS220)、算出された全樹種樹冠面積占有度に基づく全樹種樹冠占有画像データを作成する(ステップS230)。   FIG. 12 is a flowchart for explaining a forest present state information creation process in the forest present state information creation system according to the second embodiment. In the forest status information image data creation process, as shown in FIG. 12, the number of crowns of all tree species is calculated based on the tree species classification image data created by the tree species classification system 100 (step S210). Further, the total tree species crown area occupancy per unit area is calculated based on the tree species classification image data (step S220), and all tree species crown occupancy image data based on the calculated all tree species crown area occupancy is created (step S230). ).

ステップS210〜S230の処理は、全樹種樹冠占有画像データ作成部210によって行われ、全樹種樹冠占有画像データ作成部210によって全樹種樹冠占有画像データが作成されることにより、当該全樹種樹冠占有画像データに対応する全樹種樹冠占有画像(図13参照。)をディスプレイ上に表示させることができる。   The processing of steps S210 to S230 is performed by the all-tree species crown occupancy image data creation unit 210, and the all-tree species crown occupancy image data creation unit 210 creates all-tree species crown occupancy image data, thereby the all-tree species crown occupancy image. An all-tree species crown occupation image (see FIG. 13) corresponding to the data can be displayed on the display.

図13は、ディスプレイ上に表示された全樹種樹冠占有画像を示す図である。図13の調査対象地域Aは図4と同じで地域ある。図13において白色で示す領域が図4に示す森林域における全樹種の樹冠占有領域を示しており、図4に示す森林域の全体に全樹種の樹冠が存在していることがわかる。   FIG. 13 is a diagram illustrating an all-tree species crown occupation image displayed on the display. 13 is the same area as FIG. In FIG. 13, the white area indicates the crown occupation area of all tree species in the forest area shown in FIG. 4, and it can be seen that the crowns of all tree species exist in the entire forest area shown in FIG. 4.

次に、樹種分類システム100において作成された樹種分類画像データに基づいて樹種別の樹冠本数を算出する(ステップS240)。また、樹種分類画像データに基づいて単位面積当たりの樹種別樹冠面積占有度を算出し(ステップS250)、算出された樹種別樹冠面積占有度に基づく樹種別樹冠占有画像データを作成する(ステップS260)。   Next, the number of tree crowns is calculated based on the tree classification image data created by the tree classification system 100 (step S240). Further, a tree type crown area occupancy per unit area is calculated based on the tree type classification image data (step S250), and tree type crown occupancy image data based on the calculated tree type crown area occupancy is created (step S260). ).

ステップS240〜S260の処理は、樹種別樹冠占有画像データ作成部220によって行われ、樹種別樹冠占有画像データ作成部220によって樹種別樹冠占有画像データが作成されることにより、当該樹種別樹冠占有画像データに対応する樹種別樹冠占有画像(図14参照。)をディスプレイ上に表示させることができる。   The processing in steps S240 to S260 is performed by the tree-type crown occupancy image data creation unit 220, and the tree-type crown occupancy image data creation unit 220 creates the tree-type crown occupancy image data. A tree-type crown occupancy image (see FIG. 14) corresponding to the data can be displayed on the display.

図14は、ディスプレイ上に表示された樹種別樹冠占有画像を示す図である。図14の調査対象地域Aは図4と同じで地域ある。図14は樹種としてヒノキの樹冠面積占有度を示しており、図10において白色で示す領域が図4に示す森林域におけるヒノキの樹冠占有領域を示している。   FIG. 14 is a diagram illustrating a tree-type crown occupation image displayed on the display. 14 is the same area as FIG. FIG. 14 shows the occupancy degree of cypress as a tree species, and the area shown in white in FIG. 10 shows the occupancy area of cypress in the forest area shown in FIG.

次に、全樹種の樹冠本数から単位面積当たりの全樹種樹冠本数分布画像データを作成する(ステップS270)。ステップS27の処理は、全樹種樹冠本数分布画像データ作成部230によって行われ、全樹種樹冠本数分布画像データ作成部230によって全樹種樹冠本数分布画像データが作成されることにより、当該全樹種樹冠本数分布画像データに対応する全樹種樹冠本数分布画像(図15参照。)をディスプレイ上に表示させることができる。   Next, all-tree species crown number distribution image data per unit area is created from the number of crowns of all tree species (step S270). The processing of step S27 is performed by the all-tree species crown number distribution image data creation unit 230, and the all-tree species crown number distribution image data creation unit 230 creates all-tree species crown number distribution image data, whereby the number of all-tree species crowns is calculated. An all-tree species crown distribution image (see FIG. 15) corresponding to the distribution image data can be displayed on the display.

図15は、ディスプレイ上に表示された全樹種樹冠本数分布画像を示す図である。図15の調査対象地域Aは図4と同じで地域ある。図15において白色で示す領域が図4に示す森林域における全樹種の樹冠本数分布を示しており、図4に示す森林域の全体に全樹種の樹冠が高密度で存在することがわかる。   FIG. 15 is a diagram showing an image of the distribution of the number of crowns of all tree species displayed on the display. 15 is the same area as FIG. 15 indicates the crown number distribution of all tree species in the forest area shown in FIG. 4, and it can be seen that the crowns of all tree species exist at high density in the entire forest area shown in FIG.

次に、樹種別の樹冠本数から単位面積当たりの樹種別樹冠本数分布画像データを作成する(ステップS280)。ステップS280の処理は、樹種別樹冠本数分布画像データ作成部240によって行われ、樹種別樹冠本数分布画像データ作成部240によって樹種別樹冠本数分布画像データが作成されることにより、当該樹種別樹冠本数分布画像データに対応する樹種別樹冠本数分布画像(図16参照。)をディスプレイ上に表示させることができる。   Next, tree-type crown number distribution image data per unit area is created from the number of tree-type crowns (step S280). The processing in step S280 is performed by the tree-type crown number distribution image data creation unit 240, and the tree-type crown number distribution image data creation unit 240 creates the tree-type crown number distribution image data. A tree-type crown number distribution image (see FIG. 16) corresponding to the distribution image data can be displayed on the display.

図16は、ディスプレイ上に表示された樹種別樹冠本数分布画像を示す図である。図16の調査対象地域Aは図4と同じで地域ある。図16は樹種としてヒノキの樹冠本数分布を示しており、ヒノキは満遍なく存在しているが、特に、森林域B内の所定の領域B1に高密度で存在していることがわかる。   FIG. 16 is a diagram illustrating a tree-type crown number distribution image displayed on the display. The survey target area A in FIG. 16 is the same as FIG. FIG. 16 shows the crown distribution of cypress trees as tree species, and cypress is present evenly, but it can be seen that it exists in a predetermined region B1 in the forest area B at a high density.

次に、実施形態1に係る樹種分類システム100において得られた樹種分類画像データ、実施形態2に係る森林現況情報の作成システム200において得られた全樹種樹冠占有画像データ、樹種別樹冠占有画像データ、全樹種樹冠本数分布画像データ及び樹種別樹冠本数分布画像データに基づいて森林の現況区分情報を表す森林現況情報画像データを作成する(ステップS290)。ステップS290の処理は、森林現況情報画像データ作成部250によって行われ、森林現況情報画像データ作成部250によって森林現況情報画像データが作成されることにより、当該林現況情報画像データに対応する森林現況情報画像(図17参照。)をディスプレイ上に表示させることができる。   Next, tree species classification image data obtained by the tree species classification system 100 according to the first embodiment, all tree species crown occupancy image data, tree type crown occupancy image data obtained by the forest status information creation system 200 according to the second embodiment. Forest current state information image data representing the current state division information of the forest is created based on the total tree species tree crown number distribution image data and the tree type tree crown number distribution image data (step S290). The processing in step S290 is performed by the forest present state information image data creating unit 250, and the forest present state information image data is created by the forest present state information image data creating unit 250, whereby the forest present state information image data corresponding to the forest present state information image data is created. An information image (see FIG. 17) can be displayed on the display.

なお、森林現況情報画像データを作成する際、レーザープロファイラや写真判読で求めた同じ分解能の樹高画像データを用い、全樹種樹冠面積占有度画像データ、樹種別樹冠面積占有度画像データ、全樹種樹冠本数分布画像データ及び樹種別樹冠本数分布画像データに加えて樹高画像データから樹高の違いを用いるようにしてもよい。   When creating forest status information image data, use tree height image data of the same resolution obtained by laser profiler and photo interpretation, all tree species crown area area image data, tree type tree area area image data, all tree species crown The difference in tree height from the tree height image data may be used in addition to the number distribution image data and the tree type tree number distribution image data.

図17は、ディスプレイ上に表示された森林現況情報画像を示す図である。図17はモノクロ画像であるため、図17からは、樹種を色別として読み取ることは困難であるが、実際には図17の元となるカラー画像上ではヒノキCo、アカマツPd、カラマツLk、スギCjなどが色別で示されているので、これら各樹種を容易に読み取ることができる。この場合、図7の元となるカラー画像上では、ヒノキCoが赤色、アカマツPdが橙色、カラマツLkが黄緑色、スギCjが濃緑色というように樹種別に色分けされて示されている。   FIG. 17 is a diagram showing a forest present state information image displayed on the display. Since FIG. 17 is a monochrome image, it is difficult to read the tree species by color from FIG. 17, but actually, on the original color image of FIG. 17, cypress Co, red pine Pd, larch Lk, cedar Since Cj and the like are indicated by color, each of these tree species can be easily read. In this case, the cypress Co is red, the red pine Pd is orange, the larch Lk is yellow green, and the cedar Cj is dark green on the original color image in FIG.

なお、森林現況情報画像データを作成する際、ここでは、一定の面積より小さな面積を有する領域は隣接する大きな面積を有する領域に組み入れるような処理を行うものとし、図17は一定の面積より小さな面積を有する領域が、隣接する大きな面積を有する領域に組み入れられた結果を示している。図17に示す森林現況情報画像によれば、森林の現況を適切に把握することができる。   When creating the forest present state information image data, here, it is assumed that processing is performed such that a region having an area smaller than a certain area is incorporated into an adjacent region having a large area, and FIG. 17 is smaller than the certain area. The results show that the area having the area is incorporated into the adjacent area having a large area. According to the forest present state information image shown in FIG. 17, it is possible to appropriately grasp the present state of the forest.

また、図12に示すフローチャートの処理の順序は、図12に示す順序に限られるものではない。例えば、図12においては、全樹種樹冠占有画像データの作成(ステップS230)と樹種別樹冠占有画像データの作成(ステップS260)とを行ったのちに、全樹種樹冠本数分布画像データの作成(ステップS270)と樹種別樹冠本数分布画像データの作成(ステップS280)とを行うようにしたが、これに限られるものではなく、全樹種樹冠画像データの作成終了後に全樹種樹冠本数分布画像データの作成を行い、その後、樹種別樹冠画像データを作成し、樹種別樹冠本数分布画像データの作成を行うというような順序としてもよい。また、全樹種の樹冠本数の算出(ステップS210)は、全樹種樹冠本数分布画像データの作成(ステップS270)の直前であってもよく、同様に、樹種別の樹冠本数の算出(ステップS240)は、樹種別樹冠本数分布画像データの作成(ステップS280)の直前であってもよい。   Further, the order of processing in the flowchart shown in FIG. 12 is not limited to the order shown in FIG. For example, in FIG. 12, after the creation of all-tree species crown occupation image data (step S230) and the creation of tree-type crown occupation image data (step S260), creation of all-tree species crown number distribution image data (step S270) and creation of tree-type crown number distribution image data (step S280), but not limited to this, creation of all-tree species crown number distribution image data after the creation of all-tree species crown image data After that, the tree-type tree crown image data is created, and the tree-type tree crown number distribution image data is created. Further, the calculation of the number of crowns of all tree species (step S210) may be performed immediately before the generation of all-tree species crown number distribution image data (step S270), and similarly, the calculation of the number of tree types of crowns (step S240). May be immediately before the creation of tree-type crown number distribution image data (step S280).

[実施形態3]
実施形態3は、実施形態2に係る森林現況情報の作成システム200において作成された森林現況情報画像データに基づいて、間伐すべき区域(間伐対象区域)を選定するものである。なお、森林現況画像情報データは、調査対象地域における一般的な森林現況情報画像データを用いることが可能であるが、実施形態3においては、森林現況情報画像データとして、実施形態2に係る森林現況情報の作成システムにおいて作成された森林現況情報画像データを用いるものとする。
[Embodiment 3]
The third embodiment selects an area (thinning target area) to be thinned based on the forest present situation information image data created by the forest present situation information creation system 200 according to the second embodiment. In addition, although the forest present condition image information data can use the general forest present condition information image data in an investigation object area, in Embodiment 3, the forest present condition information which concerns on Embodiment 2 as forest present condition information image data is used. It is assumed that the forest present state image data created in the information creation system is used.

図18は、実施形態3に係る間伐対象区域の選定システムの構成を示す図である。実施形態2に係る間伐対象区域の選定システム300は、実施形態2に係る森林現況情報の作成システム200に加えて、森林現況小班区画画像データ作成部310と、小班区画属性情報画像データ作成部320と、間伐対象区域画像データ作成部330とを有する。   FIG. 18 is a diagram illustrating a configuration of a thinning target area selection system according to the third embodiment. The thinning target area selection system 300 according to the second embodiment includes a forest present state subcompartment image data creation unit 310 and a subcompartment segment attribute information image data creation unit 320 in addition to the forest present state information creation system 200 according to the second embodiment. And a thinning target area image data creation unit 330.

森林現況小班区画画像データ作成部310は、森林現況情報の作成システム200において作成された森林現況情報画像データに基づいて、森林を閉多角形(ポリゴン)状の小班区画に分別するための小班区画画像データを作成し、小班区画画像データを森林現況情報画像(図17参照。)に重ね合わせた森林現況小班区画画像データを作成する機能を有する。   The forest present state subcompartment image data creation unit 310 divides the forest into closed polygon (polygon) subcompartments based on the forest present state information image data created by the forest present state information creation system 200. It has a function of creating image data and creating forest present state subcompartment image data in which the subcompartment block image data is superimposed on the forest present state information image (see FIG. 17).

小班区画属性情報画像データ作成部320は、森林現況小班区画画像データを樹種分類画像データ(図8参照。)に重ね合わせることによって各小班区画の属性情報(各樹冠の樹冠本数、単位面積当たりの樹冠本数、樹冠半径、樹冠面積など各樹冠の特徴を現す情報)を抽出し、抽出された各小班区画の属性情報に関する小班区画属性情報画像データを作成する機能を有する。   The sub-compartment attribute information image data creation unit 320 superimposes the current forest sub-compartment image data on the tree species classification image data (see FIG. 8), thereby obtaining attribute information (number of crowns of each crown, per unit area) of each sub-compartment. (Information representing the characteristics of each crown, such as the number of tree crowns, crown radius, and crown area), and a function of creating subgroup attribute information image data relating to the attribute information of each extracted subsection.

間伐対象区域画像データ作成部330は、小班区画属性情報画像データとGIS500に記憶されている地理的情報(この場合、当該森林域における林道、等高線、立地因子などの情報)とに基づいて間伐対象区域を選定し、選定された間伐対象区域を表す間伐対象区域画像データを作成する機能を有する。   The thinning target area image data creation unit 330 is a thinning target based on the subcompartment attribute information image data and the geographical information stored in the GIS 500 (in this case, information on forest roads, contour lines, location factors, etc. in the forest area). It has a function of selecting an area and creating thinning target area image data representing the selected thinning target area.

なお、実施形態3に係る間伐対象区域の選定システムは、当該間伐対象区域の選定システムに含まれる上記各構成要素が有する機能がコンピュータのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピュータのソフトウエア上で実行されるものである。   Note that the thinning target area selection system according to the third embodiment is configured such that the functions of each component included in the thinning target area selection system are installed as a computer program, and predetermined data is stored in each component. By giving the above, the function of each component is executed on the software of the computer.

図19は、実施形態3に係る間伐対象区域の選定システムにおける間伐対象区域選定処理を説明するフローチャートである。間伐対象区域選定処理は、図19に示すように、まずは、実施形態2に係る森林現況情報の作成システム200によって作成された森林現況情報画像(図17参照)に基づいて、森林現況小班区画画像データを作成する(ステップS310)。この森林現況小班区画画像データは、森林域を閉多角形(ポリゴン)状の小班区画に区分するための小班区画輪郭画像を、森林現況情報画像(図17参照)に重ね合わせることによって作成することができる。ステップS310の処理は、森林現況小班区画画像データ作成部310によって行われ、森林現況小班区画画像データ作成部310によって森林現況小班区画画像データが作成さることにより、当該森林現況小班区画画像データに対応する森林現況小班区画画像(図20参照。)をディスプレイ上に表示させることができる。   FIG. 19 is a flowchart for explaining a thinning target area selection process in the thinning target area selection system according to the third embodiment. In the thinning target area selection processing, as shown in FIG. 19, first, based on the forest present state information image (see FIG. 17) created by the forest present state information creating system 200 according to the second embodiment, the forest present state subcompartment section image Data is created (step S310). The forest present state subcompartment image data is created by superimposing the subcompartment contour image for dividing the forest area into closed polygon (polygonal) subcompartments on the forest present state information image (see FIG. 17). Can do. The processing in step S310 is performed by the forest present state subcompartment image data creation unit 310, and the forest present state subcompartment image data creation unit 310 creates the forest present state subcompartment image data, thereby corresponding to the forest present state subcompartment image data. The current forest subcompartment section image (see FIG. 20) can be displayed on the display.

図20は、ディスプレイ上に表示された森林現況小班区画画像を示す図である。なお、図20はモノクロ画像であるため、図20からは、閉多角形(ポリゴン)で形成される複数の小班区画を読み取ることは困難であるが、実際には図20の元となるカラー画像上では、森林現況の1つである樹種が色別で示されており(ヒノキCoが赤色、アカマツPdが橙色、カラマツLkが黄緑色、スギCjが濃緑色など)、これら各樹種別の領域を囲む小班区画輪郭画像が白線で示されているので、図20の元となるカラー画像上では、小班区画を容易に読み取ることができる。   FIG. 20 is a diagram showing an image of the current forest subdivision section displayed on the display. Since FIG. 20 is a monochrome image, it is difficult to read a plurality of segment sections formed of closed polygons (polygons) from FIG. 20, but in reality, the color image that is the original of FIG. Above, tree species that are one of the current forest conditions are shown by color (hinoki Co is red, red pine Pd is orange, larch Lk is yellow green, cedar Cj is dark green, etc.), and the areas of each tree type Since the outline image of the segment area surrounding the segment is indicated by a white line, the segment area can be easily read on the original color image of FIG.

次に、森林現況小班区画画像データを樹種分類画像データ(図8参照)に重ね合わせることによって、各小班区画の特徴を表す属性情報(樹冠本数、単位面積当たりの樹冠本数、樹冠半径、樹冠面積などの情報)を抽出し、抽出された各小班区画の属性情報に関する小班区画属性情報画像データを作成する(ステップS320)。ステップS320の処理は、小班区画属性情報画像データ作成部320によって行われ、小班区画属性情報画像データ作成部320によって小班区画属性情報画像データが作成さることにより、当該小班区画属性情報画像データに対応する小班区画属性情報画像(図21参照。)をディスプレイ上に表示させることができる。   Next, attribute information (number of crowns, number of crowns per unit area, crown radius, crown area) representing the characteristics of each of the small plots by superimposing the forest present subcompartment image data on the tree classification image data (see FIG. 8). Etc.), and subgroup section attribute information image data relating to the extracted attribute information of each subsection section is created (step S320). The processing in step S320 is performed by the subcompartment segment attribute information image data creation unit 320, and the subcompartment segment attribute information image data creation unit 320 creates the subcompartment segment attribute information image data, thereby corresponding to the subcompartment segment attribute information image data. The segment segment attribute information image (see FIG. 21) to be displayed can be displayed on the display.

図21は、ディスプレイ上に表示された小班区画属性情報画像を示す図である。なお、図21はモノクロ画像であるため、図21からは、各小班区画輪郭画像及び各樹種を読み取ることは困難であるが、実際には図21の元となるカラー画像上では、各小班区画輪郭画像、各樹冠の樹種(ヒノキCo、アカマツPd、カラマツLk、スギCj、メタセコイアMg、広葉樹Bl)が色別で示されているので、各小班区画情報及び各樹種を容易に読み取ることができる。この場合、図21の元となるカラー画像上では、各小班区画は白線で示され、各樹冠の樹種としてのヒノキCoが赤色、アカマツPdが橙色、カラマツLkが黄緑色、スギCjが濃緑色、メタセコイアMgが黄色、広葉樹Blが灰色で示されている。   FIG. 21 is a diagram showing a segment segment attribute information image displayed on the display. Since FIG. 21 is a monochrome image, it is difficult to read the outline images and the tree species of each subcompartment from FIG. 21, but in actuality, each subcompartment section is displayed on the original color image of FIG. Since the contour image and tree species of each crown (Hinoki Co, Pinus Pd, Larch Lk, Sugi Cj, Metasequoia Mg, Broadleaf B1) are shown by color, each subcompartment information and each tree type can be easily read. . In this case, in the original color image of FIG. 21, each subdivision is indicated by a white line, cypress Co as the tree species of each crown is red, red pine Pd is orange, larch Lk is yellow green, and cedar Cj is dark green. Metasequoia Mg is shown in yellow and hardwood Bl is shown in gray.

図22は、図21における領域Dの拡大画像を示す図である。なお、図22において、領域Dに含まれる複数の小班区画のうち、ここでは、1〜5の番号(区画番号という。)が付された小班区画に注目して、区画番号1〜5が付された小班区画における全樹種の属性情報及び樹種別(ヒノキとする。)の属性情報について考える。なお、図22及び次に説明する図23においては、区画番号1〜5は、それぞれが丸で囲まれた数字で示されている。   FIG. 22 is a diagram showing an enlarged image of the region D in FIG. In FIG. 22, among the plurality of sub-compartments included in the region D, here, attention is paid to sub-compartments with numbers 1 to 5 (referred to as section numbers). Consider the attribute information of all tree species and the tree type (referred to as cypress) in the sub-compartment. In FIG. 22 and FIG. 23 to be described next, the partition numbers 1 to 5 are indicated by numbers surrounded by circles.

図23は、各小班区画から抽出された属性情報を示す図である。図23において、区画番号1の小班区画について見ると、区画番号1の小班区画は、単位面積(1ha)当たり全樹種の樹冠本数が多く、樹冠半径と樹冠面積が小さく、ヒノキの本数が多いことがわかる。区画番号2〜5の小班区画においても同様に、全樹種の属性情報及び樹種別の属性情報を抽出することができる(図23参照。)。   FIG. 23 is a diagram showing attribute information extracted from each subdivision block. In FIG. 23, when looking at the subdivision section of section number 1, the subcompartment section of section number 1 has a large number of crowns of all tree species per unit area (1 ha), a small crown radius and crown area, and a large number of cypress trees. I understand. Similarly, the attribute information of all tree species and the attribute information of the tree type can be extracted in the subdivision sections of the section numbers 2 to 5 (see FIG. 23).

なお、図23において、「total」は、樹冠本数及び樹冠本数/hについては、区画番号1〜5で求められた樹冠本数の合計及び樹冠本数/hの合計あり、樹冠直径及び樹冠面積については、区画番号1〜5で求められた樹冠直径の平均値及び樹冠面積の平均値である。 In FIG. 23, “total” is the sum of the number of crowns and the number of crowns / h, and is the sum of the number of crowns and the number of crowns / h obtained in section numbers 1 to 5, and the crown diameter and the crown area are as follows. The average value of the canopy diameter and the average value of the canopy area obtained from the division numbers 1 to 5.

次に、間伐対象区域としての評価を行うための地利条件画像データを作成する(ステップS330)。すなわち、間伐計画を立てるには、各小班区画の属性情報の他に、間伐作業のし易さ(伐採した樹木の搬出のし易さなども含む)を示す地利条件(林道からの距離、標高、傾斜など)を調べる必要がある。そこで、間伐対象区域としての評価を行うための地利条件画像データを作成する。ステップS330の処理は、間伐対象区域画像データ作成部330によって行われ、間伐対象区域画像データ作成部330によって地利条件画像データが作成されることにより、当該地利条件画像データに対応する地利条件画像(図24参照。)をディスプレイ上に表示させることができる。   Next, land condition image data for performing evaluation as a thinning target area is created (step S330). In other words, in order to make a thinning plan, in addition to the attribute information of each subcompartment, the land condition indicating the ease of thinning work (including the ease of carrying out felled trees) (distance from forest road, elevation) , Tilt etc.). Therefore, land condition image data for evaluation as a thinning target area is created. The processing in step S330 is performed by the thinning target area image data creation unit 330, and the thinning target area image data creation unit 330 creates the land condition condition image data, whereby the land condition condition image corresponding to the land condition condition image data ( FIG. 24) can be displayed on the display.

図24は、ディスプレイ上に表示された地利条件画像を示す図である。図24は地利条件の1つである林道からの距離を説明するための地利条件画像の一般例であり、本発明の各実施形態における調査対象地域A(図3参照)とは異なるものである。   FIG. 24 is a diagram showing a land condition image displayed on the display. FIG. 24 is a general example of a land condition image for explaining the distance from the forest road, which is one of the land conditions, and is different from the survey area A (see FIG. 3) in each embodiment of the present invention. .

なお、図24はモノクロ画像であるため、図24からは、林道及び当該林道からの距離を色別として読み取ることは困難であるが、実際には図19の元となるカラー画像上では、林道及び当該林道からの距離が色別で示されているので、林道及び当該林道からの距離を容易に読み取ることができる。この場合、図24の元となるカラー画像上では、林道は茶色で示され、当該林道からの距離dが「d≦50m」の領域は橙色で示され、当該林道からの距離dが「50m<d≦100m」の領域は黄色で示されている。   Note that since FIG. 24 is a monochrome image, it is difficult to read the forest road and the distance from the forest road by color from FIG. 24, but in reality, the forest road on the original color image of FIG. Since the distance from the forest road is indicated by color, the distance from the forest road and the forest road can be easily read. In this case, in the original color image of FIG. 24, the forest road is shown in brown, the region where the distance d from the forest road is “d ≦ 50 m” is shown in orange, and the distance d from the forest road is “50 m”. The region <d ≦ 100 m ”is shown in yellow.

次に、地理情報システム(GIS)500に格納されている地理情報(この場合、林道、等高線など)から、該当する小班区画の地利条件(林道からの距離、標高、傾斜など)を取得し、小班区画属性情報と併せて、間伐対象区域としての評価値の高さ(例えば、林道から近い、傾斜が緩い、樹冠面積が大きい、占有面積割合が高い、立木本数が多いなど)を数値モデル化した間伐評価モデル化画像データを作成する(ステップS340)。ステップS340の処理も、間伐対象区域画像データ作成部330によって行われ、間伐対象区域画像データ作成部330によって間伐評価モデル化画像データが作成されることにより、当該間伐評価モデル化画像データに対応する間伐評価モデル化画像(図25参照。)をディスプレイ上に表示させることができる。   Next, from the geographical information stored in the geographic information system (GIS) 500 (in this case, forest roads, contour lines, etc.), the land conditions (distance from the forest road, elevation, slope, etc.) Combined with subcompartment attribute information, numerical modeling of the height of the evaluation value as a thinning target area (for example, close to the forest road, gentle slope, large crown area, high occupied area ratio, large number of standing trees, etc.) The thinned evaluation modeled image data is created (step S340). The processing in step S340 is also performed by the thinning target area image data creation unit 330, and the thinning evaluation modeled image data is created by the thinning target area image data creation unit 330, thereby corresponding to the thinning evaluation modeled image data. A thinning evaluation modeled image (see FIG. 25) can be displayed on the display.

図25は、ディスプレイ上に表示された間伐評価モデル化画像の一例を示す図である。なお、図25は間伐評価モデル化画像の例について説明するための一般的な画像例であり、本発明の実施形態における調査対象地域A(図3参照)とは異なる区域の画像である。   FIG. 25 is a diagram illustrating an example of a thinning evaluation modeled image displayed on the display. FIG. 25 is a general image example for explaining an example of the thinning evaluation modeled image, and is an image of an area different from the survey target area A (see FIG. 3) in the embodiment of the present invention.

間伐評価モデル画像の一例として、図25に示すように、(a)標高画像、(b)傾斜画像、(c)斜面方位画像、(d)地利画像などの地利条件を示す立地因子画像と、(e)森林現況情報画像との照合からクロス集計を行い、数値モデル化することで、(f)間伐対象選定用画像を作成する。図25はモノクロ画像であるため、図25からは、(a)標高画像、(b)傾斜画像、(c)斜面方位画像、(d)地利画像、(e)森林現況情報画像及び(f)間伐対象選定用画像をそれぞれ識別することは困難であるが、実際には図25の元となるカラー画像上では、これら(a)標高画像、(b)傾斜画像、(c)斜面方位画像、(d)地利画像、(e)森林現況情報画像及び(f)間伐対象選定用画像は容易に識別することができる。   As an example of the thinning evaluation model image, as shown in FIG. 25, (a) an elevation image, (b) a tilt image, (c) a slope orientation image, (d) a location factor image indicating a land condition condition, and the like, (E) Cross-tabulation is performed from the comparison with the forest present state information image, and numerical modeling is performed, and (f) a thinning target selection image is created. Since FIG. 25 is a monochrome image, from FIG. 25, (a) an altitude image, (b) an inclined image, (c) a slope orientation image, (d) a land-use image, (e) a forest status information image, and (f) Although it is difficult to identify the thinning target selection images, in reality, on the color image that is the original of FIG. 25, these (a) elevation image, (b) tilt image, (c) slope orientation image, (D) A land image, (e) a forest status information image, and (f) a thinning target selection image can be easily identified.

なお、間伐評価モデル化画像データを作成する技術は、下記の公知文献4に記載されている。ただし、本発明の各実施形態に係る間伐対象区域の選定システムにおいては、上空から森林地区を撮影して得られた撮影画像データは、高解像度の撮影画像データであるが、下記の公知文献4は、高解像度ではない衛星データを用いている点が異なる。   In addition, the technique which produces thinning evaluation modeling image data is described in the following well-known document 4. However, in the thinning target area selection system according to each embodiment of the present invention, the photographed image data obtained by photographing the forest area from the sky is high-resolution photographed image data. Is different in that it uses non-high resolution satellite data.

公知文献4:加藤正人、“衛星リモートセンシングとGISを活用した林分の評価とモデル化”、写真測量とリモートセンシング、1991年、Vol.30,No.6   Known Document 4: Masato Kato, “Evaluation and Modeling of Forest Using Satellite Remote Sensing and GIS”, Photogrammetry and Remote Sensing, 1991, Vol. 30, no. 6

次に、間伐評価モデル画像データに基づいて間伐すべき区域(間伐対象区域)を選定するための間伐対象区域選定画像データを作成する(ステップS350)。ステップS350の処理も、間伐対象区域画像データ作成部330によって行われ、間伐対象区域画像データ作成部330によって間伐対象区域選定画像データが作成されることにより、当該間伐対象区域選定画像データに対応する間伐対象区域選定画像(図26参照。)をディスプレイ上に表示させることができる。   Next, thinning target area selection image data for selecting an area to be thinned (thinning target area) based on the thinning evaluation model image data is created (step S350). The processing of step S350 is also performed by the thinning target area image data creation unit 330, and the thinning target area image data creation unit 330 creates thinning target area selection image data, thereby corresponding to the thinning target area selection image data. A thinning target area selection image (see FIG. 26) can be displayed on the display.

図26は、ディスプレイ上に表示された間伐対象区域選定画像を示す図である。図26において、間伐対象区域として選定された小班区画の輪郭線を他の小班区画の輪郭線と区別して示している。なお、図26はモノクロ画像であるため、図26からは、間伐対象区域を示す小班区画の輪郭線を他の小班区画の輪郭線と区別しにくいが、実際には図26の元となるディスプレイ上のカラー画像上では、間伐対象区域を示す小班区画の輪郭線の色が赤色で示されているので、間伐対象区域を容易に知ることができる。   FIG. 26 is a diagram showing a thinning target area selection image displayed on the display. In FIG. 26, the outline of the small segment selected as the thinning target area is distinguished from the outline of the other small segments. Since FIG. 26 is a monochrome image, from FIG. 26, it is difficult to distinguish the outline of the subdivision section indicating the thinning target area from the outlines of the other subsection sections. In the upper color image, the outline of the subdivision section indicating the thinning target area is shown in red, so that the thinning target area can be easily known.

また、図26に示されている間伐対象区域は、図22における区画番号1の小班区画である。このように、区画番号1の小班区画が間伐対象区域として選定された理由は、間伐作業のし易さなど間伐対象区域としての評価値を高くする条件(例えば、林道から近い、傾斜が緩い、樹冠面積が大きい、占有面積割合が高い、立木本数が多いなど)を考慮した結果である。   Further, the thinning target area shown in FIG. 26 is the subdivision section of section number 1 in FIG. In this way, the reason why the subcompartment with the division number 1 is selected as the thinning target area is that the evaluation value as the thinning target area such as the ease of thinning is high (for example, close to the forest road, the slope is gentle, This is a result of taking into account the large crown area, high occupation area ratio, and the number of standing trees.

なお、図26に示すように間伐対象区域が選定されたとしても、実際に間伐すべき区域は、図26のように選定された間伐対象区域全体であるとは限らず、例えば、都道府県などで整備されている森林の管理台帳などの登録されている森林の境界(例えば、図3参照。)に基づいて、図26のように選定された間伐対象区域全体のうちの所定の区域内のみを実際の間伐対象区域とするという場合もある。   Note that even if a thinning target area is selected as shown in FIG. 26, the area to be actually thinned is not necessarily the entire thinning target area selected as shown in FIG. 26. For example, prefectures, etc. Based on the registered forest boundaries (see, for example, Fig. 3), such as the forest management ledger, etc., only within a predetermined area of the entire thinning target area selected as shown in Fig. 26 May be the actual thinning target area.

以上説明したように、実施形態3に係る間伐対象区域の選定システムによれば、従来においては、個人の労力と熟練に頼っていた間伐対象区域を、調査対象地域Aの中から客観的かつ適切に選定することできる。例えば、図4に示す調査対象地域Aについて間伐評価モデル画像を作成すれば、図4に示す調査対象地域A内においてどの区域を間伐対象区域として選定すべきかを客観的にかつ適切に選定することができる。   As described above, according to the thinning target area selection system according to the third embodiment, the thinning target area that has conventionally relied on the labor and skill of individuals is objective and appropriate from the survey target area A. Can be selected. For example, if a thinning evaluation model image is created for the survey target area A shown in FIG. 4, objectively and appropriately selecting which area should be selected as the thinning target area in the survey target area A shown in FIG. 4. Can do.

なお、本発明は上述の実施形態に限られるものではなく、本発明の要旨を逸脱しない範囲で種々変形実施可能となるものである。たとえば、実施形態3に係る間伐対象区域の選定システム300においては、森林現況情報画像データとして、実施形態2に係る森林現況作成システム200において作成された森林現況画像データを用いた場合を例示したが、森林現況画像データとしては、他の一般的な森林現況画像データを用いてもよい。   The present invention is not limited to the above-described embodiment, and various modifications can be made without departing from the spirit of the present invention. For example, in the thinning target area selection system 300 according to the third embodiment, the case where the forest present state image data created in the forest present state creation system 200 according to the second embodiment is used as the forest present state information image data is illustrated. As the forest present state image data, other general forest present state image data may be used.

実施形態1に係る樹種分類システムの構成を示す図である。It is a figure which shows the structure of the tree classification system which concerns on Embodiment 1. FIG. 実施形態1に係る樹種分類システムにおける樹種分類処理を説明するフローチャートである。It is a flowchart explaining the tree classification process in the tree classification system concerning Embodiment 1. FIG. ディスプレイ上に表示された森林域境界画像を示す図地である。It is a figure ground which shows the forest area boundary image displayed on the display. ディスプレイ上に表示された森林域画像を示す図である。It is a figure which shows the forest area image displayed on the display. ディスプレイ上に表示された樹冠画像を示す図である。It is a figure which shows the tree crown image displayed on the display. ディスプレイ上に表示された樹冠情報画像を部分的に拡大して示す図である。It is a figure which expands and shows the tree crown information image displayed on the display partially. 樹冠の日向部の分光反射特性の一例を示す図である。It is a figure which shows an example of the spectral reflection characteristic of the sunlight part of a tree crown. ディスプレイ上に表示された森林域全体の樹種分類画像を示す図である。It is a figure which shows the tree classification image of the whole forest area displayed on the display. ディスプレイ上に表示された各樹冠の画素ごとに樹種分類を行った樹種分類結果を示す図である。It is a figure which shows the tree species classification result which performed tree species classification for every pixel of each crown displayed on the display. 図9に示した画素ごとの樹種分類結果に基づいて1つの樹冠を1つ樹種として分類した樹種分類結果を示す図である。It is a figure which shows the tree species classification result which classified one crown as one tree species based on the tree species classification result for every pixel shown in FIG. 実施形態2に係る森林現況情報の作成システムの構成を示す図である。It is a figure which shows the structure of the preparation system of the forest present condition information which concerns on Embodiment 2. FIG. 実施形態2に係る森林現況情報の作成システムにおける森林現況情報作成を説明するフローチャートである。10 is a flowchart for explaining forest present state information creation in the forest present state information creating system according to the second embodiment. ディスプレイ上に表示された全樹種樹冠占有画像を示す図である。It is a figure which shows the whole tree species tree crown occupation image displayed on the display. ディスプレイ上に表示された樹種別樹冠占有画像を示す図である。It is a figure which shows the tree classification tree occupation image displayed on the display. ディスプレイ上に表示された全樹種樹冠本数分布画像を示す図である。It is a figure which shows the whole tree species tree crown number distribution image displayed on the display. ディスプレイ上に表示された樹種別樹冠本数分布画像を示す図である。It is a figure which shows the tree classification tree number distribution image displayed on the display. ディスプレイ上に表示された森林現況情報画像を示す図である。It is a figure which shows the forest present condition information image displayed on the display. 実施形態3に係る間伐対象区域の選定システムの構成を示す図である。It is a figure which shows the structure of the selection system of the thinning object area which concerns on Embodiment 3. FIG. 実施形態3に係る間伐対象区域の選定システムにおける間伐対象区域選定処理を説明するフローチャートである。It is a flowchart explaining the thinning object area selection process in the selection system of the thinning object area which concerns on Embodiment 3. FIG. ディスプレイ上に表示された森林現況小班区画画像を示す図である。It is a figure which shows the forest present condition small block division image displayed on the display. ディスプレイ上に表示された小班区画属性情報画像を示す図である。It is a figure which shows the small segment division attribute information image displayed on the display. 図21における領域Dの拡大画像を示す図である。It is a figure which shows the enlarged image of the area | region D in FIG. 各小班区画から抽出された属性情報を示す図である。It is a figure which shows the attribute information extracted from each subdivision block. ディスプレイ上に表示された地利条件画像を示す図である。It is a figure which shows the land condition image displayed on the display. ディスプレイ上に表示された間伐評価モデル化画像の一例を示す図である。It is a figure which shows an example of the thinning evaluation modeled image displayed on the display. ディスプレイ上に表示された間伐対象区域選定画像を示す図である。It is a figure which shows the thinning object area selection image displayed on the display.

符号の説明Explanation of symbols

100・・・樹種分類システム、110・・・撮影画像データ入力部、120・・・森林域画像データ作成部、130・・・樹冠画像データ作成部、140・・・樹冠情報画像データ作成部、150・・・樹種分類画像データ作成部、200・・・森林現況情報の作成システム、210・・・全樹種樹冠占有画像データ作成部、220・・・樹種別樹冠占有画像データ作成部、230・・・全樹種樹冠本数分布画像データ作成部、240・・・樹種別樹冠本数分布画像データ作成部、250・・・森林現況情報画像データ作成部、300・・・間伐対象区域の選定システム、310・・・森林現況小班区画画像データ作成部、320・・・小班区画属性情報画像データ作成部、330・・・間伐対象区域画像データ作成部、500・・・地理情報システム(GIS)、A・・・調査対象地域、B・・・・森林域、C・・・非森林域、R・・・林道、Pe・・・樹頂部、Su・・・日向部、Sh・・・日陰部、T1,T2・・・樹冠   DESCRIPTION OF SYMBOLS 100 ... Tree kind classification system, 110 ... Shooting image data input part, 120 ... Forest area image data creation part, 130 ... Tree crown image data creation part, 140 ... Tree crown information image data creation part, 150 ... Tree species classification image data creation unit, 200 ... Forest status information creation system, 210 ... All tree species crown occupation image data creation unit, 220 ... Tree type crown occupation image data creation unit, 230 .. All tree species crown number distribution image data creation unit, 240 ... Tree type crown number distribution image data creation unit, 250 ... Forest status information image data creation unit, 300 ... Thinning target area selection system, 310 ... Forest present state subcompartment section image data creation section, 320 ... Subcompartment section attribute information image data creation section, 330 ... Thinning target area image data creation section, 500 ... Geography Information system (GIS), A ... Survey area, B ... Forest area, C ... Non-forest area, R ... Forest road, Pe ... Top of tree, Su ... Hyuga part, Sh: Shade, T1, T2: Crown

Claims (11)

調査対象となる森林を含む地域を上空から撮影して得られた撮影画像データと前記森林の地理的情報とに基づいて抽出された森林域に関する森林域画像データを用いて前記森林域に含まれる各樹木を種類別に分類する樹種分類方法であって、
前記森林域画像データに基づいて前記森林域に含まれる各樹冠を抽出し、前記各樹冠に関する樹冠画像データを作成するステップと、
前記樹冠画像データに基づいて前記各樹冠の日向部を含む情報を樹冠情報として抽出し、前記樹冠情報に関する樹冠情報画像データを作成するステップと、
前記樹冠情報画像データに基づいて前記各樹冠の日向部とされた領域における分光反射特性を作成し、前記分光反射特性に基づいて当該樹冠の樹種を特定して樹種分類画像データを作成するステップと、
を有すること特徴とする樹種分類方法。
Included in the forest area using the forest area image data regarding the forest area extracted based on the captured image data obtained by photographing the area including the forest to be surveyed from the sky and the geographical information of the forest A tree species classification method for classifying each tree by type,
Extracting each canopy included in the forest area based on the forest area image data, and creating canopy image data relating to each canopy;
Extracting information including the sunny part of each canopy based on the canopy image data as canopy information, and creating canopy information image data relating to the canopy information;
Creating a spectral reflection characteristic in a region defined as a sunny part of each canopy based on the canopy information image data, identifying a tree type of the canopy based on the spectral reflection characteristic, and creating tree species classification image data; ,
A tree species classification method characterized by comprising:
請求項1に記載の樹種分類方法において、
前記各樹冠の日向部は、前記樹冠画像データにおける輝度値に基づいて抽出され、前記輝度値が当該樹冠において所定値以上である箇所を前記日向部とすることを特徴とする樹種分類方法。
The tree species classification method according to claim 1,
A tree species classification method characterized in that the sunny part of each canopy is extracted based on the brightness value in the canopy image data, and the sunny part is a place where the brightness value is equal to or greater than a predetermined value in the canopy.
請求項1又は2に記載の樹種分類方法において、
前記各樹冠から複数の樹種が抽出された場合、前記各樹冠から抽出された複数の樹種のそれぞれの占有面積を各樹種ごとに算出し、占有面積の多い樹種を当該樹冠の樹種とみなすことを特徴とする樹種分類方法。
The tree species classification method according to claim 1 or 2,
When a plurality of tree species are extracted from each canopy, the respective occupied area of the plurality of tree species extracted from each canopy is calculated for each tree species, and a tree species having a large occupied area is regarded as the tree species of the crown. Characteristic tree classification method.
調査対象となる森林を含む地域を上空から撮影して得られた撮影画像データと前記森林の地理的情報とに基づいて抽出された森林域に関する森林域画像データを用いて前記森林域に含まれる各樹木の種類を分類する樹種分類システムであって、
前記森林域画像データに基づいて前記森林域に含まれる各樹冠を抽出し、前記各樹冠に関する樹冠画像データを作成する樹冠画像データ作成部と、
前記樹冠画像データに基づいて前記各樹冠の日向部を含む情報を樹冠情報として抽出し、前記樹冠情報に関する樹冠情報画像データを作成する樹冠情報画像データ作成部と、
前記樹冠情報画像データに基づいて前記各樹冠の日向部とされた領域における分光反射特性を作成し、前記分光反射特性に基づいて当該樹冠の樹種を特定して樹種分類画像データを作成する樹種分類画像データ作成部と、
を有すること特徴とする樹種分類システム。
Included in the forest area using the forest area image data regarding the forest area extracted based on the captured image data obtained by photographing the area including the forest to be surveyed from the sky and the geographical information of the forest A tree classification system for classifying each tree type,
Extracting each canopy included in the forest area based on the forest area image data, creating a canopy image data creation unit for creating canopy image data relating to each canopy;
Information including the sunny part of each canopy based on the canopy image data is extracted as canopy information, and a canopy information image data creating unit for creating canopy information image data related to the canopy information;
A tree species classification that creates spectral reflection characteristics in a region that is the sunny part of each crown based on the canopy information image data, and identifies tree species of the crown based on the spectral reflection characteristics and creates tree classification image data An image data creation unit;
A tree species classification system characterized by comprising:
調査対象となる森林を含む地域を上空から撮影して得られた撮影画像データと前記森林の地理的情報とに基づいて抽出された森林域に関する森林域画像データを用いて前記森林の現況に関する森林現況情報を作成するための森林現況情報の作成方法であって、
前記森林域画像データに基づいて前記森林域に含まれる各樹冠を抽出し、前記各樹冠に関する樹冠画像データを作成するステップと、
前記樹冠画像データに基づいて前記各樹冠の日向部を含む情報を樹冠情報として抽出し、前記樹冠情報に関する樹冠情報画像データを作成するステップと、
前記樹冠情報画像データに基づいて前記各樹冠の日向部とされた領域における分光反射特性を作成し、前記分光反射特性に基づいて当該樹冠の樹種を特定して樹種分類画像データを作成するステップと、
前記樹種分類画像データに基づいて単位面積当たりの全樹種樹冠面積占有度を算出し、前記全樹種樹冠面積占有度に基づく全樹種樹冠占有画像データを作成するステップと、
前記樹種分類画像データに基づいて単位面積当たりの樹種別樹冠面積占有度を算出し、前記樹種別樹冠面積占有度に基づく樹種別樹冠占有画像データを作成するステップと、
前記樹種分類画像データに基づいて算出された全樹種の樹冠本数から単位面積当たりの全樹種樹冠本数分布画像データを作成するステップと、
前記樹種分類画像データに基づいて算出された樹種別の樹冠本数から単位面積当たりの樹種別樹冠本数分布画像データを作成するステップと、
前記全樹種樹冠占有画像データ、前記樹種別樹冠占有画像データ、前記全樹種樹冠本数分布画像データ、前記樹種別樹冠本数分布画像データ及び前記樹種分類画像データに基づいて森林の現況情報に関する森林現況情報画像データを作成するステップと、
を有することを特徴とする森林現況情報の作成方法。
Forests related to the current state of the forest by using forest image data relating to the forest area extracted based on the photographed image data obtained by photographing the area including the forest to be surveyed from above and the geographical information of the forest A method for creating forest status information for creating status information,
Extracting each canopy included in the forest area based on the forest area image data, and creating canopy image data relating to each canopy;
Extracting information including the sunny part of each canopy based on the canopy image data as canopy information, and creating canopy information image data relating to the canopy information;
Creating a spectral reflection characteristic in a region defined as a sunny part of each canopy based on the canopy information image data, identifying a tree type of the canopy based on the spectral reflection characteristic, and creating tree species classification image data; ,
Calculating the total tree species crown area occupancy per unit area based on the tree species classification image data, and creating all tree species crown occupancy image data based on the total tree species crown area occupancy;
Calculating a tree-type crown area occupancy per unit area based on the tree-type classification image data, and creating tree-type crown occupancy image data based on the tree-type crown area occupancy;
Creating all tree species crown number distribution image data per unit area from the number of tree crowns calculated based on the tree species classification image data;
Creating tree-type crown number distribution image data per unit area from the number of tree-type crowns calculated based on the tree-type classification image data;
Forest present state information on forest present state information based on the all tree species crown occupation image data, the tree type crown occupation image data, the all tree species crown number distribution image data, the tree type crown number distribution image data, and the tree species classification image data Creating image data;
A method for creating forest status information, characterized by comprising:
調査対象となる森林を含む地域を上空から撮影して得られた撮影画像データと前記森林の地理的情報とに基づいて抽出された森林域に関する森林域画像データを用いて前記森林の現況に関する森林現況情報を作成するための森林現況情報の作成システムであって、
前記森林域画像データに基づいて前記森林域に含まれる各樹冠を抽出し、前記各樹冠に関する樹冠画像データを作成する樹冠画像データ作成部と、
前記樹冠画像データに基づいて前記各樹冠の日向部を含む情報を樹冠情報として抽出し、前記樹冠情報に関する樹冠情報画像データを作成する樹冠情報画像データ作成部と、
前記樹冠情報画像データに基づいて前記各樹冠の日向部とされた領域における分光反射特性を作成し、前記分光反射特性に基づいて当該樹冠の樹種を特定して樹種分類画像データを作成する樹種分類画像データ作成部と、
前記樹種分類画像データに基づいて全樹種の樹冠本数を算出する機能と、単位面積当たりの全樹種樹冠面積占有度を算出し、前記全樹種樹冠面積占有度に基づく全樹種樹冠占有画像データを作成する機能とを有する全樹種樹冠占有画像データ作成部と、
前記樹種分類画像データに基づいて樹種別の樹冠本数を算出する機能と、単位面積当たりの樹種別樹冠面積占有度を算出し、前記樹種別樹冠面積占有度に基づく樹種別樹冠占有画像データを作成する機能とを有する樹種別樹冠占有画像データ作成部と、
前記全樹種の樹冠本数から単位面積当たりの全樹種樹冠本数分布画像データを作成する全樹種樹冠本数分布画像データ作成部と、
前記樹種別の樹冠本数から単位面積当たりの樹種別樹冠本数分布画像データを作成する樹種別樹冠本数分布画像データ作成部と、
前記全樹種樹冠占有画像データ、前記樹種別樹冠占有画像データ、前記全樹種樹冠本数分布画像データ、前記樹種別樹冠本数分布画像データ及び前記樹種分類画像データに基づいて森林の現況情報に関する森林現況情報画像データを作成する森林現況情報画像データ作成部と、
を有することを特徴とする森林現況情報の作成システム。
Forests related to the current state of the forest by using forest image data relating to the forest area extracted based on the photographed image data obtained by photographing the area including the forest to be surveyed from above and the geographical information of the forest A system for creating forest status information for creating current status information,
Extracting each canopy included in the forest area based on the forest area image data, creating a canopy image data creation unit for creating canopy image data relating to each canopy;
Information including the sunny part of each canopy based on the canopy image data is extracted as canopy information, and a canopy information image data creating unit for creating canopy information image data related to the canopy information;
A tree species classification that creates spectral reflection characteristics in a region that is the sunny part of each crown based on the canopy information image data, and identifies tree species of the crown based on the spectral reflection characteristics and creates tree classification image data An image data creation unit;
A function for calculating the number of crowns of all tree species based on the tree species classification image data, calculating the total tree species crown area occupancy per unit area, and creating all tree species crown occupancy image data based on the total tree species crown area occupancy An all-tree species crown occupancy image data creation unit having a function to
A function for calculating the number of crowns of a tree type based on the tree type classification image data and a tree type canopy area occupancy per unit area, and creating tree type canopy occupancy image data based on the tree type canopy area occupancy A tree-type crown occupancy image data creation unit having a function to
A total tree species crown number distribution image data creation unit for creating all tree species crown number distribution image data per unit area from the number of crowns of all tree species;
Tree-type crown number distribution image data creation unit for creating tree-type crown number distribution image data per unit area from the number of tree-type crowns;
Forest present state information on forest present state information based on the all tree species crown occupation image data, the tree type crown occupation image data, the all tree species crown number distribution image data, the tree type crown number distribution image data, and the tree species classification image data Forest status information image data creation section that creates image data;
A system for creating forest status information, characterized by comprising:
調査対象となる森林の現況に関する森林現況情報画像データに基づいて間伐すべき間伐対象区域を選定する間伐対象区域の選定方法であって、
前記森林現況情報画像データに基づいて前記森林を小班区画に分別するための小班区画画像データを作成して、前記小班区画画像データと前記森林現況情報画像データとに基づいて森林現況小班区画画像データを作成するステップと、
前記各小班区画の特徴を現す小班区画属性情報を抽出し、前記各小班区画属性情報に関する小班区画属性情報画像データを作成するステップと、
前記小班区画属性情報画像データと当該森林の地利条件とに基づいて間伐対象区域を選定し、選定された間伐対象区域を表す間伐対象区域画像データを作成するステップと、
を有することを特徴とする間伐対象区域の選定方法。
A method of selecting a thinning target area, which selects a thinning target area to be thinned based on image data of the forest present state information regarding the current state of the forest to be investigated,
Based on the forest present state information image data, a subcompartment image data for dividing the forest into subcompartments is created, and based on the subcompartment image data and the forest present state information image data, the forest present state subcompartment image data The steps of creating
Extracting the subcompartment section attribute information representing the characteristics of each subcompartment section, creating subcompartment section attribute information image data related to each subsection section attribute information;
Selecting a thinning target area based on the subdivision section attribute information image data and the land condition of the forest, and creating thinning target area image data representing the selected thinning target area;
The selection method of the thinning target area characterized by having.
請求項7に記載の間伐対象区域の選定方法において、
前記森林現況情報画像データは、請求項5に記載の森林現況情報の作成方法によって作成された森林現況情報画像データであることを特徴とする間伐対象区域の選定方法。
In the selection method of the thinning object area of Claim 7,
The method for selecting a thinning target area, wherein the forest present state information image data is forest present state information image data created by the forest present state information creating method according to claim 5.
請求項7又は8に記載の間伐対象区域の選定方法において、
前記小班区画属性情報は、当該小班区画における樹種別の樹冠本数、単位面積あたりの樹種別の樹冠本数を含むことを特徴とする間伐対象区域の選定方法。
In the selection method of the thinning object area of Claim 7 or 8,
The method of selecting a thinning target area, wherein the sub-compartment attribute information includes the number of tree-type crowns in the sub-compartment and the number of tree-type crowns per unit area.
請求項7〜9のいずれかに記載の間伐対象区域の選定方法において、
前記森林の地利条件は、当該小班区画の林道からの距離、標高、傾斜を含むことを特徴とする間伐対象区域の選定方法。
In the selection method of the thinning object area in any one of Claims 7-9,
The method for selecting a thinning target area, wherein the land condition of the forest includes a distance, an altitude, and a slope from the forest road of the subcompartment.
調査対象となる森林の現況に関する森林現況情報画像データに基づいて間伐すべき間伐対象区域を選定する間伐対象区域の選定システムであって、
前記森林現況情報画像データに基づいて前記森林を小班区画に分別するための小班区画画像データを作成して、前記小班区画画像データと前記森林現況情報画像データとに基づいて森林現況小班区画画像データを作成する森林現況小班区画画像データ作成部と、
前記各小班区画の特徴を現す小班区画属性情報を抽出し、前記各小班区画属性情報に関する小班区画属性情報画像データを作成する小班区画属性情報画像データ作成部と、
前記小班区画属性情報画像データと当該森林の地利条件とに基づいて間伐対象区域を選定し、選定された間伐対象区域を表す間伐対象区域画像データを作成する間伐対象区域画像データを作成部と、
を有することを特徴とする間伐対象区域の選定システム。
A system for selecting a thinning target area for selecting a thinning target area to be thinned based on image data of the forest present state information regarding the current state of the forest to be investigated,
Based on the forest present state information image data, a subcompartment image data for dividing the forest into subcompartments is created, and based on the subcompartment image data and the forest present state information image data, the forest present state subcompartment image data Forest current subcompartment block image data creation section to create,
Extracting subdivision section attribute information representing the characteristics of each subsection section, and generating subgroup section attribute information image data related to each subsection section attribute information.
A thinning target area image data for creating a thinning target area image data representing a selected thinning target area, selecting a thinning target area based on the subdivision section attribute information image data and the land condition of the forest,
A system for selecting areas to be thinned, characterized by comprising:
JP2008254520A 2008-09-30 2008-09-30 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 Active JP4900356B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008254520A JP4900356B2 (en) 2008-09-30 2008-09-30 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

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008254520A JP4900356B2 (en) 2008-09-30 2008-09-30 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

Publications (2)

Publication Number Publication Date
JP2010086276A true JP2010086276A (en) 2010-04-15
JP4900356B2 JP4900356B2 (en) 2012-03-21

Family

ID=42250175

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008254520A Active JP4900356B2 (en) 2008-09-30 2008-09-30 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

Country Status (1)

Country Link
JP (1) JP4900356B2 (en)

Cited By (21)

* 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
JP2012208046A (en) * 2011-03-30 2012-10-25 Fujitsu Ltd Plant species identification method and plant species identification device
WO2013002349A1 (en) 2011-06-29 2013-01-03 富士通株式会社 Plant species identification device, method, and program
JP2013054499A (en) * 2011-09-02 2013-03-21 Asia Air Survey Co Ltd Forest stereoscopic image generation method
JP2013090584A (en) * 2011-10-25 2013-05-16 Fujitsu Ltd Program for plant species specification, information processing method, and device
CN103155836A (en) * 2011-12-16 2013-06-19 东北林业大学 Method for forecasting forest pest occurrence degree
WO2014142317A1 (en) * 2013-03-15 2014-09-18 株式会社パスコ Forest type analysis device, forest type analysis method and program
KR101453965B1 (en) 2013-04-24 2014-10-28 (주)태은 Analysis of basic information for National Forest forest survey and Input system of attribute data in description of National Forest Management Planning
JP2015026369A (en) * 2013-06-20 2015-02-05 株式会社パスコ Forest physiognomy analyzing device, forest physiognomy analysis method, and program
JP2015084191A (en) * 2013-10-25 2015-04-30 株式会社パスコ Forest physiognomy analyzing apparatus, forest physiognomy analyzing method, and program
JP2015109064A (en) * 2013-10-25 2015-06-11 株式会社パスコ Forest type analyzer, forest type analysis method, and program
JP2015152340A (en) * 2014-02-12 2015-08-24 株式会社パスコ Forest type analysis apparatus, forest type analysis method and program
JP2017169520A (en) * 2016-03-25 2017-09-28 忠士 岩下 Grape cultivation management method
WO2019082519A1 (en) * 2017-10-26 2019-05-02 ソニー株式会社 Information processing device, information processing method, program, and information processing system
JP2019095910A (en) * 2017-11-20 2019-06-20 株式会社パスコ Erroneous discrimination possibility evaluation apparatus, erroneous discrimination possibility evaluation method, and program
JP2020018274A (en) * 2018-08-03 2020-02-06 一般財団法人電力中央研究所 Benthos grasping system and benthos grasping method
JP2020098442A (en) * 2018-12-18 2020-06-25 株式会社パスコ Tree width setting device and program
CN111611443A (en) * 2020-04-20 2020-09-01 东北林业大学 Tree species space isolation degree evaluation method based on neighbor relation
KR20200120403A (en) * 2019-04-12 2020-10-21 대한민국(산림청 국립산림과학원장) System and method of tree species classification using satellite image
WO2022252799A1 (en) * 2021-06-04 2022-12-08 成都数之联科技股份有限公司 Model training method, woodland change detection method, system, and apparatus, and medium
CN115862304A (en) * 2023-03-03 2023-03-28 吉林省林业科学研究院 Intelligent early warning system and method for ecological restoration of degenerated natural secondary forest

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105354534B (en) * 2015-09-29 2018-11-20 南京林业大学 A kind of tree species classification method based on multi-source same period high-definition remote sensing data
JP7169640B2 (en) * 2018-11-16 2022-11-11 国立大学法人信州大学 Thinned tree selection method and thinned tree selection device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001357380A (en) * 2000-06-13 2001-12-26 Pasuko:Kk Method for evaluating forest region by picture processing and storage medium with stored program concerning the evaluation
JP2003344048A (en) * 2002-05-22 2003-12-03 Pasuko:Kk System for processing forest information
JP2006285310A (en) * 2005-03-31 2006-10-19 Kanazawa Univ Evaluation method of canopy of forest, and its canopy evaluation program

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001357380A (en) * 2000-06-13 2001-12-26 Pasuko:Kk Method for evaluating forest region by picture processing and storage medium with stored program concerning the evaluation
JP2003344048A (en) * 2002-05-22 2003-12-03 Pasuko:Kk System for processing forest information
JP2006285310A (en) * 2005-03-31 2006-10-19 Kanazawa Univ Evaluation method of canopy of forest, and its canopy evaluation program

Cited By (33)

* 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
JP2012208046A (en) * 2011-03-30 2012-10-25 Fujitsu Ltd Plant species identification method and plant species identification device
JP5673823B2 (en) * 2011-06-29 2015-02-18 富士通株式会社 Plant species identification device, method and program
WO2013002349A1 (en) 2011-06-29 2013-01-03 富士通株式会社 Plant species identification device, method, and program
CN103620381A (en) * 2011-06-29 2014-03-05 富士通株式会社 Plant species identification device, method, and program
US20140093138A1 (en) * 2011-06-29 2014-04-03 Fujitsu Limited Plant species identification apparatus and method
US9230170B2 (en) 2011-06-29 2016-01-05 Fujitsu Limited Plant species identification apparatus and method
CN103620381B (en) * 2011-06-29 2015-09-30 富士通株式会社 Floristics recognition device and method
JP2013054499A (en) * 2011-09-02 2013-03-21 Asia Air Survey Co Ltd Forest stereoscopic image generation method
JP2013090584A (en) * 2011-10-25 2013-05-16 Fujitsu Ltd Program for plant species specification, information processing method, and device
CN103155836A (en) * 2011-12-16 2013-06-19 东北林业大学 Method for forecasting forest pest occurrence degree
CN103155836B (en) * 2011-12-16 2014-09-24 东北林业大学 Method for forecasting forest pest occurrence degree
WO2014142317A1 (en) * 2013-03-15 2014-09-18 株式会社パスコ Forest type analysis device, forest type analysis method and program
JP2014199660A (en) * 2013-03-15 2014-10-23 株式会社パスコ Forest physiognomy analysis device, forest physiognomy analysis method and program
KR101453965B1 (en) 2013-04-24 2014-10-28 (주)태은 Analysis of basic information for National Forest forest survey and Input system of attribute data in description of National Forest Management Planning
JP2015026369A (en) * 2013-06-20 2015-02-05 株式会社パスコ Forest physiognomy analyzing device, forest physiognomy analysis method, and program
JP2015084191A (en) * 2013-10-25 2015-04-30 株式会社パスコ Forest physiognomy analyzing apparatus, forest physiognomy analyzing method, and program
JP2015109064A (en) * 2013-10-25 2015-06-11 株式会社パスコ Forest type analyzer, forest type analysis method, and program
JP2015152340A (en) * 2014-02-12 2015-08-24 株式会社パスコ Forest type analysis apparatus, forest type analysis method and program
JP2017169520A (en) * 2016-03-25 2017-09-28 忠士 岩下 Grape cultivation management method
WO2019082519A1 (en) * 2017-10-26 2019-05-02 ソニー株式会社 Information processing device, information processing method, program, and information processing system
US11823447B2 (en) 2017-10-26 2023-11-21 Sony Group Corporation Information processing apparatus, information processing method, program, and information processing system
JPWO2019082519A1 (en) * 2017-10-26 2020-11-19 ソニー株式会社 Information processing equipment, information processing methods, programs, information processing systems
JP2019095910A (en) * 2017-11-20 2019-06-20 株式会社パスコ Erroneous discrimination possibility evaluation apparatus, erroneous discrimination possibility evaluation method, and program
JP2020018274A (en) * 2018-08-03 2020-02-06 一般財団法人電力中央研究所 Benthos grasping system and benthos grasping method
JP7248419B2 (en) 2018-12-18 2023-03-29 株式会社パスコ Tree width setting device and program
JP2020098442A (en) * 2018-12-18 2020-06-25 株式会社パスコ Tree width setting device and program
KR20200120403A (en) * 2019-04-12 2020-10-21 대한민국(산림청 국립산림과학원장) System and method of tree species classification using satellite image
KR102239393B1 (en) * 2019-04-12 2021-04-12 대한민국 System and method of tree species classification using satellite image
CN111611443B (en) * 2020-04-20 2023-04-14 东北林业大学 Tree species space isolation degree evaluation method based on neighbor relation
CN111611443A (en) * 2020-04-20 2020-09-01 东北林业大学 Tree species space isolation degree evaluation method based on neighbor relation
WO2022252799A1 (en) * 2021-06-04 2022-12-08 成都数之联科技股份有限公司 Model training method, woodland change detection method, system, and apparatus, and medium
CN115862304A (en) * 2023-03-03 2023-03-28 吉林省林业科学研究院 Intelligent early warning system and method for ecological restoration of degenerated natural secondary forest

Also Published As

Publication number Publication date
JP4900356B2 (en) 2012-03-21

Similar Documents

Publication Publication Date Title
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
JP4858793B2 (en) Tree number calculation method and tree number calculation device
Dorren et al. Improved Landsat-based forest mapping in steep mountainous terrain using object-based classification
Ismail et al. Satellite data classification accuracy assessment based from reference dataset
CA2656595C (en) Remote sensing and probabilistic sampling based forest inventory method
Malambo et al. Automated detection and measurement of individual sorghum panicles using density-based clustering of terrestrial lidar data
US7212670B1 (en) Method of feature identification and analysis
Zhang et al. Estimation of forest leaf area index using height and canopy cover information extracted from unmanned aerial vehicle stereo imagery
Mitri et al. A semi-automated object-oriented model for burned area mapping in the Mediterranean region using Landsat-TM imagery
Marcelino et al. Landslide inventory using image fusion techniques in Brazil
JP6544582B2 (en) Calculation method of damage classification of pine pileworm and damage classification calculation device of pine pileworm
Aliero et al. The usefulness of unmanned airborne vehicle (UAV) imagery for automated palm oil tree counting
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
JP2011024471A (en) Classification of kind of tree in forest by utilizing intermediate infrared area band
Piermattei et al. Pléiades satellite images for deriving forest metrics in the Alpine region
Mitri et al. Mapping the severity of fire using object-based classification of IKONOS imagery
Yin et al. Individual tree parameters estimation for chinese fir (cunninghamia lanceolate (lamb.) hook) plantations of south china using UAV Oblique Photography: Possibilities and Challenges
Smigaj et al. Capturing hedgerow structure and flowering abundance with UAV remote sensing
Katoh et al. Application of high-resolution airborne data using individual tree crowns in Japanese conifer plantations
Silva et al. Mapping two competing grassland species from a low-altitude Helium balloon
Arenas-Corraliza et al. Automatic mapping of tree crowns in scattered-tree woodlands using low-density LiDAR data and infrared imagery
CN115760885B (en) High-closure-degree wetland forest parameter extraction method based on consumer-level unmanned aerial vehicle image
JP2001357380A (en) Method for evaluating forest region by picture processing and storage medium with stored program concerning the evaluation
Gong et al. Vineyard identification in an oak woodland landscape with airborne digital camera imagery
van Geffen et al. SiDroForest: A comprehensive forest inventory of Siberian boreal forest investigations including drone-based point clouds, individually labelled trees, synthetically generated tree crowns and Sentinel-2 labelled image patches

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110616

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20110616

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20110726

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110823

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110914

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20111122

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20111219

R150 Certificate of patent or registration of utility model

Ref document number: 4900356

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150113

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250