RGB.rar_multivariate_rgb _多元线性回归 MATLAB_线性回归

clc,clear all; figure(1); file1='J:\ke_yan\P-grace-lag\resource_data_180\result3\part_dem.tif'; file2='J:\ke_yan\P-grace-lag\resource_data_180\result3\part_gpcc.tif'; file3='J:\ke_yan\P-grace-lag\resource_data_180\result3\part_ndvi.tif'; % [data1, R] = geotiffread(file1); [data2, ~] = geotiffread(file2); [data3, ~] = geotiffread(file3); [data, R1] = geotiffread('J:\ke_yan\P-grace-lag\resource_data_180\result3\part_tm.tif'); % % id=find(data==255);%研究区外界 % %data(isnan(data))=nan; % data1(id)=nan;%外界空值 % data2(id)=nan; % data3(id)=nan; % data1(data11==min(min(data11)))=0;%内部空值为0 % data2(data2==min(min(data2)))=0; % data3(data3==min(min(data3)))=0; % data1=round((data1-min(min(data1)))*255/(max(max(data1))-min(min(data1))));%归一化 round:四舍五入 data2=round((data2-min(min(data2)))*255/(max(max(data2))-min(min(data2)))); data3=round((data3-min(min(data3)))*255/(max(max(data3))-min(min(data3)))); data2=double(data2); [fA, fB, fC] = fractions(data1,data2,data3); fA(isnan(fA))=1; fB(isnan(fB))=1; fC(isnan(fC))=1; % fA(id)=1;%背景颜色 % fB(id)=1; % fC(id)=1; rgb=cat(3,fA, fB, fC); %按照dim连接fA,fB,fC 三个数组 转为3层 %%%%%%%%%%%%% imshow(data) hold on h=imshow(rgb); title('2000','fontsize',14,'fontweight','bold'); %fontsize：定义坐标轴标签和标题的字体大小：整数；默认值为12 %FontWeight：选择粗体或正常字体：{normal}、bold、light、demi %Title：定义坐标轴的标题 %%%%%%%%%%%%%%图中图 h2=axes('Position',[0.6 0.2 0.2 0.15]);%按照等边三角形比例设置 %Position：图形窗口的位置和大小；四维位置向量，格式为[坐标系最左下角定点的横坐标 坐标系最左下角定点的纵坐标 坐标系的总宽度 坐标系的总高度] fac = [1 2 3]; vert = [0 0;1 1*2^0.5;2 0];%边长为1 fvc = [1 0 0;0 1 0;0 0 1];%边长为1 patch('Faces',fac,'Vertices',vert,'FaceVertexCData',fvc,'FaceColor','interp') axis off %去掉坐标轴 text(0, -0.1,'dem', 'color','k', 'horizontalalignment', 'center','FontSize',14,'FontWeight','Bold');%坐标轴标签 text(2, -0.1,'gpcc', 'color','k', 'horizontalalignment', 'center','FontSize',14,'FontWeight','Bold'); text(1,1.58,'ndvi', 'color','k','horizontalalignment', 'center','FontSize',14,'FontWeight','Bold'); %%%%%%%%%%%%图中图结束 print(1,'-r500','-djpeg',strcat('use_fraction2005','.jpg'));