今天搭高鐵前往台中，我又做了一次高鐵上的觀測氣壓
原來這整段有31段隧道，顯見竹苗一帶的丘陵地形，這些隧道多到我沒辦法在圖上標示了。
編號11的新城隧道一直到編號41的神岡隧道
1.觀察氣壓變化，我認為新苗隧道、頭份一號隧道和頭份二號隧道至少都有一個豎井。氣壓在隧道內都會有一次大的變化。
2.途中氣壓變化最大的不在隧道段，是在兩個隧道之間（苑裡隧道-神岡隧道）
3.上回推測某些隧道表現出來的氣壓變化可能是跟明挖覆蓋有關，但是觀察神岡隧道（最右邊那個）卻沒看到這個現象。
今天回程的時候也進行了一次觀測，但是時間是六點多，無法藉由外界天光與隧道的光線差異去偵測隧道所在，於是想了應該也可以透過偵測聲音改變。
新竹到苗栗段的聲音是怎麼變化的呢？偵測平均音量，結果是這樣，可以觀測到隧道內的確有比較吵。我還想作的是去作FFT，把特定一個頻率段抓出來，其實在隧道裡會聽到一個聲音，幾乎是一聽就知道進隧道了。中間的低谷（8:23左右）是苗栗站，過了苗栗站，整個高鐵就吵起來了，不過不是乘客的吵，是列車的整體聲音。
回程（台中-新竹）的高鐵倒是一致性的都很吵，不過這班沒有停苗栗
作圖的時候，為了要把隧道段自動highlight，特別去找了資料來參考。上回因為隧道少，自己手動用inkscape畫一畫，這次多達31個，根本沒辦法手動。
參考
Plotting the recession - some basic data analysis with gnuplot
基本就是設定這個函數
g(x) = (x<1000?GPVAL_Y_MAX:GPVAL_Y_MIN)
然後先跑一次圖
plot "light" using 2:4
獲得GPVAL_Y_MAX和GPVAL_Y_MIN的數值
然後用下面這指令畫圖的時候，就會自動把小於1000的亮度值Highlight，用個灰色方塊畫出來（其實是個被填滿的曲線）
plot "light" using 2:(g($4)) with filledcurve lc rgb "#eeeeee" notitle,\gnuplot的命令稿
-----------------------------------
set xdata time
set datafile separator ","
set timefmt "%H:%M:%S"
set format x "%H:%M"
set xrange ["08:12:00":"08:42:00"]
set xlabel "Time"
set ylabel "light"
set y2label "pressure"
set ytics nomirror
set y2tics
g(x) = (x<1000?GPVAL_Y_MAX:GPVAL_Y_MIN)
plot "light" using 2:4
plot "light" using 2:(g($4)) with filledcurve lc rgb "#eeeeee" notitle,\
"light" using 2:4 with lines lc rgbcolor "#f4c2c2" title "light" axis x1y1,\
"pressure" using 2:4 with lines title "pressure" lc rgbcolor "blue" axis x1y2,\
pause -1