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

相关文章推荐

有腹肌的卡布奇诺 · js ...· 2 月前 ·

独立的书签 · docker ...· 昨天 ·

满身肌肉的拖把 · “30万投《水门桥》可赚百万”“易烊千玺主演 ...· 6 月前 ·

想发财的茴香 · 如何选购手机？手机处理器强弱怎么判断？看完这 ...· 1 年前 ·

听话的跑步鞋 · 合纵六国的苏秦，为何要求将自己车裂？|齐国| ...· 1 年前 ·

沉稳的脆皮肠 · 有哪些「小众高级」的圣诞礼物？ - 知乎· 2 年前 ·

仗义的夕阳 · 联想y700键盘背光新款 - ...· 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)

推荐文章

有腹肌的卡布奇诺 · js jquery-ajax/fetch请求时数据文本丢失加号和连接号的问题 - 深入学习ing

2 月前

独立的书签 · docker compose入门6—如何挂载卷_docker-compose volumes 挂载目录

昨天

满身肌肉的拖把 · “30万投《水门桥》可赚百万”“易烊千玺主演差不了”…小心电影投资也有杀猪盘|易烊千玺_新浪财经_新浪网

6 月前

想发财的茴香 · 如何选购手机？手机处理器强弱怎么判断？看完这篇文章你就知道了_百科TA说

1 年前

听话的跑步鞋 · 合纵六国的苏秦，为何要求将自己车裂？|齐国|秦国|齐宣王|燕国_网易订阅

1 年前

沉稳的脆皮肠 · 有哪些「小众高级」的圣诞礼物？ - 知乎

2 年前

仗义的夕阳 · 联想y700键盘背光新款 - 联想y700键盘背光2021年新款 - 京东

2 年前

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

2 个回答