AeroViz 0.1.6__py3-none-any.whl → 0.1.8__py3-none-any.whl

AeroViz/process/script/retrieve_RI.py

@@ -1,69 +0,0 @@
-import numpy as np
-from pandas import DataFrame
-
-from AeroViz.process.core.SizeDist import SizeDist
-from AeroViz.process.method import Mie_PESD
-
-
-def retrieve_RI(_df: DataFrame,
-                _PNSD: DataFrame,
-                nMin: float = 1.33,
-                nMax: float = 1.60,
-                kMin: float = 0.00,
-                kMax: float = 0.60,
-                spaceSize: int = 31,
-                dlogdp: float = 0.014
-                ) -> DataFrame:
-    nRange = np.linspace(nMin, nMax, num=spaceSize)
-    kRange = np.linspace(kMin, kMax, spaceSize)
-    Delta_array = np.zeros((spaceSize, spaceSize))
-    # 同一時間除了折射率其餘數據皆相同 因此在折射率的迴圈外
-    bext_mea, bsca_mea, babs_mea = _df['Extinction'], _df['Scattering'], _df['Absorption']
-
-    dp = SizeDist(data=_PNSD).dp
-    for ki, k in enumerate(kRange):
-        for ni, n in enumerate(nRange):
-            m = n + (1j * k)
-            ndp = np.array(_df[3:])
-
-            ext_dist, sca_dist, abs_dist = Mie_PESD(m, 550, dp, ndp)
-
-            bext_cal = sum(ext_dist) * dlogdp
-            bsca_cal = sum(sca_dist) * dlogdp
-            babs_cal = sum(abs_dist) * dlogdp
-
-            Delta_array[ni][ki] = ((babs_mea - babs_cal) / 18.23) ** 2 + ((bsca_mea - bsca_cal) / 83.67) ** 2
-
-    min_delta = Delta_array.argmin()
-    next_n = nRange[(min_delta // spaceSize)]
-    next_k = kRange[(min_delta % spaceSize)]
-
-    # 將網格變小
-    nMin_small = next_n - 0.02 if next_n > 1.33 else 1.33
-    nMax_small = next_n + 0.02
-    kMin_small = next_k - 0.04 if next_k > 0.04 else 0
-    kMax_small = next_k + 0.04
-    spaceSize_small = 41
-
-    nRange_small = np.linspace(nMin_small, nMax_small, spaceSize_small)
-    kRange_small = np.linspace(kMin_small, kMax_small, spaceSize_small)
-    Delta_array_small = np.zeros((spaceSize_small, spaceSize_small))
-    # 所有數據與大網格一致所以使用上方便數即可
-    for ki, k in enumerate(kRange_small):
-        for ni, n in enumerate(nRange_small):
-            m = n + (1j * k)
-            ndp = np.array(_df[3:])
-            ext_dist, sca_dist, abs_dist = Mie_PESD(m, 550, dp, ndp)
-
-            bext_cal = sum(ext_dist) * dlogdp
-            bsca_cal = sum(sca_dist) * dlogdp
-            babs_cal = sum(abs_dist) * dlogdp
-
-            Delta_array_small[ni][ki] = ((bext_mea - bext_cal) / 18.23) ** 2 + ((bsca_mea - bsca_cal) / 83.67) ** 2
-
-    min_delta_small = Delta_array_small.argmin()
-    _df['re_real'] = (nRange_small[(min_delta_small // spaceSize_small)])
-    _df['re_imaginary'] = (kRange_small[(min_delta_small % spaceSize_small)])
-
-    print(f'\t\tReal part:{_df['re_real']}\tIm part:{_df['re_imaginary']}', end='')
-    return _df['re_real':]

AeroViz/rawDataReader/script/EPA_vertical.py

@@ -1,46 +0,0 @@
-import numpy as np
-from pandas import read_csv, to_numeric
-
-from AeroViz.rawDataReader.core import AbstractReader
-
-
-class Reader(AbstractReader):
-    nam = 'EPA_vertical'
-
-    def _raw_reader(self, file):
-        with file.open('r', encoding='ascii', errors='ignore') as f:
-            # 有、無輸出有效值都可以
-            # read 查詢小時值(測項).csv
-            df = read_csv(f, encoding='ascii', encoding_errors='ignore', index_col=0, parse_dates=True,
-                          usecols=lambda col: col != 'Unnamed: 1')
-
-            df.index.name = 'Time'
-            df.rename(columns={'AMB_TEMP': 'AT', 'WIND_SPEED': 'WS', 'WIND_DIREC': 'WD'}, inplace=True)
-
-            # 欄位排序
-            desired_order = ['SO2', 'NO', 'NOx', 'NO2', 'CO', 'O3', 'THC', 'NMHC', 'CH4', 'PM10', 'PM2.5', 'WS', 'WD',
-                             'AT', 'RH']
-
-            missing_columns = []
-
-            for col in desired_order:
-                if col not in df.columns:
-                    df[col] = np.nan
-                    missing_columns.append(col)
-
-            if missing_columns:
-                self.logger.info(f"{'=' * 60}")
-                self.logger.info(f"Missing columns: {missing_columns}")
-                self.logger.info(f"{'=' * 60}")
-                print(f"Missing columns: {missing_columns}")
-
-            df = df[desired_order]
-
-            # 如果沒有將無效值拿掉就輸出 請將包含 #、L、O 的字串替換成 *
-            df.replace(to_replace=r'\d*[#LO]\b', value='*', regex=True, inplace=True)
-            df = df.apply(to_numeric, errors='coerce')
-
-        return df
-
-    def _QC(self, _df):
-        return _df

AeroViz/rawDataReader/script/Table.py

@@ -1,27 +0,0 @@
-# read meteorological data from google sheet
-
-from pandas import read_csv, to_datetime
-
-from AeroViz.rawDataReader.core import AbstractReader
-
-
-class Reader(AbstractReader):
-    nam = 'Table'
-
-    def _raw_reader(self, file):
-        with file.open('r', encoding='utf-8-sig', errors='ignore') as f:
-            _df = read_csv(f, low_memory=False, index_col=0)
-
-            _df.index = to_datetime(_df.index, errors='coerce')
-            _df.index.name = 'time'
-
-            _df.columns = _df.keys().str.strip(' ')
-
-        return _df.loc[~_df.index.duplicated() & _df.index.notna()]
-
-    def _QC(self, _df):
-        # remove negative value
-        _df = _df.mask((_df < 0).copy())
-
-        # QC data in 6h
-        return _df.resample('6h').apply(self.basic_QC).resample(self.meta.get("freq")).mean()