[顯示] [隱藏]

%function 3D brain Color
np = 128;                 
kMax = 16;      
X = zeros(kMax, np+1);       %設0矩陣x=kmax，ｙ=mp+1
Y = zeros(kMax, np+1);       %設0矩陣x=kmax，ｙ=mp+1
Z = zeros(kMax, np+1);       %設0矩陣x=kmax，ｙ=mp+1
X(:,:) = 0;                  
Y(:,:) = 0;
Z(:,:) = 0;
%surf (X,Y,Z)               
for k=1:kMax                  %for迴圈 設k=1到kmax
    bFirst  = 1;              %設bFirst=1,bFirst是個變數
    fn = sprintf ('%02d.mat', k);
    %Img = imread (fn);      
    load(fn);                 %讀取(fn)
    Img = new;                %設Img為new   
    [width,height, layers] = size(Img); %以Img去設定width,height, layers
    Img = roicolor(Img, V1, V2);      
    [C, h] = imcontour(Img,1,'b-');
    %  pause
    [num, n] = size(h);       %以h去設定num, n
    n = 1;                    %設n=1
    max = 0;                  %設max=0
    x=[];                     %設x為空矩陣
    y=[];                     %設y為空矩陣
    for i=1:num               %for迴圈 設定i=1到num
        if (C(2,n)>60 & C(1, n+1) > 10 & C(2, n+1)>10) %如果此行程式成立,將可繼續往下跑x,y
            x = [x,C(1, n+1:n+C(2,n))];
            y = [y,C(2, n+1:n+C(2,n))];
            max = max + C(2,n);
        end
        n = n + C(2,n) + 1;
    end
    subplot(1,2,1)
    plot(x,y);                %把xy連成一條線
    for i=1:np/2              %for迴圈 設定i=1到np/2
        xMin = 1000;          %設定x的最小值為1000     
        xMax = -1000;    %設定x的最大值為-1000  
        ys = height/(np/2-1) * i; %設變數為ys = height/(np/2-1) * i
        for j=1:max           %for迴圈 設j=1到最最大值
            if (abs(ys-y(j)) < 1)   %abs是絕對值或複數長度,假如程式計算值小於1,程式將會往下跑
                if (xMin > x(j))    %xMin大於x(j)
                    xMin = x(j);    %則xMin=x(j)
                end
                if (xMax < x(j))    %xMax小於x(j)
                    xMax = x(j);    %則xMax等於x(j)
                 end
            end
        en [ 瀏覽完整內容請先註冊或登入會員。]

這位大哥
你前不接後的要人修改
還真不知是要人改啥呀

呵呵~你是要寫矩陣轉換嗎?還是...說明一下囉!!

謝謝大大分享