skimage.segmentation 的mark_boundaries怎么用

时间: 2023-10-01 14:03:16 浏览: 432
`mark_boundaries` 是 `skimage.segmentation` 中的一个函数,用于在图像的边界处绘制轮廓线。它可以接受两个参数,第一个参数是原始图像,第二个参数是对应的标签图像,它们的形状应该相同。标签图像中每个像素的值表示原始图像中对应像素所属的区域编号。 下面是一个使用 `mark_boundaries` 函数的例子: ```python from skimage.segmentation import mark_boundaries from skimage.io import imread # 读入一张图像 img = imread('image.jpg') # 对图像进行分割并得到标签图像 labels = segmentation.slic(img, n_segments=100, compactness=10) # 在图像边界处绘制轮廓线 boundaries = mark_boundaries(img, labels, color=(1, 0, 0)) # 将结果显示出来 plt.imshow(boundaries) plt.show() ``` 在这个例子中,我们首先读入一张图像,并使用 SLIC 算法对其进行分割,得到标签图像。然后,使用 `mark_boundaries` 函数在图像边界处绘制轮廓线,并将结果显示出来。 需要注意的是,`mark_boundaries` 函数默认使用黑色绘制轮廓线,如果需要使用其他颜色,可以通过 `color` 参数指定。在这个例子中,我们将轮廓线的颜色设置为红色。
相关问题

skimage.segmentation的mark_boundary如何画方框计算长宽

`skimage.segmentation.mark_boundary`函数是用来在图像上标记出分割边界的,它并不会画方框或计算长宽。如果您想要画方框并计算长宽,可以使用下面的代码: ```python import numpy as np import matplotlib.pyplot as plt from skimage import segmentation # 读入图像 image = plt.imread('image.png') # 做分割 labels = segmentation.slic(image, compactness=30, n_segments=400) # 画出边界 boundary = segmentation.mark_boundaries(image, labels, color=(1, 0, 0)) # 找到最小外接矩形的左上角坐标、宽度和高度 regions = segmentation.find_boundaries(labels) nonzero_rows, nonzero_cols = np.nonzero(regions) x, y = np.min(nonzero_cols), np.min(nonzero_rows) width, height = np.max(nonzero_cols) - x, np.max(nonzero_rows) - y # 在图片上画出最小外接矩形 plt.imshow(boundary) plt.plot([x, x+width], [y, y], 'r', linewidth=2) plt.plot([x+width, x+width], [y, y+height], 'r', linewidth=2) plt.plot([x+width, x], [y+height, y+height], 'r', linewidth=2) plt.plot([x, x], [y+height, y], 'r', linewidth=2) plt.show() # 输出长宽 print('width:', width) print('height:', height) ``` 这里使用了`skimage.segmentation.find_boundaries`函数找到分割边界,然后找到最小外接矩形的左上角坐标和宽度、高度。最后使用`matplotlib`画出最小外接矩形。

有两张大小相同的图像A和B,利用代码:from skimage.segmentation import slic from skimage.segmentation import mark_boundaries from skimage.util import img_as_float import matplotlib.pyplot as plt import numpy as np import cv2 # args args = {"image": './1.png'} # load the image and apply SLIC and extract (approximately) # the supplied number of segments image = cv2.imread(args["image"]) segments = slic(img_as_float(image), n_segments=100, sigma=5) # show the output of SLIC fig = plt.figure('Superpixels') ax = fig.add_subplot(1, 1, 1) ax.imshow(mark_boundaries(img_as_float(cv2.cvtColor(image, cv2.COLOR_BGR2RGB)), segments)) plt.axis("off") plt.show() print("segments:\n", segments) print("np.unique(segments):", np.unique(segments)) # loop over the unique segment values for (i, segVal) in enumerate(np.unique(segments)): # construct a mask for the segment print("[x] inspecting segment {}, for {}".format(i, segVal)) mask = np.zeros(image.shape[:2], dtype="uint8") mask[segments == segVal] = 255 print(mask.shape) # show the masked region cv2.imshow("Mask", mask) cv2.imshow("Applied", np.multiply(image, cv2.cvtColor(mask, cv2.COLOR_GRAY2BGR) > 0)) cv2.waitKey(0),对A进行超像素分割,将A划分的每个超像素块范围进行记录,应用到B上,使B直接得到超像素图像。给出pytorch实现代码

以下是使用PyTorch实现将图像B应用超像素图像A的代码: ```python import torch import numpy as np import cv2 import matplotlib.pyplot as plt from skimage.segmentation import slic, mark_boundaries from skimage.util import img_as_float # Load the images image_a_path = './image_a.png' image_b_path = './image_b.png' image_a = cv2.imread(image_a_path) image_b = cv2.imread(image_b_path) # Generate superpixels for image A segments_a = slic(img_as_float(image_a), n_segments=100, sigma=5) # Loop over the unique segment values for image A for seg_val in np.unique(segments_a): # Construct a mask for the segment mask = np.zeros(image_a.shape[:2], dtype="uint8") mask[segments_a == seg_val] = 255 # Apply the mask to image B masked_image_b = np.multiply(image_b, cv2.cvtColor(mask, cv2.COLOR_GRAY2BGR) > 0) # Convert the masked image to a PyTorch tensor tensor_image = torch.from_numpy(masked_image_b.astype(np.float32)).permute(2, 0, 1) / 255.0 # TODO: Apply your PyTorch model to the tensor_image # Convert the output tensor to a masked image output_image = tensor_image.permute(1, 2, 0).numpy() * 255.0 output_image = output_image.astype(np.uint8) output_image = np.multiply(output_image, cv2.cvtColor(mask, cv2.COLOR_GRAY2BGR) > 0) # Show the output image fig = plt.figure('Output') ax = fig.add_subplot(1, 1, 1) ax.imshow(output_image) plt.axis("off") plt.show() ``` 在上面的代码中,我们首先加载了图像A和图像B。然后,我们使用skimage.segmentation中的slic函数生成图像A的超像素图像。对于每个超像素块,我们构造一个掩膜并将其应用到图像B上。然后,我们将掩膜后的图像转换为PyTorch张量,并将其输入到您的PyTorch模型中。最后,我们将输出张量转换回图像格式,并将其与掩膜相乘以获得最终输出图像。

相关推荐

import cv2 import numpy as np import torch import torch.nn.functional as F from skimage.segmentation import slic import matplotlib.pyplot as plt from skimage.segmentation import mark_boundaries from skimage import img_as_float # 定义超像素数量 num_segments = 100 # 加载图像 A 和 B img_a = cv2.imread('img_a.jpg') img_b = cv2.imread('img_b.jpg') # 对图像 A 进行超像素分割,并获取每个超像素块的像素范围 segments_a = slic(img_as_float(img_a), n_segments=num_segments, sigma=5) pixel_ranges = [] for i in range(num_segments): mask = (segments_a == i) indices = np.where(mask)[1] pixel_range = (np.min(indices), np.max(indices)) pixel_ranges.append(pixel_range) # 将像素范围应用到图像 B 上实现超像素分割 segments_b = np.zeros_like(segments_a) for i in range(num_segments): pixel_range = pixel_ranges[i] segment_b = img_b[:, pixel_range[0]:pixel_range[1], :] segment_b = torch.from_numpy(segment_b.transpose(2, 0, 1)).unsqueeze(0).float() segment_b = F.interpolate(segment_b, size=(img_b.shape[0], pixel_range[1] - pixel_range[0]), mode='bilinear', align_corners=True) segment_b = segment_b.squeeze(0).numpy().transpose(1, 2, 0).astype(np.uint8) gray = cv2.cvtColor(segment_b, cv2.COLOR_BGR2GRAY) _, mask = cv2.threshold(gray, 0, 255, cv2.THRESH_BINARY) segments_b[np.where(mask)] = i # 可视化超像素分割结果 fig = plt.figure('Superpixels') ax = fig.add_subplot(1, 2, 1) ax.imshow(mark_boundaries(img_as_float(cv2.cvtColor(img_a, cv2.COLOR_BGR2RGB)), segments_a)) ax = fig.add_subplot(1, 2, 2) ax.imshow(mark_boundaries(img_as_float(cv2.cvtColor(img_b, cv2.COLOR_BGR2RGB)), segments_b)) plt.axis("off") plt.show(),上述代码中segments_a = slic(img_as_float(img_a), n_segments=num_segments, sigma=5)出现错误:ValueError: Cannot convert from object to float64.

import torch import torchvision.transforms as transforms import numpy as np from skimage.segmentation import slic from skimage.segmentation import mark_boundaries from skimage.filters import sobel from skimage.color import rgb2gray from PIL import Image # 超像素数量 num_segments = 100 # 加载图像 image = Image.open('test.jpg') # 转换为 PyTorch 张量 transform = transforms.ToTensor() img_tensor = transform(image).unsqueeze(0) # 转换为 Numpy 数组 img_np = img_tensor.numpy().transpose(0, 2, 3, 1)[0] # 转换为灰度图像 gray_img = rgb2gray(img_np) # 使用 SLIC 超像素分割算法 segments = slic(img_np, n_segments=num_segments, compactness=10, sigma=1) # 绘制超像素边界线 edge_img = mark_boundaries(img_np, segments) # 转换为灰度图像 gray_edge_img = rgb2gray(edge_img) # 使用 Canny 边缘检测算法 edges = sobel(gray_edge_img) edge_map = edges > np.mean(edges) # 绘制超像素范围的线 line_map = np.zeros_like(gray_img) for i in range(num_segments): line_map[segments == i] = edge_map[segments == i].max() # 将线绘制到图像上 line_img = np.zeros_like(img_np) line_img[:, :, 0] = line_map line_img[:, :, 1] = line_map line_img[:, :, 2] = line_map result_img = img_np * (1 - line_img) + line_img * np.array([1, 0, 0]) # 显示结果 result_img = (result_img * 255).astype(np.uint8) result_img = Image.fromarray(result_img) result_img.show()上述代码出现问题:alueError: zero-size array to reduction operation maximum which has no identity

import torch import torch.nn.functional as F from skimage.segmentation import slic, mark_boundaries import torchvision.transforms as transforms import numpy as np from PIL import Image import matplotlib.pyplot as plt # 加载图像 image = Image.open('3.jpg') # 转换为 PyTorch 张量 transform = transforms.ToTensor() img_tensor = transform(image).unsqueeze(0) # 将 PyTorch 张量转换为 Numpy 数组 img_np = img_tensor.numpy().transpose(0, 2, 3, 1)[0] # 使用 SLIC 算法生成超像素标记图 segments = slic(img_np, n_segments=60, compactness=10) # 可视化超像素索引映射 plt.imshow(segments, cmap='gray') plt.show() # 将超像素索引映射可视化 segment_img = mark_boundaries(img_np, segments) # 将 Numpy 数组转换为 PIL 图像 segment_img = Image.fromarray((segment_img * 255).astype(np.uint8)) # 保存超像素索引映射可视化 segment_img.save('segment_map.jpg') # 定义超像素池化函数 def superpixel_pooling(feature_map, segments): # 获取超像素数量和特征维度 n_segments = np.unique(segments).size n_channels = feature_map.shape[0] # 初始化超像素特征 pooled_features = torch.zeros((n_segments, n_channels)) # 对每个超像素内的像素特征进行聚合 for segment_id in range(n_segments): mask = (segments == segment_id).reshape(-1, 1, 1) pooled_feature = (feature_map * mask.float()).sum(dim=(1, 2)) / mask.sum() pooled_features[segment_id] = pooled_feature return pooled_features # 进行超像素池化 pooled_features = superpixel_pooling(img_tensor, segments) # 可视化超像素特征图 plt.imshow(pooled_features.transpose(0, 1), cmap='gray') plt.show(),上述代码出现问题:AttributeError: 'numpy.ndarray' object has no attribute 'float'

rom skimage.segmentation import slic, mark_boundaries import torchvision.transforms as transforms import numpy as np from PIL import Image import matplotlib.pyplot as plt # 加载图像 image = Image.open('3.jpg') # 转换为 PyTorch 张量 transform = transforms.ToTensor() img_tensor = transform(image).unsqueeze(0) # 将 PyTorch 张量转换为 Numpy 数组 img_np = img_tensor.numpy().transpose(0, 2, 3, 1)[0] # 使用 SLIC 算法生成超像素标记图 segments = slic(img_np, n_segments=60, compactness=10) # 可视化超像素索引映射 plt.imshow(segments, cmap='gray') plt.show() # 将超像素索引映射可视化 segment_img = mark_boundaries(img_np, segments) # 将 Numpy 数组转换为 PIL 图像 segment_img = Image.fromarray((segment_img * 255).astype(np.uint8)) # 保存超像素索引映射可视化 segment_img.save('segment_map.jpg') 将上述代码中引入超像素池化代码:import cv2 import numpy as np # 读取图像 img = cv2.imread('3.jpg') # 定义超像素分割器 num_segments = 60 # 超像素数目 slic = cv2.ximgproc.createSuperpixelSLIC(img, cv2.ximgproc.SLICO, num_segments) # 进行超像素分割 slic.iterate(10) # 获取超像素标签和数量 labels = slic.getLabels() num_label = slic.getNumberOfSuperpixels() # 对每个超像素进行池化操作,这里使用平均值池化 pooled = [] for i in range(num_label): mask = labels == i region = img[mask] pooled.append(region.mean(axis=0)) # 将池化后的特征图可视化 pooled = np.array(pooled, dtype=np.uint8) pooled_features = pooled.reshape(-1) pooled_img = cv2.resize(pooled_features, (img.shape[1], img.shape[0]), interpolation=cv2.INTER_NEAREST) print(pooled_img.shape) cv2.imshow('Pooled Image', pooled_img) cv2.waitKey(0),并显示超像素池化后的特征图

最新推荐

recommend-type

grpcio-1.47.0-cp310-cp310-linux_armv7l.whl

Python库是一组预先编写的代码模块,旨在帮助开发者实现特定的编程任务,无需从零开始编写代码。这些库可以包括各种功能,如数学运算、文件操作、数据分析和网络编程等。Python社区提供了大量的第三方库,如NumPy、Pandas和Requests,极大地丰富了Python的应用领域,从数据科学到Web开发。Python库的丰富性是Python成为最受欢迎的编程语言之一的关键原因之一。这些库不仅为初学者提供了快速入门的途径,而且为经验丰富的开发者提供了强大的工具,以高效率、高质量地完成复杂任务。例如,Matplotlib和Seaborn库在数据可视化领域内非常受欢迎,它们提供了广泛的工具和技术,可以创建高度定制化的图表和图形,帮助数据科学家和分析师在数据探索和结果展示中更有效地传达信息。
recommend-type

小程序项目源码-美容预约小程序.zip

小程序项目源码-美容预约小程序小程序项目源码-美容预约小程序小程序项目源码-美容预约小程序小程序项目源码-美容预约小程序小程序项目源码-美容预约小程序小程序项目源码-美容预约小程序小程序项目源码-美容预约小程序小程序项目源码-美容预约小程序v
recommend-type

MobaXterm 工具

MobaXterm 工具
recommend-type

grpcio-1.48.0-cp37-cp37m-linux_armv7l.whl

Python库是一组预先编写的代码模块,旨在帮助开发者实现特定的编程任务,无需从零开始编写代码。这些库可以包括各种功能,如数学运算、文件操作、数据分析和网络编程等。Python社区提供了大量的第三方库,如NumPy、Pandas和Requests,极大地丰富了Python的应用领域,从数据科学到Web开发。Python库的丰富性是Python成为最受欢迎的编程语言之一的关键原因之一。这些库不仅为初学者提供了快速入门的途径,而且为经验丰富的开发者提供了强大的工具,以高效率、高质量地完成复杂任务。例如,Matplotlib和Seaborn库在数据可视化领域内非常受欢迎,它们提供了广泛的工具和技术,可以创建高度定制化的图表和图形,帮助数据科学家和分析师在数据探索和结果展示中更有效地传达信息。
recommend-type

扁平风格PPT可修改ppt下载(11).zip

扁平风格PPT可修改ppt下载(11).zip
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

【实战演练】MATLAB用遗传算法改进粒子群GA-PSO算法

![MATLAB智能算法合集](https://static.fuxi.netease.com/fuxi-official/web/20221101/83f465753fd49c41536a5640367d4340.jpg) # 2.1 遗传算法的原理和实现 遗传算法(GA)是一种受生物进化过程启发的优化算法。它通过模拟自然选择和遗传机制来搜索最优解。 **2.1.1 遗传算法的编码和解码** 编码是将问题空间中的解表示为二进制字符串或其他数据结构的过程。解码是将编码的解转换为问题空间中的实际解的过程。常见的编码方法包括二进制编码、实数编码和树形编码。 **2.1.2 遗传算法的交叉和
recommend-type

openstack的20种接口有哪些

以下是OpenStack的20种API接口: 1. Identity (Keystone) API 2. Compute (Nova) API 3. Networking (Neutron) API 4. Block Storage (Cinder) API 5. Object Storage (Swift) API 6. Image (Glance) API 7. Telemetry (Ceilometer) API 8. Orchestration (Heat) API 9. Database (Trove) API 10. Bare Metal (Ironic) API 11. DNS
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。