2016年8月7日

高鐵上的氣壓變化-新竹至台中 之二

昨天做了新竹到台中的高鐵氣壓變化之後,從苗栗站到台中站之間「沒隧道」的路段,有個變化讓我很疑惑。

出了苑裡隧道沒多久,氣壓開始降低,行車兩分後氣壓降到最低再回升進入神岡隧道。

苑裡到神岡

這到底是發生什麼事情了。於是我從「政府資料開放平臺」的高鐵資料和中研院的DEM,整合做了張高鐵路程海拔圖來看看,終於找出答案了,而且是我意料之外的事情。


如果高鐵沒隧道-苑裡神岡大甲溪

如果高鐵沒有隧道或是地下化工程,直接在地面上移動,高度變化是怎樣呢?忽高忽低的,所以高鐵在北部到中部才會有許多的隧道和地下工程,過了彰化之後,幾乎就是一路平坦。不過並不是每一個高低起伏變化大的地方都需要挖隧道,從苑裡到大甲溪這段就沒有(也就是苑裡隧道到神岡隧道)

在南下高鐵列車上,出了苑裡隧道之後,苑裡的海拔大約是70m左右,而后里、外埔的海拔約200m,然後就過大甲溪,準備進神岡隧道了。

高鐵在這一段爬升了一百多公尺,只花了兩分鐘,我有印象這段我也是耳膜不舒服。

看整個海拔資料,我想后里外埔這段應該就是高鐵在地面上的最高海拔了。而且這個點附近有個地名相當有意思,叫做「打死人崎」,這個地名還有上新聞喔。











==================================================
接下來紀錄一下作法,屬於技術性文件,給自己看的

高鐵的里程在政府資料開放平台這裡
http://data.gov.tw/node/8011 (高鐵里程座標)

下載之後是個CSV檔,紀錄了里程和經緯度,我需要把每個點的海拔也找出來。所以我需要數位高程的資料。

有了兩個資料之後,用QGIS開啟
里程資料用  圖層/加入圖層/加入分隔文字圖層  
DEM資料用  圖層/加入圖層/加入影像圖層

2016-08-07 10:29:20 的螢幕擷圖

然後QGIS要安裝一個plugin,叫做Point sampling tool
附加元件/管理與安裝附加元件
輸入Point sampling tool來搜尋,點選安裝就可以

安裝好之後,這工具會在附加元件的Analyses
點選之後在Layers containing sampling points,選擇THRC_MILE(這是里程資料)
Layers with fields/bands to get values from選擇twdtm_asterV2_30m(這是數位高程)
底下的output point vector layer是選擇要存在哪個檔案裡
按下確定之後,它就可以把每個里程點都找到海拔然後存入那個檔案
我是存成pointAltitude.csv的檔案

接下來就用gnuplot來繪圖就可以
#!/usr/local/bin/gnuplot -rv
set datafile separator ","
set xlabel "distance(100m)"
set ylabel "height(m)"
set ytics nomirror
plot  "pointAltitude.csv" using 0:3 with lines lc rgb "red" notitle

set arrow from 14,graph(0,0) to 14,graph(1,1) nohead
set arrow from 21,graph(0,0) to 21,graph(1,1) nohead
set arrow from 326,graph(0,0) to 326,graph(1,1) nohead
set arrow from 649,graph(0,0) to 649,graph(1,1) nohead
set arrow from 999,graph(0,0) to 999,graph(1,1) nohead
set arrow from 1652,graph(0,0) to 1652,graph(1,1) nohead
set arrow from 1947,graph(0,0) to 1947,graph(1,1) nohead
set arrow from 2210,graph(0,0) to 2210,graph(1,1) nohead
set arrow from 2557,graph(0,0) to 2557,graph(1,1) nohead
set arrow from 3212,graph(0,0) to 3212,graph(1,1) nohead
set arrow from 3563,graph(0,0) to 3563,graph(1,1) nohead

pause -1


像各個站的標示其實是先去找站區的里程
http://data.gov.tw/node/8010 (高鐵車站資料)
然後去對照那個車站里程在高鐵里程座標裡的第幾行,比方說台北站的里程應該是第14行,繪圖的時候就用這樣的指令畫出來
set arrow from 14,graph(0,0) to 14,graph(1,1) nohead