DEIM Forum 2016 F3-1

Storyline を適用した実数値型時系列データ可視化の一手法— 気象データを例として
八木佐也香†

伊藤

貴之†

高塚

正浩††

† お茶の水女子大学大学院人間文化創成科学研究科 〒 112–8610 東京都文京区大塚 2-1-1†† School of Information Technologies, J12, University of Sydney, NSW 2006 Australia
E-mail: †{sayaka,itot}@ itolab.is.ocha.ac.jp, ††masa.takatsuka@ sydney.edu.au
あらまし

身の回りの時系列データの分析において，可視化は有効である．本報告では実数値型の時系列データに対し，データ要素の短時間単位の特徴と全体のクラスタ変遷の双方を同時に観察することを目標に，要素間の交流や共起性の表現に用いられる storyline ベースの可視化手法を提案する．本手法では局所時間帯ごとに時系列数値情報の類
似度を算出，隣接時間帯間にてクラスタ群のマッチングを適用することで，類似した要素同士が画面内で近接するような storyline を描く．また，実数値変化の重要度にもとづいて各要素の不透明度を設定することで，重要度の高い部分を強調し，ユーザの知覚効果の高い可視化結果を得る．本報告ではアメダス気温データを適用した例を紹介し，storyline を俯瞰することで気温変化の傾向を把握する．さらに storyline 上で特徴的な気温変化を示す観測所の詳細な数値変化を，折れ線表示との連携可視化を行うことで直感的に理解可能であることを示す．
キーワード

情報可視化，時系列データ，storyline

1. は じ め にで，その相互関係に関する視認性は大きく低下する．そこでク
ラスタリングなどを適用して要素を構造化した上で可視化する

ビッグデータ時代の到来に伴い，多種多様なデータベースにことが有効である．ここで時系列データにおいては要素間の相時系列性を有する情報が蓄積されるようになった．これらの時互関係も時間とともに変遷するため，クラスタ構造の時間的変系列データの分析・予測はデータマイニングや知覚情報処理等

遷を表現することも重要となる．一方でこれらの可視化手法にの分野においては既に数多くの研究が進められている．一方

おいて，データの部分的な特徴とクラスタ構造の時間的変遷をで，時系列データの観察と理解，またそれを扱う業務での能動

同時に表現するのは簡単ではない．

的な意思決定を支援する道具として，情報可視化が有用である．

クラスタ構造の時間的遷移の表現に向いた可視化手法として，情報可視化には，目で見て全体像を把握し，対話操作によってSankey Diagram と storyline がよく利用されている．Sankey

データ中の必要な部分を能動的に選択できる点に加え，特に時

Diagram はクラスタ間の流量変化を観察するのに適している．

系列データにおいてはリアルタイムな観察のための瞬時な理解

一方，storyline は個々の要素について，クラスタから別のクにおいても有効である．一般的に情報可視化手法は，一次元，ラスタへの移動を観察するのに向いている．Storyline はもと

二次元，三次元，多次元，木構造，グラフ，時系列の 7 種類の

もと XKCD というウェブコミックに手書きイラスト “Movie

データ構造に大別される [1]．その中でも我々はこれまで，時Narrative Charts”として紹介されたものであり [4]，物語の登系列データの可視化に関する新しい手法の開発に取り組んでき場人物の共起行動を時系列のダイアグラムで表現している．図た [2]．

1 にあるように，x 軸が時刻，各曲線が映画の各登場人物を表

多量の実数値で構成される時系列データの可視化では，デーしており，複数の登場人物が同一のシーンに登場した場合に線タを構成する実数値の大小そのものだけでなく，実数値どうしが近接するような可視化結果となる．これらの可視化手法は任の関連性や変化量などを観察することも重要である．例えば多意の 2 要素間の流量や共起に関する時間変化を表現するのに用数の商品の売上を一定時刻ごとに集計して作成される時系列

いられるのが一般的であり，各要素が独立に実数値を有するよ

データでは，商品単体の売上の推移だけでなく，どの商品と同

うな時系列データに対してこれらの可視化手法を適用した事例

時期に売れるか，季節や時間帯との相関はあるか，といった関はまだ少なく，議論の余地が残っている．

連性を観察することによって，より深い分析が可能になり，売上予測や予兆発見などに貢献できると考えられる．
時系列データの可視化では折れ線グラフとヒートマップの 2
種類の視覚表現が特に多用されている [3]．しかしいずれの視覚表現においても，大量のデータ要素（以下「要素」と称する）の相互関係を全て表現するのは難しい．例えば折れ線グラフでは，要素に対応する折れ線を 1 個の画面空間に大量に描くこと

図 1 手書きによる storyline の例（ [4] より転載）．

以上を踏まえ，我々は実数値で構成される時系列データをstoryline として表現する可視化の一手法を提案する．本手法でにもとづく可視化手法と，ヒートマップにもとづく可視化手法
が特に多用されてきた [3]．

はまず，時系列データを構成する各要素に対して局所時間帯ご折れ線グラフにおいては，要素数の多いデータにおけるクとにクラスタリングを適用し，隣りあう局所時間帯間のクラス

ラッタリングがその可読性を妨げてきた．その解決策として折タ群にマッチングを適用する．その結果にもとづいて，クラス

れ線のサンプリングが有効であり，有効なサンプリングを実現

タに対応する長方形群を画面配置し，長方形群を通過する線分するために時系列データのクラスタリングがしばしば適用されとして各要素を描画する．描画時には次を工夫することで，よてきた．我々は折れ線グラフにもとづく可視化手法において，り深い分析を可能にする．

局所時間帯ごとに要素をクラスタリングしてサンプリングする•

各要素に割り当てられた数値やメタ情報から線分の色相

ことにより，可読性が高く，かつデータの特徴を見逃さないよを算出することで，要素のクラスタ移動の要因を理解しやすく

うな可視化を実現した [2]．しかしこのようなアプローチでは，する．

どの要素が局所時間帯ごとにいつクラスタを変遷したか，と•各要素の数値変化の重要度から線分の不透明度を算出することで，重要な数値変化を有する線分に視線を届きやすく

いった情報を視認することができない．
またヒートマップにもとづく可視化手法において，時間帯全する．

体にわたる類似度で要素を並べ替え，さらに局所時間帯ごとにまた，本手法では各要素のクラスタ移動の要因を分析するため見られる特徴的な数値パターンをマークすることで，データ大に，storyline 中の局所における実数値変化を別の画面領域に折域にわたるクラスタ構成と，局所に見られる特徴的な現象を同れ線グラフ表示するという連携可視化を実装している．

時に表現する手法を提案した [5]．しかしこのようなアプローチ

本手法を用いてクラスタの時間的変遷を観察することで，次の 3 種類の要素群を視覚的に分離して理解するのが容易になる．

では，局所時間帯ごとにクラスタリングを適用できない．

2. 2 Storyline を適用した可視化手法•

長時間にわたって同一クラスタに属する要素群•

外れ値を有するために小さなクラスタに属する要素群

化する，という問題設定は実数値型以外の時系列データにおい•ある時には特定のクラスタに属し，またある時には別のて活発に議論されている．人物や単語を要素とみなしたとき，クラスタに属する，というようなふらふらとした時間的遷移を

その 2 要素間の距離（例えば共起性）の時間変化を可視化する，有する要素群という問題がその最たる例である．この課題においてフレキシ

例えば，夏には商品 A と B の売上の時間的推移が類似して，冬ブルに各要素間の共通性を描く手法として，XKCD の手描きには商品 A と C の売上の時間的推移が類似している，というイラスト [4] に着想を得て storyline のメタファを用いた時系列

時系列データを構成する各要素の共通性に関する変遷を可視

商品データがあるとする．図 2（右）のように本手法でこれをデータ可視化手法が近年数多く提案されている．Ogawa らは可視化することで，商品 A が夏から冬にかけて別のクラスタにソフトウェア開発のリポジトリを題材として storyline のレイ

移動していることが storyline 上で明確に表現できる．これにアウトを自動生成する手法を提案した [6]．この手法による可視

よって，商品 A に対する店舗での陳列には注意が必要である，

化結果は，小〜中規模のプロジェクトによるソフトウェア開発というような知識の発見を期待できる．

の構造変化の観察に有用である．一方でこの可視化結果には，線分の小刻みな蛇行や線分どうしの重なりが数多く含まれており，時間軸に沿って各線分の動きを追うことは困難である．

Tanahashi らは Ogawa らの手法を拡張し，遺伝的アルゴリズ
ムを用いて線分の配置を最適化する手法を提案した [7]．この手法では，線分のずれや重なり，画面上の空白を制御することで，Ogawa らの手法の問題点を解決し，線分間の共起性に着目し
図2

商品の売上データの可視化例の比較．

やすい結果を実現している．しかしこれらの手法にはスケーラ
ビリティと計算時間に課題がある．数百以上の要素数を有する

本報告では，実数値型時系列データの一例として，全国 376
箇所で観測されたアメダス気温データを提案手法に適用した事

データにおいて線分間の重なりが非常に多くなり，各要素を目で追うことが困難になる．

例を紹介する．可視化結果において，storyline 上の色や不透明

この問題を解決するため，StoryFlow は各クラスタに詳細

度の変化から，特定の時間帯で特徴的な気温変化をしている要度制御を適用し，Storyline 全体の特徴を捉えやすくした [8]．

素に着目可能であること，また，折れ線グラフ表示によって各

StoryFlow では各要素の階層構造を考慮した配置に加え，イン要素の実数値数値変化を直接観察することで，特徴的な変化をタラクティブな線分の並べ替え，バンドリングなどの機能も実もたらした要因を推察可能であることを示す．

現している．また，Reda らは storyline のメタファを用いつ

2. 時系列データ可視化に関する関連研究

つも，数百規模の要素の描画が可能な手法を提案している [9]．

2. 1 実数値型時系列データの可視化手法

固定し，各スロット内に線分を配置している．このアプローチ

実数値型時系列データの可視化には旧来から，折れ線グラフ

はユーザが特定のクラスタ間の変化に着目したい場合，有効でReda らの手法では，ｙ軸上にコミュニティを示すスロットを

あるといえる．しかしこれらの手法は，上述のように既知の関係性やキーワード・トピック等を入力情報として要素の共通性を表現しており，実数値で構成される時系列データを対象にしたものではない．また，要素間の距離などにもとづく強調表示などは行われていない．
実数値型の時系列データにおいても，局所時間帯ごとのクラスタリング結果に焦点をおいた可視化手法は最近になっていく
つか報告されている．例として Lex ら [10] は，クラスタごと，
局所時間帯ごとに断片化されたヒートマップを Sankey Dia-gram 風に接続する可視化手法を提案している．しかし Sankey
Diagram と同様な表現を採用したこの手法では，時系列デー
タを構成する個々の要素の変遷に着目することが難しい．またTurkay ら [11] は，局所時間帯ごとにクラスタを縦に並べ，各図 3 画面上での変化量が大きく見える storyline の例（文献 [7] より転載）．

クラスタを通過するように要素を曲線で描く可視化手法を提案している．この手法は storyline の簡易な生成方法の一種と考えることもできる．本報告の提案手法における描画形式は Turkay
らの描画形式に類似しているが，次節にて議論する点において提案手法は Turkay らの手法を改良しているといえる．

2. 3 Storyline の適用に関する注意点

3. 提 案 手 法
本章では storyline を適用した実数値型時系列データ可視化手法を提案し，その各処理について論じる．本手法の処理の流れは下記の通りである．

本節では，Storyline を実数値型時系列データに適用する際（ 1 ） 入力データが対象とする時間帯を「局所時間帯」と呼の注意点と，それらを改良するためのアプローチを議論する．

ばれる短い時間帯に分割し，各局所時間帯内の時系列実数値にクラスタの位置関係: storyline による描画では，要素を表現す

クラスタリングを適用する．

る曲線の縦方向の移動を小さくすることが重要である．本手法（ 2 ） 各局所時間帯を構成するクラスタ群を長方形とみなし，では隣りあう局所時間帯において多くの要素を共有するクラス

画面空間に配置する．

タが画面上で隣接するように配置することで，曲線の移動を抑（ 3 ） 入力データを構成する各要素を，クラスタに対応する

える．上述の Turkay らの手法 [11] においても，要素の共有性

長方形群を通過する storyline として描画する．その際に，各を考慮した並び替えが行われている．一方で storyline による要素のクラスタ移動が強調されるように storyline の透明度を描画を見ただけでは，実数値型時系列データの数値自体の変遷

自動調節する．

を理解することは難しい．そこで，要素の共有性に加えて数値

図 4 に 6 個の要素集合 (A〜F) をもつ時系列データを用いた，変化も加味してクラスタを配置することが有用である．

処理の流れの概要を示す．

重要な要素の強調表示: storyline では多くの場合において，画

なお本章では，時系列データは m 個の要素の集合 D =面上で大きく移動する曲線が目立って見えるため，この曲線に{a1 , a2 , ..., am } であり，要素 ai は n 個の各時刻における実数

対応する要素が重要な要素（あるいは大きな変化を有する要値の集合 ai = {ai1 , ai2 , ..., ain } であるとする．また時系列デー素）であると連想しがちである．よって実数値の時系列データ

タが対象とする時間帯は 2 時刻を区間とする Nw 個の局所時間帯

を storyline で可視化すると，実数値の変化の小さい要素に対に分割されているとする．また，入力データの対象となる時間帯応する曲線が画面上で大きく移動しているために目立ってしまを Nw 個に区切った局所時間帯のセットを W = {w1 , ..., wNw }う，というような直感に合わない可視化結果を生むことがある．

とし，i 番目の局所時間帯 wi を次の通り記述する：

例えば図 3 において，左下部の赤の曲線は短時間で急激な変化を有しているように見える．しかし多くの既存手法において，

wi = {tstart , tend , Ci }，ただし Ci = {ci1 , ci2 ..., cik }．

曲線の画面上での移動量は，個々のクラスタに含まれる曲線のここで，tstart ，tend は wi の開始時間および終了時間であり，本数やクラスタの配置などによって決定されるため，必ずしもk は局所時間帯 wi におけるクラスタ数，cij は j 番目のクラス

実数値の変化量に比例するとは限らない．この問題を解決する

タを示す．

手段として，曲線の形状以外の視覚変数によって重要な要素を

3. 1 各局所時間帯におけるクラスタリング

強調表示することが有効である．

まず時系列データを構成する各局所時間帯に対して，要素の数値表示へのインタラクション: storyline による描画は数値自クラスタリングを適用する．ある局所時間帯が j 番目から k 番

体の変遷を表現しない．そのため，既存手法ではヒートマッ

目 (j < k) までの nw 個の時刻で構成されている場合に，本処理プを併用する数値表現 [10] が用いられているが，厳密に数値

では i 番目の要素 ai から実数値 [aij , ..., aik ] を抽出し，これをを読み取ることは難しい．この問題を解決する一手段として，

nw 次元ベクトルとして扱う．このベクトル群に対してクラス

storyline 上で興味深い局所を対話的に指定し，別の画面領域でタリングを適用することで，各局所時間帯における要素を有限数値表示することが考えられる．

個のクラスタに分類する．本手法では非階層型クラスタリング

図 4 本手法の処理の流れ．

として k-means を適用し，Davies-Bouldin Index（DBI）[12]に含まれる要素数，acij は cij を構成する実数値の平均，α はを各局所時間帯におけるクラスタ数の妥当性指標に利用する．

0<=α<= 1 を満たす定数である．
()∑(cij ∩ c(i+1)k )1α+ (1 − α)ncij + nc(i+1)kacij − ac(i+1)k

DBI は各クラスタ内の凝集性が高く，クラスタ間距離が離れているほど良いという前提にもとづき，式 1 の評価関数を最小化(4)するクラスタ数を採用してクラスタ分割する．

1∑
Si + Sj
Ri ，ただし Ri = max
j =i
|
n i=1
dij

この式で第 1 項は両クラスタ間の要素の共有性を表すものであ

n

DB =

(1)

り，第 2 項は両クラスタ間の数値の不変性を表すものである．
両者を考慮して最大マッチング問題を解くことで，要素を表現

ここで，Sn は i 番目のクラスタ Ci の中心点 ci と Ci 内の各点する曲線の移動量を抑え，かつクラスタ配置結果から数値変化との距離の平均，dij はクラスタ Ci ，Cj の重心間の距離であを説明しやすい可視化結果を得る．

り，次のように表される．

dij = ∥ci − cj ∥p =( n∑

3. 2. 2 クラスタのマッチングと配置決定)1

本手法では各局所時間帯を構成する各クラスタを長方形で表

p|ck,i − ck,j |

p

(2)

k=1

1 ∑
∥ej − ci ∥p
Sn =Ti j=1

現し，3. 2. 1 節の処理によってマッチングされたクラスタが画面上で左右に隣接するように，以下の処理によって各長方形の

Ti

(3)

Ti はクラスタ Ci の大きさ，ej はクラスタ Ci 内の要素を示す．

位置を決定する．
（ 1 ） 画面左端の局所時間帯を構成するクラスタを上から順に配置する．

上述のように，現時点での我々の実装では非階層型クラスタ（ 2 ） 左から i 番目の局所時間帯について配置が終わったら，リングを適用しているが，階層型クラスタリングを含め他のク

左から (i + 1) 番目の局所時間帯について，マッチングされた

ラスタリング手法を適用することも可能である．

クラスタができるだけ画面上で左右に並ぶように各クラスタの3. 2 クラスタの画面配置
Storyline を適用した既存の可視化手法の処理手順は，次の 2
種類に大別される．

•

XKCD のような表現を実現するために，曲線としての各要素の配置を最適化する方法 [6], [7]•

まず要素にクラスタリングを適用し，続いてクラスタの位置を決定する．

3. 3 Storyline の描画
続いて時系列データを構成する各要素を線分として描画する．
我々の実装では次の手順に沿って線分を描く．

•

i 番目の局所時間帯にて当該要素が属するクラス

タ cij の 左 右 端 お よ び 上 下 端 の 座 標 値 を [xminij , xmaxij ],位置を決定する方法 [9], [10], [11][yminij , ymaxij ] とする．このとき当該要素を表現するため本手法では数百・数千の要素を有する大規模な時系列データをに 2 点 (xminij , yij ), (xmaxij , yij ) を結ぶ線分を描く．ただし対象とするため，前者の適用は非現実である．そこで本手法でyminij <= ymaxij とする．
= yij <

は後者を適用するが，本手法では時間的に隣接する局所時間帯•

i 番目の局所時間帯と (i + 1) 番目の局所時間帯を連結す間でクラスタのマッチングを取る必要がある．そこで本手法でるために，2 点 (xmaxij , yij ), (xmin(i+1)j , y(i+1)j ) を結ぶ線

は，以下の手順でクラスタの共通度を求め，それに沿って隣接分を描く．

する局所時間帯間のクラスタをマッチングする．

このとき我々の実装では，storyline の色を以下のように定義

3. 2. 1 クラスタの共通度の算出

する．本手法では各要素に 1 次元の属性が付与されていると仮本手法では，i 番目および (i + 1) 番目の局所時間帯を構成す

定し，その属性値に対応した色相を線分に与える．1 次元の属るクラスタ集合 Ci ，Ci+1 に対して，以下の式 (4) が最大とな性がない場合には例えば，各要素の実数値 aij で線分の各頂点

るような最大マッチング問題を適用する．ここで，cij ∩ c(i+1)kに色相を与えることも可能である．また各要素の数値変化の重はマッチングされる 2 クラスタが共有する要素数，ncij は cij

要度に応じて，線分を構成する各頂点に不透明度を与える．重要度が大きい部位において不透明度を高く設定することで，観を付与することで，既に用いられている他の視覚変数（線分の察者の注意を引くことができる．現時点の我々の実装では実数

形状や色相）に干渉することなく，ユーザの注意を重要度の高

値 aij の微分値および各クラスタ内での要素の安定性にもとづい部位に向けられると考える．また，不透明度から定量的な差

いて重要度を算出している．詳細は次節にて論じる．

異を読み取ることは困難であるものの，ユーザは不透明度を見また各要素のクラスタ移動の要因を分析するために，story-line 中の局所における実数値の変化を観察することが有用である．我々の実装では，storyline を描画した画面領域でマウス操て直感的にその部位における重要度を判断可能である．

4. 実 行 結 果

作で長方形を描くと，その長方形を通過する要素群の実数値の4. 1 実行環境と処理時間時間変化を別の画面領域に表示するインタラクション機能を有

我々は JDK (Java Development Kit) 1.7.0 を用いて提案手法

する．具体的には，図 5 のように画面の上半分に storyline をを実装し，Windows 8.1 (64bit) 搭載のラップトップ PC (CPU描き，上述の長方形描画操作によって要素群を選択すると，折1.80 GHz, RAM 4.0GB) 上で実行した．表 1 は上記の環境にれ線グラフによってその要素群を画面の下半分に描く，というおける DBI を適用したクラスタリング，クラスタの画面配置

連携可視化を実装している．なお，複数回の操作を行った場合，およびインタラクションの平均処理時間を示す．なお，いずれ複数回の選択の論理積を取ることで，観察する要素を絞り込むの事例も元データ内の時間帯数は 248，局所時間帯サイズを 24ことが可能である．

時間に設定した．
表1

各処理における処理時間 (秒，5 回の平均).

要素数 クラスタリング 画面配置

インタラクション

86

0.20

0.099

0.013

376

0.64

0.083

0.16

上記の結果から，既存手法 [6]，[7] に比べ，本手法は対話的に操作する上で十分高速であるといえる．

4. 2 適用データ
本章では，全国 376 箇所の観測所について 3 時間ごとに観測図 5 本手法におけるインタラクション機能．

3. 4 Storyline の不透明度算出
前節で論じた通り，現時点での我々の実装では，次の 2 種類の基準を各時刻における重要度とみなし，storyline の不透明度算出に用いている．

されたアメダス気温データを提案手法に適用した事例を紹介する．可視化結果において，x 軸は時刻，1 本の曲線は 1 観測所を示し，近接する曲線は同じクラスタに属することを表す．各
クラスタ内の曲線の並びは観測所 ID 順，また観測所 ID から観測所を 8 つの地方に分類し，ColorBrewer 2.0 [16] の配色方
針に沿って各地方に色を割り当てた（図 6 参照）．

隣接する局所時間帯間 においては，要素の実数値 aij の微分値を基準にする．具体的には，隣接する局所時間帯間でクラスタを移動する曲線について，対応する要素の実数値の微分

dab /|dt| を求める．ここで，dab はクラスタ AB 間の距離であり，|dt| は微分値を計算する時間幅を示す．

1 つの局所時間帯内 では，要素の安定性を基準にする．具体的には，i 番目の局所時間帯について，(i − 1) 番目の対応するクラスタと共通する要素は不透明度を高く，別のクラスタから移図 6 曲線の色と地方の対応．

動してきた要素は不透明度を低く設定する．
不透明度にもとづく視覚表現を採用した理由は以下の通りで

ここで，storyline の読み取りの凡例として，図 7 に 12 観測

ある．人間の目は知覚的に変化の激しい部分に対して反応す

所の気温データを適用した実行結果の一例を示す．この例ではる [13] という Marr の理論を適用して，我々は storyline 上で局所時間帯サイズを 12 時間とした．

の要素のクラスタ移動という重要な現象が発生している部位に3 番目の局所時間帯において，各クラスタには上から順に，東

視覚変化の大きな表現を適用したいと考えた．可視化における北地方，東北から九州地方，北海道地方に属する観測所が含ま

視覚表現の手段として，Bertin [14] によって提案された視覚変れている．クラスタ 2 に着目すると，東北地方を示す青緑の線

数集合が知られている．それに対して MacEachren は，画面内分と中部地方を示す黄緑の線分がそれぞれ 1 本ずつ不透明度がの同一座標上に複数の情報を載せる場合に不透明度を用いるの

高く表示されている．これは 2 番目の局所時間ではクラスタ 1

が有効であると述べている [15]．本手法においても，storylineに属する当該線分が，3 番目の局所時間帯でクラスタ 2 に切りを構成する線分において重要度の高い部位にのみ高い不透明度

替わっているためといえる．

中国地方を中心に，北海道，東北，九州地方の一部の観測所で，特に北陸地方に属する観測所が多く見られた．これにより，当該観測所群が上述のような気温変化を示す際に，クラスタを移動しやすいといえる．

図 7 12 観測所の気温データを適用した実行例．

また本章では，3. 3 節で述べた要素選択機能を用いた折れ線
表示についても示す．折れ線の表示においては，x 軸が時刻，y軸が気温，各折れ線が各観測所を示す．

4. 3 2009 年 3 月の事例
図 8 は 2009 年 3 月の 1 カ月全体の気温変化を表す storyline
である．この事例では局所時間帯サイズを 24 時間ごと，DBI

図 9 2009 年 3 月下旬にクラスタを移動した要素群の気温変化．

の距離関数にはユークリッド距離を用いた．クラスタ数は 2 となる局所時間帯が最も多く，上旬および下旬にクラスタ数が 3から 6 個に分かれた局所時間帯が複数存在することが読み取れる．

図 8 2009 年 3 月の気温変化の storyline による表示例．
図 10

クラスタ数が変化した要因を分析するため，中〜下旬の 2 週

図 9 で選択した要素群と局所時間帯 wd , we との論理積を取った結果．

間分を拡大した結果を図 9（上）に示す．図 9（上）において，桃色の円で囲んだ部分は，クラスタ数が 2 から 3 へ変化して
おり，しかもクラスタ移動部分の不透明度が高く表示されていることから，重要度が高いと考えられる．当該要素の数値変化を観察するため，図 9（上) の白い破線で示した 2 クラスタを選択し，論理積を取ったものを図 9（下）に示す．青い破線で囲まれた時間帯 wa および青でハイライトされた時間帯 wb がstoryline 上での各局所時間帯に相当する．折れ線グラフの形状に着目すると，局所時間帯 wa では朝晩の気温が低く，昼に気温が上昇するという変化をしているのに対し，局所時間帯 wb
では昼になっても気温がほとんど上がらず，気温が急降下するという変化をしていることが読み取れる．また，図 10（上）において，桃色の円で囲んだクラスタ移動部分も同様に不透明度が高く表示されており，重要度が高いと推察できる．図 9 で選択した局所時間帯に加え，さらに 2 回の選択を行い，論理積を取った要素群の気温変化を図 10（下）に示す．当該要素群の折れ線の形状に着目すると，局所時間帯 wd では昼に気温が上昇しているのに対し，局所時間帯 we では昼に気温があまり上がらず，夜に向けて気温が低下するという，上記と同様の変化をしていることが観察できる．図 11 に上述の 2 箇所のクラスタ
移動部分を拡大した結果を示す．いずれの時間帯においてもクラスタを移動した要素群は，storylinen の色から中部，近畿・

図 11

図 9，10 のクラスタ切替部分を拡大した結果．

また，図 12（上）の storyline において，1 番上のクラスタは
青の北海道の観測所を多く含み，安定して推移しているが，桃色の円で囲んだ部分のみ不透明度が高くなっており，重要度が高いと考えられる．当該期間の気温変化を観察するため，wc ，wd の 2 箇所を選択した結果を図 12（下）に示す．折れ線の形状から，2 箇所の局所時間帯での気温変化のしかたは共通しているものの，wd′ では wc に比べ高い気温で推移していることがわかる．この平均気温の差により，当該時間の重要度が高いと判定されたと考えられる．以上のように，storyline 上でのクラスタ数の推移や不透明度の変化を観察することにより，興味深
い実数値変化を発見しやすくなると考えられる．

[7]

[8]

[9]
図 12 安定したクラスタ内において，storyline の不透明度が高い部
分を選択した結果．

[10]

5. まとめと今後の課題
本報告では，storyline を用いた実数値型時系列データの可[11]視化手法を提案した．本手法では局所時間ごとに時系列実数値情報の類似度を算出し，隣接する局所時間帯間のクラスタ群にマッチングを適用することでクラスタの配置を決定し，類似し[12]た要素同士が画面内で近接するような storyline を描く．また，本手法では各要素の数値変化の重要度にもとづいて要素の不透明度を設定することで，重要度の高い部分を強調し，ユーザの[13]知覚に近い可視化結果を得る．さらに本報告では，アメダス気温データを storyline に適用した例を紹介し，大規模なデータ[14]の中から特徴的な気温変化を示している観測所群を発見し，折れ線表示により細かな数値変化を観察するというアプローチの[15]有効性を示した．
今後の課題として，気象以外のデータでの適用事例を示し，[16]

ユーザテスト等によって可視化結果の妥当性を検証することを計画中である．また手法部分の拡張として，各要素の実数値間[17]の類似度を算出する単位となる局所時間帯を可変にすることで，よりデータの局所的特徴に応じたクラスタ切り替えの数や時間を検出したいと考えている．この点に関して，SAX 法 [17] に対し，不等間隔離散化を適用したパターン抽出手法 [18] などを用いることが考えられる．

文

献[1] B. Shneiderman. The eyes have it: a task by data typetaxonomy for information visualizations. Proceedings 1996
IEEE Symposium on Visual Languages, pp. 336–343, 1996.
[2] Sayaka Yagi, Yumiko Uchida, and Takayuki Itoh. A
Polyline-based Visualization Technique for Tagged Timevarying Data. 2012 16th International Conference on Information Visualisation, pp. 106–111, 2012.
[3] 伊藤貴之. 視覚協創学 (5): マッピング技術〜特に時系列データ可視化技術の体系化に向けて〜. 第 17 回計算工学講演会論文集,
Vol. 17, , 2012.
[4] Randall Munroe. Movie narrative charts. http://xkcd.
com/657/, accessed June, 2015.
[5] Maiko Imoto and Takayuki Itoh. A 3d visualization technique for large scale time-varying data. In Information Visualisation (IV), 2010 14th International Conference, pp.
17–22. IEEE, 2010.
[6] Michael Ogawa and Kwan-Liu Kl Ma. Software evolution[18]

storylines. SOFTVIS ’10 Proceedings of the 5th international symposium on Software visualization, pp. 35–42,
2010.
Yuzuru Tanahashi and Kwan Liu Ma. Design considerations
for optimizing storyline visualizations. IEEE Transactions
on Visualization and Computer Graphics, Vol. 18, No. 12,
pp. 2679–2688, 2012.
Shixia Liu, Yingcai Wu, Enxun Wei, Mengchen Liu, andYang Liu. StoryFlow: Tracking the evolution of stories.
IEEE Transactions on Visualization and Computer Graphics, Vol. 19, No. 12, pp. 2436–2445, 2013.
Khairi Reda, Chayant Tantipathananandh, Andrew Johnson, Jason Leigh, and Tanya Berger-Wolf. Visualizing theevolution of community structures in dynamic social networks. Computer Graphics Forum, Vol. 30, No. 3, pp. 1061–
1070, 2011.
A. Lex, M. Streit, C. Partl, K. Kashofer, and D. Schmalstieg. Comparative Analysis of Multidimensional, Quantitative Data. IEEE Transactions on Visualization and Computer Graphics, Vol. 16, No. 6, pp. 1027–1035, 2010.
C. Turkay, P. Filzmoser, and H. Hauser. Brushing Dimensions - A Dual Visual Analysis Model for High-Dimensional
Data. IEEE Transactions on Visualization and Computer
Graphics, Vol. 17, No. 12, pp. 2591–2599, 2011.
David L. Davies and Donald W. Bouldin. A Cluster Separation Measure. IEEE Transactions on Pattern Analysis
and Machine Intelligence, Vol. PAMI-1, No. 2, pp. 224–227,
1979.
David Marr. VISION - A Computational Investigation into
the Human Representation and Processing of Visual Information. Henry Holt and Co., Inc. New York, 1982.
Jacques Bertin. Semiology of Graphics: Diagrams, Networks, Maps. University of Wisconsin Press, Madison, Wisconsin, 1983.
Alan M MacEachren. How Maps Work: Representation,
Visualization, and Design. The Guilford Press, New York,
1995.
Mark Harrower and Cynthia A Brewer. Colorbrewer 2.0.
http://colorbrewer2.org/, accessed February, 2015.
Jessica Lin, Eamonn Keogh, Stefano Lonardi, and Bill Chiu.
A Symbolic Representation of Time Series, with Implications for Streaming Algorithms. Proceedings of the 8th ACM
SIGMOD Workshop on Research Issues in Data Mining
and Knowledge Discovery, pp. 2–11, 2003.
Miho Ohsaki, Hidenao Abe, and Takahira Yamaguchi. Numerical Time-Series Pattern Extraction Based on Irregular
Piecewise Aggregate Approximation and Gradient Specification. New Generation Computing, Vol. 25, No. 3, pp.
213–222, 2007.