======= distance2-east-japan-earthquake-meiji-sanriku-cumulative-human-damage === ==== 東日本大震災と明治三陸津波の ==== 【震源からの距離】と【人的被害の累和】との合成比較グラフ ================================================================================== $$u // ユーザデータ・セクション ---------------------　東日本大震災分 ------------------------ $c 0001.00,0043.00,aaa // ケース始点,終点番号 ,★第１系列名 緯度 ,bbb // 空白で同一ケース範囲,★第２系列名 経度 ,ccc // 空白で同一ケース範囲,★第３系列名 浸水域人口 ,ddd // 空白で同一ケース範囲,★第４系列名 人的被害 $d // データ入力指示コマンド ctype // ケース毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ 40.5122 141.4883 5229 2 40.6833 141.3688 1924 2 40.5992 141.3978 3820 0 40.4525 141.6211 1189 0 39.6414 141.9572 18378 901 39.0823 141.7083 19073 463 40.1903 141.7753 7171 4 39.0150 141.6295 16640 2203 39.2759 141.8856 13164 1355 39.3599 141.9066 11915 1711 39.4675 141.9489 11418 938 39.8431 141.7968 1137 7 39.9303 141.8889 1582 36 40.0053 141.8861 1115 1 40.1103 141.8179 3177 38 40.4086 141.7187 2733 0 38.2681 140.8697 29962 860 38.4342 141.3028 112276 5734 38.3145 141.0222 18718 22 38.9083 141.5700 40331 1533 38.1717 140.8919 12155 1052 38.2940 141.0042 17144 188 38.1045 140.8699 8051 186 38.4264 141.2104 34014 1769 38.0378 140.8527 14080 273 37.9625 140.8778 8990 749 38.3814 141.0692 4053 4 38.3047 141.0594 9149 76 38.3303 140.9756 542 4 38.4455 141.4444 8048 1156 38.6778 141.4464 14389 1173 37.0506 140.8878 32520 383 37.7967 140.9197 10436 469 37.6422 140.9571 13377 938 37.2145 140.9947 1385 3 37.2825 140.9935 1746 13 37.3456 141.0086 1401 11 37.4044 140.9835 1127 44 37.4492 141.0125 1278 35 37.4947 141.0008 3356 186 37.8763 140.9195 4666 115 36.8020 140.7511 7212 6 35.7203 140.6467 8303 15 ---------------------　明治三陸津波分 ------------------------ $c // クロスセクションデータ属性コマンド 0001.00,0043.00,CCC // ケース始点,終点番号 ,★第５系列名　人口 ,DDD // 空白で同一ケース範囲,★第６系列名　人的被害 $d // データ入力指示コマンド ctype // 変量毎に読むタイプ ------- ユーザ自身が文字・数値データを ------- テキストファイルまたはExcelシートからコピー＆ペーストする． ------- ユーザデータの各行の末尾にも「//」を挿入してコメント文を記述できる． ------- ★ユーザ文字・数値データをこの行直後にペーストする ★↓↓★ #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A 18093 3071 11995 3174 4404 487 4257 818 12665 6477 6983 600 9844 2119 820 364 529 232 925 302 1352 260 4247 254 #N/A #N/A 7875 271 #N/A #N/A 12178 1887 #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A 2521 1 7148 1234 #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A ====================== $$v // 変量分析セクション $a // 変量記号の割り当て a,aaa // 緯度 b,bbb // 経度 c,ccc // 浸水域人口 東日本大震災 d,ddd // 人的被害 東日本大震災 C,CCC // 人口 明治三陸津波 D,DDD // 人的被害 明治三陸津波 ---------------- $d // 表示範囲 all // 全範囲 ---------------- $t // 変数変換 r=(d/c*100) // 人的被害人口比（％）東日本大震災 R=(D/C*100) // 人的被害人口比（％）明治三陸津波 p=:ci(a) // 個体識別文字列p作成 P=:ci(a) // 個体識別文字列P作成 =pr*(a,b,c,d,r,C,D,R,P) // 数値プリント ------------------------------------------------------------ ---- 起点（震源）から各終点までの距離 東日本大震災 ---- U=(38.1) // ★東日本大震災の震源の緯度のスカラー V=(142.9) // ★東日本大震災の震源の経度のスカラー ((R(( // 繰り返しルーティン R の開始 A=abs(a-U) // 各市の緯度の震源からの乖離 B=abs(b-V) // 各市の経度の震源からの乖離 ..................... 起点（震源）から各終点(各市町村庁）までの距離の計測処理 n=cos(U/180*3.14159) // 起点緯度（ラジアン）のcos l=cos(a/180*3.14159) // 終点緯度（ラジアン）のcos m=cos((U+a)/2/180*3.14159) // 起点終点の平均緯度（ラジアン）のcos　 M=(1+((6378.137**2-6356.752**2)/6356.752**2)*(m**2))**0.5 M=(6378.137**2/6356.752)/(M**3) // 起点終点の平均緯度の子午線曲率半径 N=(1+((6378.137**2-6356.752**2)/6356.752**2)*(n**2))**0.5 N=(6378.137**2/6356.752)/N // 起点緯度の卯酉線曲率半径 L=(1+((6378.137**2-6356.752**2)/6356.752**2)*(l**2))**0.5 L=(6378.137**2/6356.752)/L // 終点緯度の卯酉線曲率半径 =pr*(M,N,L) // 数値プリント h=(A/180*3.14159*M) // 南北の緯度差（ラジアン）の子午線弧長 i=(B/180*3.14159*N*n) // 起点緯度の東西の経度差（ラジアン）弧長 j=(B/180*3.14159*L*l) // 終点緯度の東西の経度差（ラジアン）弧長 T=(h**2+i*j)**0.5 // 起点（震源）から各終点までの距離　km =pr*(h,i,j,P,T) // 数値プリント ))R)) // 繰り返しルーティン R の終了 t=(T) // ルーティン Rで求めた「震源から各市町村までの距離」変量T を（東日本大震災）の距離変量 tとする ------------------------------------------------------------ ---- 起点（震源）から各終点までの距離 明治三陸津波 ---- U=(39.4) // ★明治三陸津波の震源の緯度のスカラー V=(144) // ★明治三陸津波の震源の経度のスカラー ((R)) // 繰り返しルーティン R の実行 .......　ルーティン Rで求めた「震源から各市町村までの距離」変量T がそのまま（明治三陸津波）の距離変量T ------------------------------------------------------------ ---- 震源からの距離順に並び替えと累和 東日本大震災 ---- j=r.l(t)blank // tの小さい順(定数項blankで欠測値にも末尾の順位)の順位変量j ....j=r.g(t)blank // 大きい順の場合は先頭の....を取る c=pmt(c,j) // 並び替え（順序数変量jによる） 浸水域人口 d=pmt(d,j) // 並び替え（順序数変量jによる） 人的被害 r=pmt(r,j) // 並び替え（順序数変量jによる） 人的被害浸水域人口比 t=pmt(t,j) // 並び替え（順序数変量jによる） 震源から市町村庁までの距離 p=pmt(p,j) // 並び替え（順序数変量jによる） 個体識別（印字）文字列 p,nam,:ci,p=pmt(p,j) // 印字系列pが文字系列であることを示す変量名に変更 ..................... ?c=cum(c) // 浸水域人口の累和?c ?d=cum(d) // 人的被害の累和?d ?r=(?d/?c*100)) // 累和の人的被害浸水域人口比?r =pr*(?c,?d,?r,t,p) // 数値プリント ------------------------------------------------------------ ---- 震源からの距離順に並び替えと累和 明治三陸津波 ---- J=r.l(T)blank // Tの小さい順(定数項blankで欠測値にも末尾の順位)の順位変量j ....j=r.g(T)blank // 大きい順の場合は先頭の....を取る C=pmt(C,J) // 並び替え（順序数変量Jによる） 人口 D=pmt(D,J) // 並び替え（順序数変量Jによる） 人的被害 R=pmt(R,J) // 並び替え（順序数変量Jによる） 人的被害人口比 T=pmt(T,J) // 並び替え（順序数変量Jによる） 震源から市町村庁までの距離 P=pmt(P,J) // 並び替え（順序数変量Jによる） 個体識別（印字）文字列 P,nam,:ci,P=pmt(P,J) // 印字系列Pが文字系列であることを示す変量名に変更 ..................... ?C=cum(C) // 人口の累和?C ?D=cum(D) // 人的被害の累和?D ?R=(?D/?C*100)) // 累和の人的被害人口比?R =pr*(?C,?D,?R,T,P) // 数値プリント =============== $$g // グラフセクション $g // スケールの目盛り指示コマンド（標準10ポイント） ?c,002 // 変量?cの目盛りを細かく2ポイントごとに ?d,002 ?r,002 $z // ゼロ軸表示 ?r?R // 変量?r,?Rのゼロ軸表示 $3 // ３次元図　 縦軸に【（浸水域)人口の累和】 横軸に【震源からの距離】 ?c,t, ,p=c,* // 縦軸?c,横軸t,奥行軸なし,個体識別p=バブル変量c,合成用保存* 東日本大震災 ?C,T, ,P=C,* // 縦軸?C,横軸T,奥行軸なし,個体識別P=バブル変量C,合成用保存* 明治三陸津波 // 合成 （浸水域)人口の累和と【震源からの距離】のグラフ　　　　　　　　 $3 // ３次元図　 縦軸に【人的被害の累和】 横軸に【震源からの距離】 ?d,t, ,p=d,* // 縦軸?d,横軸t,奥行軸なし,個体識別p=バブル変量d,合成用保存* 東日本大震災　 ?D,T, ,P=D,* // 縦軸?D,横軸T,奥行軸なし,個体識別P=バブル変量D,合成用保存* 明治三陸津波　 // 合成 人的被害の累和と【震源からの距離】のグラフ $3 // ３次元図　 縦軸に【累和の人的被害（浸水域）人口比】 横軸に【震源からの距離】 ?r,t, ,p=d,* // 縦軸?r,横軸t,奥行軸なし,個体識別p=バブル変量d,合成用保存* 東日本大震災 ?R,T, ,P=D,* // 縦軸?R,横軸T,奥行軸なし,個体識別P=バブル変量D,合成用保存* 明治三陸津波 // 合成 累和の人的被害（浸水域）人口比と【震源からの距離】のグラフ ======================= $$ // 終了セクション