花蓮地震 20080812

今天花蓮連續發生了幾次地震 第一次大概四級 接下來的都比較小 地震發生前都有一些低頻的聲音 不是很確定是不是建築物搖晃的聲音 可能不是 因為在感覺道建物搖晃之前好像就聽到聲音。 第一次聽到這種聲音，但是似乎1022嘉義地震據說有這樣的現象，這樣的聲音也被許多文獻紀錄著！

Install GMT (Generic Mapping Tools) under Windows XP (and Windows VISTA)

1. Go to SOEST (University of Hawaii), click on "Mirrors" in the left, and pick an appropriate FTP server;
2. Download netcdf-3.6.1-win32.zip, unzip the file and put it in a folder and rename the folder to "NETCDF"
3. Click on "4", to enter GMT version 4's directory, and download the files listed below:

GMT_exe.zip All executables
GMT_share.zip Support run-time files (patterns, cpt)
GSHHS_coast.zip Basic GSHHS coastlines

GMT_suppl_exe.zip Supplemental executables
GMT_pdf.zip PDF documentation and Unix man pages
GMT_man.zip Unix man pages for GMT
GMT_web.zip HTML versions of all documentation
GMT_tut.zip Tutorial data sets
GMT_src.zip All source code
GMT_scripts.zip Example scripts and data
GMT_suppl.zip Supplemental program source and docs
GSHHS_high.zip High resolution GSHHS coastline files
GSHHS_full.zip Full resolution GSHHS coastline files

Unzip all the files whose name started with "GMT", this will create ONE folder named "GMT", in this folder there will be many sub-folders:

bin (gmtprograms)
examples (gmtexamples)
lib
man
netcdf
share ( rivers, borders, shorelines, patterns etc)
src
tutorial
www

4. move the "GMT" folder to C:\
5. move the "NETCDF" folder to C:\ (this folder was from step 2)
6. Copy C:\NETCDF\netcdf.dll to C:\windows\system32
7. In control panel -> system -> advance -> environmental variables, find "system variables", and add one item:
variable name is "GMTROOT", variable is "C:\GMT"; (all without the double quotes)

In PATH, add C:\GMT\bin

and save

在 Windows XP 下裝設 GMT (Generic Mapping Tools)

1. 到
夏威夷大學SOEST的網站 點至左方Mirrors， 選取適當的FTP站 （在台灣大概都是日本的清水 Shimizu）

2. 下載netcdf-3.6.1-win32.zip， 解壓縮後，置於一新增檔案夾中，更改檔案夾名稱為NETCDF

3. 進入"4"， 即GMT 4版的目錄，下載以下的檔案

GMT_exe.zip All executables
GMT_share.zip Support run-time files (patterns, cpt)
GSHHS_coast.zip Basic GSHHS coastlines

GMT_suppl_exe.zip Supplemental executables
GMT_pdf.zip PDF documentation and Unix man pages
GMT_man.zip Unix man pages for GMT
GMT_web.zip HTML versions of all documentation
GMT_tut.zip Tutorial data sets
GMT_src.zip All source code
GMT_scripts.zip Example scripts and data
GMT_suppl.zip Supplemental program source and docs
GSHHS_high.zip High resolution GSHHS coastline files
GSHHS_full.zip Full resolution GSHHS coastline files

將GMT開頭的一起解壓縮（這樣會出現一個檔案夾叫做GMT，其下有

bin (gmtprograms)
examples (gmtexamples)
lib
man
netcdf
share ( rivers, borders, shorelines, patterns etc)
src
tutorial
www
等的檔案夾 裡面又各有許多檔案。

4. 將這個GMT檔案夾移至 C:\

5. 將netcdf-3.6.1-win32.zip 解壓縮置於 C:\NETCDF

6. Copy C:\netcdf\netcdf.dll 到 C:\windows\system32

7. 到控制台->系統設定->進階->環境變數 的 “系統變數“中

增加一項

變數名稱為 GMTROOT  而變數值為 C:\GMT 

PATH 中 加入 C:\GMT\bin

儲存

8. 在DOS窗框中 打入 pscoast 看看是否有錯誤訊息; 若是有一長串說明檔出現 就表示gmt本身安裝成功

9. 下載Ghostscript以及GSView 並安裝於系統中

ROI-PAC 產生檔案

Description of SAR images

Image Types
cpx = alternating real and imaginary pixels
rmg = alternating lines of amplitude and phase
• slc: cpx
• int: cpx
• cor: rmg
• unw: rmg
• hgt: rmg
• amp: cpx

一些有用的 Unix 連結

General Unix Links:

Unix tutorial for beginners

Unix for beginners by Gerard J. Kleywegt (written by an X-ray crystallographer). Comprehensive.

Unix for beginners by Ian Tickle (written by an X-ray crystallographer). Some Dos/unix conversions.

Unixhelp for Users Edinburgh Univ. (a bit more advanced)

A note about Unix Shells William Scott's page. (brief notes about [t]csh vs. sh-type shells).

Apple's Command Line Primer covers basic shell concepts, frequently used commands, environment variables, and running programs.

使用StaMPS以及Doris

1. copy 所有SLC到 ~/StaMPS/SLC/
2. 將檔案夾名稱改成日期 (yyyymmdd)
3. 選取master影像
4. 進入master檔案夾
5. 執行 “step_master_read_whole"
6. 決定crop的區域
7. 在master檔案夾 ccp $MY_SCR/master_crop.in .
8. 修改其中的起始與結束位置
9. 回到SLC檔案夾
\ls -d [1,2]* > make_slcs.list
去掉其中的master日期
10. 執行 make_read

暫時結束

使用awk來處理USGS地震資料

USGS 的地震資料格式如下：

Year Month Day Time(hhmmss.mm)UTC Latitude Longitude Magnitude Depth
1973 01 30 210112.50 18.48 -103.00 7.5 43
1973 02 06 103710.10 31.40 100.58 7.7 33
.
.

如果想要利用GMT畫地震分布時，awk會是一個很方便的工具（相對於Micro$oft的東西來說， 利用很簡單的shell script就可以轉換檔案：

awk '{if (NR>1) print $6, $5, $7, sqrt($8)}' infile > outfile

其中要輸出的 field 也可以做運算 例如上面的'sqrt($8)'