代码之家 › 专栏 › 技术社区 › corvus

R: 从density()获取数字输出

scientific-notation raster gis r

corvus · 技术社区 · 8 年前

spatialgrid <- as(density(mypattern,5000,eps=50),'SpatialGridDataFrame')
rastergrid <- raster(spatialgrid)
writeRaster(rastergrid, filename=â/file.tifâ,format=âGTiffâ)

然而,当我将生成的光栅加载到QGIS中时,由于单元格值是用科学符号而不是简单的数字书写的,所以我有一些问题。

基于 this question ,我试过了 format(density(),scientific=FALSE) 但这导致CPU出现严重峰值,运行时间很长,因此我最终终止了该进程。

我想找到一种方法,让density()函数输出整数值。或者,也许有一种方法可以将数据帧转换为整数数据类型?

2 回复 | 直到 8 年前

jgadoury 8 年前

你可以使用 options(scipen=99) 按照@dash2的建议,删除科学符号。
将你的值乘以一个因子,例如1000。尝试一些方法,直到你的值不再是科学符号。优点是光栅的前导零不会像选项1那样多,并且占用的硬盘空间更少,但缺点是您必须将QGIS中的值乘以该因子。

Robert Hijmans 8 年前

“科学记数法”与此无关。这是用于文本表示,但使用writeHost,您可以写入数字格式。它们只是真实的数字。

我认为您可以简化:

d <- density(mypattern,5000,eps=50)
r <- raster(d)

问题可能是密度很低。如果您想要更大的数字(正如jgadoury建议的那样)

r <- r * 1000000

可能为圆形

r <- round(r)

然后

writeRaster(r, 'file.tif')

推荐文章

dwiz · R中具有焦函数和parlappy的光栅并行处理列表

7 年前

krishnab · 谷歌地球引擎:寻找正确的图像尺寸以获得卫星图像的最大分辨率

7 年前

Hugo · 具有随机森林的图像分类(光栅堆栈)(package ranger)

7 年前

B. Davis · 创建包含特定光栅值的多边形

7 年前

MapUser · 如何从光栅中提取多边形?

8 年前

corvus · R: 从density()获取数字输出

8 年前

Joke O. · 使用R中的savitzky golay/sgolayflt(信号)平滑rasterstack

8 年前

N.Varela · R: 如何通过在基本R图中用光栅覆盖地图来获得正确的分辨率(光栅)?

9 年前

Geo-sp · 在水平图面板区域外添加文本

9 年前

user32882 · 将光栅文件作为阵列处理?

10 年前