从图像中选择像素位置并计算反射率，最后将结果保存为反射率影像，并能框选范围计算范围内反射率平均值。MATLAB

首先，可以使用MATLAB中的`imtool`函数打开图像，然后使用鼠标在图像中选择像素位置，获取像素的位置信息。然后使用上一个问题中的方法计算该位置的反射率，并使用`imwrite`函数保存为反射率影像。

对于框选范围计算范围内反射率平均值，可以使用`imrect`函数在图像上绘制矩形框，然后计算该框选范围内的所有像素的反射率平均值。具体步骤如下：

1. 打开图像并选择像素位置：

im = imread('image.tif');
imshow(im);
[x, y] = ginput(1); % 获取一个像素位置

2. 计算该位置的反射率：

% 设置反射率的计算参数，例如太阳入射角、大气校正等
% 这里假设已经得到了这些参数
solar_angle = 30;
atmospheric_correction = 0.5;

% 计算反射率
reflectance = (im(y, x) * pi) / (solar_angle * atmospheric_correction);

3. 保存反射率影像：

imwrite(reflectance, 'reflectance.tif');

4. 绘制矩形框并计算反射率平均值：

rect = imrect(); % 绘制矩形框并返回其句柄
mask = rect.createMask(); % 创建矩形框范围的掩膜
reflectance_roi = reflectance .* mask; % 在反射率图像上应用掩膜
reflectance_mean = mean(reflectance_roi(:)); % 计算反射率平均值

5. 显示结果：

imshow(im); % 显示原始图像
rect.draw(); % 绘制矩形框
text(x, y, num2str(reflectance), 'Color', 'r'); % 在选择的像素位置处显示反射率
text(rect.X + rect.Width/2, rect.Y + rect.Height/2, num2str(reflectance_mean), 'Color', 'r'); % 在矩形框中心位置处显示反射率平均值

这样就可以选择图像中的像素位置，并计算反射率。同时，还可以框选范围计算范围内反射率平均值，并在图像上显示结果。