開啟輔助訪問
 找回密碼
 立即加入

合成向日葵8號衛星全彩雲圖教學

簽到天數: 588 天

[LV.9]以壇為家II

2015-5-13 21:13 | 顯示全部樓層
本帖最後由 Meow 於 2015-7-29 20:26 編輯

JMA向日葵8號衛星出類拔萃,不僅解析度翻倍,16個波道每10分鐘就截取1張,其中3個可見光波道還可合成全彩照片。但是你們想按照JMA官方指定方式,下載3個超大檔案、軟體解碼轉換、裁切出小部份、再合成出全彩照片嗎?一般要合成RGB照片起先想到的是Adobe Photoshop,但除非是工作或學業用途,你們會為了照片三合一而每月花費9.99美元嗎?

Noul 2015-05-10 0530Z.jpg
颱風紅霞巔峰照

Meow在此要讓大家簡簡單單合成向日葵8號衛星全彩雲圖,而且零成本、零負擔、不馬虎。美國國立衛生研究院基於Java開發了一套叫做ImageJ的跨平台公有領域軟體,可讓用戶操作一系列影響處理,包含這邊會用到的合成技術。如果沒有Adobe Photoshop,ImageJ是最好的選擇。

安裝好所需軟體後請下載由RealEarth™提供的圖檔,但只會保留3天:

  1. https://re.ssec.wisc.edu/proxy/image.php?products=HIMAWARI-B03_20150707_020000&width=720&height=480&center=15,140&zoom=7&labels=lines

  2. https://re.ssec.wisc.edu/proxy/image.php?products=HIMAWARI-B波道_年月日_時分秒&width=寬&height=高&center=北緯,東經&zoom=縮放&labels=lines
複製代碼


波道:可指定01到16,可見光RGB)分別對應030201,13是常用的紅外線波道。
時間:請使用UTC時間。
寬高:最高可達1136。
經緯:可達小數後1位,前面加負號即變成南緯或西經。
縮放:可見光波道可達7,其餘13波道可達6。
移除&labels=lines即移除邊線。

請務必確認當時系統中心位置再下載,最好分別命名rgb。如果覺得長寬1136像素不符需求,可下載4張位置不一樣的圖拼接,中心向外延伸5.5度成效最好。下載後請先比對3張圖有無偏差,如果偏差很明顯,合成後會有很多雜點甚至會出現彩虹,就得嘗試下載另一時間的圖檔。

下載後執行ImageJ,開啟那3個可見光雲圖的檔案,最好按照RGB順序開啟。點選Image > Color > Merge Channels...會出現以下視窗:

Merge Channels.png

請為C1、C2、C3分別選擇已開啟的RGB圖檔,下方3個選項勾不勾選都不影響合成效果,但最好勾選第2個。點選OK以後就會看到合成後的全彩圖像了,而且真的是彩色的。

RGB.png
帝汶海

存檔請點選File > Save As,會有一系列檔案格式可選擇,建議使用PNG,這是所有無損格式之中最小的。有時輸出後的圖檔看起來有層白色壓住,如果覺得不舒服可斟酌修改,Meow自己會加黑10%,如果太暗會增加曝光10%。

進階版

中國颱風論壇會員Nasdaq提供了基於Python寫成的腳本,可自動生成向日葵8號衛星全彩可見光雲圖。下載附件解壓縮後請放Generator.py到想要生成圖檔的檔案夾,然後以Python的IDLE執行Generator.py並點擊Run Module,就會開始生成圖檔。

Generator.py.zip (736 Bytes, 下載次數: 12)

請確認已安裝Python最新正式版並也安裝了Pillow,如果執行from PIL import Image語法不出錯就理應成功了。

以下語法可修改的部份以中文註釋,除時間前面不添加0以外,其他和基礎版的說明一致。

  1. width=720, height=480

  2. width=寬, height=高
複製代碼

  
  1. start_time=datetime(2015,7,7,2,0,0)
  2. end_time=datetime(2015,7,7,2,0,0)
  3. num = (end_time - start_time) // timedelta(minutes=10) + 1
  4. time = start_time
  5. for i in range(num):
  6.     #make_RGB(time.strftime('%Y-%m-%dT%H%M%SZ.png'), (15,140), 7, time)
  7.     #ssec_download(time.strftime('%Y-%m-%dT%H%M%SZ.png'), 'HIMAWARI-B13', (15,140), 6, time)
  8.     time = time + timedelta(minutes=10)

  9. start_time=datetime(年,月,日,起始時,起始分,0)
  10. end_time=datetime(年,月,日,終止時,終止分,0)
  11. num = (end_time - start_time) // timedelta(minutes=10) + 1
  12. time = start_time
  13. for i in range(num):
  14.     #make_RGB(time.strftime('%Y-%m-%dT%H%M%SZ.png'), (北緯,東經), 縮放, time)
  15.     #ssec_download(time.strftime('%Y-%m-%dT%H%M%SZ.png'), 'HIMAWARI-B##', (北緯,東經), 縮放, time)
  16.     time = time + timedelta(minutes=10)
複製代碼

移除其中一個#可選擇合成全彩圖像或從單一波道下載;起始時間和終止時間一致就只會生成單一圖像。

評分

參與人數 2水氣能量 +250 風力強度 +1 收起 理由
krichard2011 + 150 贊一個!
ben811018 + 100 + 1 感謝分享!!!

查看全部評分

回復

使用道具 舉報

簽到天數: 588 天

[LV.9]以壇為家II

Meow|2015-5-14 22:05 | 顯示全部樓層
還有進階版了,用Python自動生成的,不過要像第一張圖那麼大還是得自己拼接。:lol
回復

使用道具 舉報

簽到天數: 588 天

[LV.9]以壇為家II

Meow|2015-5-27 21:17 | 顯示全部樓層
向日葵8號衛星將在7月7日上午10點(02Z)正式啟用。
回復

使用道具 舉報

簽到天數: 588 天

[LV.9]以壇為家II

Meow|2015-6-6 23:30 | 顯示全部樓層
小幅調整Python語法,現在輸出圖片檔名是日期、時間了,不再是RGB或00開始。
回復

使用道具 舉報

簽到天數: 588 天

[LV.9]以壇為家II

Meow|2015-7-28 23:26 | 顯示全部樓層
透過和RealEarth團隊信件溝通得知波道圖像不重合的癥結在連結格式,現在微調後已不再有這種問題,請務必使用目前最新的辦法獲得圖像。
回復

使用道具 舉報

簽到天數: 241 天

[LV.8]以壇為家I

krichard2011|2015-7-29 09:56 | 顯示全部樓層
改了程式語法後 雖然影像較連貫但仍有差不多的問題
以下嘗試用python下載了昨天 下午兩點到三點的雲圖
但只有 14:10 與 15:00 能夠合成順利
其餘仍會缺圖檔 導致合成失敗
anigif.gif

點評

缺圖的問題仍無解,這只是解決不重合問題。  發表於 2015-7-29 20:35
回復

使用道具 舉報

快速回覆
您需要登錄後才可以回帖 登錄 | 立即加入

本版積分規則

Meow

    主題

    帖子

    28萬

    積分

    11級[狂風Violent storm]

    Rank: 11Rank: 11Rank: 11Rank: 11

11級[狂風Violent storm]

發訊息

本平台僅供學術討論之用,預報應以氣象局為準

威普網站虛擬主機贊助公司

臺灣第一個天氣類型社群平台 即時天氣資訊、精準颱風動態

線上客服
FB傳送訊息
廣告行銷
精準行銷 物超所值
官方粉專
發佈 快速回復 返回頂部 返回列表