第16号

国立防災科学技術セソター研究報告

1976年12月

508．7，578．7：681．3

降ひょ う記録計記録の電子計算機による白動検測矢崎忍＊

国立防災科学技術セソター

Automatic皿easurement
Hai1−Padbyof

Records011

ComputerBy

S．YazakiNα肋伽1肋5θ肌乃Cθ〃θ7仰刎5α∫伽P〃閉

4o〃，τo妙o

Abst工act
NRCDP
of

se七s

abou亡300hai1・pads

hai1fa11distribu士ion．The

hai1fa1l
bai1fa1lasaredepressionsis（1）depressionh泓i1−pa（1consis士son士heasautomate士his

foi1．From七heof

Kan士oahard

depression

djstric七士oa1uminum

paせem±heobtaindata

foi1and−records

in士ensi七yoftheis士aken

on

（2）

Nex士士his（3）士hen（4）pla七e

thebyisonpaper士oob士ain

white

circu1arspot

digitized

wi亡h

paper

pa七ternsomecompu士er

isandsuchasprin士edasisascannedandproperties

of

each

direc士ion

arec01nputed，aretransforn1ed−in士oand士hen士hese

hai1s士one

es亡ima七ion

pエesen七ed．An

ou七1ineof士heapa。士tern．

fo11ows．

At丘rs亡士hisarea，eccen七工一ici士yeach

nortbern

es士ima七ed一．

A士ria1士oprocesses

in士heva1ues

sp〇七（whjch

diame士er，i刀cidenceang1eand

isina

c1ear

binaryb1ack

ground．

igurescanner，regarded

thekine七ic

asanEache11ipse）suchcor工esponding

Properties

asofenergy．

1．はじめに
国立防災科学技術セソター第1研究部では，ひょう雲の構造と行動についての調査研究のため，降ひょう発生頻度の高い北関東平野部に約300台の簡易型降ひょう記録計を設置L，降ひょう分布調査を行っている（清野ほか1961）．

この簡易型降ひょう記録計とは，発泡スチロールの台の上に硬いアルミ箔（厚さ0．1mm

叉は0．2mm；大きさ（受ひょう部）220皿m×220m皿）を密着してのせたものであり，降＊第3研究部計測研究室

一47一

降ひょう記録計記録の電子計算機による自動検測一矢崎

ひょうはこのアルミ箔の凹みとして記録される（写真1）．

局所的変化の激しい降ひょう分布のようすを適確につかむためには，平面的に密なデータ得ることが重要である．そのためにはを多数の記録計を設置しなげればならず，したがって個々の記録計については，簡易性と低廉性は最も要求されるものである．この設録計自体ほ，記録の質や精度の点で十分なものとはいえないが，上記の要求によく合致したものであり，今のところこれにとって代るもの
写真1は考えにくい．

Photo．1

アルミ箔上のひょう痕
Hai1records

on

ana1uminum

foil．

さて，この記録から我々が求めたいのは農作物に対する被害率をよく説明しうる量である．小元ほか（1976）にょると，被害率を決定る最も重要たす物理量は，降ひょうの衝突エネルギーであることが知らせている．しかし，衝突エネルギーをある程度正確に求めるには，多くのアルミ箔について個々のひょう痕の大きさや数を逐一計測しなければならず，大へんな労力と時間を要する．
そこで現在は，簡便な方法として，一人の研究者がアルミ箔全体を一見して判断し，「目測強度」といった値を定め代用している．小元ほか（1976）によると，このr目測強度」と衝突エネルギーの問にはかなり良い対応がつくことがわかった．アルミ箔上の痕跡記録そのものの精度を考えると，この程度の評価でも一応の用はなすと考えられるが，この方法は読み取りがいくぶん主観的であることは避けられない．

したがって，今後さらに大量の資料を解析するにあたっては，これを客観的かつ迅速に処理する方法を開発しなければならない．本論文はこの目的のための電算機を用いた処理の自動化について述べられる．

2・

図形読取装置によるアルミ面上のひょう痕の読み取リ

平面上に記録された情報を計算機で処理するにはまずこれを走査Lて格子点上の値の列を得なげれぱならない．

このように平面記録を読み坂る機械として我々が持っているのは，図形読坂装置（諾星1976）である・この機械は平面形を光学的に走査するものであり，本来の使用法としては，平面土の情報が濃度色の濃淡の形で与えられたものを読み坂るのに用いるものである．

ところでアルミ面上のひょう痕の場合，その情報の形は面の凹み（変形）であるから，走査して得たいのは面上の各点におげる面の変位の値である．しかし今述べたようにわれわれ

一48一

国立防災科学技術セソター研究報告第16号

1976年12月

受光器

W
光源

SN
みW
W

資料面

図1光Fig．1

学

Detec士ing

系
part

図2面の傾きと光学系の読み
of

the丘gurescanner．

Fig．2

Relationandbetweenreadingof

gradient

theofasurfacescanner．の図形読取装置では，これを直接得ることはできたい．

そこでいかに記録を読み取るかが問題になるが，これには二つの考え方がある．一つは，アルミ板をともかくこの機械で直接読んで，あとは計算機プログラムで処理しようというものである．もう一つは，これを機械にかける前に何らかの方法で処理して，その変位を紙面の濃度に変換して，それからこの機械で読み坂るという考え方である．
初めに，前者の方法を検討してみる．

この場合の読み（出力）は，表面の光学的状態（濃度や散乱度）による成分と，面の変形による成分の和である．したがってまず，これから光学的な変動によるものを坂り除かなげればならないが，これは資料によって異なり，また同一資料上でも，表面のサビやよごれによって不均一が生じている場合があるといった具合で，なかなか面倒である．

次に，表面の状態が一様であった（あるいは一様にした）として，変形と読みの関係を調べてみる．この機械の検出部（光学系）は図1のよう変 ⊥止になっている．これによると読みは主に面の傾き（変読み位の勾配）による（正確には変位自体にもよるが，これは傾きによる変動に較べて小さい）．傾きて読みの関係はたとえば図2のようになるだろう．このように一読みの形が非対称だから，面の傾きは2次元量として扱わなけれぱならない（対称ならぼ，傾きの絶対値だけなら1次元量である）．したがって，傾きから読みへの変換が求まったとしても，逆の変換は1つの変換式では求まらない．原理的には，偏徴分方程式の解とし

一49一

図3凹みの形とそれに対する読みFig．3

Shapeofareadingof

tha士．

depre

ssionand降ひょう記録計記録の電子計算機による自動検測一矢崎て求まるはずだが実際にはほとんど不可能である．すたわち，このときの読みから機械的な計算だげによって変位を求めることはできないということである．

そこでこの場合に可能な方法として考えられるのは，対象に制限を加えることによって（ひょう痕だけを考えるということ）いわゆるバターソ認識を行うことである．凹みの形とその読みの関係はたとえぼ図3のようになるだろうから，これを念頭において逆に読みの形から凹みの形を推算することになる・

しかし実際には，必ずしもこのようなきれいなパターソカミ得られるとは限らず，またひょう痕以外の凹凸がある場合もあり，そのプログラムは非常にむずかしいものになることが予想される1しかも，そのような複雑なプログラムを作ったとしても，それから得られる情報は，せいぜい変位に関する二値情報（凹んでいるかいないかの区別）であろうと思われる．

このように，この直接読むという方法は機械にかけるまでの手問が全くか

写真2写真のアルミ箔を紙に写し取ったもの

Photo．2からないという点では最も望ましいもPrinted

pa士ernof

asurfaceof

the・1uminumf・i1・fPhot・・1・のであるが，その計算機処理が面倒であり，したがって計算機の使用時問を多く費すことになり，今回の目的には通さないと思われた．

今回行ったのは，アルミ箔の表面にイソクを塗って，これを紙に写し取るという方法である
実際に用いたのは謄写用イソクであるが，これは油性で乾きにくいので薄く均一に塗れ，またにじみにくいので小さなものもよく採れて好都合である（写真2）．

この方法は，イソクの塗り方や紙のあて方にかなり依存するように思われるが，実際に行ってみると十分安定した（試行ごとのバラツキが少ない）ものが得られることがわかった・

これは，板面の変位をひょうが当ったか当らないかによって，二値化したことにほかたらない．変位を二値化するということは，ひょう痕による凹みを平面的な形（切り口の形ということにする）と深さに分けるとする）と深さのほうの情報を捨ててしまうことである．したがって，これが妥当であるかどうかを検討しておかなければならない・消極的な理由としては，図形読取装置を用いる方法では，どんな方法でも深さを求めることは一般に困難であるということがある．しかし物理的にも後に述べるように，切り口の形と深さは近似的には独立に変化しえない．すなわちこれらの間には関数関係があるということである．したがって，この記録自体の精度を考えると切り口の形を調べるだけで一応十分であるということが一50一

国立防災科学技術セソター研究報告第16号

1976年12月

し・える．

次に，この方法の利点をいくつかあげてみる．

まず，これは面の変位が紙面の濃度に変換されているから，これを図形読取装置で読んだデータは，そのまま変位のデータとして使えるということがある．しかもこの場合，白黒の差が十分大きいから，これを数値的に二値化するのが大へん楽である．

それから，読み坂り図形の情報の形が濃度であるから，修正がしやすいということがある・実際の降ひょう計アルミ箔には，ひょう痕と区別つけにくい凹みなどがある場合がある
が・これを計算機で判断するのは容易ではない．したがってこういうものは人問が見てひょう痕以外は塗りつぶしてしまえば，計算機プログラムは非常に簡単なものになる．要はこれらの処置がオフライソでできるということである．

この方法は・機械にかけるまでの手問がいくらかかるが，あとの機械による処理は非常に速く・結局全体の所要時問を考えると，いちばん速いものの一つであると思われる．
3．

計算機プログラム

粒子（ひょう痕）の識別と，個々の粒子について形と大きさの算出を行なう．

（1）粒子の識別対象とする図形は二値図形であるから，粒子を記述するには例えばその境界の座標を粒子ごとに全部書き上げられれぼよい・粒子とは一つのつながった（連結な）領域であるから，粒子の識別は領域のつたがりを見ていけばよい．一方，計算機内部（主記憶装置）へのデータの読み込みは，走査線単位でしかも走査と同じ順序で行っていくのが最も望ましい．したがって領域のつながりは隣接する走査線上でのつながりを次々に追っていくという方法で行なわれる．

1本の走査線が粒子を切る部分を，ここでは「弦」ということにする．各粒子は弦の集りと見なせるから，粒子の識別は，弦を粒子ごとに分類することによってなされる．ところで各粒子は，その性質上凸図形と考えられる．このことに注意すると，次のようなプログラムが可能である．

走査線を順に第1列，第2列，…とよぶことにする．隣接する2本の走査線第ゴー1列と第ゴ列に対して，第ゴー1列上のある弦にはすべて異った番号がついているとする．そこで，これら2本の走査線上の弦を比較することによって，第ゴ列の弦に次の方法で番号をつける
図4）．

（i）第ゴ列の弦が第ゴー1列の弦aに接していたら，bにはaと同じ番号をつける．
（ii）第ゴ列の弦にc接する弦が第ゴー1列になかったら，cには新しい番号をつける．

この操作のくり返しによって，すべての弦の分類，すなわち粒子の識別ができる．しかしこのままでは，このデータは新しいデータの読み込みによって失なわれてしまうから，これ一51一降ひょう記録計記録の電子計算機による自動検測一矢崎

cfP

9hqr

↓Ph

図4q

粒子と「弦」

Fig．4

Spots

and

apair0土adjacentSCa．nning1ineS．

をたとえば弦の両端の座標に変換して，番号ごとにあらかじめ（主記憶設置の中に）用意した表に書き込んでおく．しかし一般には，全部の粒子のデータをしまうほどの容量を主記憶装置はもっていないから，そこで（iii）第ゴー1列の弦dに接する弦が第ゴ列になかったら，この弦の属する粒子は走査し終ったと判断し，そのデータ（粒子の境界の座標）をこの時点で処理してしまうか，または磁気ディスク等に書いておく．

こうすれぼ，1個の表をくり返し使うことができ，用意しておく表の数はr1列上にある粒子の数の最大値」以上であれぼよいことになる．

以上の方法は，初めに述べたように，粒子が凸であるときにだけ有効なものである．ところで実際の読み坂りデータ上の各粒子は，全体的な形としては凸であるとしても，細かく見れぱ境界の形に小さい凹凸（ぎざぎざ）があるかも知れない。今回，後に述べるように図形
読坂装置で読んだデータを，そのままこのプログラムのデータとして使用したいから，このような場合も考慮しておかなげれぱならたい．この小さい凹凸は意味のないものであるから，これによって結果に大きな誤差を生じる可能性さえ取り除いておけばよい．それには，上の（i），（ii），（iii）の手順の前に次の修正を行っておけばよい・すなわち（o）第ク列の弦hに接する弦が第ゴー1列に2つ以上あったら，これらのうち1番左のも（e）のを除いて他（f，9）は消してLまう．また第i−1列の弦pに接する弦が第i列に2個以上あったら，これら（9，r）をつないでしまう（図5）．

この処置によって，図6のような修正が行なわれたことになる．これによってこのプログラムは小さなノイズに対して安定になる．

なお，この方法は各粒子が凸からのはずれが小さくて，それを無視できるような場合にのみ有効なものであるが，一般の（凸とは限らない）粒子を対象としたときにも，これを一部修正することによって同様のプログラムが可能である・（2）形状の計算

まず何を求めるかを決めなげれぼならない，最も重要なものは面積である．これだげを求

一52一

国立防災科学技術セソター研究報告

第1−1例

一5＼、第1例

b第16号

1976年12月

d／

C

図6修

図5r弦」の

修

正

Fig．5

Correction

of

da亡a．

Fig6

An

正例

examp1eof

correctionof

asp〇七一めるのならぼ．プログラムはもっと簡単になる．もう少し詳しく記述するにはたとえば，これをだ円と見て面積長軸の向き離心率を求めることができる．

ここでは，これらを与えられたデータから求める計算法を結果だけあげておく，走査方向をκ軸，送り方向をツ軸にとって連続座標で表わすことにする．与えられたデータは弦の両端の座標，すなわち（κエ（ツ），κ。（ツ），ツ）（κ。＜。）という形をしているとしたが，ここではむLろ
2
幻＝一（κ。十κ。）
1∠κ＝κrκ。
で定義される（κ。け），∠π（ツ），ツ）という形で求めておいたほうが便利である（図7）．
まず面積（∫）は，y
X＝X。（y）

長軸？

・】（y）・。（y）・。（y）

図7広Fg，7標系

Deinitionofvariables．

X・一／1：〃（ツ）・ツで求められることはすぐにわかる．

9

次に長軸の向き（φ）を求めるには，まず（κ。

（ツ，ツ）を結ぶ線の向きθを求める（図8）・実
際にこれを求めるには，（κ。（ツ），ツ）のデータから最小自乗法で図8長軸と直線κ＝耽（ツ）Fig．8

κ。（ツ）＝tanθ・ツ十定数

Majoraxis

andα1inex。（y）ofan・11ip・e．

一53一x＝

降ひょう記録計記録の電子計算機による自動検測一矢扉として求めてもよいしもっと簡単にはt。。θ＝κ・（ツ・）一κ・（ツ・）ツ2一ツ1としてもよい．するとgはこのθと前に求めた面積Kを用いて・一；…一・（。．1黒n、㌧）（ただし一号≦θ・ρ≦云で・・はθと同符号）で求められる．ここに，

A＝1gS2π2（ツ。一ツ。）4とおし・た．

離心率（θ）は，上で求めたθとgを用いて
L2tanθ

θ （t・・2p−t・・θ）〉lr干t…2φ十（1−t・・22ρ）t・・θ
一で求まる．

（3）システムと処理時間について
このプログラムは，その処理速度が非常に大きいのが特徴である．

まずデータであるが，これは図形読坂装置で走査した生のデータをそのまま用いることができる．

それは，このデータ白体がほとんど二値データになっており，改めて二値化するという手続きは実質的に要らないこと，求めるものが面積とか概形であり，小さなノイズによる影響は平均化されてしまって結果にはほとんどきかないこと，さらにプログラム自体にもノイズを無害化する処置がしてあることによる．

次に処理時問であるが・実際に調べてみると走査線当りだいたい100msであった1走査線のデータ数2，200の場合）．これは図形読坂装置のドラムの回転周期（200ms）と同じオーダーの値である．データの処理が入出力より速くできれば，いわゆる実時問処理ができる．

図形走査は必ずしもドラムの各回転ごとに行う必要はなく，たとえば2回転に1回走査すれば入カデータの転送時問問隔は400msとなる．したがってこれを適当に選べば，このプログラムは図形走査に対して実時問で処理を行えるはずである．

この方法によれば，データを磁気テープなどに往復させる時問の節約になり，全体の処理時問は約半分になる．

今日行ったものでは，画面の大きさ220mmx220mlnで，メヅシュの大きさ0．1mmxC・1mmならぼ約16分，0・2mmx0・2mmならば約4分，O．4叫m×0．4mmのときに約2分で読み取り，演算，出力が可能である．

図11は・写真2のパターソを0・2mm×0．2mmのメッシュで読んで処理した結果を，モー54一

国立防災科学技術セソター研究報告第16号

1976年12月

ニターとしてライソプリソタに打たせたものである・ただしここでは面積だけを求めた．数字は各粒子の面積を円の直径に換算して，mm単位で測った値である．下の表は大きさ毎の度数を数えたもの．

4．

アルミ面上のひょう痕からひょう粒子への変換

ひょう痕の測定ができたら，次にこれからひょう粒子自体を求めなけれぱならない．それには，衝突するひょう粒子のもつ独立な物理量を調べる必要がある．なぜならぼ，衝突によってできる凹みの形状には，これらすべての量が関与するはずであるから．

一般にこれらの量として考えられるのは，形・大きさ（体積）・硬さ・質量・速度・入射角および衝突時の粒子の向きである．しかし，今の目的のためにはこれらすべてを考える必要はない．このうち変動の大きないくつかの量（成分）だけを考えれば十分である．

まず，硬さは一定としてよいだろう．形による影響はあまり大きくないと思われるから形も一定（球形）とする・次に，密度（質量1体積）は一定であり，速度はほぽ終端速度（ヒれも質量と大きさだけで決まる）であるから，大きさ，質量，速度のうち独立なものは1個だけになる・この量として，ここでは運動エネルギー（1／2×質量×速度2）を用いることにすると，結局独立なものとしては，運動エネルギーと入射角だけになる．

ところで，入射角もほぼ0すなわち面に直角に当たるとすれば，残りは1個だけであるから，これによって作られる凹みについても，この近似では独立な変量は1個ということになる・これが，ひょう痕を測定するとき，切り口の大きさと深さを別々に測る必要はないといった理由である．

しかし，ひょう痕については3個の量：面積（∫），離心率（θ），向き（g）を求めたから
ひょう粒子についても3個の量：運動エネルギー（E），入射角（θ，g）をもって記述することにする（図9）．

zy③

ExyS，e予x
図・9一ひょう粒子とひょう痕の記述
Fig．9

Description

depressjonof

on亡heahai1s亡one（righ七）andfoi1（orpatern）．

一55一aspot

onathe，降ひょう記録計記録の電子計算機による自動検測一矢崎

Lたがって求める変換は（∫，θ，ρ）から（亙，θ，ψ）へのそれである．

まずψと〆は他と独立で，g＝〆としてよいだろう．そこでこの変換は一般に，
E＝E（∫，θ）
θ二θ（∫，2）〆＝Pと書ける．しかし実際にはこれを，
E＝E（∫）
θ＝θ（2）〆＝9と近似して十分であると思われる．すなわち，ひょう痕の切り口の面積がひょう粒子の建物エネルギーに，形（離心率）が入射角に対応すると考えるわけである．

このうち最も重要なのは最初の関係である．そこで，これを実験によって求めた．実験結果は図10のとおりである．

この実験で衝突粒子として用いたのは鉄製のボールである．鉄とひょう粒では密度が違うから，上にあげたすべての量を再現することはできない．合わせたのは大きさと運動エネルギーである（したがって質量，速度，運動量などは同じではない）・それは，この2個が凹みの形状に最も大きく効くと考えたからである．また硬さも同じではない．これら（同じでな
ひょう粒gよう痕㍗子の直径

直径（mm一（mm〕

25107

205×l06

152×l06106

105xlOヨ12345678

ひよう粒子の運動エネルギー（erg〕

図10 ひょう痕の大きさとそれに対応するひょう粒子の大きさおよび運動エネルギーの関係
Fig．10Transfo工mation
andkineticofenergy

size
of

the一56一ofaspot

hails士one．

intosize

国立防災科学技術セソター研究報告第16号

1976年12月い分）が結果にどの程度影響を及ぽすかは検討する必要があるかも知れない．

5．おわリに今回行った読み取り方法は，その計算機の使用時間が非常に短くてすむという点で望まし
■

2．

■

12

一1

ヨ22

一二︒

13221R＾Nκ

l2
ヨ01＾MFTF嗜

【〕

O−

12〕．…

1−2：：：：：＾5．

COuNTS｛02甲19j■ ＾ ．． ．
＾一5 ．． ．

6｛

一9罧i11高．．．・・｛o．．．■

1

図11写真2のパターソの処理結果
Fig．11 Aresu1tofprocessingofapatemofPho士o．2．Eachigureindecates

areain士erms

ofdiame士er（uni士mm）of

一57一

eachspot．

2■

降ひょう記録計記録の電子計算機による自動検測一矢崎

いものである．また測定精度も，この簡易型記録計から期待される精度に対して一応満足できるものである．

ところで，この方法ではひょう痕の深さに関する情報は捨ててしまった．これには一応合理的な理由があったが，中には凹みの形を次元的に見なければ正しい判断ができないようなものもある．その一つは，2個以上のひょう痕が近接して重なったものであるもう．一つは，ひょう粒子がアノレミ板の上を跳ねて2回以上面にあたる場合があるが，この識別である．これは2回目以後の衝突は速度がちがうから，切り口の大きさと深さの関係を見れば判断がつくはずである．しかしこのような判断を計算機で行なうには，非常に面倒なプログラムを要し，所要時問も長くなってしまう．今回の場合のようにあまり精度が要求されないものに対しては，むしろ少し雑でも計算機処理の速いものを作って，それで得た結果に何らかの一たとえぼ粒子数が少ないときには人問が見て判断し，粒子数が多いときには統計的に処理をするといった一補正を加えることで十分であると思われる．

なお今回は，ひょう粒子の運動ニネルギーを求める関係だげを求めたが，入射角についても必要があれば同様に求めることができる．

最後に，今目行ったアルミ箔の表面を紙に写し取るという方法は，菅原正已当セソター前所長の御教示によるものである．また本報告を書くにあたって，第一研究部小元敬男異常気侯防災研究室長をはじめ同研究室の諸氏にいろいろ御助言をいただいた．記して謝意を表します．

参考文献

1）諸星敏一（1976）：電子計算機による強震記録の読み取り（第2報）．国立防災科学技術セソター研
究報告，第16号，29−45．

2）小元敬男，清野害谷，八木鶴平，米谷恒春（1976）：農作物のひょう害と降ひょう強度との関係．国立防災科学技術セソター昭和51年度研究発表会講演要旨，2−5．
3）清野離，小元敬男（1976）：1972年8月3日のひょう雲の構造と行動，国立防災科学技術セソター
研究報告，第14号，53−63．

（1976年8月11目

一58一

原稿受理）