2014-10-02

QGIS欄位更新語法

QGIS欄位更新語法

使用現有欄位字串資訊更新欄位內容



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


2014-09-04

雨量站坐標處理與清查

同樣雨量測站資料表面上可能有多種不同坐標及投影系統的資料來源

氣象局雨量站採用經緯度方式呈現的資料 

宣稱內容為GRS67 的資料但部份測站實際使用為WGS84

水利署雨量站採用二度分帶方式呈現的資料

但宣稱是 TWD97 但資料內容中可能有TWD67的資料

延續去年進行全台雨量站圖面坐標清查工作後今年再加入

水利署雨量站進行清查

處理完之後感覺頭一直呈現暈眩狀況

但這一個問題在國內至少有一半以上需要處理雨量站資料的人

在使用氣象局雨量站資料時可能都不曉得

需要將GRS67經緯度的資料轉換成為WGS84經緯度才能將測站

正確套疊在可能的坐標位置上

最快檢查的方式優先檢查台灣西海岸沿海的雨量站 如 枋寮  枋山 旗津 等雨量

如果上述雨量站套疊底圖位於海上表示目前使用的經緯度資料為GRS67的坐標系統

需要經過坐標轉換才能正確顯示所在位置

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

2014-08-30

二仁溪 累積降雨量網格呈現

利用20130828~20130901 QPESUMS降雨量資料進行
二仁溪累積降雨量空間 分布輸出展示


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

2014-06-10

農地空間資訊


農委會的農地空間資訊整合內容


有很多東西還沒完成但至少整合在一起還是值得參考

https://talis.coa.gov.tw/ASSO/


雖然台灣政府單位的首長通常不重視觀測資料與基礎的資料的建置

目前有些資料還是只有1996年的資料或更早之前的資料感覺有點..

是不是可以花多一點精神向長官說明基礎資料建置的重要性

少一辦些跳舞的活動 歌唱大賽 或者是煙火式的活動

就可以更新很多基礎資料內容

2014-06-02

台中市有沙坑的公園狀況



在台中市公園內可以讓小孩玩的沙坑屈指可數

北屯兒童公園 24.169985, 120.694483

(維護不佳)


英才兒童公園 24.154373, 120.671066

(目前是乎是維護較佳的沙坑 但有小黑蚊出沒 )

文英兒童公園 24.156369, 120.689629

(規劃最理想但 環境持續惡化中)

葫蘆墩公園 24.256707, 120.714549  有噴水池沙坑

( 有噴水池沙坑 但 周邊設施開始出現毀損的狀況 )


福安兒童公園 24.179677, 120.624662
(目前 維護狀況是乎是最佳的沙坑 )

但今天去了一個被台中市府廢棄的砂坑

地點在調查局台中市調查站對面綠園道的沙坑

坐標位置24.145967, 120.664710


沙坑出現積水(下面是水泥鋪面) 與長草的現象

這幾年臺中市政府對於市政建設真是"另人刮目相看"

首長喜歡的建設都是 預算是十億來百億去的建設

如 臺灣塔 媽祖園區 企鵝館 新社-大坑纜車

但對於花小錢就能維護好的與市民息息相關的公園

卻是疏於維護與管理

臺中市政府社會局兒童發展資訊網


提供的小朋友可以去的遊憩的地點 但是乎內容也是 錯誤百出


http://www.tc-rc.org/service_map.php?cate=37

上述相關沙坑地點空間分布請參考

https://mapsengine.google.com/map/edit?mid=z0Bw0l98s87U.kyiveEZvxges

2014-05-15

GIS檔案格式介紹GeoEAS

在GIS軟體中由於各家廠商由於對於軟體的銷售與保護措施或演算法考量

在GIS檔案格式的在使用上往往會有些保護與限制條件

例如如早年shp只能在esri軟體開啟隨著使用量大到一定程度esri開放shp相關解碼說明

使得目前各GIS軟體大多數都可以支援shp檔案的格式

目前在在QGIS軟體中即可支援大多數的GIS軟體的檔案格式

最近在處理水利署的檔案格式時

廠商檔案內容
如下

氣象局w0模式降雨預報: yyyymmddhh預報yyyymmddhh  解析度(度):  0.0125
3
Longitude
Latitude
intensity(mm/hr)
121.5125  25.2875    0.0
121.5375  25.2875    0.0
121.5500  25.2875    0.0
121.5625  25.2875    0.0
121.5750  25.2875    0.1
121.5875  25.2875    0.1
121.5000  25.2750    0.0
121.5125  25.2750    0.0

上述內容為水利署廠商供應的一組資料 依據廠商的說法該資料為pcraster檔案格式

但熟悉的pcraster預設檔案格式為*.map 而非文字檔

在網路上搜尋了一下 該檔案格式為pcraster支援的ascii檔案格式
Point data column file format
資料格式的名稱實際為GeoEAS file
實際為史丹佛大學發展的GSLIB使用的檔案格式
並非pcraster原本預設的輸出檔案格式


GeoEAS file相關網頁介紹與說明連結

http://www.gslib.com/gslib_help/format.html

https://www.gammadesign.com/GSWinHelp/importing_data_from_external_files/text_input_file_formats/text_input_file_formats__geoeas_format.htm

http://ag.arizona.edu/plp/GIS/Case_Study_Af/GeoEAS/fileformat.html

http://pcraster.geo.uu.nl/pcraster/4.0.1/doc/manual/secdatbase.html#secdatbasepointform




在FEWS_Taiwan時序列輸出時增加副檔名csv

在FEWS_Taiwan輸出時序列檔案時可以在檔案名稱內容的tag中
設定如 xxxxx.csv 的內容

詳細設定內容如圖所示










在輸出檔案時系統就可以自動將csv的副檔名寫入

ACTPdata_2010092013.csv

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

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

QGIS Raster to Vector use Python console step2

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