2013-06-27

QGIS add Grass Create Watershed

3 step 

1. Import Raster
r.in.gdal -o input=G:/giswork/7thrmo/1850_b5km.tif output=1850

2. Create Watershed
r.watershed elevation=1850@wgs84 threshold=125 accumulation=9 drainage=10 stream=11 basin=12

3.raster to vector

r.to.vect input=12@wgs84 output=v1850 feature=area -s


 

Department of Civil Engineering 
National Taipei University of Technology
Jhih-Cyuan Shen

臺北科技大學土木工程系
沈志全
 

2013-06-25

氣象局五分山雷達維修保養時間

2013/06/24
 氣象局五分山氣象雷達站因清洗天線罩,於6月27/28兩日08時00分起至18時00分停止觀測,如提早完成,則提前恢復觀測。
2013/06/24
氣象局五分山氣象雷達站因月保養,自25日13時00分起至25日18時00分停止觀測。

2013-06-13

QGIS Intersect (相交) 切圖

在 QGIS中利用 Intersect 功能可以進行兩個圖層切割功能 (如圖1)

input 圖層為 基本主要圖層

 Intersect圖層 為要用來切割的圖層

例如將全台水位站上游集水區搭配全台氣象局徐昇多邊形進行相交套疊切割

找出各水位站上游集水區的所包含的氣象局徐昇多邊型範圍與邊界





    













圖1 QGIS  Intersect  功能


水位站上游集水區邊界
全台氣象局雨量站徐昇多邊形

兩個圖層套疊
切割後的成果


Department of Civil Engineering
National Taipei University of Technology
Jhih-Cyuan Shen


2013-06-01

GMTED2010數值高程資料

GMTED2010

Global Multi-resolution Terrain Elevation Data 2010 

美國 U.S. Geological Survey (USGS) 及 National Geospatial-Intelligence Agency (NGA)

提供的新的全球數值高程資料取代原先的GTOPO30 資料

目前資料解析度如果使用7.5秒 的資料在台灣大約是180~200公尺一個點

適合進行初步的水文模式分析與基本誌突的需求

在官方網站中有提到 全球資料在垂直精度於控制點上可達到6M的精度

資料內容區分為30秒 15秒 7.5秒等三種解析度
 整體而言 
30秒資料 高程誤差之RMSE 可能介於 25 ~ 42 meters
15秒資料 高程誤差之RMSE 可能介於 29 ~ 32 meters
7.5秒資料 高程誤差之RMSE 可能介於 26 ~ 30 meters

詳細相關資訊可以由https://lta.cr.usgs.gov/GMTED2010 取得更完整的說明資訊


台灣鄰近區域資料查詢


              流域圖層套疊成果

QGIS Raster to Vector use Python console step2

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