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?:）?