- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
我有一张图像,想为每个非零值像素计算它到最近的零值像素的距离。我尝试这样做的方式如下:
import numpy as np
from scipy.spatial.distance import cdist
from skimage import io
im=io.imread('imagepath')
#getting array where elements are 0
a,b = np.where(im == 0)
# create a list with (row,column)
x = list(zip(a, b))
#getting array where elements are non zero
r, c =np.where(im!=0)
#create a list with (row, column) of all non 0 values
#note coordinates are in y, x format
obj = list(zip(r,c))
dist_dict={}
#calculating for each pixel of the object
for o in obj:
d = (cdist(np.array([o]), x, metric='euclidean')).min()
dist_dict.update({o:d})
我相信这应该可行,但速度很慢。对于单个像素,计算距离大约需要 0.2 秒。对于大约 50.000 像素大的物体,因此每张图像需要大约三个小时的计算时间,这是根本不可行的。我在这里看到的一个问题是我只是遍历所有非零像素。有没有一种方法不是从数组的开头而是从当前坐标开始搜索,直到找到零值?或者有任何其他建议可以加快这个过程吗?
最佳答案
您可以使用 scipy.ndimage.morphology.distance_transform_edt它找到距离输入像素欧氏距离最小的最近背景点(值为 0)。
from scipy import ndimage
import pprint
def nearest_zero(image):
" Finds closest background (zero) element for each element in image "
# Find closest zero elements in the inverted image (same as closest non-zero for image)
edt = ndimage.distance_transform_edt(image, return_indices=False)
# Create dictionary of indexes
return {(r,c):edt[r][c] for r in range(image.shape[0]) for c in range(image.shape[1]) if image[r][c]}
使用示例
image = np.array(([0,50,75,15,0],
[2,0,111,10,15],
[0,112,25,110,115],
[0,10,110,115,0],
[15,12,115,0,0]))
d = nearest_zero(image)
pp = pprint.PrettyPrinter(indent=4)
print('Original Image')
print(image)
print('\nDictionary of Distances to closest background pixel for each non-background pixel')
pp.pprint(sorted(d.items(), key=lambda x: x[0]))
输出
Original Image
[[ 0 50 75 15 0]
[ 2 0 111 10 15]
[ 0 112 25 110 115]
[ 0 10 110 115 0]
[ 15 12 115 0 0]]
Dictionary of Distances to closest background pixel for each non-background pixel
[ ((0, 1), 1.0),
((0, 2), 1.4142135623730951),
((0, 3), 1.0),
((1, 0), 1.0),
((1, 2), 1.0),
((1, 3), 1.4142135623730951),
((1, 4), 1.0),
((2, 1), 1.0),
((2, 2), 1.4142135623730951),
((2, 3), 1.4142135623730951),
((2, 4), 1.0),
((3, 1), 1.0),
((3, 2), 1.4142135623730951),
((3, 3), 1.0),
((4, 0), 1.0),
((4, 1), 1.4142135623730951),
((4, 2), 1.0)]
性能测试
结果:SciPy 快约 100 倍
测试数据生成 -- 随机图像(大小 250x250 = 62、500 像素)
import random
size = 250
z = [random.randrange(0, 255) for r in range(size) for c in range(size)]
image = np.array(z).reshape(size, size)
测试图像中的数字零
print(np.count_nonzero(image==0)) # 62262
timeit 使用原始帖子中的方法
11.6 s ± 89.1 ms per loop (mean ± std. dev. of 7 runs, 1 loop each)
which is an average of 1.8e-04 seconds per non-zero point
timeit 使用 SciPy 方法
119 ms ± 17.1 ms per loop (mean ± std. dev. of 7 runs, 10 loops each)
Which is an average of 1.9e-06
因此 SciPy 快了约 100 倍
关于python - 在python中计算矩阵中值距离的有效方法,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/59126316/
我遇到了一个奇怪的问题。我有这个: $(document).ready(function () {
我正在编写一个程序,它从列表中读取一些 ID,从中找出不同的 URL,然后将图像保存到我的 C: 驱动器中。 如果我在浏览器中导航到图像 URL,它们就会起作用。此外,如果我尝试从不同的服务器获取图像
我编写了一个 REST WCF RIA Silverlight 4.0 兼容服务,我可以从 javascript + jQuery.1.4.2.js + JSON2.js(当然,还可以从 .NET 4
我很确定这个网站实际上还没有得到回答。一劳永逸地,与 32 位有符号整数范围内的数字字符串匹配的最小正则表达式是什么,范围是 -2147483648至 2147483647 . 我必须使用正则表达式进
我有两个data.table;我想从那些与键匹配的元素中随机分配一个元素。我现在这样做的方式相当慢。 让我们具体点;这是一些示例数据: dt1<-data.table(id=sample(letter
我已经安装了 celery 、RabitMQ 和花。我可以浏览到花港。我有以下简单的工作人员,我可以将其附加到 celery 并从 python 程序调用: # -*- coding: utf-8 -
我正在使用 ScalaCheck 在 ScalaTest 中进行一些基于属性的测试。假设我想测试一个函数,f(x: Double): Double仅针对 x >= 0.0 定义的, 并返回 NaN对于
我想检查文件是否具有有效的 IMAGE_DOS_SIGNATURE (MZ) function isMZ(FileName : String) : boolean; var Signature: W
在 Herbert Schildt 的“Java:完整引用,第 9 版”中,有一个让我有点困惑的例子。它的关键点我无法理解可以概括为以下代码: class Test { public stat
我在工作中查看了一些代码,发现了一些我以前没有遇到过的东西: for (; ;) { // Some code here break; } 我们一直调用包含这个的函数,我最近才进去看看它是
在 Herbert Schildt 的“Java:完整引用,第 9 版”中,有一个让我有点困惑的例子。它的关键点我无法理解可以概括为以下代码: class Test { public stat
我试图编写一个函数,获取 2D 点矩阵和概率 p 并以概率 p 更改或交换每个点坐标 所以我问了一个question我试图使用二进制序列作为特定矩阵 swap_matrix=[[0,1],[1,0]]
这个问题在这里已经有了答案: Using / or \\ for folder paths in C# (5 个答案) 关闭 7 年前。 我在某个Class1中有这个功能: public v
PostgreSQL 10.4 我有一张 table : Column | Type ------------------------- id | integer| title
我正在 Postgresql 中编写一个函数,它将返回一些针对特定时区(输入)计算的指标。 示例结果: 主要问题是这只是一个指标。我需要从其他表中获取其他 9 个指标。 对于实现此目标的更简洁的方法有
我需要在 python 中模拟超几何分布(用于不替换采样元素的花哨词)。 设置:有一个装满人口许多弹珠的袋子。弹珠有两种类型,红色和绿色(在以下实现中,弹珠表示为 True 和 False)。从袋子中
我正在使用 MaterializeCSS 框架并动态填充文本输入。我遇到的一个问题是,在我关注该字段之前,valid 和 invalid css 类不会添加到我的字段中。 即使我调用 M.update
是否有重叠 2 个 div 的有效方法。 我有以下内容,但无法让它们重叠。 #top-border{width:100%; height:60px; background:url(image.jpg)
我希望你们中的一位能向我解释为什么编译器要求我在编译单元中重新定义一个静态固定长度数组,尽管我已经在头文件中这样做了。这是一个例子: 我的类.h: #ifndef MYCLASS_H #define
我正在使用旧线程发布试图解决相同问题的新代码。什么是安全 pickle ? this? socks .py from socket import socket from socket import A
我是一名优秀的程序员,十分优秀!