k.m

来自「几个计算」· M 代码 · 共 26 行

% % A：1GAIN、2OFFSET、3中心波长、4植被反射率、5土壤反射率、6直射辐射通密、7散射辐射通密、8比值
A(3,3)=0.545;%green
A(3,2)=0.645;%red
A(3,1)=0.840;%ir
A(4,1)=0.108615;
A(4,2)=0.055911;
A(4,3)=0.493660;
% A(5,1)=0.154629%*2;
% A(5,2)=0.142932%*2;
% A(5,3)=0.189277%*2;
% A(5,1)=0.154629;
% A(5,2)=0.269561;
% A(5,3)=0.189277;
A(5,3)=0.357328;
A(5,2)=0.269667;
A(5,1)=0.158707;
A(6,1)=1858;
A(6,2)=1575;
A(6,3)=1047;

Es_0561 = 0.1*1858;
for i = 1:3
    A(7,i) = ((0.545/A(3,i))^4)*Es_0561;
    A(8,i) = A(7,i)/(A(6,i)*cos(44.448334*pi/180)+A(7,i));
end
save A A;