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

Meow · 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天：

https://re.ssec.wisc.edu/proxy/image.php?products=HIMAWARI-B03_20150707_020000&width=720&height=480&center=15,140&zoom=7&labels=lines
https://re.ssec.wisc.edu/proxy/image.php?products=HIMAWARI-B波道_年月日_時分秒&width=寬&height=高&center=北緯,東經&zoom=縮放&labels=lines

複製代碼

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

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

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

Merge Channels.png

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

帝汶海

存檔請點選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以外，其他和基礎版的說明一致。

width=720, height=480
width=寬, height=高

複製代碼

start_time=datetime(2015,7,7,2,0,0)
end_time=datetime(2015,7,7,2,0,0)
num = (end_time - start_time) // timedelta(minutes=10) + 1
time = start_time
for i in range(num):
#make_RGB(time.strftime('%Y-%m-%dT%H%M%SZ.png'), (15,140), 7, time)
#ssec_download(time.strftime('%Y-%m-%dT%H%M%SZ.png'), 'HIMAWARI-B13', (15,140), 6, time)
time = time + timedelta(minutes=10)
start_time=datetime(年,月,日,起始時,起始分,0)
end_time=datetime(年,月,日,終止時,終止分,0)
num = (end_time - start_time) // timedelta(minutes=10) + 1
time = start_time
for i in range(num):
#make_RGB(time.strftime('%Y-%m-%dT%H%M%SZ.png'), (北緯,東經), 縮放, time)
#ssec_download(time.strftime('%Y-%m-%dT%H%M%SZ.png'), 'HIMAWARI-B##', (北緯,東經), 縮放, time)
time = time + timedelta(minutes=10)