SAR C矩阵转T矩阵

function T = C2T（C）

% check argument
if nargin ~= 1
    error（‘Invalid argument.‘）;
end

if ~isnumeric（C）
    error（‘Invalid argument.‘）;
end

if ndims（C） ~= 3
    error（‘Invalid argument.‘）;
end

[row col pag] = size（C）;

if （row <= 0 || col <=0） || pag ~= 9
    error（‘Invalid argument.‘）;
end

% convert to coherency matrix
T = zeros（row col 9）;

T（: : 1） = （C（: : 1） + C（: