NETSTARS輸砂模式資料轉換工具

NETSTARS為謝慧民老師與李鴻源老師共同開發之擬似二維的沖淤模式，可模擬支流、複雜河系、陡坡、緩坡、水躍、定量流及變量流之水理及相對應底床沖淤特性，更深了解NETSTARS可詳見謝慧民老師的網站 http://www.hmhsieh.idv.tw/hmhsieh/netstars.htm。
但是 NETSTARS在資料輸入上，仍需藉由資料卡輸入的方式建置，最近剛好需要跑NETSTARS模式，就藉此機會利用 Matlab與Excel建立”NETSTARS輸砂模式資料轉換工具”，目前為1.00版，此工具特色包括:

• 使用者能利用Excel介面輸入資料，操作簡單
• 不須額外安裝Matlab (但需安裝MCRinstaller.exe)
• 能夠一次轉換全部的斷面資料(GR卡和NS卡)
• 能夠一次轉換全部參數(BR與BS卡)
• 自動在最後一筆資粒加上負號(資料結束)
• CQ卡的功能
• 目前僅現單一河段

下載位置:  NETSTARS輸砂模式資料轉換工具 V1.40

Google Earth 6.0版推出了

這次新版的Google Earth整合Google Map的街景服務，在導覽工具上出現街景人形圖示，操作方法跟Google Map一樣把圖案拖進即可。

Google Earth 6.0 新功能

與Google Map不同的是，在Google Earth瀏覽街景更為順暢、也更直覺。

新版本個另一項賣點就是3D樹景，但是目前似乎僅於在東京(Tokyo)等大都市郊區有。

3D樹景

當然新版本還加強其他原有的功能，詳情見 Google Earth http://www.google.com/intl/zh-TW/earth/index.html
Google Earth 6.0 下載位置  http://www.google.com/intl/zh-TW/earth/download/ge/agree.html

中央氣象局大雨和豪雨的定義

目前中央氣象局大雨和豪雨的定義為:
大雨： 指24小時累積雨量達50毫米以上，且其中至少有1小時雨量達15毫米以上之降雨現象。
豪雨： 指24小時累積雨量達130毫米以上之降雨現象。
如24小時累積雨量達200毫米以上之降雨現象，則為大豪雨。
如24小時累積雨量達350毫米以上之降雨現象，則為超大豪雨。

1.豪雨、大豪雨、超大豪雨之間都只差一個字，很難分辨，建議比照颱風強度的描述方法，較為一班大眾接受，另外，最大等級僅24小時350mm完全不敷使用。
我建議的雨量分級標準:
急雨:1小時雨量超過50mm以上之降雨現象 (超過一般都市排水極限)
輕度豪雨:累積雨量超過200mm。
中度豪雨:累積雨量超過500mm。(可能會在地質脆弱區引起土石流)
強烈豪雨:累積雨量超過800mm。(達到一般地區山崩土石流的臨界值)
超強豪雨:累積雨量超過1200mm。(將引起嚴重的淹水或山崩等天然災害)

利用ArcGIS Model Builder建立差異化植生指數(NDVI)模組

差異化植生指數(Normalized Difference Vegetation Index；NDVI)是一個常見的環境因子，利用多光譜之衛星影像用來評估地表上的植生覆蓋情形，用來計算植物生長與生產力，另外也常被用來應用在崩塌地判釋或土壤沖蝕量估算或是土地利用變化之監測，原理是利用植生在不同波段之反射率的差異，將此差異標準化，如圖A所示，植生(綠線)在近紅外光段反射率約50%；而在紅光段的反射率卻只有10%左右。

圖 A 土壤、水和植生的典型光譜反射曲線
利用上述之特性，NDVI方程式表為:  NDVI = (NIR – RED) / (NIR + RED) ，其中NIR為光譜近紅外光段；RED為光譜紅光段，計算的值範圍會在 -1 至 1 之間。

賈伯斯在史丹佛的演講

Steve Jobs' 2005 Stanford Commencement Address

7:30秒之後才是演講內容  

走過大崩壞 為台灣土地把脈

走過大崩壞 為台灣土地把脈01

我們的島 宜蘭做大水(2010/10/26)

我們的島  宜蘭做大水  片頭


線上影音
http://hichannel.hinet.net/event/pts/player/ptsplayer.html?eid=74&id=247601

統計面資料內polygon的網格

我們在使用 ArcGIS時，常用到需要向量式的shapefile和網格式資料要交互分析，例如:一個集水區的平均坡度或平均高程，由於集水區通常為shapefile，坡度圖通常圍網格資料，因此一個集水區通常會包含許多網格資料，要如何去計算最大值、最小值、加總或平均值，這時就需要用到zonal stastics，以下就計算坡單元平均坡度為例，來說明如何使用zonal stastics


 
1. 開DEM算slope

2.把坡單元也加進來

 
3.check坡單元欄位中，必須有一個欄位是可以分別每個feature（FID不算），如果沒有，新增一個long integer的欄位，再用field calculator，欄位等於FID即可。

 
4.zonal

 
按OK


 
5.EXPORT

完成


 
開啟 Table欄位中的mean則代表，平均坡度。


 
坡向的算法也是類似，但是坡向比較麻煩，因為1度和359應該相似，但是計算出的結果=(1+359)/2 = 180，導致北向坡與坡向變化大的坡有的錯誤，如果懶的時候也是可以這樣算。正確做法要先算出坡向圖後，把角度換成向量ex:(0,-1)，去做zonal，再用sin與cos去分辨像現，則可以計算出平均坡向180。


 

Matlab讀取 txt文字檔, csv, excel

path = 'C:\\';  %設定絕對路徑

%read txt
cwbdata = load(strcat(path, 'cwb_466920_2006.txt'));

%read CSV
cwbdata = csvread(strcat(path, 'cwb_466920_2006.csv'),2,0);

% Matlab 7.0以後才有的功能
%read excel
cwbdata = xlsread(strcat(path, 'cwb_466920_2006.xls'));

[VB6] 讀區資料夾中所有檔案名稱

Private Sub CommandButton1_Click()
ComboBox1.Clear
Dim s As String
  s = Dir("C:\*.*")
  Do Until s = ""
    ComboBox1.AddItem (s)
    s = Dir()
  Loop
End Sub

Matlab常用的算術函數

Matlab常用的算術函數(Arithmetic functions)


sin(x), cos(x), tan(x)


exp(x), log(x), log10(x)

- round(x) 四捨五入(Rounds a number to the nearest integer)


- ceil(x) 無條件進位(Rounds a number up to the nearest integer)


- floor(x) 無條件捨去(Rounds a number down to the nearest integer)


- fix(x) 往0的方向取最接近的整數，ex: fix(-3.3) = 3.3 (Rounds a number to the nearest integer towards zero)


- rem(x,y) x/y的餘數(The remainder left after division)


- mod(x,y) 也是餘數，但是和rem()的差別在於處理負數有差(The signed remainder left after division)


- abs(x) 絕對值(The absolute value of x)


- sign(x) 判釋正負The sign of x

iRobot 打掃機器人

最近新買一台iRobot 打掃機器人，操作介面非常簡單，只要按下Clean鍵他就會開始打掃，本來擔心他沒辦法房間的角落，經過約50分鐘之後，他已經跑遍我們家的每個地方，掃完之後比我自己跪在地板上擦地還要乾淨。

打掃完後會自己回去休息。

整個家裡都乾乾淨淨，心情也就好了。

Optimally Interacting Minds 臭皮匠勝諸葛亮 但得實力相當

過去我們常聽到三個臭皮匠勝過一個諸葛亮，也就是說集合眾人的智慧可以得到更好的想法，而Science這一期就刊出由University College London團的的研究指出，廣為諮詢眾人的智慧確實可以做出更好的決定，但是必須是眾人必須有相當的智慧，否則會得到更為糟糕的決定。

老照片-玉山主峰

2003年8月13日，山社最後一個大眾活動隊 玉山群峰 。

玉山登山口之早晨

往北峰路上回望主峰

8/11 台大→塔塔加→C0

8/12 C0→登山口→西峰觀景台→排雲山莊→圓峰山屋→C1

8/13 C1→玉山主峰→東峰→北峰→圓峰山屋→C2

8/14 C2→南峰→東小南山→圓峰山屋→C3

8/15 C3→西峰→前峰→登山口→台北

舊照片-金瓜石草山

2004年10月30日，那時摔車摔傷的腿都還沒好，就忍不住想去爬山，一拐一拐地從東北角海邊爬上500m的草山戰備道。
迎面著東北季風的金瓜石，總是第一個迎接到秋天濕冷的風，水氣化作霧氣徘徊於山峰間，與根根芒草結為露水。

如何在ArcGIS中設定TWD97與TWD67座標系統(2)

繼上一篇  如何在ArcMAP設定TWD97與TWD67座標 ，因為這個方法實在有點麻煩本身在設定座標，而且實在沒把握到底參數設定的正確性，建議可以採用別的高手已經用設好的座標檔。推薦中研院計算中心的一個網站，裡面有許多很好用的GIS小工具。

1. 開啟 研究支援工具(http://www.ascc.sinica.edu.tw/gis/ISTIS/tools.html)

[caption id="attachment_142" align="alignnone" width="503" caption="研究支援工具"][/caption]

1. 開啟壓縮檔，將檔案解壓縮到 C:\Program File\ArcGIS\Coordinate Systems

[caption id="attachment_144" align="alignnone" width="608" caption="解壓縮至Coordinate Systems"][/caption]

1. 開啟 ArcMAP >> Layers  >> Properties >>  Coordinate Systems >> Predefined就可以看到剛剛課訂好的坐標系系統。
   
   [caption id="attachment_145" align="alignnone" width="638" caption="開啟Properties "][/caption]
   
   [caption id="attachment_148" align="alignnone" width="486" caption="看到剛剛加入的座標檔"][/caption]

如何在ArcMAP設定TWD97與TWD67座標 1

TWD97和TWD67是台灣常用座標系統，在ArcGIS上面如何設定？

各大瀏覽器佔有率

Source: StatCounter Global Stats - Browser Market Share

資料來源:http://gs.statcounter.com/

Chrome在八月初時已經突破10%，Firefox則持平，IE為52%，已經快要跌破一半

北二高往基隆段發生三十年來第一場嚴重山崩事件

在較大的地圖上查看山崩事件

Google 街景 崩塌前 http://0rz.tw/mEqma



圖片來源：路透社



圖片來源：路透社

這突如其來的崩山真是讓人嚇了一跳，在沒有豪雨或地震的情況下，會產生如此重大的崩山，實在令人有點匪夷所思。

根據中央氣象局雨量觀測資料

縣市名稱	雨量(毫米)	測站名稱	測站代碼	所在地
基隆市	1.0	五堵	01B03	基隆市七堵區(十河局)
台北縣	5.0	太平	L1A84	台北縣雙溪鄉(翡翠水庫)
台北市	7.5	格致國中	A1AA1	台北市士林區(水利處)

基隆市過去2天以來的最大降雨量僅僅1.0mm，等於沒有下雨。

但是從Google Earth 2006年3月18日 影像上，似乎可以看到一條V字形的溝，不像是一般的沖蝕溝，是不是張裂縫還很難說，而且崩塌發生的範圍是大於VV字形的溝。



資歷來源：Google Earth

從經建版地形地形圖來看，上邊坡坡度大約是10～15度，從直覺上來判斷，既不像是個舊崩塌，更不覺它是個有危險的崩塌地坡面。



紅色星號與虛線為崩塌區

但是若是從地質圖來看，可以發現這邊是位於大寮層的順向坡，傾角約10度。



地質圖，黃色線為北二高，紅色路段為隧道

(註：地圖會說話 提到與林肯大郡在同一岩層上)

大寮層地層特性為砂頁岩互層，當雨水滲入頁岩層上方，就會大幅降低岩層間的內摩擦角，可以說是最適合順向坡破壞的地質條件。

四個條件構成大規模崩塌的要素：

1. 10度的地形坡度


2. 10度的地層頃角


3. 坡角開挖造成的自由面


4. 砂頁岩互層


5. 雨水滲入

以上四個要素都是先天地質不良的要素，但是缺乏雨水滲入的trigger factor而產生崩塌，還是令人感到相當意外，只能說這種地質條件發生崩塌災害機率是很低的，但是一旦發生，其災害的規模通常相當大，可以說是最難以防範的崩塌災害。

如何讓GridView刪除資料時，出現確認對話框

如何讓GridView刪除資料時，出現確認對話框：

網路上有一推範例，找到一種最簡單的方式。

<asp:LinkButton ID="EditButton" runat="server" CausesValidation="False"

CommandName="Edit" Text="編輯" />

&nbsp;<asp:LinkButton ID="DeleteButton" runat="server" CausesValidation="False"

CommandName="Delete" Text="刪除" OnClientClick="if(!window.confirm('確定要刪除此項目嗎？'))return false;" />

&nbsp;<asp:LinkButton ID="NewButton" runat="server" CausesValidation="False"

CommandName="New" Text="新增" />