0

我需要在 15 分钟内比较 Meteosat-8 的 2、3、6 和 9 波段。我已经计算了亮度温度,并设法将所有 Tif 读入 IDL。

现在我做了一些比率。像 B6-B4、B9-3... 现在我想创建一个新的 Tif,当达到某些条件时,它会在其中写入(然后加起来)1 或 0。

结果应该是一个 Tif,其中我有每个像素中“1”的数量。我想我必须用 cols/rows 创建一个新数组,并在每个“for循环”之后写入其中的 1 和 0 并每次将其相加。

谢谢!布米

4

1 回答 1

0

我设法让它读取 tiff 并在新的 tif 中比较它们。我做了一个 for 循环来加载所有路径。

该程序一直运行,直到缺少文件(路径不存在)。所以我想跳过(继续)所有丢失的“路径”。

我想,我可以用 file_info 或 file_test 来做。

x=file_info(filepath(path), exists) 如果 x eq 0 然后继续

存在:如果文件存在,则为 True(1)。如果不存在,则为假 (0)。

谢谢

于 2018-06-15T10:26:36.187 回答