我需要比较Meteosat-8的波段2,3,6和9全部15分钟。我已经完成了亮度温度的计算,我设法将所有Tif读入IDL。
现在我做了一些比率。就像B6-B4,B9-3 ......现在我想创建一个新的Tif,当达到某些条件时,它会写入(后来加起来)1或0。
结果应该是Tif,其中我有" 1"在每个像素。 我想我必须创建一个带有cols / rows的新数组,并在每个" for循环后写入"它中的1和0并且每次都加起来。
谢谢! 布米
答案 0 :(得分:0)
我设法读取tiff并在新的tif中比较它们。 我做了一个for循环来加载所有路径。
程序一直有效,直到丢失文件(路径不存在)。所以我想跳过(继续)所有缺失的“路径”。
我想,我可以用file_info或file_test来做。
x = file_info(文件路径(路径),存在) 如果x eq 0则继续
存在:如果文件存在,则为True(1)。如果不存在则为假(0)。
感谢