使所有峰在Matlab中清晰可见

Question

我终于在这里用lennon310解决了我的问题。 我在时频图片中有数千个细峰的图片。 我无法在一张图片中同时看到所有图片。 根据我时间窗口的物理宽度，一些窗口会出现，而有些窗口会变得可见。

我的数据的图片，我通过绘制imagesc

所有图片均来自相同的数据点T，F，B 。

如何在Matlab的图片中一次绘制所有峰？

Answer 1

您需要使用重新采样来调整图像的大小，以防止出现混叠效果（这种重叠被描述为不可避免的）。

例如，MATLAB imresize函数可以执行抗锯齿。 不要使用“最近”调整大小方法，这就是您现在所拥有的。

Answer 2

@BenVoigt答案的扩展

我的尝试

B = abs(B);
F1 = filter2(B,T); % you need a different filter kernel because resolution is lower    
T = filter2(B,F);
F = F1;
image1 = imagesc(B);
display1 = imresize(image1, [600 600], 'bilinear');   
imshow(T*t, F*fs, display1); 

哪里有问题。

我再次得到分辨率不够的图片

Answer 3

BenVoigt答案的第二次扩展

我对一个内核滤波器的建议是相对随机误差的卷积

data(find(data ~= 0)) = sin(pi .* data(find(data ~= 0))) ./ (pi*data(find(data ~= 0)));        
data(find(data == 0)) = 1; % removing lastly the discontinuity
data1 = data + 0.0000001 * mean(abs(data(:))) * randn(size(data));
data = conv(data, data1);

这是BenVoigt内核过滤器进行分发的意思吗？

这给出了类似的结果

分辨率仍然是一个问题。 如果调整窗口大小，则中心峰往往会轻易繁殖。

我在上图中激活了旧代码，但没有改变结果。 上面的代码对于显示的内核过滤器仍然不够。 可能某些滤波器功能仍必须分别应用于时间和频率轴，例如：

F1 = filter2(B,T); % you need a different kernel filter because resolution is lower    
T = filter2(B,F);
F = F1;

这些过滤器将两个轴上的值弄乱了。 我需要更好地了解他们才能解决此问题。 但是首先要了解它们是否是正确的方法。

该图仍已调整大小。 数据的大小是5001x1的两倍，而F和T的大小是13635x1的两倍。 所以我认为我应该在设置轴，标签和标题后最后调整大小

imresize(image, [13635 13635], 'bilinear');   

因为分布是双线性的。

Answer 4

BenVoigt答案的第三次扩展

我现在画图

imagesc([0 numel(T)], [0 numel(F)], B);              

我的图片中出现严重的锯齿问题。 可能这样的事情应该是操纵时间频率表示

T = filter2(B,t); % you need a different filter kernel because resolution is lower    
F = filter2(B,fs);

Answer 5

BenVoigt答案和评论的第四个扩展

我删除了过滤器并添加了随机相对误差。 我设置T，F，B，和运行的大小

imagesc([0 numel(T)], [0 numel(F)], B, [0 numel(B)])

我仍然出现明显的锯齿，但图片不同

Answer 6

这并不是轻率的，但我认为唯一的方法是获得一个具有更高像素密度的更宽的显示器。 MATLAB和您的视频卡只能在屏幕上显示这么多像素，并且必须决定显示哪个像素，而忽略哪个像素。 数据仍然存在，只是没有显示。 由于您有很多非常狭窄的线条，因此在决定要点亮的像素时，其中一些会被跳过。 更改时间窗口的大小会更改哪些行被混淆，哪些行足够幸运地出现在屏幕上。

我的建议没有特别的顺序：沿时间轴用高斯对线进行卷积以加宽它们，从而增加部分峰出现在屏幕上的可能性。 在600 dpi打印机上将它们打印出来，看它们是否出现。 绘制多个图，每个图都放大一个单独的时间窗口。

Answer 7

BenVoigt答案的第三次扩展

我试图告诉imagesc制作图像的大小，而不是让监视器使用大小，因此它不会丢失数据。

imagesc(T*t, F*fs, B, [50 50]);           
% imagesc(T*t, F*fs, B');

% This must be last
imresize(image, 'bilinear', [64 30],);   

我不知道我应该在哪里指定x轴和y轴的像素量，无论是在命令imagesc或imresize 。

在这里调整图像大小的正确方法是什么？