使用R更新下垫面资料为FROM-GLC 30m数据

aslamic

最早发于知乎，因为当时论坛崩了。
受该帖子的启发，早在2017年底就实现了这个方法，当时用的是中科院1km的土地利用类型数据，但是接壤的地方没有数据是个大问题，去年更新了FROM-GLC 2015，但但是这套数据虽然分辨率高，有些地方分类又是错的，总之不管怎样，共享给大家供讨论。
FROM-GLC 2017分类倒是没啥大问题了，但是只分了10类，比较局限。程序只适用于modis数据。

堆码水平有限，见笑了。
建议1：如果没有办法订正FROM-GLC的分类问题，只更新城市和水体；
建议2：在Windows下面运行。

详情在此：
https://zhuanlan.zhihu.com/p/118510946

核心的思路是：WRF静态数据（Static Data）中很重要的就是下垫面资料，包括土地利用类型、土壤类型、叶面积指数等，模型自带了这些数据，但就土地利用类型而言，我国城市化速度快，下垫面变化较为明显，模型自带的下垫面资料实效性较差，USGS为1992年，MODIS则为2002年。AMLUpdater则是通过读取现有土地利用类型数据，建立分类与模型分类的关系，利用网格经纬度范围获取现有下垫面数据，计算各网格内不同土地利用类型之间的关系，最后写入geo_em文件。
更新前：

                               
登录/注册后可看大图

更新后：

                               
登录/注册后可看大图

过几天空了上传到github，多路更新代码太折腾了，这里就不贴代码了，代码量不多，但是实现比较暴力。

在这个代码的基础上，可以通过修改实现土壤类型等数据的更新。
至于为什么不直接生成geog二进制文件，这块目前还有点技术难点没有解决，而且我暂时没法脱离gis纯靠R代码实现，等后续成熟了再共享给大家。

shubabu

好帖子 学习了学习了

传说小小星

{:5_235:}{:5_235:}{:5_235:}{:5_235:}{:5_235:}牛啊