1. 項目概述從一道數學建模題說起幾年前我接手指導一個學生團隊參加全國大學生數學建模競賽CUMCM他們抽到的題目是2011年的A題關于城市表層土壤重金屬污染分析。拿到數據的第一刻團隊里編程最好的同學就懵了——Excel里密密麻麻的采樣點數據夾雜著大量明顯不合理甚至為負數的重金屬濃度值。直接拿這樣的數據去擬合、做空間插值結果必然是失真的甚至可能得出“某區域土壤重金屬含量為負”這種荒謬結論。這就是異常值數據分析路上第一只也是最常見的“攔路虎”。這個項目我們就以這個真實的數學建模賽題為背景用Python這把“手術刀”來系統性地處理數據中的異常值。這不僅僅是刪除幾個數字那么簡單它關乎你后續所有模型的可信度。無論是數學建模、商業分析還是科研數據處理清洗數據尤其是識別并妥善處理異常值是決定項目成敗的第一步。如果你正在為數據中的“噪點”頭疼或者想系統學習Python在數據清洗中的實戰技巧那么這篇基于真實項目復盤的經驗分享就是為你準備的。我們將從原理到實操一步步拆解如何用Pandas、NumPy、Matplotlib等工具像偵探一樣找出數據中的“異類”并科學地決定它們的去留。2. 核心思路為什么不能簡單“一刪了之”面對異常值新手最容易犯的錯誤就是直接刪除。在2011年A題的數據里看到“砷(As)”濃度出現一個“-999”或者“0.001”這種與其他樣本相差幾個數量級的值第一反應可能就是把它當成錯誤數據扔掉。但且慢在動刀之前我們必須先回答三個問題它真的是“異常”嗎如果是它為什么會出現我們應該如何處理它2.1 異常值的本質與分類異常值或稱離群點是指與數據集中其他觀測值顯著不同的數據點。這種“不同”可能源于數據錄入或測量錯誤這是最理想的情況比如儀器故障、記錄筆誤多輸了一個0、單位混淆把mg/kg記成了g/kg。2011年A題中部分負值濃度很可能就屬于此類。采樣誤差或實驗變異例如采樣點恰好位于一個污染源排放口其濃度自然遠高于背景值。這雖然是“異常”的但卻是真實的、有價值的信息直接刪除會掩蓋重要的污染特征。數據本身的自然分布某些數據本身就服從重尾分布如金融收益、某些環境污染物濃度存在極端值是正常的概率事件。因此我們的處理策略絕不能是“一刀切”。處理異常值的核心思路是一個決策流程檢測 - 診斷 - 處理。檢測是用統計或可視化方法找出嫌疑點診斷是結合業務背景在數學建模中就是題目背景和專業知識判斷嫌疑點的性質處理則是根據診斷結果選擇保留、修正、刪除或用穩健方法替代。2.2 數學建模場景下的特殊考量在數學建模競賽中處理異常值還需額外考慮兩點可解釋性你的處理方法必須在論文中清晰闡述并給出令人信服的理由。直接說“我們刪除了大于3倍標準差的數據”是蒼白的必須結合題目背景如土壤環境質量標準、重金屬的地球化學背景值來論證閾值設定的合理性。數據完整性競賽數據通常有限每一個樣本都可能攜帶關鍵信息。盲目刪除可能導致樣本量過小影響后續建模如插值生成污染分布圖時邊界出現空洞。因此修正或插補往往是比刪除更優的選擇。基于以上思路我們為本次實戰設計了一個四步流程首先進行探索性數據分析EDA可視化異常其次運用統計方法進行定量檢測然后結合背景知識進行診斷與閾值確定最后根據診斷結果選擇適當的處理策略并實現。下面我們就用Python一步步實現它。3. 環境準備與數據加載工欲善其事必先利其器。我們首先搭建一個用于數據清洗的Python環境。我個人強烈推薦使用Anaconda來管理環境它能很好地解決包依賴問題。3.1 創建專屬分析環境打開你的終端或Anaconda Prompt執行以下命令來創建一個純凈的、包含所有必要庫的環境。# 創建一個名為data_cleaning的新環境指定Python版本為3.9兼容性好 conda create -n data_cleaning python3.9 # 激活該環境 conda activate data_cleaning # 安裝核心數據分析庫 conda install pandas numpy matplotlib seaborn scipy jupyter -y # 安裝用于更高級統計和插值的庫可選但推薦 conda install scikit-learn statsmodels -y注意如果你沒有安裝Anaconda也可以直接使用pip install命令安裝上述包。但使用虛擬環境venv或conda env是一個好習慣可以避免不同項目間的包版本沖突。3.2 加載與初窺數據由于2011年A題的原始數據需要從賽題頁面下載這里我模擬生成了一份結構、特征相似的合成數據集以便完整演示流程。假設我們有一個soil_heavy_metal.csv文件包含采樣點ID、坐標X Y以及8種重金屬As, Cd, Cr, Cu, Hg, Ni, Pb, Zn的濃度數據。import pandas as pd import numpy as np import matplotlib.pyplot as plt import seaborn as sns # 設置繪圖風格讓圖表更好看 plt.style.use(seaborn-v0_8-darkgrid) sns.set_palette(husl) # 1. 加載數據 df pd.read_csv(soil_heavy_metal.csv) print(數據形狀行列:, df.shape) print(\n前5行數據) print(df.head()) print(\n數據基本信息) print(df.info()) print(\n描述性統計) print(df.describe())運行這段代碼你會立刻對數據有個整體印象有多少行、多少列、是否有缺失值、各重金屬濃度的均值、標準差、最小最大值等。df.describe()輸出的最大值、最小值往往是發現異常值的第一條線索。比如你看到某元素的最小值是負數或者最大值是均值的上千倍紅燈就亮起了。3.3 數據質量初步檢查在深入分析前先做一次快速體檢。# 檢查缺失值 missing_sum df.isnull().sum() print(各列缺失值數量\n, missing_sum[missing_sum 0]) if missing_sum.sum() 0: print(數據集中暫無缺失值。) # 檢查重復行 duplicate_rows df.duplicated().sum() print(f\n重復行數量{duplicate_rows}) # 檢查明顯不合理的值例如濃度為負 negative_cols [] for col in df.columns: if df[col].dtype in [int64, float64]: # 只檢查數值列 if (df[col] 0).any(): negative_cols.append(col) print(f\n存在負值的列{negative_cols}) if negative_cols: for col in negative_cols: neg_count (df[col] 0).sum() print(f {col}: {neg_count} 個負值例如{df.loc[df[col] 0, col].iloc[:3].values})這一步能幫你快速抓住最明顯的“硬傷”比如錄入錯誤導致的負值。對于這些明顯的錯誤我們可以在后續步驟中優先處理。4. 異常值檢測可視化與統計方法雙管齊下檢測異常值我習慣“先看圖再算數”。可視化能給你直觀感受統計方法則提供定量依據。4.1 可視化檢測讓異常值無所遁形4.1.1 箱線圖Boxplot箱線圖是識別異常值的經典工具。它顯示了數據的中位數、上下四分位數Q1, Q3以及“須”的范圍。通常落在Q1 - 1.5IQR四分位距以下或Q3 1.5IQR以上的點被視為潛在的異常值在箱線圖中常以圓點形式標出。# 繪制所有重金屬濃度的箱線圖 heavy_metal_cols [As, Cd, Cr, Cu, Hg, Ni, Pb, Zn] plt.figure(figsize(14, 8)) df[heavy_metal_cols].boxplot() plt.title(重金屬濃度箱線圖初步識別異常值) plt.ylabel(濃度 (mg/kg)) plt.xticks(rotation45) plt.tight_layout() plt.show()通過箱線圖你可以一眼看出哪些元素的濃度分布存在大量“飛點”以及這些異常值是偏于高值還是低值。例如Cd鎘可能顯示出較多的高位異常點這很可能對應著局部的強污染源。4.1.2 直方圖與Q-Q圖直方圖看分布形態Q-Q圖分位數-分位數圖檢驗數據是否服從正態分布。許多統計檢測方法如3σ原則基于正態假設如果數據嚴重偏離正態這些方法可能失效。fig, axes plt.subplots(2, 4, figsize(16, 10)) axes axes.ravel() for idx, col in enumerate(heavy_metal_cols): # 直方圖 axes[idx].hist(df[col].dropna(), bins30, edgecolorblack, alpha0.7) axes[idx].set_title(f{col} 分布) axes[idx].set_xlabel(濃度) axes[idx].set_ylabel(頻數) plt.suptitle(重金屬濃度直方圖, fontsize16) plt.tight_layout() plt.show() # 以As為例繪制Q-Q圖 import scipy.stats as stats plt.figure(figsize(6,6)) stats.probplot(df[As].dropna(), distnorm, plotplt) plt.title(As濃度Q-Q圖檢驗正態性) plt.show()如果直方圖嚴重右偏長尾在右且Q-Q圖上的點明顯偏離對角線說明數據非正態。這時你需要考慮使用基于中位數和四分位數的穩健方法如箱線圖法而不是基于均值和標準差的方法。4.1.3 散點圖矩陣觀察多個變量之間的關系以及是否存在遠離主體群的觀測點。# 選擇幾個關鍵元素繪制散點圖矩陣 sns.pairplot(df[[As, Cd, Pb, Zn]], diag_kindkde, cornerTrue) plt.suptitle(關鍵重金屬元素散點圖矩陣, y1.02) plt.show()有時單變量看不是異常但在多變量空間里卻是異常比如As和Cd濃度都異常高。散點圖能幫你發現這種多維異常。4.2 統計方法檢測設定定量閾值可視化給了我們直覺接下來需要用數字說話。4.2.1 標準差法3σ原則適用于近似正態分布的數據。認為99.7%的數據落在均值上下3個標準差的范圍內之外的即為異常值。def detect_outliers_std(df, column, n_std3): 使用標準差法檢測異常值 mean df[column].mean() std df[column].std() lower_bound mean - n_std * std upper_bound mean n_std * std outliers df[(df[column] lower_bound) | (df[column] upper_bound)] return outliers, lower_bound, upper_bound # 檢測As的異常值 as_outliers_std, low_std, up_std detect_outliers_std(df, As, n_std3) print(fAs濃度 - 標準差法±3σ:) print(f 邊界[{low_std:.2f}, {up_std:.2f}]) print(f 異常值數量{len(as_outliers_std)}) print(f 異常值索引示例{as_outliers_std.index[:5].tolist()})4.2.2 箱線圖法IQR法這是一種更穩健的方法不依賴于正態分布假設對極端值不敏感。def detect_outliers_iqr(df, column): 使用IQR法檢測異常值 Q1 df[column].quantile(0.25) Q3 df[column].quantile(0.75) IQR Q3 - Q1 lower_bound Q1 - 1.5 * IQR upper_bound Q3 1.5 * IQR outliers df[(df[column] lower_bound) | (df[column] upper_bound)] return outliers, lower_bound, upper_bound # 檢測As的異常值IQR法 as_outliers_iqr, low_iqr, up_iqr detect_outliers_iqr(df, As) print(f\nAs濃度 - IQR法1.5IQR:) print(f 邊界[{low_iqr:.2f}, {up_iqr:.2f}]) print(f 異常值數量{len(as_outliers_iqr)})實操心得對比兩種方法的結果。對于嚴重偏態的數據標準差法可能會誤將許多正常數據判為異常因為均值被極端值拉偏而IQR法通常更可靠。在數學建模論文中我建議同時報告兩種方法的結果并解釋為什么最終選擇其中一種例如“由于數據呈右偏分布我們采用對極端值不敏感的IQR法進行異常值識別”。4.2.3 基于模型的方法孤立森林對于更復雜的數據集可以使用機器學習算法。孤立森林通過隨機劃分特征空間來隔離樣本異常點因為“與眾不同”而容易被快速隔離。from sklearn.ensemble import IsolationForest # 假設我們針對所有重金屬濃度特征進行檢測 X df[heavy_metal_cols].fillna(df[heavy_metal_cols].median()) # 用中位數填充缺失值 # 初始化并訓練孤立森林模型 # contamination參數估計異常值比例可根據箱線圖結果大致設定 iso_forest IsolationForest(contamination0.05, random_state42) outlier_labels iso_forest.fit_predict(X) # 標記異常值-1表示異常 df[iso_forest_outlier] outlier_labels outlier_count (df[iso_forest_outlier] -1).sum() print(f\n孤立森林檢測出的異常樣本數{outlier_count})孤立森林能捕捉多維異常但解釋性稍差。在數學建模中可以作為輔助驗證手段不建議作為主要判定依據因為論文中需要清晰可解釋的閾值。5. 異常值診斷與處理策略制定檢測出的“嫌疑點”名單有了現在進入關鍵的診斷與決策階段。這一步需要結合題目背景知識。5.1 診斷它是“錯誤”還是“特征”以2011年A題為例我們需要查閱《土壤環境質量 農用地土壤污染風險管控標準試行》GB 15618-2018等資料了解各重金屬的背景值范圍和風險篩選值。例如某區域土壤砷的背景值可能在5-15 mg/kg如果檢測到一個200 mg/kg的點它很可能是一個真實的高污染點特征而非錯誤。反之一個-1 mg/kg的值顯然是錯誤。診斷流程對照物理/化學極限濃度不能為負通常有檢測下限。負值和低于檢測下限的異常低值可初步判定為錯誤。對照背景值或標準值查閱文獻或國家標準了解該區域該元素的正常范圍。遠超正常范圍上限的可能是污染點遠低于下限的可能是分析誤差或特殊地質背景。空間關聯分析在地圖上標出異常點。如果高值點聚集在工業區、交通干線附近則很可能是真實污染如果孤立且遠離任何潛在源則錯誤可能性增大。元素比值分析某些重金屬之間存在穩定的地球化學比值。計算異常點的元素比值如Pb/Zn Cu/Ni看是否偏離區域背景比值輔助判斷。# 示例結合空間位置分析異常值 # 假設df中有X, Y坐標列 as_outliers, _, _ detect_outliers_iqr(df, As) plt.figure(figsize(10, 8)) # 繪制所有采樣點 plt.scatter(df[X], df[Y], cblue, alpha0.5, label正常點, s20) # 高亮標注As異常點 plt.scatter(as_outliers[X], as_outliers[Y], cred, marker^, s100, labelAs異常點 (IQR法)) plt.xlabel(X坐標) plt.ylabel(Y坐標) plt.title(As濃度異常點空間分布) plt.legend() plt.grid(True, alpha0.3) plt.tight_layout() plt.show()5.2 處理策略選擇根據診斷結果選擇處理方式診斷結果可能原因推薦處理策略Python實現思路明顯錯誤錄入錯誤、儀器故障、負值、低于檢測限的零值刪除或修正直接刪除行或用np.nan標記后插補真實異常污染點局部強污染源保留但單獨分析保留數據在后續建模中可作為重點研究對象或使用穩健統計方法如中位數真實異常自然極端數據本身重尾分布保留或穩健處理使用對數變換、Box-Cox變換使數據更接近正態或使用非參數方法不確定信息不足無法判斷多重插補或標記用多重插補法生成多個可能值或創建布爾列標記為“可疑點”在數學建模中的建議對于明顯錯誤在論文中明確說明“我們發現了N個負值/零值數據點這些在物理意義上不可能判定為錄入錯誤予以刪除。”對于疑似真實高值不要輕易刪除可以這樣處理“識別出M個顯著高于背景值的采樣點如As X mg/kg。經查這些點均位于歷史工業區附近我們認為是真實污染予以保留。在后續空間插值時我們采用了穩健的克里金方法以降低這些極端值對整體插值面的過度影響。”穩健性處理在計算區域平均濃度、進行相關性分析時使用中位數和Spearman秩相關代替均值和Pearson相關它們對異常值不敏感。6. 異常值處理實操Python代碼實現現在我們將診斷決策轉化為具體的Python操作。假設我們對As列的診斷結果是負值均為錯誤予以刪除部分極高值 Q3 3IQR為真實污染保留但標記其他IQR法檢測出的高值用中位數插補。6.1 處理明顯錯誤如負值# 記錄原始數據形狀 original_shape df.shape print(f原始數據形狀{original_shape}) # 1. 處理負值直接刪除或設為NaN # 假設我們判定所有負值為錯誤 for col in heavy_metal_cols: # 首先查看負值情況 neg_mask df[col] 0 if neg_mask.any(): print(f在 {col} 列中發現 {neg_mask.sum()} 個負值正在刪除對應行...) # 刪除任何重金屬濃度為負的整行數據因為一個點采樣所有元素數據應同進退 df df[~neg_mask] print(f刪除負值后數據形狀{df.shape}) print(f共刪除 {original_shape[0] - df.shape[0]} 行。)6.2 處理疑似異常高值分情況處理# 2. 識別并分情況處理高異常值 # 2.1 使用更嚴格的IQR邊界如3倍IQR識別“極端高值”我們認為是真實污染保留但標記。 def get_extreme_high_outliers(df, column, iqr_multiplier3): Q1 df[column].quantile(0.25) Q3 df[column].quantile(0.75) IQR Q3 - Q1 upper_bound_extreme Q3 iqr_multiplier * IQR extreme_outliers df[df[column] upper_bound_extreme] return extreme_outliers, upper_bound_extreme # 以As為例 as_extreme, as_extreme_bound get_extreme_high_outliers(df, As, iqr_multiplier3) print(f\nAs濃度極端高值 Q33IQR {as_extreme_bound:.2f}數量{len(as_extreme)}) print(這些點將被保留并標記為‘真實污染點’。) # 在數據框中添加標記列 df[As_extreme_high] df[As] as_extreme_bound # 2.2 對于普通高異常值Q31.5IQR value Q33IQR我們認為是可能的誤差或輕度異常用中位數插補。 as_normal_outliers, low_bound, up_bound detect_outliers_iqr(df, As]) # 1.5IQR邊界 # 注意這里的as_normal_outliers包含了極端高值我們需要排除它們 as_normal_outliers as_normal_outliers[~as_normal_outliers[As_extreme_high]] if not as_normal_outliers.empty: print(f\nAs濃度普通高異常值1.5IQR value 3IQR數量{len(as_normal_outliers)}) print(這些點將用該列中位數進行插補。) # 計算中位數排除極端高值避免被它們拉高 median_as df.loc[~df[As_extreme_high], As].median() # 記錄原始值便于追溯 df[As_original] df[As].copy() # 進行插補 df.loc[as_normal_outliers.index, As] median_as print(f用于插補的中位數為{median_as:.4f}) else: print(\n未發現需要插補的普通高異常值。)6.3 處理其他元素的異常值我們可以將上述邏輯封裝成函數批量處理所有重金屬列。def process_column_outliers(df, column_name, extreme_multiplier3, normal_multiplier1.5): 處理單列異常值 1. 負值已在前一步刪除此處不考慮。 2. 識別極端高值 Q3extreme_multiplier*IQR保留并標記。 3. 識別普通高異常值Q3normal_multiplier*IQR value Q3extreme_multiplier*IQR用中位數插補。 # 創建標記列 extreme_flag_col f{column_name}_extreme_high original_col f{column_name}_original # 計算邊界 Q1 df[column_name].quantile(0.25) Q3 df[column_name].quantile(0.75) IQR Q3 - Q1 upper_bound_extreme Q3 extreme_multiplier * IQR upper_bound_normal Q3 normal_multiplier * IQR # 標記極端高值 df[extreme_flag_col] df[column_name] upper_bound_extreme # 識別普通高異常值在正常邊界和極端邊界之間且未被標記為極端 normal_high_mask (df[column_name] upper_bound_normal) (df[column_name] upper_bound_extreme) (~df[extreme_flag_col]) # 保存原始值 if original_col not in df.columns: df[original_col] df[column_name] # 插補普通高異常值 if normal_high_mask.any(): # 使用非極端值的中位數進行插補 median_val df.loc[~df[extreme_flag_col], column_name].median() df.loc[normal_high_mask, column_name] median_val print(f {column_name}: 插補了 {normal_high_mask.sum()} 個普通高異常值使用中位數 {median_val:.4f}) else: print(f {column_name}: 無普通高異常值需插補。) # 返回極端高值數量和信息 extreme_count df[extreme_flag_col].sum() return extreme_count, upper_bound_extreme print(\n--- 開始批量處理各重金屬列 ---) for col in heavy_metal_cols: if col in df.columns: # 確保列還存在未被刪除 extreme_cnt, bound process_column_outliers(df, col, extreme_multiplier3, normal_multiplier1.5) print(f {col}: 保留了 {extreme_cnt} 個極端高值 {bound:.2f}視為潛在污染點。)6.4 處理后的數據驗證與保存處理完成后務必驗證數據質量并保存清洗后的數據。# 驗證再次繪制箱線圖觀察異常值處理效果 plt.figure(figsize(14, 8)) df[heavy_metal_cols].boxplot() plt.title(異常值處理后重金屬濃度箱線圖) plt.ylabel(濃度 (mg/kg)) plt.xticks(rotation45) plt.tight_layout() plt.show() # 查看處理摘要 print(\n 數據處理摘要 ) print(f最終數據形狀{df.shape}) print(\n各列極端高值真實污染點數量) extreme_cols [col for col in df.columns if _extreme_high in col] for col in extreme_cols: metal col.replace(_extreme_high, ) count df[col].sum() print(f {metal}: {count}) # 保存清洗后的數據 output_file soil_heavy_metal_cleaned.csv df.to_csv(output_file, indexFalse) print(f\n清洗后的數據已保存至{output_file})7. 高級話題與常見問題排查在實際操作中你肯定會遇到比教程更復雜的情況。下面分享一些進階技巧和踩坑經驗。7.1 缺失值、異常值與插補的協同處理很多時候數據集中同時存在缺失值NaN和異常值。處理順序很重要。我推薦的流程是標記異常值先用IQR或業務規則找出異常值但不急于修改或刪除只是新增一個布爾列進行標記。處理缺失值對缺失值進行插補。關鍵點插補時應排除那些被標記為異常值的樣本避免用異常值去影響插補模型。例如用KNN插補時特征矩陣中不應包含異常值樣本。處理異常值根據診斷對標記的異常值進行刪除、修正或保留操作。# 示例在插補缺失值時排除異常值 from sklearn.impute import KNNImputer # 假設我們已經有了‘As_is_extreme’標記列 # 1. 創建用于插補的數據副本并將極端異常值暫時設為NaN不參與建模 df_for_impute df[heavy_metal_cols].copy() for col in heavy_metal_cols: extreme_flag f{col}_extreme_high if extreme_flag in df.columns: df_for_impute.loc[df[extreme_flag], col] np.nan # 2. 使用KNN插補缺失值現在包括原始的缺失值和被我們暫時設為NaN的異常值位置 imputer KNNImputer(n_neighbors5) df_imputed_array imputer.fit_transform(df_for_impute) df_imputed pd.DataFrame(df_imputed_array, columnsheavy_metal_cols, indexdf.index) # 3. 將插補后的值填回原數據集僅填充原為NaN或我們暫時設為NaN的位置 for col in heavy_metal_cols: mask_original_missing df[col].isna() mask_extreme_set_nan df_for_impute[col].isna() ~mask_original_missing fill_mask mask_original_missing | mask_extreme_set_nan df.loc[fill_mask, col] df_imputed.loc[fill_mask, col]7.2 基于業務規則的定制化檢測統計方法有其局限性必須結合領域知識。例如在環境領域常用“背景值幾何標準差”的方法。def detect_outliers_geometric(df, column, log_base10, n_std2): 基于對數正態分布假設的異常值檢測。 適用于濃度等通常呈對數正態分布的環境數據。 # 取對數處理前確保沒有0或負值 data_log np.log10(df[column].clip(lower1e-10)) # clip避免0 geometric_mean np.power(10, data_log.mean()) geometric_std np.power(10, data_log.std()) lower_bound geometric_mean / (geometric_std ** n_std) upper_bound geometric_mean * (geometric_std ** n_std) outliers df[(df[column] lower_bound) | (df[column] upper_bound)] return outliers, lower_bound, upper_bound # 使用示例 as_outliers_geo, low_geo, up_geo detect_outliers_geometric(df, As, n_std2) print(f幾何法檢測2倍幾何標準差邊界[{low_geo:.2f}, {up_geo:.2f}]) print(f異常值數量{len(as_outliers_geo)})7.3 常見問題與排查技巧實錄問題1處理后的數據分布嚴重變形現象箱線圖顯示數據被“壓扁”在一個很窄的范圍內失去了變異信息。原因可能過于激進地刪除了大量“異常值”或者用中位數/均值插補了太多點。排查檢查異常值檢測的閾值是否過嚴如用了2σ而非3σ。回顧診斷步驟確認被處理的是否大部分是“真實錯誤”而非“真實特征”。解決放寬閾值或對“疑似真實異常”采用更保守的處理方式如Winsorizing縮尾法僅將極端值替換為指定分位數而非全部替換為中位數。def winsorize_column(series, limits(0.05, 0.05)): 縮尾處理將兩端極端值替換為指定分位數 from scipy.stats.mstats import winsorize # limits(0.05, 0.05) 表示將最低5%和最高5%的值縮尾 winsorized_data winsorize(series, limitslimits) return pd.Series(winsorized_data, indexseries.index) # 對As列進行5%縮尾處理 df[As_winsorized] winsorize_column(df[As], limits(0.05, 0.05))問題2多變量異常檢測結果與單變量不一致現象用孤立森林找出的異常點在單個變量的箱線圖上看起來并不極端。原因該樣本在所有單變量維度上都正常但多個變量的組合關系異常例如Cd濃度中等Pb濃度中等但兩者的比值異常高。排查檢查這些樣本在多維散點圖或降維圖如PCA得分圖上的位置。解決這類異常點可能揭示了特殊污染類型或數據質量問題需要結合空間位置和背景知識深入分析不能簡單忽略。問題3處理后后續建模如克里金插值在異常點位置出現“空洞”或“尖峰”現象刪除異常點后插值圖在該區域出現空白保留極端高值則產生不合理的“污染島”。原因空間插值方法對數據密度和極值敏感。解決對于刪除的點考慮使用考慮變量相關性的多重插補方法在刪除點位置生成合理的估計值再參與插值。對于保留的極端點在插值時使用穩健變異函數Robust Variogram估計它受極端值影響小。或者在插值前對數據進行對數轉換以緩和極端值的影響。# 使用Ordinary Kriging插值前進行對數轉換 df[As_log] np.log1p(df[As]) # log1p log(1x)避免x0的問題 # 或者使用穩健變異函數需要專業地統計學庫如pykrige或gstools # 這里以概念為例問題4自動化腳本在處理新數據集時效果不佳現象為2011年數據寫的處理流程用在2012年數據上要么漏掉很多異常值要么把正常值誤殺了。原因不同數據集的數據分布、量級、背景值可能不同。硬編碼的閾值如固定倍數IQR不具備普適性。解決將閾值參數化并編寫一個數據質量報告生成函數在新數據上先運行報告人工審查建議的閾值再進行調整。def generate_data_quality_report(df, numeric_cols): 生成數據質量報告輔助設定閾值 report {} for col in numeric_cols: desc df[col].describe(percentiles[.01, .05, .25, .5, .75, .95, .99]) Q1, Q3 desc[25%], desc[75%] IQR Q3 - Q1 report[col] { mean: desc[mean], std: desc[std], min: desc[min], 1%: desc[1%], 5%: desc[5%], median: desc[50%], 95%: desc[95%], 99%: desc[99%], max: desc[max], IQR: IQR, 建議箱線圖法上限 (Q31.5IQR): Q3 1.5 * IQR, 建議箱線圖法下限 (Q1-1.5IQR): Q1 - 1.5 * IQR, 偏度: df[col].skew(), 峰度: df[col].kurtosis() } report_df pd.DataFrame(report).T return report_df # 生成報告 quality_report generate_data_quality_report(df, heavy_metal_cols) print(quality_report[[min, 1%, 建議箱線圖法下限, median, 建議箱線圖法上限, 99%, max, 偏度]])這份報告能幫你快速了解每個變量的分布特別是“偏度”指標。如果偏度遠大于0右偏那么基于均值的標準差法就可能不適用報告中的百分位數1% 99%和箱線圖法建議邊界能給你更合理的調整起點。數據處理尤其是異常值處理從來不是一道有標準答案的數學題。它混合了統計學、領域知識和實際經驗。在數學建模中清晰記錄你的處理流程、每一步的理由以及參數選擇的依據比追求一個“完美”的結果更重要。因為評審專家看的正是你面對真實、混亂數據時所展現出的科學決策能力和嚴謹態度。希望這篇長文能為你提供一套可復現、可解釋的方法論工具箱。