coop: 2011

2011-12-28

FEWS_Taiwan工作日誌2011.12.28

完成氣象局測站更新後觀測資料下載測試工作

並且同步進行報表模組使用測站id調整工作

修正測站id調整後報表模組無法取得資料之設定問題

9thRMO_Report 1.00 default.xml,Location with id C1Z020 does not exist;

5thRMO_Report 1.00 default.xml,Location with id C1M550 does not exist;

4thRMO_Report 1.00 default.xml,Location with id C0I090 does not exist;

3rdRMO_Report 1.00 default.xml,Location with id C1F870 does not exist;

2ndRMO_Report 1.00 default.xml,Location with id C0E530 does not exist;

1stRMO_Report 1.00 default.xml,Location with id C1U500 does not exist;

National Taipei University of Technology

Jhih-Cyuan Shen

2011-12-27

Maptoolkit

http://www.bikemap.net
使用Maptoolkit 進行google map 資料呈現的加值

Maptoolkit 網址
http://www.toursprung.com/products/maptoolkit/

2011-12-22

C1E460 --> C1E461
C1E450--> C1E451
C1E710-->C1E711
C1E690-->C1E691
C1E520 -->C0E520
C1E680 -->C1E681
C1E720 -->C1E721
C1E700 -->C1E701
C1E560 -->廢站
C1E550 -->C0E550
C1E570 --> 廢站本站與苗栗農改場(K2E360)設置地點相同
C1E540 -->C0E540
C1D430 -->C0D430
C1E440 -->C0E440
C1E500-->廢站
C1U680-->C0U680
C1A650-->C0A650
C1A640-->C1A640
C1G630-->C1G631
C1H970-->C1H971
C1G620-->C0G620
C1G660-->C0G660
C1G680-->廢站
C0G650-->C0G651
C0G640-->C0G641
C1M550-->C1M551
C1T900-->C0T900
C0T840--> C0T841
C1Z020-->C0Z020
C0K240-->C0K241
C1I130-->C1I131
C1A700-->廢站
C1A690-->廢站
C1A680-->C0A680
C1A670-->廢站
C1A660-->C0A660
C1A650-->C0A650
C1A640-->C0A640
C0A9I0-->C0A9I1
C0A9H0-->廢站
C0A590-->廢站
C0A510-->廢站
C1T900-->C0T900
C0T840-->C0T841
C1Z020-->C0Z020
C1M550-->C1M551
C1I200-->C1I201
C1I130-->C1I131
C1I120-->C1I121
C1I100-->C1I101
C0K240-->C0K241
C0I110-->C0I111
C0I090-->C0I091
C1H940-->C1H941
C1H9B0-->C1H9B1
C1I210-->C1I211
C1H960-->C0H960
C1G620-->C0G620
C1F9J0-->C1F9J1
C0F990-->C0F991
C1F9B0-->C1F9B1
C1F000-->C0F000
C1F970-->C0F970
C0H9C0-->C0H9C1
C1H970-->C1H971
C1F9A0-->C0F9A0
C1H950-->C0H950
C1F9E0-->C1F9E1
C1F9E0-->C1F9E1
C1F940-->C1F941
C1H9B0-->C1H9B1
C1F890-->C1F891
C1F920-->C0F920
C1F9D0-->C1F9D1
C1F9C0-->C1F9C1
C1F9H0-->C1F9H1
C0F860-->C0F861
C1F850-->C0F850
C1F910-->C1F911
C1F9G0-->C1F9G1
C0T790-->C0T791
C1F930-->C0F930
C0H9C0-->C0H9C1
C1E500-->廢站
C1F9F0-->C1F9F1
C1F870-->C1F871
C1C620-->C0C620
C0E530-->C0E531
C0E580-->廢站
C1E590-->C0E590
C1E600-->C1E601
C1D370-->C0D370
C1C500-->C0C500
C0C450-->廢站
C1C490-->C0C490
C1C540-->C0C540
C1C480-->C0C480
C1A710 -->C0A710
C1U620-->C0U620
C1U680-->C0U680
C1U500-->C1U501
C1U510-->C1U511
C1U620 -->C0U620

2011-12-21

2011濁水溪流域平均降雨量

一月 48.5mm二月 21.6mm三月 46.4mm四月 14.6mm五月 184.3mm六月 201.9mm1~6月累積降平均降雨量528.2mm

2011-11-19

REFOR模式與SOBEK模式計算成果比對

在FEWS_Taiwan平台下

利用十河局REFOR模式與SOBEK模式針對 1003豪雨進行計算測試

如圖所示為社后橋計算成果

粉紅色為REFOR模式

綠色為SOBEK模式

模式在計算成果各有其優缺點

目前模式在水文架構與使用資料上仍有些差異

後續可以針對細部設定近一步比對

使模式應用更加完整

National Taipei University of Technology
Jhih-Cyuan Shen

2011-11-17

QGIS Raster to Point

本文將說明在QGIS中如何將RASTER資料轉換為POINT資料 使用的套件包含RASTER （ GDAL Tools）Delimited Text Layer Plugin
在開始前務必確認相關套件已經啟用

Qgis raster to point

View more presentations from J.C Shen

2011-11-14

OGC WaterML2.0

OGC WaterML2.0
http://www.opengeospatial.org/projects/groups/waterml2.0swg
結合OGC GML與 Observations and Measurements (O&M) standard 的產物
主要由美國NOAA 澳洲 CSIRO 與荷蘭 Deltares 加拿大等單位合力推動的標準

Delft-FEWS 就上次使用這會議的方向也會支援此標準
National Taipei University of Technology
Jhih-Cyuan Shen

QGIS 使用Spatial Query 進行圖層空間分析

如何利用QGIS 進行圖形空間分析
在本文將介紹利用Spatial Query 進行空間分析之流程
首先安裝Spatial Query plugin
開啟兩個預計要進行分析之圖層
例如哪些點位在哪些多邊形之圖層之內
先選擇source 圖層為點位 point
接著下條件 where the feature "within"
在選擇Reference 圖層 polygon
再點選 apply進行空間分析
QGIS 會將選擇出來的圖層利用不同的顏色加以區分

但目前在分析流程上不會相商業軟體自動將選擇出來之成果自動分群組
（假設有同時多個多邊形同時分析時）
需要個別另存或分析資料
另外查詢出來之成果可以利用
Selected features 下的create new selection 將選擇的圖檔範圍另外存為一個新圖層方便使用者進行分析
但如果要後續可以持續使用可能要記得將圖層輸出成為shp檔
National Taipei University of Technology
Jhih-Cyuan Shen

2011-11-04

FEWS_Taiwan 工作日誌2011.11.04

1.水位站上游集水區預報降雨量無法正常計算
檢查原因為水位站上游集水區shp檔案更新（修正部份shp檔無法正常載入問題）
修正烏溪大安溪之水位站上游集水區shp
2.修正非等距資料展示時char圖會將最後一筆時序列資料會持續向後延伸之狀況
修該增加path 檔修正該問題
3.labels字型展示問題
原先字型需要大於13 顯示才會正常更新版本後已經可以將字型設定小於13
但建議字型展示設定大小不要小於12，否則使用者在進行辨識時會出現酷困難
4.新版本fews 可以將所有congfig檔利用壓縮檔加以存放在區域資料夾中方便後續複製使用
5.增加資料快速瀏覽功能，類似地圖之鷹眼功能。可以查詢長時間段時序列資料找到想要的區段進行顯示

2011 Delft-FEWS User meeting

2011-10-26

FEWS_Taiwan工作日誌2011.10.26

海堤構造物增加顯示資訊中包含海堤名稱、海堤最高高程、最低高程
（台灣本島、澎湖、馬祖、金門、綠島、蘭嶼等地區之海堤構造物資訊）
測站資料檢核城林橋 HP73 資料庫時間 2011.09.28 17:00 開始

National Taipei University of Technology
臺北科技大學土木工程系

2011-10-20

NTUT 100校慶 poster

2011-10-18

淡水河流域測站對應之斷面編號

淡水河
TE10 土地公鼻 HP02
TE13 獅子頭 HP03
TE24.A 臺北橋 HP04
TE36.A 新海大橋 HP06
新店溪
H01 華江橋
H10 中正橋 HP08
H02 萬板大橋
H17.A 秀朗 HP18
基隆河
KE11.A 百齡橋 HP41
KE043 南湖大橋 HP24
KE50 社后橋 HP23

2011-10-14

FEWS_Taiwan工作日誌2011.10.14

氣象局雨量站更新
<location id="C0U780" name="五結">
<location id="C1U830" name="烏石鼻">
<location id="C1U840" name="東澳嶺">
<location id="C1U850" name="觀音海岸">
<location id="C0AG90" name="中和">
<location id="C0AH00" name="汐止">

2011-10-07

FEWS_Taiwan 即時二維演算功能測試

在水利署利用FEWS_Taiwan銜接 sobek 2D模式進行即時二維淹水模式之演算
測試案例高屏溪支流美濃溪流域
目前使用80m 數值高程模型進行測試
每小時輸出一筆資料
19筆時序列資料計算時間約需3~5分鐘
如圖所示分別為水深與流速與流向之二維網格展示介面
後續需要持續針對流量與水深持續進行檢較工作
並且持續針對較高解析度之料進行測試

2011-10-03

在日本播出的犀利人妻節目

在icfm5 結束後回到飯店晚上看到的電視節目

畫質在數位電視的效果下

顆粒較細緻有立體的感覺

FEWS_Taiwan 河川局2011-10-03客製化自動水情報表輸出內容

如圖所示分別為第一河川局第八河川局與第九河川局
利用FEWS_Taiwan 自動產出之水情資訊
水位站上游集水區平均降雨資訊與河川局管轄之鄉鎮平均降雨資訊

荒川旁的超級堤防

在日本荒川旁的超級堤防

提供了在超級堤防上之住戶淹水問題改善之狀況

但在超級堤防旁的低窪地區可能因為

超級堤防之興建造成更嚴重之淹水現象

但要完成該項工程民眾政府與建築規劃單位需要有共識

且需要一段很長時間的溝通與協調工作

才能將超級堤防順利完工

日本首都圈外郭放水路

在日本江戶川河川事務所的管轄下的
首都圈外郭放水路
目前最大流量為平成20年8月（2008年）調節約1200萬立方米之水量

大約是是 25000個日本國小標準之游泳池

春日部市

2011.09.30 在icfm5 參訪行程中抵達了

春日部市

進行東京相關防洪設施與防洪應變單位之參訪

過程中看到春日部市利用蠟筆小新作為資訊查詢告示板主角的內容

2011-09-28

ICFM5 參考網站2

http://www.rijksoverheid.nl/documenten-en-publicaties/rapporten/2009/12/01/nationaal-waterplan-2009-2015%5B2%5D.html
http://www.ewi.tudelft.nl/live/binaries/a4751543-3c3e-4787-b3c2-ae34f717f351/doc/citatie149.pdf
http://citg.tudelft.nl/fileadmin/Faculteit/CiTG/Over_de_faculteit/Afdelingen/Afdeling_Waterbouwkunde/sectie_waterbouwkunde/people/personal/jonkman/doc/Jonkmanetal_flood_risk_south_holland.pdf
http://www.deltacommissaris.nl/
http://www.nfrmp.us/IFRMA/docs/presentations/vanAlphen_Presentation.pdf
Flood risk management in Europe: innovation in policy and practice
http://books.google.co.jp/books?id=A7THGJVD74EC&pg=PA481&lpg=PA481&dq=cost+effectiveness+dike++ring+level&source=bl&ots=LZVS63Iw52&sig=UqcBXEgdNNm3x4iYjlh3WeiZeNw&hl=zh-TW&ei=RemCTq7xLOmgmQWAjsFD&sa=X&oi=book_result&ct=result&resnum=2&ved=0CCoQ6AEwAQ

2011-09-27

ICFM5 note01 2011.09.27

100% reliable flood early warning system does not exit

http://www.siwi.org/documents/Resources/Water_Front_Articles/2009/WRMASystemsView.pdf

http://www.deltacommissaris.nl/english/

Flood Ensemble forecasting are using in Australia ,Canada, and Japan

2011-09-22

牛鬥橋地面光達 3D 資料

利用地面光達資料重新產出之3d成果

並且輸出5m dem資料

本次成果較前一次測試成果增加一個測站資料

National Taipei University of Technology
Jhih-Cyuan Shen
臺北科技大學土木工程系
沈志全

2011-09-07

FEWS_Taiwan讀取ESRI SHP屬性展示設定

若要顯示中文字串需要新增<charset>使FEWS_Taiwan知道 SHP檔的文字為哪一種編碼

在<geoDatum >標籤之後
<visible>標籤之前

<esriShapeLayer id="CWB_Meteo67" name="CWB_Meteo67">
<file>%REGION_HOME%/MapLayerFiles/CWBMeteo_67.shp</file>
<geoDatum>TWD 1967</geoDatum>
<charset>Cp950</charset> <--! 中文編碼設定-- >
<visible>false</visible>
<label>%NAME%</label>
<labelFontSize>20</labelFontSize>
<toolTip>%ID%</toolTip>
<pointSize>20</pointSize>
</esriShapeLayer>

2011-09-01

FEWS_Taiwan在二維網格圖層之上新增測站位置與資訊

FEWS_Taiwan可以將中央氣象局提供雷達降雨二維網格資訊進行計算與展示功能

但如何在二維網格圖層之上新增測站位置與資訊

需利用新增圖層方式將測站圖層新增於網格展示中

xml 設定方法如下
<esriShapeLayer id="CWB_Meteo" name="CWB_Meteo">
<file>%REGION_HOME%/MapLayerFiles/CWBMeteo_84.shp</file> 
<geoDatum>WGS 1984</geoDatum> 
<visible>false</visible>
<label>%ID%</label> 
<labelFontSize>20</labelFontSize>
<pointSize>20</pointSize>
</esriShapeLayer>

其成果如圖所示

2011-08-30

2011-08-30 NANMADOL颱風恆春淹水事件雨量站空間位置關係圖

南瑪都颱風

恆春淹水事件

雨量站空間位置關係圖

淹水事件

2011-08-27 00:00 GMT+8時間開始至 08-30 15:00截止
佳樂水最大降雨 72mm/hr 累積雨量1018 mm
墾丁最大降雨127mm/hr 累積雨量827 mm
恆春最大降雨67 mm/hr 累積雨量975.5 mm

2011-08-24

0820研究室聚餐

2011-08-22

FEWS_Taiwan工作日誌20110822

1.台灣颱風洪水研究中心（颱洪中心）颱風定量降雨數值模式系集預報實驗預報降雨量銜接測試完成
在Delos（柏愷）檔案轉換流程的協助下完成颱洪中心預報降雨量銜接測試
由於氣象投影轉換至規則網格之問題仍無法有效解決，目前仍舊採用強制空間對應方式
將雨量網格指定在對應之規則網格中，目前已經完成單一組預報成果之銜接測試
後續將持續針對各系集預報成員成果內容加以轉換以完成洪水系集預報之測試工作

2.利用氣象局qpesums雷達降雨資料進行水利署雨量測站成果輸出之測試工作
利用北區水資源局石門水庫之雨量站進行測試
雷達降雨量資料使用最接近該雨量測站之網格方式加以指定
利用凡那比颱風為例
初步可以發現在整場颱風豪雨事件兩種觀測量誤差最小2.9mm 最大誤差108.4mm
霞雲34.4mm與高義108.4mm之雨量出現較大的差異
其他雨量站鎮西堡石門白石玉峰池端巴陵嘎拉賀總降雨量誤差在2.9~12.8mm之間
由於兩種資料觀測方式在觀測方式與理論上屬於不相同之方式
需要近一步探討相關空間與資料之差異與確認儀器之正確性

國立臺北科技大學,土木工程系,光達與高效能試算防災中心,經濟部水利署

2011-08-21

心態與想法

這星期在工作時聽到一位在政府部門的資訊前輩的談話內容
將談話內容彙整如下

在工作上
-每天開心來上班工作
-不斷提升自己本職學能
-知道自己在作什麼要完成什麼

在心態上
-不要用我認為我以為等抗拒的心態來面對新事物
-不管怎麼樣使用者都是對的（不論使用者的操作或流程是否有問題）
-試著使用多個面相看待事情
-試著改變自己的看法

在學習上
-每一天至少花30分鐘學習新東西

提供給大家共勉之

另外如果公務人員都有這樣的心態國家的進步應該不只是現在這種情況

2011-08-18

FEWS_Taiwan20110818工作日誌

1.完成 soebk模式新增計算點位於FEWS_Taiwan平台之教學文件製作
2.完成後龍溪 soebk模新增測站輸出程序 1350H020;彼岸橋 1350H019;玉清大橋 1350H018;社寮角橋輸出
中港溪 soebk模新增測站輸出程序 1340H010;尖山大橋 1340H011;東興橋輸出

3.教學文件內容規劃

標題：教學文件標題
說明：針對教學文件內容之概念敘述
更新流程概述：說明主要更新流程
相關檔案：可能調整之相關檔案
詳細流程：針對流程進行細部說明
相關文件：相關可參考之文件

後需相關教學文件內容將參照上述架構新增並且針對既有文件進行調整與配置

2011-08-14

降雨量資訊問題

雷達降雨量資訊與測站降雨量資訊出現差異

凡那比事件資料在高屏溪流域月眉大橋上出現計算成果之差異

使用雷達平均降雨量資訊之成果出現小於測站內查計算之觀測降雨量

有兩張集水區降雨資訊中可以看出山區之集水區降雨量可能出現較大之差異性

在時間許可下應該可以近一步分析出差異量與改善或修正計算之流程

2011-07-28

FEWS_Taiwan工作日誌20110728

上傳新版展示filter於mc主機

針對海象資訊、氣象與降雨量資訊、區域雨量資訊等內容進行調整
詳細資訊內容如下

海象資訊
--潮位
--資料浮標
--河口模擬潮位
氣象與降雨量資訊(測站)
--氣象局氣候站
--氣象局雨量站
區域雨量資訊（面）
--台灣各鄉鎮分區降雨
--水位站上游集水區
--水庫集水區

並且將地下水分區完成分區設定工作其內容如下

地下水分區
--臺北盆地
--桃園中壢臺地
--新苗地區
--臺中地區
--濁水溪沖積扇
--嘉南平原
--屏東平原
--蘭陽平原
--花東縱谷
--澎湖地區

CWB WRF vs TTFRI WRF Grid Location Different

最外圍網格為中央氣象局所提供之WRF資料網格資料
較小網格為颱洪中心所提供之WRF資料網格資料

資料解析度約為 5 KM

TTFRI WRF Grid point

配合TTFRI 供應的WRF系集預報資料

重新進行網格點位確認工作

目前在資料中包含3354個網格點位

空間解析度約5~5.1公里

時間解析度1小時

預報資料長度 78HR

每天固定輸出四次(00,06,12,18)每次間隔6HR (時間標準GMT+0 )

資料開始輸出至輸出完成約需耗費1.5~3HR

如附圖所示為TTFRI WRF Grid Point 資訊

2011-07-15

FEWS_Taiwan工作日誌20110714

1.完成網格型態資料轉換成為點位資料測試
可將網格型態的資料內容利用空間內插方式取得特定點位之數值

預計將針對

2004艾利 AERE 08/23∼08/26（只有測站資訊）
2005馬莎 MATSA 08/03∼08/06
2005海棠 HAITANG 07/16∼07/20
2005泰利 TALIM 08/30∼09/01
2007柯羅莎 KROSA 10/04∼10/07
2008 鳳凰 FUNG-WONG 07/26∼07/29
2008辛樂克 SINLAKU 09/11∼09/16
2008薔蜜 JANGMI 09/26∼09/29
2009 莫拉克 MORAKOT 08/05∼08/10
等颱風進行測試

2.取得颱洪中心系集預報降雨資料
開始進行資料匯入轉換程式之撰寫工作

3.後續資料傳輸工作之方式
協助外單位進行資料轉換傳輸位置只限定於所屬機關位置網段
協力單位自行向所屬機關取得

2011-07-07

FEWS_Taiwan工作日誌 20110707

完成水利署浮標資料更新工作

將浮標資料區分成為水利署與氣象局兩個分群

氣象局浮標資訊目前可以取得
觀測浪高
波浪週期
波浪方向

水利署浮標資訊目前可以取得
觀測流速
觀測流向
觀測水溫
乾球觀測溫度
觀測氣壓
觀測風速
觀測三秒陣風風速
觀測浪高
波浪週期
波浪方向

附圖為將水利署浮標資訊之觀測數據同時展示之畫面

2011-07-06

0626 三重淹水

水深約10~15cm

QGIS 1.7 版釋出

QGIS 1.7 版在官方網站上已經釋出

http://www.qgis.org/wiki/Download

QGIS 1.7.0 was released on 19 June 2011

檔案容量91.8mb

專案名稱 Wrocław 弗羅茨瓦夫（波蘭地名）

在圖層屬性調整的介面進行大幅度的變動

將原先選項置於左側垂直方式選擇調整為置於上方水平選擇

如圖所示分別為1.6 與1.7 版介面之差異

由說明文件中指出修正277個bug 和新增許多新功能

有興趣的人可以試試看

2011-06-24

CWB QPE Rainfall with different basemap

在FEWS_Taiwan網格式資料展示功能中

將 wms 圖層加入顯示

如圖所示分別為利用微軟 VE 與open street map 之圖層作為基本底圖

將可以有效將降雨資訊與地理空間產生高度關聯性

訂閱：文章 (Atom)

C0G880	二水	209335.3	2634550	92
C0AD10	八里	289021.5	2782903	27
C0F9K0	大安	206551.9	2693773	30
C0F9N0	大里	217877.1	2665717	75
C0G840	北斗	198917.1	2641941	32
C1D640	外坪	257564.6	2726796	718
C0F9Q0	外埔	218475	2693976	149
C0G850	田中	206438	2639646	62
C0G820	田尾	198403.2	2644622	29
C0F9L0	后里	220904.7	2689824	253
C0E810	西湖	222860.3	2716346	53
C0G890	伸港	195937	2671999	24
C0G780	秀水	197878.7	2659299	27
C0K291	宜梧	163516.5	2604260	5
C0R430	東港	191711.1	2485355	31
C0G860	社頭	206210.9	2644154	45
C0G790	花壇	202788.5	2657273	27
C0G870	芳苑	179116.3	2646983	12
C0F9U0	南屯	211578.2	2670657	74
C0E750	苗栗	230578.5	2718001	62
C1D630	飛鳳山	260368.3	2740326	462
C0G810	埔心	200072.1	2649631	26
C0G800	埔鹽	190525.6	2655589	15
C0F9S0	烏日	210846.1	2667183	43
C0G830	埤頭	192617.9	2640700	27
C0AD20	深坑	310659.8	2766460	70
C0F9P0	清水	203913.1	2690080	26
C0E590	通霄	218400.8	2707719	40
C0E740	造橋	234652.5	2725980	20
C0V490	新興	176245.5	2503600	20
C0E820	獅潭	240252.8	2715151	220
C0G770	福興	191148.8	2660106	15
C0E780	銅鑼	227141.1	2709684	160
C1F9W0	德基	267654.1	2683670	1970
C0F9O0	潭子	218276.1	2679046	175
C0F9R0	龍井	200489.8	2675942	20
C0F9M0	豐原	219961.5	2683645	230
C0AD30	蘆洲	295990.9	2775879	20
C0AC80	文山	306453.4	2766565	40