double to table, create new tabel and add data to existing table with different properties

16 ビュー (過去 30 日間)
Frederik Reese
Frederik Reese 2022 年 5 月 17 日
回答済み: Abhishek Chakram 2023 年 9 月 21 日
Hello,
I have a Matlab table for which I carry out many individual calculations as can be seen in the code. the results are columns and tables in double format. I'm having a hard time putting the data back together in a table. I would like a separate table for BHQ, HQextrem, HQ5000 and HQ 10000 with the following data
  • colums 1:6; PP_D colums 1:6; with column names OID, FKM, Lage, x, y, z
  • colums 7; ...MWH name: Maximalewasserhöhe
  • colums 8; ...MWA name:Maximaler Wasseranstieg
  • colums 9; ...MWR name: Maximaler Wasserrückgang
  • colums 10;...FD name: Maximale Überflutungsdauer
  • colums 11: 152; ...WH names: WH Time 0:142
  • colums 153: 295: ... WA names WA Time 0:142
Thanks for your help
here the code and the data:
%% Bearbeitung Punkte
clc
clear all
clear workspace
load ('PP_D.mat')
PP_D_new = removevars(PP_D, {'OID_','FKM','Lage'});
PP_D_string = table2array(PP_D_new);
PP_D_string(PP_D_string==0)=nan
%% Wasserhöhen
for i=1:151
PP_D_Wasserhoehe (i,:)= PP_D_string([i],[4:513])- PP_D_string(i,3);
end
% PP_D_Wasserhoehe=str2double(PP_D_Wasserhoehe);
D_BHQ_WH= PP_D_Wasserhoehe([1:151],[1:83])
D_HQextrem_WH= PP_D_Wasserhoehe([1:151],[84:225])
D_HQ5000_WH= PP_D_Wasserhoehe([1:151],[226:367])
D_HQ10000_WH= PP_D_Wasserhoehe([1:151],[368:509])
%Maximale Wasserhöhen
for i=1:151
D_HQ10000_MWH (i,1)= max(D_HQ10000_WH(i,[1:142]));
D_HQ5000_MWH (i,1)= max(D_HQ5000_WH(i,[1:142]));
D_HQextrem_MWH (i,1)= max(D_HQextrem_WH(i,[1:142]));
D_BHQ_MWH (i,1)= max(D_BHQ_WH(i,[1:83]));
end
%% Anstieg pro Zeitschritt [m/h]
PP_D_Wasserhoehe(isnan(PP_D_Wasserhoehe))=0;
for i=2:510
PP_D_Wasserhoehe_Wasseranstieg(:,i)= PP_D_Wasserhoehe([1:151],i)- PP_D_Wasserhoehe([1:151],i-1);
end
D_BHQ_WA= PP_D_Wasserhoehe_Wasseranstieg([1:151],[1:83])
D_HQextrem_WA= PP_D_Wasserhoehe_Wasseranstieg([1:151],[84:225])
D_HQ5000_WA= PP_D_Wasserhoehe_Wasseranstieg([1:151],[226:367])
D_HQ10000_WA= PP_D_Wasserhoehe_Wasseranstieg([1:151],[368:509])
%Maximaler Wasseranstieg
for i=1:151
D_HQ10000_MWA (i,1)= max(D_HQ10000_WA(i,[1:142]));
D_HQ5000_MWA (i,1)= max(D_HQ5000_WA(i,[1:142]));
D_HQextrem_MWA (i,1)= max(D_HQextrem_WA(i,[1:142]));
D_BHQ_MWA (i,1)= max(D_BHQ_WA(i,[1:83]));
end
%Maximaler Wasserrückgang
for i=1:151
D_HQ10000_MWR (i,1)= min(D_HQ10000_WA(i,[1:142]));
D_HQ5000_MWR (i,1)= min(D_HQ5000_WA(i,[1:142]));
D_HQextrem_MWR (i,1)= min(D_HQextrem_WA(i,[1:142]));
D_BHQ_MWR (i,1)= min(D_BHQ_WA(i,[1:83]));
end
%% Überflutungsdauer
PP_D_Wasserhoehe(isnan(PP_D_Wasserhoehe))=0;
D_BHQ_x= PP_D_Wasserhoehe([1:151],[1:83]);
D_HQextrem_x= PP_D_Wasserhoehe([1:151],[84:225]);
D_HQ5000_x= PP_D_Wasserhoehe([1:151],[226:367]);
D_HQ10000_x= PP_D_Wasserhoehe([1:151],[368:509]);
for i=1:151
D_HQ10000_FD (i,1)= nnz(D_HQ10000_x(i,[1:142]));
D_HQ5000_FD (i,1)= nnz(D_HQ5000_x(i,[1:142]));
D_HQextrem_FD (i,1)= nnz(D_HQextrem_x(i,[1:142]));
D_BHQ_FD (i,1)= nnz(D_BHQ_x(i,[1:83]));
end

回答 (1 件)

Abhishek Chakram
Abhishek Chakram 2023 年 9 月 21 日
Hi Frederik Reese,
I understand that you are facing difficulty in creating separate tables for BHQ, HQextrem, HQ5000 and HQ 10000 with your the calculated data.
Here’s the example for doing it:
%% Bearbeitung Punkte
clc
clear all
clear workspace
load ('PP_D.mat')
PP_D_new = removevars(PP_D, {'OID_','FKM','Lage'});
PP_D_string = table2array(PP_D_new);
PP_D_string(PP_D_string==0)=nan
%% Wasserhöhen
for i=1:151
PP_D_Wasserhoehe (i,:)= PP_D_string([i],[4:513])- PP_D_string(i,3);
end
% PP_D_Wasserhoehe=str2double(PP_D_Wasserhoehe);
D_BHQ_WH= PP_D_Wasserhoehe([1:151],[1:83])
D_HQextrem_WH= PP_D_Wasserhoehe([1:151],[84:225])
D_HQ5000_WH= PP_D_Wasserhoehe([1:151],[226:367])
D_HQ10000_WH= PP_D_Wasserhoehe([1:151],[368:509])
%Maximale Wasserhöhen
for i=1:151
D_HQ10000_MWH (i,1)= max(D_HQ10000_WH(i,[1:142]));
D_HQ5000_MWH (i,1)= max(D_HQ5000_WH(i,[1:142]));
D_HQextrem_MWH (i,1)= max(D_HQextrem_WH(i,[1:142]));
D_BHQ_MWH (i,1)= max(D_BHQ_WH(i,[1:83]));
end
%% Anstieg pro Zeitschritt [m/h]
PP_D_Wasserhoehe(isnan(PP_D_Wasserhoehe))=0;
for i=2:510
PP_D_Wasserhoehe_Wasseranstieg(:,i)= PP_D_Wasserhoehe([1:151],i)- PP_D_Wasserhoehe([1:151],i-1);
end
D_BHQ_WA= PP_D_Wasserhoehe_Wasseranstieg([1:151],[1:83]);
D_HQextrem_WA= PP_D_Wasserhoehe_Wasseranstieg([1:151],[84:225]);
D_HQ5000_WA= PP_D_Wasserhoehe_Wasseranstieg([1:151],[226:367]);
D_HQ10000_WA= PP_D_Wasserhoehe_Wasseranstieg([1:151],[368:509]);
%Maximaler Wasseranstieg
for i=1:151
D_HQ10000_MWA (i,1)= max(D_HQ10000_WA(i,[1:142]));
D_HQ5000_MWA (i,1)= max(D_HQ5000_WA(i,[1:142]));
D_HQextrem_MWA (i,1)= max(D_HQextrem_WA(i,[1:142]));
D_BHQ_MWA (i,1)= max(D_BHQ_WA(i,[1:83]));
end
%Maximaler Wasserrückgang
for i=1:151
D_HQ10000_MWR (i,1)= min(D_HQ10000_WA(i,[1:142]));
D_HQ5000_MWR (i,1)= min(D_HQ5000_WA(i,[1:142]));
D_HQextrem_MWR (i,1)= min(D_HQextrem_WA(i,[1:142]));
D_BHQ_MWR (i,1)= min(D_BHQ_WA(i,[1:83]));
end
%% Überflutungsdauer
PP_D_Wasserhoehe(isnan(PP_D_Wasserhoehe))=0;
D_BHQ_x= PP_D_Wasserhoehe([1:151],[1:83]);
D_HQextrem_x= PP_D_Wasserhoehe([1:151],[84:225]);
D_HQ5000_x= PP_D_Wasserhoehe([1:151],[226:367]);
D_HQ10000_x= PP_D_Wasserhoehe([1:151],[368:509]);
for i=1:151
D_HQ10000_FD (i,1)= nnz(D_HQ10000_x(i,[1:142]));
D_HQ5000_FD (i,1)= nnz(D_HQ5000_x(i,[1:142]));
D_HQextrem_FD (i,1)= nnz(D_HQextrem_x(i,[1:142]));
D_BHQ_FD (i,1)= nnz(D_BHQ_x(i,[1:83]));
end
% Create separate tables for BHQ, HQextrem, HQ5000, and HQ10000
BHQ_table = table(PP_D.OID_, PP_D.FKM, PP_D.Lage, PP_D.X, PP_D.Y, PP_D.z, D_BHQ_MWH, D_BHQ_MWA, D_BHQ_MWR, D_BHQ_FD);
BHQ_table.Properties.VariableNames = {'OID', 'FKM', 'Lage', 'x', 'y', 'z', 'Maximalewasserhöhe', 'Maximaler Wasseranstieg', 'Maximaler Wasserrückgang', 'Maximale Überflutungsdauer'};
HQextrem_table = table(PP_D.OID_, PP_D.FKM, PP_D.Lage, PP_D.X, PP_D.Y, PP_D.z, D_HQextrem_MWH, D_HQextrem_MWA, D_HQextrem_MWR, D_HQextrem_FD);
HQextrem_table.Properties.VariableNames = {'OID', 'FKM', 'Lage', 'x', 'y', 'z', 'Maximalewasserhöhe', 'Maximaler Wasseranstieg', 'Maximaler Wasserrückgang', 'Maximale Überflutungsdauer'};
HQ5000_table = table(PP_D.OID_, PP_D.FKM, PP_D.Lage, PP_D.X, PP_D.Y, PP_D.z, D_HQ5000_MWH, D_HQ5000_MWA, D_HQ5000_MWR, D_HQ5000_FD);
HQ5000_table.Properties.VariableNames = {'OID', 'FKM', 'Lage', 'x', 'y', 'z', 'Maximalewasserhöhe', 'Maximaler Wasseranstieg', 'Maximaler Wasserrückgang', 'Maximale Überflutungsdauer'};
HQ10000_table = table(PP_D.OID_, PP_D.FKM, PP_D.Lage, PP_D.X, PP_D.Y, PP_D.z, D_HQ10000_MWH, D_HQ10000_MWA, D_HQ10000_MWR, D_HQ10000_FD);
HQ10000_table.Properties.VariableNames = {'OID', 'FKM', 'Lage', 'x', 'y', 'z', 'Maximalewasserhöhe', 'Maximaler Wasseranstieg', 'Maximaler Wasserrückgang', 'Maximale Überflutungsdauer'};
for i = 1:83
columnName = sprintf('WH Time %d', i-1);
BHQ_table.(columnName) = D_BHQ_WH(:, i);
end
for i = 1:83
columnName2 = sprintf('WA Time %d', i-1);
BHQ_table.(columnName2) = D_BHQ_WA(:, i);
end
for i = 1:142
columnName = sprintf('WH Time %d', i-1);
HQextrem_table.(columnName) = D_HQextrem_WH(:, i);
HQ5000_table.(columnName) = D_HQ5000_WH(:, i);
HQ10000_table.(columnName) = D_HQ10000_WH(:, i);
end
for i = 1:142
columnName2 = sprintf('WA Time %d', i-1);
HQextrem_table.(columnName2) = D_HQextrem_WA(:, i);
HQ5000_table.(columnName2) = D_HQ5000_WA(:, i);
HQ10000_table.(columnName2) = D_HQ10000_WA(:, i);
end
You can refer to the following documentation to know more about the functions used:
Best Regards,
Abhishek Chakram

カテゴリ

Help Center および File ExchangeData Type Conversion についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by