【问题标题】:Get an accurate Y value from a raster on R从 R 上的栅格中获取准确的 Y 值
【发布时间】:2018-12-14 15:48:50
【问题描述】:

我正在尝试获取栅格在空间中的确切位置,以便精确确定与栅格元素重叠的线参数。但是,在我使用的存储库中,ymin 和 ymax 坐标总是四舍五入(总是显示为整数)。

> ymin(vegetation_raster_1)
[1] 4894166

我们可以通过图形看到,ymin 不完全是 4894166,而是接近 4894165.7。

> vegetation_raster_1
class       : RasterLayer 
dimensions  : 94, 95, 8930  (nrow, ncol, ncell)
resolution  : 0.05744, 0.05744  (x, y)
extent      : 284823.7, 284829.2, 4894166, 4894171  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=31 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : in memory
names       : band1 
values      : 0, 0.3753507  (min, max)

我尝试过使用 extent() 函数,但结果相同。

> extent(vegetation_raster_1)
class       : Extent 
xmin        : 284823.7 
xmax        : 284829.2 
ymin        : 4894166 
ymax        : 4894171 

有谁知道小数点后至少有一位数字的 ymin 和 ymax 的方法?

问候,

弗朗索瓦

【问题讨论】:

  • 您能否提供一些可以复制和粘贴问题的数据,以便我们自己重现问题?
  • 抱歉耽搁了。在我发送示例之前,Robert 已经提供了该问题的解决方案...

标签: r coordinates r-raster


【解决方案1】:

这些值没有四舍五入。打印时可能不会显示所有小数。你可以用不同的方式控制它。例如

x <- 12.12300000100000123   
x
#[1] 12.123
sprintf("%.10f", x)
#[1] "12.1230000010"

【讨论】:

    猜你喜欢
    • 2019-11-08
    • 2019-06-27
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2021-11-30
    • 1970-01-01
    相关资源
    最近更新 更多