代码最终运行结果为什么是全黑的，没有显示结果。麻烦帮忙改正代码。
代码：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 = # 将沙脊线部分标记为红色
# 配置输出的空间参考信息
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, image.shape, 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 于 2023-8-13 18:49 编辑
问题出在这行代码上：
image_path = 'F:\duneline\dune\dune.tif'
在Python中，\是一个转义字符，用于表示特殊字符。因此，\d在字符串中被解释为一个特殊字符，而不是两个普通的字符。这会导致路径无法正确解析，从而无法读取图像数据。
为了解决这个问题，你可以使用原始字符串（raw string）来表示路径。原始字符串以r开头，可以让字符串中的转义字符失效。修改代码如下：
image_path = r'F:\duneline\dune\dune.tif'
这样，\d将被解释为两个普通的字符，而不是一个特殊字符。这样就可以正确读取图像数据了。 本帖最后由 陶远航 于 2023-8-13 17:30 编辑
加个r 陶远航 发表于 2023-8-13 17:28
问题可能出在代码的最后一行，print("保存成功！")。根据提供的代码，你使用了中文的引号字符（U+3003）而 ...
陶远航 发表于 2023-8-13 17:28
加个r
Snipaste 会记录截图记录的哦~
歌者文明清理员 发表于 2023-8-13 17:34
Snipaste 会记录截图记录的哦~
举报又没奖励，还没人管。
要不然，我靠举报都要发财了。
页:
[1]