【栅格数据计算分析】分类后的栅格提取增加和减少区域

经过一些头铁尝试，现在觉得“遥感处理一般用栅格数据，少用矢量”，是挺有道理的

---

提取属性中特定字段

之前看的是用栅格计算器【Pick函数】。啥也不是。对提取特定字段根本不好用。
想了解的可以参考：【教程】ArcGIS中栅格计算器常用函数的使用
https://www.bilibili.com/video/BV1Dc41187Pd/?p=3&vd_source=849c6a56de214c389a4ae60afc
讲得还是比较清楚。

经过又一番搜索，使用下面的工具。
栅格数据的空间分析——提取分析https://www.jianshu.com/p/565d69186ec6
工具：按属性提取

下一步就是：重分类
原因：知乎：Arcgis栅格数据重分类——实例操作https://zhuanlan.zhihu.com/p/385435549

我只想要一类特定数据。最后的结果为：
上面是重分类的，下面是刚提取出来的

另，如果栅格有多帧，建议在工具处，右击，选择【批处理】

---

注意，上面的提取有一点点问题。就是提取的结果没有加tif后缀。虽然是栅格，但不知道是什么栅格格式，换一台电脑或者从map改道pro，如上图所示的Reclass_18就不稳定。故，请加上tif后缀。

或者，直接用【extract by attribute】，提取特定栅格。

栅格计算器

这一步是为了获得，两帧影像相同的地方。方便后面计算增加还是减少。
【其实是走了很多弯路暴躁*n之后，才发现的】

用栅格计算器。con(a.tif != b.tif,0,1)
我是这个公式。算出来的数据是两个栅格相同的地方。这里因为我数据有点，奇怪，所以这一步请有需要的同学，自己搜下。我的有，特殊性和巧合性。

栅格相减

首先，栅格计算器直接相减，pass。出来的结果还是两个重叠的地方。原因暂不知。
分析工具-减法，pass。理由同上。
交集取反【】SymDifff，只适合要素，pass
Erase，擦除，只适合要素，Pass

最后：
栅格擦除(反向裁剪)https://blog.csdn.net/u010885059/article/details/82348514
ArcMap中如何用栅格裁剪矢量（栅格擦除）https://www.likecs.com/show-203310216.html
我主要参考的是，第一个链接。步骤少一点

注意我划红线的地方，第一个是需要擦除的影像，第二个是包含想要留下部分的影像。
以及，加一个tif后缀，稳定一点。

下一步没什么，属性提取、因为你擦除后的影像会新增一个字段。就是你想要的结果，把这个字段提取出来。

下图红圈是原来的。

图中的红圈里面，就是我栅格计算后，由原来1个字段，新增了1个。新增的这个就是我想要的。
可以上下两个图对比一下。

最后，提醒下增加和减少的问题。
比如，18年和19年。
用18减去18_19重合的，就是18年减少的。
用19减去18_19重合的，就是19年增加的。

完毕。

---

其实还有GEE或者arcpy方法可以做。后期可以搞一下。