设为首页 加入收藏

TOP

R中的空间数据分析
2019-08-15 00:10:05 】 浏览:36
Tags:空间 数据分析


> library(sp)
> library(maptools)
> library(raster)
> library(rgeos)

> maxd3 = readAsciiGrid("max.d3.txt")    # 加载数据
> rasterlayer.maxd3 = raster(maxd3)      # 转化为Raster
> maxd1 = readAsciiGrid("max.d1.txt")
> rasterlayer.maxd1 = raster(maxd1)     
> rasterlayer.max.d3.d1 = rasterlayer.maxd3 - rasterlayer.maxd1   # 比较3天与1天的数据
> rc <- reclassify(rasterlayer.max.d3.d1, c(-Inf,0,1, 0,Inf,NA))  # 对比较结果进行分类,提取出 有问题的数据

> rasterlayer.maxRange = rasterlayer.maxd1>0   #提取范围
> maxRange = rasterToPolygons(rasterlayer.maxRange, dissolve=TRUE)

> cols = grey.colors(20, 0.1, 0.9, 2.2)
> image(rc, col=cols, useRaster=TRUE)    #显示有问题的数据
> plot(maxRange, add=TRUE)               #显示整个范围


R中的空间数据分析 https://www.cppentry.com/bencandy.php?fid=91&id=228916

】【打印繁体】【投稿】【收藏】 【推荐】【举报】【评论】 【关闭】 【返回顶部
上一篇RxJS 实现摩斯密码(Morse) 【内附.. 下一篇时间序列分析工具箱——tibbletime