代码最终运行结果为什么是全黑的，没有显示结果。麻烦帮忙改正代码。

Justheartyoung

代码：

from tifffile import imread, imwrite
from skimage import filters, feature, color
import matplotlib.pyplot as plt
import numpy as np
from osgeo import gdal, osr
# 读取TIFF格式无人机影像数据
image_path = 'F:\duneline\dune\dune.tif'
output_path = 'overlay_image5.tif'
image = imread(image_path)

# 转换为灰度图像
gray_image = color.rgb2gray(image)
# 边缘检测，提取沙丘脊线
edges = filters.sobel(gray_image)
# 设置合适的阈值来确定沙脊线的二值化图像
threshold = 0.03
binary = edges > threshold
# 叠加沙脊线在原始影像上
overlay = np.copy(image)
overlay[binary] = [255, 0, 0]  # 将沙脊线部分标记为红色
# 配置输出的空间参考信息
in_ds = gdal.Open(image_path)
projection = in_ds.GetProjection()
geotransform = in_ds.GetGeoTransform()

# 创建新的输出图像
driver = gdal.GetDriverByName('GTiff')
out_ds = driver.Create(output_path, image.shape[1], image.shape[0], 1, gdal.GDT_Byte)  # 修改数据类型为gdal.GDT_Byte

# 设置输出图像的空间参考信息
out_ds.SetProjection(projection)
out_ds.SetGeoTransform(geotransform)

# 将沙脊线部分写入输出图像
out_band = out_ds.GetRasterBand(1)
out_band.WriteArray(binary.astype(np.uint8))  # 修改为正确的数据类型

# 设置颜色表
color_table = gdal.ColorTable()
color_table.SetColorEntry(0, (0, 0, 0, 0))  # 黑色，对应值为0
color_table.SetColorEntry(255, (255, 255, 255, 255))  # 白色，对应值为255
out_band.SetColorTable(color_table)

# 释放资源
out_band.FlushCache()
out_ds = None
in_ds = None

print("保存成功！")

isdkz

问题出在这行代码上：

image_path = 'F:\duneline\dune\dune.tif'



在Python中，\是一个转义字符，用于表示特殊字符。因此，\d在字符串中被解释为一个特殊字符，而不是两个普通的字符。这会导致路径无法正确解析，从而无法读取图像数据。



为了解决这个问题，你可以使用原始字符串（raw string）来表示路径。原始字符串以r开头，可以让字符串中的转义字符失效。修改代码如下：



image_path = r'F:\duneline\dune\dune.tif'



这样，\d将被解释为两个普通的字符，而不是一个特殊字符。这样就可以正确读取图像数据了。

陶远航

加个r

歌者文明清理员

歌者文明清理员

陶远航 发表于 2023-8-13 17:28
问题可能出在代码的最后一行，print("保存成功！")。根据提供的代码，你使用了中文的引号字符（U+3003）而 ...

歌者文明清理员

陶远航 发表于 2023-8-13 17:28
加个r

Snipaste 会记录截图记录的哦~

ba21

歌者文明清理员 发表于 2023-8-13 17:34
Snipaste 会记录截图记录的哦~

举报又没奖励，还没人管。
要不然，我靠举报都要发财了。