[分享] 自己寫的程式 (TCFS V1.00 beta)

看板TY_Research作者 (小多啦)時間13年前 (2012/06/01 03:05), 編輯推噓14(14023)
留言37則, 16人參與, 最新討論串1/1
TC生成模擬程式 (TC Formation Simulator) V1.00 測試版 這是以去年開始寫的TCFS系列程式(之前寫到V0.09)為基礎 經過了半年多以來所改寫而成的程式,簡稱為TCFS V1.00 beta beta的意思就是現在只是先放出來"公開測試"的 裡面的參數或許還沒有調整到最好 等多蒐集一點資料後 可能今年暑假會發布TCFS V1.00 正式版 (參數設定會經過優化) 以下是對這個程式的基本介紹 (之前雖然就貼過了,但後面會新增"有更新"的地方) ------------------------------------------------------------------------- 這個軟體跟小程式 (颱風路徑預測程式)是不同的軟體 裡面所使用的模型也幾乎不同 這個程式顧名思義是偏重在預測熱帶氣旋(TC)的生成 而且使用者不用輸入任何參數 完全讀取圖形 判斷高低壓分佈 甚至就算現在西太平洋上沒有任何擾動或TC 這個程式也可以預測出幾個小時後會在哪邊出現新的擾動 然後會不會接著發展起來 先介紹一下使用方法 請先點以下網址 把這幾張圖按右鍵存成BMP檔案 (檔名不變 只是從GIF變成BMP檔案) 應該就是256色的BMP圖片檔案了 只需要下載最新的圖就好 渦度 http://tropic.ssec.wisc.edu/real-time/westpac/winds/wgmsvor.GIF
海溫 http://www.oceanweather.com/data/NPAC-Western/sst.html 風切 http://tropic.ssec.wisc.edu/real-time/westpac/winds/wgmssht.GIF
幅散 http://tropic.ssec.wisc.edu/real-time/westpac/winds/wgmsdvg.GIF
幅合 http://tropic.ssec.wisc.edu/real-time/westpac/winds/wgmsconv.GIF
最低層的駛流圖 (目前) http://tropic.ssec.wisc.edu/real-time/westpac/winds/wgmsdlm1.GIF
最低層的駛流圖 (12小時前) http://tropic.ssec.wisc.edu/real-time/westpac/winds/wgmsdlm1-4.GIF
PS:如果主要想分析的對象是已經有一定強度的TC的話 可以搭配適合該系統強度 不同層面的駛流圖 然後檔名再改回 wgmsdlm1.bmp 和 wgmsdlm1-4.bmp 就好 然後再上CWB的網站 下載你所要分析時間點的全球衛星雲圖 http://www.cwb.gov.tw/V7/observe/satellite/Sat_EA.htm# 要選擇色調強化圖才行 按右鍵另存成BMP檔案 檔名請改成 satellite (加上副檔名全名為 satellite.bmp) 但要注意的是 從這裡另存下來的BMP檔案會是24位元的BMP檔案 而這個程式要讀的是256色的BMP檔案 所以請再用類似小畫家的軟體 改存成256色的BMP檔案 然後此時小畫家會跳出 "存成這種格式會造成色彩遺失 是否繼續" 請選"是" 其他軟體也是一樣 反正就是要硬存成256色的BMP檔案就對了 然後這個程式的模擬範圍為 3N~58N 100~170E 強度方面是以氣壓hPA的方式表示的 預測的時間點最多可以從+0H跑到+240H 也就是未來10天內 把上面的那些圖(總共8個)都丟進去跟這個程式同一個資料夾後 只要點這個程式就會開始執行 當目前進度進展到第17個方塊時 就代表完成模擬了 一般來說大約要跑5~15分鐘 如果說執行到一半程式出錯的話 最有可能的原因就是那8張圖有問題 一定要確定都是BMP 256色執行檔案 執行結果大約長成這樣: http://ppt.cc/Vz8f 其中有些低壓擾動的起始時間不是 +0H 而是其他時間 例如+78H 代表說這個低壓擾動 可能會在78小時後 才在這邊出現 載點如下: TC Formation Simulator V1.00 beta http://www.mediafire.com/?xggchhz3d82yud9 另外補充說明一下 這個程式主要是用來預測哪裡可能會有擾動會TC生成 所以在路徑預測的準確度上可能是輸小程式 (颱風路徑預測程式V4.03) 所以如果說想預測的是目前已知的TC (TD強度以上)的話 用小程式反而比較好 而且TCFS這個程式還正在測試中 所以強度預測僅供參考 相信隨著未來多蒐集點數據之後 借由這些數據去微調程式裡的模型 強度預測的準確度應該會跟著提高 ----------------------------------------------------------------- 前面部分是以前就貼過的操作說明 接下來的部分是V1.00 beta版所更新的部分 1.強度模式有經過調整 現在比較不會像以前那樣常常高估TC強度 2.利用渦度圖來判斷短波槽的位置 (這點是舊版本所忽略的) 對判斷副高勢力的變化有幫助 3.我有做2000~2011年所有颱風的統計資料 共200多個颱風 上萬筆的各時間點詳細資訊 包括 颱風所在經緯度 強度 行進方向 速度等 並且以此為基準 以我自己的方法建立了一整套的颱風資料庫 所以TCFS再進行路徑預測時 會有部分的比例是參考這套內建資料庫的資料 對預測的路徑進行修正 進而提升路徑預報的準確度 4.加入了新的模式 "TCFS S+ mode" 其中S+代表Statistic plus的意思 所以以後TCFS會同時跑兩個模式:TCFS V1.00 mode 和 TCFS S+ mode 兩者的差別在於 TCFS S+ mode 比較注重過去的統計資料 (颱風資料庫) 會以此對路徑進行大幅度的修正 而TCFS V1.00 mode 雖然也是會參考資料庫 不過只是微調而已 修正的程度會比TCFS S+ mode還要少很多 因此對於大氣環境較單純 且颱風路徑較單純的情況下 TCFS S+ mode的準確度會比較高 (畢竟很注重過去的統計數據) 但是大氣環境只要一複雜 或許就會變成TCFS V1.00 mode比較準 5.把結果自動存成txt文字檔 (TCFS V1.00 beta.txt) 而且還有寫成kml檔形式的預測結果 可以用GOOGLE Earth開啟kml檔 就會有已經畫好的路徑圖 關於kml檔的詳細操作說明可以參考載點裡面的附件 (KML manual.doc) 如果說大家對使用上有任何疑問的話 可以寄站內信給我 或在以下推文 那就先這樣了~祝大家玩得愉快 ~^^ -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.112.218.91

06/01 05:36, , 1F
推 謝謝分享
06/01 05:36, 1F

06/01 08:57, , 2F
推認真~而且可以寫出這個真的很厲害 ^^
06/01 08:57, 2F

06/01 09:36, , 3F
06/01 09:36, 3F

06/01 09:36, , 4F
GOOD
06/01 09:36, 4F

06/01 12:27, , 5F
另一個「小程式」
06/01 12:27, 5F

06/01 12:50, , 6F
推,靠圖就能預測TD走向!太強了!!
06/01 12:50, 6F

06/01 22:15, , 7F
這麼熱心不給個推太說不過去了
06/01 22:15, 7F

06/01 23:26, , 8F
推強者!
06/01 23:26, 8F

06/02 00:28, , 9F
06/02 00:28, 9F

06/02 23:33, , 10F
能想到把圖轉成bmp當作bin檔來讀也蠻厲害的,不知道
06/02 23:33, 10F

06/02 23:34, , 11F
是不是資工系影像處理的背景~不過我更好奇的是
06/02 23:34, 11F

06/02 23:34, , 12F
Governing Equation的架構是什麼XD因為熱帶擾動除了
06/02 23:34, 12F

06/02 23:35, , 13F
除了高頻的波動會影響以外,從模式邊界吃進來的低頻
06/02 23:35, 13F

06/02 23:35, , 14F
波動也會影響擾動的生成位置
06/02 23:35, 14F

06/02 23:57, , 15F
我其實不是資工系的XD 只是有修過一學期的C++語言
06/02 23:57, 15F

06/02 23:58, , 16F
當然要做數值模式 一定會用到很多複雜的公式和微分
06/02 23:58, 16F

06/02 23:59, , 17F
方程 如果每個式子都要求出精確的解的話 那可能就
06/02 23:59, 17F

06/03 00:00, , 18F
不是憑個人的力量寫的出來的 可能還要用到超級電腦
06/03 00:00, 18F

06/03 00:00, , 19F
然後必須取得大量的初始場資料 (通常個人無法取得)
06/03 00:00, 19F

06/03 00:01, , 20F
因此基於EE背景 最喜歡做的事情就是 Approximation
06/03 00:01, 20F

06/03 00:02, , 21F
也就是做一大堆假設來化簡複雜的計算 當然這些假設
06/03 00:02, 21F

06/03 00:03, , 22F
必須要跟實際的情況不會相差太遠
06/03 00:03, 22F

06/03 00:04, , 23F
所以這個程式雖然準確度不會太高 但是適用個人電腦
06/03 00:04, 23F

06/03 00:05, , 24F
就可以跑的 (所需時間5~15分鐘) 而且使用的初始場資
06/03 00:05, 24F

06/03 00:05, , 25F
料也都是網路上就可以取得的 這才是我寫這些程式的
06/03 00:05, 25F

06/03 00:06, , 26F
目的 要讓對大氣有興趣的一般人 也能享受跑數值的
06/03 00:06, 26F

06/03 00:07, , 27F
樂趣 ^^
06/03 00:07, 27F

06/03 00:35, , 28F
其實現在全世界的天氣模式也都是參數化的結果啦XD除
06/03 00:35, 28F

06/03 00:35, , 29F
了日本的NICAM model以外
06/03 00:35, 29F

06/04 00:09, , 30F
真強,推一個
06/04 00:09, 30F

06/04 00:12, , 31F
推!
06/04 00:12, 31F
載點更新了 修正輸出時不同系統間的空行問題 有空行會比較方便閱讀 ※ 編輯: DoraBoy 來自: 140.112.218.91 (06/04 22:06)

06/04 22:07, , 32F
如果想分析的對象是有一定強度的TC 則可以嘗試看看
06/04 22:07, 32F

06/04 22:07, , 33F
搭配適合該系統強度 那個層面的駛流圖 然後再把檔名
06/04 22:07, 33F

06/04 22:08, , 34F
改回 wgmsdlm1.bmp 和 wgmsdlm1-4.bmp 這樣可能較準
06/04 22:08, 34F

06/06 23:50, , 35F
大推!
06/06 23:50, 35F

07/31 01:33, , 36F
07/31 01:33, 36F

08/13 23:04, , 37F
目的 要讓對大氣有興趣 https://noxiv.com
08/13 23:04, 37F
文章代碼(AID): #1Fny3f4- (TY_Research)