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

相关文章推荐

神勇威武的红酒 · Python dat文件读入 ...· 2 周前 ·

豁达的红烧肉 · Ant Design Vue Tree ...· 2 周前 ·

挂过科的人字拖 · 数据分析与可视化项目技术参考-腾讯云开发者社 ...· 1 周前 ·

温柔的豆浆 · TypeScript加载json数组_mob ...· 1 周前 ·

眉毛粗的草稿纸 · Python -更新现有json文件中的项 ...· 1 周前 ·

老实的牛肉面 · 在PHP中，为什么imagecreatefr ...· 5 月前 ·

胡子拉碴的绿豆 · maven打包忽略单元测试-掘金· 10 月前 ·

淡定的蚂蚁 · select case 语句中case ...· 1 年前 ·

暴躁的香烟 · matlab msgbox ...· 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 dat文件读入 python读取dat文件内容_mob64ca1406d617的技术博客_51CTO博客

2 周前

豁达的红烧肉 · Ant Design Vue Tree 树形获取父节点信息_a-tree-select 怎么拿到父value-CSDN博客

2 周前

挂过科的人字拖 · 数据分析与可视化项目技术参考-腾讯云开发者社区-腾讯云

1 周前

温柔的豆浆 · TypeScript加载json数组_mob649e815e9bc9的技术博客_51CTO博客

1 周前

眉毛粗的草稿纸 · Python -更新现有json文件中的项 - 腾讯云开发者社区 - 腾讯云

1 周前

老实的牛肉面 · 在PHP中，为什么imagecreatefromjpeg“不是jpeg文件”更改为致命错误？-腾讯云开发者社区-腾讯云

5 月前

胡子拉碴的绿豆 · maven打包忽略单元测试-掘金

10 月前

淡定的蚂蚁 · select case 语句中case 后跟的判断语句为什么不能用like进行字符串的配比，有别的替代方法吗？十分感谢。_百度知道

1 年前

暴躁的香烟 · matlab msgbox 警告,Excel VBA关闭消息提示警告窗口_晚安呼的博客-CSDN博客

1 年前

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

2 个回答