AeroViz 0.1.3__py3-none-any.whl → 0.1.4__py3-none-any.whl

AeroViz-0.1.4.dist-info/RECORD

@@ -0,0 +1,112 @@
+AeroViz/__init__.py,sha256=A5W6SR71uY_eW44Sh-Yk6blJQ_G1aHrkSzeP2YTPQc4,371
+AeroViz/data/DEFAULT_DATA.csv,sha256=eeeyeh8vSLKkE5tAF0TYnUNOyQIH98VA41bJaAP204Y,2248526
+AeroViz/data/DEFAULT_PNSD_DATA.csv,sha256=imLvLA80oYwo_jzXZtlQn5hZ76d47HUIlK2jp0tZPrg,2636511
+AeroViz/dataProcess/__init__.py,sha256=D3rTVUiGfs_daGuaotVtbijOgLAp6HaRWchj-zoEnHw,828
+AeroViz/dataProcess/Chemistry/__init__.py,sha256=vWlQVQqyX2R_qUQ8cLrIYWmdqLBaxbYwK3PaO7zRqIw,1986
+AeroViz/dataProcess/Chemistry/_calculate.py,sha256=q7ojTFPok0vg8k_1PMECNdP5CPanR9NWQ4Rx5iTcHew,599
+AeroViz/dataProcess/Chemistry/_isoropia.py,sha256=PJPdqErssUPnrksR6BocmgQch1DB7lULyge7diOq9uc,3304
+AeroViz/dataProcess/Chemistry/_mass_volume.py,sha256=oqFjRDYZ1om0_kG4g8XNi2QpUPcHsukpnWKydQoNEVg,5394
+AeroViz/dataProcess/Chemistry/_ocec.py,sha256=75r2EroxVrcp2tSHgNC9TAif_MeJeRlrFyo6DvCLemo,5957
+AeroViz/dataProcess/Chemistry/_partition.py,sha256=tKhb6BJns46UiUlEq6Zq7ahYnvUJ_whY3tWE54C3bqU,1023
+AeroViz/dataProcess/Chemistry/_teom.py,sha256=IiM-TrifWpQLTbKllG-4k4c3mvQulfcmjswWu6muCXA,486
+AeroViz/dataProcess/Chemistry/isrpia.cnf,sha256=iWXTqsOZFmNrJxAI9nYuilZ9h6ru1icdPFVim7YKc_k,566
+AeroViz/dataProcess/Optical/Angstrom_exponent.py,sha256=KSGM_hf2oUPsLoLuhYSRDK5zdekwce9KgqhzRnf2Y4g,706
+AeroViz/dataProcess/Optical/_IMPROVE.py,sha256=_u0uyKC9ZQpB4uLMx3x8Dm5hSUMcdGUifBcC2bGaM9c,1675
+AeroViz/dataProcess/Optical/__init__.py,sha256=h_IYhaum0DhK7Qw6lWFvjAEv7uAqLVVPEwU0V4hJnX4,1371
+AeroViz/dataProcess/Optical/_absorption.py,sha256=pYHIoj_R23wKvxR4acUqqU9JXKY40wgsFvZH78rETrY,1234
+AeroViz/dataProcess/Optical/_extinction.py,sha256=anUemSXmLvJi2CSxqTYgQCgUW4pSUqgZ0q0xZ4zs0zY,1189
+AeroViz/dataProcess/Optical/_mie.py,sha256=K_I8NqRStCdExRVxq2rX5hQdvYmfKuu9JCkNtL1OzfE,352
+AeroViz/dataProcess/Optical/_mie_sd.py,sha256=cbWSorfG-l92csaXsZJ7005f8ij2H_VZopIhNPCwx70,4812
+AeroViz/dataProcess/Optical/_scattering.py,sha256=w5xXXOfafV9a_dXY60YqTyzFDZppjBJeLcKl5PeCQe4,961
+AeroViz/dataProcess/SizeDistr/__init__.py,sha256=PJVAdb-v1q8h38JXIkpVcr6q4mg38V26sk4xD_2Mjlc,2186
+AeroViz/dataProcess/SizeDistr/__merge.py,sha256=A-ZmduUA9a0nVM8szhG1AWhsfUPu2UhGIaqFLcmg7Lk,9422
+AeroViz/dataProcess/SizeDistr/_merge.py,sha256=CyZwUc2M4lL-OE74e8mXq41nC6X0LJTAV8TwpO3fS2o,9176
+AeroViz/dataProcess/SizeDistr/_merge_v1.py,sha256=6Anb8DszoatK66tc9ccA6ZApbqtL7pZRPR5TKHj-Tbs,9526
+AeroViz/dataProcess/SizeDistr/_merge_v2.py,sha256=8OzUKw7hTg-yuQBipuFKgBS_7c7zbApN_BNr00G8q9c,9046
+AeroViz/dataProcess/SizeDistr/_merge_v3.py,sha256=HN2ARFmeWOawOWRPPv_pHEGBBZNgXVbH4dDTxcN7rdY,18749
+AeroViz/dataProcess/SizeDistr/_merge_v4.py,sha256=b8RVAievGIOLrmJHJXRsKXQ1tkMkm6rx43S7XAfeXE4,16228
+AeroViz/dataProcess/SizeDistr/_size_distr.py,sha256=o5fTwLH7j9j5129e3uciSJrOR8AGCN3tkAyGffoPENg,3127
+AeroViz/dataProcess/VOC/__init__.py,sha256=3W5G28UowKVfwsj-tFlxdsCEzXSMxrlrBjpiqUfneaA,290
+AeroViz/dataProcess/VOC/_potential_par.py,sha256=H_9EoDpbz6dZ193FuJEsoaOQ7jCDz9LqPqhEce_nOFQ,3065
+AeroViz/dataProcess/VOC/voc_par.json,sha256=cJF0ynKCMccIruHp7YJ19fZmXZV6ASlLszCJCqduDJQ,7077
+AeroViz/dataProcess/core/__init__.py,sha256=FhGwqlP8ZwXZQS0LjH0YOz2g1Q2yJniUUhRZ_jDOQ6k,2707
+AeroViz/plot/__init__.py,sha256=BxoarxAeDHNdMHizS1bCTwV8NEuUCgeDB75GbSEsNWs,394
+AeroViz/plot/bar.py,sha256=cSPKJLWGifCzICCZVPcbZtrDrea_AStlwSli9-b-lsA,4246
+AeroViz/plot/box.py,sha256=8zDd8jYndSv5f0P1lFeX0m7VWYh9okYgJTClh2Lj7Dg,2115
+AeroViz/plot/pie.py,sha256=PVOqc_tAMlfjojTggP7HPL9Nk3ImigETo_EUFD6K2BY,7421
+AeroViz/plot/regression.py,sha256=mTCNrCAYvH7QjIb_5WlltCR1Kqwer1anpGEsxmO2S3o,7292
+AeroViz/plot/scatter.py,sha256=2Q_Um7Ai-pym2nUiMk6tY_zEXwbKE-BpkgrDkMyNbYA,6293
+AeroViz/plot/violin.py,sha256=pU2Z2yTWocEtImmCAmbtn0WvXtUOrnCGOdDOrLxjooU,2689
+AeroViz/plot/distribution/__init__.py,sha256=nhbIegWczkuEfWsE7-2jfF0dnpmPDzJJzjq8Fuh6q5k,28
+AeroViz/plot/distribution/distribution.py,sha256=sAjqtqKavFwQqI8PGPFnpvZFSU-w2UKjcTTC5L91f4E,20595
+AeroViz/plot/meteorology/__init__.py,sha256=hhGfQE3IUzS3Eaju_nO7LomPPHJnd-zAAZZweXOXs2M,27
+AeroViz/plot/meteorology/meteorology.py,sha256=h5KiaReatt238FXZyGjP9CAMkzTF1aVjOaFAGq3CKRY,11197
+AeroViz/plot/optical/__init__.py,sha256=5S2WPq8NouPYgc46c1XBveSL8Cxs3w8wNYN7vWIpB-s,23
+AeroViz/plot/optical/optical.py,sha256=xtuJuaonyM4OUS_BlkxKudLZKEteCS1u_-Svj3lQPco,13728
+AeroViz/plot/templates/__init__.py,sha256=VXp73SqrSF1WR9WYhHpHtbhvFcb1CR4JeV1tPQvCEEI,227
+AeroViz/plot/templates/ammonium_rich.py,sha256=a46rfwhSZ3LIj71LxdOnklAFyy6Xe_VQ-kh1komnZ5g,1225
+AeroViz/plot/templates/contour.py,sha256=fAM6Wn_EmENQM1IW6KJpwmLva6D9P2dEAaK6cNQBtdQ,1678
+AeroViz/plot/templates/corr_matrix.py,sha256=VeiwBwmg_v6LYslJvhSqweiLBheR3vGfPYxrtHoqJRI,3507
+AeroViz/plot/templates/diurnal_pattern.py,sha256=Vp8nDzSe95ppZBCDkFXnpUEh1cb-19yDS3H0qftoPaw,1434
+AeroViz/plot/templates/koschmieder.py,sha256=XmsYn8x9lu2BgN1FXKHb0quZdRw7Etv1dT5LBjBB0ok,3199
+AeroViz/plot/templates/metal_heatmap.py,sha256=ngqJFiOPCoExPKv-lHvgB1XCSr9fF_jSKmx4VwjUUXQ,4961
+AeroViz/plot/timeseries/__init__.py,sha256=LQ_aE32pvp5ouC--BRa0NwSJ48DYO3rFaOt4AKGmn_s,50
+AeroViz/plot/timeseries/template.py,sha256=E9dxaH_Y5QxnkIHK4VH59IVShnf4i1-TT_QE4NqVqrQ,1460
+AeroViz/plot/timeseries/timeseries.py,sha256=EitNnjleHj3Qu1PBHClWNi9kCgJ8tULGZSicS3xhoZw,14626
+AeroViz/plot/utils/__init__.py,sha256=sxp3nfUnyVIIXYd-F874gEHkhQiDnMPDE1nJ_Z98270,104
+AeroViz/plot/utils/_color.py,sha256=Vhjo-v4UWFME1JbbBR_Ro8FRj6dHB2O5Sud_9lpdu_U,2797
+AeroViz/plot/utils/_unit.py,sha256=Abum0WDnq08DzkSoxkKvBc0mYKo8izIDmgduPo9Q27A,1671
+AeroViz/plot/utils/fRH.json,sha256=t-2ux4TLOYAB-4jJ72LSM4jv1jk9XkaxKYNMDepMHIg,6038
+AeroViz/plot/utils/plt_utils.py,sha256=7Au3r2-7AZQmzrO2OfcyTFomJRLHgu1Npb7wxQqUvzY,3438
+AeroViz/plot/utils/sklearn_utils.py,sha256=hKFfkVkYLRxkIDKvO9COHXwhjD_UWqQigdT3mDToni4,2098
+AeroViz/plot/utils/units.json,sha256=JKcqvLA6dkc8REV_NhX85Jl9LA4iAZxzw9RJp6JTla0,2965
+AeroViz/process/__init__.py,sha256=vqXO0vhowU9yNW5tnVYNg3ceBxLt96PgBNkZ2KUxKXw,1309
+AeroViz/process/core/DataProc.py,sha256=6MBDcHyQmVU16Z0oDo0GItvhU5f7fWLXFMGqUyoW1eQ,419
+AeroViz/process/core/SizeDist.py,sha256=V6CpxTRfWjsZ8GwIEOqyKJ8_qcz7dO0ZmfHQvR6SgW4,2097
+AeroViz/process/core/__init__.py,sha256=WmVSFDx97urX0WPA8dtRVoQF04IDaU-IgxHK979-kM8,98
+AeroViz/process/method/PyMieScatt_update.py,sha256=g3vlzATjzYSYZd3LwmwxEmdkpo4cHJ3KY4rePY4jwZk,21065
+AeroViz/process/method/__init__.py,sha256=_yxoNAbbI82ORhrVsi-6x8vPnWQNHavjZUXy-ioZ2wg,123
+AeroViz/process/method/mie_theory.py,sha256=EC3SsUaq5G3SdQyPMQxcosbEokqAgJw-4Gz6moJiOsk,8666
+AeroViz/process/method/prop.py,sha256=7f18SOkkLgDiY_g1c1We8g-3WDJJIFcxkfR4F99jUgc,1845
+AeroViz/process/script/AbstractDistCalc.py,sha256=wXF7OA634m_i-Y79lU4THj0c0yZx5NUtLNvORH1FWZg,5351
+AeroViz/process/script/Chemical.py,sha256=jjBIdILe0fO5Rqidb4pg3CP6qmMLqFvst-ZSz8mbb7Q,6413
+AeroViz/process/script/IMPACT.py,sha256=Zz1GRQLJopNQGw5qUTdWl7Zxy_CBnFBsL1B-CCemaDQ,1520
+AeroViz/process/script/IMPROVE.py,sha256=5FNl33S_785J7CsSue_CvZP58TWevpwGr2ePQM6qScw,5709
+AeroViz/process/script/Others.py,sha256=TOtSKVZ9Z-zX1LFfE4p36H0hMjfmyyQw_Xaeq8HgMBk,2325
+AeroViz/process/script/PSD.py,sha256=unSvyHp1RqgQecY-Lab-wRFGP2u45blSO4NCNTaxcRQ,3721
+AeroViz/process/script/PSD_dry.py,sha256=P9rSQVA78jHXrY2bDiC1iQvG0w9Js8QiBNfooPOXzfY,3269
+AeroViz/process/script/__init__.py,sha256=QExUHHK79oiwZXvk8u00P8rLwWhGOc1fpZY2Fx8Vkxg,188
+AeroViz/process/script/retrieve_RI.py,sha256=9OQJU1w6VzVqKz00e9GdWswRAnaETe_1meAHp5N_s54,2711
+AeroViz/rawDataReader/__init__.py,sha256=ytqVoXXYk1xeuyqdo2VJkC3JbIA5RmBz20uzgBnMKW0,1881
+AeroViz/rawDataReader/config/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+AeroViz/rawDataReader/config/supported_instruments.py,sha256=3h25eUGysciCzvDeqReDr7YtYLC8fKAcU903FGiGs9g,4728
+AeroViz/rawDataReader/core/__init__.py,sha256=Ne2_eAzasAcwrWmG-NU35cJy1QMkLU1-iWEPLGZdcqU,10304
+AeroViz/rawDataReader/script/AE33.py,sha256=uwadLYg3RsxW1Yhdtb83qy5lWTTBpliYgfYv4GXV-Js,1110
+AeroViz/rawDataReader/script/AE43.py,sha256=-gdpPJg6EHnU7mKNNCsDe-aytbDlR_qyRWo_kgg0Ss8,1058
+AeroViz/rawDataReader/script/APS_3321.py,sha256=6fYwMPLG0w2RduNqIelwQ3Q2rsiJUZKQUetkyaDg0BU,1696
+AeroViz/rawDataReader/script/Aurora.py,sha256=I2xGA-HwE-mlSqZZQ_GHRM3OxE27L3uvbEXfK0p58ig,1238
+AeroViz/rawDataReader/script/BC1054.py,sha256=XD3CZjg9rd4JzQEF4MCNqyRP1931JhMi-pGkgE7B-Ig,1488
+AeroViz/rawDataReader/script/EPA_vertical.py,sha256=AUXbRlnO2T-CCcBdvyMS4TxYGSp3RI679aZH90thh6k,1671
+AeroViz/rawDataReader/script/GRIMM.py,sha256=wvWDksAewCqYNNMR3IuxilGrzooQWzzWFK2b1CbxWjo,950
+AeroViz/rawDataReader/script/IGAC.py,sha256=3F0gLz0-z5_Sj6JQmO4wklrjfRHkX-EnHfBiioBpcvI,2950
+AeroViz/rawDataReader/script/MA350.py,sha256=yjg4WYGpPiZbtlF1gKfrWF1WS-nbGPTpbKZ3ItixOPs,1289
+AeroViz/rawDataReader/script/Minion.py,sha256=x29UQvjHZbWQigSJuNuGGOyQQxP9O4kPSJ2S1ixR5XM,3942
+AeroViz/rawDataReader/script/NEPH.py,sha256=hTWzZmoBeTm7qKkVFYuyC07iDYOmAexamcA7gfdd2Ak,2806
+AeroViz/rawDataReader/script/SMPS_TH.py,sha256=-ME3MxXv2jWBSwsO3uDcBe28sWcrJ9eJpqcS5VMva2A,1470
+AeroViz/rawDataReader/script/SMPS_aim11.py,sha256=zYE482FQICuFVDo0KtpKL6lBOfdXaxObL7ep2I_ufeE,1634
+AeroViz/rawDataReader/script/SMPS_genr.py,sha256=vFNJXsXpUmZhfuobaFsfypVU98iA4bWE3Tk32LTIh1Q,1611
+AeroViz/rawDataReader/script/Sunset_OCEC.py,sha256=O-o8It8XtVDc50zVT043OqdG_eHkPS229-zEt5WR9cs,1953
+AeroViz/rawDataReader/script/TEOM.py,sha256=pFxks2To9y7n_OlpbrGPcRJEmzK-Y_oVoccHvQqQMg8,1798
+AeroViz/rawDataReader/script/Table.py,sha256=piLBiBhLeqR5Wmd9LLma6ubEbGa--3rdrvRCrQa8NPM,854
+AeroViz/rawDataReader/script/VOC.py,sha256=tl6rzqfVU5ofa2R3Ax-_1i37Lno9rUuMhlylPObYFFY,716
+AeroViz/rawDataReader/script/__init__.py,sha256=UndjtWjvJHMma7WtArL5sF9UHiQsanx5sT2otQDqzng,270
+AeroViz/tools/__init__.py,sha256=tPUmCI9Fi1LgE-QSVnzan5jRdYQBIENH9PRDXB_DhTw,109
+AeroViz/tools/database.py,sha256=05VzjJyhlRrhsZdhfFQ__7CxGm4MdFekLjz3_Is5h9U,3430
+AeroViz/tools/dataclassifier.py,sha256=_wpv0PlZ5EGkcNqHxfFtdEsYvHP5FVE8sMZXikhm_YE,4492
+AeroViz/tools/dataprinter.py,sha256=Jq2Yztpa9YCOeLDVTrRs7PhSdNIPhEAexVj1YSuJ7hY,2249
+AeroViz/tools/datareader.py,sha256=iTQ0U8hdNMjCdbiH7EiKW10UEoxzxXRHc4s5_1IikJo,1933
+AeroViz-0.1.4.dist-info/LICENSE,sha256=E-679GpGGkp3irmtuJXiT7R4cNUA4cmsH6Q7QUgPf5U,1069
+AeroViz-0.1.4.dist-info/METADATA,sha256=RQchK9CdQftok3LH-Cc3gfbzC5g-R-7UWzVQl7VagMA,5323
+AeroViz-0.1.4.dist-info/WHEEL,sha256=HiCZjzuy6Dw0hdX5R3LCFPDmFS4BWl8H-8W39XfmgX4,91
+AeroViz-0.1.4.dist-info/top_level.txt,sha256=BYsmTst_o4FZOKRP1XIvIMlN6mMTTXNfnSToL2_nVbQ,8
+AeroViz-0.1.4.dist-info/RECORD,,

AeroViz/plot/improve/__init__.py

@@ -1 +0,0 @@
-from .improve import *

AeroViz/plot/improve/improve.py

@@ -1,240 +0,0 @@
-import matplotlib.pyplot as plt
-import numpy as np
-import pandas as pd
-from matplotlib.pyplot import Figure, Axes
-from pandas import DataFrame, read_json
-from scipy.optimize import curve_fit
-from pathlib import Path
-
-from AeroViz import plot
-from AeroViz.plot.utils import *
-from AeroViz.tools import DataBase, DataReader, DataClassifier
-
-# TODO: this file has to be reorganized
-
-__all__ = ['chemical_enhancement',
-		   'ammonium_rich',
-		   'pie_IMPROVE',
-		   'MLR_IMPROVE',
-		   'fRH_plot',
-		   ]
-
-
-@set_figure
-def chemical_enhancement(data_set: DataFrame = None,
-						 data_std: DataFrame = None,
-						 ax: Axes | None = None,
-						 **kwargs
-						 ) -> tuple[Figure, Axes]:
-	fig, ax = plt.subplots() if ax is None else (ax.get_figure(), ax)
-
-	ser_grp_sta, ser_grp_sta_std = DataClassifier(DataBase('/Users/chanchihyu/NTU/2020能見度計畫/data/All_data.csv'),
-												  by='State')
-	species = ['AS', 'AN', 'POC', 'SOC', 'Soil', 'SS', 'EC', 'ALWC']
-	data_set, data_std = ser_grp_sta.loc[:, species], ser_grp_sta_std.loc[:, species]
-
-	width = 0.20
-	block = width / 4
-
-	x = np.array([1, 2, 3, 4, 5, 6, 7])
-	for i, state in enumerate(['Clean', 'Transition', 'Event']):
-		val = np.array(data_set.iloc[i, :-1])
-		std = (0,) * 7, np.array(data_std.iloc[i, :-1])
-
-		plt.bar(x + (i + 1) * (width + block), val, yerr=std, width=width, color=Color.colors3[:-1],
-				alpha=0.6 + (0.2 * i),
-				edgecolor=None, capsize=None, label=state)
-
-	ax.set(xlabel=r'$\bf Chemical\ species$',
-		   ylabel=r'$\bf Mass\ concentration\ ({\mu}g/m^3)$',
-		   xticks=x + 2 * (width + block),
-		   xticklabels=species,
-		   ylim=(0, 25),
-		   title=r'$\bf Chemical\ enhancement$')
-
-	ax.vlines(8, 0, 25, linestyles='--', colors='k')
-
-	ax2 = ax.twinx()
-	for i, state in enumerate(['Clean', 'Transition', 'Event']):
-		val = np.array(data_set.iloc[i, -1])
-		std = np.array([[0], [data_std.iloc[i, -1]]])
-		plt.bar(8 + (i + 1) * (width + block), val, yerr=std, width=width, color='#96c8e6',
-				alpha=0.6 + (0.2 * i), edgecolor=None, capsize=None, label=state)
-
-	ax2.set(ylabel=r'$\bf Mass\ concentration\ ({\mu}g/m^3)$',
-			ylim=(0, 100),
-			xticks=x + 2 * (width + block),
-			xticklabels=species
-			)
-
-	a = (np.array(data_set.loc['Event']) + np.array(data_set.loc['Transition'])) / 2
-	b = (np.array(data_set.loc['Transition']) + np.array(data_set.loc['Clean'])) / 2
-	c = np.array(data_set.loc['Event']) / np.array(data_set.loc['Transition'])
-	d = np.array(data_set.loc['Transition']) / np.array(data_set.loc['Clean'])
-
-	for i, (posa, posb, vala, valb) in enumerate(zip(a, b, c, d)):
-		if i < 7:
-			ax.text(i + 1.5, posa, '{:.2f}'.format(vala), fontsize=6, weight='bold', zorder=1)
-			ax.text(i + 1.25, posb, '{:.2f}'.format(valb), fontsize=6, weight='bold', zorder=1)
-		else:
-			ax2.text(i + 1.5, posa, '{:.2f}'.format(vala), fontsize=6, weight='bold', zorder=1)
-			ax2.text(i + 1.25, posb, '{:.2f}'.format(valb), fontsize=6, weight='bold', zorder=1)
-
-	plt.show()
-
-	return fig, ax
-
-
-@set_figure
-def ammonium_rich(df: DataFrame,
-				  **kwargs
-				  ) -> tuple[Figure, Axes]:
-	df = df[['NH4+', 'SO42-', 'NO3-', 'PM25']].dropna().copy().div([18, 96, 62, 1])
-	df['required_ammonium'] = df['NO3-'] + 2 * df['SO42-']
-
-	fig, ax = plt.subplots()
-
-	scatter = ax.scatter(df['required_ammonium'].to_numpy(), df['NH4+'].to_numpy(), c=df['PM25'].to_numpy(),
-						 vmin=0, vmax=70, cmap='jet', marker='o', s=10, alpha=1)
-
-	ax.axline((0, 0), slope=1., color='k', lw=2, ls='--', alpha=0.5, label='1:1')
-	plt.text(0.97, 0.97, r'$\bf 1:1\ Line$', color='k', ha='right', va='top', transform=ax.transAxes)
-
-	ax.set(xlim=(0, 1.2),
-		   ylim=(0, 1.2),
-		   xlabel=r'$\bf NO_{3}^{-}\ +\ 2\ \times\ SO_{4}^{2-}\ (mole\ m^{-3})$',
-		   ylabel=r'$\bf NH_{4}^{+}\ (mole\ m^{-3})$',
-		   title=kwargs.get('title', ''))
-
-	color_bar = plt.colorbar(scatter, label=Unit('PM25'), extend='both')
-
-	# fig.savefig(f'Ammonium_rich_{title}')
-	plt.show()
-
-	return fig, ax
-
-
-def pie_IMPROVE():
-	Species1 = ['AS_ext_dry', 'AN_ext_dry', 'OM_ext_dry', 'Soil_ext_dry', 'SS_ext_dry', 'EC_ext_dry']
-	Species2 = ['AS_ext_dry', 'AN_ext_dry', 'OM_ext_dry', 'Soil_ext_dry', 'SS_ext_dry', 'EC_ext_dry', 'ALWC_ext']
-	Species3 = ['AS_ext', 'AN_ext', 'OM_ext', 'Soil_ext', 'SS_ext', 'EC_ext']
-
-	ser_grp_sta, _ = DataClassifier(DataBase(), by='State')
-
-	ext_dry_dict = ser_grp_sta.loc[:, Species1]
-	ext_amb_dict = ser_grp_sta.loc[:, Species2]
-	ext_mix_dict = ser_grp_sta.loc[:, Species3]
-
-	plot.donuts(data_set=ext_dry_dict, labels=['AS', 'AN', 'OM', 'Soil', 'SS', 'BC'], unit='Extinction')
-	plot.donuts(data_set=ext_mix_dict, labels=['AS', 'AN', 'OM', 'Soil', 'SS', 'BC'], unit='Extinction')
-	plot.donuts(data_set=ext_amb_dict, labels=['AS', 'AN', 'OM', 'Soil', 'SS', 'BC', 'ALWC'],
-				unit='Extinction', colors=Color.colors2)
-
-
-def MLR_IMPROVE(**kwargs):
-	"""
-	Perform multiple linear regression analysis and generate plots based on IMPROVE dataset.
-
-	Parameters
-	----------
-	**kwargs : dict
-		Additional keyword arguments for customization.
-
-	Returns
-	-------
-	None
-
-	Examples
-	--------
-	Example usage of MLR_IMPROVE function:
-
-	>>> MLR_IMPROVE()
-
-	Notes
-	-----
-	This function performs multiple linear regression analysis on the IMPROVE dataset and generates plots for analysis.
-
-	- The function first selects specific species from the dataset and drops NaN values.
-	- It calculates a 'Localized' value based on a multiplier and the sum of selected species.
-	- Data from 'modified_IMPROVE.csv' and 'revised_IMPROVE.csv' are read and concatenated with the dataset.
-	- Statistical analysis is performed using DataClassifier to calculate mean and standard deviation.
-	- Plots are generated using linear_regression for Extinction vs. Revised/Modified/Localized and Pie.donuts for a
-	  pie chart showing the distribution of species based on Extinction.
-
-	"""
-	species = ['Extinction', 'Scattering', 'Absorption',
-			   'total_ext_dry', 'AS_ext_dry', 'AN_ext_dry', 'OM_ext_dry', 'Soil_ext_dry', 'SS_ext_dry', 'EC_ext_dry',
-			   'AS', 'AN', 'POC', 'SOC', 'Soil', 'SS', 'EC', 'OM']
-
-	df = DataBase('/Users/chanchihyu/NTU/2020能見度計畫/data/All_data.csv')[species].dropna().copy()
-
-	# multiple_linear_regression(df, x=['AS', 'AN', 'POC', 'SOC', 'Soil', 'SS'], y='Scattering', add_constant=True)
-	# multiple_linear_regression(df, x=['POC', 'SOC', 'EC'], y='Absorption', add_constant=True)
-	# multiple_linear_regression(df, x=['AS', 'AN', 'POC', 'SOC', 'Soil', 'SS', 'EC'], y='Extinction', add_constant=False)
-
-	multiplier = [2.675, 4.707, 11.6, 7.272, 0, 0.131, 10.638]
-	df['Localized'] = df[['AS', 'AN', 'POC', 'SOC', 'Soil', 'SS', 'EC']].mul(multiplier).sum(axis=1)
-	# TODO: remove name
-	modify_IMPROVE = DataReader('modified_IMPROVE.csv')['total_ext_dry'].rename('Modified')
-	revised_IMPROVE = DataReader('revised_IMPROVE.csv')['total_ext_dry'].rename('Revised')
-
-	df = pd.concat([df, revised_IMPROVE, modify_IMPROVE], axis=1)
-
-	n_df = df[['AS', 'AN', 'POC', 'SOC', 'Soil', 'SS', 'EC']].mul(multiplier)
-	mean, std = DataClassifier(n_df, 'State')
-
-	ser_grp_sta, _ = DataClassifier(DataBase(), by='State')
-	mass_comp = ser_grp_sta.loc[:, ['AS', 'AN', 'POC', 'SOC', 'Soil', 'SS', 'EC']]
-
-	# plot
-	plot.linear_regression(df, x='Extinction', y=['Revised', 'Modified', 'Localized'], xlim=[0, 400], ylim=[0, 400],
-						   regression=True, diagonal=True)
-	plot.donuts(data_set=mass_comp, labels=['AS', 'AN', 'POC', 'SOC', 'Soil', 'SS', 'EC'],
-				unit='PM25', colors=Color.colors3)
-	plot.donuts(mean, labels=['AS', 'AN', 'POC', 'SOC', 'Soil', 'SS', 'EC'], unit='Extinction', colors=Color.colors3)
-
-
-@set_figure
-def fRH_plot(**kwargs) -> tuple[Figure, Axes]:
-	frh = read_json(Path(__file__).parent.parent / 'utils' / 'fRH.json')
-
-	def fitting_func(RH, a, b, c):
-		f = a + b * (RH / 100) ** c
-		return f
-
-	x = frh.index.to_numpy()
-	y = frh['fRHs'].to_numpy()
-
-	result = curve_fit(fitting_func, x, y)
-	params = result[0].tolist()
-	val_fit = fitting_func(x, *params)
-
-	fig, ax = plt.subplots(figsize=(3, 3))
-
-	ax.plot(frh.index, frh['fRH'], 'k-o', ms=2, label='$f(RH)_{original}$')
-	ax.plot(frh.index, frh['fRHs'], 'g-o', ms=2, label='$f(RH)_{small\\ mode}$')
-	ax.plot(frh.index, frh['fRHl'], 'r-o', ms=2, label='$f(RH)_{large\\ mode}$')
-	ax.plot(frh.index, frh['fRHSS'], 'b-o', ms=2, label='$f(RH)_{sea\\ salt}$')
-
-	ax.set(xlim=(0, 100),
-		   ylim=(1, None),
-		   xlabel='$RH (\\%)$',
-		   ylabel='$f(RH)$',
-		   title=f'$Hygroscopic\\ growth\\ factor$'
-		   )
-
-	ax.grid(axis='y', color='gray', linestyle='dashed', linewidth=0.4, alpha=0.4)
-
-	ax.legend()
-
-	plt.show()
-	# fig.savefig('fRH_plot')
-
-	return fig, ax
-
-
-if __name__ == '__main__':
-	# chemical_enhancement()
-	# MLR_IMPROVE()
-	# ammonium_rich()
-	fRH_plot()

AeroViz/plot/optical/aethalometer.py

@@ -1,77 +0,0 @@
-import matplotlib.pyplot as plt
-import numpy as np
-from pandas import date_range
-
-from AeroViz.plot.utils import *
-
-__all__ = ['plot_MA350',
-		   'plot_MA3502',
-		   'plot_day_night']
-
-
-@set_figure(figsize=(15, 5))
-def plot_MA350(df, **kwargs):
-	fig, ax = plt.subplots()
-
-	# ax.scatter(df.index, df['UV BCc'], marker='o', c='purple', alpha=0.5, label='UV BCc')
-	# ax.scatter(df.index, df['Blue BCc'], c='b', alpha=0.5, label='Blue BCc')
-	# ax.scatter(df.index, df['Green BCc'], c='g', alpha=0.5, label='Green BCc')
-	# ax.scatter(df.index, df['Red BCc'], c='r', alpha=0.5, label='Red BCc')
-	mean, std = round(df.mean(), 2), round(df.std(), 2)
-
-	label1 = rf'$MA350-0171\ :\;{mean["MA350_0171 IR BCc"]}\;\pm\;{std["MA350_0171 IR BCc"]}\;(ng/m^3)$'
-	label2 = rf'$MA350-0176\ :\;{mean["MA350_0176 IR BCc"]}\;\pm\;{std["MA350_0176 IR BCc"]}\;(ng/m^3)$'
-	label3 = rf'$BC-1054\ :\;{mean["BC1054 IR BCc"]}\;\pm\;{std["BC1054 IR BCc"]}\;(ng/m^3)$'
-	ax.scatter(df.index, df['MA350_0171 IR BCc'], s=10, ls='-', marker='o', c='#a3b18a', alpha=0.5, label=label1)
-	ax.scatter(df.index, df['MA350_0176 IR BCc'], s=10, ls='-', marker='o', c='#3a5a40', alpha=0.5, label=label2)
-	ax.scatter(df.index, df['BC1054 IR BCc'], s=10, ls='-', marker='o', c='g', alpha=0.5, label=label3)
-	ax.legend(prop={'weight': 'bold'}, loc='upper left')
-
-	st_tm, fn_tm = df.index[0], df.index[-1]
-	tick_time = date_range(st_tm, fn_tm, freq=kwargs.get('freq', '10d'))
-
-	ax.set(xlabel=kwargs.get('xlabel', ''),
-		   ylabel=kwargs.get('ylabel', r'$BC\ (ng/m^3)$'),
-		   xticks=kwargs.get('xticks', tick_time),
-		   xticklabels=kwargs.get('xticklabels', [_tm.strftime("%F") for _tm in tick_time]),
-		   xlim=kwargs.get('xlim', (st_tm, fn_tm)),
-		   ylim=kwargs.get('ylim', (0, None)),
-		   )
-
-
-@set_figure
-def plot_MA3502(df):
-	fig, ax = plt.subplots()
-
-	bins = np.array([375, 470, 528, 625, 880])
-	vals = df.dropna().iloc[:, -5:].values
-
-	ax.boxplot(vals, positions=bins, widths=20,
-			   showfliers=False, showmeans=True, meanline=True, patch_artist=True,
-			   boxprops=dict(facecolor='#f2c872', alpha=.7),
-			   meanprops=dict(color='#000000', ls='none'),
-			   medianprops=dict(ls='-', color='#000000'))
-
-	ax.set(xlim=(355, 900),
-		   ylim=(0, None),
-		   xlabel=r'$\lambda\ (nm)$',
-		   ylabel=r'$Absorption\ (1/Mm)$', )
-
-
-@set_figure(figsize=(6, 5))
-def plot_day_night(df):
-	# Group by hour of day and calculate mean
-	df_grouped = df.groupby(df.index.hour).mean()
-
-	# Create figure and plot
-	fig, ax = plt.subplots()
-	ax.plot(df_grouped.index, df_grouped['MA350_0171 IR BCc'], marker='o', c='k', alpha=0.5, label='MA350-0171')
-	ax.plot(df_grouped.index, df_grouped['MA350_0176 IR BCc'], marker='o', c='r', alpha=0.5, label='MA350-0176')
-	ax.plot(df_grouped.index, df_grouped['BC1054 IR BCc'], marker='o', c='b', alpha=0.5, label='BC-1054')
-
-	ax.set(xlim=(0, 23),
-		   xlabel='Hour of Day',
-		   ylabel=r'$BC\ (ng/m^3)$',
-		   title=f'Diurnal pattern', )
-
-	ax.legend()

AeroViz/plot/templates/event_evolution.py

@@ -1,65 +0,0 @@
-from os.path import join as pth
-
-import matplotlib.pyplot as plt
-import numpy as np
-from pandas import date_range, read_csv
-
-from AeroViz.plot.utils import *
-
-# TODO:
-
-# read csv file
-blh = read_csv(pth('事件分析.csv'), parse_dates=['Time'], index_col='Time')
-
-
-@set_figure(figsize=(12, 5))
-def event_evolution(_df, **kwargs):
-	print(f'Plot : {_df.month[0]}')
-
-	st_tm, fn_tm = _df.index[0], _df.index[-1]
-	tick_time = date_range(st_tm, fn_tm, freq='1d')  # set tick
-
-	# seperate day and night
-	_df_day = _df.between_time('6:00', '17:00').reindex(date_range(st_tm, fn_tm, freq='1h'))
-	_df_night = _df.between_time('18:00', '5:00').reindex(date_range(st_tm, fn_tm, freq='1h'))
-
-	## plot
-	fig, ax = plt.subplots()
-
-	## plot background
-	shade_value, _ = np.meshgrid(_df['PM2.5'], np.arange((1., 2500), 100))
-	ax.pcolormesh(_df.index, np.arange((1., 2500), 100), shade_value, cmap='binary', vmin=0, vmax=300,
-				  shading='auto')
-
-	## plot day and night
-	ld = ax.scatter(_df.index[0:], _df_day['Ext'], s=50, c='#73b9ff', label='Day Ext', marker='o', alpha=.7)
-	ln = ax.scatter(_df.index[0:], _df_night['Ext'], s=50, c='#00238c', label='Night Ext', marker='o', alpha=.7)
-
-	ax2 = ax.twinx()
-	# ld, = ax2.plot(_df_day['VC'],c='#FF9797',label='day 06:00~18:00')
-	# ln, = ax2.plot(_df_night['VC'],c='#FF0000',label='night 18:00~06:00')
-	ld2 = ax2.scatter(_df.index, _df_day['VC'], s=50, c='#FF9797', label='Day VC', marker='o', alpha=.5)
-	ln2 = ax2.scatter(_df.index, _df_night['VC'], s=50, c='#FF0000', label='Night VC', marker='o', alpha=.5)
-
-	# add legend on the first axes
-	ax.legend(handles=[ld, ln, ld2, ln2], framealpha=0, prop={'weight': 'bold'}, loc='upper left')
-
-	# add xlabel, ylabel, suptitle
-	ax.set(xlabel='Date',
-		   ylabel='Ext (1/Mm)',
-		   xlim=(st_tm, fn_tm),
-		   ylim=(1., 600),
-		   xticks=tick_time,
-		   xticklabels=[_tm.strftime("%F %H:00") for _tm in tick_time])
-
-	ax2.set(ylabel=r'$VC (m^{2}/s)$',
-			ylim=(1., 2500))
-
-	fig.suptitle(f'Event evolution ({st_tm.strftime("%F")}_{fn_tm.strftime("%F")})')
-
-	# save figure
-	fig.savefig(pth(f"event_evolution_{st_tm.strftime("%F")}_{fn_tm.strftime("%F")}"))
-
-
-if __name__ == '__main__':
-	event_evolution(blh)