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

相关文章推荐

另类的蚂蚁 · 使用pandas读取csv文件指定的前几行_ ...· 3 月前 ·

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

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

潇洒的煎鸡蛋 · 天津市博物馆、图书馆、出版社、文艺院团等推出 ...· 11 月前 ·

体贴的黄瓜 · 经济学原理(第8版) 微观经济学分册曼昆 ...· 1 年前 ·

文武双全的板栗 · 民盟静安区委赴昆明开展主题教育活动· 1 年前 ·

胆小的毛衣 · 国家税务总局浙江省税务局总局文件财政部 ...· 1 年前 ·

满身肌肉的水桶 · 阿拉善盟行政公署旗区新闻 ...· 2 年前 ·

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)

推荐文章

另类的蚂蚁 · 使用pandas读取csv文件指定的前几行_pandas读取csv前几行

3 月前

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

4 周前

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

6 天前

潇洒的煎鸡蛋 · 天津市博物馆、图书馆、出版社、文艺院团等推出丰富的“六一”系列活动

11 月前

体贴的黄瓜 · 经济学原理(第8版) 微观经济学分册曼昆课后习题答案 - 哔哩哔哩

1 年前

文武双全的板栗 · 民盟静安区委赴昆明开展主题教育活动

1 年前

胆小的毛衣 · 国家税务总局浙江省税务局总局文件财政部税务总局关于延续实施宣传文化增值税优惠政策的公告

1 年前

满身肌肉的水桶 · 阿拉善盟行政公署旗区新闻 “美在阿拉善——岩画与居延汉简艺术展”在宝鸡青铜博物院开展

2 年前

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

2 个回答