如何使用 scipy.ndimage.filters.gereric_filter?

2023-11-21

我正在尝试使用 scipy.ndimage.filters.generic_filter 来计算邻域的加权和。邻域在某些时候会发生变化,但目前我正在努力实现 3x3。 到目前为止,这就是我所在的地方:

    def Func(a):
         a = np.reshape((3,3))
         weights = np.array([[0.5,.05,0.5],[0.5,1,0.5],[0.5,0.5,0.5]])
         a = np.multiply(a,weights)
         a = np.sum(a)
         return a

ndimage.filters.generic_filter(Array,Func,footprint=np.ones((3,3)),mode='constant',cval=0.0,origin=0.0)

我从 ndimage 收到一个错误,提示“TypeError: a float is required”,但我不知道它指的是什么参数,它看起来与我见过的其他示例基本相同。


这对我有用。代码有几个小问题:

import scipy.ndimage.filters
import numpy as np

Array = rand( 100,100 )

def Func(a):
    a = a.reshape((3,3))
    weights = np.array([[0.5,.05,0.5],[0.5,1,0.5],[0.5,0.5,0.5]])
    a = np.multiply(a,weights)
    a = np.sum(a)
    return a

out = scipy.ndimage.filters.generic_filter(Array,Func,footprint=np.ones((3,3)),mode='constant',cval=0.0,origin=0.0)

You had a = np.reshape( (3,3) )这是不正确的。那是你要的吗?

[update]

根据我们的讨论对此进行一些清理:

import scipy.ndimage.filters
import numpy as np

Array = rand( 100,100 )

def Func(a):
    return np.sum( a * r_[0.5,.05,0.5, 0.5,1,0.5, 0.5,0.5,0.5] )

out = scipy.ndimage.filters.generic_filter(Array,Func,footprint=np.ones((3,3)),mode='constant',cval=0.0,origin=0.0)
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

如何使用 scipy.ndimage.filters.gereric_filter? 的相关文章

  • 无法在 PyCharm 版本 9.3.3 中安装 NumPy。 Python版本3.8.2

    在 PyCharm 中安装 NumPy 时出错 尝试安装 Microsoft Visual C 14 0 还是行不通 NumPy 正在通过命令安装pip3 install numpy在 cmd 终端中 但是当尝试将其安装在 PyCharm
  • numpy polyfit 中使用的权重值是多少以及拟合误差是多少

    我正在尝试对 numpy 中的某些数据进行线性拟合 Ex 其中 w 是该值的样本数 即对于点 x 0 y 0 我只有 1 个测量值 该测量值是2 2 但对于这一点 1 1 我有 2 个测量值 值为3 5 x np array 0 1 2 3
  • numpy python 中的“AttributeError:'matrix'对象没有属性'strftime'”错误

    我有一个维度为 72000 1 的矩阵 该矩阵涉及时间戳 我想使用 strftime 如下所示 strftime d m y 为了得到像这样的输出 11 03 02 我有这样一个矩阵 M np matrix timestamps 我使用了
  • 如何在Python中高效地添加稀疏矩阵

    我想知道如何在Python中有效地添加稀疏矩阵 我有一个程序 可以将大任务分解为子任务 并将它们分配到多个 CPU 上 每个子任务都会产生一个结果 一个 scipy 稀疏矩阵 格式为 lil matrix 稀疏矩阵尺寸为 100000x50
  • 如何通过减少请求来改进 AJAX 实时搜索

    我正在构建一个 AJAX 实时搜索页面 到目前为止 一切都按预期运行 但我注意到我正在进行大量的 AJAX 调用 我知道发生这种情况的地点和原因 但我找不到阻止这些 AJAX 调用发生的方法 我将尝试给出快速解释 然后粘贴下面的代码 在页面
  • Matplotlib 与多处理冻结计算机

    我对 matplotlib 和多重处理有疑问 我启动第一个进程 在其中显示图像并选择一个区域 然后关闭图形 然后我启动另一个进程 在其中调用定期更新的图形函数 至此 一切正常 然后 当我尝试使用相同的图形功能启动另一个进程时 它冻结了我的整
  • 张量流服务错误:参数无效:JSON 对象:没有命名输入

    我正在尝试使用 Amazon Sagemaker 训练模型 并且希望使用 Tensorflow 服务来为其提供服务 为了实现这一目标 我将模型下载到 Tensorflow 服务 docker 并尝试从那里提供服务 Sagemaker 的训练
  • Python 2:SMTPServerDisconnected:连接意外关闭

    我在用 Python 发送电子邮件时遇到一个小问题 me my email address you recipient s email address me email protected cdn cgi l email protectio
  • 如何通过 TLS 1.2 运行 django runserver

    我正在本地 Mac OS X 机器上测试 Stripe 订单 我正在实现这段代码 stripe api key settings STRIPE SECRET order stripe Order create currency usd em
  • 如何在 Windows 命令行中使用参数运行 Python 脚本

    这是我的蟒蛇hello py script def hello a b print hello and that s your sum sum a b print sum import sys if name main hello sys
  • 格式化多维数组 Python

    如何编写代码来更改多维数组中每个单独数组的值a出现负值后立即归零 所以里面的第二个数组a有负值 12 34 5 6 88 10 30 75 10 会将其所有值及其后的值转为零 将数组变成 12 34 5 6 88 0 0 0 我怎样才能获得
  • 在 python 中使用 pandas 计算行的出现次数

    我有一个包含数千行和 4 列的 pandas 数据框 IE A B C D 1 1 2 0 3 3 2 1 3 1 1 0 有没有办法统计某一行出现了多少次 例如 可以找到多少次 3 1 1 0 并返回这些行的索引 如果你只寻找一行 那么我
  • numpy NPV 和 Excel NPV 有区别吗?

    我的 Excel 中有一行包含 11 个值 TotalSavings 0 8000 8000 8000 8000 8000 8000 8000 8000 8000 8000 贴现率为 0 08 我在 Excel 中使用 计算 NPVNPV
  • Sklearn 随机 Logistic 回归给出错误“ValueError:类的数量必须大于 1”

    我发现了 sklearn Randomized Logistic 中的一个错误 由于我花了很长时间才解决它 所以我会将其发布在这里 以防其他人遇到同样的问题 发生的情况是 在格式完美的数据上 sklearn RandomizedLogist
  • scipy.optimize on pandas dataframe

    我试图搜索它 但结果很差 有人可以向我解释一下如何在 Pandas DataFrame 上执行 optimize minimize 以便最小化 DataFrame 中的类别和结果列之间的错误 考虑这个例子 import pandas as
  • 通过 python 中的另外两个修改数组[重复]

    这个问题在这里已经有答案了 假设我们有三个一维数组 A 长度为 5 B 长度相同 示例中为5 C 更长 比如长度为 100 C最初用零填充 A给出索引C应更改的元素 它们可能会重复 以及B给出应添加到初始零的值C 例如 如果A 1 3 3
  • ValueError:数据必须为正(boxcox scipy)

    我正在尝试将我的数据集转换为正态分布 0 8 298511e 03 1 3 055319e 01 2 6 938647e 02 3 2 904091e 02 4 7 422441e 02 5 6 074046e 02 6 9 265747e
  • 使用 NumPy 将非均匀数据从文件读取到数组中

    假设我有一个如下所示的文本文件 33 346 1223 10 23 11 23 12 23 13 23 14 23 15 23 16 24 10 24 11 24 12 24 13 24 14 24 15 24 16 25 14 25 15
  • 如何在 OSX 上安装 numpy 和 scipy?

    我是 Mac 新手 请耐心等待 我现在使用的是雪豹 10 6 4 我想安装numpy和scipy 所以我从他们的官方网站下载了python2 6 numpy和scipy dmg文件 但是 我在导入 numpy 时遇到问题 Library F
  • Ubuntu 上的 Python 2.7

    我是 Python 新手 正在 Linux 机器 Ubuntu 10 10 上工作 它正在运行 python 2 6 但我想运行 2 7 因为它有我想使用的功能 有人敦促我不要安装 2 7 并将其设置为我的默认 python 我的问题是 如

随机推荐

  • 一旦 std:thread 进入 C++Ox,pthreads 就会过时吗?

    Closed 这个问题是基于意见的 目前不接受答案 显然我们仍然会维护它 但是一旦 C 标准保证了 它会有多大用处 随着新标准的出现 同步原语 互斥体 条件变量 怎么样 您是否认为 pthread 比 std thread 更难掌握 C不会
  • SIMD 代码在 Debug 中有效,但在 Release 中无效

    该代码在调试模式下工作 但由于发布模式下的断言而出现恐慌 use std arch x86 64 fn main unsafe let a vec 2 0f32 0 0 0 0 0 0 0 0 0 0 0 0 0 0 let b 1 0f3
  • Google Speech API 返回 403 PERMISSION_DENIED

    几个月来 我一直在使用 Google Speech API 将我的 PHP 应用程序 使用 Google Cloud PHP 客户端 中的音频转录为文本 没有出现任何问题 但我的呼叫现在开始返回 403 错误 状态为 PERMISSION
  • 无法将 XML 包安装到 R

    我需要将 XML 包安装到 R 它总是报告 XML 不适用于 R 版本 3 1 1 我尝试下载 tar 但仍然无法正常工作 谁能帮我解决这个问题吗 gt install packages XML repos http www omegaha
  • 更改 Lumen 或 Laravel 5 中的时区

    我正在使用流明框架 如何将时区更改为欧洲 巴黎 CEST 我在我的中添加了一个变量 env file APP TIMEZONE Europe Paris 但这行不通 更新时区的正确方法是什么 您可以设置您的应用时区通过配置app php文件
  • io.BufferedReader peek 函数返回缓冲区中的所有文本

    我在 Windows 8 上使用 Python 3 4 1 我想读取一个带有缓冲接口的文件 该接口允许我提前查看一定数量的字节以及读取字节 io BufferedReader似乎是正确的选择 很遗憾 io BufferReader peek
  • 三种 .htaccess 过期规则之间的区别

    以下三个 htaccess 规则之间有什么区别以及何时使用每个规则 是否有任何特殊用例更喜欢其中一个 Header set Cache Control max age 290304000 Header set Expires Thu 15
  • 安装plpython时Postgres数据库崩溃

    我试图在 Postgres 9 1 中安装 plpython 但它使服务器崩溃 postgres dataserver1 gt opt postgres 9 1 bin psql d mydb psql bin 9 1 4 Type hel
  • Android Studio Gradle 未能完成 Gradle 执行

    我正在尝试运行 appengineUpdate 来部署我的后端 但收到以下错误消息 Failed to complete Gradle execution Cause Gradle build daemon disappeared unex
  • 如何在sitemesh装饰器中获取模型属性或spring的bean?

    我正在使用 Spring 3 和 sitemesh 我想参考 sitemesh 中定义的装饰器页面中的 spring 上下文 bean 问题是 SiteMesh 过滤器在 Spring 上下文之外工作 因此 sitemesh 装饰器 jsp
  • 从多元正态分布中有效地随机抽取

    只是想知道是否有人遇到过这样的问题 他 她需要从非常高维的多元正态分布 假设维度 10 000 中随机抽取 因为rmvnorm的功能mvtnorm包对此来说是不切实际的 我知道这个article has an Rcpp实施dmvnorm的功
  • R 可以读取 html 编码的表情符号字符吗?

    Question 我的问题如下 如何使用 R 读取包含 HTML 表情符号代码的字符串 例如 55358 56599 我想 1 表示表情符号 例如 作为 unicode 符号 在解析后的字符串中 OR 2 将其转换为等效文本 hugging
  • 如何在 Java 中将分数格式的字符串转换为 Decimal 或 Float?

    我从数据库中获取的字符串值很少 例如 1 4 2 3 但是在显示为 Android ListView 内容时 我需要将其显示为0 25 0 66 现在我不想分割字符串 然后将单个字符串转换为数字 然后 然后将它们分开即可得到结果 有谁知道
  • 如何确保我的 bash 脚本尚未运行?

    我有一个 bash 脚本 我想从 cron 每 5 分钟运行一次 但是该脚本的上一次运行可能尚未完成 在这种情况下 我希望新的运行退出 我不想只依赖 tmp 中的锁定文件 我想确保在执行锁定文件 或其他文件 之前该进程实际上正在运行 这是到
  • 为什么 [-1] 不返回文件中该行的最后一个字符?

    我想打印 python 中从文件中读取的字符串的最后一个字符 我打电话作为str 1 但它没有按预期工作 t txt 包含 Do not laugh please 9 Are you kidding me 4 我的代码是 with open
  • 如何使用 FLAG_GRANT_READ_URI_PERMISSION 授予对自定义内容提供程序的临时访问权限

    我正在尝试从另一个应用程序 应用程序 B 查询自定义内容提供程序 应用程序 A 当内容提供商没有权限保护时我可以这样做 具体来说 我在应用程序 A 上构建自定义内容提供程序 并将包含 URI 的意图发送到应用程序 B 这是App A中的意图
  • GCM 无效发件人错误

    我被一个神秘的问题困扰了一段时间 所以我在我的 Android 应用程序中实现了 GCM 但收到 INVALID SENDER 错误 我的SENDER ID是来自google项目的url的id Android 推送通知已开启 奇怪的是我已经
  • 如何卸载Python和所有包

    我希望卸载 Python 2 7 以及与其连接的所有包 我最初从官方网站安装了Python 并使用 pip install 命令安装了所有软件包 从控制面板卸载Python也会自动卸载所有包吗 我之所以要卸载Python 是因为我想使用An
  • 避免在 Kubernetes 中为一个 cron 执行点运行多个 cron 作业

    编辑 问题已解决 这是我的错误 我只是使用了错误的 cron 设置 我假设 2 每天只会在 2 点运行一次 但实际上它在 2 点之后每分钟都会运行一次 因此 Kubernetes 的行为是正确的 我一直在一个 cron 执行点运行多个作业
  • 如何使用 scipy.ndimage.filters.gereric_filter?

    我正在尝试使用 scipy ndimage filters generic filter 来计算邻域的加权和 邻域在某些时候会发生变化 但目前我正在努力实现 3x3 到目前为止 这就是我所在的地方 def Func a a np resha