Python 3D图像分割在分水岭距离图中寻找局部峰值
。
首先，让我们了解一下相关的概念和背景知识。
- Python：Python是一种高级编程语言，具有简单易学、可读性强、拥有丰富的第三方库等特点，广泛应用于数据分析、机器学习、人工智能等领域。
- 3D图像分割：3D图像分割是指将三维图像中的不同物体或区域分离出来的过程。它在医学影像处理、计算机视觉等领域有着重要的应用。
- 分水岭距离图：分水岭距离图是一种图像处理技术，用于表示图像中每个像素点到最近的分水岭的距离。它常用于图像分割和边缘检测。
- 局部峰值：局部峰值是指在一定范围内具有最大或最小值的数据点。在图像处理中，局部峰值常用于检测图像中的显著特征或边缘。
接下来，我们将讨论如何使用Python进行3D图像分割，并在分水岭距离图中寻找局部峰值。
- 导入必要的库和模块：
import numpy as np
import scipy.ndimage as ndimage
import scipy.ndimage.filters as filters
import matplotlib.pyplot as plt- 加载3D图像数据：
# 读取3D图像数据
image = np.load('image.npy')- 计算分水岭距离图：
# 计算梯度幅值图
gradient = filters.sobel(image)
# 计算分水岭距离图
distance = ndimage.distance_transform_edt(gradient)- 寻找局部峰值：
# 使用最大值滤波器寻找局部峰值
local_maxima = filters.maximum_filter(distance, size=3) == distance
# 使用掩码过滤非局部峰值
background = (distance == 0)
local_maxima &= ~background
# 获取局部峰值的坐标
coordinates = np.transpose(np.nonzero(local_maxima))- 可视化结果：
# 可视化分水岭距离图和局部峰值
fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(10, 5))
ax1.imshow(distance, cmap='gray')
ax1.set_title('Distance Map')
ax2.imshow(image, cmap='gray')
ax2.scatter(coordinates[:, 1], coordinates[:, 0], c='r', s=5)
ax2.set_title('Local Maxima')
plt.show()以上代码演示了如何使用Python进行3D图像分割，并在分水岭距离图中寻找局部峰值。你可以根据实际需求进行调整和优化。
推荐的腾讯云相关产品和产品介绍链接地址：
- 腾讯云图像处理：https://cloud.tencent.com/product/imagex
- 腾讯云人工智能：https://cloud.tencent.com/product/ai
- 腾讯云云服务器：https://cloud.tencent.com/product/cvm
- 腾讯云数据库：https://cloud.tencent.com/product/cdb
- 腾讯云存储：https://cloud.tencent.com/product/cos
- 腾讯云区块链：https://cloud.tencent.com/product/baas
- 腾讯云物联网：https://cloud.tencent.com/product/iot
- 腾讯云移动开发：https://cloud.tencent.com/product/mpp
- 腾讯云音视频：https://cloud.tencent.com/product/vod
- 腾讯云网络安全：https://cloud.tencent.com/product/ddos
相关·内容
1回答
用粗“线”分割三维形状
、、、
我正在寻找一种方法，寻找形状在三维图像中的matlab。我现在没有真正的3D样本图像，事实上，我的3D图像实际上是一组量化的2D图像。
下面的数字是我正在努力实现的目标：
虽然上面的例子是一个2D图像，但请理解我正试图在3D中这样做。输入的形状有这些“触角”，我必须在其中寻找不规则的形状。触角从一个点到另一个点的大小可以改变，但以“一致和平稳”的速度-也就是说，它可以是大的开始，然后逐渐变小。但是如果突然间，形状变得不那么大了，就像上图中右下角的红色区域，那么这就是兴趣的体积之一。请注意，这些形状有更多的倾向是圆形和球形，但其中有些是完全任意和随机的。
到目前为止，我已经尝试了以下方
浏览 0提问于2011-08-15得票数 2
回答已采纳
3回答
我很难找到这张图中的4个主要峰值
信号是值非常不稳定，因为它们先升后降，很难找到最大值和它的索引。
function [peaks, locations] = findMaxs (mag, threshold)
len = length(mag);
prev = 1;
cur = 2;
next = 3;
k = 1; %number of peaks
while next < len
if mag(cur) - mag(prev) > threshold
if mag(cur) > mag(next)
peaks(
浏览 5提问于2014-04-30得票数 2
回答已采纳
这个图 (fig_1)是三维中的单变量正态分布，对应的方程位于左上角。
📷
看来x从-5到5不等；西格玛从0到4不等。
这个Python代码试图再现这个图形。
ax = plt.figure().gca(projection='3d')
xx, yy = np.meshgrid(np.arange(-5,5,.1), np.arange(0.1,4,.1))
zz = norm.pdf(xx,0,yy)
ax.plot_surface(xx, yy, zz, alpha=.5)
我得到了这个数字
📷
而这家伙
zz = norm.pdf(xx,0,np.sqrt(yy))
给出
浏览 0提问于2019-10-15得票数 0
回答已采纳
2回答
我想要检测下图中显示的点：
到目前为止，我已经这样做了：
[X,map] = rgb2ind(img,0.0);
img = ind2gray(X,map); % Convert indexed to grayscale
level = graythresh(img); % Compute an appropriate threshold
img_bw = im2bw(img,level);% Convert grayscale to binary
mask = zeros(size(img_bw));
mask(2:end-2,2:end-2) = 1;
img_bw(m
浏览 3提问于2017-07-20得票数 2
1回答
我正在尝试创建一个像这样的3d曲面图，链接在这里：
但问题是，我只有有限的数据可用，我只有峰位置和峰高的数据，而其余的数据都丢失了。在示例中，z-data需要25X25值625个数据点才能生成有效的曲面图。
我的数据看起来像这样：
所以我的问题是，有没有可能使用一些带有峰值位置值的多项式函数作为约束，以基于我所拥有的信息来生成Z数据？
欢迎任何讨论。任何形式的建议都是受欢迎的。
浏览 8提问于2021-11-23得票数 1
2回答
我正在尝试实现荧光成像中神经元/轴突数据的3D分割和跟踪算法。我研究了分水岭算法，活动轮廓，OTSU阈值，但无法决定采用哪种算法。此外，在不同的帧之间进行跟踪的最佳方式是什么，尤其是这些细胞体移动，有时甚至相互重叠？
我将使用python库来完成这项工作。我正在考虑使用Mahotas来实现，因为它有一个高效的C++实现。此外，我将使用vispy来进行可视化。非常感谢!
浏览 1提问于2015-02-23得票数 1
2回答
科学色谱的图像分析
、、、
我正在寻找一种从GC或HPLC色谱图中提取信息的方法。色谱图如下所示：
我对图像处理/分析不是很感兴趣，所以我正在寻找一种工具/算法来从这些色谱图中提取峰的长度(如果可能的话，还有峰下的表面)。解决方案可以是Python，也可以是C#。
提前谢谢。
浏览 0提问于2012-04-29得票数 1
回答已采纳
2回答
我试着把连接的物体分开。Python和分水岭算法(scipy实现)似乎非常适合处理这一问题。
这是我的图像和自动生成的分水岭种子点(阈值和距离转换图像的局部最大值)：
seeds = myGenSeeds( image_grey )
到目前为止，一切都很好，每个物体都有一种种子。
不过，当我管理分水岭的时候，事情就坏了：
segmented = ndimage.measurements.watershed_ift( 255 - image_grey, seeds)`
中上部星系团和中央星系团的分离程度很低。在最上面的星系团中，有一个物体围绕着另两个物体。在中心星系团中，虽然
浏览 10提问于2013-10-15得票数 6
1回答
我正在寻找Python绘图中的最小点。是否有相反版本的“从scipy.signal导入find_peaks"？如何找到python上图的局部极小值的x点？
浏览 3提问于2020-04-22得票数 2
回答已采纳
1回答
我有一个信号，我想用MATLAB找到它的局部最大值，找到它们的振幅，从最小到最大排序，然后选择集合中的第75个百分位数。
我该怎么做呢？
** 75%百分位数:表示集合中75%的数据等于或小于此值的具体元素。例如:如果我们有从0到100的值，75是第75个百分位数。
浏览 0提问于2012-09-12得票数 0
回答已采纳
我对蟒蛇不熟悉，所以我就..。
我想要制作一个python脚本，它允许我在这样的图像上分离相邻的粒子：
进入像这样的独立区域：
有人建议我使用分水岭法，据我所知，这种方法会给我带来这样的结果：
编辑实际上发现这是距离变换而不是分水岭。
然后我可以用一个门槛把它们分开。跟随这个，但它只起作用，去除粒子。无法“转换”代码来做我想做的事。
然后我又采取了另一种方法。尝试使用，它给了我很好的粒子轮廓。然后，我一直在紧张地寻找一种简单的方法来执行多边形偏移，以便像这样缩小边缘：
使用偏移等高线(多边形)的中心应该给出粒子的数目。但是我只是找到了一种简单的方法来用p
浏览 4提问于2012-10-07得票数 8
回答已采纳
2回答
我是从这里来的：
我还在麻省理工学院开放式频道观看了youtube 。
根据我的理解，该算法是寻找一个局部峰值。
这篇文章的顶部评论说，峰值不一定是全球性的maximum...but，是不是全球最大峰值？这不是真的吗？
如果我从这个算法开始，那么初始的中间元素将决定如何找到峰值。另一边(本地或全球)的任何峰值都将被忽略。
也是在讲座中说，“如果它存在”。那是什么意思？它如何依赖于>=标志？
请指出我错的地方。
浏览 2提问于2016-05-06得票数 0
回答已采纳
我正在试着组合一个简单的程序，让我可视化一系列连续的切割在木板上使用一个特定的切割头的路由器。我正在尝试寻找一个几何库，它能够通过扫描其他3D实体来对3D形状或实体进行减法修改。
具体地说，我希望能够定义矩形实体(木板)，然后定义钻头轮廓形状，并通过矩形实体进行切割(有时在直线上，有时在圆弧上)。有谁知道有什么可以做到这一点的吗？
浏览 0提问于2010-04-10得票数 2
3回答
鞍顶算法
、
我有许多建筑物的地图(openstreetmap项目)。每座建筑都是一个多边形。如何为每个建筑轮廓创建鞍形屋顶零件多边形？
算法应将2D中的一个多边形转换为2D (或3D)中的一组多边形。
这种转换的原因是可视化-更好地渲染等轴测视图。
例如(着色并不重要)：
谢谢
浏览 0提问于2010-03-28得票数 9
回答已采纳
1回答
我有一个3D图，如下所示：
📷
参考文献:谷歌图片
它有两个角度，如X和Y，Z轴是振幅值(每个3D图代表一个像素)。我想把它建模成一些有用的数据结构，比如图表或向量，考虑从上面的3D图中提取的一些参数，这样我就能够将它输入到分类算法中。但是，我无法提取所有的局部最小值/最大值，或斜坡。我如何在python中使用它？我并不是在要求代码(库和代码肯定会受到赞赏)，而是所使用的方法。
我可以使用机器学习从图形中提取某些参数吗？
我对这个问题非常陌生，我甚至不知道如何准确地构建这个问题，所以我为低质量的问题道歉。请指给我看什么东西，这样我就可以从那里看和读。
浏览 0提问于2021-07-07得票数 0
我正在尝试在MATLAB中估计信号的振幅：
可以看出，我不能简单地使用max()和min()来获得振幅。信号不断地游荡，而且在任何情况下都不是恒定的。然而，我所寻找的只是一个平均值--一个单一的数字，用来估计上图中主要分量的平均振幅(重复8次)。从视觉上看，我想要找出的东西是相当明显的。然而，从程序上讲...
我听说过‘希尔伯特变换’--希尔伯特()--和‘快速傅立叶变换’-- fft() --但对信号分析知之甚少--不知道有没有人能把我引向正确的方向，或者解释一下如何利用这两个函数。我的数据只是一个以1000 at存储在一维数组中的浮点数。
感谢您的帮助！
浏览 2提问于2012-02-13得票数 1
回答已采纳
3回答
在2d矩阵中寻找脉冲峰值(dirac增量)的最佳方法是什么？
更具体地说，我想找出给定图像的谐波频率，因此我需要在图像的绝对值DFT中找到脉冲峰值。
我曾想过使用findpeaks，但没有2d版本。我也看过之前关于使用imdilate和/或imextendedmax查找普通峰值的帖子，但这些帖子可以找到2d矩阵中的所有峰值，而我只对脉冲峰值感兴趣。我相信DSP的人对此有一个通用的方法…
请帮帮我
谢谢
浏览 13提问于2010-11-22得票数 4
回答已采纳
1回答
我有一个传感器，它每分钟输出一个电压输出，就会不断地收集数据(以蓝色显示)。我有一个参考传感器收集数据(以红色显示)，输出在我感兴趣的单位。我感兴趣的是确定一个缩放因子，这样我就可以缩放蓝色传感器的数据来匹配红色传感器的数据。
通常，我会在任何给定的时间对两个传感器的值做一个简单的线性回归，这会给我一个基于回归斜率的缩放因子。然而，我注意到，红色传感器在感知环境变化方面速度较慢，而且可以在6到15分钟后的任何地方出现--这使得回归变得困难，因为在任何特定时间，这两个传感器可能测量的是不同的东西。
我想知道是否有任何类型的曲线拟合，可以进行这样我可以提取一个缩放因子，以便我可以缩放蓝色传感器的数
浏览 2提问于2017-11-15得票数 1
回答已采纳
1回答
假设我们有下面的图片，它代表功率谱图。
我的目标是：
1.检测功率谱图的峰值。
2.检测其频率
首先，我从下面的命令中得到了这张图片
[Pxx,f]=periodogram(B,[],[],100);
plot(f,Pxx);
当B是输入信号，100是采样频率时，我尝试在matlab中使用find峰命令，如下所示
[pxx_peaks,location]=findpeaks(Pxx);
然后找到
f(location)
但是它似乎不适合实际的频率，所以请告诉我如何从给定的峰值找到频率?谢谢。
例子如下：
高峰如下：
0.417543614272817
0.38992218758101
浏览 3提问于2014-03-18得票数 1
回答已采纳
3回答
据我所知，Python没有内置的多边形功能。我想创建一个3D地图，并认为多边形将是最好的方法。
不想重新发明轮子，我做了一些谷歌搜索，发现有很多Python的东西，但我找不到我想要的东西。因此，在我重新发明轮子(或整个轮子)之前，有没有人知道Python的多边形系统？
请注意，它确实需要是3D的(我找到了相当多的2D )。还要注意的是，我对它们的显示不感兴趣，而是对在Python中存储它们和数据结构感兴趣。
谢谢
浏览 0提问于2009-01-11得票数 2
回答已采纳
我正在尝试在几个mp3 (音频文件)中找到单词之间的间隔时间，然后再进一步解释一下我的音频文件，如下所示：
我要求我的实验对象在60秒内生成尽可能多的动物名称。我录下了他们的演讲。因此，所有的mp3文件都是动物名称的序列。(单词，不是连续的句子)
首先，我读取其中一个文件，只需输入以下内容即可绘制图形：
plot(data);
图表如下：
我认为一些局部最大值是动物命名的候选者。然而，我不能计算和得到动物名字之间的间隔时间(可能是准确的)。从峰顶到峰顶就行了？或者其他计算单词之间间隔的方法？谢谢:)
浏览 4提问于2014-01-20得票数 1
在Matlab中实现标记控制分水岭时，我遇到了一个问题。
输入图像是具有两个聚类对象的二值掩模。另一个图像是另一个二进制图像，其中白色区域表示标记。
然后，我尝试使用标记控制分水岭来分割聚类对象。代码如下：
bw=imread('Im.jpg'); % read image
bwMarker=imread('Marker.jpg');% read marker
bw(bwMarker)=0; % impose the marker on the image
L=watershed(bw);% do the watershed
%% plot
rgb
浏览 5提问于2011-09-15得票数 8
回答已采纳
2回答
我正在尝试理解MINPEAKDISTANCE是如何工作的。我返回到文档，但是不太清楚这个参数是如何工作的。
你能稍微澄清一下吗？
谢谢。
浏览 15提问于2014-09-30得票数 0
1回答
我可以检测出彼此分离的矩形。但是，我在接触矩形时遇到了如下问题：
我应该在图像中检测到两个矩形。我像预期的那样使用了findContours，并且我已经尝试了各种modes:CV_RETR_TREE、CV_RETR_LIST。我总是得到最外面的一条轮廓，如下所示：
我已经尝试过使用或不使用精明的边缘检测。我做的事情如下：
cv::Mat element = cv::getStructuringElement(cv::MORPH_RECT, cv::Size(3,3));
cv::erode(__mat,__mat, element);
cv::dilate(__mat,__mat, ele
浏览 20提问于2016-08-17得票数 0
回答已采纳
1回答
我正在寻找Python函数，它可以计算从3D中的一个点(x_0、y_0、z_0)到它的端点(x_1、y_1、z_1)和(x_2、y_2、z_2)定义的线段之间的距离。
对于这个问题，我只找到了2D的解决方案。
有一些解决方案可以在3d中找到从点到线的距离，但不能找到线段，如以下所示：
(照片取自)
浏览 2提问于2019-06-05得票数 4
回答已采纳
我正在寻找一个名为3d doctor的软件的替代品，它可以将一堆2d图像(下图中的步骤1和2)转换为具有良好渲染效果的3d图像(步骤3)。在python中有没有一个包可以做同样的事情？谢谢你的建议！
浏览 0提问于2013-05-24得票数 1
回答已采纳
4回答
我感兴趣的是在直方图中找到大致类似的局部最小值
我希望找到109.258的局部最小值，最简单的方法是确定109.258的计数数是否低于周围某个区间(包括109.258)的平均计数数。对我来说，识别这个区间是最困难的部分。
至于这些数据的来源，它是一个具有100个非均匀宽度的直方图。每个存储箱都有一个值(显示在x轴上)，以及落入该存储箱的采样计数(显示在y轴上)。我想要做的是找到分割直方图的“最佳”位置。作为分类算法的一部分，分割的每一边向下传播二叉树。
我认为我最好的做法是尝试将一条曲线拟合到这个直方图中，使用类似于的东西，然后比较局部最小值以找到“最佳”的一个。适当的“最佳”度量将
浏览 0提问于2011-07-31得票数 12
回答已采纳
在Jeff Edmonds的“如何思考算法”一文中，有一节解释了网络流和线性规划章节中的原始-对偶爬山。我很难想象屋顶的指数数量，以及为什么“最低的、因此最理想的屋顶高于最高的、因此最适宜站立的地方”。
浏览 4提问于2012-07-27得票数 1
1回答
基于颜色的图像分割
、、、、
我一直在寻找使用openCV和python对电阻器频带中存在的不同颜色进行分类的方法。
什么算法可以用来分割图像到不同的波段等。我尝试使用分水岭算法，但我不能得到正确的标记，也没有得到想要的结果。
我使用的示例图像：
浏览 1提问于2017-03-05得票数 1
我有一个场景，我有2个或更多的系列和单个图表区域( 3D样式)。当我为两个系列添加超过1000个数据点时，我无法看到两个系列数据点之间的距离。注意:我尝试了系列的点宽，也添加了边框厚度和颜色，但整个网格挤压。
在这张图中，它有1000个数据点，两个系列看起来都是挤压/重叠的，我需要不同系列的点之间的距离：
此图为两个系列的100个点，它显示了系列之间的数据点之间的一些空间：
浏览 0提问于2018-01-04得票数 0
我有一些代码可以制作图像的直方图。我目前正在使用scipy方法find_peaks来寻找峰值，但它也在检测一些冗余值。下面是我使用的代码： from scipy.signal import find_peaks
...
hist = cv2.reduce(threshed,1, cv2.REDUCE_AVG).reshape(-1)
peaks, _ = find_peaks(hist)
plt.figure(figsize=(20,20))
plt.plot(hist)
plt.plot(peaks, hist[peaks], "x")
plt.show()
... 这是
浏览 15提问于2019-06-26得票数 0
1回答
我有以下问题。我想做一个基于区域的分割，就像在图像分割教程中提供的。本教程非常清楚。但在我看来，他们使用灰色值直方图中最极端部分的部分似乎不清楚。它们使用以下直方图中的值30和150，而不解释如何获得它们：
现在我的问题是如何得到这些外部极端？
浏览 0提问于2018-07-26得票数 1
回答已采纳
2回答
使用体素生成球体
、、
我一直在使用，我是3D编程的新手，它在中说这段代码生成了“球体世界”：
generate: function(x,y,z) {
return x*x+y*y+z*z <= 20*20 ? 1 : 0 // sphere world
},
这实际上是如何生成球体的？从我的简单理解来看，我认为它基本上是在3D世界中的每个“块”中“循环”？任何进一步的解释或指向一个好的教程将是一个巨大的帮助！
浏览 5提问于2013-01-24得票数 3
回答已采纳
2回答
一种新的blob检测算法
、、、、
这是我的实验的一个样本图像，将被处理，以找到所有的气泡(中心，尺寸等)的信息。基本上，下面的图像包含了一些噪音，几条粗体线条，当然也有几个斑点。
我的目标是用一个图像处理软件包"Sapera“编程一个FPGA。Sapera APF是一个用于编程Teledyne DALSA帧采集器来执行图像处理的应用程序。不幸的是，它不允许我们用常规的方式编写代码，也不允许使用基本的编码特性(如果其他语句、循环等)。但是只支持一些内置的代码块形式的函数，我们只能对这些块进行排列。FYI，下面是Sapera支持的函数列表的。
现在，我的第一个目标是在matlab上测试我的算法，然后在Sapera
浏览 6提问于2014-06-24得票数 1
1回答
MATLAB网格图中的峰值分析
、、、、
我正在对下面的图像进行曲率分析，需要能够测量峰值之间的距离。有人能说明一下这是怎么做的吗？我只对图的正部分感兴趣，如果每个峰值的高度也能被提取出来，它会很有用。
我试过使用find峰值，但它要求一个向量，但我没有一个可以给它。
我正在使用来自online存储库的surfature函数。为了解决这个问题，“X”是我输入到高斯滤波器中的图像：
G = fspecial('gaussian',[19 19],6);
X = imfilter(X,G,'same');
[XI YI]=meshgrid(1:size(X,2), 1:size(X
浏览 1提问于2015-02-11得票数 0
回答已采纳
1回答
分水岭分割总是返回黑色图像
、、、、
我最近一直在研究角膜内皮细胞的分割过程，我发现了一篇相当不错的论文，它描述了如何以良好的效果来实现这一目标。我一直试图跟随那篇论文，并使用scikit-image和openCV来实现它，但是我已经陷入了分水岭分割的泥潭。
我将简要介绍这一进程是如何进行的：
首先，您有原始的内皮细胞图像。
然后，他们指导你进行形态学灰度重建，以达到图像灰度的水平(然而，他们并没有解释如何获得灰度的标记，所以我一直在胡闹，并试图在自己的道路上得到一些)。
这就是重建的图像应该是什么样子：
这就是我的重建图像(让我们将其标记为r)的样子：
目的是利用重建图像获取分水岭分割的标志，如何实现？！我们得到原始图像(让我们
浏览 3提问于2017-11-28得票数 2
回答已采纳
2回答
我想过滤一个真正有噪声的信号，以便获得超过某个阈值的峰值数量。这是我使用的matlab代码：
thresh = 3; % Specify threshold
x = VarName1;
% Essentially need to pick out values which exceed threshold with the condition that the previous value
% needs to be below the threshold
idxl = x>=thresh;
idxl(1) = 0;
idx = find(idxl);
yest = x(idx-
浏览 2提问于2018-05-07得票数 0
1回答
我被赋予了一项任务，给出了离散的数据，如下所示
我需要将它分割成5块，由它创建的模板决定。
我不允许猜测模板，因为每个输入看起来都不一样。
我的方法是在数据中找到峰值(高于或低于零)，然后使用这种峰值模式来分割数据。以下是我得到的信息：(以上数据不适用)
上面的图是图中的峰值，因为我知道我有5个部分，15个点，我可以说每个部分都有3个点，然后切片，这是图中的第二个图。
在40个输入中，我成功地完成了其中的5个，因为我的“峰值检测”算法非常基本。
我应该使用什么峰值检测算法，它也可以找到局部极小值，并且有PHP实现/简单的psudo代码？我是这个领域的数据分析的初学者，所以我
浏览 5提问于2016-08-03得票数 0
在混合2D和3D子图时，我遇到了删除多余空格的问题。对于纯3D子图，我可以使用fig.subplots_adjust()调整正在绘制的区域以删除空白，请参见。
然而，如果这个3D图像在2D子图中，同样的技巧就不起作用了。我创建了混合的子图，如下所示：
import matplotlib.pyplot as plt
from matplotlib import cm
from mpl_toolkits.mplot3d import axes3d
fig,axes = plt.subplots(2,2)
ax = axes.flat
for a in range(3):
ax[
浏览 95提问于2017-07-05得票数 7
1回答
我是倒谱分析的新手。所以这就是问题所在。
我有一个长度为4096，采样率为8000 Hz的信号。我做了快速傅立叶变换，得到了长度为4096*2的数组(2*i位置是余弦系数，2*i+1位置是正弦系数)。频率步长为(采样率/信号长度== 8000/4096)。因此，我可以这样计算i位置的频率: i*sampleRate/signalLength。
然后，我进行倒谱变换。我不明白如何找到频率步长，以及如何找到给定频率的频率。
浏览 3提问于2013-08-28得票数 1
现在我使用标准相位相关进行图像拼接。它给出了正常的结果，但在硬图像上给出了错误的结果，但ImageJ(斐济)中的Stitch 2D插件在大多数情况下都给出了很好的结果。本文中描述的插件中使用的算法是，但我不能理解它。然而，在真实图像中，F−1(Q)包含多个峰值，标记出具有高相关性的不同平移。此外，由于傅立叶空间的周期性，每个峰值描述8种不同的可能平移(在3D中)。为了确定正确的平移，我们从F−1(Q)中选择n个最高的局部最大值(3×3×3邻域)，并通过图像A，B的重叠区域上的互相关来评估其8种可能的平移。具有最高相关性的峰值被选为两个图像之间的平移。如果没有一个峰值超过一定的限制，则认为这些平
浏览 0提问于2012-02-17得票数 1
回答已采纳
我试图在Matlab中找到一些峰值，但是函数findpeaks.m没有宽度选项。我想要被检测到的峰值在球里。所有检测到的都在红色方块中。如你所见，它们的宽度很小。有什么帮助吗？
下面是我使用的代码：
[pk,lo] = findpeaks(ecg);
lo2 = zeros(size(lo));
for m = 1:length(lo) - 1
if (ecg(m) - ecg(m+1)) > 0.025
lo2(m) = lo(m);
end
end
p = find(lo2 == 0);
lo2(p) = [];
figure, p
浏览 0提问于2013-07-02得票数 2
回答已采纳
我在python中使用下面的代码绘制了一个3D图。这里使用三个数组x、y和z来绘制这幅图。我想显示图中数组的最后一点(或3D线的结束点)。我使用了在2d绘图中使用的方法，也就是说，我只要求使用这个命令ax.plot(x[-1],y[-1],z[-1],'o')绘制每个数组的最后一点。但不起作用。
import numpy as np
from matplotlib import pyplot as plt
from mpl_toolkits.mplot3d import Axes3D # noqa: F401 unused import
x=np.linspace(0,2*
浏览 5提问于2020-11-23得票数 0
回答已采纳
1回答
我在上礼貌课。我尽了最大努力，得到了40%的分数，然后尝试寻找答案。
能不能请得到以下代码的人向我解释一下maxFlags行发生了什么(为什么maxFlags是那样的，从我理解的最大标志来看，我们可以有等于峰值的数目，然后根据距离进行过滤)
这是我找到的100%的代码
function solution(A) {
if (A.length <= 2) return 0;
var peaks = [];
var size = 0;
for (var i = 1; i < A.length-1; ++i) {
if (A[i] >
浏览 0提问于2020-11-16得票数 1
回答已采纳
1回答
我试图在tensorflow中创建一个3d距离图(大小:W*H* D)，用于训练的损失函数。我有一个地面实况图(大小为W*H*D的二进制体积)，我将使用它来创建距离图，即距离图中每个像素的值将是该像素到地面实况图中正值(即pixel=1)形状的最小距离。对于3d形状问题，L2NORM将轴减少到2D形状，并使此问题完全可微。任何建议或指示都将不胜感激。
浏览 10提问于2019-11-27得票数 2
回答已采纳
1回答
下面的kmeans聚类有什么不同吗？
a)将图像转换为灰度，并对1D特征向量执行kmeans
b)保持3个通道RGB，并对3D特征向量执行kmeans
c)获取图像直方图并对分布执行kmeans
第一种解决方案肯定更快，但会有区别吗？也许某些像素具有不同的RGB，但具有相同的强度？
浏览 0提问于2012-09-25得票数 0
1回答
有没有一个python模块可以集成简单的色谱/痕迹分析算法？我正在寻找基线校正，峰值检测和峰值积分功能的简单时间进程(与数据存储在numpy数组)。
我现在花了相当多的时间寻找，似乎没有任何真正让我惊讶的东西。
浏览 3提问于2017-05-09得票数 0
1回答
我正在尝试分割python中多孔网络的3d断层图像。我可以用ndimage.distance_transform_edt计算距离图，用feature.peak_local_max计算峰值图。当我应用分水岭算法时，得到了一个可以接受的结果，但是峰的标记不在可见峰上，见距离图的图像。
提前感谢
这里的代码a是图像
D = ndimage.distance_transform_edt(a)
localMax = feature.peak_local_max(D, indices=False, min_distance=50,
labels=a)
localMax2
浏览 50提问于2018-08-15得票数 2
我有一个连通部分的图像(用圆圈填充)，我想分割它们，我可以使用分水岭算法。我更喜欢为分水岭编写自己的函数，而不是使用opencv中的内置函数。我已经成功了，我如何使用.If找到对象的区域最大值？
浏览 2提问于2012-06-26得票数 3
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
