matlabでの計算について

Question

matlabで計算をしたいのですが、使い方がいまいちわかっておりません
計算というのは例えば
13
1
15
7
≈
58
576
のように１行６１４４００列の行列のファイルがあり、偶数列目に256をかけてから１列目と２列目、３列目と４列目、５列目と６列目というように機数列目と偶数列目の和を１行３０７２００列の行列にしたいのです。
c言語ならだいたい
[x] = a[2n]+a[2n-1]*256
n = n + 1
的な感じになりそうなんですが、そちらも使い方があまりよく分かっておりません…(´;ω;｀)

ちなみにそこから480×680にして画像データにしたいのですが、どちらか一方でもアドバイスをくだされば幸いです

Answer 1

前の方と同じく１行６１４４００列のデータとして回答します。

x = 1:614400;
% x =
%     1 2 3 4 5 6 7 8 9 10 ...

reshape を使って、隣り合うデータをペアにして列に並べた行列を作ります。

A = reshape(x, 2, length(x)/2);
% A =
%     1  3  5  7  9
%     2  4  6  8 10 ...

行列計算を行い、上の列を256倍して、下の列に足します。

B = [256 1] * A;
% B =
%     258      772      1286     1800     2314 ...
% つまり
%     256*1+2  256*3+4  256*5+6  256*7+8  256*9+10 ...

もう一度 reshape を使って、最後の答えを求めます。

reshape(B,480,640);

Comments

[kosuke1126]

ご回答ありがとうございます。
何とか画像データにもってくことができました。
ご丁寧にありがとうございました

Answer 2

１行６１４４００列

とありますが、数字が縦に並んでいるところをみると
６１４４００行１列
の間違いではないでしょうか？
どちらでもいいのですが、データが１行６１４４００列で与えられているとして話をすすめます。

すでに回答があがっていますが、よりMATLABにふさわしい書き方を紹介します。
データを d に格納します。
しかし、プログラムが意図したとおりに動作するかを簡単にチェックしたいので、
まずは
d = 1:614400;
という自然数が並んだデータを用意しましょう。上記の回答ではd=rand(1,614400);としていますが、これだと結果が正しいかチェックしずらいのでおすすめしません。
（本当は d= 1:10;くらいに設定してテストするのがちょうどよいのですが。。。）
さて、dの偶数番目に256をかけ,dの奇数番目と和をとるには次のようにします。

dd = 256*d(2:2:end) + d(1:2:end);

なぜこれでよいかはマニュアルを読んで自習してください。
d(1:2:end-1)では？と思われるかもしれませんが、d(1:2:end)で十分です。
ちなみにMATLABはd(1)が最初の要素なので上の回答では奇数番目を256倍することになってしまいます。

これを480x640にするにはreshapeで整形します。

ddd = reshape(dd,480,640);

＞480x680
とありますが、サイズが合いません。おそらく480 x 640の間違いでしょう。
もし本当に680にしたいのであればデータ長さが足りないので0をうめることにします。
その場合は
ddd = zeros(480,680);
ddd(1:length(dd))=dd;
でOKです。

C:row-major order
MATLAB:column-major order
にも注意してください。
参考

Comments

[kosuke1126]

睡魔と闘いながら質問したのですごいミスがありましたねｗｗ
ご指摘ありがとうございます。
10万枚くらい画像データがあるので、参考にさしていただきます

Answer 3

配列の変換については、例えば以下のようにすればよいかと思います。

x = randn(1,614400);
newx = zeros(1,307200);
for i=1:25
    newx(1,i)=x(2*i)+x(2*i-1)*256;
end

上記では、xに乱数を入れていますが、ここには、現在お手持ちのデータを入れていただければと思います。
excelデータであれば、例えば以下のようにxlsread関数でデータをそのまま読み込めます。
xxx.xlsxのsheet1を読み込み、dataに格納する場合、

[data]=xlsread('data.xlsx','sheet1')

dataに格納したデータの1行目のデータをxdataに格納したい場合、

xdata=data[1,:]