三次元(XYZ)を1つの関数に置く
Show older comments
図1
図2
図1の縦方向をY軸、横方向をX軸とし、CXVファイルをZ軸方向としています。
この(XYZ)を使ってA=(XYZ)のような１つの関数に置きたいです。
宜しくお願いします。
Accepted Answer
下記にサンプルコードを作成しました。全てのCSVファイルのX,Yのそれぞれのデータ数は全部等しいことを前提にしています。
% サンプルとなるCSVファイルを作成
fileName = 'file_';
for i = 1 : 10
data = zeros(11,11);
data(1,2:end) = 0:9;
data(2:end,1) = 0:9;
data(2:end, 2:end) = reshape(1:100, 10, 10); % サンプルデータはとりあえず連番データを挿入
fileData = num2cell(data);
fileData{1,1} = 'y/x';
writecell(fileData, [fileName, num2str(i-1), '.CSV']);
end
fileData % CSVに保存したデータのひとつを表示
list = dir('*.CSV'); % CSVファイルの取り出し
filePath = fullfile(list(1).folder, list(1).name); % 一つ目のCSVファイルパスを取得
file1Data = readmatrix(filePath, 'Range', [2,2]); % データを示すセル（2行2列目）からデータの読み込み
csvData = zeros([size(file1Data), length(list)]); % データの初期化
csvData(:,:,1) = file1Data; % ひとつめのファイルのデータを挿入
for i = 2 : length(list) % 2ファイル目以降のファイル読み込みのためのループ文
filePath = fullfile(list(i).folder, list(i).name);
csvData(:,:,i) = readmatrix(filePath, 'Range', [2,2]); % CSVファイルのデータの代入
end
csvData(:,:,1) % 一つ目のファイルのデータを表示
9 Comments
皓
on 24 Jul 2023
ご回答ありがとうございます。
2つ目のコードを書くと以下に添付している代入を実行できませんと出てきます。
似たような質問内容があったのですが、良く分からないです。
よろしくお願いいたします。
交感神経優位なあかべぇ
on 24 Jul 2023
エラーが出ているのは「csvData(:,:,i) = readmatrix(filePath, 'Range', [2,2]); % CSVファイルのデータの代入」の行だと思いますが、これは、readmatrixで読み取ったファイルにデータがなく（サイズ0x0）、左辺(csvData(:,:,i) (csvData(1:10,1:10,i)と同意))で要求しているデータサイズ10x10が満たせていませんと言っています。
おそらく、エラーが出たcsvファイルの中身が空だったりすると予想しますので、コードをスクリプトファイル(.mファイル)上に記述して、実行して、ブレークポイントを貼ったりして、エラー原因となったcsvファイルを特定してみるとよいと思います。
皓
on 25 Jul 2023
上記のサンプルだと、csvData(:,:,1)の1つ目の:がY軸、2つ目の:をX軸になっていますので、サンプルに転置処理を加えて、csvData(:,:,1)の1つ目の:がX軸、2つ目の:をY軸に変換させた変数csvDataを作成しました。
また、csvData(:,:,1)のコンター図を描画しました。
% サンプルとなるCSVファイルを作成
fileName = 'file_';
for i = 1 : 10
data = zeros(11,11);
data(1,2:end) = 0:9;
data(2:end,1) = 0:9;
data(2:end, 2:end) = reshape(1:100, 10, 10); % サンプルデータはとりあえず連番データを挿入
fileData = num2cell(data);
fileData{1,1} = 'y/x';
writecell(fileData, [fileName, num2str(i-1), '.CSV']);
end
fileData % CSVに保存したデータのひとつを表示
list = dir('*.CSV'); % CSVファイルの取り出し
filePath = fullfile(list(1).folder, list(1).name); % 一つ目のCSVファイルパスを取得
fileData = readmatrix(filePath, 'Range', [2,2]); % データを示すセル（2行2列目）からデータの読み込み
file1Data = fileData'; % 転置をして、(:,:,1)のひとつめをX軸、２つめをY軸となるように変換する。
csvData = zeros([size(file1Data), length(list)]); % データの初期化
csvData(:,:,1) = file1Data; % ひとつめのファイルのデータを挿入
for i = 2 : length(list) % 2ファイル目以降のファイル読み込みのためのループ文
filePath = fullfile(list(i).folder, list(i).name);
fileData = readmatrix(filePath, 'Range', [2,2])'; % CSVファイルの読み取り
csvData(:,:,i) = fileData'; % 転置をして、(:,:,1)のひとつめをX軸、２つめをY軸となるように変換する。
end
csvData(:,:,1) % 一つ目のファイルのデータを表示
contourf(0:size(csvData, 1)-1, 0:size(csvData, 2) - 1, csvData(:,:,1)); % コンター図の描画 （X,YデータはCSVのデータにあわせて0スタート）
colorbar; % カラーバーも一緒に描画
皓
on 27 Jul 2023
コンター図を描くことができました。
ありがとうございます。
①listを1,2,3,100と小さい順にしたいのですが、下図のように100を超えると1の後ろに来てしまいます。何か方法はありますでしょうか？
➁csvData = zeros([size(file1Data), length(list)]); % データの初期化
csvData(:,:,1) = file1Data; % ひとつめのファイルのデータを挿入
初心者なので、あまり分かっていないのですが、一回データを初期化してからファイルのデータを挿入しないといけないのでしょうか？プログラミングが上手く行かないのかとは思いますが、なぜうまくいかないのかが分からないです。
➂ list(1).folder, list(1).nameの中のlist(1)の番号は何に反映されるのでしょうか？
番号変えた場合でも同じ結果になったので、どうなっているのか分からなかったです。
① ファイル名の数値を抜け出し、昇順に入れ替えるといった方法になると思います。
下記にlist構造体を数値の昇順になるように入れ替えるサンプルを書きました。
% サンプルとなるCSVファイルを作成
fileName = 'file_';
for i = 1 : 10
data = zeros(11,11);
data(1,2:end) = 0:9;
data(2:end,1) = 0:9;
data(2:end, 2:end) = reshape(1:100, 10, 10); % サンプルデータはとりあえず連番データを挿入
fileData = num2cell(data);
fileData{1,1} = 'y/x';
writecell(fileData, [fileName, num2str(i), '.CSV']);
end
list = dir('*.CSV');
{list.name}' % 全ファイル名を表示。2番目のファイルがfile_10.CSVになってしまっている。
allnames = {list.name};% list構造体から全ファイル名を取得
numberStrCell = regexp(allnames, '^file_(\d+)\.', 'once', 'tokens'); % ファイル名のfile_〇〇. の〇〇（数字）を抽出
numberStr = [numberStrCell{:}];% セル配列の中にセル配列の中に文字列というデータ型になっているので、最初のセル配列の中身を出す。
number = str2double(numberStr);% セル配列の文字列を数値配列に変換する。
[~, sortIdx] = sort(number);% 数値を昇順に入れ替えを行ったそのインデックスを取得
list = list(sortIdx); % list構造体を数値の昇順に入れ替え
{list.name}' % 全ファイル名を表示。ファイル名が昇順に入れ替わった。
② 初期化せずに実行することも可能です。ただし変数の配列数がひとつずつ追加されている形になるので、実行速度が遅くなります。（MATLABで下記サンプルのコードを書くと警告が出ます。）
data = rand(100,100);
% 変数の初期化なしの実行速度
tic;
for i = 1 : 1000
test1(:,:,i) = data;
end
toc
% 変数を初期化した場合の実行速度
test2 = zeros(100, 100, 1000);
tic;
for i = 1 : 1000
test2(:,:,i) = data;
end
toc
③ listは構造体配列になっており、list(1)は構造体配列の1番目の要素を差しています。今回はたまたまfolderメンバーは全て同一の値のため、差がないように見えてしまったようです。下記のようになっています。
配列の中身は、ワークスペースの変数名をダブルクリックするか、コマンドで「open('list')」のように打つと、下記図のようなダイアログが現れて変数の詳細を確認できます。
皓
on 28 Jul 2023
①昇順に入れ替えることができました。
➁➂も分かりやすく解説してくださり、ありがとうございます。
④ ①のnumberStrCell = regexp(allnames, '^file_(\d+)\.', 'once', 'tokens'); % ファイル名のfile_〇〇. の〇〇（数字）を抽出でのコードで分からないところがあります。% ファイル名のfile_〇〇. の〇〇（数字）を抽出と書いてくださっているので、何をこのコードで実行しているのかは分かるのですが、'^file_(\d+)\.', 'once', 'tokens'が何を表しているのかが分からないです。
⑤ ①の[~, sortIdx] = sort(number);% 数値を昇順に入れ替えを行ったそのインデックスを取得のコードでは何を行っているのかが分からないです。配列を変換して昇順に順番を入れ替えているということでしょうか？
⑥colorbarのラベルと色がcsvファイルによって図１のようにラベルの数値と色が変わってしまうので、毎回同じ色と数値にしていただきたいです。数値は0,10,20,30,40,50,60,70,80,90,100で色は図2が理想形です。
図１
図2
④ '^file_(\d+)\.'は正規表現です。^は文字列の先頭を指し、\dは数字、\d+は一つ以上連続した数値の文字列、\.は、.を指します。（正規表現で.は任意の文字を表すので、文字の.を表現するには\.になります。）
また、(\d+)の()はトークンを表します。
a = [3,2,1,5,4];
[data, sortIdx] = sort(a) % 出力１は昇順に並び替えたデータで、出力2はソート結果の並び替えインデックス
また、配列データのインデックスに並び替えたい数値配列を指定すると、その配列順にデータを並び替えることができます。
orderIdx = [3,2,1];
data = [1,2,3];
orderData = data(orderIdx)% orderIdxで指定した順序に並び替え
この２つを組み合わせると、sortで指定した並び替え順序を他のデータに反映させることができます。
data = [3,2;1,2;2,7;6,8;5,11;4,100]%１列目の数値が昇順になるように行全体を並べ替えたい。
[~, sortIdx] = sort(data(:,1));% 1列目の数値配列のソート結果のインデックス取得
sortData = data(sortIdx,:) %ソート結果のインデックスをもとに行の並び替え
⑥ カラーのラベルについて、climでカラーの上下限値を設定できるので、全てのグラフでclim([0,100])を指定すれば、全て0~100の幅を持ったカラーバーに統一できます。(MATLAB R2022aより前のバージョンの場合はcaxisです。)
サンプルを下記に書きました。
sample1 = peaks(0:9) .* 10; % サンプルデータ1の作成
[min(sample1(:)), max(sample1(:))] % data1の上下限の確認
sample2 = sample1 .* 2;% サンプルデータ2の作成。サンプル1より2倍大きい
colormap(turbo()); % カラーマップの色の指定 理想形に近そうな色を選択
contourf(0:9,0:9,sample1); % sample1のコンター図描画
clim([0,100]); % カラーマップの上下限値の指定
colorbar; % カラーバー表示
contourf(0:9,0:9,sample2);% sample2のコンター図描画
clim([0,100]);% カラーマップの上下限値の指定
colorbar;% カラーバー表示
皓
on 2 Aug 2023
More Answers (0)
Categories
Find more on cell 配列 in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
