R语言高分辨率遥感卫星数据处理

目的：由于自己的需求，需要对原始30m分辨率的土地利用数据的分类ID进行重写。该工作并不复杂，但由于数据分辨率太高且需要处理全国尺度，处理时遇到了很多问题，其中最主要的问题是数据读写溢出。对于此类问题，自己主要通过将原始数据进行瓦片分割、处理和再合并。期间踩过了很多坑，这些坑主要标注在下面的代码中，与大家分享。

数据来源：http://data.ess.tsinghua.edu.cn/fromglc2017v1.html

处理流程：批下载、瓦片分割、分类ID重新写入写出和数据合并

处理代码：

###第一步:网格数据转投影和重采样成30m分辨率(img Byte)
###第二步:根据第一步的数据创建合适的渔网
###第三步:根据渔网裁剪栅格并修改土地利用类型编号
###第四步:合并数据
###第五步:输出domain