iso_s_smooth.sci

来自「用来实现三维阻抗及光学断层成像重建的matlab程序」· SCI 代码 · 共 106 行

function [sm2]=iso_s_smooth(%simp,w)
sm2=[];
[nargout,nargin] = argn(0)
//function [sm2] = iso_f_smooth(simp,w);
// 
//Calculates a second order discrete Gaussian smoothing operator
// 
// 
// 
//simp = The simplices matrix.
//w    = smoothing weight, w=1...k, default value = 1
//sm2  = The second order smoothing regulariser. 
 
 
if nargin<2 then
  w = 1;
end
 
if w<0 then
  error('Weight must be possitive');
end
 
 
ndsrch = [%simp,(1:size(%simp,1))'];
 
//!! Unknown function spalloc ,the original calling sequence is used
sm2 = spalloc(size(%simp,1),size(%simp,1),80*size(%simp,1));
 
 
for i = 1:size(ndsrch,1)
   
  t_id = ndsrch(i,1:size(%simp,2));
   
  X2simp = [];
  Xsimp = [];
   
   
  //! unknown arg type, using mtlb_length 
  for j = 1:mtlb_length(t_id)
     
    t_nd = t_id(j);
     
    X1 = find(%simp(:,1)==t_nd)';
    X2 = find(%simp(:,2)==t_nd)';
    X3 = find(%simp(:,3)==t_nd)';
    X4 = find(%simp(:,3)==t_nd)';
     
     
    //!! Unknown function unique ,the original calling sequence is used
    XX = unique([X1;X2;X3;X4]);
     
    Xsimp = [Xsimp;XX];
  end
   
  Joined = %simp(Xsimp,:);
   
   
  //!! Unknown function unique ,the original calling sequence is used
  Jr = unique(Joined);
   
   
  //! unknown arg type, using mtlb_length 
  for k = 1:mtlb_length(Jr)
     
    t_2nd = Jr(k);
     
    XX1 = find(%simp(:,1)==t_2nd)';
    XX2 = find(%simp(:,2)==t_2nd)';
    XX3 = find(%simp(:,3)==t_2nd)';
    XX4 = find(%simp(:,4)==t_2nd)';
     
     
    //!! Unknown function unique ,the original calling sequence is used
    XXX = unique([XX1;XX2;XX3;XX4]);
     
    X2simp = [X2simp;XXX];
     
  end
   
   
  //!! Unknown function setdiff ,the original calling sequence is used
  Xdif = setdiff(Xsimp,i);
   
  //!! Unknown function setdiff ,the original calling sequence is used
  Xdif2 = setdiff(X2simp,i);
   
  sm2(i,Xdif) = -0.5*w;
  sm2(i,Xdif2) = -1/w;
   
  //! mtlb_sum(sm2(i,:)) may be replaced by 
  //!     sum(sm2(i,:)) if sm2(i,:)is a vector
  //!     sum(sm2(i,:),1) if sm2(i,:)is a matrix
  sm2(i,i) = abs(mtlb_sum(sm2(i,:)));
end
 
 
 
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
// This is part of the EIDORS suite.
// Copyright (c) N. Polydorides 2001
// Copying permitted under terms of GNU GPL
// See enclosed file gpl.html for details.
// EIDORS 3D version 1.0
// MATLAB version 5.3 R11
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%