使用Astropy将数组转为无量纲遇阻，调用np.log报错求助

你的问题出在lum_in_W*u.dimensionless_unscaled这一步：该操作只是给原带单位的量乘了一个无量纲缩放因子，并不会移除原有单位（本质等价于(数值 × W/Hz) × 1 = 数值 × W/Hz），因此最终的量依然带有W/Hz单位，触发np.log()的单位错误。

以下是两种可行的解决方法：

方法一：提取纯数值直接运算 #

如果不需要保留Astropy Quantity类型，直接提取数值即可完成对数计算：

log_lum = np.log(lum_in_W.value)

方法二：转换为无量纲Quantity #

如果需要保留Quantity类型（仅为无量纲），可以通过两种方式实现：

1. 用自身单位做除法，得到无量纲量：

lum_dimensionless = lum_in_W / lum_in_W.unit
# 验证单位，输出应为 dimensionless_unscaled
print(lum_dimensionless.unit)
log_lum = np.log(lum_dimensionless)

2. 基于原数值创建无量纲Quantity：

lum_dimensionless = u.Quantity(lum_in_W.value, u.dimensionless_unscaled)
log_lum = np.log(lum_dimensionless)

另外，若你需要计算的是相对参考值的对数光度（比如以某特征光度为基准），可直接计算比值（自动转为无量纲）：

# 示例：定义参考光度
ref_luminosity = u.Quantity(1e20, u.Watt/u.Hertz)
# 计算比值，自动成为无量纲量
lum_ratio = lum_in_W / ref_luminosity
log_lum = np.log(lum_ratio)

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