2008-06-30

fews表單內容更新與建置工作

FEWS_Taiwan資料內容整並與調整
儘速完成下列表單之更新工作
1.苗栗縣、新竹縣水情通報單
2.苗栗縣、新竹縣洪水警報單
3.氣象局雨量站降雨量現況表
4.水位站水位現況表單
5.水文分析表單樣版建立

2008-06-29

Using_Arcgis_Command_line_syntax_to_get_slope_value

1.匯入asc檔並轉換成raster檔案
使用ASCIIToRaster_conversion H:\NCTU\5m_asc\0_5m.asc H:\NCTU\5m_asc\16_5m FLOAT
語法:ASCIIToRaster_conversion <in_ascii_file> <out_raster> {INTEGER | FLOAT}
2. 計算格網點坡度
使用Slope_sa 16_5m H:\NCTU\5m_asc\16_slope2 DEGREE 1
語法Slope_sa <in_raster> <out_raster> {DEGREE | PERCENT_RISE} {z_factor}
3.將格網資料輸出成為點資型態
使用RasterToPoint_conversion 16_slope2 H:\NCTU\5m_asc\p_16_slope2.shp Value
語法 : RasterToPoint_conversion <in_raster> <out_point_features> {raster_field}
4.提供xyz格式資料給其它模式使用
使用ExportXYv p_16_slope2 GRID_CODE Comma H:\NCTU\5m_asc\16p_s.txt
語法ExportXYv <Input_Feature_Class> <Value_Field> <Space | Comma |
Semi-Colon> <Output_Ascii_File>

2008-06-10

FEWS Taiwan Modules WAR Rainfall 小時雨量計算模組更新完成

將水利署非等距之觀測降雨量轉換成為等距之每小時降雨量資訊

ruleBasedTransformation rule採用 nonequitoequidistantwithzero

將非等距轉換為等距資料

在利用aggregate rule的accumulative方式將資料累計成為每小時雨量資訊

成果如圖所示

相關ruleBasedTransformation rule 語法如下
sampleforecastvalues,
selectpeakvalues,
selectlowvalues,
selectpeakvalueswithincertaingap,
selectlowvalueswithincertaingap,
equitononequidistant,
equitononequidistantforinstantaneousseries,
equitononequidistantforaccumulativeseries,
datahierarchy,
typicalprofiletotimeseries,
nonequitoequidistant,
nonequitoequidistantforinstantaneousseries,
nonequitoequidistantforaccumulativeseries,
nonequitoequidistantwithzero,
nonequitoequidistantwithmissing,
zerodegreealtitudelevel,
datatotimeseries,
shiftdata,
timeshift,
timedifference,
soilmoisture,
opticalDepth,
harmonictotimeseries,
selectmaxvalue,
selectminvalue,
convertDatumLocalToGlobal,
convertDatumGlobalToLocal,
elevationlapse,
samplehistorical

FEWS Taiwan Modules 計算累積降雨量

使用FEWS Taiwan 累積降雨計算

aggregate rule採用accumulative

<outputVariable variableId="1hrCWB">
<outputVariable variableId="3hrCWB">
<outputVariable variableId="6hrCWB">
<outputVariable variableId="12hrCWB">
<outputVariable variableId="24hrCWB">
<outputVariable variableId="48hrCWB">

workflowDescriptor 中 forecast設定為true

完成後即可取得如圖之累計雨量資訊

FEWS_Taiwan Modules problem in sobek batch

在第一次安裝系統時
若沒有重新註冊dll檔將無法藉由fews啟動sobek
在\diagnostics\Diagnostics.xml檔案中會出現以下訊息
Error while loading file in routine Load_XML_File: ActiveX component
can't create object Operation Cancelled
解決方式
先將
SBKAdvBatchLib.dll 589824 byte
LtXmlComLib5.dll 937984 byte
LtXmlComHelp5.dll 266240 byte
等三個檔案找到原先舊版的檔案位置加以置換
並且利用 regsvr32 將dll重新註冊
regsvr32 SBKAdvBatchLib.dll
regsvr32 LtXmlComLib5.dll
regsvr32 LtXmlComHelp5.dll
完成後即可經由fews啟動sobek

注意執行上述流程時時要使用系統管理員身份的權限執行



2008-06-09

工作日誌20080608

預計更新天文潮位資料測站名稱
1596;大武
1386;小琉球
1196;後壁湖
123;梗枋
1786;永安
1356;澎湖
111;竹圍
1676;綠島
1246;蘇澳
1396;蘭嶼

二河局所需防汛配合資訊
1.目前施工地點(廠商名稱 主辦協辦工程司 預計完工日期 完成百分比)
2.高風險度區 聯繫資訊
3.警戒水位資訊 (水位站警戒區域)
4.警戒雨量資訊
5.斷面資訊
6.防汛道路資訊
7.重現期距 水位雨量資訊
8.堤防位置 (名稱 坐標 )
9.防汛物資堆置場所 (名稱 品名 數量 )

防汛作業時 取得之相關資訊與數據如何因應與如何處理 ( 降雨量、水位 )
作業時間規劃 執行與研判所需時間需考量

QGIS Raster to Vector use Python console step2

資料處理步驟: 1.讀取 Raster 並計算 25 等級等值線級距 2.用 GDAL 產生等值線(Contour Lines) 3.Polygonize:Raster 轉等值面多邊形 4.將多邊形圖層載入 QGIS 5.為每個 polygon 計算其 VALUE 所屬級距(LE...