为计算方便，按照图3所示，为无人机拍摄的每张标定版影像建立一个行列数相同的掩模影像。计算掩模影像中灰度为255的区域在标定版影像上对应区域的灰度平均值，按照式（1）Fi=pi/Σi∈P Li/n所示计算该波段的辐射定标系数。式中，Ri为波段i的标准反射率;P是由灰板掩模中灰度值为255的像素构成的集合，n为集合中像素的数目;Li为标定版影像经辐射校正后的波谱辐射亮度。MATLAB

代码如下：

% 读取标定版影像和掩模影像
calib_img = imread('calib_img.tif');
mask_img = imread('mask_img.tif');

% 计算掩模影像中灰度为255的区域在标定版影像上对应区域的灰度平均值
gray_vals = calib_img(mask_img == 255);
gray_avg = mean(gray_vals);

% 计算辐射定标系数
Ri = 0.5; % 假设波段i的标准反射率为0.5
n = numel(gray_vals);
Li = gray_avg;
Fi = Ri * Li / (n * Ri);

需要注意的是，这里的波段i的标准反射率Ri是根据实际情况进行设定的，需要根据具体数据进行调整。