粒度分析端元分析方法及代码分享

很多小伙伴问粒度常用的端元分析如何做？Matlab中Analysize程序包是最常用、最经典的方法。由于我电脑之前卸载了matlab安装包，实在不想安装16G的软件，所以一直没有机会给大家录视频。今天正好整理资料，看到之前的一篇学习记录，分享给有需要做端元分析的朋友。以下是简要的分析步骤，如果有需要Matlab和Analysize程序包的小伙伴，也麻烦一键三连后给我私信。

Matlab加载上AnalySize程序包：

Matlab窗口，输入命令AnalySize，运行程序(注意大小写)

加载数据

File，Load data，excel数据做成下面的格式

根据数据的格式选取，列和行序列。

然后执行⑧perform， 并选择无参还是有参的计算方法。有参中可选择SSG、webull函数等

根据R2和拐点的情况判断端元选取情况，选端元个数（一般2-4个）。可以在这里选择save plot，保存为eps文件，后期在coreldraw中矢量化调整。

Ok之后可以在图3中根据数据情况进行分析，选取tool中进行单个异常数据剔除和某范围内异常数据删除。之后，可以在①中选取对应样品看拟合情况，④代表整个序列中提取出的端元情况。然后在⑥中导出数据，All data代表原始粒度数据，不用到处。Export plot data,代表导出当前样品的数据。Export EM data代表导出端元的数据，包括贡献量（Abundance）和概率分布（density）。其中贡献量最重要，根据density可作出图④，当然也可以直接下载③和④为eps文件，后期矢量修改。 注：在导出Abundance的时候，样品号specimen这一列的数据可能导出为异常值，可以直接用复制粘贴来获取abundance的数据。

此外还可以通过Other Tools对端元参数进行分析，根据实际情况选取几何还是对数等分析方法：

同时，该程序包还提供了CM图像和三端元分析，也可以导出为EPS矢量图。

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删