|
马上注册,结交更多好友,享用更多功能^_^
您需要 登录 才可以下载或查看,没有账号?立即注册
x
代码:
# canny检测算子
from tifffile import imread, imwrite
from skimage import filters, feature, color
import numpy as np
# 读取TIFF格式无人机影像数据
image = imread('F:\\duneline\\dune\\dune2.tif')
# 转换为灰度图像
gray_image = color.rgb2gray(image)
# 边缘检测,提取沙丘脊线
edges = feature.canny(gray_image)
# 将阈值应用于图像,得到二值化图像
binary = edges
# 叠加沙脊线在原始影像上
overlay = np.copy(image)
overlay[binary] = [255, 0, 0] # 将沙脊线部分标记为红色
# 导出结果
imwrite('can_2.tif', overlay)
要将代码中的检测算子更改为拉普拉斯算子检测算子,您可以使用skimage库中的laplace函数来实现。以下是修改后的代码:
- from tifffile import imread, imwrite
- from skimage import filters, feature, color
- import numpy as np
- # 读取TIFF格式无人机影像数据
- image = imread('F:\\duneline\\dune\\dune2.tif')
- # 转换为灰度图像
- gray_image = color.rgb2gray(image)
- # 使用拉普拉斯算子进行边缘检测
- edges = filters.laplace(gray_image)
- # 将阈值应用于图像,得到二值化图像
- binary = edges > 0
- # 叠加沙脊线在原始影像上
- overlay = np.copy(image)
- overlay[binary] = [255, 0, 0] # 将沙脊线部分标记为红色
- # 导出结果
- imwrite('laplace_2.tif', overlay)
复制代码
在修改后的代码中,我们使用了 filters.laplace 函数来进行拉普拉斯算子边缘检测。然后,将得到的边缘图像进行阈值处理,得到二值化图像。最后,我们将二值化图像叠加到原始影像上,并将沙脊线部分标记为红色。最终结果将保存为"laplace_2.tif"文件。
|
|