asteroid_spinprops 0.2.30__tar.gz → 0.2.32__tar.gz

{asteroid_spinprops-0.2.30 → asteroid_spinprops-0.2.32}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: asteroid_spinprops
-Version: 0.2.30
+Version: 0.2.32
 Summary: Collection of tools used for fitting sHG1G2 and ssHG1G2 photometric models to sparse asteroid photometry
 License: MIT
 Author: Odysseas

{asteroid_spinprops-0.2.30 → asteroid_spinprops-0.2.32}/asteroid_spinprops/ssolib/periodest.py

@@ -194,9 +194,29 @@ def get_multiband_period_estimate(
 
     period_in = 2 * (1 / f_chosen)
 
+    model_final = LombScargleMultiband(
+        residuals_dataframe["jd"].values,
+        residuals_dataframe["residuals"].values,
+        residuals_dataframe["filters"].values,
+        residuals_dataframe["sigma"].values,
+        normalization="standard",
+        fit_mean=True,
+        nterms_base=k_val,
+        nterms_band=1,
+    )
+    frequency_final, power_final = model_final.autopower(
+        method="fast",
+        sb_method="fastnifty_chi2",
+        minimum_frequency=1 / period_range[1],
+        maximum_frequency=1 / period_range[0],
+        samples_per_peak=5,
+    )
+
     ls = LombScargle(
         residuals_dataframe["jd"].values,
         np.ones(len(residuals_dataframe["jd"])),
+        fit_mean=False,
+        center_data=False,
     )
     frequencyW, powerW = ls.autopower(
         minimum_frequency=1 / period_range[1],
@@ -205,7 +225,7 @@ def get_multiband_period_estimate(
     )
 
     pindex, heights = find_peaks(
-        power,
+        power_final,
         height=0.1,
         threshold=None,
         distance=260,
@@ -219,7 +239,7 @@ def get_multiband_period_estimate(
     hindex = np.argsort(heights["peak_heights"])[::-1][:5]
     hindexW = np.argsort(heightsW["peak_heights"])[::-1][:5]
 
-    signal_peaks = frequency[pindex[hindex]]
+    signal_peaks = frequency_final[pindex[hindex]]
     window_peaks = frequencyW[pindexW[hindexW]]
 
     return period_in, k_val, p_rms, signal_peaks, window_peaks

{asteroid_spinprops-0.2.30 → asteroid_spinprops-0.2.32}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "asteroid_spinprops"
-version = "0.2.30"
+version = "0.2.32"
 description = "Collection of tools used for fitting sHG1G2 and ssHG1G2 photometric models to sparse asteroid photometry"
 authors = [
     {name = "Odysseas",email = "odysseas.xenos@proton.me"}