获取二维数组中高于某一值的局部最大值的坐标

相关文章推荐

光明磊落的花生 · 退出页面前发送请求实现方案 - 公瑾当年 · 2 月前 ·

飘逸的口罩 · 如何将CSV导入的嵌套列表转换为指定列数据类 ...· 1 月前 ·

气宇轩昂的豆腐 · Apache ...· 1 年前 ·

暴走的火锅 · tuple Class | ...· 1 年前 ·

稳重的企鹅 · C++ STL ...· 2 年前 ·

文武双全的鸵鸟 · C语言字符串详解-腾讯云开发者社区-腾讯云· 2 年前 ·

飘逸的蘑菇 · 【图像处理】python实现图像标准化_py ...· 3 年前 ·

from PIL import Image
import numpy as np
from scipy.ndimage.filters import maximum_filter
import pylab
# the picture (256 * 256 pixels) contains bright spots of which I wanna get positions
# problem: data has high background around value 900 - 1000
im = Image.open('slice0000.png')
data = np.array(im)
# as far as I understand, data == maximum_filter gives True-value for pixels
# being the brightest in their neighborhood (here 10 * 10 pixels)
maxima = (data == maximum_filter(data,10))
# How can I get only maxima, outstanding the background a certain value, let's say 500 ?

import numpy as np
import scipy
import scipy.ndimage as ndimage
import scipy.ndimage.filters as filters
import matplotlib.pyplot as plt
fname = '/tmp/slice0000.png'
neighborhood_size = 5
threshold = 1500
data = scipy.misc.imread(fname)
data_max = filters.maximum_filter(data, neighborhood_size)
maxima = (data == data_max)
data_min = filters.minimum_filter(data, neighborhood_size)
diff = ((data_max - data_min) > threshold)
maxima[diff == 0] = 0
labeled, num_objects = ndimage.label(maxima)
xy = np.array(ndimage.center_of_mass(data, labeled, range(1, num_objects+1)))
plt.imshow(data)

推荐文章

光明磊落的花生 · 退出页面前发送请求实现方案 - 公瑾当年

2 月前

飘逸的口罩 · 如何将CSV导入的嵌套列表转换为指定列数据类型的NumPy数组？

1 月前

气宇轩昂的豆腐 · Apache Tomcat修改默认400/404/500页面 - 银河系派件员 - 博客园

1 年前

暴走的火锅 · tuple Class | Microsoft Learn

1 年前

稳重的企鹅 · C++ STL vector删除元素的几种方式（超级详细）

2 年前

文武双全的鸵鸟 · C语言字符串详解-腾讯云开发者社区-腾讯云

2 年前

飘逸的蘑菇 · 【图像处理】python实现图像标准化_python图像标准化_BIT可达鸭的博客-CSDN博客

3 年前

获取二维数组中高于某一值的局部最大值的坐标 内容来源于 Stack Overflow，遵循 CC BY-SA 4.0 许可协议进行翻译与使用。IT领域专用引擎提供翻译支持 腾讯云小微IT领域专用引擎提供翻译支持

2 个回答