clear all; close all; clc
num_rows = 256;%行数64
num_cols = 200;%列数64
N = 1;
M = num_rows * num_cols;%每一波段平面,像素个数
X2 =imread('D:\Matlab\RXDDatagroup\飞机5\100.bmp')
X1 = reshape(X2,N,M);
X = im2double(X1);
X_mean = mean(X,2); %按照行求和取平均值(mean本来是按照列求和,转置求均值再转置)
X = X - repmat(X_mean,1,M);%将X_mean扩增[1 M]倍(repmat 复制平铺),减去均值
DataTest = reshape(X', num_rows, num_cols, N);%转置,排列成三维64*64*169
DataTest_ori = DataTest(:, :, :);
win_out = 11; % 外窗口
win_in = 3; % 内窗口
%% Local RX算法开始
Data = DataTest_ori;
[a b c] = size(Data);%返回数据几行几列
result = zeros(a, b);%创建一个a行,b列的零矩阵
t = fix(win_out/2);%朝零方向取整,Z_OUT
t1 = fix(win_in/2);%Z_IN
M1 = 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, :) = Data(:, b:-1:1, :);%原始数据左右翻转,填入西
DataTest(a+1:2*a, 2*b+1:3*b, :) = Data(:, b:-1:1, :);%原始数据左右翻转,填入东
DataTest(1:a, :, :) = DataTest(2*a:-1:(a+1), :, :);%中间的行镜像填入上面的行
DataTest(2*a+1:3*a, :, :) = DataTest(2*a:-1:(a+1), :, :);%中间的行镜像填入下面的行
for i = 1+b: 2*b
for j = 1+a: 2*a
block = DataTest(j-t: j+t, i-t: i+t, :);%规定外窗口大小(11*11)*3
y = squeeze(DataTest(j, i, :)).';%将1*1*3的三维矩阵转成1*3的二维矩阵。y依次等于每个像素点的二维矩阵(squeeze函数)
block(t-t1+1:t+t1+1, t-t1+1:t+t1+1, :) = NaN;%将内窗口处(3*3)*3区域内赋为非数值元素NaN
block = reshape(block, M1, c);%将三维转为二维*3,每一列为一个波段
block(isnan(block(:, 1)), :) = [];%将空矩阵赋给nan值(去掉每一列的nan值,现在每一列为121-9=112个数)(isnan判断是否为非数值元素,如果是返回1)
H = block.'; % 转置,每一行为一个波段
Sigma = (H * H');%同global的计算,自相关矩阵
Sigma_pinv = pinv(Sigma); %同global的计算,计算逆矩阵
% Sigma_inv = inv(Sigma);
result(j-a, i-b) = y * Sigma_pinv * y';%Rx算子
% result(j-a, i-b) = y * Sigma_inv * y';
end
end
%% Local RX算法结束,得出RX算子result
r2 = result;
r2 = reshape(r2, 1, M);%将每个像素点的Rx结果值排成行矩阵
tmp = reshape(r2, 256, 200);%重组成64*64
colormap('gray'),imagesc(tmp),colorbar; axis image;%根据值的大小颜色不同画图
title('双窗口 RX')
没有合适的资源?快使用搜索试试~ 我知道了~
Matlab.zip_dim small_expectwvx_individual57v_rxddetector_异常检测
共616个文件
jpg:500个
bmp:110个
m:3个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
5星 · 超过95%的资源 1 下载量 159 浏览量
2022-09-22
14:33:41
上传
评论
收藏 12.07MB ZIP 举报
温馨提示
针对弱小目标信号提出的RXD异常检测算法,实现对目标的检测
资源推荐
资源详情
资源评论
收起资源包目录
Matlab.zip_dim small_expectwvx_individual57v_rxddetector_异常检测 (616个子文件)
RXDDoubleWins.asv 2KB
RXDGLOBAL.asv 865B
43.bmp 51KB
13.bmp 51KB
19.bmp 51KB
84.bmp 51KB
48.bmp 51KB
37.bmp 51KB
36.bmp 51KB
28.bmp 51KB
56.bmp 51KB
29.bmp 51KB
71.bmp 51KB
34.bmp 51KB
1.bmp 51KB
110.bmp 51KB
3.bmp 51KB
98.bmp 51KB
60.bmp 51KB
47.bmp 51KB
108.bmp 51KB
66.bmp 51KB
20.bmp 51KB
44.bmp 51KB
83.bmp 51KB
100.bmp 51KB
90.bmp 51KB
58.bmp 51KB
93.bmp 51KB
14.bmp 51KB
68.bmp 51KB
59.bmp 51KB
75.bmp 51KB
65.bmp 51KB
80.bmp 51KB
72.bmp 51KB
55.bmp 51KB
40.bmp 51KB
106.bmp 51KB
81.bmp 51KB
53.bmp 51KB
54.bmp 51KB
88.bmp 51KB
6.bmp 51KB
11.bmp 51KB
89.bmp 51KB
7.bmp 51KB
67.bmp 51KB
25.bmp 51KB
109.bmp 51KB
62.bmp 51KB
27.bmp 51KB
92.bmp 51KB
77.bmp 51KB
101.bmp 51KB
70.bmp 51KB
32.bmp 51KB
12.bmp 51KB
105.bmp 51KB
31.bmp 51KB
86.bmp 51KB
74.bmp 51KB
4.bmp 51KB
78.bmp 51KB
17.bmp 51KB
82.bmp 51KB
30.bmp 51KB
33.bmp 51KB
8.bmp 51KB
57.bmp 51KB
26.bmp 51KB
85.bmp 51KB
103.bmp 51KB
39.bmp 51KB
63.bmp 51KB
22.bmp 51KB
95.bmp 51KB
16.bmp 51KB
51.bmp 51KB
23.bmp 51KB
79.bmp 51KB
52.bmp 51KB
64.bmp 51KB
9.bmp 51KB
107.bmp 51KB
76.bmp 51KB
91.bmp 51KB
50.bmp 51KB
61.bmp 51KB
18.bmp 51KB
5.bmp 51KB
42.bmp 51KB
99.bmp 51KB
38.bmp 51KB
69.bmp 51KB
10.bmp 51KB
21.bmp 51KB
94.bmp 51KB
2.bmp 51KB
49.bmp 51KB
共 616 条
- 1
- 2
- 3
- 4
- 5
- 6
- 7
资源评论
- zhizhi2050302024-03-05非常有用的资源,可以直接使用,对我很有用,果断支持!
小波思基
- 粉丝: 74
- 资源: 1万+
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功