我用ENVI对图像做了波段处理后合成了一张图像，怎么把其中每个大于0.2的数值求和？

bigliyang12 · 发表于 2018-12-23 20:35

我用ENVI对图像做了波段处理后合成了一张图像（可以看作算了一个NDVI)，请问现在想把数值大于0.2的所有值求和，应该怎么操作呢。然后有没有办法把这些值所占的区域的面积求出来呢？

jimi21 · 发表于 2018-12-24 09:30

提示: 作者被禁止或删除内容自动屏蔽

lk137095576 · 发表于 2018-12-24 10:12

用matlab
[row,column] = size(image);
sum = 0;
count = 0;
for m = 1:column
for n = 1:row
if image(m,n) > 0.2
sum = image(m,n) + sum;
count = count +1;
end
end
end

bigliyang12 · 发表于 2018-12-24 21:37

lk137095576 发表于 2018-12-24 10:12
用matlab
[row,column] = size(image);
sum = 0;

你好，matlab我也不大懂，要怎么操作呢？比如怎么用matlab读取图像信息，才能用这个代码实现

lk137095576 · 发表于 2018-12-25 09:43

bigliyang12 发表于 2018-12-24 21:37
你好，matlab我也不大懂，要怎么操作呢？比如怎么用matlab读取图像信息，才能用这个代码实现

用envi也能做，用快速统计工具，得到表格，然后自己根据需求算一下，如果matlab不太懂，推荐这种方法。
如果用matlab，就酱：
首先把图导出tif,虽然envi导出的tif看着是三个波段，实际上三个波段都一样。比如导出来文件叫'a.tif',当前文件夹调到文件所在文件夹,你把下面这段复制到命令行里面。得到结果里面，sum就是大于0.2的数值的和，count是大于0.2的数值的像元个数，面积用它乘以分辨率的平方。

a = imread('a.tif');
b = a(:,:,1);
[row,column] = size(b);
count = 0;
sum = 0;
for c = 1:column
for r = 1:row
      if b(r,c) > 0.2
         count = count + 1;
         sum = b(r,c) + sum;
      end
end
end