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

相关文章推荐

性感的猴子 · Chapter 8 Basic: read ...· 1 周前 ·

纯真的白开水 · 顶顶通呼叫中心中间件(mod_cti基于Fr ...· 6 天前 ·

没读研的大葱 · Query and transform ...· 2 天前 ·

重感情的斑马 · 3DMax打不开代码错误无法初始化怎么办？ ...· 1 月前 ·

一身肌肉的烤红薯 · 北京市体育局2018年政府信息公开工作年度报 ...· 2 月前 ·

打盹的板栗 · 社会艺术水平考级管理办法(文化部令第31号)· 2 月前 ·

魁梧的饭盒 · 上海外滩踩踏当晚官员“吃喝门”真相追踪-新华网· 8 月前 ·

直爽的黄花菜 · 一个方法自测血管是否硬化！这8个指标控制好， ...· 9 月前 ·

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)

推荐文章

性感的猴子 · Chapter 8 Basic: read and write | R 语言：从基础到进阶

1 周前

纯真的白开水 · 顶顶通呼叫中心中间件(mod_cti基于FreeSWITCH)-拨号方案和路由配置 | 顶顶通 - 呼叫中心二次开发接口,smartivr,mod_vad,FreeSWITCH,语音识别,呼叫中心中间

6 天前

没读研的大葱 · Query and transform data | Grafana documentation

2 天前

重感情的斑马 · 3DMax打不开代码错误无法初始化怎么办？ -BIM免费教程_腿腿教学网

1 月前

一身肌肉的烤红薯 · 北京市体育局2018年政府信息公开工作年度报告_市政府部门_首都之窗_北京市人民政府门户网站

2 月前

打盹的板栗 · 社会艺术水平考级管理办法(文化部令第31号)

2 月前

魁梧的饭盒 · 上海外滩踩踏当晚官员“吃喝门”真相追踪-新华网

8 月前

直爽的黄花菜 · 一个方法自测血管是否硬化！这8个指标控制好，老来生病少_手机新浪网

9 月前

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

2 个回答