How do I create multiple figure of cwt in for loop and save them on folder

2 次查看（过去 30 天）
saratt Ratana-Ubol
评论： saratt Ratana-Ubol 2022-4-8
Hi everyone!
I record data from emg sensor and segment the signal. I want to convert signal in cwt figure and save them as a picture my example file contain 2560*46 (46 figures).
This is my code and file
load Data_segment.mat
%% Normalized frequency (*pi rad/sample)
frequencyLimits = [0 1];
voicesPerOctave = 8;
%% Convert to cycles/samples
frequencyLimits =frequencyLimits/2;
path = 'C:\Users\payusan\Matlab\BB\cwt\AI\Dataset_creation\CWT_creation' ;
myfolder = 'output_CH1' ;
folder = mkdir([path,filesep,myfolder]) ;
path = [path,filesep,myfolder] ;
for k = 1:length(Data_segment(end,:))
for i = 1 : length(Data_segment(:,end))
for j = 1 : length(Data_segment(end,:))
figure(j)
frequencyLimits(1) = max(frequencyLimits(1),...
cwtfreqbounds(numel(Data_sgement(:,i))));
[WT,F]= cwt(Data_segment(:,i), ...
'VoicesPerOctave',voicesPerOctave, ...
'FrequencyLimits',frequencyLimits);
pcolor(1:numel(Data_segment(:,i),log2(F),abs(WT)));
shading interp;
end
end
temp=[path,filesep,'seg_',num2str(k),'.jpg'];
saveas(figure(k),temp);
end
How do I fix this code.
  4 个评论
显示 2更早的评论
KSSV
KSSV 2022-4-8
We cannot run the code as all the variables are not given,
Unrecognized function or variable 'frequencyLimits'.
Error in Junk (line 14)
frequencyLimits(1) = max(frequencyLimits(1),...
Note that pcolor needs matrix as input. In your case it is not, so error popped out.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

KSSV
KSSV 2022-4-8
May be you need th add the below shown line after loading the data from mat file.
load Data_segment
Data_segment = New_data ;
  1 个评论
saratt Ratana-Ubol
saratt Ratana-Ubol 2022-4-8
编辑：saratt Ratana-Ubol 2022-4-8
I already try that one this is my full code of this part.
D = load('C:\Users\payusan\Matlab\BB\cwt\AI\Mat_file_dataset\label_0\2.mat');
% Signal_segmentation
data = D.CH1;
N_max=length(data);
n_window=2560;
overlap_win=-1;
max_count = ceil((length(data)-n_window)/(n_window-overlap_win))+1;
new_data = zeros(max_count,n_window);
n_start=1;
for count = 1:max_count
n_end = n_start + n_window - 1;
if n_end > N_max
new_data(count,1:N_max-n_start+1) = data(n_start:N_max);
else
new_data(count,:) = data(n_start:n_end);
end
n_start = n_end-overlap_win;
end
%% Normalized frequency (*pi rad/sample)
frequencyLimits = [0 1];
voicesPerOctave = 8;
%% Convert to cycles/samples
frequencyLimits =frequencyLimits/2;
%%
Data_segment = new_data';
path = 'C:\Users\payusan\Matlab\BB\cwt\AI\Dataset_creation\CWT_creation' ;
myfolder = 'output_CH2' ;
folder = mkdir([path,filesep,myfolder]) ;
path = [path,filesep,myfolder] ;
for k = 1:length(Data_segment(end,:))
for i = 1 : length(Data_segment(:,end))
for j = 1 : length(Data_segment(end,:))
figure(j)
frequencyLimits(1) = max(frequencyLimits(1),...
cwtfreqbounds(numel(Data_segment(:,i))));
[WT,F]= cwt(Data_segment(:,i), ...
'VoicesPerOctave',voicesPerOctave, ...
'FrequencyLimits',frequencyLimits);
pcolor(1:numel(Data_segment(:,i),log2(F),abs(WT)));
shading interp;
end
end
temp=[path,filesep,'seg_',num2str(k),'.jpg'];
saveas(figure(k),temp);
end

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Continuous Wavelet Transforms 的更多信息

标签

产品


版本

R2020a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by