clc;
clear;
close all;
warning off;
addpath(genpath(pwd));
size_x = 200;
size_y = 200;
size_lambda = 100;
y = zeros(size_x,size_y,size_lambda);
D = zeros(size_x*size_y, size_lambda);
load('Data/celegans_spec_init.mat');
for i = 0:(size_lambda - 1)
if i<10
temp = ['Stack000',num2str(i),'.tif'];
else
temp = ['Stack00',num2str(i),'.tif'];
end
img = imread(temp);
temp_column = reshape(img,size_x*size_y,1);
y(:,:,i+1) = img;
D(:,i+1) = temp_column;
end
clear size_x size_y size_lambda temp temp_column i img;
rmpath('Data/CElegans');
lambda = 0.6e1;
copt = C_ML_Lasso( y, celegans_spec_init',lambda);
sopt = celegans_spec_init;
Comp1 = copt(:,:,1);
Comp2 = copt(:,:,2);
Comp3 = copt(:,:,3);
Comp4 = copt(:,:,4);
Comp1 = imadjust(uint8(255.*(Comp1 - min(Comp1(:) ))/(max(Comp1(:))-min(Comp1(:) ))));
Comp2 = imadjust(uint8(255.*(Comp2 - min(Comp2(:) ))/(max(Comp2(:))-min(Comp2(:) ))));
Comp3 = imadjust(uint8(255.*(Comp3 - min(Comp3(:) ))/(max(Comp3(:))-min(Comp3(:) ))));
Comp4 = imadjust(uint8(255.*(Comp4 - min(Comp4(:) ))/(max(Comp4(:))-min(Comp4(:) ))));
rgbComp1 = cat(3, Comp1*0.5, Comp1*0, Comp1*0.5); % magneta for oxidized lipid
rgbComp2 = cat(3, Comp2*1, Comp2*0, Comp2*0); % red for nucleus
rgbComp3 = cat(3, Comp3*0, Comp3*0.8, Comp3*0); % green for protein
rgbComp4 = cat(3, Comp4*0, Comp4*0.3, Comp4*0.3);% cyan for water
C_disp = rgbComp1+rgbComp2+rgbComp3+rgbComp4;
figure;
imshow(C_disp,'InitialMagnification',150);
Spec1 = sopt(1,:);
Spec2 = sopt(2,:);
Spec3 = sopt(3,:);
Spec4 = sopt(4,:);
figure;
plot(Spec1,'m','Linewidth',1.5);
hold on;
plot(Spec2,'r','Linewidth',1.5);
plot(Spec3,'g','Linewidth',1.5);
plot(Spec4,'c','Linewidth',1.5);
xlabel 'Channel';
ylabel 'Int. (a.u.)';
legend('Oxidized lipid','Fat droplets','lysosome-related organelles','Protein')
pbaspect([2 1 1]);
没有合适的资源?快使用搜索试试~ 我知道了~
温馨提示
1.版本:matlab2021a,我录制了仿真操作录像,可以跟着操作出仿真结果 2.领域:LASSO算法 3.内容:基于LASSO算法的稀疏正则化高光谱图像的光谱解混合算法matlab仿真+仿真录像
资源推荐
资源详情
资源评论
收起资源包目录
基于LASSO算法的稀疏正则化高光谱图像的光谱解混合算法matlab仿真+仿真录像 (106个子文件)
操作录像0035.avi 48.8MB
312.jpg 20KB
12.jpg 18KB
Runme.m 2KB
C_ML_Lasso.m 1KB
celegans_spec_init.mat 2KB
Stack0079.tif 39KB
Stack0032.tif 39KB
Stack0045.tif 39KB
Stack0051.tif 39KB
Stack0064.tif 39KB
Stack0053.tif 39KB
Stack0034.tif 39KB
Stack0092.tif 39KB
Stack0085.tif 39KB
Stack0068.tif 39KB
Stack0016.tif 39KB
Stack0056.tif 39KB
Stack0060.tif 39KB
Stack0049.tif 39KB
Stack0038.tif 39KB
Stack0006.tif 39KB
Stack0043.tif 39KB
Stack0039.tif 39KB
Stack0058.tif 39KB
Stack0089.tif 39KB
Stack0095.tif 39KB
Stack0097.tif 39KB
Stack0002.tif 39KB
Stack0082.tif 39KB
Stack0047.tif 39KB
Stack0033.tif 39KB
Stack0063.tif 39KB
Stack0044.tif 39KB
Stack0030.tif 39KB
Stack0088.tif 39KB
Stack0022.tif 39KB
Stack0012.tif 39KB
Stack0005.tif 39KB
Stack0059.tif 39KB
Stack0019.tif 39KB
Stack0081.tif 39KB
Stack0009.tif 39KB
Stack0037.tif 39KB
Stack0028.tif 39KB
Stack0048.tif 39KB
Stack0000.tif 39KB
Stack0036.tif 39KB
Stack0074.tif 39KB
Stack0007.tif 39KB
Stack0041.tif 39KB
Stack0098.tif 39KB
Stack0026.tif 39KB
Stack0076.tif 39KB
Stack0014.tif 39KB
Stack0021.tif 39KB
Stack0094.tif 39KB
Stack0065.tif 39KB
Stack0084.tif 39KB
Stack0073.tif 39KB
Stack0010.tif 39KB
Stack0080.tif 39KB
Stack0062.tif 39KB
Stack0057.tif 39KB
Stack0067.tif 39KB
Stack0099.tif 39KB
Stack0018.tif 39KB
Stack0078.tif 39KB
Stack0050.tif 39KB
Stack0086.tif 39KB
Stack0046.tif 39KB
Stack0054.tif 39KB
Stack0015.tif 39KB
Stack0001.tif 39KB
Stack0066.tif 39KB
Stack0008.tif 39KB
Stack0020.tif 39KB
Stack0075.tif 39KB
Stack0025.tif 39KB
Stack0023.tif 39KB
Stack0091.tif 39KB
Stack0070.tif 39KB
Stack0011.tif 39KB
Stack0017.tif 39KB
Stack0042.tif 39KB
Stack0087.tif 39KB
Stack0096.tif 39KB
Stack0003.tif 39KB
Stack0027.tif 39KB
Stack0031.tif 39KB
Stack0077.tif 39KB
Stack0093.tif 39KB
Stack0090.tif 39KB
Stack0029.tif 39KB
Stack0055.tif 39KB
Stack0071.tif 39KB
Stack0083.tif 39KB
Stack0072.tif 39KB
Stack0052.tif 39KB
Stack0013.tif 39KB
共 106 条
- 1
- 2
资源评论
- weixin_458610472023-03-14这个资源内容超赞,对我来说很有价值,很实用,感谢大佬分享~
- 博德安2023-10-02发现一个宝藏资源,资源有很高的参考价值,赶紧学起来~
fpga和matlab
- 粉丝: 15w+
- 资源: 2548
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功