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

AeroViz/process/script/Chemical.py

@@ -1,177 +0,0 @@
-from pathlib import Path
-
-import numpy as np
-from pandas import read_csv, concat, notna, DataFrame, to_numeric
-
-from AeroViz.process.core import DataProc
-from AeroViz.tools.datareader import DataReader
-
-
-class ChemicalProc(DataProc):
-    """
-    A class for process chemical data.
-
-    Parameters:
-    -----------
-    reset : bool, optional
-        If True, resets the process. Default is False.
-    filename : str, optional
-        The name of the file to process. Default is None.
-
-    Methods:
-    --------
-    mass(_df):
-        Calculate mass-related parameters.
-
-    volume(_df):
-        Calculate volume-related parameters.
-
-    volume_average_mixing(_df):
-        Calculate volume average mixing parameters.
-
-    process_data():
-        Process data and save the result.
-
-    Attributes:
-    -----------
-    DEFAULT_PATH : Path
-        The default path for data files.
-
-    Examples:
-    ---------
-
-    """
-
-    def __init__(self, file_paths: list[Path | str] = None):
-        super().__init__()
-        self.file_paths = [Path(fp) for fp in file_paths]
-
-    @staticmethod
-    def mass(_df):  # Series like
-        Ammonium, Sulfate, Nitrate, OC, Soil, SS, EC, PM25 = _df
-        status = (Ammonium / 18) / (2 * (Sulfate / 96) + (Nitrate / 62))
-
-        if status >= 1:
-            _df['NH4_status'] = 'Enough'
-            _df['AS'] = 1.375 * Sulfate
-            _df['AN'] = 1.29 * Nitrate
-
-        if status < 1:
-            _df['NH4_status'] = 'Deficiency'
-            mol_A = Ammonium / 18
-            mol_S = Sulfate / 96
-            mol_N = Nitrate / 62
-            residual = mol_A - 2 * mol_S
-
-            if residual > 0:
-                _df['AS'] = 1.375 * Sulfate
-                _df['AN'] = residual * 80 if residual <= mol_N else mol_N * 80
-
-            else:
-                _df['AS'] = mol_A / 2 * 132 if mol_A <= 2 * mol_S else mol_S * 132
-                _df['AN'] = 0
-
-        _df['OM'] = 1.8 * OC
-        _df['Soil'] = 28.57 * Soil / 1000
-        _df['SS'] = 2.54 * SS
-        _df['EC'] = EC
-        _df['SIA'] = _df['AS'] + _df['AN']
-        _df['total_mass'] = _df[['AS', 'AN', 'OM', 'Soil', 'SS', 'EC']].sum()
-        species_lst = ['AS', 'AN', 'OM', 'Soil', 'SS', 'EC', 'SIA', 'unknown_mass']
-
-        _df['unknown_mass'] = PM25 - _df['total_mass'] if PM25 >= _df['total_mass'] else 0
-        for _species, _val in _df[species_lst].items():
-            _df[f'{_species}_ratio'] = _val / PM25 if PM25 >= _df['total_mass'] else _val / _df['total_mass']
-
-        return _df['NH4_status':]
-
-    @staticmethod
-    def volume(_df):
-        _df['AS_volume'] = (_df['AS'] / 1.76)
-        _df['AN_volume'] = (_df['AN'] / 1.73)
-        _df['OM_volume'] = (_df['OM'] / 1.4)
-        _df['Soil_volume'] = (_df['Soil'] / 2.6)
-        _df['SS_volume'] = (_df['SS'] / 2.16)
-        _df['EC_volume'] = (_df['EC'] / 1.5)
-        _df['ALWC_volume'] = _df['ALWC']
-        _df['total_volume'] = sum(_df['AS_volume':'EC_volume'])
-
-        for _species, _val in _df['AS_volume':'ALWC_volume'].items():
-            _df[f'{_species}_ratio'] = _val / _df['total_volume']
-
-        _df['density'] = _df['total_mass'] / _df['total_volume']
-        return _df['AS_volume':]
-
-    @staticmethod
-    def volume_average_mixing(_df):
-        _df['n_dry'] = (1.53 * _df['AS_volume_ratio'] +
-                        1.55 * _df['AN_volume_ratio'] +
-                        1.55 * _df['OM_volume_ratio'] +
-                        1.56 * _df['Soil_volume_ratio'] +
-                        1.54 * _df['SS_volume_ratio'] +
-                        1.80 * _df['EC_volume_ratio'])
-
-        _df['k_dry'] = (0.00 * _df['OM_volume_ratio'] +
-                        0.01 * _df['Soil_volume_ratio'] +
-                        0.54 * _df["EC_volume_ratio"])
-
-        # 檢查_df['ALWC']是否缺失 -> 有值才計算ambient的折射率
-        if notna(_df['ALWC']):
-            v_dry = _df['total_volume']
-            v_wet = _df['total_volume'] + _df['ALWC']
-
-            multiplier = v_dry / v_wet
-            _df['ALWC_volume_ratio'] = (1 - multiplier)
-
-            _df['n_amb'] = (1.53 * _df['AS_volume_ratio'] +
-                            1.55 * _df['AN_volume_ratio'] +
-                            1.55 * _df['OM_volume_ratio'] +
-                            1.56 * _df['Soil_volume_ratio'] +
-                            1.54 * _df['SS_volume_ratio'] +
-                            1.80 * _df['EC_volume_ratio']) * multiplier + \
-                           (1.33 * _df['ALWC_volume_ratio'])
-
-            _df['k_amb'] = (0.00 * _df['OM_volume_ratio'] +
-                            0.01 * _df['Soil_volume_ratio'] +
-                            0.54 * _df['EC_volume_ratio']) * multiplier
-
-            _df['gRH'] = (v_wet / v_dry) ** (1 / 3)
-
-            return _df[['n_dry', 'k_dry', 'n_amb', 'k_amb', 'gRH']]
-
-    @staticmethod
-    def kappa(_df, diameter=0.5):
-        surface_tension, Mw, density, universal_gas_constant = 0.072, 18, 1, 8.314  # J/mole*K
-
-        A = 4 * (surface_tension * Mw) / (density * universal_gas_constant * (_df['AT'] + 273))
-        power = A / diameter
-        a_w = (_df['RH'] / 100) * (np.exp(-power))
-
-        _df['kappa_chem'] = (_df['gRH'] ** 3 - 1) * (1 - a_w) / a_w
-        _df['kappa_vam'] = np.nan
-
-    @staticmethod
-    def ISORROPIA():
-        pass
-
-    def process_data(self, reset: bool = False, save_file: Path | str = None) -> DataFrame:
-        save_file = Path(save_file)
-        if save_file.exists() and not reset:
-            return read_csv(save_file, parse_dates=['Time'], index_col='Time')
-        else:
-            df = concat([DataReader(file) for file in self.file_paths], axis=1).apply(to_numeric, errors='coerce')
-
-            df_mass = df[['NH4+', 'SO42-', 'NO3-', 'Optical_OC', 'Fe', 'Na+', 'Optical_EC', 'PM2.5']].dropna().apply(
-                self.mass,
-                                                                                                        axis=1)
-            df_mass['ALWC'] = df['ALWC']
-            df_volume = df_mass[['AS', 'AN', 'OM', 'Soil', 'SS', 'EC', 'total_mass', 'ALWC']].dropna().apply(
-                self.volume,
-                axis=1)
-            df_volume['ALWC'] = df['ALWC']
-            df_vam = df_volume.dropna().apply(self.volume_average_mixing, axis=1)
-
-            _df = concat([df_mass, df_volume.drop(['ALWC'], axis=1), df_vam], axis=1).reindex(df.index.copy())
-            _df.to_csv(save_file)
-
-            return _df

AeroViz/process/script/IMPACT.py

@@ -1,49 +0,0 @@
-from pathlib import Path
-
-from pandas import DataFrame, read_csv, concat
-
-from AeroViz.process.core import DataProc
-from AeroViz.tools.datareader import DataReader
-
-
-class ImpactProc(DataProc):
-    """
-    A class for processing impact data.
-
-    Parameters:
-    -----------
-    reset : bool, optional
-        If True, resets the processing. Default is False.
-    save_filename : str or Path, optional
-        The name or path to save the processed data. Default is 'IMPACT.csv'.
-
-    Methods:
-    --------
-    process_data(reset: bool = False, save_filename: str | Path = 'IMPACT.csv') -> DataFrame:
-        Process data and save the result.
-
-    save_data(data: DataFrame, save_filename: str | Path):
-        Save processed data to a file.
-
-    Attributes:
-    -----------
-    DEFAULT_PATH : Path
-        The default path for data files.
-
-    Examples:
-    ---------
-    >>> df_custom = ImpactProc().process_data(reset=True, save_filename='custom_file.csv')
-    """
-
-    def __init__(self, file_paths: list[Path | str] = None):
-        super().__init__()
-        self.file_paths = [Path(fp) for fp in file_paths]
-
-    def process_data(self, reset: bool = False, save_file: Path | str = None) -> DataFrame:
-        save_file = Path(save_file)
-        if save_file.exists() and not reset:
-            return read_csv(save_file, parse_dates=['Time'], index_col='Time')
-        else:
-            _df = concat([DataReader(file) for file in self.file_paths], axis=1)
-            _df.to_csv(save_file)
-            return _df

AeroViz/process/script/IMPROVE.py

@@ -1,161 +0,0 @@
-from pathlib import Path
-from typing import Literal
-
-from pandas import read_csv, concat, read_json
-
-from AeroViz.process.core import DataProc
-from AeroViz.tools.datareader import DataReader
-
-
-class ImproveProc(DataProc):
-    """
-    A class for process improved chemical data.
-
-    Parameters:
-    -----------
-    reset : bool, optional
-        If True, resets the process. Default is False.
-    filename : str, optional
-        The name of the file to process. Default is None.
-    version : str, optional
-        The version of the data process. Should be one of 'revised' or 'modified'.
-        Default is None.
-
-    Methods:
-    --------
-    revised(_df):
-        Calculate revised version of particle contribution.
-
-    modified(_df):
-        Calculate modified version of particle contribution.
-
-    gas(_df):
-        Calculate gas contribution.
-
-    frh(_RH, version=None):
-        Helper function to get frh values based on relative humidity (RH) and version.
-
-    process_data():
-        Process data and save the result.
-
-    Attributes:
-    -----------
-    DEFAULT_PATH : Path
-        The default path for data files.
-
-    Examples:
-    ---------
-    >>> df = ImproveProc(reset=True, filename='revised_IMPROVE.csv', version='revised').process_data()
-
-    """
-
-    def __init__(self, file_paths: list[Path | str] = None):
-        super().__init__()
-        self.file_paths = [Path(fp) for fp in file_paths]
-
-    @staticmethod
-    def frh(_RH):
-        _frh = read_json(Path(__file__).parent.parent.parent / 'plot' / 'utils' / 'fRH.json')
-        if _RH is not None:
-            if _RH > 95:
-                _RH = 95
-            _RH = round(_RH)
-            return _frh.loc[_RH].values.T
-
-        return 1, 1, 1, 1
-
-    def revised(self, _df):
-        def mode(Mass):
-            L_mode = Mass ** 2 / 20 if Mass < 20 else Mass
-            S_mode = Mass - L_mode if Mass < 20 else 0
-
-            return L_mode, S_mode
-
-        _frh, _frhss, _frhs, _frhl = self.frh(_df['RH'])
-
-        L_AS, S_AS = mode(_df['AS'])
-        L_AN, S_AN = mode(_df['AN'])
-        L_OM, S_OM = mode(_df['OM'])
-
-        _df['AS_ext_dry'] = 2.2 * 1 * S_AS + 4.8 * 1 * L_AS
-        _df['AN_ext_dry'] = 2.4 * 1 * S_AN + 5.1 * 1 * L_AN
-        _df['OM_ext_dry'] = 2.8 * S_OM + 6.1 * L_OM
-        _df['Soil_ext_dry'] = 1 * _df['Soil']
-        _df['SS_ext_dry'] = 1.7 * 1 * _df['SS']
-        _df['EC_ext_dry'] = 10 * _df['EC']
-        _df['total_ext_dry'] = sum(_df['AS_ext_dry':'EC_ext_dry'])
-
-        _df['AS_ext'] = (2.2 * _frhs * S_AS) + (4.8 * _frhl * L_AS)
-        _df['AN_ext'] = (2.4 * _frhs * S_AN) + (5.1 * _frhl * L_AN)
-        _df['OM_ext'] = (2.8 * S_OM) + (6.1 * L_OM)
-        _df['Soil_ext'] = (1 * _df['Soil'])
-        _df['SS_ext'] = (1.7 * _frhss * _df['SS'])
-        _df['EC_ext'] = (10 * _df['EC'])
-        _df['total_ext'] = sum(_df['AS_ext':'EC_ext'])
-
-        _df['ALWC_AS_ext'] = _df['AS_ext'] - _df['AS_ext_dry']
-        _df['ALWC_AN_ext'] = _df['AN_ext'] - _df['AN_ext_dry']
-        _df['ALWC_SS_ext'] = _df['SS_ext'] - _df['SS_ext_dry']
-        _df['ALWC_ext'] = _df['total_ext'] - _df['total_ext_dry']
-
-        _df['fRH_IMPR'] = _df['total_ext'] / _df['total_ext_dry']
-
-        return _df['AS_ext_dry':]
-
-    def modified(self, _df):
-        _frh, _frhss, _frhs, _frhl = self.frh(_df['RH'])
-
-        _df['AS_ext_dry'] = 3 * 1 * _df['AS']
-        _df['AN_ext_dry'] = 3 * 1 * _df['AN']
-        _df['OM_ext_dry'] = 4 * _df['OM']
-        _df['Soil_ext_dry'] = 1 * _df['Soil']
-        _df['SS_ext_dry'] = 1.7 * 1 * _df['SS']
-        _df['EC_ext_dry'] = 10 * _df['EC']
-        _df['total_ext_dry'] = sum(_df['AS_ext_dry':'EC_ext_dry'])
-
-        _df['AS_ext'] = (3 * _frh * _df['AS'])
-        _df['AN_ext'] = (3 * _frh * _df['AN'])
-        _df['OM_ext'] = (4 * _df['OM'])
-        _df['Soil_ext'] = (1 * _df['Soil'])
-        _df['SS_ext'] = (1.7 * _frhss * _df['SS'])
-        _df['EC_ext'] = (10 * _df['EC'])
-        _df['total_ext'] = sum(_df['AS_ext':'EC_ext'])
-
-        _df['ALWC_AS_ext'] = _df['AS_ext'] - _df['AS_ext_dry']
-        _df['ALWC_AN_ext'] = _df['AN_ext'] - _df['AN_ext_dry']
-        _df['ALWC_SS_ext'] = _df['SS_ext'] - _df['SS_ext_dry']
-        _df['ALWC_ext'] = _df['total_ext'] - _df['total_ext_dry']
-
-        _df['fRH_IMPR'] = _df['total_ext'] / _df['total_ext_dry']
-
-        return _df['AS_ext_dry':]
-
-    @staticmethod
-    def gas(_df):
-        _df['ScatteringByGas'] = (11.4 * 293 / (273 + _df['AT']))
-        _df['AbsorptionByGas'] = (0.33 * _df['NO2'])
-        _df['ExtinctionByGas'] = _df['ScatteringByGas'] + _df['AbsorptionByGas']
-        return _df['ScatteringByGas':]
-
-    def process_data(self, reset: bool = False, save_file: Path | str = None,
-                     version: Literal["revised", "modified"] = "revised"):
-        save_file = Path(save_file)
-        if save_file.exists() and not reset:
-            return read_csv(save_file, parse_dates=['Time'], index_col='Time')
-        else:
-            # data_files = ['EPB.csv', 'IMPACT.csv', 'chemical.csv']
-            df = concat([DataReader(file) for file in self.file_paths], axis=1)
-
-            # particle contribution '銨不足不納入計算'
-            improve_input_df = df.loc[df['NH4_status'] != 'Deficiency', ['AS', 'AN', 'OM', 'Soil', 'SS', 'EC', 'RH']]
-
-            df_improve = improve_input_df.dropna().copy().apply(self.revised if version == 'revised' else self.modified,
-                                                                axis=1)
-
-            # gas contribution
-            df_ext_gas = df[['NO2', 'AT']].dropna().copy().apply(self.gas, axis=1)
-
-            _df = concat([df_improve, df_ext_gas], axis=1).reindex(df.index.copy())
-            _df.to_csv(save_file)
-
-            return _df

AeroViz/process/script/Others.py

@@ -1,65 +0,0 @@
-from pathlib import Path
-
-import numpy as np
-from pandas import read_csv, concat, DataFrame
-
-from AeroViz.process.core import DataProc
-from AeroViz.tools.datareader import DataReader
-
-
-class OthersProc(DataProc):
-    """
-    A class for process impact data.
-
-    Parameters:
-    -----------
-    reset : bool, optional
-        If True, resets the process. Default is False.
-    filename : str, optional
-        The name of the file to process. Default is None.
-
-    Methods:
-    --------
-    process_data():
-        Process data and save the result.
-
-    Attributes:
-    -----------
-    DEFAULT_PATH : Path
-        The default path for data files.
-
-    Examples:
-    ---------
-    >>> df = OthersProc().process_data(reset=True, filename=None)
-
-    """
-
-    def __init__(self, file_paths: Path | list[Path | str] = None):
-        super().__init__()
-        self.file_paths = [Path(fp) for fp in file_paths]
-
-    def process_data(self, reset: bool = False, save_file: Path | str = None) -> DataFrame:
-        save_file = Path(save_file)
-        if save_file.exists() and not reset:
-            return read_csv(save_file, parse_dates=['Time'], index_col='Time')
-        else:
-            df = concat([DataReader(file) for file in self.file_paths], axis=1)
-
-            results = DataFrame(index=df.index)
-
-            results['PG'] = df[
-                ['Scattering', 'Absorption', 'ScatteringByGas', 'AbsorptionByGas']].dropna().copy().apply(np.sum,
-                                                                                                          axis=1)
-            results['MAC'] = df['Absorption'] / df['T_EC']
-            results['Ox'] = df['NO2'] + df['O3']
-            results['N2O5_tracer'] = df['NO2'] * df['O3']
-            results['Vis_cal'] = 1096 / df['Extinction']
-            # results['fRH_Mix'] = df['Bext'] / df['Extinction']
-            # results['fRH_PNSD'] = df['Bext_internal'] / df['Bext_dry']
-            results['fRH_IMPR'] = df['total_ext'] / df['total_ext_dry']
-            results['OCEC_ratio'] = df['O_OC'] / df['O_EC']
-            results['PM1/PM25'] = np.where(df['PM1'] / df['PM2.5'] < 1, df['PM1'] / df['PM2.5'], np.nan)
-            # results['MEE_PNSD'] = df['Bext_internal'] / df['PM25']
-            # results['MEE_dry_PNSD'] = df['Bext_dry'] / df['PM25']
-
-            return results

AeroViz/process/script/PSD.py

@@ -1,103 +0,0 @@
-from pathlib import Path
-
-from pandas import concat, read_csv, DataFrame
-
-from AeroViz.process.core import DataProc
-from AeroViz.process.core.SizeDist import SizeDist
-from AeroViz.process.script.AbstractDistCalc import DistributionCalculator
-
-
-class ParticleSizeDistProc(DataProc):
-	"""
-    A class for process particle size distribution (PSD) data.
-
-    Parameters
-    ----------
-    filename : str, optional
-        The name of the PSD data file.
-        Defaults to 'PNSD_dNdlogdp.csv' in the default path.
-
-    Attributes
-    ----------
-    file_path : Path
-        The directory path where the PSD data file is located.
-
-    psd : SizeDist
-        The SizeDist object.
-
-    Methods
-    -------
-    process_data(filename='PSD.csv')
-        Process and save overall PSD properties.
-
-    Examples
-    --------
-    Example 1: Use default path and filename
-    >>> psd_data = ParticleSizeDistProc(filename='PNSD_dNdlogdp.csv').process_data(reset=True)
-    """
-
-	def __init__(self, file_path: Path | str = None):
-		super().__init__()
-		self.file_path = Path(file_path)
-
-		self.psd = SizeDist(read_csv(file_path, parse_dates=['Time'], index_col='Time'))
-
-	def process_data(self, reset: bool = False, save_file: Path | str = None) -> DataFrame:
-		save_file = Path(save_file)
-		if save_file.exists() and not reset:
-			return read_csv(save_file, parse_dates=['Time'], index_col='Time')
-
-		number = DistributionCalculator('number', self.psd).useApply()
-		surface = DistributionCalculator('surface', self.psd).useApply()
-		volume = DistributionCalculator('volume', self.psd).useApply()
-
-		surface.to_csv(save_file.parent / 'PSSD_dSdlogdp.csv')
-		volume.to_csv(save_file.parent / 'PVSD_dVdlogdp.csv')
-
-		result_df = concat(
-			[DistributionCalculator('property', SizeDist(data=number, weighting='n')).useApply(),
-			 DistributionCalculator('property', SizeDist(data=surface, weighting='s')).useApply(),
-			 DistributionCalculator('property', SizeDist(data=volume, weighting='v')).useApply()
-			 ], axis=1)
-
-		result_df.to_csv(save_file)
-		return result_df
-
-
-class ExtinctionDistProc(DataProc):
-
-	def __init__(self, file_path: Path | str = 'PNSD_dNdlogdp.csv', file_path_chem: Path | str = 'chemical.csv'):
-		super().__init__()
-		self.file_path = Path(file_path)
-		self.file_path_chem = Path(file_path_chem)
-
-		self.psd = SizeDist(read_csv(file_path, parse_dates=['Time'], index_col='Time'))
-		self.RI = read_csv(file_path_chem, parse_dates=['Time'], index_col='Time')[['n_dry', 'n_amb', 'k_dry', 'k_amb',
-																					'AS_volume_ratio',
-																					'AN_volume_ratio',
-																					'OM_volume_ratio',
-																					'Soil_volume_ratio',
-																					'SS_volume_ratio',
-																					'EC_volume_ratio',
-																					'ALWC_volume_ratio']]
-
-	def process_data(self, reset: bool = False, save_file: Path | str = 'PESD.csv'):
-		save_file = Path(save_file)
-		if save_file.exists() and not reset:
-			return read_csv(save_file, parse_dates=['Time']).set_index('Time')
-
-		ext_internal = DistributionCalculator('extinction', self.psd, self.RI, method='internal',
-											  result_type='extinction').useApply()
-		ext_external = DistributionCalculator('extinction', self.psd, self.RI, method='external',
-											  result_type='extinction').useApply()
-
-		ext_internal.to_csv(save_file.parent / 'PESD_dextdlogdp_internal.csv')
-		ext_external.to_csv(save_file.parent / 'PESD_dextdlogdp_external.csv')
-
-		result_df = concat([
-			DistributionCalculator('property', SizeDist(data=ext_internal, weighting='ext_in')).useApply(),
-			DistributionCalculator('property', SizeDist(data=ext_internal, weighting='ext_ex')).useApply(),
-		], axis=1)
-
-		result_df.to_csv(save_file)
-		return result_df

AeroViz/process/script/PSD_dry.py

@@ -1,93 +0,0 @@
-from pathlib import Path
-
-import numpy as np
-from pandas import DataFrame, read_csv, concat
-
-from AeroViz.process.core import DataProc
-from AeroViz.process.core.SizeDist import SizeDist
-
-
-class DryPSDProc(DataProc):
-    """
-    A class for process impact data.
-
-    Parameters
-    ----------
-    reset : bool, optional
-        If True, resets the process. Default is False.
-    filename : str, optional
-        The name of the file to process. Default is None.
-
-    Methods
-    -------
-    process_data():
-        Process data and save the result.
-
-    Attributes
-    ----------
-    DEFAULT_PATH : Path
-        The default path for data files.
-
-
-    Examples
-    --------
-    >>> df = DryPSDProc(reset=True, filename='PNSD_dNdlogdp_dry.csv').process_data()
-    """
-
-    def __init__(self, file_path: Path | str = 'PNSD_dNdlogdp.csv', file_path_chem: Path | str = 'chemical.csv'):
-        super().__init__()
-        self.file_path = Path(file_path)
-        self.file_path_chem = Path(file_path_chem)
-
-        self.psd = SizeDist(read_csv(file_path, parse_dates=['Time'], index_col='Time'))
-        self.RI = read_csv(file_path_chem, parse_dates=['Time'], index_col='Time')[['n_dry', 'n_amb', 'k_dry', 'k_amb',
-                                                                                    'AS_volume_ratio',
-                                                                                    'AN_volume_ratio',
-                                                                                    'OM_volume_ratio',
-                                                                                    'Soil_volume_ratio',
-                                                                                    'SS_volume_ratio',
-                                                                                    'EC_volume_ratio',
-                                                                                    'ALWC_volume_ratio']]
-
-    def process_data(self, reset: bool = False, save_filename: Path | str = None) -> DataFrame:
-        save_filename = Path(save_filename)
-        if save_filename.exists() and not reset:
-            return read_csv(save_filename, parse_dates=['Time']).set_index('Time')
-        _df = concat([self.psd, self.RI], axis=1)
-        _df.to_csv(save_filename)
-        return _df
-
-
-def dry_PNSD_process(dist, dp, **kwargs):
-    ndp = np.array(dist[:np.size(dp)])
-    gRH = resolved_gRH(dp, dist['gRH'], uniform=True)
-
-    dry_dp = dp / gRH
-    belong_which_ibin = np.digitize(dry_dp, dp) - 1
-
-    result = {}
-    for i, (ibin, dn) in enumerate(zip(belong_which_ibin, ndp)):
-        if dp[ibin] not in result:
-            result[dp[ibin]] = []
-        result[dp[ibin]].append(ndp[i])
-
-    dry_ndp = []
-    for key, val in result.items():
-        dry_ndp.append(sum(val) / len(val))
-
-    return np.array(dry_ndp)
-
-
-def resolved_gRH(dp, gRH=1.31, uniform=True):
-    if uniform:
-        return np.array([gRH] * dp.size)
-
-    else:
-        lognorm_dist = lambda x, geoMean, geoStd: (gRH / (np.log10(geoStd) * np.sqrt(2 * np.pi))) * np.exp(
-            -(x - np.log10(geoMean)) ** 2 / (2 * np.log10(geoStd) ** 2))
-        abc = lognorm_dist(np.log10(dp), 200, 2.0)
-        return np.where(abc < 1, 1, abc)
-
-
-if __name__ == '__main__':
-    pass

AeroViz/process/script/__init__.py

@@ -1,5 +0,0 @@
-from .Chemical import ChemicalProc
-from .IMPACT import ImpactProc
-from .IMPROVE import ImproveProc
-from .Others import OthersProc
-from .PSD import ParticleSizeDistProc, ExtinctionDistProc