芸術科学会論文誌 Vol. 14, No. 5, pp. 188 – 197

Storyline を適用した実数値型時系列データ可視化の一手法
八木佐也香 1) (学生会員)

伊藤貴之 1) (正会員)

1) お茶の水女子大学大学院人間文化創成科学研究科

高塚正浩 2)

2) シドニー大学情報理工学部

A Storyline-based Visualization Technique for Consecutive
Numerical Time-varying Data
Sayaka Yagi1) Takayuki Itoh1) Masahiro Takatsuka2)
1) The Graduate School of Humanities and Sciences, Ochanomizu University
2) Faculty of Engineering & Information Technologies, The University of Sydney{sayaka, itot} ＠ itolab.is.ocha.ac.jp, masa.takatsuka ＠ sydney.edu.au
概要
我々の身の回りには数多くの時系列データが存在しており，その分析において，可視化は有効である．代表的な時系列データの可視化手法として，折れ線グラフとヒートマップがあるが，いずれの手法もデータの部分的な特徴と全体のクラスタの変遷の両方を同時に観察するのは難しい．本論文では実数値型の時系列データを対象に，要素間の交流や共起性の表現に用いられる storyline による可視化手法を提案する．本手法ではまず，実数値の類似度にもとづいて局所時間帯ごとに要素をクラスタリングする．続いて隣接する時間帯間にてクラスタ群のマッチングを適用し，その結果に沿ってクラスタの配置を決定する．
この配置結果にしたがって各要素を描くことで，類似した要素同士が画面内で近接するような storyline を描く．また，各要素の実数値変化の重要度にもとづいて要素の不透明度を設定することで，重要度の高い部分を強調し，ユーザの知覚効果の高い可視化結果を得る．さらに本手法では，storyline 上で興味深い局所を対話的に指定することで，選択した要素群の実数値
変化を別の画面領域に表示するというインタラクション機能を持つ．本手法を用いてクラスタの時間的変遷を観察することにより，長時間にわたって類似した値/外れ値を有する要素だけでなく，所属するクラスタが頻繁に切り替わるような要素に着目しやすくなる．また，インタラクション機能を用いて元の数値情報を参照することで，クラスタの切り替えの要因を分析することが可能である．

Abstract
Information visualization is an eﬀective approach to analyze time-varying data in our daily lives. We commonly
represent time-varying values applying polyline charts or heatmaps; however, it is diﬃcult to simultaneously observe
short-term features of time-varying values and cluster transitions while applying either polyline charts or heatmaps.
This paper proposes a storyline-based visualization technique for consecutive numerical time-varying data. Storyline
is a visualization technique to show associative features among elements over time. Our technique ﬁrst measuressimilarity of elements in each time-step, and divides the elements into clusters. The technique then deﬁnes thecluster layout by matching corresponding clusters between two adjacent time-steps, and draws similar elements asproximity storyline. Reﬂecting transparency on storyline as a visual variable, the technique also emphasizes theamount of line changes. Moreover, the technique provides a user interface so that users can interactively selectinteresting parts on storyline, and explore the numerical values by observing a polyline-based visualization. We
believe it is important to focus on elements which switch their clusters frequently. We suppose that by making theappearances of numerical changes prominent based on the amount of changes, a user would be able to eﬀortlessly pay
his/her attention to where those changes occurred. This easy recognition of numerical changes would lead to further
focused investigation on the causality through examination of the original numerical values and other associatedinformation.

– 188 –

芸術科学会論文誌 Vol. 14, No. 5, pp. 188 – 197ている．Sankey Diagram はクラスタ間の流量変化を観

1 はじめに察するのに適している．一方，storyline は個々の要素

ビッグデータ時代の到来に伴い，多種多様なデータについて，クラスタから別のクラスタへの移動を観察す

ベースに時系列性を有する情報が蓄積されるようになっるのに向いている．Storyline はもともと XKCD といた．これらの時系列データの分析・予測はデータマイニ

うウェブコミックに手書きイラスト “Movie Narrative

ングや知覚情報処理等の分野においては既に数多くの研Charts”として紹介されたものであり [5]，物語の登場人

究が進められている．一方で，時系列データの観察と理物の共起行動を時系列のダイアグラムで表現している．

解，またそれを扱う業務での能動的な意思決定を支援す

図 1 にあるように，x 軸が時刻，各曲線が映画の各登場る道具として，情報可視化が有用である．情報可視化に人物を表しており，複数の登場人物が同一のシーンに登は，目で見て全体像を把握し，対話操作によってデータ

場した場合に線が近接するような可視化結果となる．こ中の必要な部分を能動的に選択できる点に加え，特に時れらの可視化手法は任意の 2 要素間の流量や共起に関系列データにおいてはリアルタイムな観察のための瞬時する時間変化を表現するのに用いられるのが一般的であな理解においても有効である．一般的に情報可視化手法

り，各要素が独立に実数値を有するような時系列データ

は，一次元，二次元，三次元，多次元，木構造，グラフ，に対してこれらの可視化手法を適用した事例はまだ少な

時系列の 7 種類のデータ構造に大別される [1]．その中

く，議論の余地が残っている．

でも我々はこれまで，時系列データの可視化に関する新
しい手法の開発に取り組んできた [2]．
多量の実数値で構成される時系列データの可視化では，データを構成する実数値の大小そのものだけでなく，実数値どうしの関連性や変化量などを観察することも重
図 1 手書きによる storyline の例（[5] より転載）．

要である．例えば多数の商品の売上を一定時刻ごとに集計して作成される時系列データでは，商品単体の売上の推移だけでなく，どの商品と同時期に売れるか，季節や以上を踏まえて，我々は実数値で構成される時系列

時間帯との相関はあるか，といった関連性を観察するこ

データを storyline として表現する可視化の一手法を提とによって，より深い分析が可能になり，売上予測や予案する．本手法ではまず，時系列データを構成する各要兆発見などに貢献できると考えられる．

素に対して局所時間帯ごとにクラスタリングを適用し，時系列データの可視化では折れ線グラフとヒートマッ

隣りあう局所時間帯間のクラスタ群にマッチングを適用

プの 2 種類の視覚表現が特に多用されている [3]．しかする．その結果にもとづいて，クラスタに対応する長方

しいずれの視覚表現においても，大量のデータ要素（以形群を画面配置し，長方形群を通過する線分として各要

下「要素」と称する）の相互関係を全て表現するのは難素を描画する．描画時には次を工夫することで，より深しい．例えば折れ線グラフでは，要素に対応する折れ線

い分析を可能にする．

を 1 個の画面空間に大量に描くことで，その相互関係に関する視認性は大きく低下する．そこでクラスタリング• 各要素に割り当てられた数値やメタ情報から線分のなどを適用して要素を構造化した上で可視化することが色相を算出することで，要素のクラスタ移動の要因有効である．ここで時系列データにおいては要素間の相を理解しやすくする．

互関係も時間とともに変遷するため，クラスタ構造の時• 各要素の数値変化の重要度から線分の不透明度を算

間的変遷を表現することも重要となる．一方でこれらの出することで，重要な数値変化を有する線分に視線可視化手法において，データの部分的な特徴とクラスタを届きやすくする．

構造の時間的変遷を同時に表現するのは簡単ではない．

本手法を用いてクラスタの時間的変遷を観察すること

クラスタ構造の時間的遷移の表現に向いた可視化手法

で，次の 3 種類の要素群を視覚的に分離して理解するのとして，Sankey Diagram と storyline がよく利用されが容易になる．

– 189 –

芸術科学会論文誌 Vol. 14, No. 5, pp. 188 – 197• 長時間にわたって同一クラスタに属する要素群ることで，データ大域にわたるクラスタ構成と，局所に• 外れ値を有するために小さなクラスタに属する要見られる特徴的な現象を同時に表現する手法を提案した[4]．しかしこのようなアプローチでは，局所時間帯ごと

素群• ある時には特定のクラスタに属し，またある時には別のクラスタに属する，というようなふらふらとしにクラスタリングを適用できない．

2.2

た時間的遷移を有する要素群

Storyline を適用した可視化手法

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

例えば，夏には商品 A と B の売上の時間的推移が類似
して，冬には商品 A と C の売上の時間的推移が類似している，という商品データがあるとする．図 2（右）のように本手法でこれを可視化することで，商品 A が夏から冬にかけて別のクラスタに移動していることが storyline
上で明確に表現できる．これによって，商品 A に対する
店舗での陳列には注意が必要である，というような知識の発見を期待できる．

を可視化する，という問題設定は実数値型以外の時系列
データにおいて活発に議論されている．人物や単語を要素とみなしたとき，その 2 要素間の距離（例えば共起性）の時間変化を可視化する，という問題がその最たる例である．この課題においてフレキシブルに各要素間の共通性を描く手法として，XKCD の手描きイラスト [5] に着想を得て storyline のメタファを用いた時系列データ
可視化手法が近年数多く提案されている．Ogawa らは
ソフトウェア開発のリポジトリを題材として storylineのレイアウトを自動生成する手法を提案した [6]．この手法による可視化結果は，小〜中規模のプロジェクトによるソフトウェア開発の構造変化の観察に有用である．
一方でこの可視化結果には，線分の小刻みな蛇行や線分

図2

どうしの重なりが数多く含まれており，時間軸に沿って商品の売上データの可視化例の比較．

各線分の動きを追うことは困難である．Tanahashi らは

Ogawa らの手法を拡張し，遺伝的アルゴリズムを用いて

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

線分の配置を最適化する手法を提案した [7]．この手法

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

では，線分のずれや重なり，画面上の空白を制御するこ

実数値型時系列データの可視化には旧来から，折れ線とで，Ogawa らの手法の問題点を解決し，線分間の共起グラフにもとづく可視化手法と，ヒートマップにもとづ性に着目しやすい結果を実現している．しかしこれらのく可視化手法が特に多用されてきた [3]．

手法にはスケーラビリティと計算時間に課題がある．数折れ線グラフにおいては，要素数の多いデータにおけ百以上の要素数を有するデータにおいて線分間の重なりるクラッタリングがその可読性を妨げてきた．その解決が非常に多くなり，各要素を目で追うことが困難になる．

策として折れ線のサンプリングが有効であり，有効なサ

この問題を解決するため，StoryFlow は各クラスタにンプリングを実現するために時系列データのクラスタリ

詳細度制御を適用し，Storyline 全体の特徴を捉えやす

ングがしばしば適用されてきた．我々は折れ線グラフにくした [8]．StoryFlow では各要素の階層構造を考慮しもとづく可視化手法において，局所時間帯ごとに要素をた配置に加え，インタラクティブな線分の並べ替え，バ

クラスタリングしてサンプリングすることにより，可読

ンドリングなどの機能も実現している．また，Reda ら性が高く，かつデータの特徴を見逃さないような可視化

は storyline のメタファを用いつつも，数百規模の要素を実現した [2]．しかしこのようなアプローチでは，どのの描画が可能な手法を提案している [9]．Reda らの手法

要素が局所時間帯ごとにいつクラスタを変遷したか，とでは，ｙ軸上にコミュニティを示すスロットを固定し，いった情報を視認することができない．

各スロット内に線分を配置している．このアプローチはまたヒートマップにもとづく可視化手法において，時ユーザが特定のクラスタ間の変化に着目したい場合，有間帯全体にわたる類似度で要素を並べ替え，さらに局所効であるといえる．しかしこれらの手法は，上述のよう時間帯ごとに見られる特徴的な数値パターンをマークすに既知の関係性やキーワード・トピック等を入力情報と– 190 –

芸術科学会論文誌 Vol. 14, No. 5, pp. 188 – 197

して要素の共通性を表現しており，実数値で構成される有しているように見える．しかし多くの既存手法におい時系列データを対象にしたものではない．また，要素間て，曲線の画面上での移動量は，個々のクラスタに含まの距離などにもとづく強調表示などは行われていない．

れる曲線の本数やクラスタの配置などによって決定され実数値型の時系列データにおいても，局所時間帯ごとるため，必ずしも実数値の変化量に比例するとは限らなのクラスタリング結果に焦点をおいた可視化手法は最近

い．この問題を解決する手段として，曲線の形状以外のになっていくつか報告されている．例として Lex ら [10]視覚変数によって重要な要素を強調表示することが有効

は，クラスタごと，局所時間帯ごとに断片化されたヒーである．

トマップを Sankey Diagram 風に接続する可視化手法を提案している．しかし Sankey Diagram と同様な表現を採用したこの手法では，時系列データを構成する個々の要素の変遷に着目することが難しい．また Turkay ら[11] は，局所時間帯ごとにクラスタを縦に並べ，各クラ
スタを通過するように要素を曲線で描く可視化手法を提案している．この手法は storyline の簡易な生成方法の一種と考えることもできる．本論文の提案手法における描画形式は Turkay らの描画形式に類似しているが，次節にて議論する点において提案手法は Turkay らの手法
図3を改良しているといえる．

2.3

画面上での変化量が大きく見える storyline の例（文献 [7] より転載）．

Storyline の適用に関する注意点

本節では，Storyline を実数値型時系列データに適用する際の注意点と，それらを改良するためのアプローチを

数値表示へのインタラクション: storyline による描画は議論する．

数値自体の変遷を表現しない．そのため，既存手法ではクラスタの位置関係: storyline による描画では，要素をヒートマップを併用する数値表現 [10] が用いられている表現する曲線の縦方向の移動を小さくすることが重要でが，厳密に数値を読み取ることは難しい．この問題を解ある．本手法では隣りあう局所時間帯において多くの要決する一手段として，storyline 上で興味深い局所を対話

素を共有するクラスタが画面上で隣接するように配置す的に指定し，別の画面領域で数値表示することが考えらることで，曲線の移動を抑える．上述の Turkay らの手れる．

法 [11] においても，要素の共有性を考慮した並び替えが3 提案手法行われている．一方で storyline による描画を見ただけ
では，実数値型時系列データの数値自体の変遷を理解す

本章では storyline を適用した実数値型時系列データ

ることは難しい．そこで，要素の共有性に加えて数値変可視化手法を提案し，その各処理について論じる．本手化も加味してクラスタを配置することが有用である．

法の処理の流れは下記の通りである．

重要な要素の強調表示: storyline では多くの場合におい1. 入力データが対象とする時間帯を「局所時間帯」とて，画面上で大きく移動する曲線が目立って見えるため，呼ばれる短い時間帯に分割し，各局所時間帯内の時この曲線に対応する要素が重要な要素（あるいは大きな

系列実数値にクラスタリングを適用する．

変化を有する要素）であると連想しがちである．よって

2. 各局所時間帯を構成するクラスタ群を長方形とみな実数値の時系列データを storyline で可視化すると，実し，画面空間に配置する．

数値の変化の小さい要素に対応する曲線が画面上で大3. 入力データを構成する各要素を，クラスタに対応すきく移動しているために目立ってしまう，というようなる長方形群を通過する storyline として描画する．

直感に合わない可視化結果を生むことがある．例えば図その際に，各要素のクラスタ移動が強調されるよう3 において，左下部の赤の曲線は短時間で急激な変化をに storyline の透明度を自動調節する．

– 191 –

芸術科学会論文誌 Vol. 14, No. 5, pp. 188 – 197

図4

本手法の処理の流れ．

• XKCD のような表現を実現するために，曲線として

図 4 に 6 個の要素集合 (A〜F) をもつ時系列データを用の各要素の配置を最適化する方法 [6, 7]

いた，処理の流れの概要を示す．

• まず要素にクラスタリングを適用し，続いてクラス

なお本章では，時系列データは m 個の要素の集合

D = {a1 , a2 , ..., am } であり，要素 ai は n 個の各時刻における実数値の集合 ai = {ai1 , ai2 , ..., ain } であるとする．また時系列データが対象とする時間帯は 2 時刻を区間とする Nw 個の局所時間帯に分割されているとする．
また，入力データの対象となる時間帯を Nw 個に区切った局所時間帯のセットを W = {w1 , ..., wNw } とし，i 番目の局所時間帯 wi を次の通り記述する：

wi = {tstart , tend , Ci }， ただし Ci = {ci1 , ci2 ..., cik }．
ここで，tstart ，tend は wi の開始時間および終了時間であり，k は局所時間帯 wi におけるクラスタ数，cij は j
番目のクラスタを示す．

タの位置を決定する方法 [9, 10, 11]
本手法では数百・数千の要素を有する大規模な時系列
データを対象とするため，前者の適用は非現実である．
そこで本手法では後者を適用するが，本手法では時間的に隣接する局所時間帯間でクラスタのマッチングを取る必要がある．そこで本手法では，以下の手順でクラスタの共通度を求め，それに沿って隣接する局所時間帯間のクラスタをマッチングする．

3.2.1 クラスタの共通度の算出
本手法では，i 番目および (i + 1) 番目の局所時間帯を構成するクラスタ集合 Ci ，Ci+1 に対して，以下の式 (1)3.1 各局所時間帯におけるクラスタリング
まず時系列データを構成する各局所時間帯に対して，
要素のクラスタリングを適用する．ある局所時間帯が j
番目から k 番目 (j < k) までの nw 個の時刻で構成されている場合に，本処理では i 番目の要素 ai から実数
値 [aij , ..., aik ] を抽出し，これを nw 次元ベクトルとしが最大となるような最大マッチング問題を適用する．ここで，cij ∩ c(i+1)k はマッチングされる 2 クラスタが共有する要素数，ncij は cij に含まれる要素数，acij は cijを構成する実数値の平均，α は 0 ≤ α ≤ 1 を満たす定数である．

て扱う．このベクトル群に対してクラスタリングを適用することで，各局所時間帯における要素を有限個のクラスタに分類する．現時点での我々の実装では非階層型ク∑(

ラスタリングを適用しているが，階層型クラスタリング
を含め他のクラスタリング手法を適用することも可能で
ある．

)

(1)
この式で第 1 項は両クラスタ間の要素の共有性を表すものであり，第 2 項は両クラスタ間の数値の不変性を表す

3.2 クラスタの画面配置ものである．両者を考慮して最大マッチング問題を解く

Storyline を適用した既存の可視化手法の処理手順は，次の 2 種類に大別される．

(cij ∩ c(i+1)k )1+ (1 − α)αncij + nc(i+1)kacij − ac(i+1)kことで，要素を表現する曲線の移動量を抑え，かつクラ
スタ配置結果から数値変化を説明しやすい可視化結果を得る．

– 192 –

芸術科学会論文誌 Vol. 14, No. 5, pp. 188 – 197

3.2.2 クラスタのマッチングと配置決定

通過する要素群の実数値の時間変化を別の画面領域に表本手法では各局所時間帯を構成する各クラスタを長方示するインタラクション機能を有する．具体的には，図形で表現し，3.2.1 節の処理によってマッチングされたク

5 のように画面の上半分に storyline を描き，上述の長方ラスタが画面上で左右に隣接するように，以下の処理に

形描画操作によって要素群を選択すると，折れ線グラフ

よって各長方形の位置を決定する．

によってその要素群を画面の下半分に描く，という連携

1. 画面左端の局所時間帯を構成するクラスタを上から可視化を実装している．

順に配置する．

2. 左から i 番目の局所時間帯について配置が終わったら，左から (i + 1) 番目の局所時間帯について，マッ
チングされたクラスタができるだけ画面上で左右に並ぶように各クラスタの位置を決定する．

3.3

Storyline の描画

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

図5

本手法におけるインタラクション機能．

• i 番目の局所時間帯にて当該要素が属するク
ラ ス タ cij の 左 右 端 お よ び 上 下 端 の 座 標 値 を

3.4

Storyline の不透明度算出[xminij , xmaxij ], [yminij , ymaxij ] とする．この

前節で論じた通り，現時点での我々の実装では，次の 2とき当該要素を表現するために 2 点 (xminij , yij ),種類の基準を各時刻における重要度とみなし，storyline(xmaxij , yij ) を結ぶ線分を描く．ただし yminij ≤の不透明度算出に用いている．

yij ≤ ymaxij とする．

隣接する局所時間帯間 においては，要素の実数値 aij の• i 番 目 の 局 所 時 間 帯 と (i + 1) 番 目 の 局 所

微分値を基準にする．具体的には，隣接する局所時間帯

時 間 帯 を 連 結 す る た め に ，2 点 (xmaxij , yij ),間でクラスタを移動する曲線について，対応する要素の(xmin(i+1)j , y(i+1)j ) を結ぶ線分を描く．

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

このとき我々の実装では，storyline の色を以下のよを示す．

うに定義する．本手法では各要素に 1 次元の属性が付

1 つの局所時間帯内 では，要素の安定性を基準にする．

与されていると仮定し，その属性値に対応した色相を線

具体的には，i 番目の局所時間帯について，(i − 1) 番目分に与える．1 次元の属性がない場合には例えば，各要の対応するクラスタと共通する要素は不透明度を高く，素の実数値 aij で線分の各頂点に色相を与えることも可

別のクラスタから移動してきた要素は不透明度を低く設能である．また各要素の数値変化の重要度に応じて，線定する．

分を構成する各頂点に不透明度を与える．重要度が大き

不透明度にもとづく視覚表現を採用した理由は以下のい部位において不透明度を高く設定することで，観察者通りである．人間の目は知覚的に変化の激しい部分に対の注意を引くことができる．現時点の我々の実装では実して反応する [12] という Marr の理論を適用して，我々

数値 aij の微分値および各クラスタ内での要素の安定性

は storyline 上での要素のクラスタ移動という重要な現にもとづいて重要度を算出している．詳細は次節にて論象が発生している部位に視覚変化の大きな表現を適用しじる．

たいと考えた．可視化における視覚表現の手段として，また各要素のクラスタ移動の要因を分析するために，

Bertin [13] によって提案された視覚変数集合が知られstoryline 中の局所における実数値の変化を観察するこている．それに対して MacEachren は，画面内の同一とが有用である．我々の実装では，storyline を描画した

座標上に複数の情報を載せる場合に不透明度を用いる

画面領域でマウス操作で長方形を描くと，その長方形をのが有効であると述べている [14]．本手法においても，– 193 –

芸術科学会論文誌 Vol. 14, No. 5, pp. 188 – 197

storyline を構成する線分において重要度の高い部位にのみ高い不透明度を付与することで，既に用いられている他の視覚変数（線分の形状や色相）に干渉することなく，ユーザの注意を重要度の高い部位に向けられると考える．また，不透明度から定量的な差異を読み取ることは困難であるものの，ユーザは不透明度を見て直感的に
図6

その部位における重要度を判断可能である．

曲線の色と地方の対応．

4 実行結果4.1 実行環境と処理時間
我々は JDK (Java Development Kit) 1.7.0 を用いて提案手法を実装し，Windows 8.1 (64bit) 搭載のラップ
トップ PC (CPU 1.80 GHz, RAM 4.0GB) 上で実行した．表 1 は上記の環境におけるクラスタリング，クラ
図 7 12 観測所の気温データを適用した実行例．

スタの画面配置およびインタラクションの平均処理時間を示す．なお，いずれの事例も元データ内の時間帯数は240，また，クラスタ数を 3，局所時間帯サイズを 12 時

れぞれ 1 本ずつ不透明度が高く表示されている．これは間に設定した．

2 番目の局所時間ではクラスタ 1 に属する当該線分が，3

表1

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

番目の局所時間帯でクラスタ 2 に切り替わっているため

要素数

クラスタリング

画面配置

インタラクション

といえる．

86

0.40

0.08

0.03

また本章では，3.3 節で述べた要素選択機能を用いた

376

1.83

0.19

0.18

折れ線表示についても示す．折れ線の表示においては，

x 軸が時刻，y 軸が気温，各折れ線が各観測所を示す．
上記の結果から，既存手法 [6]，[7] に比べ，本手法は対4.3

2009 年 4 月の事例

図 8（上）は 2009 年 4 月の 1 カ月全体の気温変化を表話的に操作する上で十分高速であるといえる．

4.2 適用データ

す storyline である．上のクラスタから順に，青で示さ

本章では，全国 376 箇所の観測所について 3 時間ごとれた北海道，青緑から黄緑で示される東北から中部，黄に観測された 2 種類のアメダス気温データを提案手法に色から赤で示される近畿・中国から沖縄・南西諸島に属適用した事例を紹介する．可視化結果において，x 軸はする観測所を多く含んでおり，ところどころ変動が見ら

時刻，1 本の曲線は 1 観測所を示し，近接する曲線は同じれるものの，全体的に似たような地方が安定して同じククラスタに属することを表す．この事例ではクラスタ数

ラスタに配置されている．北海道の観測所を多く含むクを 3，局所時間帯サイズを 12 時間ごと，各クラスタ内のラスタについて，前半（wa ）と月末（wb ）にそれぞれ破曲線の並びは観測所 ID 順とした．また観測所 ID から線で囲んだ 1 局所時間帯分のみ，青から橙で示させる北観測所を 8 つの地方に分類し，ColorBrewer 2.0[15] の海道から九州にかけての観測所が含まれており，全国的配色方針に沿って各地方に色を割り当てた（図 6 参照）．

に同じような気温の変化をしていたと推察できる．これ

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

らの局所時間帯について，長方形による要素選択機能を観測所の気温データを適用した実行結果の一例を示す．

用いて，各クラスタに属する要素群を折れ線で描いた結3 番目の局所時間帯において，各クラスタには上から順

果を図 8（中）（下）に示す．wa ，wb いずれにおいてもに，東北地方，東北から九州地方，北海道地方に属する

他の時間帯に比べ，全国的に似たような気温変化を示し観測所が含まれている．クラスタ 2 に着目すると，東北ており，気温差が小さくなっていることが確認できる．

地方を示す青緑の線分と中部地方を示す黄緑の線分がそ

局所時間帯 wb についてさらに観察するため，月末の– 194 –

芸術科学会論文誌 Vol. 14, No. 5, pp. 188 – 197

図82009 年 4 月の気温変化：（上）storyline による表示．局所時間帯 wa （中）/wb （下）における，選択

図 9 2009 年 4 月末 1 週間分の気温変化：（上）sto-

したクラスタに属する要素の折れ線表示．

ryline による表示．（中）クラスタ cb1 /（下）クラスタ
cb3 に属する要素の折れ線表示．

1 週間分を拡大した結果を図 9 に示す．storyline 上の破線で囲まれた局所時間帯 wb について，青から橙で示された要素が含まれるクラスタが 2 種類存在することが観察できる．クラスタ cb1 ，cb3 それぞれについて要素選択機能を用いて折れ線を表示した例を図 9（中）（下）に示

図 10 2009 年 4 月の気温変化を折れ線によって可視す．各クラスタに含まれる折れ線群は似たような気温変化した結果．

化を示しているものの，平均気温や気温差が異なってい
ることがわかる．これは日本海側と太平洋側で異なる気

4.4

2009 年 11 月の事例

温変化を示したためと推察できる．選択した折れ線群に

図 11（上）は 2009 年 11 月の上旬 2 週間の気温変化を着目すると，いずれのクラスタにおいても，選択した局表す storyline である．y 軸上中央に配置されたクラス

所時間帯以前は低温で推移していた折れ線群が，この時タは比較的安定して推移しているが，図 11（上）の下部

間帯だけ気温が高くなったため，全国的に同じような気に配置されたクラスタに着目すると，期間の前半は，青温変化となったと考えられる．図 9（下）に示されるク

緑から黄色で示される東北から中部地方の観測所を多く

ラスタ cb3 内の要素については特にその傾向が顕著であ含んでいるが，破線で囲まれた局所時間帯において，青り，storyline 上でも青や青緑で表された北海道や東北地で示された北海道の観測所ばかりを含む表示に切り替わ

方の観測所の要素が不透明度が高く表示されている．こることが観察できる．この局所時間帯の前後で他のクラのことから，直接数値変化を示していない storyline にスタも観察してみると，前半は北海道の観測所が複数のよる表示からでも，これらの地方で他の時間帯とは違うクラスタに分散しているが，切り替え後は一番下のクラ特徴的な気温変化があったと推察できる．

スタに集約されており，その後の時間も不透明度が低い

同じデータについて，折れ線によって直接全観測所のことから，そのまま安定して推移していることがわかる．

気温を表示した例を図 10 に示す．期間全体を通して似

このことから，11 月の 2 週目以降，北海道で他の地方とたような気温変化を示しているため，どの部分に着目すは極端に異なる気温変化を示していたと推察できる．

べきかが判断し難く，この図のみから上記のような知見

前述の局所時間帯の前後 2 箇所を選択し，該当する折を得ることは困難である．この結果から，storyline にれ線を表示させた結果を図 11（下），12（下）に示す．図よって全体の傾向を俯瞰し，細かな数値変化を見たい場

11（下）は切り替えの前後 2 局所時間帯分，図 12（下）合は折れ線グラフによる結果と比較するというアプロー

は切り替え後，青の要素で安定して推移している 9 局所チが有効であるといえる．

時間帯分を選択したものである．図 11（下）について，– 195 –

芸術科学会論文誌 Vol. 14, No. 5, pp. 188 – 197

選択した局所時間帯の前数日間は折れ線群全体が似たよい可視化結果を得る．

うな気温変化を示しているのに対し，選択部分以降は数

今後の課題として，気象以外のデータでの適用事例を日間にわたって気温の開きが顕著であることが観察でき

示し，ユーザテスト等によって可視化結果の妥当性をる．図 12（下）の安定したクラスタを選択した結果と比検証することを計画中である．また手法部分の拡張とし較してみると，このクラスタに含まれる折れ線群（北海て，各要素の実数値間の類似度を算出する単位となる局

道の観測所群）が当該期間において，極端に低い気温変所時間帯を可変にすることで，よりデータの局所的特徴化を示していることがわかる．

に応じたクラスタ切り替えの数や時間を検出したいと考これらの結果から，storyline 上のクラスタの安定性かえている．この点に関して，不等間隔離散化を適用した

ら，気温変化が他と異なる時間帯を推測できるといえる．

パターン抽出手法 [16] などを用いることが考えられる．
さらに，3.1 節でも示したように，階層型クラスタリン
グ手法やクラスタ数可変型のクラスタリング手法などの適用を検討したい．

参考文献[1] B. Shneiderman. The eyes have it: A task by data
type taxonomy for information visualizations. Visual Languages, 1996. Proceedings., IEEE Sym図 11 2009 年 11 月上旬 2 週間分の気温変化：（上）posium on, 336–343, 1996.

storyline による表示．（下）選択した 2 局所時間帯に[2] S. Yagi, Y. Uchida, T. Itoh. A polyline-based

含まれる要素の折れ線表示．

visualization technique for tagged time-varying
data. 16th International Conference on Information Visualisation, 106–111, 2012.
[3] 伊藤貴之. 視覚協創学 (5) マッピング技術: 特に時
系列データ可視化技術の体系化に向けて. 日本計算工学会計算工学講演会, 17, 4, 2012.

[4] M. Imoto, T. Itoh. A 3D visualization techniquefor large scale time-varying data. 14th International Conference on Information Visualisation,
図 12 2009 年 11 月上旬 2 週間分の気温変化：（上）17-22, 2010.

storyline による表示．（下）選択した 9 局所時間帯に[5] R. Munroe. Movie narrative charts. http://

含まれる要素の折れ線表示．

xkcd.com/657/, accessed 4 June, 2015.
[6] M. Ogawa, K.-L. Ma. Software evolution story-

5 まとめと今後の課題

lines. ACM Iternational Symposium on Software
visualization, 35–42. 2010.

本論文では，storyline を用いた実数値型時系列データの可視化手法を提案した．本手法では局所時間ごとに時[7] Y. Tanahashi, K.-L. Ma. Design considerations
for optimizing storyline visualizations.

系列実数値情報の類似度を算出し，隣接する局所時間帯

Transactions on Visualization and Computer間のクラスタ群にマッチングを適用することでクラスタ

Graphics, 18(12), 2679–2688, 2012.

の配置を決定し，類似した要素同士が画面内で近接する
ような storyline を描く．また，本手法では各要素の数

IEEE[8] S. Liu, Y. Wu, E. Wei, M. Liu, Y. Liu. Storyflow:

値変化の重要度にもとづいて要素の不透明度を設定する
ことで，重要度の高い部分を強調し，ユーザの知覚に近– 196 –

Tracking the evolution of stories. IEEE Transactions on Visualization and Computer Graphics,

芸術科学会論文誌 Vol. 14, No. 5, pp. 188 – 197

19(12), 2436–2445, 2013.

年お茶の水女子大学大学院人間文化創成科学研究科理[9] K. Reda, C. Tantipathananandh, A. Johnson, J.
Leigh, T. Berger-Wolf.

Visualizing the evolu-tion of community structures in dynamic social
networks. Computer Graphics Forum, 30, 1061–
1070, 2011.

学専攻博士前期課程修了．現在，お茶の水女子大学大学
院人間文化創成科学研究科理学専攻博士後期課程に在
学中．
伊藤

貴之

[10] A. Lex, C. Partl, D. Schmalstieg. Comparative
analysis of multidimensional, quantitative data.
IEEE Transactions on Visualization and Computer Graphics, 16(6), 1027-1035, 2010.
[11] C. Turkay, J. Parulek, N. Reuter, H. Hauser. Interactive visual analysis of temporal cluster structures. Eurographics/IEEE Symposium on Visu-

1990 年早稲田大学理工学部電子通信学科卒業．1992 年

alization (EuroVis2011), 711-720, 2011.

早稲田大学大学院理工学研究科電気工学専攻修士課程修[12] D. Marr. Vision: A computational investigation

了．同年日本アイ・ビー・エム (株) 入社．1997 年博士

into the human representation and processing of(工学)．2000 年米国カーネギーメロン大学客員研究員．

visual information. New York: Henry Holt and2003 年から 2005 年まで京都大学大学院情報学研究科

Co., Inc, 1982.

COE 研究員（客員助教授担当）．2005 年日本アイ・ビー・[13] J. Bertin. Semiology of graphics: diagrams, net-

エム (株) 退職，お茶の水女子大学理学部情報科学科助教

works, maps. University of Wisconsin press, 1983.

授．2011 年同大学教授，同大学シミュレーション科学教[14] A. M. McEachren. How maps work: represen-

育研究センター長，ACM, IEEE Computer Society, 芸

tation, visualization, and design. The Guilford
Press, New York, 1995.

術科学会，他会員．
高塚正浩[15] M. Harrower, C. A. Brewer. Colorbrewer 2.0.
http://colorbrewer2.org/, accessed February,
2015.
[16] M. Ohsaki, H. Abe, T. Yamaguchi. Numericaltime-series pattern extraction based on irregularpiecewise aggregate approximation and gradientspecification. New Generation Computing, 25(3),213–222, Springer, 2007.

Masahiro Takatsuka received his Ph.D in Elecrical andComputer Systems Engineering from Monash University in 1997. He is Director of the High-Performance

八木 佐也香Visualization Research Laboratory (ViSLAB) at the
University of Sydney. He has been active researcher
in the area of Scientific Visualization and Visual Analytics since 1998. He was a Senior Research Faculty at
GeoVISTA Centre, the Pennsylvania State University
between 1998 and 2002. He joined the University of2011 年お茶の水女子大学理学部情報科学科卒業．2013

Sydney in 2002 and was also a researcher at National
ICT Australia between 2004 and 2007.

– 197 –