Python Pint库如何正确自定义温度梯度单位实现转换

• Pint默认将.解析为单位乘法连接符，因此deg.C会被拆分为deg（平面角单位“度”）和C（电荷量单位“库仑”）做乘法运算，根本不会匹配你自定义的别名规则，最终才会输出coulomb * degree / kilometer的错误解析结果。
• 原自定义配置文件还存在两个问题：一是温度梯度换算系数错误，二是没有正确区分绝对温度单位和温差单位，温度梯度计算必须使用温差单位degC/degF，而非带偏移量的绝对温度单位celsius/fahrenheit。

第一步：修改自定义单位定义文件 #

将GEOPHIRES3_newunits.txt的内容替换为如下配置：

# 温度梯度相关单位定义
# 给内置温差单位绑定带点格式的别名，别名优先级高于点号乘法解析规则
@alias degC = deg.C
@alias degF = deg.F
@alias mile = mi

# 可选：定义梯度量纲，用于后续量纲合法性校验
[gradient] = [delta_temperature] / [length]

# 梯度单位定义，使用精确换算系数
celsius_per_km = degC / km = degC/km = °C/km
fahrenheit_per_mile = celsius_per_km * (9/5) / 1.609344 = degF/mi = deg.F/mile = °F/mile
fahrenheit_per_km = celsius_per_km * 9/5 = degF/km = deg.F/km = °F/km
celsius_per_mile = celsius_per_km / 1.609344 = degC/mi = deg.C/mile = °C/mile

第二步：修正Python代码逻辑 #

注意修正原代码中的变量名拼写错误（Myquanity漏写字符t），参考如下可直接运行的代码：

from pint import UnitRegistry, Quantity
ureg = UnitRegistry()
ureg.load_definitions('GEOPHIRES3_newunits.txt')

my_quantity = Quantity("55 deg.C/km")
print(my_quantity)  # 输出：55 degree_Celsius / kilometer

# 转换为deg.F/mile
my_quantity_f_mile = my_quantity.to("deg.F/mile")
print(my_quantity_f_mile)  # 输出约为156.23 degree_Fahrenheit / mile

# 转换为deg.C/mile
my_quantity_c_mile = my_quantity.to("deg.C/mile")
print(my_quantity_c_mile)  # 输出约为88.51 degree_Celsius / mile

• 使用@alias指令绑定别名的优先级高于默认的.乘法解析规则，从根源解决deg.C被错误拆分识别的问题。
• 温度梯度计算必须使用温差单位degC/degF，这类单位代表1摄氏度/华氏度的温度差值，没有零点偏移，直接按比例换算即可；如果误用绝对温度单位celsius/fahrenheit，会因为温度单位自带的零点偏移量导致换算结果完全错误。
• 换算时使用精确换算比例：摄氏度和华氏度的温差换算比例为9:5，1英里精确等于1.609344公里，比手动写近似值的换算精度更高。

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