美国国家航空航天局开发了探测森林扰动的早期预警系统

为了开发一个近实时系统来检测和监测48个相邻州的植被变化，NASA的工程师需要分析Aqua和Terra卫星捕获的多光谱图像，建立收集数据的多时间尺度基线，然后将传入的数据与基线进行比较，以识别异常。

每天，MODIS数据生成14个网格瓷砖，跨越48个州，每个瓷砖包含4800 x 4800像素。在更新用于识别植被变化的年度基线时，该系统每年两次处理15至25 tb的数据。

NASA的工程师需要处理和分析存储在多维数组中的大量时间序列图像数据集。他们需要一个开发环境来交互开发和快速评估新的分析算法。在演示了他们的算法之后，他们想在不重新编码的情况下将其部署到生产环境中。

NASA工程师使用MATLAB开发了ForWarn的两个关键组件:时间序列产品工具(TSPT)，它对MODIS数据进行时间处理，以及物候参数估计工具(PPET)，它使用处理过的数据来计算绿色程度和生长季节的高峰期以及美国林务局感兴趣的其他物候参数。

对于TSPT, NASA的团队编写了一个MATLAB脚本，用于检索NASA存档中存储在分级数据格式(HDF)文件中的MODIS数据，并将数据导入MATLAB。TSPT调用Mapping Toolbox™函数将导入的经纬度数据转换为投影的地图坐标系。

TSPT波段处理模块也是在MATLAB中开发的，它从时间序列数据中生成归一化植被指数(NDVI)，以及土壤、水分、水分等指标。

在MATLAB中，该团队开发了TSPT算法，以消除因云、阴影、视图顶角和其他影响而扭曲的像素。

在将来自Aqua和Terra卫星的数据合并为单个时间序列后，TSPT使用信号处理工具箱™功能来识别和删除峰值和其他外围数据点。

一旦TSPT算法删除了异常值，它们将使用优化工具箱™和图像处理工具箱™进一步过滤和重新采样时间序列。TSPT算法应用信号处理工具箱中的Savitzky-Golay滤波器来插值任何缺失像素的值。

工程师们使用MATLAB、优化工具箱和图像处理工具箱开发了PPET，该工具对时间序列数据进行曲线拟合，以识别与每年周期性生长季节相关的植物状态，如发绿、成熟、衰老和休眠。后来，他们通过识别每日卫星数据和时间序列基线之间的差异，增强了PPET，以探测森林扰动。

使用MATLAB Compiler™，该团队创建了基于MATLAB算法的独立可执行版本，可以由未安装MATLAB的用户运行。

ForWarn团队因他们的努力获得了多个奖项，包括跨部门合作伙伴奖，该奖项表彰来自至少两个不同机构的联邦雇员，他们“在技术转移方面合作完成了出色的工作”。ForWarn团队由美国林务局、美国国家航空航天局、美国能源部橡树岭国家实验室、美国地质调查局和北卡罗来纳大学阿什维尔分校组成。

• 新想法在数小时内实现和测试。洛克希德·马丁公司的计算机工程师杰瑞·加瑟说:“信号和图像处理算法内置在MATLAB中，所以我不必从头开始编写它们。”“我没有编写和调试几百行C代码来测试一个新的算法思想，而是使用MATLAB命令交互式地验证它，然后编写一个脚本来测试多个参数值。我在几分钟或几小时内就能得到结果，而不是几天或几周。”
• 节省了数年的开发时间。“MATLAB处理多维数据的方式非常适合分析卫星图像，”阿姆斯特朗说。“如果我们使用传统的编程语言，比如C语言，我们会花更多的时间。”
• 向不断增长的用户社区交付可重用的生产软件。世界杯预选赛小组名单“ForWarn是美国大陆第一个近实时森林威胁早期预警系统，”阿姆斯特朗指出。它不断增长的用户群包括7000多名林业专家、学生和土地资源管理人员。我们也在重用MATLAB中开发的算法和模块，用于其他客户和目的，与世界各地的合作伙伴合作。”