python编译环境配置的实现步骤

# -*- coding: utf-8 -*-
"""
Created on Sun May 25 19:13:37 2025

@author: galax
"""

# 导入必要的库
import cv2
import numpy as np
from skimage.feature import graycomatrix, graycoprops
from matplotlib import pyplot as plt


# 读取PNG图像
#"C:\Users\galax\Desktop\数据Data\p\Vesta03.png"
#image_path = '/kaggle/input/asteroids-pngtest/Bennu02.png'
image_path = r'C:\Users\galax\Desktop\Data\p\Vesta03.png'
image = cv2.imread(image_path, cv2.IMREAD_UNCHANGED)

# 检查图像是否加载成功
if image is None:
    raise ValueError("无法加载图像！请检查文件路径")

# 显示原始图像信息
print("[调试信息] 原始图像形状:", image.shape)  # 应显示 (1024, 1024)
print("[调试信息] 图像数据类型:", image.dtype)  # 应显示 uint8

# 灰度图像专用处理流程 ---------------------------------------------------
# 直接使用原始灰度图像（无需转换）
if len(image.shape) == 2:
    gray = image  # 直接使用原始灰度数据
    print("[调试信息] 检测到灰度图像，跳过通道转换")
else:
    # 如果意外出现彩色图像则转换为灰度
    gray = cv2.cvtColor(image, cv2.COLOR_BGR2GRAY)

# 01: 提取颜色特征（灰度直方图）
color_hist = cv2.calcHist([gray], [0], None, [256], [0, 256])  # 单通道直方图
color_hist = cv2.normalize(color_hist, None).flatten()  # 修正归一化参数

# 02：提取纹理特征（GLCM）
glcm = graycomatrix(gray, 
                   distances=[1], 
                   angles=[0, np.pi/4, np.pi/2],  # 增加角度多样性
                   levels=256,
                   symmetric=True, 
                   normed=True)
contrast = np.mean(graycoprops(glcm, 'contrast'))
homogeneity = np.mean(graycoprops(glcm, 'homogeneity'))
energy = np.mean(graycoprops(glcm, 'energy'))

# 03: 优化边缘检测参数
edges = cv2.Canny(gray, threshold1=50, threshold2=150)  # 调整阈值适应太空图像

# 可视化设置
plt.figure(figsize=(15, 5))

# 原始灰度图像
plt.subplot(131)
plt.imshow(gray, cmap='gray')
plt.title('原始灰度图像')

# 颜色直方图（灰度分布）
plt.subplot(132)
plt.bar(range(256), color_hist.ravel()[:256], width=1.0)  # 显示前256个bin
plt.title('灰度直方图分布')
plt.xlim([0, 256])

# 边缘检测结果
plt.subplot(133)
plt.imshow(edges, cmap='gray')
plt.title('优化边缘检测')

plt.tight_layout()
plt.show()

# 特征输出
print("\n特征提取结果：")
print(f"直方图维度: {len(color_hist)} (0-255灰度级分布)")
print(f"纹理特征 - 平均对比度: {contrast:.2f}")
print(f"         平均同质性: {homogeneity:.2f}")
print(f"         平均能量: {energy:.2f}")
print(f"边缘特征 - 边缘像素占比: {np.mean(edges > 0)*100:.1f}%")