遥感数据处理

bil = csvread('BIL.txt');
[h,w] = size(bil);
C = zeros(h/3,w,3);
C(:,:,1) = bil(1:3:3348,:);
C(:,:,2) = bil(2:3:3348,:);
C(:,:,3) = bil(3:3:3348,:);
C1 = cast(C, 'uint8');

bip = csvread('BIP.txt');
[h,w] = size(bip);
C = zeros(h,w/3,3);
C(:,:,1) = bip(:,1:3:w);
C(:,:,2) = bip(:,2:3:w);
C(:,:,3) = bip(:,3:3:w);
C2 = cast(C, 'uint8');

bsq = csvread('BSQ.txt');
[h,w] = size(bsq);
h = h/3;
C = zeros(h,w,3);
C(:,:,1) = bsq(1:h,:);
C(:,:,2) = bsq(h+1:2*h,:);
C(:,:,3) = bsq(2*h+1:3*h,:);
C3 = cast(C, 'uint8');

C = [C3,C1,C2];
[h,w] = size(C);
image([1 w],[1 h],C)