2016年7月17日

地震模型(2)-用QGIS作地震模型

還好我有寫紀錄的習慣啊,為了這個模型,我可是自己看了自己的影片呢

要作的東西其實就跟這差不多《QGIS作透明的立體等高線地形模型》,只是把等高線換成震央分佈罷了。

所以我就只要把前一篇做出來的earthquake.csv匯入QGIS當一個圖層就好。

匯入圖層的方式很簡單
圖層(L)/加入圖層/加入分隔文字圖層..
選擇好檔案之後,記得把X座標設成第三欄,Y座標設成第四欄

接著再按照這篇《QGIS製作行政區圖》,到政府開放平臺,找直轄市、縣市界線(TWD97經緯度)的SHP檔案,拿去匯入QGIS就可以了。

有縣市界線,又有地震資料,就可以做出像這樣的圖


TaiwanEarthquake03


可是還沒完啊,不是說要作立體模型嗎?那就依照前篇的觀察資料結果,每10公里的震源放一張圖,反正台灣地震大多數都在100公里以內的深度。

至於切的方法,就照著這篇的作法來弄就可以了《QGIS作透明的立體等高線地形模型》,要記得深度是第5個欄位就是了。

earthquake


製作好的PDF檔就在此,有分兩個,是把震源標示的點用不同大小來作的,就看看哪個比較清楚了。

https://drive.google.com/file/d/0Bzwhi7Oh9a5ZMFRNREh3TjloZzg/view?usp=sharing

https://drive.google.com/file/d/0Bzwhi7Oh9a5ZVHFCR0FXQ0FVdEE/view?usp=sharing


https://drive.google.com/file/d/0Bzwhi7Oh9a5ZVzZWckl0cmlBLTA/view?usp=sharing
呈現規模4以上的地震,並且把震央的點圖示調小,台灣的縣市邊線加粗

地震模型(1)-地震數據分析


前幾天看到yenwen貼了一個連結,是台灣地震活動分布立體模型

覺得很有意思,想到去年去車籠埔斷層保存園區看到的這個模型:
「1973-2012年台灣規模四以上地震震源及活動斷層分布」

我猜是壓克力作的,然後用雷射雕刻在不同深度,這樣就可以藉由震源看到板塊交界的位置呢!既然興趣來了,我就想來自作一下。








但是數據打哪來呢?在中央氣象局的「地震活動彙整」中,可以下載到歷年的資料

http://www.cwb.gov.tw/V7/earthquake/rtd_eq.htm

本來我想一口氣把所有資料都下載,不過那樣作,反而伺服器不給下載,後來我只好500筆500筆的抓,抓了19個csv檔案,再把這些檔案用cat 和資料流導向的>符號,把所有的檔案結合在一起
cat *.csv > all.csv

統計了一下從1995年1月到抓到資料的當下,合計9433個地震,顯著有感地震3063個, 小區域地震6370個。

然後,就是整理資料的部份了,首先是每個csv檔都有標題,所以我就直接用複製取代的方式把所有的標題都刪掉了。
編號,地震時間,經度,緯度,深度,規模,震央位置

再來是時間的部份,把預設的12時制轉換成24時制,這樣才能送給gnuplot去處理。
12時制實在有點難搞
上午12點是凌晨12點,也就是24時制的0時
下午12點就是中午12點,也就是24時制的12時

這部份一樣用複製取代的方式去作,得到的檔案就在這
earthquake.csv


先用gnuplot看看地震規模的直方圖和地震深度的直方圖吧,gnuplot用的命令稿我就寫在圖片底下。


這是地震深度和次數的比較圖,台灣的地震深度少數可以到300公里之下,多數都是在50公里以內的。
地震深度次數


set datafile separator ","
plot "earthquake.csv" using 5:(1) smooth frequency with boxes




這個是特別把10公里之內的深度放大來看

地震深度次數10公里




地震規模大約是3到5之間是最多的
地震規模次數

set datafile separator ","
plot "earthquake.csv" using 6:(1) smooth frequency with boxes




這是把深度和規模放在一起看,深度深的,地震規模都蠻大的。淺源地震的規模有大有小

地震規模深度

 plot "earthquake.csv" using 5:6 with dots



既然有了一些資料,就想來看看很多年前有學生問的問題,他問地震是不是都在晚上。我認為不會和時間有關係,不過那時候沒證據,現在就拿已有的數據來分析。圖在此。
地震和時間

set datafile separator ","
gnuplot> set timefmt "%Y/%m/%d %H:%M:%S"
plot "earthquake.csv" using (tm_hour(timecolumn(2))):(1) smooth frequency with boxes fill solid 0.25




「那麼哪一天地震多呢?」既然資料都有了,就分析一下囉。如預期一般,31號果然很少(廢話),但有幾天的次數看得出來低些(但不一定有顯著差異)
地震和day

set datafile separator ","
gnuplot> set timefmt "%Y/%m/%d %H:%M:%S"
gnuplot> plot "all-cuttitle-adjustdate.csv" using (tm_mday(timecolumn(2))):(1) smooth frequency with boxes fill solid 0.25



哪一年地震多呢?這答案不一定能從數據中得到,也許測站變多,或是精密度提高,都會讓地震的數據增加。不過單就圖上來看,1999年之後,數據是比較多了
地震和year

set datafile separator ","
gnuplot> set timefmt "%Y/%m/%d %H:%M:%S"
gnuplot> plot "all-cuttitle-adjustdate.csv" using (tm_year(timecolumn(2))):(1) smooth frequency with boxes fill solid 0.25


ㄟ...不是說好要作地震模型嗎?哎呦,總是得把數據收集好,先分析觀察一下囉。下一篇就來講怎麼製作成模型

2016年7月13日

修復看片箱

幾年前,一位獸醫朋友送我一個壞掉的看片箱,我高高興興拿回家之後,就放進櫃子裡,過了好幾年才想起「啊,有這個東西啊」

其實當初是有拆開來看哪裡有問題,不過當時年紀小,什麼都看不懂。阿現在懂得多一些了,看一看之後,交叉檢驗就發現是燈管出問題了。
DSC05649


這麼細的燈管叫做冷陰極螢光燈管(CCFL),以前真的不知道這個東西,是後來逛電子材料行才知道有這個的。這種管子常用在電腦螢幕裡頭,不過現在很多都是改用LED了。知道材料之後,要修復就簡單了,上網拍找一下,看到尺寸合適的買回來換上去,就成功了。
DSC05652


前面幾篇文章裡,拍照的背景亮亮的就釋放在這看片箱上拍照的,像是以下幾張照片。拍照的時候,EV值要再調高一點。


IMAG2173


DSC05654


DIY玻片放置架-第二版

六年前寫這篇「DIY玻片放置架」之後,其實沒有怎麼常使用,還是都讓學生把載玻片洗乾淨插在講桌前面的架子。

不過早上我在整理東西的時候,看到之前這些隔板,突然想到另一種的組合方式。試了一下,覺得非常好用,好拿好放,底部還能通風滴水。這種尺寸的隔板剛好可以拆成7格的兩片和12格的兩片,組合成這樣,就可以放9片玻片了。算一算手邊的隔板,每小組放個兩組放置架是足夠的呢。
DSC05654

間隔有長的和短的,我建議是插在短的那邊,這樣才好放置和取出。玻片有沒有缺少,有沒有洗乾淨,都可以一目了然。俯看可以算玻片,側看可以看有沒有洗乾淨,有洗乾淨,9片玻片就可以透過去。
DSC05655

還有另外一種放置方式,跟之前那種類似,也是用插的,不過比起上面那種,其實不是太穩固
DSC05657


拿掉玻片的樣子就是這樣
DSC05658

DSC05659

這個高是6公分左右
DSC05661

每塊方格,大約是1公分寬。不過也不用拘泥於尺寸,反正就買小型的隔板回來試一試,應該都能組合出這種玻片放置架。
DSC05662

2016年7月8日

自製印章研習

雖然海報上寫的是marker,但是這場研習真是的maker,作的是印章。
教務主任接了這個研習計畫,請我協助規劃,其中一場我就找了季瑾一起來作這個研習。

IMAG2125

因為之前她有一個課程就是讓學生畫圖片之後,請老王幫忙用印章機做出實體印章,我覺得太有趣了,而我自己作印章的課程則是這個《氮循環遊戲-氮原子的旅行


分享完我倆的課程之後,就是用老王的印章機幫大家作印章。大家在白紙上先畫出自己想要作印章的圖案,然後影印在描圖紙上,接著再把描圖紙墊在光敏印章墊上面(中間夾一張透明片,避免碳粉沾粘),送到印章機裡頭。

光敏墊是一種可吸墨的材質,但是在強光照射下,材質特性會改變,變成無法吸墨。而印章機是一個會發出強烈閃光的機器。在光敏墊上,有圖案的部份,因為有遮光,所以維持原有特性,可以吸墨,而沒圖案的地方就不能吸墨。

「啪」的一次大閃光,印章就完成了
  IMAG2123

這個就是這次研習中,諸位主任的印章成品就是這樣啦,這些可都是大家親筆畫出的喔,看起來還真像是用電腦印出來的吧!
  IMAG2122


你可能對這有興趣

Related Posts Plugin for WordPress, Blogger...