Link Search Menu Expand Document
Focus on Air Quality

LGHAP數據之下載及轉檔

Table of contents

下載與解壓縮

下載

  • 由[zenodo][zenodo]網站取得zip檔案之網址後，以規則化檔名進行下載。規則為
    • 目錄：https://zenodo.org/record/5652265/files/
    • 檔名：LGHAP.PM25.D001.Y$y.zip
for y in 20{00..20};do 
wget -q https://zenodo.org/record/5652265/files/LGHAP.PM25.D001.Y$y.zip;done &

解壓縮

  • 內容為日均值
  • 每天存一個檔案
  • 檔案規則：LGHAP.${spec}.D001.A${YYYYMMDD}.nc
for i in 20{00..20};do unzip *${i}*.zip;done

轉檔

檔案格式

  • 為自行定義之ncf檔案。(不是CF、ioaip、USGS等 convention)
  • 2個維度：先是經度、再來是緯度(需轉置)
  • 因使用正值作為nan(海上)之替代值，且該值似乎不同年度有異，需特別注意。
netcdf LGHAP.PM25.D001.A20151018 {
dimensions:
        lat = 4100 ;
        lon = 6800 ;
variables:
        float lat(lat) ;
        float lon(lon) ;
        ushort PM25(lon, lat) ;
                PM25:_FillValue = 65535US ;
                PM25:units = "micrograms per cubic meter (μg m**-3)" ;
                PM25:scale_factor = 1. ;
                PM25:add_offset = 0. ;

// global attributes:
                :_NCProperties = "version=2,netcdf=4.7.3,hdf5=1.8.12," ;
}

將LGHAP轉成ioapi濃度檔

  • 使用rdPM25.py(個別轉檔)或tmPM25.py(求取平均值並轉檔)，轉成CMAQ濃度檔(ioapi convention)；
  • 全年合併成1個檔案
  • 引數：年代
  • 詳見[程式說明][rdPM25.md]
for i in 20{00..20};do sub python rdPM25.py $i ;done

for i in 20{00..20};do sub python tmPM25.py $i ;done

程式碼下載