Matlab Local RX 探测算法代码

function result = Local_RX（Data win_out win_in）
% Hyperspectral Local RX anomaly detector
%   Compiled by Zephyr Hou
%
% Usage
%   [result] = hyperRxDetector（Data window lambda）
% Inputs
%   Data - 3D data matrix （num_row x num_col x num_dim）
%   window - spatial size window （e.g. 3 5 7 9...）
%   lambda - regularization parameter
% Outputs
%   result - Detector output （num_row x num_col）
%  

[a b c] = size（Data）;
result = zeros（a b）;
t = fix（win_out/2）;
t1 = fix（win_in/2）;
M = win_out^2;

% padding avoid edges
DataTest = zeros（3*a 3*b c）;
DataTest（a+1:2*a b+1:2*b :） = Data;
DataTest（a+1:2*a 1:b :）