Robin-Rhee
diff --git a/‎v1/src/sensing/feature_extractor.py‎
Lines changed: 21 additions & 2 deletions b/‎v1/src/sensing/feature_extractor.py‎
Lines changed: 21 additions & 2 deletions
@@ -103,6 +103,8 @@ def extract(self, samples: List[WifiSample]) -> RssiFeatures:
 
         # Trim to window
         samples = self._trim_to_window(samples)
+        if len(samples) < 4:
+            return RssiFeatures(n_samples=len(samples))
         rssi = np.array([s.rssi_dbm for s in samples], dtype=np.float64)
         timestamps = np.array([s.timestamp for s in samples], dtype=np.float64)
 
@@ -158,17 +160,34 @@ def extract_from_array(
 
         return features
 
+    # -- window trimming -----------------------------------------------------
+
+    def _trim_to_window(self, samples: List[WifiSample]) -> List[WifiSample]:
+        """Keep only samples within the most recent ``window_seconds``."""
+        if not samples:
+            return samples
+        latest_ts = samples[-1].timestamp
+        cutoff = latest_ts - self._window_seconds
+        trimmed = [s for s in samples if s.timestamp >= cutoff]
+        return trimmed
+
     # -- time-domain ---------------------------------------------------------
 
     @staticmethod
     def _compute_time_domain(rssi: NDArray[np.float64], features: RssiFeatures) -> None:
         features.mean = float(np.mean(rssi))
         features.variance = float(np.var(rssi, ddof=1)) if len(rssi) > 1 else 0.0
         features.std = float(np.std(rssi, ddof=1)) if len(rssi) > 1 else 0.0
-        features.skewness = float(scipy_stats.skew(rssi, bias=False)) if len(rssi) > 2 else 0.0
-        features.kurtosis = float(scipy_stats.kurtosis(rssi, bias=False)) if len(rssi) > 3 else 0.0
         features.range = float(np.ptp(rssi))
 
+        # Guard against constant signals where higher moments are undefined
+        if features.std < 1e-12:
+            features.skewness = 0.0
+            features.kurtosis = 0.0
+        else:
+            features.skewness = float(scipy_stats.skew(rssi, bias=False)) if len(rssi) > 2 else 0.0
+            features.kurtosis = float(scipy_stats.kurtosis(rssi, bias=False)) if len(rssi) > 3 else 0.0
+
         q75, q25 = np.percentile(rssi, [75, 25])
         features.iqr = float(q75 - q25)