讲师中心 微信公众号
AI工具推荐 视频效率加速

在Python中怎样通过NumPy掩码高效剔除离群异常值

梦敏吖_6387

梦敏吖_6387

发布时间:2026-10-29 06:56:32

|

671人浏览过

|

来源于php中文网

原创

最常用但需注意:zscore默认按列计算,一维数组须指定axis=0;其均值和标准差易受离群值污染,建议改用中位数与MAD构造稳健阈值;IQR法更鲁棒,应使用np.quantile(arr, [0.25, 0.75])计算Q1、Q3,并用nanpercentile处理含NaN数据。

在python中怎样通过numpy掩码高效剔除离群异常值

用 np.abs(zscore) 做标准差掩码最常用但要注意中心化

直接用 scipy.stats.zscore 计算 z 值再设阈值(如 3)是最直观的离群值剔除方式,但容易忽略两点:一是 zscore 默认按列计算,一维数组要加 axis=0 防止维度错乱;二是它内部用均值和标准差中心化,而均值本身会被离群值污染——所以实际中建议先用中位数和 MAD(中位数绝对偏差)替代,或先粗筛再重算统计量。

实操建议:

  • 对一维数据,务必传 axis=0,否则 zscore(arr) 可能返回标量或形状异常
  • 若数据含明显偏态或重尾,改用 scipy.stats.median_abs_deviation + 中位数构造稳健阈值
  • zscore 掩码后用布尔索引取子集,别用 np.delete,后者需传入索引数组,效率低且易出错

布尔索引配合 np.quantile 做 IQR 掩码更鲁棒

IQR(四分位距)法不依赖正态假设,对非高斯分布更安全。核心是计算 Q1 = np.quantile(arr, 0.25) 和 Q3 = np.quantile(arr, 0.75),然后定义上下界:low_bound = Q1 - 1.5 * (Q3 - Q1),up_bound = Q3 + 1.5 * (Q3 - Q1)。注意 np.quantile 默认插值方式是 linear,小样本时可考虑 method='midpoint' 减少边界抖动。

常见错误现象:

立即学习“Python免费学习笔记(深入)”;

macOS UI 自动化技能包 —— 用 Python 控制你的 Mac
macOS UI 自动化技能包 —— 用 Python 控制你的 Mac

让 macOS 自动化像说话一样简单 —— 用代码控制你的 Mac,无需学习 AppleScript 系统要求 macOS 14+ (Sonoma) Python 3.9+ 需要开启:辅助功能 + 屏幕录制权限

下载
  • 直接写 arr[(arr > low_bound) & (arr 忘了括号,导致 <code>& 运算符优先级高于比较,报 ValueError: The truth value of an array with more than one element is ambiguous
  • 对空数组或全 NaN 数组调用 np.quantile 会返回 nan,后续边界计算失效;应提前用 np.isfinite(arr).sum() > 0 校验
  • 多维数组未指定 axis,默认在展平后计算,结果不符合逐列/逐行剔除预期

用 np.where + np.nanpercentile 处理含 NaN 的真实数据

传感器或日志数据常含 nan,而 np.quantile 遇到 nan 直接返回 nan。此时必须换用 np.nanpercentile,它自动跳过 nan 计算分位数。但注意:它不支持 method 参数(仅 numpy>=1.22 版本起支持 interpolation,已弃用),老版本只能接受默认线性插值。

使用场景:

  • 原始数据含缺失值且不能简单删除整行(如时间序列需保持续性)
  • 需保留原始数组 shape,只把离群位置置为 nan,后续用 pd.Series.interpolate() 或 scipy.ndimage.uniform_filter1d 平滑
  • 想定位离群索引而非仅取子集,用 np.where((arr up_bound)) 返回元组,解包为 outlier_idx = np.where(...)[0]

性能关键:避免重复计算统计量与掩码链式调用

在循环中反复对同一数组调用 np.quantile 或 zscore 是典型性能陷阱。NumPy 统计量函数本身不缓存,每次都是全量扫描。若需多次用相同阈值过滤不同字段(如 DataFrame 多列),应预先算好 Q1/Q3 或 mean/std,再广播到布尔表达式。

实操建议:

  • 对二维数组做逐列 IQR,用 np.quantile(arr, [0.25, 0.75], axis=0) 一次获取两组分位数,比两次调用快 40%+
  • 掩码生成后立即用于索引,不要拆成 mask = ... 和 arr[mask] 两步——除非需复用该 mask,否则中间变量无益于性能
  • 超大数据(>1e7 元素)慎用 zscore,因它内部会显式构造均值广播数组,内存翻倍;改用原地减均值再除标准差更省空间

离群值剔除没有银弹:zscore 快但怕偏态,IQR 稳但对小样本分位数估计不准,真正复杂的数据往往得组合用——比如先 IQR 粗筛,再对剩余数据用 zscore 细筛。别忘了检查掩码后是否还剩足够样本量,尤其在分组操作中,空组会静默破坏下游逻辑。

热门AI工具

更多
WorkBuddy

一款AI办公效率工具,主要用于腾讯云推出的AI原生桌面智能体工作台,适合需要提升相关任务效率的用户。

DeepSeek

DeepSeek是一款面向对话、写作、编程和推理场景的AI大模型工具。

超级简历WonderCV

一款AI办公效率工具,主要用于免费求职简历模版下载制作,应届生职场人必备简历制作神器,适合需要提升相关任务效率的用户。

切问学术

切问学术是一款AI论文写作工具,复旦大学NLP团队推出的AI学术智能体。

Seko
Seko Hot

一款AI视频创作工具,主要用于商汤科技推出的创编一体的AI短视频创作Agent,适合需要提升相关任务效率的用户。

火山引擎

火山引擎是一款面向企业的云计算与AI服务平台。

豆包大模型

豆包大模型是一款由字节跳动推出的企业级大语言模型服务平台。

立刻MV
立刻MV Hot

立刻MV是一款AI文本写作工具,AI 音乐视频(MV)创作工具。

超级麦吉

一款通过数字员工协作完成调研、数据分析与文档制作的AI办公平台。

相关专题

更多
python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

1831

2023.07.20

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

4884

2023.07.25

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

1789

2023.07.31

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

29137

2023.08.03

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

3487

2023.08.04

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

3487

2023.08.04

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

1223

2023.08.11

python合并两个列表
python合并两个列表

Python是一种强大的编程语言,具有许多方便的功能和工具。在Python中,有多种方法可以合并两个列表。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

616

2023.08.10

Vite 项目如何打包和部署
Vite 项目如何打包和部署

Vite项目打包部署专题汇集生产构建、产物预览、输出目录、资源路径、Nginx与GitHub Pages部署及常见白屏和路由404问题的已审核文章,帮助读者完成上线流程。

40

2026.10.28

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
Uvicorn官方安装与入门手册
Uvicorn官方安装与入门手册

共0课时 | 0人学习

Python uv官方项目管理指南
Python uv官方项目管理指南

共0课时 | 0人学习

Python uv官方安装与入门文档
Python uv官方安装与入门文档

共0课时 | 0人学习

关于我们 免责申明 举报中心 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号
PHP中文网订阅号
每天精选资源文章推送

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn