JP7169640B2 - 間伐木の選定方法及び間伐木の選定装置 - Google Patents

Description

本発明は、間伐対象となる森林域から間伐木を選定する間伐木の選定方法及び間伐木の選定装置に関する。

植栽によって形成された森林においては、植栽した木が大きくなるにつれ、木と木の間が混み合ってきて日光が当たらなくなり、木が成長できなくなる。そこで間伐を行い、残存木の成長を促すことが必要である。このような間伐は植栽されて２０年程度経過したところで、形質の悪い木を伐採し、優良木を残す定性間伐が一般的である。

定性間伐を行う際には、間伐木を選定する必要があるが、現地での間伐木の選定は熟練性と多くの人手を要するとともに、蜂や熊など危険な生物との遭遇、また、斜面の笹薮など足場の悪い中での選定作業など種々の危険が伴う。また、森林全体を見て間伐木を選定するものではないために、樹木の密度が高い範囲や密度が低い範囲も同じ割合で選木されることとなり、樹木の密度の違いが解消できず、残存木が適正配置とならない場合も多い。

また、間伐木選定結果は個人差によるバラツキも多く、間伐すべき樹木の見落としや間伐すべき樹木が集中してしまうなど課題もある。また、現場経験を重視するため、若い技術者や女性の新規参入の障害にもなっている。

間伐木の選定は、広大な森林域において設定された間伐対象となる森林域（間伐対象区域という。）内で行われるが、このような間伐対象区域を選定する方法は、すでに提案されている（例えば、特許文献１参照。）。

特許文献１に記載されている間伐対象区域の選定方法は、当該特許文献１に記載されている森林現況情報の作成方法によって作成された森林現況情報画像データに基づいて、間伐対象区域を選定するものである。

すなわち、特許文献１に記載されている間伐対象区域の選定方法は、森林現況情報画像データに基づいて分別された各小班区画における各樹冠の特徴を現す小班区画属性情報を抽出し、抽出された各小班区画属性情報と当該森林の地理条件とに基づいて間伐対象区域を選定するようにしている。これにより、従来においては、個人の労力と熟練に頼っていた間伐対象区域の選定を客観的かつ適切に選定することできる。

また、特許文献１に記載されている間伐対象区域の選定方法は、現時点における実際の森林の状況として、例えば、樹冠面積が大きい、占有面積割合が高い、立木本数が多いなどが考慮され、かつ、間伐作業のし易さとして、例えば、林道から近い、標高があまり高くない、傾斜が緩いなどが考慮された間伐対象区域を選定することができる。なお、「樹冠」とは森林を構成する各樹木を上空から見たときに個々の樹木の枝と葉で構成された樹木の部分をいう。

このように、特許文献１に記載されている間伐対象区域の選定方法は、広大な森林域から、間伐すべき間伐対象区域を選定する場合には、優れた間伐対象区域の選定方法であるといえる。

特許４９００３５６号公報

しかしながら、特許文献１に記載されている間伐対象区域の選定方法は、間伐対象区域を適正かつ効率よく選定可能となるが、実際に間伐すべき間伐木の選定をどのようにして行うかについては示されていない。

前述したように、現地での間伐木の選定は熟練性と多くの人手を要し、危険も多く、また、残存木が適正配置とならない場合も多いといった課題や、間伐木選定結果は個人差によるバラツキも多く、間伐すべき樹木の見落としや間伐すべき樹木が集中してしまうなど課題が多い。このため、現地での人手による間伐木の選定に頼ることなく、間伐木の選定を効率的にかつ森林の現況に適応するように行うことが要求されている。

そこで本発明は、現地での人手による間伐木の選定に頼ることなく、間伐木の選定を効率的にかつ森林の現況に適応するように行うことができる間伐木の選定方法及び間伐木の選定装置を提供することを目的とする。

［１］本発明の間伐木の選定方法は、間伐対象区域に存在する樹木に関する情報が個々の樹木ごとに記録されている森林資源情報に基づいて、前記間伐対象区域に存在する樹木のうち、間伐すべき樹木を間伐木として選定する間伐木の選定方法であって、前記間伐対象区域における間伐率を設定する間伐率設定ステップと、前記間伐対象区域においてどのような間伐を行うかを示す間伐条件を設定する間伐条件設定ステップと、前記間伐対象区域に存在する樹木の個々の樹木間の距離を算定する樹木間距離算定ステップと、前記樹木間距離算定ステップで算定された距離のうち、最も近接している２本の樹木を最短距離樹木として選択する最短距離樹木選択ステップと、前記間伐条件設定ステップによって設定されている前記間伐条件に基づいて、前記最短距離樹木選択ステップによって選択されている２本の樹木のうちの１本の樹木を前記間伐木として選定する間伐木選定ステップと、を有し、前記樹木間距離算定ステップ、前記最短距離樹木選択ステップ及び前記間伐木選定ステップは、前記間伐率設定ステップにおいて設定されている前記間伐率に達するまで順次行うことを特徴とする。

本発明の間伐木選定方法によれば、このような手順によって間伐木の選定を行うことにより、現地での人手による間伐木の選定に頼ることなく、間伐木の選定を効率的にかつ森林の現況に適応するように行うことができる。

［２］本発明の間伐木の選定方法においては、前記間伐率は、前記間伐対象区域に存在する樹木の本数に基づく間伐率であることが好ましい。

これにより、間伐率を樹木の本数によって管理できる。例えば、間伐対象区域に存在する樹木の本数が３００本であって、当該間伐対象区域に存在する樹木の本数に基づく間伐率を３０％と設定した場合には、間伐木として選定すべき樹木の本数は９０本となる。この場合、間伐木として選定すべき樹木の本数を、設定された間伐率未満に抑えるということも可能である。

［３］本発明の間伐木の選定方法においては、前記間伐率は、前記間伐対象区域に存在する樹木の合計材積に基づく間伐率であることも好ましい。

これにより、間伐率を樹木の材積によって管理できる。例えば、間伐対象区域に存在する樹木の合計材積が１００ｍ^３であって、当該間伐対象区域に存在する樹木の合計材積に基づく間伐率を３５％と設定した場合には、間伐木として選定すべき樹木の合計材積は３５ｍ^３となる。この場合、間伐木として選定すべき樹木の合計材積を、設定された間伐率未満に抑えるということも可能である。

なお、樹木の材積は、例えば、各樹木の胸高直径と樹高の２変数材積式から求めることができる。ここで、胸高直径というのは、樹木の幹の太さを表すものであり、各樹木において人間のほぼ胸の高さの直径を各樹木の胸高直径としている。但し、樹木の幹は実際には真円ではないため、厳密には「直径」とは言えないが、この明細書においては、幹の太さを「直径」と表記することとする。

［４］本発明の間伐木の選定方法においては、前記間伐条件は、上層間伐、下層間伐又は全層間伐のいずれかであり、前記間伐条件として前記上層間伐が設定されている場合には、前記間伐木選定ステップは、前記最短距離樹木選択ステップによって選択されている２本の樹木のうちの大きい方の樹木を前記間伐木として選択し、前記間伐条件として前記下層間伐が設定されている場合には、前記間伐木選定ステップは、前記最短距離樹木選択ステップによって選択されている２本の樹木のうちの小さい方の樹木を前記間伐木として選択し、前記間伐条件として前記全層間伐が設定されている場合には、前記間伐木選定ステップは、前記最短距離樹木選択ステップによって選択されている２本の樹木のうちの１本の樹木を所定の規則に基づいて前記間伐木として選択することが好ましい。

間伐条件として、上層間伐、下層間伐及び全層間伐のいずれかを設定することによって、それぞれの間伐条件に応じた間伐木の選定が可能となる。例えば、上層間伐による間伐木の選定を行うことによって、最短距離として選択された２本の樹木のうち、大きい樹木が間伐木として選定されることとなる。このため、間伐計画の立案者は、木材価値の高い樹木から自動選択することができ、高い木材収入を得ることができる。また間伐後の中層木と下層木の残存木が均等配置されることで、広い生育空間があることから旺盛な生育が期待できる森林となる。

また、下層間伐による間伐木の選定を行うことによって、最短距離として選択された２本の樹木のうち、小さい樹木が間伐木として選定されることとなる。すなわち、生育の障害となる混んだ箇所に生えている木材価値の低い下層木が間伐木として選定されることとなるため、生育の良い樹木が残存木としてそれぞれ所定距離を置いて配置されることとなり、間伐した後の木材価値の向上が期待できる森林になる。

また、全層間伐による間伐木の選定を行うことによって、最短距離として選択されている２本の樹木のうちの１本の樹木が所定の規則に基づいて選択されることとなるため、所定の規則に応じた間伐木の選定が可能となる。

［５］本発明の間伐木の選定方法においては、前記所定の規則は、現時点の処理回数を表す数値が奇数か偶数かによって、前記最短距離樹木選択ステップによって選択されている２本の樹木のうちの大きい方の樹木又は小さい方の樹木のいずれかを前記間伐木として選択することである。

これにより、最短距離として選択された２本の樹木のうち、小さい方の樹木と大きい方の樹木とを交互に選定することができる。すなわち、処理回数を、１，２，３，・・・というように自然数で表した場合、例えば、処理回数を表す数値が偶数の場合には大きい方の樹木を選定し、処理回数を表す数値が奇数の場合には小さい方の樹木を選定するというような規則を設定しておけば、最短距離として選択された２本の樹木のうちの大きい方の樹木又は小さい方の樹木が交互に間伐木として選定されることとなる。このため、樹木の大小に関わらず、生育の障害となる混んだ箇所の樹木を間伐することができ、現在の上層木、中層木、下層木が配置されている森林状態の構成を持続的に維持しながら、生育の障害となる混んだ箇所の樹木を間伐することができる。これにより、全層（上層、中層及び下層）からの木材価値と森林の価値を長期に渡り継続的に得ることができる。

［６］本発明の間伐木の選定方法においては、前記樹木の大きさは、当該樹木の胸高直径で表わされることが好ましい。

このように樹木の大きさを胸高直径で表すことによって、樹木の大きさの大きい方の樹木というのは、幹の太い樹木ということとなり、また、樹木の大きさの小さい方の樹木というのは、幹の細い樹木ということとなる。

［７］本発明の間伐木の選定方法においては、前記樹木間距離算定ステップは、ディープラーニングのニューラルネットワークを用いて行うことが好ましい。

樹木間距離算定ステップが行う樹木間距離算定処理は、間伐対象区域に存在するすべての樹木の個々の樹木間の距離を算定する処理であり、しかも、当該個々の樹木間の距離を算定する処理は、間伐率に達するまで繰り返し行うこととなるため、演算量が膨大なものとなる。このため、樹木間距離算定ステップが行う樹木間算定処理を、ディープラーニング（深層学習）のニューラルネットワークを用いることで、短時間で高精度に間伐対象区域に存在する個々の樹木間の距離を算定することができる。なお、樹木間距離算定ステップで算定された距離のうち、最も近接している２本の樹木を最短距離樹木として選択する最短距離樹木選択ステップまでを含めてディープラーニング（深層学習）のニューラルネットワークを用いて行ってもよい。

［８］本発明の間伐木の選定方法においては、前記間伐対象区域を複数に区分した各区分ごとの傾斜角度を表す傾斜区分画像データを作成する傾斜区分画像データ作成ステップと、前記傾斜区分画像データに基づいて所定の傾斜角度未満の傾斜を有する区分を間伐可能区分とした間伐可能区分画像データを作成する間伐可能区分画像データ作成ステップと、前記間伐可能区分画像データと前記森林資源情報とに基づいて、前記間伐可能区分に存在する樹木を間伐可能樹木として選定する間伐可能樹木選定ステップと、をさらに有し、前記間伐可能樹木選定ステップにより選定された前記間伐可能樹木に対して、前記間伐率設定ステップにより設定されている前記間伐率と、前記間伐条件設定ステップにより設定されている前記間伐条件とに基づいて、前記樹木間距離算定ステップが行う樹木間距離を算定する処理、前記最短距離樹木選択ステップが行う最短距離樹木を選択する処理及び前記間伐木選定ステップが行う間伐木を選定する処理を行うことが好ましい。

このような処理を行うのは、主に、林業用ハーベスタが走行可能な傾斜地に存在している樹木を間伐木として選定できるようにするためである。なお、林業用ハーベスタというのは、一台で伐採、枝払い、玉切りなどの作業をすべて行える高性能収穫機を指している。従って、間伐可能区分画像データを作成して、当該間伐可能区分画像データと森林資源情報から得られる個々の樹木の樹冠位置とを重ね合わせることにより、間伐対象区域内において林業用ハーベスタが走行可能な区域に存在する樹木を選定することができる。このようにして選定した樹木（間伐可能樹木）に対して、樹木間距離算定ステップ、最短距離樹木選択ステップ及び前記間伐木選定ステップを行うことにより、間伐木として選定された樹木は、林業用ハーベスタが走行可能な位置に存在することとなり、当該樹木の伐採、枝払い、玉切りなどの作業の機械化が可能となる。

［９］本発明の間伐木の選定方法においては、前記間伐木選定ステップによって選定した間伐木に関する間伐木資源情報を作成するとともに、当該間伐木資源情報に基づいて間伐木の概要に関する間伐木概要情報を作成する間伐木集計ステップをさらに有し、前記間伐木資源情報には、間伐木として選定された個々の樹木の番号、個々の樹木の位置情報、個々の樹高、個々の樹木の胸高直径、個々の樹木の材積が含まれており、前記間伐木概要情報には、間伐木選定前の樹木に関する情報及び間伐木選定後の選定間伐木に関する情報として、本数、所定面積当たりの本数、平均胸高直径、平均樹高、合計材積、所定面積当たりの合計材積及び選定間伐木の間伐率が含まれていることが好ましい。

このように作成された間伐木資源情報は、間伐木として選定された個々の樹木に対応した高精度な資源情報となる。また、当該間伐木資源情報に基づいて間伐木の概要に関する間伐木概要情報を作成することによって、選定された間伐木の概要を容易に取得できる。

なお、間伐対象区域における森林資源情報（間伐木選定前の森林資源情報）、間伐木選定後の間伐木資源情報及び間伐木概要情報は、それぞれデータベース化してコンピューターに保存させることによって、随時、一覧表としてディスプレイ上に表示させることができる。これにより、間伐木と残存木の分布状況を容易に把握することができる。また、これら各情報をプリントアウトしたものを持参して現地に入ったり、データベース化した情報の表示が可能なモバイル端末を持参して現地に入ったりすることによって、熟練者でなくてもＧＰＳを併用することにより、現地で容易に間伐木として選定した樹木のマーキングが可能となる。

［１０］本発明の間伐木の選定方装置は、間伐対象区域に存在する樹木に関する情報が個々の樹木ごとに記録されている森林資源情報に基づいて、前記間伐対象区域に存在する樹木のうち、間伐すべき樹木を間伐木として選定する間伐木の選定装置であって、前記間伐対象区域における間伐率を設定する間伐率設定部と、前記間伐対象区域においてどのような間伐を行うかを示す間伐条件を設定する間伐条件設定部と、前記間伐対象区域に存在する樹木の個々の樹木間の距離を算定する樹木間距離算定部と、前記樹木間距離算定部で算定された距離のうち、最も近接している２本の樹木を最短距離樹木として選択する最短距離樹木選択部と、前記間伐条件設定部によって設定されている前記間伐条件に基づいて、前記最短距離樹木選択部によって選択されている２本の樹木のうちの１本の樹木を前記間伐木として選定する間伐木選定部と、を有し、前記樹木間距離算定部が行う樹木間距離を算定する処理、前記最短距離樹木選択部が行う最短距離樹木を選択する処理及び前記間伐木選定部が行う間伐木を選定する処理を、前記間伐率設定部によって設定されている前記間伐率に達するまで順次行うことを特徴とする。

本発明の間伐木の選定装置がこのような構成を有していることから、前記［１］に記載の本発明の間伐木の選定方法によって得られる効果と同様の効果が得られる。

［１１］本発明の間伐木の選定方装置においては、前記間伐対象区域を複数に区分した各区分ごとの傾斜角度を表す傾斜区分画像データを作成する傾斜区分画像データ作成部と、前記傾斜区分画像データに基づいて所定の傾斜角度未満の傾斜を有する区分を間伐可能区分とした間伐可能区分画像データを作成する間伐可能区分画像データ作成部と、前記間伐可能区分画像データと前記森林資源情報とに基づいて、前記間伐可能区分に存在する樹木を間伐可能樹木として選定する間伐可能樹木選定部と、をさらに有し、前記間伐可能樹木選定部により選定された前記間伐可能樹木に対して、前記間伐率設定部により設定されている前記間伐率及び前記間伐条件設定部により設定されている前記間伐条件に基づいて、前記樹木間距離算定部が行う樹木間距離を算定する処理、前記最短距離樹木選択部が行う最短距離樹木を選択する処理及び前記間伐木選定部が行う間伐木を選定する処理を行うことが好ましい。

本発明の間伐木の選定装置がこのような構成を有することにより、前記［８］に記載の本発明の間伐木の選定方法によって得られる効果と同様の効果が得られる。

［１２］本発明の間伐木の選定方装置においては、前記間伐木選定部によって選定した間伐木に関する間伐木資源情報を作成するとともに、作成された間伐木資源情報に基づいて間伐木の概要に関する間伐木概要情報を作成する間伐木集計部をさらに有し、前記間伐木資源情報には、間伐木として選定された個々の樹木の番号、個々の樹木の位置情報、個々の樹高、個々の樹木の胸高直径、個々の樹木の材積が含まれており、前記間伐木概要情報には、間伐木選定前における樹木に関する情報及び間伐木選定後のおける選定間伐木に関する情報として、それぞれの本数、それぞれの所定面積当たりの本数、それぞれの平均胸高直径、それぞれの平均樹高、それぞれの合計材積、それぞれの所定面積当たりの合計材積及び選定間伐木の間伐率が含まれていることが好ましい。

本発明の間伐木の選定装置がこのような構成を有することにより、前記［９］に記載の本発明の間伐木の選定方法によって得られる効果と同様の効果が得られる。

なお、本発明の間伐の選定装置は、前記［２］～［７］に記載の本発明の間伐木の選定方法の特徴を有することが好ましい。

実施形態１に係る間伐木の選定方法における間伐木の選定処理を説明するために示すフローチャートである。 間伐対象区域における森林資源情報の一例を示す図である。 ディスプレイ上に表示された間伐対象区域を含む森林域画像を示す図である。 ディスプレイ上に表示された上層間伐による間伐対象区域画像を示す図である。 ディスプレイ上に表示された下層間伐による間伐対象区域画像を示す図である。 ディスプレイ上に表示された全層間伐による間伐対象区域画像を示す図である。 選定された間伐木に関する間伐木資源情報の一例を示す図である。 図７に示す間伐木資源情報に基づいて作成された間伐木概要情報を示す図である。 全層間伐による間伐木選定を行った結果の残存木の胸高直径（ＤＢＨ）の分布を間伐木の選定前と間伐木選定後との比較で示すヒストグラムである。 実施形態２に係る間伐木の選定方法における間伐木の選定処理を説明するためのフローチャートである。 実施形態２に係る間伐木選定方法において作成された間伐木概要情報を示す図である。 材積に基づく間伐率と本数に基づく間伐率との両方を満足する場合の間伐木概要情報の一例を示す図である。 実施形態３に係る間伐木の選定方法における間伐木の選定処理を説明するために示すフローチャートである。 ディスプレイに表示された傾斜区分画像を示す図である。 実施形態１に係る間伐木選定装置１及び実施形態２に係る間伐木選定装置２を説明するために示す図である。 実施形態３に係る間伐木選定装置３を説明するために示す図である。

以下、本発明の実施形態について説明する。実施形態に係る間伐木の選定方法は、間伐対象区域に存在する多数の樹木の個々の樹木に関する情報が記録されている森林資源情報に基づいて、当該間伐対象区域に存在する多数の樹木のうち間伐すべき樹木を選定する間伐木の選定方法である。

［実施形態１］
図１は、実施形態１に係る間伐木の選定方法における間伐木の選定処理を説明するために示すフローチャートである。

実施形態１に係る間伐木の選定方法について図１を参照して説明する。実施形態１に係る間伐木の選定方法は、図１に示すように、コンピューターなどの情報処理装置にデータベース化した状態で保存されている森林資源情報から間伐対象区域内の森林資源情報を取得する森林資源情報取得処理（ステップＳ１０）と、間伐対象区域内における間伐率（Ｒとする。）を設定する間伐率（Ｒ）設定処理（ステップＳ２０）と、間伐対象区域に存在する多数の樹木についてどのような間伐条件（詳細は後述する。）で間伐を行うかを設定する間伐条件設定処理（ステップＳ３０）と、間伐対象区域に存在する樹木の個々の樹木間の距離を算定する樹木間距離算定処理（ステップＳ４０）と、樹木間距離算定処理で算定された距離のうち、最も近接している２本の樹木を最短距離樹木として選択する最短距離樹木選択処理（ステップ５０）と、間伐条件設定処理（ステップＳ３０）によって設定されている間伐条件に基づいて、最短距離樹木選択処理（ステップＳ５０）によって選択されている２本の樹木のうちの１本の樹木を間伐木として選定する間伐木選定処理（ステップＳ６０）と、を有している。

上述の各処理のうち、樹木間距離算定処理（ステップＳ４０）から間伐木選定処理（ステップＳ６０）の処理は、間伐率設定処理（ステップＳ２０）において設定されている間伐率（設定間伐率と表記する場合もある。）Ｒに達するまで順次行う。すなわち、現時点の間伐率をＲｎとしたとき、当該現時点の間伐率Ｒｎが設定間伐率Ｒに達するまで行う。

具体的には、図１に示すように、間伐木選択処理（ステップＳ６０）を行った後に、現時点の間伐率Ｒｎが設定間伐率Ｒ未満か否か（Ｒｎ＜Ｒか否か）を判定する処理（ステップＳ７０）を行い、Ｒｎ＜Ｒであれば、現時点の処理回数ｎに１を加算する処理（ｎ＝ｎ＋１）を行って（ステップＳ８０）、樹木間距離算定処理（ステップＳ４０）に戻る。一方、現時点の間伐率Ｒｎが設定間伐率Ｒに達していれば（Ｒｎ≧Ｒ）、これまでに選定された間伐木について間伐木集計処理を行う（ステップＳ９０）。なお、実施形態１に係る間伐木の選定方法における間伐率は、間伐対象区域に存在する樹木の本数に基づく間伐率である。

また、間伐条件設定処理（ステップＳ３０）は、当該間伐対象区域においてどのような間伐を実施するかを示すものであり、「上層間伐」、「下層間伐」、「全層間伐」の中からいずれかを選択する。

ここで、「上層間伐」は、森林資源情報データベースに記録されている間伐対象区域内の各樹木の大小に基づいて、相対的に大きい樹木を間伐する間伐方法である。また、「下層間伐」は、森林資源情報データベースに記録されている間伐対象区域内の各樹木の大小に基づいて、相対的に小さい樹木を間伐する間伐方法であり、「全層間伐」は、間伐対象区域における森林資源情報に記載されている間伐対象区域内のすべての樹木から所定の規則（後述する。）に基づいて満遍なく間伐する間伐方法である。なお、ここでは、樹木の大小は、胸高直径（ＤＢＨ）の大小で表わされるものとする。すなわち、大きい樹木というのは、胸高直径（ＤＢＨ）の大きい（幹が太い）樹木であり、小さい樹木というのは、胸高直径（ＤＢＨ）の小さい（幹が細い）樹木を意味している。

なお、間伐率設定処理（ステップＳ２０）及び間伐条件設定処理（ステップＳ３０）を実施する順番は、図１においては、間伐率設定処理（ステップＳ２０）を行った後に、間伐条件設定処理（ステップＳ３０）を行うといった順番となっているが逆でもよい。すなわち、間伐件設定処理を行った後に、間伐率設定処理を行うようにしてもよい。また、間伐率設定処理（ステップＳ２０）における間伐率Ｒの設定及び間伐条件設定処理（ステップＳ３０）における間伐条件の設定は、それぞれ予め設定しておくことができる。例えば、間伐率Ｒとしては３０％を設定しておき、間伐条件としては「上層間伐」、「下層間伐」、「全層間伐」のいずれかを設定しておけば、図１に示すフローチャートにおいて、間伐率設定処理（ステップＳ２０）及び間伐条件設定処理（ステップＳ３０）は省略できる。

続いて、図１に示すフローチャートの各処理について順を追って説明する。ここでは、初期値として、設定間伐率Ｒとして３０（％）、処理回数ｎ＝１が設定されているとする。また、第１回目（ｎ＝１）の処理を開始する前の段階においては、現時点の間伐率ＲｎはＲｎ＝０である。なお、現時点における間伐率Ｒｎは、Ｒｎ＝（現時点の選定間伐木数／間伐対象区域に存在する樹木の本数）×１００で表わされる。

図１に示すフローチャートにおいて、森林資源情報取得処理（ステップＳ１０）は、既に作成されている森林資源情報を森林資源情報データベースから取得する処理であり、ここで用いる森林資源情報は公知の森林資源情報作成技術によって作成することができる。

図２は、間伐対象区域における森林資源情報の一例を示す図である。なお、図２に示す森林資源情報は、例えば、下記公知文献に示されている森林資源情報算定方法及び森林資源算定装置に基づいて作成することができる。当該公知文献は、本発明の出願人が平成２８年１１月２２日に出願したものである。

公知文献：「森林資源情報算定方法及び森林資源算定装置」、特願２０１６－２２７２０７（特開２０１８－８４４７２号公報）。

実施形態１に係る間伐木の選定方法を説明するために設定した間伐対象区域（実施形態１に係る間伐対象区域と表記する場合もある。）は、上記公知文献に記載されている「森林資源情報算定方法及び森林資源算定装置」において森林資源算定のために設定されている森林域とは異なるため、図２に示す森林資源情報は、特願２０１６－２２７２０７（特開２０１８－８４４７２号公報）に記載されている森林資源情報算定方法及び森林資源算定装置において作成された森林資源情報とは個々の樹木に関する情報が異なったものとなっている。

実施形態１に係る間伐対象区域において作成された森林資源情報には、図２に示すように、個々の樹木の番号（ナンバー）、各樹木（各樹木の樹冠の中心位置）の位置情報（ｘ座標及びｙ座標）、樹種（カラマツ）、樹高（ｍ）、胸高直径（ｃｍ）、材積（ｍ^３）が個々の樹木に対応して示されている。なお、上記公知文献には、樹種については記載されていないが、樹種の分類については、前述の「先行技術文献」として用いた特許文献１（特許４９００３５６号公報）に記載されている。

また、図２においては、間伐対象区域に存在する樹木の一部のみが示されているが、実際には、百本単位、千本単位の樹木が存在する場合もあり、又はそれ以上の樹木が存在する場合もある。なお、「胸高直径」は、前述したように、樹木の幹の太さを表すものであり、「ＤＢＨ」と表記する場合もある。

図２に示すような間伐対象区域における森林資源情報を森林域資源情報データベースから取得して（ステップＳ１０）、間伐率設定処理（ステップＳ２０）以降の処理を行う。間伐率Ｒ設定処理（ステップＳ２０）においては、間伐率Ｒ（設定間伐率Ｒ）がＲ＝３０％と設定されたとする。

続いて、樹木間距離算定処理（ステップＳ４０）を行う。樹木間距離算定処理は、間伐対象区域に存在するすべての樹木において樹木間距離を算定し、最も近接している２本の樹木を最短距離樹木として選択する処理である。例えば、間伐対処区域に存在する樹木が１００本であったとすると、当該１００本の樹木間でそれぞれ樹木間距離を算定する。このときの樹木間距離は、例えば、ピタゴラスの定理を用いて行うことができる。すなわち、図１における破線枠Ｓ４０’内に示すように、樹木間距離Ｄｎ（ｃｍ）は、
Ｄｎ^２＝（（Ｘｎ＋１）－Ｘｎ）^２＋（（Ｙｎ＋１）－Ｙｎ）^２ ・・・（１）
から求めることができる。

なお、（１）式において、Ｘｎは、処理回数が第ｎ回目（例えば、ｎ＝１とする。）のときの樹木間距離算定対象となる２本のうちの一方の樹木のｘ座標、Ｘｎ＋１は当該処理回数（ｎ＝１）のときの樹木間距離算定対象となる２本のうちの他方の樹木のｘ座標、Ｙｎは、当該処理回数（ｎ＝１）のときの樹木間距離算定対象となる２本のうちの一方の樹木のｙ座標、Ｙｎ＋１は当該処理回数（ｎ＝１）のときの樹木間距離算定対象となる２本のうちの他方の樹木のｙ座標である。

ここで、間伐対象区域における樹木の総数を例えば１００本とした場合、第１回目の処理（ｎ＝１）においては、１００本の樹木それぞれについて（１）式の計算を行って、Ｄｎが最小（最短距離）の２本の樹木を選択する（ステップＳ５０）。そして、当該処理回数（ｎ＝１）において選択された最短距離の２本の樹木のうちの１本の樹木を、間伐条件設定処理（ステップＳ３０）において設定されている間伐条件に基づいて間伐木として選定する（ステップＳ６０）。

間伐木選定処理（ステップＳ６０）は、ステップＳ３０において設定されている間伐条件に基づいて行われる。すなわち、図１における破線枠Ｓ６０’内に示すように、間伐条件設定処理（ステップＳ３０）において「上層間伐」が設定されている場合には、処理Ａ１を行う。処理Ａ１は、最短距離の２本の樹木のうち、胸高直径（ＤＢＨ）が大きい方の樹木（ＤＢＨｍａｘ）を間伐木として選定する。また、「下層間伐」が設定されている場合には処理Ａ２を行う。処理Ａ２は、最短距離の２本の樹木のうち、胸高直径（ＤＢＨ）が小さい方の樹木（ＤＢＨｍｉｎ）を間伐木として選定する。また、「全層間伐」が設定されている場合には処理Ａ３を行う。

処理Ａ３は、最短距離樹木として選択された２本の樹木のうちの１本を所定の規則に基づいて選択する。ここで、「所定の規則」は、現時点の処理回数ｎを表す数値が奇数か偶数かによって、最短距離樹木選択処理（ステップＳ５０）によって選択されている２本の樹木のうちの胸高直径（ＤＢＨ）が大きい方の樹木（ＤＢＨｍａｘ）の樹木又は胸高直径（ＤＢＨ）が小さい方の樹木（ＤＢＨｍｉｎ）のいずれかを間伐木として選択することである。

具体的には、現時点の処理回数ｎが偶数である場合には、最短距離樹木として選択された２本の樹木のうち胸高直径が大きい方の樹木（ＤＢＨｍａｘ）を間伐木として選定し、現時点の処理回数ｎが奇数であれば、胸高直径が小さい方の樹木（ＤＢＨｍｉｎ）を間伐木として選定する。このような処理は、モジュロ演算によって得られた余りに基づいて行うことができる。すなわち、ｎ÷２を行ったときの余りが「０」であれば偶数、ｎ÷２を行ったときの余りが「１」であれば奇数とする。なお、第１回（ｎ＝１）の場合は、奇数として設定しておけばよい。

間伐木選定処理（ステップＳ４０）により、現時点の処理回数（例えば、ｎ＝１）において、１本の間伐木が選定されたら、当該処理回数（ｎ＝１）における間伐率ＲｎがステップＳ２０において設定されている間伐率Ｒ（設定間伐率Ｒ）未満であるか否か（Ｒｎ＜Ｒであるか否か）を判定し（ステップＳ７０）、Ｒｎ＜Ｒであれば（設定間伐率Ｒ未満であれば）、ｎ＝ｎ＋１（ステップＳ８０）を行って、樹木間距離算定処理（ステップＳ４０）に戻って、当該樹木間距離算定処理（ステップＳ４０）以降の処理（第２回目以降の処理）を行う。

第２回目の処理（ｎ＝２）においては、残りの９９本の樹木それぞれについて（１）式により樹木間距離算定処理（ステップＳ４０）を行い、９９本の樹木のうちＤｎが最小の２本の樹木（最短距離樹木）を選択する処理（ステップＳ５０）を行う。そして、当該処理回数（ｎ＝２）において最短距離樹木とされた２本の樹木のうちの１本の樹木を間伐条件設定処理（ステップＳ３０）において設定されている間伐条件に基づいて間伐木として選定する。

このような処理をステップＳ２０において設定した間伐率Ｒ（設定間伐率Ｒ）に達するまで行う。一方、Ｒｎ＜Ｒであるか否かの判定処理（ステップＳ７０）において、Ｒｎ＜Ｒでなければ（設定間伐率Ｒ（Ｒ＝３０％）に達していれば）、間伐木集計処理を行う（ステップＳ９０）。なお、間伐木集計処理を行う際には、現時点の間伐率Ｒｎが設定間伐率Ｒに達した直前の樹木までを、間伐木として選定された樹木として間伐木集計を行う。

ところで、ステップＳ７０における「Ｒｎ＜Ｒ」の判定処理は、例えば、間伐対象区域に存在する樹木の総数が３００本である場合には、当該処理回数ｎがｎ＝１（第１回）のときは、選定された間伐木の本数は１本であるため、現時点（第１回）における間伐率Ｒｎは、Ｒｎ＝（１÷３００）×１００≒０．３３％であり、また、当該処理回数ｎがｎ＝２のときは、選定された間伐木の本数は２本であるため、現時点（第２回）における間伐率Ｒｎは、Ｒｎ＝（２÷３００）×１００≒０．６７％である。

また、当該処理回数ｎがｎ＝９０のときは、選定された間伐木の本数は９０本であるため、現時点（第９０回）における間伐率Ｒｎは、Ｒｎ＝（９０÷３００）×１００＝３０．０％である。このため、当該処理回数ｎがｎ＝９０となると、Ｒ（ｎ）＜Ｒを満たさなくなる。従って、樹木総数が例えば３００本の場合においては、当該処理回数ｎがｎ＝８９までは、Ｒｎ＜Ｒを満たすこととなり、ステップＳ４０に処理が戻る。一方、当該処理回数ｎがｎ＝９０となると、現時点の間伐率Ｒｎは３０％となり、Ｒｎ＜Ｒを満たさなくなるため、間伐木集計処理（ステップＳ９０））に移る。

以上、図１に示すフローチャートを参照して実施形態１に係る間伐木選定方法の全体的な処理を概略的に説明したが、続いて、実施形態１に係る間伐木選定方法について、図１のフローチャート及び図２～図１２を参照して具体的に説明する。

図３は、ディスプレイ上に表示された間伐対象区域を含む森林域画像を示す図である。図３（ａ）は間伐対象区域を含む森林域を上空からドローンなどによって撮影して得られた撮像画像をディスプレイ上に表示した画像（森林域画像）であり、図３（ｂ）は図３（ａ）に示す森林域の画像から間伐対象区域及びその周辺の画像を拡大してディスプレイ上に表示した画像である。

図３（ａ）及び図３（ｂ）において、黒枠線Ｌ１で囲った区域が間伐対象区域（定性間伐区ともいう。）であるとする。図３（ｂ）においては、１つ１つの樹木に対応する樹冠が樹木ごとの樹冠として判別することができる。なお、図３はモノクロ画像であるが、図３の元となるカラー画像上では、各樹木の樹冠をさらに容易に識別することができる。

ここで、図３（ｂ）において、個々の樹冠に示されている白抜き丸印（なお、正確には丸印ではないが丸印と表記する。）は、当該白抜き丸印の大きさが胸高直径（ＤＢＨ）の大きさ（ｃｍ）を表している。ここでは、ＤＢＨ＝３０．０ｃｍ以下、ＤＢＨ＝３０．１ｃｍ～３５．０ｃｍ、ＤＢＨ＝３５．１ｃｍ～４０．０ｃｍ、ＤＢＨ＝４０．１ｃｍ～４５．０ｃｍ、ＤＢＨ＝４５．１ｃｍ以上の５段階に分けている。なお、各樹木の胸高直径（ＤＢＨ）は、図２に示した森林資源情報を参照することによって取得することができる。

続いて、実施形態１に係る間伐木選定方法における具体的な間伐木選定方法について説明する。ここで、具体的な間伐木選定方法として、「上層間伐」による間伐木選定処理、「下層間伐」による間伐木選定処理、「全層間伐」による間伐木選定処理を例示して説明する。なお、これらの３つの間伐木選定処理においては、いずれも、間伐率設定処理（ステップＳ２０）においては、間伐率Ｒとして３０％が設定されているとする。

［上層間伐による間伐木選定処理］
この場合、間伐条件設定処理（ステップＳ３０）においては、「上層間伐」が設定されている。間伐条件として「上層間伐」が設定された後に、樹木間距離算定処理を行う（ステップＳ４０）。この樹木間距離算定処理においては、処理回数が第１回目（ｎ＝１）の処理として、間伐対象区域に存在する樹木すべてにおいて、前述の（１）式によって、それぞれ２本の樹木間の樹木間距離Ｄｎを計算し（ステップＳ４０）、Ｄｎが最小の２本の樹木（最短距離樹木）を選択する（ステップＳ６０）。

そして、当該第１回目（ｎ＝１）の処理において最短距離樹木とされた２本の樹木から１本の樹木を所定の間伐条件（この場合、「上層間伐」）に基づいて１本の樹木を間伐木として選定する（ステップＳ６０）。具体的には、間伐方法として、「上層間伐」が設定されている場合には、図１に示すフローチャートのステップＳ６０における破線枠Ｓ６０’に示すように、処理Ａ１が行われる。すなわち、最短距離として選択された２本の樹木のうち、ＤＢＨが大きい方（ＤＢＨｍａｘ）の樹木を間伐木として選定する。

ところで、樹木間距離算定処理（ステップＳ４０）は、例えば、樹木総数が３００本であったとすると、第１回目（ｎ＝１）においては、３００本の樹木の個々の樹木が（１）式を自分の樹木以外の２９９本の樹木に対して行うこととなるため、第１回目（ｎ＝１）においては、（１）式を３００×（３００－１）＝８９７００回計算することとなる。なお、３００×（３００－１）＝８９７００は、２本の樹木の組み合わせにおいて相互に距離計算がなされることとなるため、２本の樹木の組み合わせにおいて同じ距離計算が重複して行われることとなる。ここで、重複する距離計算を省くとしても、初回（ｎ＝１）における（１）式の距離計算回数は、（３００×（３００－１））÷２＝４４８５０回となる。

なお、間伐対象区域に存在する樹木は、実際には、数千本以上にも達する場合もあるため、樹木間距離算定処理（ステップＳ４０）は、ディープランニング（深層学習）のニューラルネットワークを用いることが好ましい。また、樹木間距離算定処理（ステップＳ４０）と、当該樹木間距離算定処理によって得られた各樹木間距離Ｄｎに基づいて、Ｄｎが最小の２本の樹木（最短距離樹木）を選択する処理（ステップＳ５０）までをディープランニング（深層学習）のニューラルネットワークを用いるようにしてもよい。これにより、間伐対象区域に存在する樹木が、数千本存在したとしても、最短距離の２本の樹木を短時間で、かつ、高精度に選択することができる。

上述の間伐木選定処理（ステップＳ６０）が終わると、この時点（第１回目（ｎ＝１）の間伐率Ｒｎが設定間伐率Ｒ（この場合、Ｒ＝３０％）に達しているか否か（Ｒｎ＜Ｒ）を判定する（ステップＳ７０）。この時点では、間伐木の選定は１本であるため、この時点における間伐率Ｒ（ｎ）は、Ｒｎ＜Ｒである。なお、この時点の間伐率Ｒｎというのは、間伐対象区域に存在する樹木の総数が、仮に３００本であるとすれば、第１回（ｎ＝１）における間伐率Ｒｎ（ｎ＝１）は、（１÷３００）×１００≒０．３３％となる。

このようにして、第１回目（ｎ＝１）の処理において、間伐すべき１本の樹木が選定されたこととなり、かつ、現時点における間伐率が設定間伐率（３０％）に達していないため、処理回数ｎに１を加算する処理（ｎ＝ｎ＋１）を行い（ステップＳ８０）、ステップＳ４０（樹木間距離算定処理）に戻って、第２回目（ｎ＝２）の処理を行う。当該第２回目（ｎ＝２）の処理においては、残りの２９９本の樹木それぞれについて、上述同様に（１）式を行って、２９９本の樹木の中からＤｎが最小の２本の樹木（最短距離樹木）を選択する（ステップＳ５０）。

そして、当該第２回目（ｎ＝２）の処理において選択されたＤｎが最小の２本の樹木（最短距離樹木）について間伐木の選定処理（ステップＳ６０）を行う。この場合も、既に、ステップＳ３０において、間伐方法として「上層間伐」が設定されているため、第１回目（ｎ＝１）の処理と同様に２本の樹木のうち、ＤＢＨの大きい方の樹木を間伐木として選定する。これによって、ここまでに２本の間伐木が選定されたこととなる。

そして、現時点の間伐率Ｒｎ、すなわち、第２回目（ｎ＝２）の間伐率Ｒｎが設定間伐率Ｒ（この場合、Ｒ＝３０％）未満であるか否か（Ｒｎ＜Ｒであるか否か）を判定する（ステップＳ７０）。現時点では、間伐木の選定は２本であるため、間伐対象区域に存在する樹木の総数が、仮に３００本であるとすれば、第２回（ｎ＝２）における間伐率Ｒｎ（ｎ＝１）は、（２÷３００）×１００≒０．６７％となる。

このようにして、第２回目（ｎ＝１）までの処理において、間伐すべき２本の樹木が選定されたこととなり、かつ、現時点における間伐率が設定間伐率（３０％）に達していないため、処理回数ｎに１を加算する処理（ｎ＝ｎ＋１）を行い、ステップＳ４０（樹木間距離算定処理）に戻って、第３回目（ｎ＝３）の処理を行う。当該第３回目（ｎ＝３）の処理においては、残りの２９８本の樹木それぞれについて、上述同様に（１）式の計算を行って、２９９本の樹木の中からＤｎが最小の２本の樹木（最短距離樹木）を選択する（ステップＳ５０）。

このように、ステップＳ２０において設定されている間伐率（設定間伐率）Ｒと、ステップＳ３０において設定されている間伐方法（この場合「上層間伐」）とに基づいて、ステップＳ４０～Ｓ７０を繰り返し行う。そして、ステップＳ７０において、現時点の間伐率Ｒｎが設定間伐率（Ｒ＝３０％）未満であるか否かの判定（Ｒｎ＜Ｒ）を行った結果、Ｒｎ＜Ｒを満たさなくなると、これまでに選定した間伐木に関する資源情報を集計する間伐木集計処理（ステップＳ９０）を行う。

例えば、間伐対象区域に存在する樹木の総数が３００本であって、設定間伐率が３０％であるとすれば、Ｒｎ＜Ｒを満たさなくなる直前の本数（間伐率が３０％未満の最大の本数）は、８９本であるため、選定された８９本の樹木に関する間伐木資源情報が集計される。

図４は、ディスプレイ上に表示された上層間伐による間伐対象区域画像を示す図である。なお、「上層間伐による間伐対象区域画像」というのは、「上層間伐」の間伐条件に基づいて間伐木の選定がなされた後の間伐対象区域画像である。また、図４に示す間伐対象区域画像は図３（ｂ）において示した間伐対象区域と同じ森林域であり、黒枠線Ｌ１で囲まれた範囲が間伐対象区域である。

また、図４に示す各樹木に対応する樹冠に示されている丸印は、図３（ｂ）において示した白抜き丸印に対応するものであり、図４においては、「上層間伐」の間伐条件に基づいて選定された間伐木については、丸印が塗りつぶされている。また、図４においては図３（ｂ）と同様に、丸印の大きさがＤＢＨの大きさを表している。また、図４において、白抜き丸印の樹木は、「上層間伐」による間伐条件により間伐木として選定されなかった残存木を表している。

なお、図４はモノクロ画像であるため、選定された間伐木の判別がしにくいが、図４の元となるカラー画像上では、選定された間伐木の丸印は赤色で示されているため、選定された間伐木を容易に判別することができる。なお、図４においては図３（ｂ）と同様に、丸印の大きさがＤＢＨの大きさを表している。

「上層間伐」により間伐においては、最短距離として選定された２本の樹木のうち、胸高直径（ＤＢＨ）の大きい方の樹木を間伐木として選定しているため、図４に示すように、全体的に、ＤＢＨの大きい樹木が間伐木として選定されていることが分かる。

このように「上層間伐」を設定した場合の間伐においては、全体的に、ＤＢＨの大きい樹木が間伐木として選定されることとなるため、間伐計画の立案者は、木材価値の高い樹木から自動選択することができ、高い木材収入を得ることができる。また間伐後の中層木と下層木の残存木が均等配置されることで、広い生育空間があることから旺盛な生育が期待できる森林となる。

［下層間伐による間伐木選定処理］
「下層間伐」による間伐木選定処理における基本的な処理フローは、図１に示すフローチャートと同じである。「下層間伐」による間伐木選定処理が「上層間伐」による間伐木選定処理と異なる主な点は、図１に示す間伐方法設定処理（ステップＳ３０）において「下層間伐」が設定されている点及び間伐木選定処理（ステップＳ６０）において最短距離として選択された２本の樹木のうち、ＤＢＨの小さい方（ＤＢＨｍｉｎ）の樹木を間伐木として選定する点である。これら以外の各ステップにおける処理は、実施形態１に係る間伐木の選定方法と同じであるため、説明は省略する。

図５は、ディスプレイ上に表示された下層間伐による間伐対象区域画像を示す図である。なお、「下層間伐による間伐対象区域画像」というのは、「下層間伐」の間伐条件に基づいて間伐木の選定がなされた後の間伐対象区域画像である。また、図５は図３（ａ）において示した間伐対象区域と同じ森林域であり、黒枠線Ｌ１で囲まれた範囲が間伐対象区域である。

また、図５に示す各樹木に対応する樹冠に示されている丸印は、図３（ｂ）において示した丸印（白抜き丸印）に対応するものであり、図５において、当該「下層間伐」の間伐方法に基づいて選定された間伐木については、丸印が塗りつぶされている。図５においては図３（ｂ）と同様に、丸印の大きさがＤＢＨの大きさを表している。また、図５において、白抜き丸印の樹木は、「下層間伐」により間伐木として選定されなかった残存木を表している。

なお、図５はモノクロ画像であるため、選定された間伐木の判別がしにくいが、図５の元となるカラー画像上では、選定された間伐木の丸印は赤色で示されているため、選定された間伐木を容易に判別することができる。

「下層間伐」による間伐においては、最短距離として選定された２本の樹木のうち、胸高直径（ＤＢＨ）の小さい方の樹木を間伐木として選定しているため、図５に示すように、全体的に、ＤＢＨの小さい樹木が間伐木として選定されていることが分かる。

このように、「下層間伐」を設定した場合の間伐においては、全体的に、ＤＢＨの小さい樹木が間伐木として選定されることとなるため、間伐計画の立案者は、間伐した後の木材価値が向上する森林になるように、生育の良い樹木（例えば、ＤＢＨの大きい樹木）を残存木として均等に残し、生育の障害となる混んだ箇所の木材価値の低い下層木（ＤＢＨの小さい樹木）を間伐することができる。

［全層間伐による間伐木選定処理］
「全層間伐」による間伐木選定処理における基本的な処理フローも、図１に示すフローチャートと同じである。「全層間伐」による間伐木選定処理が「上層間伐」及び「下層間伐」による間伐木選定処理実と異なる主な点は、図１に示す間伐条件設定処理（ステップＳ３０）において、「全層間伐」が設定されている点及び間伐木の選定処理（ステップＳ６０）において、最短距離として選択された２本の樹木のうち１本を間伐木として選定する際に、最短距離樹木として選択された２本の樹木のうちの１本を所定の規則に基づいて選択する点である。

すなわち、間伐木選定処理（ステップＳ６０）は、現時点の処理回数ｎが偶数であれば、最短距離樹木として選択された２本の樹木のうち、胸高直径が大きい方の樹木（ＤＢＨｍａｘ）を間伐木として選定し、現時点の処理回数ｎが奇数であれば、胸高直径が小さい方の樹木（ＤＢＨｍｉｎ）を間伐木として選定する。これら以外の各ステップにおける処理は実施形態１に係る間伐木の選定方法と同じであるため、説明は省略する。

このように、「全層間伐」の場合においては、最短距離樹木として選択された２本の樹木のうちの１本を選択する際に、現時点の処理回数（ｎ）が奇数か偶数かによって、ＤＢＨの大小を交互に切り替えて間伐木を選定するようにしている。このため、「全層間伐」による間伐木の選定は、全層から満遍なく行うことができる。

図６は、ディスプレイ上に表示された「全層間伐」による間伐対象区域画像を示す図である。なお、「全層間伐による間伐対象区域画像」というのは、「全層間伐」による間伐木の選定がなされた後の間伐対象区域画像である。また、図６は図３（ｂ）において示した間伐対象区域と同じ森林域であり、黒枠線Ｌ１で囲まれた範囲が間伐対象区域である。

また、図６に示す各樹木に対応する樹冠に示されている丸印は、図３（ｂ）において示した丸印（白抜き丸印）に対応するものであり、図６においては、「全層間伐」の間伐方法に基づいて選定された間伐木に付されている丸印が塗りつぶされている。図６においては図３（ｂ）と同様に、丸印の大きさがＤＢＨの大きさを表している。また、図６において、白抜き丸印の樹木は、「全層間伐」により間伐木として選定されなかった残存木を表している。

なお、図６はモノクロ画像であるため、選定された間伐木の判別がしにくいが、図６の元となるカラー画像上では、選定された間伐木の丸印は赤色で示されているため、選定された間伐木を容易に判別することができる。

「全層間伐」により間伐においては、最短距離として選定された２本の樹木のうち、胸高直径（ＤＢＨ）の小さい方の樹木と大きい方の樹木とが交互に間伐木として選定されるため、図６に示すように、全体的に、胸高直径（ＤＢＨ）の小さい方の樹木と大きい方の樹木とが満遍なく間伐木として選定されていることが分かる。

このように、「全層間伐」を設定した場合の間伐においては、胸高直径（ＤＢＨ）の小さい方の樹木と大きい方の樹木とが満遍なく間伐木として選定されているため、間伐計画の立案者は、現在の上層木、中層木、下層木が配置されている森林状態の構成を持続的に維持できる森林になるように、胸高直径の大小に関わらず、生育の障害となる混んだ箇所の樹木を間伐することができる。これにより、全層からの木材収入と森林の価値を長期に渡り継続的に得ることができる。

ところで、上述した３種類の間伐木の選定処理、すなわち、「上層間伐」による間伐木選定処理、「下層間伐」による間伐木の選定処理、「全層間伐」による間伐木の選定処理のそれぞれにおいて、現時点の間伐率Ｒｎが設定間伐率（Ｒ＝３０％）未満であるか否かの判定Ｒｎ＜Ｒを行った結果、Ｒｎ＜Ｒを満たさなくなると、間伐木集計処理（ステップＳ９０）にて、間伐木集計処理を行う。この間伐木集計処理には、選定されたすべての間伐木の個々の間伐木に関する間伐木資源情報を作成する処理と、当該間伐木資源情報に基づいて間伐木概要情報を作成する処理とが含まれている。
続いて、この間伐木集計処理（ステップＳ９０）について説明する。ここでは、「全層間伐」による間伐木集計処理について説明する。

図７は、間伐木集計処理によって作成された間伐木資源情報の一例を示す図である。図７に示す間伐木資源情報は、図２に示す間伐対象区域に存在する樹木に関する森林資源情報の中から「全層間伐」によって選定された間伐木に関する間伐木資源情報の一部である。当該間伐木資源情報は、データベース化されてパーソナルコンピューターなどの情報処理装置の記憶部に保存される。なお、当該間伐木資源情報には、間伐木として選定された樹木のナンバー、各樹木の樹冠位置（樹冠の中心位置）のｘ座標及びｙ座標、樹種、樹高（ｍ）、胸高直径（ＤＢＨ：ｃｍ）、材積（ｍ^３）が記録される。

図８は、図７に示す間伐木資源情報に基づいて作成された間伐木概要情報を示す図である。図８に示す間伐木概要情報は、０．５４ｈａ（正確には５３５６ｍ^２＝０．５３５６ｈａであるとする。）の間伐対象区域についての間伐木概要情報を示している。当該間伐木概要情報には、当該間伐対象区域における間伐木選定前及び間伐木選定後についての樹木に関する情報として、「本数」、「所定面積（１ｈａとする。）当たりの本数」、「平均ＤＢＨ（ｃｍ）」、「平均樹高（ｍ）」、「合計材積（ｍ^３）」、「所定面積（１ｈａとする。）当たりの合計材積（ｍ^３）」、「間伐率」が記録されている。なお、「本数」は「間伐木選定前」においては間伐対象区域に存在する樹木の本数（総数）であり、「間伐木選定後」においては選定された間伐木数（選定間伐木数）である。

図８に示す間伐木概要情報によれば、「間伐率」は、選定間伐木数（６３本とする。）÷総数（２１１本とする。）×１００≒２９．８６％となっており、設定間伐率Ｒ（３０％）未満に収まっている。

なお、１ｈａ当たりの本数のうち、間伐木選定前の１ｈａ当たりの本数（３９４本）は、「２１１÷０．５３５６≒３９４」によって得られた値であり、間伐木選定後の１ｈａ当たりの本数（１１８本）は、「６３÷０．５３５６≒１１８」によって得られた値である。

また、１ｈａ当たりの合計材積（ｍ^３）のうち、間伐木選定前の１ｈａ当たりの合計材積（４７０．７ｍ^３）は、２５２．１÷０．５３５６≒４７０．７によって得られた値であり、間伐木選定後の１ｈａの合計材積（１４７．３ｍ^３）は、８７．５÷０．５３５６≒１４７．３によって得られた値である。

なお、図７に示す間伐木資源情報及び図８に示す間伐木概要情報はそれぞれデータベース化されてパーソナルコンピューターなどの情報処理装置に保存され、随時、ディスプレイ上に表示させることができる。

なお、ここでは、「全層間伐」による間伐木選定処理についての間伐木資源情報及び間伐木概要情報について説明したが、「上層間伐」及び「下層間伐」についても、数値は異なるが形式としては図７及び図８と同様の形式の間伐木資源情報及び間伐木概要情報を作成することができるが、これらの説明及び図示は省略する。

なお、前述した３種類の間伐方法（「上層間伐」、「下層間伐」及び「全層間伐」）は、それぞれの間伐方法固有の効果が得られる。例えば、「上層間伐」は、全体的にＤＢＨの大きい樹木が間伐木として選定されることとなるため、木材価値の高い樹木を間伐木として選定でき、高い木材収入を得ることができる。また、間伐後に残された中層木及び下層木の残存木がそれぞれ所定の距離を置いて配置されることとなるため、これら中層木及び下層木は広い生育空間が得られることから旺盛な生育が期待できる森林となる。また、「下層間伐」は、全体的にＤＢＨの小さい樹木が間伐木として選定されることとなる。すなわち、生育の障害となる混んだ箇所に生えている木材価値の低い下層木（ＤＢＨの小さい樹木）が間伐木として選定されることとなるため、生育の良い樹木（例えば、ＤＢＨの大きい樹木）が残存木としてそれぞれ所定距離を置いて配置されることとなり、間伐した後の木材価値の向上が期待できる森林になる。

また、「全層間伐」は、全体的に胸高直径（ＤＢＨ）の小さい方の樹木と大きい方の樹木とが満遍なく間伐木として選定されるため、胸高直径の大小に関わらず、生育の障害となる混んだ箇所の樹木を間伐することができ、現在の上層木、中層木、下層木が配置されている森林状態の構成を持続的に維持しながら、生育の障害となる混んだ箇所の樹木を間伐することができる。これにより、全層（上層、中層及び下層）からの木材収入と森林の価値を長期に渡り継続的に得ることができる。

特に、「全層間伐」による間伐方法においては、上述したように、現在の上層木、中層木、下層木が配置されている森林状態の構成を持続的に維持できることが大きな特徴でもある。

図９は、全層間伐による間伐木選定を行った結果の残存木の胸高直径（ＤＢＨ）の分布を間伐木の選定前と間伐木選定後との比較で示すヒストグラムである。図９（ａ）は間伐木選定前における残存木の胸高直径（ＤＢＨ）の分布を示すヒストグラムであり、図９（ｂ）は間伐木選定後における残存木の胸高直径（ＤＢＨ）の分布を示すヒストグラムである。

図９に示すように、図２に示す間伐対象区域内において、「全層間伐」による間伐木の選定を行ったところ、間伐木選定後のＤＢＨのヒストグラムは、間伐木選定前のＤＢＨのヒストグラムに比べて、ヒストグラムの形を大きく変えることなく、樹木の本数を減らす（間伐する）ことができることがわかる。このように、「全層間伐」による間伐木の選定は、現在の上層木、中層木、下層木が配置されている森林状態の構成を持続的に維持しながら、生育の障害となる混んだ箇所の樹木を間伐することができる間伐木の設定方法であると言える。

［実施形態２］
上述の実施形態１に係る間伐木の選定方法においては、間伐率設定処理（ステップＳ２０）における間伐率Ｒは、樹木の本数に基づく間伐率を設定したが、実施形態２に係る間伐木の選定方法においては、樹木の材積に基づく間伐率を設定する。

図１０は、実施形態２に係る間伐木の選定方法における間伐木の選定処理を説明するためのフローチャートである。図１０に示すフローチャートは、全体の処理の流れとしては、図１に示すフローチャートとほぼ同様である。

図１０に示すフローチャートが図１に示すフローチャートと異なる点は、間伐率設定処理（ステップＳ２０）において、間伐率Ｒを樹木の材積に基づく間伐率で設定する点である。ここでは、間伐率Ｒとして３５％を設定するものとする。すなわち、間伐率設定処理（ステップＳ２０）においては、間伐率Ｒ（Ｒ＝３５（％））を設定する。なお、ここでも、間伐率Ｒを「設定間伐率Ｒ」と表記する場合もある。

そして、ステップＳ４０～Ｓ６０においては、上述の実施形態１に係る間伐木の選定方法と同様の処理を行い、最短距離の２本の樹木について、ステップＳ３０で設定した間伐条件に基づいて２本のうちの１本を間伐木として選定する処理を行う。なお、ステップＳ６０の間伐木の選定処理は、実施形態２に係る間伐木選定方法においても、実施形態１に係る間伐木選定方法と同様に、間伐条件として「上層間伐」、「下層間伐」、「全層間伐」のいずれかが設定された状態で、間伐木の選定処理を行う。

その後、現時点における材積による間伐率Ｒｎが、設定間伐率Ｒ（Ｒ＝３５％）より小さいか否か（Ｒｎ＜Ｒであるか否か）を判定し（ステップＳ７０）、Ｒｎ＜Ｒであれば、処理回数ｎに１を加算（ｎ＋１）する処理を行って（ステップＳ８０）、ステップＳ４０に戻り、当該ステップＳ４０以降の処理を行う。なお、この場合、現時点の材積による間伐率Ｒｎは、Ｒｎ＝（現時点の選定間伐木合計材積／間伐対象区域に存在する樹木合計材積）×１００で表わされる。

なお、ステップＳ７０において「Ｒｎ＜Ｒであるか否か」を判定する処理を行う際のＲｎは、現時点の処理回数ｎが例えば第１回目（ｎ＝１）であって、当該第１回目（ｎ＝１）の処理において選定された１本の樹木の材積が例えば１．４５ｍ^３であり、間伐対象区域に存在する樹木合計材積が例えば２５２．１ｍ^３であったとすると、現時点の処理回数ｎが第１回目（ｎ＝１）の処理における間伐率Ｒｎ（ｎ＝１）は、（１．４５÷２５２．１）×１００≒０．５７％となる。また、現時点の処理回数ｎが例えば第４０回目（ｎ＝４０）であって、当該第４０回目までにおいて選定された４０本の樹木の合計材積が例えば６０．７ｍ^３であったとすると、現時点の処理回数ｎ（ｎ＝４０）における間伐率Ｒｎ（ｎ＝４０）は、（６０．７÷２５２．１）×１００≒２４．１％となる。なお、各樹木の材積は、森林資源情報（間伐対象区域における森林資源情報）から取得できる。

このような処理を順次行い、Ｒｎ＜Ｒ（＝３５％）を満たさなくなると（Ｒｎ≧Ｒとなると）、間伐木資源情報集計処理を行う。例えば、第６３回目までにおいて選定された６３本の樹木の合計材積が例えば８７．５ｍ^３であったとすると、現時点の処理回数ｎ（ｎ＝６３）における間伐率Ｒｎ（ｎ＝６３）は、（８７．５÷２５２．１）×１００≒３４．７％となり、設定間伐率Ｒ（Ｒ＝３５％）直前の間伐率となる。そして、次の処理（第６４回目）の処理においける合計材積が設定間伐率Ｒ（３５％）に達したとする。この場合も、図７に示すような間伐木資源情報を得ることができ、当該間伐木資源情報に基づいて間伐木概要情報を作成することができる。

図１１は、実施形態２に係る間伐木選定方法において作成された間伐木概要情報を示す図である。図１１に示す間伐木概要情報は、図８に示した実施形態１に係る間伐木概要情報と基本的には同じ内容であるが、材積における間伐率が記録されている点が図８と異なる。図１１に示す間伐木概要情報によれば、間伐率は、間伐木選定後の合計材積（８７．５ｍ^３）÷間伐対象区域の樹木総材積（２５２．１ｍ^３）×１００＝３４．７％となっており、設定間伐率Ｒ（Ｒ＝３５％）未満に収まっている。

以上説明したように、実施形態２に係る間伐木選定方法によれば、材積による間伐率を設定して、設定した材積による間伐率（Ｒ＝３５％）未満に抑えた間伐率となるように間伐木選定を行うことができる。なお、材積に基づく間伐率と本数に基づく間伐率とを組み合わせるようにしてもよい。この場合、材積に基づく間伐率と本数に基づく間伐率との両方を満足するように間伐木を選定するようにして、図１２に示すような間伐木概要情報を作成するようにしてもよい。

図１２は、材積に基づく間伐率と本数に基づく間伐率との両方を満足する場合の間伐木概要情報の一例を示す図である。図１２によれば、実施形態１に係る間伐木の選定方法において得られた本数に基づく間伐率（２９．８６％）と、実施形態２に係る間伐木の選定方法において得られた材積による間伐率（３４．７％）とがそれぞれ記録されている。なお、図１２に示す間伐木概要情報においては、本数に基づく間伐率及び材積による間伐率がともに、設定間伐率Ｒ未満となっているが、材積に基づく間伐率が、設定間伐率Ｒに収まっていればよい場合もあり、逆に、材積による間伐率が設定間伐率Ｒ未満であればよい場合もある。

ここでは、「全層間伐」についての間伐木資源情報及び間伐木集計表について説明したが、「上層間伐」による間伐処理及び「下層間伐」による間伐処理を行った場合にも、数値は異なるが、同じ形式の間伐木資源情報及び間伐木集計表を得ることができる。このため、「上層間伐」による間伐処理及び「下層間伐」による間伐処理を行った場合に作成される間伐木資源情報及び間伐木集計表について図示は省略する。

［実施形態３］
図１３は、実施形態３に係る間伐木の選定方法における間伐木の選定処理を説明するために示すフローチャートである。実施形態３に係る間伐木の選定方法における間伐木の選定処理が実施形態１に係る間伐木の選定方法における間伐木の選定処理と異なるのは、間伐対象区域の傾斜角度を判定して、間伐対象区域内において所定の傾斜角度未満の区域の画像データを作成する処理が含まれている点である。このように、間伐対象区域内において所定の傾斜角度未満の区域の画像データを作成する処理を行うのは、主に、林業用ハーベスタが走行可能であるか否かを判断するために行われるものである。なお、林業用ハーベスタが走行可能である傾斜角度をここでは「間伐可能傾斜角度」と呼ぶことにする。

図１３に示すフローチャートが図１に示すフローチャートと異なる点は、間伐率設定処理（ステップＳ２０）の前に、間伐対象区域を複数に区分した各区分ごとの傾斜角度を表す傾斜区分画像データを作成する傾斜区分画像データ作成処理（ステップＳ１１）と、傾斜区分画像データに基づいて所定の傾斜角度未満の傾斜を有する区分を間伐可能区分とした間伐可能区分画像データを作成する間伐可能区分画像データ作成処理（ステップ１２）と、間伐可能区分画像データと森林資源情報とに基づいて、間伐可能区分に存在する樹木を間伐可能樹木として選定する間伐可能樹木選定処理（ステップＳ１３）とを有する点である。なお、所定の傾斜角度未満というのは、ここでは３５度未満の傾斜角度とする。この角度は、林業用ハーベスタが走行可能な角度（間伐可能傾斜角度）である。

傾斜区分画像データ作成処理（ステップＳ１１）は、間伐対象区域に対応する森林域画像データにおけるメッシュ化（例えば、５０ｃｍ×５０ｃｍでメッシュ化）された５０ｃｍの画素ごとのデジタル標高モデルデータ（地表の標高を表すデータ）を作成して、当該メッシュ化されたデジタル標高画像データに基づいて傾斜区分画像データを作成する。具体的には、メッシュ化されたデジタル標高画像データに基づいて、隣接する画素を直角三角形の底辺（５０ｃｍ）として、当該当該直角三角形の高さを隣接する画素の標高差として求める処理を各隣接する画素において行って、隣接する画素ごとの傾斜角度を求めることにより、傾斜区分画像データを作成することができる。なお、メッシュ化されたデジタル標高画像データは、上空から調査対象となる森林域にレーザー光を照射することよって得られたレーザー計測データから作成された３次元の点群データに基づいて作成することができる。

図１４は、ディスプレイ上に表示された傾斜区分画像を示す図である。図１４に示す傾斜区分画像は、図１３のフローチャートにおける傾斜区分画像データ作成処理（ステップＳ１１）によって作成された傾斜区分画像である。なお、図１４は各メッシュにおける傾斜角度が各樹冠に対応付けられて示されている。図１４においては、各樹冠が薄い灰色から濃い灰色のグラデーションで示されているが、灰色が濃くなるほど傾斜が大きいことを示している。

なお、図１４はモノクロ画像であるため、色の濃さによる傾斜の大きさの区別が判別しにくいが、図１４の元となるカラー画像上では、傾斜の大きさは黄色（傾斜角度が小）から赤（傾斜角度が大）のグラデーションで表わされており、赤色に近いほど傾斜が大きい。なお、ここで、最も薄い色の黄色は傾斜角度が０度（平坦）であることを示しており、赤色は傾斜角度が８０度以上であることを示している。なお、図１４に示す傾斜区分画像における間伐対象区域内の傾斜角度は、全体的に間伐可能傾斜角度に収まっている。

ところで、傾斜区分画像データ作成処理（ステップＳ１１）において傾斜区分画像データが作成されると、当該傾斜区分画像データに基づいて林業用ハーベスタが走行可能な傾斜角度（この場合、３５度未満）となっている間伐可能区分画像データを作成する（ステップＳ１２）。そして、間伐可能区分画像データと森林資源情報とに基づいて、間伐可能区分に存在する樹木を間伐可能樹木として選定する（ステップＳ１３）。なお、ステップＳ１３において、間伐可能区分に存在する樹木を間伐可能樹木として選定する処理は、間伐可能区分画像データと森林資源情報とから得られる個々の樹木の樹冠位置とを重ね合わせることにより、間伐可能区分に存在する樹木を間伐可能樹木として選定することができる。このようにして、間伐可能樹木（間伐対象区域内において林業用ハーベスタが走行可能な区域に存在する樹木）を選定した後の処理は、実施形態１又は実施形態２において説明した処理を行う。

すなわち、実施形態１又は実施形態２において説明した間伐率設定処理（ステップＳ２０）及び間伐条件設定処理（ステップＳ３０）の後に行われる樹木間距離算定処理（ステップＳ４０）は、間伐対象区域内において間伐可能傾斜角度となっている区域（間伐可能区域とする。）に存在する樹木（間伐可能樹木）について行われる。なお、間伐率設定処理（ステップＳ３０）及び間伐条件設定処理（ステップＳ４０）は、実施形態１又は実施形態２において説明した処理と同様であり、また、樹木間距離算定処理（ステップＳ４０）及び当該樹木間距離算定処理以降の処理は実施形態１又は実施形態２において説明した処理と同様であるため、ここではその説明は省略する。

このように、実施形態３に係る間伐木選定方法においては、間伐対象区域内において林業用ハーベスタが走行可能な間伐可能区域に存在する樹木（間伐可能樹木）について、実施形態１又は実施形態２において説明した間伐木選定を行うため、３５度未満の傾斜に存在する樹木だけが間伐対象となる。このため、実施形態３に係る間伐木選定方法において、図７に示すような間伐木資源情報を作成すると、実施形態３に係る間伐木選定方法において作成される間伐木資源情報は、３５度未満の傾斜に存在する樹木の中から間伐木として選定された個々の樹木についての資源情報となる。これにより、実施形態３に係る間伐木選定方法において作成される間伐木資源情報に基づいて、間伐作業を行えば、間伐木として選定された樹木は、林業作業の効率化のために導入されている林業用ハーベスタが走行可能な位置に存在することとなり、当該樹木の伐採、枝払い、玉切りなどの作業の機械化が可能となる。

以上説明したように、各実施形態に係る間伐木の選定方法によれば、現地での人手による間伐木の選定に頼ることなく、間伐木の選定を効率的にかつ森林の現況に適応するように行うことができる。また、間伐条件として、「上層間伐」、「下層間伐」又は「全層間伐」のいずれかの間伐条件を設定して間伐木の選定を行うことにより、どのような間伐を行うかといった間伐目的に応じた間伐木の選定が可能となる。

また、いずれかの間伐条件（例えば、「全層間伐」）によって間伐木の選定を行った場合、図６に示すような間伐対象区域画像、図７に示すような間伐木資源情報、図８に示すような間伐木概要情報及び図９に示すようなヒストグラムなどをディスプレイ上に表示させることによって、間伐対象区域における間伐木の選定結果を机上にて容易に把握することができる。また、実際に現地に行くときには、これら各情報をプリントアウトしたもの、または、これらを現地で表示可能なモバイル端末を持参し、かつ、ＧＰＳを併用することで、熟練者でなくても容易に間伐木として選定された樹木をマーキングすることができる。なお、これは層間伐による間伐木の選定に限られることなく、「上層間伐」又は「下層間伐」においても同様である。

以上、実施形態１～３に係る間伐木の選定定方法について説明したが、続いて、実施形態１～３に係る間伐木の選定定方法に対応する間伐木の選定装置について説明する。

図１５は、実施形態１に係る間伐木の選定方法に対応する間伐木の選定装置１及び実施形態２に係る間伐木の選定方法に対応する間伐木の選定装置２を説明するために示す図である。間伐木の選定装置１及び間伐木の選定装置２は、図１に示すフローチャート及び図１０に示すフローチャートの各ステップの処理を行うための装置であり、図１におけるステップＳ１０において取得すべき森林資源情報を記憶している森林資源情報記憶部１０と、図１におけるステップＳ２０における間伐率の設定を行う間伐率設定部２０と、図１におけるステップＳ３０の間伐条件の設定を行う間伐条件定部３０と、図１におけるステップＳ４０の樹木間距離算定処理を行う樹木間距離算定部４０と、図１におけるステップＳ４０の樹木間距離算定処理を行う樹木間距離算定部４０と、図１におけるステップＳ５０の最短距離の２本の樹木を選択する最短距離樹木選択部５０、図１におけるステップＳ６０の間伐木の選定を行う間伐木選定部６０と、図１におけるステップＳ７０の間伐率の判定を行う間伐率判定部７０と、図１におけるステップＳ９０の処理（間伐木に関する間伐木資源情報を作成するとともに、当該間伐木資源情報に基づいて間伐木の概要に関する間伐木概要情報を作成する処理）を行う間伐木集計部９０とを有している。

なお、間伐率設定部２０は、実施形態１に係る間伐木の選定方法における樹木の本数に基づく間伐率の設定が可能であるとともに、実施形態２に係る間伐木の選定方法における樹木の材積に基づく間伐率の設定が可能である。また、間伐率判定部７０は、間伐率設定部２０において樹木の本数に基づく間伐率Ｒの設定がなされた場合には、Ｒｎ＝（現時点の選定間伐木数／間伐対象区域に存在する樹木の本数）×１００の演算結果が間伐率Ｒ（設定間伐率Ｒ）未満であるか否かの判定（Ｒｎ＜Ｒか否かの判定）を行う。一方、間伐率設定部２０において樹木の材積に基づく間伐率Ｒの設定がなされた場合には、Ｒｎ＝（現時点の選定間伐木合計材積／間伐対象区域に存在する樹木合計材積）×１００の演算結果が間伐率Ｒ（設定間伐率Ｒ）未満であるか否かの判定（Ｒｎ＜Ｒか否かの判定）を行う。

このように構成されている間伐木の選定装置１及び間伐木の選定装置２を構成する各構成要素が、図１に示すフローチャート又は図１０に示すフローチャートに示す手順に沿って処理を行うことにより、実施形態１に示す間伐木の選定方法及び実施形態２に係る間伐木の設定方法において説明した間伐木の選定が可能となる。

図１６は、実施形態３に係る間伐木の選定方法に対応する間伐木の選定装置３を説明するために示す図である。間伐木の選定装置３が図１５に示した間伐木の選定装置１及び間伐木の選定装置２と異なるのは、間伐対象区域を複数に区分した各区分ごとの傾斜角度を表す傾斜区分画像データを作成する傾斜区分画像データ作成部１１と、当該傾斜区分画像データに基づいて所定の傾斜角度未満の傾斜を有する区分を間伐可能区分とした間伐可能区分画像データを作成する間伐可能区分画像データ作成部１２と、間伐可能区分画像データと森林資源情報とに基づいて、間伐可能区分に存在する樹木を間伐可能樹木として選定する間伐可能樹木選定部１３とを有する点である。
傾斜区分画像データ作成部１１は、図１３に示すフローチャートにおける傾斜区分画像データ作成処理（ステップＳ１１）を行うものであり、間伐可能区分画像データ作成部１２は、図１３に示すフローチャートにおける傾斜区分画像データ作成処理（ステップＳ１２）を行うものであり、間伐可能樹木選定部１３は、図１３に示すフローチャートにおける傾斜区分画像データ作成処理（ステップＳ１２）を行うものである。
そして、間伐可能樹木選定部１３により選定された間伐可能樹木に対して、間伐率設定部２０により設定されている間伐率及び間伐条件設定部３０により設定されている間伐条件に基づいて、樹木間距離算定部４０が行う樹木間距離を算定する処理、最短距離樹木選択部５０が行う最短距離樹木を選択する処理及び間伐木選定部６０が行う間伐木を選定する処理を行う。

このように構成されている間伐木の選定装置３を構成する各構成要素が、図１３に示すフローチャートに示す手順に沿って処理を行うことにより、実施形態３に示す間伐木の選定方法において説明した間伐木の選定が可能となる。

また、間伐木の選定装置１、間伐木の選定装置２及び間伐木の選定装置３は、当該各間伐木の選定装置１～３に含まれる上記各構成要素（図１５及び図１６参照。）が有する機能がコンピューターのプログラムとしてインストールされており、上記各構成要素に所定のデータを与えることによって、それぞれの構成要素が持つ機能がコンピューターのソフトウエア上で実行されるものである。

なお、本発明は上述の実施形態に限られるものではなく、本発明の要旨を逸脱しない範囲で種々変形実施可能となるものである。たとえば、下記に示すような変形実施も可能である。

（１）上記各実施形態においては、森林を形成する樹種としては、カラマツを例示したが、カラマツに限られるものではなく、アカマツ、スギ、ヒノキなど種々の樹木を対象とすることができる。なお、本発明は、主に、針葉樹の植栽林の間伐を想定しているが、針葉樹の植栽林の間伐に限られるものではない。

（２）上記各実施形態においては、樹木の本数に基づく間伐率を設定する場合には、間伐率Ｒ（設定間伐率Ｒ）を３０％とし、樹木の材積に基づく間伐率を設定する場合には間伐率Ｒ（設定間伐率Ｒ）を３５％としたが、これは一例であって、間伐率Ｒは任意に設定可能である。また、実施形態３において、間伐可能傾斜角度は、林業用ハーバスタが走行可能な傾斜角度として、３５度を設定したが、これに限られるものではなく、より大きな傾斜角度においても走行可能な林業用ハーベスタが存在する場合には、より大きな傾斜角度を設定可能である。また、林業用ハーバスタが走行可能であるか否かとは関係なく、任意の角度を間伐可能傾斜角度として設定可能である。例えば、間伐可能傾斜角度として、３５度よりも小さい２０度を設定したり、さらに小さい１５度を設定したりすることも可能である。このように、間伐可能傾斜角度として任意の角度を設定可能とすることによって、林業用ハーバスタが入り込めないような区域で人力による伐採作業を行わざるを得ない場合、人力による伐採作業の危険性を考慮した間伐木の選定を行うことが可能となる。

（３）上記各実施形態においては、間伐木選定処理（ステップＳ６０が「上層間伐」、「下層間伐」又は「全層間伐」に基づいて間伐木を選定する際には、胸高直径（ＤＢＨ）で間伐木を選定する場合を例示したが、胸高直径（ＤＢＨ）ではなく、樹高によって間伐木を選定するようにしてもよい。例えば、間伐条件として「上層間伐」が設定されているとすれば、最短距離の２の本の樹木のうち、樹高の高い方の樹木を選定し、間伐条件として「下層間伐」が設定されているとすれば、最短距離の２の本の樹木のうち樹高の低い方の樹木を選定するようにしてもよく、また、間伐条件として「全層間伐」が設定されているとすれば、最短距離の２の本の樹木のうち、樹高の高い樹木と樹木とを交互に選定するようにしてもよい。

１，２，３・・・間伐木選定装置、１０・・・森林資源情報記憶部、１１・・・傾斜区分画像データ作成部、１２・・・間伐可能区分画像データ作成部、１３・・・間伐可能樹木選定部、２０・・・間伐率設定部、３０・・・間伐条件設定部、４０・・・樹木間距離算定部、５０・・・最短距離樹木選択部、６０・・・間伐木選定部、７０・・・間伐率判定部、９０・・・間伐木集計部、ｎ・・・処理回数、Ｒ・・・間伐率（設定間伐率）、Ｒｎ・・・現時点の処理回数の間伐率、Ｓ１０・・・森林資源情報取得処理、Ｓ１１・・・傾斜区分画像データ作成処理（傾斜区分画像データ作成ステップ）、Ｓ１２・・・間伐可能区分画像データ作成処理（間伐可能区分画像データ作成ステップ）、Ｓ１３・・・間伐可能樹木選定処理（間伐可能樹木選定ステップ）、Ｓ２０・・・間伐率設定処理（間伐率設定ステップ）、Ｓ３０・・・間伐条件設定処理（間伐条件設定ステップ）、４０・・・樹木間距離算定処理（樹木間算定ステップ）、Ｓ５０・・・最短距離樹木選択処理（最短距離選択ステップ）、Ｓ６０・・・間伐木選定処理（間伐木選定ステップ）、Ｓ７０・・・間伐率判定処理（間伐率判定ステップ）、Ｓ９０・・・間伐木集計処理（間伐木集計ステップ）

Claims (11)

1. 間伐対象区域に存在する樹木に関する情報が個々の樹木ごとに記録されている森林資源情報に基づいて、前記間伐対象区域に存在する樹木のうち、間伐すべき樹木を間伐木として選定する間伐木の選定方法であって、
   前記間伐対象区域における間伐率を設定する間伐率設定ステップと、
   前記間伐対象区域においてどのような間伐を行うかを示す間伐条件を設定する間伐条件設定ステップと、
   前記間伐対象区域に存在する樹木の個々の樹木間の距離を算定する樹木間距離算定ステップと、
   前記樹木間距離算定ステップで算定された距離のうち、最も近接している２本の樹木を最短距離樹木として選択する最短距離樹木選択ステップと、
   前記間伐条件設定ステップによって設定されている前記間伐条件に基づいて、前記最短距離樹木選択ステップによって選択されている２本の樹木のうちの１本の樹木を前記間伐木として選定する間伐木選定ステップと、
   を有し、
   前記樹木間距離算定ステップ、前記最短距離樹木選択ステップ及び前記間伐木選定ステップは、前記間伐率設定ステップにおいて設定されている前記間伐率に達するまで順次行うことを特徴とする間伐木の選定方法。
2. 請求項１に記載の間伐木の選定方法において、
   前記間伐率は、前記間伐対象区域に存在する樹木の本数に基づく間伐率であることを特徴とする間伐木の選定方法。
3. 請求項１に記載の間伐木の選定方法において、
   前記間伐率は、前記間伐対象区域に存在する樹木の合計材積に基づく間伐率であることを特徴とする間伐木の選定方法。
4. 請求項１～３のいずれかに記載の間伐木の選定方法において、
   前記間伐条件は、上層間伐、下層間伐又は全層間伐のいずれかであり、
   前記間伐条件として前記上層間伐が設定されている場合には、前記間伐木選定ステップは、前記最短距離樹木選択ステップによって選択されている２本の樹木のうちの大きい方の樹木を前記間伐木として選択し、
   前記間伐条件として前記下層間伐が設定されている場合には、前記間伐木選定ステップは、前記最短距離樹木選択ステップによって選択されている２本の樹木のうちの小さい方の樹木を前記間伐木として選択し、
   前記間伐条件として前記全層間伐が設定されている場合には、前記間伐木選定ステップは、前記最短距離樹木選択ステップによって選択されている２本の樹木のうちの１本の樹木を所定の規則に基づいて前記間伐木として選択することを特徴とする間伐木の選定方法。
5. 請求項４に記載の間伐木の選定方法において、
   前記所定の規則は、現時点の処理回数を表す数値が奇数か偶数かによって、前記最短距離樹木選択ステップによって選択されている２本の樹木のうちの大きい方の樹木又は小さい方の樹木のいずれかを前記間伐木として選択することを特徴とする間伐木の選定方法。
6. 請求項４又は５に記載の間伐木の選定方法において、
   前記樹木の大きさは、当該樹木の胸高直径で表わされることを特徴とする間伐木の選定方法。
7. 請求項１～６のいずれかに記載の間伐木の選定方法において、
   前記間伐対象区域を複数に区分した各区分ごとの傾斜角度を表す傾斜区分画像データを作成する傾斜区分画像データ作成ステップと、
   前記傾斜区分画像データに基づいて所定の傾斜角度未満の傾斜を有する区分を間伐可能区分とした間伐可能区分画像データを作成する間伐可能区分画像データ作成ステップと、
   前記間伐可能区分画像データと前記森林資源情報とに基づいて、前記間伐可能区分に存在する樹木を間伐可能樹木として選定する間伐可能樹木選定ステップと、
   をさらに有し、
   前記間伐可能樹木選定ステップにより選定された前記間伐可能樹木に対して、前記間伐率設定ステップにより設定されている前記間伐率と、前記間伐条件設定ステップにより設定されている前記間伐条件とに基づいて、前記樹木間距離算定ステップが行う樹木間距離を算定する処理、前記最短距離樹木選択ステップが行う最短距離樹木を選択する処理及び前記間伐木選定ステップが行う間伐木を選定する処理を行うことを特徴とする間伐木の選定方法。
8. 請求項１～７のいずれかに記載の間伐木の選定方法において、
   前記間伐木選定ステップによって選定した間伐木に関する間伐木資源情報を作成するとともに、当該間伐木資源情報に基づいて間伐木の概要に関する間伐木概要情報を作成する間伐木集計ステップをさらに有し、
   前記間伐木資源情報には、間伐木として選定された個々の樹木の番号、個々の樹木の位置情報、個々の樹高、個々の樹木の胸高直径、個々の樹木の材積が含まれており、
   前記間伐木概要情報には、間伐木選定前の樹木に関する情報及び間伐木選定後の選定間伐木に関する情報として、本数、所定面積当たりの本数、平均胸高直径、平均樹高、合計材積、所定面積当たりの合計材積及び選定間伐木の間伐率が含まれていることを特徴とする間伐木の選定方法。
9. 間伐対象区域に存在する樹木に関する情報が個々の樹木ごとに記録されている森林資源情報に基づいて、前記間伐対象区域に存在する樹木のうち、間伐すべき樹木を間伐木として選定する間伐木の選定装置であって、
   前記間伐対象区域における間伐率を設定する間伐率設定部と、
   前記間伐対象区域においてどのような間伐を行うかを示す間伐条件を設定する間伐条件設定部と、
   前記間伐対象区域に存在する樹木の個々の樹木間の距離を算定する樹木間距離算定部と、
   前記樹木間距離算定部で算定された距離のうち、最も近接している２本の樹木を最短距離樹木として選択する最短距離樹木選択部と、
   前記間伐条件設定部によって設定されている前記間伐条件に基づいて、前記最短距離樹木選択部によって選択されている２本の樹木のうちの１本の樹木を前記間伐木として選定する間伐木選定部と、
   を有し、
   前記樹木間距離算定部が行う樹木間距離を算定する処理、前記最短距離樹木選択部が行う最短距離樹木を選択する処理及び前記間伐木選定部が行う間伐木を選定する処理を、前記間伐率設定部によって設定されている前記間伐率に達するまで順次行うことを特徴とする間伐木の選定装置。
10. 請求項９に記載の間伐木の選定装置において、
    前記間伐対象区域を複数に区分した各区分ごとの傾斜角度を表す傾斜区分画像データを作成する傾斜区分画像データ作成部と、
    前記傾斜区分画像データに基づいて所定の傾斜角度未満の傾斜を有する区分を間伐可能区分とした間伐可能区分画像データを作成する間伐可能区分画像データ作成部と、
    前記間伐可能区分画像データと前記森林資源情報とに基づいて、前記間伐可能区分に存在する樹木を間伐可能樹木として選定する間伐可能樹木選定部と、
    をさらに有し、
    前記間伐可能樹木選定部により選定された前記間伐可能樹木に対して、前記間伐率設定部により設定されている前記間伐率及び前記間伐条件設定部により設定されている前記間伐条件に基づいて、前記樹木間距離算定部が行う樹木間距離を算定する処理、前記最短距離樹木選択部が行う最短距離樹木を選択する処理及び前記間伐木選定部が行う間伐木を選定する処理を行うことを特徴とする間伐木の選定装置。
11. 請求項９又は１０に記載の間伐木の選定装置において、
    前記間伐木選定部によって選定した間伐木に関する間伐木資源情報を作成するとともに、作成された間伐木資源情報に基づいて間伐木の概要に関する間伐木概要情報を作成する間伐木集計部をさらに有し、
    前記間伐木資源情報には、間伐木として選定された個々の樹木の番号、個々の樹木の位置情報、個々の樹高、個々の樹木の胸高直径、個々の樹木の材積が含まれており、
    前記間伐木概要情報には、間伐木選定前における樹木に関する情報及び間伐木選定後のおける選定間伐木に関する情報として、それぞれの本数、それぞれの所定面積当たりの本数、それぞれの平均胸高直径、それぞれの平均樹高、それぞれの合計材積、それぞれの所定面積当たりの合計材積及び選定間伐木の間伐率が含まれていることを特徴とする間伐木の選定装置。

Images