是否有更高效的方法读取DEC格式二进制文件的4字节浮点数？

问题背景 #

我有一批采用DEC格式存储的旧二进制文件，目前读取其中4字节浮点数的步骤为：

• 读取字节；
• 交换前后两个字（即交换第1、2字节与第3、4字节）；
• 使用readBin()函数将字节转换为数值；
• 将结果除以4。

原本以为readBin()的endian参数（可选值：'little'、'big'、'swap'）能处理该格式，但实际无法满足需求。当前实现的示例代码如下：

# 示例目标值为1.290，对应4字节DEC格式数据
# 实际场景中dec_bytes通过readBin(con, raw(), n=4)读取
dec_bytes <- writeBin(1.290, raw(), size=4)
# 交换前后两个字的字节顺序
pc_bytes <- c(dec_bytes[3], dec_bytes[4], dec_bytes[1], dec_bytes[2])
# 转换为数值
pc_float <- readBin(pc_bytes, numeric(), n=1, size=4)
pc_float 
# [1] 0.5161456
# 除以4得到正确结果
pc_float <- pc_float / 4
pc_float 
# [1] 0.1290364

我可以将这些步骤封装为函数，但希望找到更简便高效的实现方式。另外，我有一段30年前的C代码，通过调整字节的方式无需除以4即可得到正确值，核心逻辑是重排字节后对最后一个字节减1（调整指数）：

float ConvertDecToFloat(char bytes[4])
{
    char p[4];
    p[0] = bytes[2];
    p[1] = bytes[3];
    p[2] = bytes[0];
    p[3] = bytes[1];
    if (p[0] || p[1] || p[2] || p[3])
        --p[3];          // 调整指数

    return *(float*)p;
}

想知道在R中能否不通过整数转换再转回字节的方式实现该逻辑。

解决方案 #

1. 直接实现C代码的字节调整逻辑（R版本） #

在R中可以直接操作raw类型的字节，无需绕路整数转换，直接复刻C代码的逻辑：

convert_dec_float <- function(dec_bytes) {
  # 交换前后两个字的字节顺序
  pc_bytes <- c(dec_bytes[3:4], dec_bytes[1:2])
  # 非零值时调整最后一个字节（修正指数偏移）
  if (any(pc_bytes != as.raw(0))) {
    pc_bytes[4] <- pc_bytes[4] - as.raw(1)
  }
  # 转换为标准浮点数
  readBin(pc_bytes, numeric(), n=1, size=4)
}

# 测试示例
dec_bytes <- writeBin(1.290, raw(), size=4)
convert_dec_float(dec_bytes)
# [1] 0.1290364

该方法直接对应DEC格式到IEEE单精度浮点数的底层转换逻辑：DEC浮点数的指数偏移量为128，而IEEE单精度为127，对最后一个字节（指数位所在字节）减1，等价于将指数偏移修正为127，对应数值缩小为原来的1/4（即2^-2），省去了额外的浮点除法步骤。

2. 效率对比 #

两种方法的性能差异极小，但直接调整字节的方式在语义上更贴合格式转换原理，极端大数据量场景下会略快。可以用microbenchmark包验证：

library(microbenchmark)

dec_bytes <- writeBin(1.290, raw(), size=4)
original_method <- function() {
  pc_bytes <- c(dec_bytes[3:4], dec_bytes[1:2])
  readBin(pc_bytes, numeric(), n=1, size=4) / 4
}

# 运行10000次对比
microbenchmark(
  original_method(),
  convert_dec_float(dec_bytes),
  times = 10000
)

3. 关于endian参数的说明 #

readBin()的endian='swap'仅负责交换字节顺序，但DEC格式与IEEE单精度浮点数的差异不仅在于字节顺序，还包括指数偏移量（DEC为128，IEEE为127），因此单纯调整字节顺序无法完成正确转换，必须额外修正指数。

内容的提问来源于stack exchange，提问作者bmacwilliams