ararpy 0.1.11__py3-none-any.whl → 0.1.12__py3-none-any.whl

ararpy/smp/basic.py

@@ -534,10 +534,11 @@ def set_params(smp: Sample, params: Union[List, str], flag: Optional[str] = None
         smp.TotalParam[58:67] = remove_none(smp.TotalParam[58:67], params[4:13], n, 67 - 58)
         smp.TotalParam[97:100] = remove_none(smp.TotalParam[97:100], params[13:16], n, 100 - 97)
         smp.TotalParam[115:120] = remove_none(smp.TotalParam[115:120], params[16:21], n, 120 - 115)
-        smp.TotalParam[120:123] = remove_none(smp.TotalParam[120:123], params[21:24], n, 123 - 120)
+        smp.TotalParam[126:136] = remove_none(smp.TotalParam[126:136], params[21:31], n, 136 - 126)
+        smp.TotalParam[120:123] = remove_none(smp.TotalParam[120:123], params[31:34], n, 123 - 120)
         smp.TotalParam[100:114] = remove_none(
             smp.TotalParam[100:114],
-            [['Linear', 'Exponential', 'Power'][params[24:27].index(True)] if True in params[24:27] else '', *params[27:]], n, 114 - 100)
+            [['Linear', 'Exponential', 'Power'][params[34:37].index(True)] if True in params[34:37] else '', *params[37:]], n, 114 - 100)
     else:
         raise KeyError(f"{flag = } is not supported. It must be 'calc' for Calc Params, "
                        f"'irra' for Irradiation Params, or 'smp' for Sample Params.")

ararpy/smp/calculation.py

@@ -47,6 +47,7 @@ def recalculate(
     re_plot_style
     re_set_table
     re_table_style
+    re_corr_gain
     kwargs
 
     Returns
@@ -59,11 +60,14 @@ def recalculate(
     # print(f"{sample.SelectedSequence1 = }")
     # print(f"{sample.SelectedSequence2 = }")
     # --- initializing ---
+    re_corr_gain = re_corr_blank
     if re_initial:  # 1
         initial.re_set_smp(sample)
     # --- calculating ---
     if re_corr_blank:  # 2
         corr.corr_blank(sample)
+    if re_corr_gain:  # 2 2024-10-04 add
+        corr.corr_gain(sample)
     if re_corr_massdiscr:  # 3
         corr.corr_massdiscr(sample)
     if re_corr_decay:  # 4

ararpy/smp/corr.py

@@ -45,6 +45,32 @@ def corr_blank(sample: Sample):
     sample.CorrectedValues = copy.deepcopy(sample.BlankCorrected)
 
 
+# =======================
+# Corr Blank
+# =======================
+def corr_gain(sample: Sample):
+    """Blank Correction"""
+    corrGain = True
+    # if not corrBlank:
+    #     sample.BlankCorrected = copy.deepcopy(sample.SampleIntercept)
+    #     sample.CorrectedValues = copy.deepcopy(sample.BlankCorrected)
+    #     return
+    gain_corrected = [[]] * 10
+    try:
+        for i in range(5):
+            gain_corrected[i * 2:2 + i * 2] = calc.corr.gain(
+                *sample.BlankCorrected[i * 2:2 + i * 2], *sample.TotalParam[126 + i * 2:128 + i * 2])
+    except Exception as e:
+        print(f"Gain correction failed")
+        print(traceback.format_exc())
+        return
+    for i in range(0, 10, 2):
+        gain_corrected[i] = [gain_corrected[i][index] if corrGain else j for index, j in enumerate(sample.BlankCorrected[i])]
+        gain_corrected[i + 1] = [gain_corrected[i + 1][index] if corrGain else j for index, j in enumerate(sample.BlankCorrected[i + 1])]
+    sample.BlankCorrected = copy.deepcopy(gain_corrected)
+    sample.CorrectedValues = copy.deepcopy(sample.BlankCorrected)
+
+
 # =======================
 # Mass Discrimination
 # =======================

ararpy/smp/diffusion_funcs.py

@@ -167,7 +167,7 @@ class DiffArrmultiFunc(DiffSample):
         print(f"{self.telab = }")
         print(f"{self.tilab = }")
 
-    def main(self):
+    def main(self, **params):
 
         print(f"\n======================================")
         print(f"Run Arrmulti Main | ")
@@ -1436,7 +1436,7 @@ class DiffArrmultiFunc(DiffSample):
         return xran, yran
 
     def arr(self, e, ord, dzx, fmod, f, telab, tilab, xlogd,
-            ns=200, r=1.987e-3, pi=3.141592654, ee=0.4342944879, acut=0, b=1, imp=1):
+            ns=200, r=1.987e-3, pi=3.141592654, ee=0.4342944879, acut=0.5, b=8, imp=2):
 
         # print(f"{e = }")
         # print(f"{ord = }")
@@ -3677,13 +3677,10 @@ class InsideTemperatureCalibration:
 
 
 class SmpTemperatureCalibration:
-    def __init__(self, arr_path="", helix_log_path="", libano_log_path=None,
+    def __init__(self, helix_log_path="", libano_log_path=None,
                  name="smp_example", loc=r'D:\DjangoProjects\webarar\static\settings'):
 
         self.name = name
-        self.smp = ap.from_arr(file_path=arr_path)
-        if isinstance(self.smp, Sample):
-            self.name = self.smp.name()
 
         self.loc = loc
 
@@ -3858,27 +3855,28 @@ class SmpTemperatureCalibration:
 
         file_out.close()
 
+        np.savetxt(os.path.join(self.loc, f"{self.name}-temp.txt"), libano_log, delimiter=',')
+        heating_out = open(os.path.join(self.loc, f"{self.name}-heated-index.txt"), "w")
+        heating_out.writelines([f'{min(i)},{max(i)}\n' for i in yellow_data_index])
+        heating_out.close()
 
-        fig, ax = plt.subplots()
-        fig.set_size_inches(w=12, h=10)
+        plt_plot = False
+        if plt_plot:
 
-        ax.plot(libano_log[0], libano_log[1], c='green')
-        ax.plot(libano_log[0], libano_log[2], c='blue')
-        ax.plot(libano_log[0], libano_log[4], c='red')
-        ax.plot(libano_log[0], libano_log[5], c='red')
-        for i in range(nstep):
-            _index = yellow_data_index[i]
-            ax.plot(libano_log[0, _index], libano_log[4, _index], c='yellow')
-            ax.plot(libano_log[0, _index], libano_log[5, _index], c='yellow')
+            fig, ax = plt.subplots()
+            fig.set_size_inches(w=12, h=10)
 
-        fig.tight_layout()
-        plt.show()
+            ax.plot(libano_log[0], libano_log[1], c='green')
+            ax.plot(libano_log[0], libano_log[2], c='blue')
+            ax.plot(libano_log[0], libano_log[4], c='red')
+            ax.plot(libano_log[0], libano_log[5], c='red')
+            for i in range(nstep):
+                _index = yellow_data_index[i]
+                ax.plot(libano_log[0, _index], libano_log[4, _index], c='yellow')
+                ax.plot(libano_log[0, _index], libano_log[5, _index], c='yellow')
 
-
-        np.savetxt(os.path.join(self.loc, f"{self.name}-temp.txt"), libano_log, delimiter=',')
-        heating_out = open(os.path.join(self.loc, f"{self.name}-heated-index.txt"), "w")
-        heating_out.writelines([','.join([str(j) for j in i]) for i in yellow_data_index])
-        heating_out.close()
+            fig.tight_layout()
+            plt.show()
 
 
         #
@@ -4613,15 +4611,15 @@ def dr2_lovera(f, ti, ar, sar):
     # Lovera
     f = np.array(f)
     ti = np.array(ti)
-    ar = np.array(ar)
-    sar = np.array(sar)
+    # ar = np.array(ar)
+    # sar = np.array(sar)
     ti = ti * 60  # in seconds
     pi = math.pi
     pi = 3.141592654
 
     # sf = sqrt(b ** 2 * siga ** 2 + a ** 2 * sigb ** 2) / (a + b) ** 2
-    sf = [math.sqrt(sar[i + 1:].sum() ** 2 * (sar[:i + 1] ** 2).sum() + ar[:i + 1].sum() ** 2 * (sar[i + 1:] ** 2).sum()) / ar.sum() ** 2 for i in range(len(ar) - 1)]
-    sf.append(0)
+    # sf = [math.sqrt(sar[i + 1:].sum() ** 2 * (sar[:i + 1] ** 2).sum() + ar[:i + 1].sum() ** 2 * (sar[i + 1:] ** 2).sum()) / ar.sum() ** 2 for i in range(len(ar) - 1)]
+    # sf.append(0)
 
     f = np.where(f >= 1, 0.9999999999999999, f)
 
@@ -4631,41 +4629,67 @@ def dr2_lovera(f, ti, ar, sar):
     dr2 = [(dtr2[i] - (dtr2[i - 1] if i > 0 else 0)) / ti[i] * imp ** 2 for i in range(len(dtr2))]
     xlogd = [np.log10(i) for i in dr2]
 
-    wt = errcal(f, ti, a39=ar, sig39=sar)
+    wt = errcal(f, ti, ar, sar)
 
     return dr2, xlogd, wt
 
 
-def dr2_yang(f, ti):
+def dr2_yang(f, ti, ar, sar):
+    """
+
+        if i == 0:
+            dr2[i] = _dr2(f[i]) / ti[i]
+        else:
+            dr2[i] = _dr2((f[i] - f[i - 1]) / (1 - f[i - 1])) / ti[i]
+
+    Parameters
+    ----------
+    f: array NOTE：这里的 f 既不是释放百分数也不是累积百分数，而是各个阶段释放量占其及后续所有气体的分数，即占释放前的比例
+    ti
+
+    Returns
+    -------
+
+    """
     f = np.array(f)
     ti = np.array(ti)
     n = min(len(f), len(ti))
     ti = ti * 60  # in seconds
+    sti = 5.
     pi = math.pi
     pi = 3.141592654
+    ee = 0.4342944819
     dr2 = np.zeros(n)
+    sdr2 = np.zeros(n)
 
-    if f[-1] >= 1:
-        f[-1] = 0.99999999
+    f = np.where(f >= 1, 0.9999999999999999, f)
 
     def _dr2(_fi):
+        # return: dr2, dy/df
         if _fi <= 0.85:
-            return (1 - math.sqrt(1 - pi * _fi / 3)) ** 2 / pi
+            return (1 - math.sqrt(1 - pi * _fi / 3)) ** 2 / pi, 1 / (3 * math.sqrt(1 - pi * _fi / 3)) - 1 / 3
         else:
-            return math.log((1 - _fi) / (6 / pi ** 2)) / - (pi ** 2)
+            return math.log((1 - _fi) / (6 / pi ** 2)) / - (pi ** 2), 1 / (pi ** 2 * (1 - _fi))
 
-    for i in range(len(f)):
-        if i == 0:
-            dr2[i] = _dr2(f[i]) / ti[i]
-        else:
-            dr2[i] = _dr2((f[i] - f[i - 1]) / (1 - f[i - 1])) / ti[i]
+    for i in range(n):
+        dtr2, d = _dr2((f[i] - (f[i-1] if i > 0 else 0)) / (1 - (f[i-1] if i > 0 else 0)))
+        # sf 中 f = (f[i] - f[i - 1]) / (1 - f[i - 1])
+        sf = sum([ar[i] ** 2 * sar[j] ** 2 / sum(ar[i:]) ** 4 for j in range(i+1, n)]) + sum(ar[i+1:]) ** 2 * sar[i] ** 2 / sum(ar[i:]) ** 4
+        sdtr2 = math.sqrt(d ** 2 * sf)
+        dr2[i] = dtr2 / ti[i]
+        sdr2[i] = ap.calc.err.div((dtr2, sdtr2), (ti[i], sti))
+
+    xlogd = np.log(dr2)
+    sxlogd = np.abs(sdr2 / dr2)
 
     print(f"yang {dr2 = }")
+    print(f"yang {xlogd = }")
+    print(f"yang {sxlogd = }")
 
-    return dr2, np.log(dr2)
+    return dr2, xlogd, sxlogd
 
 
-def errcal(f, ti, a39, sig39):
+def errcal(f, ti, ar, sar):
 
     # ns = 200
     # r = 1.987E-3
@@ -4681,9 +4705,9 @@ def errcal(f, ti, a39, sig39):
     sigt0 = 90.
 
     ni = len(f)
-    sumat = sum(a39)
-    sigsm = [i / sumat for i in sig39]
-    siga = [sig39[i] / a39[i] for i in range(ni)]
+    sumat = sum(ar)
+    sigsm = [i / sumat for i in sar]
+    siga = [sar[i] / ar[i] for i in range(ni)]
 
     an1 = math.pi ** 2
     sigat = 0.

ararpy/smp/export.py

@@ -41,6 +41,52 @@ def to_pdf(file_path: str, figure: str, smp: Sample):
     pdf.save(figure=figure)
 
 
+def export_chart_to_pdf(data: dict):
+    # create a canvas
+    cv = pm.Canvas(width=17, height=12, unit="cm", show_frame=True, clip_outside_plot_areas=False)
+    # change frame outline style
+    cv.show_frame(color="grey", line_width=0.5)
+    sc = (*data['xAxis'][0]['extent'], *data['yAxis'][0]['extent'])
+    pt = cv.add_plot_area(name="Plot1", plot_area=(0.15, 0.15, 0.8, 0.8), plot_scale=sc, show_frame=True)
+    # draw axis
+    for stick in data['xAxis'][0]['interval']:
+        start = pt.scale_to_points(stick, sc[2])
+        end = pt.scale_to_points(stick, sc[2])
+        end = (end[0], end[1] - 5)
+        pt.line(start=start, end=end, width=1, line_style="solid", y_clip=False, coordinate="pt", z_index=100)
+        pt.text(x=start[0], y=end[1] - 15, text=f"{stick}", clip=False,
+                coordinate="pt", h_align="middle", z_index=150)
+    for stick in data['yAxis'][0]['interval']:
+        start = pt.scale_to_points(sc[0], stick)
+        end = pt.scale_to_points(sc[0], stick)
+        end = (end[0] - 5, end[1])
+        pt.line(start=start, end=end, width=1, line_style="solid", x_clip=False, coordinate="pt", z_index=100)
+        pt.text(x=end[0] - 5, y=end[1], text=f"{stick}", clip=False,
+                coordinate="pt", h_align="right", v_align="center", z_index=150)
+    # axis titles
+    nameLocation = pt.scale_to_points(sum(sc[:2]) / 2, sc[2])
+    pt.text(x=nameLocation[0], y=nameLocation[1] - 30, text=data['xAxis'][0]['title'], clip=False, coordinate="pt",
+            h_align="middle", v_align="top", z_index=150)
+    nameLocation = pt.scale_to_points(sc[0], sum(sc[2:4]) / 2)
+    pt.text(x=nameLocation[0] - 50, y=nameLocation[1], text=data['yAxis'][0]['title'], clip=False, coordinate="pt",
+            h_align="middle", v_align="bottom", rotate=90, z_index=150)
+    # draw series
+    for se in data['series']:
+        data = se.get('data', [])
+        if 'line' in se['type']:
+            for index in range(1, len(data)):
+                pt.line(start=data[index - 1], end=data[index], width=1, line_style='solid', name=se['name'],
+                        color=se.get('color', 'black'), clip=True, line_caps=se.get('line_caps', 'none'), z_index=9)
+        if 'scatter' in se['type'] and se['name'] != 'Text':
+            for each in data:
+                pt.scatter(each[0], each[1], fill_color=se.get('color', 'black'), size=2)
+        if 'scatter' in se['type'] and se['name'] == 'Text' or 'text' in se['type']:
+            for each in data:
+                pt.text(*each[:2], **se)
+
+    return cv
+
+
 class ExcelTemplate:
     def __init__(self, **kwargs):
         self.name = ""

ararpy/smp/json.py

@@ -22,7 +22,8 @@ def dumps(a):
 
 
 def loads(a):
-    return json.loads(a)
+    # null will be converted to None by default, replace None with np.nan
+    return json.loads(a, object_hook=myHook)
 
 
 class MyEncoder(json.JSONEncoder):
@@ -51,3 +52,13 @@ class MyEncoder(json.JSONEncoder):
         if not isinstance(obj, (int, str, list, dict, tuple, float)):
             print(f"Special type, {type(obj) = }, {obj = }")
         return super(MyEncoder, self).default(obj)
+
+
+def myHook(obj):
+    if isinstance(obj, dict):
+        return {k: myHook(v) for k, v in obj.items()}
+    elif isinstance(obj, list):
+        return [myHook(v) for v in obj]
+    elif obj is None:
+        return np.nan
+    return obj

ararpy/smp/plots.py

@@ -61,16 +61,20 @@ def set_plot_data(sample: Sample, isInit: bool = True, isIsochron: bool = True,
         try:
             initial_plot_data(sample)
         except (Exception, BaseException):
-            print(traceback.format_exc())
+            # print(traceback.format_exc())
             pass
 
     # Recalculate isochron lines
     if isIsochron:
         try:
             recalc_isochrons(sample, **kwargs)
+        except (Exception, BaseException):
+            # print(traceback.format_exc())
+            pass
+        try:
             reset_isochron_line_data(sample)
         except (Exception, BaseException):
-            print(traceback.format_exc())
+            # print(traceback.format_exc())
             pass
 
     # Recalculate plateaus
@@ -78,7 +82,7 @@ def set_plot_data(sample: Sample, isInit: bool = True, isIsochron: bool = True,
         try:
             recalc_plateaus(sample)
         except (Exception, BaseException):
-            print(traceback.format_exc())
+            # print(traceback.format_exc())
             pass
 
 
@@ -237,7 +241,7 @@ def get_isochron_results(data: list, smp: Sample, sequence, figure_type: int = 0
     try:
         regression_res = regression_method(*data[:5])
     except (Exception, BaseException):
-        print(f"Warning: {traceback.format_exc()}")
+        # print(f"Warning: {traceback.format_exc()}")
         return iso_res
     else:
         iso_res.update(dict(zip(reg_res_index, regression_res)))
@@ -305,11 +309,11 @@ def get_3D_results(data: list, sequence: list, sample: Sample):
             Q = 1 - np.exp(-1 * sample.TotalParam[46][0] * sum(sample.TotalParam[32]) / len(sample.TotalParam[32]))
             P = PQ / Q
         except:
-            print(f"Warning: {traceback.format_exc()}")
+            # print(f"Warning: {traceback.format_exc()}")
             P = 0
         age = basic.calc_age([f, sf], smp=sample)
     except:
-        print(f"Warning: {traceback.format_exc()}")
+        # print(f"Warning: {traceback.format_exc()}")
         k = [0] * 15
         age = [0] * 4
         ar40ar36, sar40ar36, P = 0, 0, 0

ararpy/smp/sample.py

@@ -119,8 +119,8 @@ TOTAL_PARAMS_HEADERS = [
     '\u2074\u2070Ar/\u00B3\u2076Ar air', '%1\u03C3',
     '\u00B3\u2078Ar/\u00B3\u2076Ar air', '%1\u03C3',  # 95-98
     'Isochron Fitting', 'Convergence', 'Iteration', 'Discrimination',  # 99-102
-    'Not Zero', 'Corr Blank', 'Corr Discr', 'Corr \u00B3\u2077Ar Decay',
-    'Corr \u00B3\u2079Ar Decay',  # 103-107
+    'Not Zero', 'Corr Blank', 'Corr Discr', 'Corr \u00B3\u2077Ar Decay',  # 103-106
+    'Corr \u00B3\u2079Ar Decay',  # 107
     'Ca Degassing', 'K Degassing', 'Cl Degassing', 'Trap Degassing',  # 108-111
     'Using Min Equation',
     # 'Recalibration', 'Using Std Age', 'Use Std Ratio',  # 112-115  to be completed
@@ -137,6 +137,11 @@ TOTAL_PARAMS_HEADERS = [
     'Heating Time (s)',  # 125
     'Heating Actual Temp (C)',  # 126
     'Heating AT 1\u03C3',  # 127
+    '36Ar Gain', '1\u03C3',  # 128-129
+    '37Ar Gain', '1\u03C3',  # 130-131
+    '38Ar Gain', '1\u03C3',  # 132-133
+    '39Ar Gain', '1\u03C3',  # 134-135
+    '40Ar Gain', '1\u03C3',  # 136-137
 ]
 
 SAMPLE_INTERCEPT_SHORT_HEADERS = [
@@ -245,6 +250,11 @@ TOTAL_PARAMS_SHORT_HEADERS = [
     'HeatingTime',  # 125
     'HeatingActualTemp',  # 126
     'HeatingActualTempError',  # 127
+    '36Gain', '1s',  # 128-129
+    '37Gain', '1s',  # 130-131
+    '38Gain', '1s',  # 132-133
+    '39Gain', '1s',  # 134-135
+    '40Gain', '1s',  # 136-137
 ]
 
 DEFAULT_PLOT_STYLES = {
@@ -731,7 +741,7 @@ DEFAULT_PLOT_STYLES = {
     },
 }
 
-VERSION = '20240730'
+VERSION = '20240930'
 
 NAMED_DICT = {
     "unknown": {"header": SAMPLE_INTERCEPT_HEADERS.copy()},
@@ -854,7 +864,7 @@ class Sample:
 
     def sequence(self) -> ArArBasic: ...
 
-    def recalculate(self, *args): ...
+    def recalculate(self, *args, **kwargs): ...
 
     def plot_age_plateau(self): ...
 

ararpy/smp/table.py

@@ -9,6 +9,7 @@
 #
 #
 """
+import copy
 from .. import calc
 from . import (sample as samples, basic)
 
@@ -47,6 +48,13 @@ def update_table_data(smp: Sample, only_table: str = None):
             data = calc.arr.merge(
                 smp.SequenceName, smp.SequenceValue, *smp.DegasValues)
         elif key == '5':
+            smp.PublishValues[0] = copy.deepcopy(smp.DegasValues[ 0])
+            smp.PublishValues[1] = copy.deepcopy(smp.DegasValues[ 8])
+            smp.PublishValues[2] = copy.deepcopy(smp.DegasValues[10])
+            smp.PublishValues[3] = copy.deepcopy(smp.DegasValues[20])
+            smp.PublishValues[4] = copy.deepcopy(smp.DegasValues[24])
+            smp.PublishValues[5:7] = copy.deepcopy(smp.ApparentAgeValues[2:4])
+            smp.PublishValues[7:9] = copy.deepcopy(smp.ApparentAgeValues[6:8])
             data = calc.arr.merge(
                 smp.SequenceName, smp.SequenceValue, *smp.PublishValues)
         elif key == '6':

ararpy/test.py

@@ -33,5 +33,68 @@ def test():
     print(sample.blank().to_df().iloc[:, [1, 2, 3]])
 
 
+def export_pdf_demo():
+    import numpy as np
+    import pdf_maker as pm
+
+    example_dir = os.path.join(os.path.dirname(os.path.abspath(__file__)), r'examples')
+
+    # ------ 将以下五个样品的年龄谱图组合到一起 -------
+    arr_files = [
+        os.path.join(example_dir, r'22WHA0433.arr'),
+        os.path.join(example_dir, r'20WHA0103.age'),
+    ]
+    colors = ['#1f3c40', '#e35000', '#e1ae0f', '#3d8ebf']
+    series = []
+
+    # ------ 构建数据 -------
+    for index, file in enumerate(arr_files):
+        smp = ap.from_arr(file_path=file) if file.endswith('.arr') else ap.from_age(file_path=file)
+        age = smp.ApparentAgeValues[2:4]
+        ar = smp.DegasValues[20]
+        data = ap.calc.spectra.get_data(*age, ar, cumulative=False)
+        series.append({
+            'type': 'series.line', 'id': f'line{index * 2 + 0}', 'name': f'line{index * 2 + 0}', 'color': colors[index],
+            'data': np.transpose([data[0], data[1]]).tolist(), 'line_caps': 'square',
+        })
+        series.append({
+            'type': 'series.line', 'id': f'line{index * 2 + 1}', 'name': f'line{index * 2 + 1}', 'color': colors[index],
+            'data': np.transpose([data[0], data[2]]).tolist(), 'line_caps': 'square',
+        })
+        series.append({
+            'type': 'text', 'id': f'text{index * 2 + 0}', 'name': f'text{index * 2 + 0}', 'color': colors[index],
+            'text': f'{smp.name()}<r>{round(smp.Info.results.age_plateau[0]["age"], 2)}',
+            'size': 10, 'data': [[index * 15 + 5, 23]],
+        })
+    data = {
+        "data": [
+            {
+                'xAxis': [{'extent': [0, 100], 'interval': [0, 20, 40, 60, 80, 100],
+                           'title': 'Cumulative <sup>39</sup>Ar Released (%)', 'nameLocation': 'middle', }],
+                'yAxis': [{'extent': [0, 250], 'interval': [0, 50, 100, 150, 200, 250],
+                           'title': 'Apparent Age (Ma)', 'nameLocation': 'middle', }],
+                'series': series
+            }
+        ],
+        "file_name": "WHA",
+        "plot_names": ["all age plateaus"],
+    }
+
+    # write pdf
+    file = pm.NewPDF(filepath=os.path.join(example_dir, f"{data['file_name']}.pdf"))
+    for index, each in enumerate(data['data']):
+        # rich text tags should follow this priority: color > script > break
+        file.text(page=index, x=50, y=780, line_space=1.2, size=12, base=0, h_align="left",
+                  text=f"The PDF can be edited with Adobe Acrobat, Illustrator and CorelDRAW")
+        cv = ap.smp.export.export_chart_to_pdf(each)
+        file.canvas(page=index, base=0, margin_top=5, canvas=cv, unit="cm", h_align="middle")
+        if index + 1 < len(data['data']):
+            file.add_page()
+
+    # save pdf
+    file.save()
+
+
 if __name__ == "__main__":
     test()
+    # export_pdf_demo()

{ararpy-0.1.11.dist-info → ararpy-0.1.12.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: ararpy
-Version: 0.1.11
+Version: 0.1.12
 Summary: A project for Ar-Ar geochronology
 Home-page: https://github.com/wuyangchn/ararpy.git
 Author: Yang Wu

{ararpy-0.1.11.dist-info → ararpy-0.1.12.dist-info}/RECORD

@@ -1,13 +1,13 @@
 ararpy/__init__.py,sha256=XKlsgklkzGfvqQkt59y-UST9aVhgCf-eXpEDiaZaynQ,6725
-ararpy/test.py,sha256=jdHhBvMbvDlxi0o2eRVG99nsUBNDA-0s2bbGbnXweo8,1130
+ararpy/test.py,sha256=4F46-JJ1Ge12HGae0qO44Qc6kiEMHBgn2MsY_5LlHDo,3973
 ararpy/calc/__init__.py,sha256=kUjRuLE8TLuKOv3i976RnGJoEMj23QBZDu37LWs81U4,322
 ararpy/calc/age.py,sha256=CUex9UxZbYinl9ypd7o86R3VRQFzahku3gW__fBeoss,5811
-ararpy/calc/arr.py,sha256=ntu6p3lUcF654SxLnFNrToLubXRWl101LRNFVJFulTQ,14871
+ararpy/calc/arr.py,sha256=P0DdYMXXnN69tNGL9NpB0hftNN03EEB_KFmfmZdIy-A,14922
 ararpy/calc/basic.py,sha256=wPzUXmrHAXdFkvd53JsH9Vae-o7O3LxYWztGQ5OpDhw,2871
-ararpy/calc/corr.py,sha256=n-mz-Uw-oC7ig6kSdbcq1Uca22NY21PeXgHY9Z-MZY4,18686
+ararpy/calc/corr.py,sha256=mQd5_CfDP2aBQrkQHKTV5R-zKbNOjXh8a21zHNrRNmU,19102
 ararpy/calc/err.py,sha256=63LtprqjemlIb1QGDst4Ggcv5KMSDHdlAIL-nyQs1eA,2691
 ararpy/calc/histogram.py,sha256=0GVbDdsjd91KQ1sa2B7NtZ4KGo0XpRIJapgIrzAwQUo,5777
-ararpy/calc/isochron.py,sha256=wJOiNKftwi8ipvNVanEcTfEgHV6slsZ1RL-y4v3rvq4,5610
+ararpy/calc/isochron.py,sha256=ej9G2e68k6yszonWHsLcEubh3TA7eh1upTJP_X0ttAA,5726
 ararpy/calc/jvalue.py,sha256=zHUhJ1iYe5mPrY95mGYxoUPAp7hwu4coUgiHKiruKfM,1138
 ararpy/calc/plot.py,sha256=v-vJwciJIhKrdXKWAUZbYslez0zBrIa6awI-DopmRZg,1956
 ararpy/calc/raw_funcs.py,sha256=fhAmqpvTvWUegbLi2PNbFRIbKgbDmycDIWNic3z25hU,2536
@@ -29,6 +29,7 @@ ararpy/examples/NGX-XLS.input-filter,sha256=iyfTLsLHf0h58ePOsNQoBtG6UBLIf7V-zTGi
 ararpy/examples/Qtegra-exported-xls.input-filter,sha256=KuD3Dey82ecuKnqNMQqOcY3VXHgKbpE9Yxfl1ZC4_so,441
 ararpy/examples/S01-239.csv,sha256=J_PHT85XCxPd3TEGzBARGwKSkiAzz1nzNbrKUf00GBU,39566
 ararpy/examples/WH01.irra,sha256=Ws78m5p3fD4oybtkOWoAtvFyuLPxpL_A_ueIg2tnhnw,365
+ararpy/examples/WHA.pdf,sha256=iopbptHkWPmYdgihEVDDfv2nm2XE-Q-e7fFnkR44Xh0,178692
 ararpy/examples/raw_example.xls,sha256=ftcSiXRx_7nYnbqJVma1Yl3Yr_iuceAWlEjhJwlAvFM,1929895
 ararpy/examples/sample-default.smp,sha256=YNkoQGgPrsL_fXS7ZHxfRtLQWekCDqT9czS6vBScImk,432
 ararpy/files/__init__.py,sha256=l5B5ZQ01WdtvjjN0aMkyAFNgpwANdM_1I0tQbqnRuEY,69
@@ -39,22 +40,22 @@ ararpy/files/new_file.py,sha256=efblARIBROVLWS2w3-98BxLX5VZ8grRpiTkJFtf_rAk,214
 ararpy/files/raw_file.py,sha256=B2k4RjAp4EViwHKv6CbI5400Liaa7OW_eBjGSPlN9Ek,22768
 ararpy/files/xls.py,sha256=8ibT4ZRY2grK34ikKm1pDmlWFlVTgDL7vSMO6mphzrY,702
 ararpy/smp/__init__.py,sha256=k6_fa27UJsQK7K7oC5GYlwMo6l0Xd8af3QtOrZz2XJk,478
-ararpy/smp/basic.py,sha256=1OgPXk-2OysJOZ9e_HuqB__fp3TNm844y0KKSlF5Wz0,21405
-ararpy/smp/calculation.py,sha256=JXo6vS-HcBhOe7b8OMLds-Sn9vofFjV6maY5u3E-X3c,2793
+ararpy/smp/basic.py,sha256=DrwU4ApXOkxYNMwax7B-ErBndJnUwG5n64fLU6zR63Y,21506
+ararpy/smp/calculation.py,sha256=1dUKYbq83j0YpZ9Nbb5rM0wbaEXezDKqcbUCc9TRKng,2919
 ararpy/smp/consts.py,sha256=XIdjdz8cYxspG2jMnoItdlUsxr3hKbNFJjMZJh1bpzw,393
-ararpy/smp/corr.py,sha256=3lvs-2A3WjK_yfRnQRCG1Tb_QNhI4j82HU6a5OAN4Q8,23402
-ararpy/smp/diffusion_funcs.py,sha256=u7cEH7sac5tw0Z9Yvh7AuH93ZYEENJVtZ3Npni8dWUE,167616
-ararpy/smp/export.py,sha256=MBcisVfPwmdMyYqE7Tm3oLAMsVA5Rk0fl3Tqn0LxouI,78231
+ararpy/smp/corr.py,sha256=MXF04dLYGcZN7JHFq8OJqawkZ-yUjN45qG0NJX_kBB0,24528
+ararpy/smp/diffusion_funcs.py,sha256=5yDlyVVcaKQDMPLd7ZvyCetyiQNBR57XtrJkLT_Hsww,168544
+ararpy/smp/export.py,sha256=SEUySwBnDm6ao0Ujc93y98FamUx_UmEsqnd3cmFgKPA,80817
 ararpy/smp/info.py,sha256=iKUELm-BuUduDlJKC1d8tKKNHbwwbNmhUg2pi6bcBvA,489
 ararpy/smp/initial.py,sha256=nZNrV4fv_PlNXYfqqIOOOkcQjnSOTLymVpJyJ9d_QtM,15359
-ararpy/smp/json.py,sha256=87zMMopD5eL3VsemDYTEAfhSWXO_1Z7Y782S4GR2UUs,1740
-ararpy/smp/plots.py,sha256=s6xDZyxeGl512Wn4knquhzFeECucZCWNAsfI70OInqI,31362
+ararpy/smp/json.py,sha256=Y7xoTZczSTrvVHlhfLmpu-zoxgLw9MFn2lgRYkjgHTY,2079
+ararpy/smp/plots.py,sha256=bM4ofS12NedBUXXleEy7-xcK9C8svVuSrTXYmeiMY3k,31495
 ararpy/smp/raw.py,sha256=KtAmzTAwCDK4x1dRgvwj2TbOzHRv7LFwaU2hHCknF3c,6173
-ararpy/smp/sample.py,sha256=S1vgt0LRihcAiqpXhuLLJnkft3kVp0QwzjH_JQkMyc0,54611
+ararpy/smp/sample.py,sha256=KPysezU7avqUBu6ykIEtP3_2wa5a16n0gC0iF_kqSoU,54988
 ararpy/smp/style.py,sha256=NZwTeodhtHBBdCyjvupF64nh2B5EDdMHzlT8oE360zg,6746
-ararpy/smp/table.py,sha256=fyITVbj0JrPoB4SPW4dHoJl2jliKEZ_HpS2FSF9FQx4,6274
-ararpy-0.1.11.dist-info/LICENSE,sha256=cvG5t_C1qY_zUyJI7sNOa7gCArdngNPaOrfujl2LYuc,1085
-ararpy-0.1.11.dist-info/METADATA,sha256=uLONl17cQulcXUYmJ4BvhpxHyYrhIOs3i8Z9uOb6zD0,24331
-ararpy-0.1.11.dist-info/WHEEL,sha256=GV9aMThwP_4oNCtvEC2ec3qUYutgWeAzklro_0m4WJQ,91
-ararpy-0.1.11.dist-info/top_level.txt,sha256=9iTpsPCYuRYq09yQTk9d2lqB8JtTEOmbN-IcGB-K3vY,7
-ararpy-0.1.11.dist-info/RECORD,,
+ararpy/smp/table.py,sha256=yNr4PLxbH4ezeUXpKBxLjRfq4RuYiW-G3sMWFDtjlVM,6802
+ararpy-0.1.12.dist-info/LICENSE,sha256=cvG5t_C1qY_zUyJI7sNOa7gCArdngNPaOrfujl2LYuc,1085
+ararpy-0.1.12.dist-info/METADATA,sha256=BiPSBAkyiLpDd-5VsJm14TqeH58VkFf1G1mpQh3Xx4E,24331
+ararpy-0.1.12.dist-info/WHEEL,sha256=GV9aMThwP_4oNCtvEC2ec3qUYutgWeAzklro_0m4WJQ,91
+ararpy-0.1.12.dist-info/top_level.txt,sha256=9iTpsPCYuRYq09yQTk9d2lqB8JtTEOmbN-IcGB-K3vY,7
+ararpy-0.1.12.dist-info/RECORD,,