计算不同分类特征指标对9种土壤类型的J-M距离值python代码

这里提供一份示例代码，希望能帮到你：

import pandas as pd
from scipy.spatial.distance import pdist
from scipy.spatial.distance import squareform

# 读取数据
data = pd.read_csv('soil_data.csv')

# 提取分类特征和J-M距离值所需的数据
features = data.iloc[:, 1:11]   # 假设分类特征在第1到第10列
distances = data.iloc[:, 11:]   # 假设J-M距离值在第11列及以后

# 计算不同分类特征指标对J-M距离值的距离
dist_matrix = pdist(distances, metric='jaccard')   # 使用Jaccard距离
dist_matrix = squareform(dist_matrix)

for i in range(features.shape[1]):
    for j in range(i+1, features.shape[1]):
        feature_pair = [features.iloc[:, i], features.iloc[:, j]]
        feature_pair = pd.concat(feature_pair, axis=1)
        feature_pair = feature_pair.dropna()
        feature_dist_matrix = pdist(feature_pair, metric='jaccard')
        feature_dist_matrix = squareform(feature_dist_matrix)
        print(f"分类特征{i+1}和{j+1}的J-M距离值：")
        print(feature_dist_matrix)

这份代码假设数据文件`soil_data.csv`中第1到第10列是分类特征，第11列及以后是J-M距离值。

代码中使用的是Jaccard距离，你也可以根据需要选择其他距离度量方法。最后输出的是每对分类特征指标对应的J-M距离值，共有45对（C(10,2)）不同的组合。