【Python语义分割】Segment Anything(SAM)模型交互式分割+掩膜保存(三)

news2024/11/25 5:42:12

我之前分享了Segment Anything(SAM)模型的基本操作,这篇给大家分享下交互式语义分割代码,可以通过鼠标点击目标物生成对应的掩膜,同时我还加入了掩膜保存的代码。

1 Segment Anything介绍

1.1 概况

        Meta AI 公司的 Segment Anything 模型是一项革命性的技术,该模型能够根据文本指令或图像识别,实现对任意物体的识别和分割。这一模型的推出,将极大地推动计算机视觉领域的发展,并使得图像分割技术进一步普及化。

        论文地址:https://arxiv.org/abs/2304.02643

        项目地址:Segment Anything

1.2 使用方法

        具体使用方法上,Segment Anything 提供了简单易用的接口,用户只需要通过提示,即可进行物体识别和分割操作。例如在图片处理中,用户可以通过 Hover & Click 或 Box 等方式来选取物体。值得一提的是,SAM 还支持通过上传自己的图片进行物体分割操作,提取物体用时仅需数秒。

        总的来说,Meta AI 的 Segment Anything 模型为我们提供了一种全新的物体识别和分割方式,其强大的泛化能力和广泛的应用前景将极大地推动计算机视觉领域的发展。未来,我们期待看到更多基于 Segment Anything 的创新应用,以及在科学图像分析、照片编辑等领域的广泛应用。

2 交互式分割代码实现

2.1 用于生成显示掩膜的函数(初始化)

        里面包含三个封装好的函数,一个是生成掩膜(分割的轮廓)的函数,一个是显示标记点(自己选择需要分割的目标)的函数。

def show_mask(mask, ax, random_color=False):
    if random_color:
        color = np.concatenate([np.random.random(3), np.array([0.6])], axis=0)
    else:
        color = np.array([30 / 255, 144 / 255, 255 / 255, 0.6])
    h, w = mask.shape[-2:]
    mask_image = mask.reshape(h, w, 1) * color.reshape(1, 1, -1)
    ax.imshow(mask_image)


def show_points(coords, labels, ax, marker_size=25):
    pos_points = coords[labels == 1]
    neg_points = coords[labels == 0]
    ax.scatter(pos_points[:, 0], pos_points[:, 1], color='green', marker='o', s=marker_size, edgecolor='white',
               linewidth=1.25)
    ax.scatter(neg_points[:, 0], neg_points[:, 1], color='red', marker='o', s=marker_size, edgecolor='white',
               linewidth=1.25)

2.2 交互事件&分割掩膜

        这里在原始单点生成掩膜的基础上加入了鼠标点击的事件,可以实现点击后自动分割该点的目标,同时右键可以清楚当前窗口所有的掩膜。

def select_point(event):
    # 鼠标点击事件
    global predictor, img_raster, image  # 声明全局变量
    if event.button == 1:
        x, y = event.xdata, event.ydata
        input_point = np.array([[int(x),int(y)]])
        # 单点 prompt  输入格式为(x, y)和并表示出点所带有的标签1(前景点)或0(背景点)。
        input_label = np.array([1])  # 点所对应的标签
        masks, scores, logit = predictor.predict(
            point_coords=input_point,
            point_labels=input_label,
            multimask_output=False,  # 为False时,它将返回一个掩码
        )
        # print(masks.shape)  # (3, 2160, 3840)波段,高,宽
        for i, (mask, score) in enumerate(zip(masks, scores)):
            # 三个置信度不同的图
            show_mask(mask, plt.gca())
            img_raster[mask] = 1
            show_points(input_point, input_label, plt.gca())
            print("[%s]坐标:%s,%s已分割完成......" % (datetime.now().strftime('%Y-%m-%d %H:%M:%S'), int(x), int(y)))
    if event.button == 3:  # 如果点击了鼠标右键
        print("[%s]正在清除内容......" % datetime.now().strftime('%Y-%m-%d %H:%M:%S'))
        img_raster = np.zeros((image.shape[0], image.shape[1]))  # 覆盖原来的掩膜,全新的0数组
        plt.clf()  # 清楚所有图形
        plt.imshow(image)  # 显示原始图片
        plt.axis('off')  # 不显示坐标系
    plt.draw()  # 更新窗口中的图片

2.3 基础函数&模型加载

        这里是分割的图片加载、模型加载、参数设置,里面链接了之前的鼠标点击事件,可以整体实现左键点击分割目标,右键删除所有目标。

global predictor, img_raster, image  # 声明全局变量
print("【程序准备阶段】")
print("[%s]正在读取图片......" % datetime.now().strftime('%Y-%m-%d %H:%M:%S'))
image = cv2.imread(image_path)  # 读取的图像以NumPy数组的形式存储在变量image中
print("[%s]正在转换图片格式......" % datetime.now().strftime('%Y-%m-%d %H:%M:%S'))
image = cv2.cvtColor(image, cv2.COLOR_BGR2RGB)  # 将图像从BGR颜色空间转换为RGB颜色空间,还原图片色彩(图像处理库所认同的格式)
print("[%s]正在初始化模型参数......" % datetime.now().strftime('%Y-%m-%d %H:%M:%S'))
img_raster = np.zeros((image.shape[0], image.shape[1]))
# 创建一个二维数组,用于保存掩膜做栅格转面
fit = plt.figure(figsize=(15, 15))
plt.imshow(image)
sys.path.append("..")  # 将当前路径上一级目录添加到sys.path列表,这里模型使用绝对路径所以这行没啥用
sam_checkpoint = model_path  # 定义模型路径
model_type = model_type  # 定义模型类型
device = device  # "cpu"  or  "cuda"
sam = sam_model_registry[model_type](checkpoint=sam_checkpoint)
sam.to(device=device)  # 定义模型参数
predictor = SamPredictor(sam)  # 调用预测模型
predictor.set_image(image)
# 调用`SamPredictor.set_image`来处理图像以产生一个图像嵌入。`SamPredictor`会记住这个嵌入,并将其用于随后的掩码预测
print("【单点分割阶段】")
print("[%s]正在分割图片......" % datetime.now().strftime('%Y-%m-%d %H:%M:%S'))
fit.canvas.mpl_connect('button_press_event', select_point)
# 窗口点击事件
plt.axis('off')
plt.show()
print("【结果保存阶段】")
print("[%s]正在保存分割数据......" % datetime.now().strftime('%Y-%m-%d %H:%M:%S'))
driver = gdal.GetDriverByName('GTiff')  # 载入数据驱动,用于存储内存中的数组
ds_result_raster = driver.Create(out_path, image.shape[1], image.shape[0], bands=1, eType=gdal.GDT_Float64)
ds_result_raster.GetRasterBand(1).SetNoDataValue(0)  # 将无效值设为0
ds_result_raster.GetRasterBand(1).WriteArray(img_raster)  # 将结果写入数组

3 实现效果

        这里我选择了多个点,每个点都是通过鼠标左键点击生成的掩膜,单点鼠标右键可以清空整个窗口的所有掩膜。事实上,实现多点、正负点、矩形+负点的交互式操作会大大减少我们的工作量,但我没什么精力所以就没往后继续写了。大家如果有兴趣可以自己去实现,然后分享给我(手动偷笑)。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1101560.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

HarmonyOS 音视频开发概述

在音视频开发指导中,将介绍各种涉及音频、视频播放或录制功能场景的开发方式,指导开发者如何使用系统提供的音视频 API 实现对应功能。比如使用 TonePlayer 实现简单的提示音,当设备接收到新消息时,会发出短促的“滴滴”声&#x…

【API篇】三、转换算子API(上)

文章目录 0、demo数据1、基本转换算子:映射map2、基本转换算子:过滤filter3、基本转换算子:扁平映射flatMap4、聚合算子:按键分区keyBy5、聚合算子:简单聚合sum/min/max/minBy/maxBy6、聚合算子:归约聚合re…

第三章 内存管理 七、具有快表的地址变换结构

目录 一、什么是快表 二、快表有什么用? 例子: 三、快表和慢表同时查询 四、局部性原理 五、总结 一、什么是快表 快表,又称联想寄存器(TLB,translation lookaside buffer),是一种访问速度比内存快很…

教程更新 | 持续开源 RK3568驱动指南-驱动基础进阶篇

《iTOP-RK3568开发板驱动开发指南》手册文档更新,手册内容对应视频教程,后续资料会不断更新,不断完善,帮助用户快速入门,大大提升研发速度。 ✦ 第一篇 驱动基础 第1章 前言 第2章 你好!内核源码 第3章 …

如何部署WampServer并结合cpolar内网穿透工具实现公网访问本地服务?

文章目录 前言1.WampServer下载安装2.WampServer启动3.安装cpolar内网穿透3.1 注册账号3.2 下载cpolar客户端3.3 登录cpolar web ui管理界面3.4 创建公网地址 4.固定公网地址访问 前言 Wamp 是一个 Windows系统下的 Apache PHP Mysql 集成安装环境,是一组常用来…

苏东坡在元丰五年

北宋神宗元丰二年(1079年)十二月二十八日,震惊朝野上下的“乌台诗案” 正式结案:“祠部员外郎、直史馆苏轼责授检校水部员外郎、黄州团练副使,本州安置、不得签书公事,令御史台差人转押前去。” 元丰三年&…

PAM从入门到精通(三)

接前一篇文章:PAM从入门到精通(二) 本文参考以下博文: PAM模块详解及sudo命令 PAM 的应用开发和内部实现源码分析 PAM详解(一)PAM介绍 Linux进阶_PAM认证机制 Oracle Solaris 10 开发者安全性指南 ——…

[Jdk版本不一致问题 ]终端查看jdk版本不一致

目录 前言: 原因 是因为环境变量配置在读取的时候 顺序不一致,或者有多个jdk版本 安装 去对应环境变量 里面查找对应 jdk配置的路径信息: 前言: cmd 终端 在查看java版本的时候, 发现 使用java -version 显示的是 :1.8 版本, javac -version 显示的是 javac11 版本 原因 是因…

Linux | 一篇文章带你深刻理解粘滞位

目录 前言 一、为什么需要粘滞位 二、粘滞位是什么? 三、怎么设置粘滞位 前言 关于粘滞位的学习,首先你得了解Linux下的权限,当然,文章会稍稍带着复习关键知识点,具体关于权限,请移步 Linux权限 一、为…

关于ts的keyof

type props_type {name: string,age: number }const props: props_type {name: tjq,age: 18 }for (const key in props) { //props[key]出现红色波浪线const value props[key]; }why? 经过我查阅多方资料,在网上看到一个比较合适的例子 地址&#xf…

如何在Python中更新代码但还想保留原有代码

Python作为一门功能强大的编程语言,为开发者提供了许多方便的方法来更新代码并且还能保留原有代码。在本文中,我们将从多个方面来详细阐述如何在Python中更新代码但还想保留原有代码。 一、使用函数参数 许多Python程序员通过将函数的参数作为字典或者…

Selenium工作原理详解

前言 Selenium是一个用于Web应用程序自动化测试工具。Selenium测试直接运行在浏览器中,就像真正的用户在操作一样。支持的浏览器包括IE(7, 8, 9, 10, 11),Mozilla Firefox,Safari,Google Chrome&#xff0…

css-雷达扫描附近好友

效果图 html <!DOCTYPE html> <html><head><meta charset"utf-8"><title></title><link type"text/css" rel"stylesheet" href"dy.css"/></head><body><div class"r…

pg嵌套子查询

1.概念 查询里面还有查询 进阶版&#xff1a;关联子查询 2.相关运算符补充 in/all/any all&#xff1a;表中的所有内容遍历一边&#xff0c;等价与max some/any&#xff1a;表中任何一个&#xff0c;等价与min

第三章 内存管理 六、基本地址变换结构

目录 一、定义 二、物理地址计算步骤 1、设页面大小为L&#xff0c;逻辑地址A到物理地址E的变换过程如下: 2、例子 三、总结 一、定义 基本地址变换结构是指计算机系统中用于实现虚拟地址到物理地址转换的基本硬件结构。 基本地址变换结构通常包括以下几个部分&#xff1…

【Linux内核】内存映射原理

【Linux内核】内存映射原理 物理地址空间 物理地址是处理器在总线上能看到的地址,使用RISC(Reduced Instruction Set Computing精简指令集)的处理器通常只实现一个物理地址空间,外围设备和物理内存使用统一的物理空间, 有些架构的处理器把分配给外围设备的物理地址称为设备内存…

基于matlab/simulink永磁直驱风力发电低电压穿越控制仿真模型

很久没更新了&#xff0c;下面继续更新电能质量治理这块的内容&#xff0c;今天分享的直驱式风力发电LVRT控制。 随着风力发电规模不断扩大&#xff0c;风电在电网系统中所占份额逐渐增加。当风电并网运行时&#xff0c;会发生很多种并网问题&#xff0c;其中最常见的问题之一就…

macos 12 支持机型 macOS Monterey 更新中新增的功能

macOS Monterey 能让你以全然一新的方式与他人沟通联络、共享内容和挥洒创意。尽享 FaceTime 通话新增的音频和视频增强功能&#xff0c;包括空间音频和人像模式。通过功能强大的效率类工具&#xff08;例如专注模式、快速备忘录和 Safari 浏览器中的标签页组&#xff09;完成更…

软件产品质量如何保障?找对软件测试公司是关键

软件产品的质量对于企业来说至关重要&#xff0c;它直接关系到用户的满意度和企业的声誉。然而&#xff0c;由于软件开发的复杂性和多样性&#xff0c;确保软件产品质量成为一项挑战。为了解决这一问题&#xff0c;找对合适的软件测试公司是至关重要的。 软件测试公司具备丰富…

考试成绩怎么查?

如何创建一个学生自助查询考试成绩的系统&#xff1f;还包括一个网页界面&#xff1f;能够让学生们通过多种方式轻松查询他们的成绩&#xff1f; 我先把最简单的方法告诉老师们&#xff0c;使用易查分&#xff0c;它可以帮助老师快速发布成绩&#xff0c;学生轻松查询&#xff…