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

相关文章推荐

挂过科的机器猫 · [Python] Pandas ...· 3 周前 ·

温暖的弓箭 · MATLAB直方图全解析_matlab ...· 1 周前 ·

火爆的电影票 · 关于公开征求《蚌埠市城市安全和防灾减灾体系建 ...· 9 月前 ·

淡定的香烟 · 江苏省科技厅科技动态 ...· 1 年前 ·

道上混的电影票 · 亿万富翁袁宝璟遗孀和儿子被判还债220万元- ...· 1 年前 ·

乖乖的芹菜 · 千年前南京人自称“侬”_手机新浪网· 1 年前 ·

善良的红茶 · 【五育并举促全面发展】2022级研究生体育教 ...· 1 年前 ·

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)

推荐文章

挂过科的机器猫 · [Python] Pandas 对数据进行查找、替换、筛选、排序、重复值和缺失值处理 - 哆啦梦乐园

3 周前

温暖的弓箭 · MATLAB直方图全解析_matlab histogram

1 周前

火爆的电影票 · 关于公开征求《蚌埠市城市安全和防灾减灾体系建设国土空间专项规划（2021—2035年）》草案意见建议的公告_蚌埠市应急管理局

9 月前

淡定的香烟 · 江苏省科技厅科技动态南京麒麟科创园“两区三城”建设展新貌——在麒麟，专“新”致“质”打造创新高地

1 年前

道上混的电影票 · 亿万富翁袁宝璟遗孀和儿子被判还债220万元-搜狐新闻

1 年前

乖乖的芹菜 · 千年前南京人自称“侬”_手机新浪网

1 年前

善良的红茶 · 【五育并举促全面发展】2022级研究生体育教育侧记-西安交通大学新闻网

1 年前

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

2 个回答