AeroViz 0.1.6__tar.gz → 0.1.8__tar.gz

aeroviz-0.1.8/AeroViz/data/240228_00.txt

@@ -0,0 +1,101 @@
+     1     1
+    AWRF    24     1    25     0     0
+     4 BACKWARD OMEGA   
+    24     2    27    16   22.630  120.346    10.0
+    24     2    27    16   22.630  120.346   100.0
+    24     2    27    16   22.630  120.346   500.0
+    24     2    27    16   22.630  120.346  1000.0
+     1 PRESSURE
+     1     1    24     2    27    16     0   808     0.0   22.630  120.346     10.0   1018.1
+     2     1    24     2    27    16     0   808     0.0   22.630  120.346    100.0   1007.1
+     3     1    24     2    27    16     0   808     0.0   22.630  120.346    500.0    961.2
+     4     1    24     2    27    16     0   808     0.0   22.630  120.346   1000.0    906.2
+     1     1    24     2    27    15     0   807    -1.0   22.700  120.337     13.4   1016.4
+     2     1    24     2    27    15     0   807    -1.0   22.779  120.325    115.8   1004.4
+     3     1    24     2    27    15     0   807    -1.0   22.615  120.313    471.0    964.8
+     4     1    24     2    27    15     0   807    -1.0   22.535  120.469   1035.2    902.9
+     1     1    24     2    27    14     0   806    -2.0   22.795  120.323     24.3   1015.4
+     2     1    24     2    27    14     0   806    -2.0   23.026  120.322    128.8   1001.6
+     3     1    24     2    27    14     0   806    -2.0   22.606  120.315    471.5    964.7
+     4     1    24     2    27    14     0   806    -2.0   22.484  120.626   1025.3    893.1
+     1     1    24     2    27    13     0   805    -3.0   22.936  120.308     56.6   1010.4
+     2     1    24     2    27    13     0   805    -3.0   23.275  120.355    142.2   1000.7
+     3     1    24     2    27    13     0   805    -3.0   22.612  120.341    492.9    961.7
+     4     1    24     2    27    13     0   805    -3.0   22.490  120.729    641.9    860.0
+     1     1    24     2    27    12     0   804    -4.0   23.133  120.283     74.9   1010.2
+     2     1    24     2    27    12     0   804    -4.0   23.491  120.413    187.7    994.2
+     3     1    24     2    27    12     0   804    -4.0   22.670  120.382    482.0    960.7
+     4     1    24     2    27    12     0   804    -4.0   22.538  120.774    830.3    801.8
+     1     1    24     2    27    11     0   803    -5.0   23.363  120.276     94.7   1007.0
+     2     1    24     2    27    11     0   803    -5.0   23.702  120.464    202.2    992.8
+     3     1    24     2    27    11     0   803    -5.0   22.708  120.455    327.6    977.5
+     4     1    24     2    27    11     0   803    -5.0   22.576  120.736    629.4    831.5
+     1     1    24     2    27    10     0   802    -6.0   23.611  120.281    123.4   1003.6
+     2     1    24     2    27    10     0   802    -6.0   23.924  120.507    235.0    988.9
+     3     1    24     2    27    10     0   802    -6.0   22.691  120.521     95.4   1004.1
+     4     1    24     2    27    10     0   802    -6.0   22.631  120.727    385.5    834.4
+     1     1    24     2    27     9     0   801    -7.0   23.910  120.338    185.0    996.9
+     2     1    24     2    27     9     0   801    -7.0   24.136  120.570    314.7    974.2
+     3     1    24     2    27     9     0   801    -7.0   22.702  120.585    125.1    999.7
+     4     1    24     2    27     9     0   801    -7.0   22.669  120.749    200.9    817.0
+     1     1    24     2    27     8     0   800    -8.0   24.287  120.521    321.6    982.5
+     2     1    24     2    27     8     0   800    -8.0   24.195  120.654    531.7    942.2
+     3     1    24     2    27     8     0   800    -8.0   22.707  120.608    170.4    993.9
+     4     1    24     2    27     8     0   800    -8.0   22.689  120.756    313.0    822.5
+     1     1    24     2    27     7     0   799    -9.0   24.612  120.810    428.8    969.2
+     2     1    24     2    27     7     0   799    -9.0   24.193  120.736    697.7    917.4
+     3     1    24     2    27     7     0   799    -9.0   22.719  120.610    229.7    986.5
+     4     1    24     2    27     7     0   799    -9.0   22.689  120.751    248.2    833.6
+     1     1    24     2    27     6     0   798   -10.0   24.859  121.146    538.3    932.4
+     2     1    24     2    27     6     0   798   -10.0   24.145  120.757    664.1    919.1
+     3     1    24     2    27     6     0   798   -10.0   22.738  120.605    305.9    978.9
+     4     1    24     2    27     6     0   798   -10.0   22.711  120.761      0.0    888.3
+     1     1    24     2    27     5     0   797   -11.0   24.984  121.430    392.9    975.2
+     2     1    24     2    27     5     0   797   -11.0   24.074  120.749    654.4    923.5
+     3     1    24     2    27     5     0   797   -11.0   22.731  120.603    495.6    958.3
+     4     1    24     2    27     5     0   797   -11.0   22.700  120.785    175.9    829.1
+     1     1    24     2    27     4     0   796   -12.0   25.158  121.638    503.0    954.6
+     2     1    24     2    27     4     0   796   -12.0   23.989  120.737    641.7    927.8
+     3     1    24     2    27     4     0   796   -12.0   22.693  120.592    648.1    943.6
+     4     1    24     2    27     4     0   796   -12.0   22.685  120.810    145.4    853.0
+     1     1    24     2    27     3     0   795   -13.0   25.451  121.825    386.1    978.5
+     2     1    24     2    27     3     0   795   -13.0   23.906  120.739    645.3    929.6
+     3     1    24     2    27     3     0   795   -13.0   22.657  120.583    682.1    940.8
+     4     1    24     2    27     3     0   795   -13.0   22.655  120.810      8.7    903.0
+     1     1    24     2    27     2     0   794   -14.0   25.789  121.989    450.1    971.0
+     2     1    24     2    27     2     0   794   -14.0   23.854  120.722    640.7    924.2
+     3     1    24     2    27     2     0   794   -14.0   22.618  120.552    686.5    941.2
+     4     1    24     2    27     2     0   794   -14.0   22.628  120.817    147.7    890.5
+     1     1    24     2    27     1     0   793   -15.0   26.127  122.136    554.7    959.4
+     2     1    24     2    27     1     0   793   -15.0   23.815  120.684    526.0    936.4
+     3     1    24     2    27     1     0   793   -15.0   22.537  120.538    763.0    933.1
+     4     1    24     2    27     1     0   793   -15.0   22.630  120.825    156.0    898.7
+     1     1    24     2    27     0     0   792   -16.0   26.425  122.289    711.1    940.7
+     2     1    24     2    27     0     0   792   -16.0   23.793  120.665    527.9    942.8
+     3     1    24     2    27     0     0   792   -16.0   22.498  120.642    901.2    897.5
+     4     1    24     2    27     0     0   792   -16.0   22.633  120.822    130.8    898.3
+     1     1    24     2    26    23     0   791   -17.0   26.705  122.454    872.6    921.4
+     2     1    24     2    26    23     0   791   -17.0   23.798  120.671    524.1    940.3
+     3     1    24     2    26    23     0   791   -17.0   22.514  120.689    734.8    863.9
+     4     1    24     2    26    23     0   791   -17.0   22.628  120.825     80.9    904.3
+     2     1    24     2    26    22     0   790   -18.0   23.816  120.688    493.1    939.0
+     3     1    24     2    26    22     0   790   -18.0   22.559  120.685    771.9    865.4
+     4     1    24     2    26    22     0   790   -18.0   22.623  120.829     56.0    907.8
+     2     1    24     2    26    21     0   789   -19.0   23.836  120.711    427.5    946.0
+     3     1    24     2    26    21     0   789   -19.0   22.557  120.666    660.4    899.0
+     4     1    24     2    26    21     0   789   -19.0   22.623  120.833     82.4    907.7
+     2     1    24     2    26    20     0   788   -20.0   23.865  120.751    489.8    941.5
+     3     1    24     2    26    20     0   788   -20.0   22.525  120.624    588.1    942.1
+     4     1    24     2    26    20     0   788   -20.0   22.636  120.848     98.3    914.3
+     2     1    24     2    26    19     0   787   -21.0   23.874  120.791    693.5    909.2
+     3     1    24     2    26    19     0   787   -21.0   22.492  120.608    624.9    943.0
+     4     1    24     2    26    19     0   787   -21.0   22.667  120.865    193.7    891.2
+     2     1    24     2    26    18     0   786   -22.0   23.874  120.793    743.4    902.1
+     3     1    24     2    26    18     0   786   -22.0   22.482  120.601    703.4    938.1
+     4     1    24     2    26    18     0   786   -22.0   22.747  120.984    287.0    906.5
+     2     1    24     2    26    17     0   785   -23.0   23.885  120.793    766.4    902.4
+     3     1    24     2    26    17     0   785   -23.0   22.471  120.582    842.0    922.5
+     4     1    24     2    26    17     0   785   -23.0   22.867  121.070    462.0    920.8
+     2     1    24     2    26    16     0   784   -24.0   23.895  120.778    859.2    898.6
+     3     1    24     2    26    16     0   784   -24.0   22.446  120.618   1090.1    893.0
+     4     1    24     2    26    16     0   784   -24.0   22.973  121.135    790.6    906.9

{aeroviz-0.1.6 → aeroviz-0.1.8}/AeroViz/dataProcess/Chemistry/_ocec.py

@@ -1,6 +1,8 @@
+import warnings
+
 import numpy as np
 from pandas import concat, DataFrame
-from scipy.optimize import curve_fit
+from scipy.optimize import curve_fit, least_squares, OptimizeWarning
 
 from AeroViz.dataProcess.core import union_index
 
@@ -22,15 +24,26 @@ def _min_Rsq(_oc, _ec, _rng):
     for _ocec, _out in _out_table.items():
         _df = DataFrame([_out.values, _ec.values]).T.dropna()
 
-        _x, _y = _df[0], _df[1]
-        _opt, _ = curve_fit(_func, _x, _y)
+        _x, _y = _df[0].values, _df[1].values
+
+        # 初始參數估計
+        slope_guess = (_y[-1] - _y[0]) / (_x[-1] - _x[0])
+        intercept_guess = _y[0] - slope_guess * _x[0]
+
+        try:
+            with warnings.catch_warnings():
+                warnings.filterwarnings('error')
+                _opt, _ = curve_fit(_func, _x, _y, p0=[slope_guess, intercept_guess], maxfev=5000)
+        except (RuntimeWarning, OptimizeWarning):
+            # 如果 curve_fit 失敗，嘗試使用 least_squares
+            residuals = lambda p: _func(_x, *p) - _y
+            _opt = least_squares(residuals, [slope_guess, intercept_guess]).x
 
-        _tss = np.sum((_y - _y.mean()) ** 2.)
-        _rss = np.sum((_y - _func(_x, *_opt)) ** 2.)
+        _tss = np.sum((_y - np.mean(_y)) ** 2)
+        _rss = np.sum((_y - _func(_x, *_opt)) ** 2)
 
-        _r2_dic[round(_ocec, 3)] = 1. - _rss / _tss
+        _r2_dic[round(_ocec, 3)] = 1 - _rss / _tss
 
-    # get the min R2
     _ratio = DataFrame(_r2_dic, index=[0]).idxmin(axis=1).values[0]
 
     return _ratio, _out_table[_ratio]

{aeroviz-0.1.6 → aeroviz-0.1.8}/AeroViz/plot/__init__.py

@@ -3,7 +3,9 @@ from . import meteorology
 from . import optical
 from .bar import bar
 from .box import box
+from .hysplit import hysplit
 from .pie import pie, donuts
+from .radar import radar
 from .regression import linear_regression, multiple_linear_regression
 from .scatter import scatter
 from .templates import *

aeroviz-0.1.8/AeroViz/plot/hysplit/__init__.py

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

aeroviz-0.1.8/AeroViz/plot/hysplit/hysplit.py

@@ -0,0 +1,79 @@
+from pathlib import Path
+
+import cartopy.crs as ccrs
+import cartopy.feature as cfeature
+import matplotlib.pyplot as plt
+import pandas as pd
+
+from AeroViz.plot.utils import set_figure
+
+# TODO: Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model
+
+
+__all__ = ['hysplit']
+
+
+def read_hysplit_data(file: Path):
+    data = pd.read_csv(file, skiprows=8, sep=r'\s+', names=range(0, 12), engine='python')
+    data = data.reset_index(drop=False)
+    data.columns = ['category', 'name', 'abc', 'year', 'month', 'hour', 'min', 'cont', 'backward', 'lat', 'lon',
+                    'height', 'pressure']
+
+    return data
+
+
+@set_figure
+def hysplit(file: Path):
+    data = read_hysplit_data(file)
+
+    # 創建地圖
+    fig, ax = plt.subplots(figsize=(4, 5), subplot_kw={'projection': ccrs.PlateCarree()})
+
+    # 設置地圖範圍
+    ax.set_extent([116, 126, 17, 30], crs=ccrs.PlateCarree())
+
+    # 添加自然地理特徵
+    ax.add_feature(cfeature.LAND)
+    ax.add_feature(cfeature.OCEAN)
+    ax.add_feature(cfeature.COASTLINE)
+    ax.add_feature(cfeature.BORDERS, linestyle=':')
+
+    # 添加經緯度網格
+    ax.gridlines(draw_labels=True, dms=True, x_inline=False, y_inline=False)
+
+    # 定義四種顏色
+    colors = ['red', 'blue', 'green', 'purple']
+
+    # 繪製四條軌跡線
+    group = data.groupby('category')
+    for i, (name, _data) in enumerate(group):
+        trajectory = _data
+        ax.plot(trajectory['lon'], trajectory['lat'], color=colors[i],
+                linewidth=2, transform=ccrs.Geodetic(),
+                label=f'Trajectory {name}')
+
+        # 添加起點和終點標記
+        # ax.plot(trajectory['lon'].iloc[-1], trajectory['lat'].iloc[-1], 'o',
+        #         color=colors[i], markersize=4, transform=ccrs.Geodetic())
+        # ax.plot(trajectory['lon'].iloc[0], trajectory['lat'].iloc[0], 's',
+        #         color=colors[i], markersize=4, transform=ccrs.Geodetic())
+
+    ax.legend(loc='upper right')
+    # 添加色標
+    # cbar = plt.colorbar(scatter, ax=ax, shrink=0.6, pad=0.12)
+    # cbar.set_label('Height (m)')
+
+    # 添加標題
+    plt.title("HYSPLIT model", pad=12)
+
+    plt.tight_layout()
+
+    # 保存地圖
+    plt.savefig('backward_hysplit.png', dpi=300, bbox_inches='tight')
+
+    # 顯示地圖（可選）
+    plt.show()
+
+
+if __name__ == "__main__":
+    hysplit(Path("../../data/240228_00.txt"))  # 請替換為您的實際檔案路徑

{aeroviz-0.1.6 → aeroviz-0.1.8}/AeroViz/plot/meteorology/meteorology.py

@@ -70,6 +70,7 @@ def wind_rose(df: DataFrame,
             rlabel_position=rlabel_pos,
             theta_direction=-1,
             theta_zero_location='N',
+            title=kwargs.get('title', None)
         )
         ax.set_thetagrids(angles=[0, 45, 90, 135, 180, 225, 270, 315],
                           labels=["N", "NE", "E", "SE", "S", "SW", "W", "NW"])
@@ -81,6 +82,7 @@ def wind_rose(df: DataFrame,
     return fig, ax
 
 
+# TODO: fix the bug of the CBPF function
 @set_figure(figsize=(4.3, 4))
 def CBPF(df: DataFrame,
          WS: Series | str,

{aeroviz-0.1.6 → aeroviz-0.1.8}/AeroViz/plot/optical/optical.py

@@ -6,13 +6,14 @@ import numpy as np
 # from PyMieScatt import ScatteringFunction
 from matplotlib.pyplot import Figure, Axes
 
+from AeroViz.plot.optical.PyMieScatt_update import ScatteringFunction
+from AeroViz.plot.optical.mie_theory import Mie_Q, Mie_MEE, Mie_PESD
 from AeroViz.plot.utils import *
-from AeroViz.process.method.mie_theory import Mie_Q, Mie_MEE, Mie_PESD
 
 __all__ = ['Q_plot',
            'RI_couple',
            'RRI_2D',
-           # 'scattering_phase',
+           'scattering_phase',
            'response_surface',
            ]
 
@@ -252,59 +253,59 @@ def RRI_2D(mode: Literal["ext", "sca", "abs"] = 'ext',
     return fig, ax
 
 
-# @set_figure
-# def scattering_phase(m: complex = 1.55 + 0.01j,
-#                      wave: float = 600,
-#                      dp: float = 200) -> tuple[Figure, Axes]:
-#     """
-#     Generate a polar plot to visualize the scattering phase function.
-#
-#     Parameters
-#     ----------
-#     m : complex, optional
-#         The complex refractive index of the scattering medium. Default is 1.55 + 0.01j.
-#     wave : float, optional
-#         The wavelength of the incident light in nanometers. Default is 600 nm.
-#     dp : float, optional
-#         The particle diameter in nanometers. Default is 200 nm.
-#
-#     Returns
-#     -------
-#     ax : Axes
-#         Matplotlib Axes object containing the generated polar plot.
-#
-#     Examples
-#     --------
-#     Example usage of the scattering_phase function:
-#
-#     >>> ax = scattering_phase(m=1.55 + 0.01j, wave=600, dp=200)
-#     """
-#     theta, _SL, _SR, _SU = ScatteringFunction(m, wave, dp)
-#
-#     SL = np.append(_SL, _SL[::-1])
-#     SR = np.append(_SR, _SR[::-1])
-#     SU = np.append(_SU, _SU[::-1])
-#
-#     angles = ['0', '60', '120', '180', '240', '300']
-#
-#     fig, ax = plt.subplots(subplot_kw={'projection': 'polar'})
-#
-#     theta = np.linspace(0, 2 * np.pi, len(SL))
-#
-#     plt.thetagrids(range(0, 360, int(360 / len(angles))), angles)
-#
-#     plt.plot(theta, SL, '-', linewidth=2, color='#115162', label='SL')
-#     plt.fill(theta, SL, '#afe0f5', alpha=0.5)
-#     plt.plot(theta, SR, '-', linewidth=2, color='#7FAE80', label='SR')
-#     plt.fill(theta, SR, '#b5e6c5', alpha=0.5)
-#     plt.plot(theta, SU, '-', linewidth=2, color='#621129', label='SU')
-#     plt.fill(theta, SU, '#f5afbd', alpha=0.5)
-#
-#     plt.legend(loc='best', bbox_to_anchor=(1, 0, 0.2, 1), prop={'weight': 'bold'})
-#     plt.title(r'$\bf Scattering\ phase\ function$')
-#
-#     plt.show()
-#     return fig, ax
+@set_figure
+def scattering_phase(m: complex = 1.55 + 0.01j,
+                     wave: float = 600,
+                     dp: float = 200) -> tuple[Figure, Axes]:
+    """
+    Generate a polar plot to visualize the scattering phase function.
+
+    Parameters
+    ----------
+    m : complex, optional
+        The complex refractive index of the scattering medium. Default is 1.55 + 0.01j.
+    wave : float, optional
+        The wavelength of the incident light in nanometers. Default is 600 nm.
+    dp : float, optional
+        The particle diameter in nanometers. Default is 200 nm.
+
+    Returns
+    -------
+    ax : Axes
+        Matplotlib Axes object containing the generated polar plot.
+
+    Examples
+    --------
+    Example usage of the scattering_phase function:
+
+    >>> ax = scattering_phase(m=1.55 + 0.01j, wave=600, dp=200)
+    """
+    theta, _SL, _SR, _SU = ScatteringFunction(m, wave, dp)
+
+    SL = np.append(_SL, _SL[::-1])
+    SR = np.append(_SR, _SR[::-1])
+    SU = np.append(_SU, _SU[::-1])
+
+    angles = ['0', '60', '120', '180', '240', '300']
+
+    fig, ax = plt.subplots(subplot_kw={'projection': 'polar'})
+
+    theta = np.linspace(0, 2 * np.pi, len(SL))
+
+    plt.thetagrids(range(0, 360, int(360 / len(angles))), angles)
+
+    plt.plot(theta, SL, '-', linewidth=2, color='#115162', label='SL')
+    plt.fill(theta, SL, '#afe0f5', alpha=0.5)
+    plt.plot(theta, SR, '-', linewidth=2, color='#7FAE80', label='SR')
+    plt.fill(theta, SR, '#b5e6c5', alpha=0.5)
+    plt.plot(theta, SU, '-', linewidth=2, color='#621129', label='SU')
+    plt.fill(theta, SU, '#f5afbd', alpha=0.5)
+
+    plt.legend(loc='best', bbox_to_anchor=(1, 0, 0.2, 1), prop={'weight': 'bold'})
+    plt.title(r'$\bf Scattering\ phase\ function$')
+
+    plt.show()
+    return fig, ax
 
 
 @set_figure
@@ -368,7 +369,7 @@ def response_surface(real_range=(1.33, 1.7),
     ax.plot_surface(real, gmd, ext, rstride=1, cstride=1, cmap=plt.get_cmap('jet'), edgecolor='none')
 
     ax.set(xlabel='Real part (n)', ylabel='GMD (nm)', zlabel=Unit('Extinction'),
-           title='Sensitive tests of Extinction')
+           title='Sensitive tests of extinction')
 
     ax.zaxis.get_offset_text().set_visible(False)
     exponent = math.floor(math.log10(np.max(ext)))
@@ -381,8 +382,8 @@ def response_surface(real_range=(1.33, 1.7),
 
 
 if __name__ == '__main__':
-    # Q_plot(['AS', 'AN', 'OM', 'Soil', 'SS', 'BC'], x='dp', y='MEE')
-    # Q_plot(['AS', 'AN', 'OM', 'Soil', 'SS', 'BC'], x='dp', y='Q')
+    Q_plot(['AS', 'AN', 'OM', 'Soil', 'SS', 'BC'], x='dp', y='MEE')
+    Q_plot(['AS', 'AN', 'OM', 'Soil', 'SS', 'BC'], x='dp', y='Q')
 
-    # RI_couple()
+    RI_couple()
     response_surface()

{aeroviz-0.1.6 → aeroviz-0.1.8}/AeroViz/plot/pie.py

@@ -81,7 +81,9 @@ def pie(data_set: DataFrame | dict,
     radius = 4
     width = 4 if style == 'pie' else 1
 
-    text = [''] * pies if style == 'pie' else [Unit(unit) + '\n\n' + '{:.2f}'.format(x) for x in data.sum(axis=1)]
+    text = [''] * pies if style == 'pie' else [Unit(unit) + '\n\n' +
+                                               '{:.2f} ± {:.2f}'.format(x, s)
+                                               for x, s in zip(data.sum(axis=1), data.std(axis=1))]
     pct_distance = 0.6 if style == 'pie' else 0.88
 
     fig, ax = plt.subplots(1, pies, figsize=((pies * 2) + 1, 2)) if ax is None else (ax.get_figure(), ax)
@@ -99,7 +101,17 @@ def pie(data_set: DataFrame | dict,
                   pctdistance=1.3, radius=radius, wedgeprops=dict(width=width, edgecolor='w'))
         ax[i].axis('equal')
         ax[i].text(0, 0, text[i], ha='center', va='center')
-        ax[i].set_title(category_names[i])
+
+        if kwargs.get('title') is None:
+            ax[i].set_title(category_names[i])
+
+        else:
+            if len(kwargs.get('title')) == pies:
+                title = kwargs.get('title')
+            else:
+                raise ValueError('The length of the title list must match the number of pies.')
+
+            ax[i].set_title(title[i])
 
     ax[-1].legend(labels, loc='center left', prop={'size': 8, 'weight': 'normal'}, bbox_to_anchor=(1, 0, 1.15, 1))
 

aeroviz-0.1.8/AeroViz/plot/radar.py

@@ -0,0 +1,184 @@
+import matplotlib.pyplot as plt
+import numpy as np
+from matplotlib.patches import Circle, RegularPolygon
+from matplotlib.path import Path
+from matplotlib.projections import register_projection
+from matplotlib.projections.polar import PolarAxes
+from matplotlib.spines import Spine
+from matplotlib.transforms import Affine2D
+
+from AeroViz.plot.utils import *
+
+__all__ = ['radar']
+
+
+def radar_factory(num_vars, frame='circle'):
+    """
+    Create a radar chart with `num_vars` axes.
+
+    This function creates a RadarAxes projection and registers it.
+
+    Parameters
+    ----------
+    num_vars : int
+        Number of variables for radar chart.
+    frame : {'circle', 'polygon'}
+        Shape of frame surrounding axes.
+
+    """
+    # calculate evenly-spaced axis angles
+    theta = np.linspace(0, 2 * np.pi, num_vars, endpoint=False)
+
+    class RadarTransform(PolarAxes.PolarTransform):
+
+        def transform_path_non_affine(self, path):
+            # Paths with non-unit interpolation steps correspond to gridlines,
+            # in which case we force interpolation (to defeat PolarTransform's
+            # autoconversion to circular arcs).
+            if path._interpolation_steps > 1:
+                path = path.interpolated(num_vars)
+            return Path(self.transform(path.vertices), path.codes)
+
+    class RadarAxes(PolarAxes):
+
+        name = 'radar'
+        PolarTransform = RadarTransform
+
+        def __init__(self, *args, **kwargs):
+            super().__init__(*args, **kwargs)
+            # rotate plot such that the first axis is at the top
+            self.set_theta_zero_location('N')
+
+        def fill(self, *args, closed=True, **kwargs):
+            """Override fill so that line is closed by default"""
+            return super().fill(closed=closed, *args, **kwargs)
+
+        def plot(self, *args, **kwargs):
+            """Override plot so that line is closed by default"""
+            lines = super().plot(*args, **kwargs)
+            for line in lines:
+                self._close_line(line)
+
+        def _close_line(self, line):
+            x, y = line.get_data()
+            if x[0] != x[-1]:
+                x = np.append(x, x[0])
+                y = np.append(y, y[0])
+                line.set_data(x, y)
+
+        def set_varlabels(self, labels):
+            self.set_thetagrids(np.degrees(theta), labels)
+
+        @staticmethod
+        def _gen_axes_patch():
+            # The Axes patch must be centered at (0.5, 0.5) and of radius 0.5
+            # in axes coordinates.
+            if frame == 'circle':
+                return Circle((0.5, 0.5), 0.5)
+            elif frame == 'polygon':
+                return RegularPolygon((0.5, 0.5), num_vars,
+                                      radius=.5, edgecolor="k")
+            else:
+                raise ValueError("Unknown value for 'frame': %s" % frame)
+
+        def _gen_axes_spines(self):
+            if frame == 'circle':
+                return super()._gen_axes_spines()
+            elif frame == 'polygon':
+                # spine_type must be 'left'/'right'/'top'/'bottom'/'circle'.
+                spine = Spine(axes=self,
+                              spine_type='circle',
+                              path=Path.unit_regular_polygon(num_vars))
+                # unit_regular_polygon gives a polygon of radius 1 centered at
+                # (0, 0) but we want a polygon of radius 0.5 centered at (0.5,
+                # 0.5) in axes coordinates.
+                spine.set_transform(Affine2D().scale(.5).translate(.5, .5)
+                                    + self.transAxes)
+                return {'polar': spine}
+            else:
+                raise ValueError("Unknown value for 'frame': %s" % frame)
+
+    register_projection(RadarAxes)
+    return theta
+
+
+@set_figure(figsize=(3, 3))
+def radar(data, labels=None, legend_labels=None, **kwargs) -> tuple[plt.Figure, plt.Axes]:
+    """
+    Creates a radar chart based on the provided data.
+
+    Parameters
+    ----------
+    data : list of list
+        A 2D list where each inner list represents a factor, and each element
+        within the inner lists represents a value for a species.
+        Shape: (n_factors, n_species)
+        Example: [[0.88, 0.01, 0.03, ...], [0.07, 0.95, 0.04, ...], ...]
+
+    labels : list, optional
+        A list of strings representing the names of species (variables).
+        If provided, it should have the same length as the number of elements
+        in each inner list of `data`.
+        Example: ['Sulfate', 'Nitrate', 'EC', 'OC1', 'OC2', 'OP', 'CO', 'O3']
+
+    legend_labels : list, optional
+        A list of strings for labeling each factor in the legend.
+        If provided, it should have the same length as the number of inner lists in `data`.
+
+    **kwargs : dict
+        Additional keyword arguments to be passed to the plotting function.
+        This may include 'title' for setting the chart title.
+
+    Returns
+    -------
+    tuple[plt.Figure, plt.Axes]
+        A tuple containing the Figure and Axes objects of the created plot.
+
+    Example
+    -------
+    >>> data = [[0.88, 0.01, 0.03, 0.03, 0.00, 0.06, 0.01, 0.00],
+    >>>         [0.07, 0.95, 0.04, 0.05, 0.00, 0.02, 0.01, 0.00],
+    >>>         [0.01, 0.02, 0.85, 0.19, 0.05, 0.10, 0.00, 0.00],
+    >>>         [0.02, 0.01, 0.07, 0.01, 0.21, 0.12, 0.98, 0.00],
+    >>>         [0.01, 0.01, 0.02, 0.71, 0.74, 0.70, 0.30, 0.20]]
+    >>> labels = ['Sulfate', 'Nitrate', 'EC', 'OC1', 'OC2', 'OP', 'CO', 'O3']
+    >>> fig, ax = radar(data, labels=labels, title='Basecase')
+
+    Note
+    ----
+    The first dimension of `data` represents each factor, while the second
+    dimension represents each species.
+    """
+    theta = radar_factory(np.array(data).shape[1], frame='polygon')
+
+    fig, ax = plt.subplots(subplot_kw=dict(projection='radar'))
+    fig.subplots_adjust(wspace=0.25, hspace=0.20, top=0.80, bottom=0.05, right=0.80)
+
+    colors = ['b', 'r', 'g', 'm', 'y']
+
+    # Plot the four cases from the example data on separate axes
+    for d, color in zip(data, colors):
+        ax.plot(theta, d, color=color)
+        ax.fill(theta, d, facecolor=color, alpha=0.25, label='_nolegend_')
+
+    ax.set_varlabels(labels)
+    ax.set_rgrids([0.2, 0.4, 0.6, 0.8])
+    ax.set(title=kwargs.get('title', ''))
+
+    # add legend relative to top-left plot
+    legend_labels = legend_labels or ('Factor 1', 'Factor 2', 'Factor 3', 'Factor 4', 'Factor 5')
+    legend = ax.legend(legend_labels, loc=(0.95, 0.95), labelspacing=0.1)
+
+    plt.show()
+
+    return fig, ax
+
+
+if __name__ == '__main__':
+    data = [[0.88, 0.01, 0.03, 0.03, 0.00, 0.06, 0.01, 0.00],
+            [0.07, 0.95, 0.04, 0.05, 0.00, 0.02, 0.01, 0.00],
+            [0.01, 0.02, 0.85, 0.19, 0.05, 0.10, 0.00, 0.00],
+            [0.02, 0.01, 0.07, 0.01, 0.21, 0.12, 0.98, 0.00],
+            [0.01, 0.01, 0.02, 0.71, 0.74, 0.70, 0.30, 0.20]]
+
+    fig, ax = radar(data=data, labels=['Sulfate', 'Nitrate', 'EC', 'OC1', 'OC2', 'OP', 'CO', 'O3'], title='Basecase')