Python代码 对数比变换（汇总）

       代码包括“中心对数比”、“等角对数比”、“对称对数比”，其他的对数比的代码，后边有时间再补充。

#不存在0值时处理的代码
"""import numpy as np
#定义中心对数比函数
def center_log_transform(data):
    mean = np.mean(data)
    transformed_data = np.log(data) - mean
    return transformed_data
#示例数据
data = np.array([1,2,3,4,5])
#对数据进行中心对数比变换
transformed_data = center_log_transform(data)
print(transformed_data)"""

#存在零值时代码的处理
import numpy as np

'''# 定义中心对数比函数
def center_log_transform(data): 
    data = data + 1  # 加1操作
    mean = np.mean(data)
    transformed_data = np.log(data) - mean
    return transformed_data

# 示例数据
data = np.array([1, 2, 3, 4, 5])

# 对数据进行中心对数比变换
transformed_data = center_log_transform(data)
print(transformed_data)'''

#以下是数据不存在0值的情况
#中心对数比
import numpy as np
import pandas as pd

# 定义中心对数比函数
def center_log_transform(data):
    data = data
    mean = np.mean(data)
    transformed_data = np.log(data) - mean
    return transformed_data

# 读取Excel文件
excel_file_path = r'your_excel_file.xlsx'  #将 your_excel_file.xlsx替换为你的Excel文件路径
df = pd.read_excel(excel_file_path)

# 选择需要进行中心对数比变换的列
selected_column = 'your_column_name'  # 将 'your_column_name' 替换为你要处理的列名
data = df[selected_column].to_numpy()

# 对数据进行中心对数比变换
transformed_data = center_log_transform(data)

# 将变换后的数据添加到 DataFrame
df['TransformedData'] = transformed_data

# 将 DataFrame 存储为 Excel 文件
output_excel_path = r'输出路径'  # 输出 Excel 文件路径
df.to_excel(output_excel_path, index=False)

print(f"Transformed data saved to {output_excel_path}")

'''
#等角对数比
import numpy as np
import pandas as pd

# 定义等角对数比函数
def ilr_transform(data):
    data = data
    geometric_mean = np.exp(np.mean(np.log(data), axis=0))
    transformed_data = np.log(data / geometric_mean)
    return transformed_data

# 读取Excel文件
excel_file_path = r'文件输入路径'
df = pd.read_excel(excel_file_path)

# 选择需要进行等角对数比变换的列
selected_columns = ['x','y','z']  # 替换为你要处理的列名列表
data = df[selected_columns].to_numpy()

# 对数据进行等角对数比变换
transformed_data = ilr_transform(data)

# 将变换后的数据添加到 DataFrame
for i, col in enumerate(selected_columns):
    df[f'TransformedData_{i+1}'] = transformed_data[:, i]

# 将 DataFrame 存储为 Excel 文件
output_excel_path = r'输出路径'
df.to_excel(output_excel_path, index=False)

print(f"Transformed data saved to {output_excel_path}")'''


'''#对称对数比
import numpy as np
import pandas as pd

# 定义对称对数比函数
def slr_transform(data):
    data = data
    geometric_mean = np.exp(np.mean(np.log(data), axis=0))
    transformed_data = np.log(data / geometric_mean)
    return transformed_data

# 读取Excel文件
excel_file_path = r'输入路径'
df = pd.read_excel(excel_file_path)

# 选择需要进行对称对数比变换的列
selected_columns = ['x','y','z']  # 替换为你要处理的列名列表
data = df[selected_columns].to_numpy()

# 对数据进行对称对数比变换
transformed_data = slr_transform(data)

# 将变换后的数据添加到 DataFrame
for i, col in enumerate(selected_columns):
    df[f'TransformedData_{i+1}'] = transformed_data[:, i]

# 将 DataFrame 存储为 Excel 文件
output_excel_path = r'输出路径'
df.to_excel(output_excel_path, index=False)

print(f"Transformed data saved to {output_excel_path}")'''