Python(TensorFlow和PyTorch)两种显微镜成像重建算法模型(显微镜学)

news2024/9/21 3:29:16

🎯要点

🎯受激发射损耗显微镜算法模型:🖊恢复嘈杂二维和三维图像 | 🖊模型架构:恢复上下文信息和超分辨率图像 | 🖊使用嘈杂和高信噪比的图像训练模型 | 🖊准备半合成训练集 | 🖊优化沙邦尼尔损失和边缘损失 | 🖊使用峰值信噪比、归一化均方误差和多尺度结构相似性指数量化结果 | 🎯训练荧光显微镜模型和对抗网络图形转换模型

🍪语言内容分比

在这里插入图片描述
在这里插入图片描述

🍇Python图像归一化

在图像处理中,归一化是改变像素强度值范围的过程。例如,应用包括由于眩光而对比度较差的照片。归一化有时称为对比度拉伸或直方图拉伸。在更一般的数据处理领域(例如数字信号处理),它被称为动态范围扩展。

在各种应用中,动态范围扩展的目的通常是将图像或其他类型的信号带入感官更熟悉或正常的范围,因此称为归一化。通常,其动机是使一组数据、信号或图像的动态范围保持一致,以避免精神分散或疲劳。例如,报纸会努力使一期中的所有图像都具有相似的灰度范围。

归一化对 n 维灰度图像 I : { X ⊆ R n } → { I:\left\{ X \subseteq R ^n\right\} \rightarrow\{ I:{XRn}{ Min, …, Max } \} } 进行变换,强度值在 (Min, Max) 范围内),转换为新图像 I N : { X ⊆ R n } → { I_N:\left\{ X \subseteq R ^n\right\} \rightarrow\{ IN:{XRn}{ newMin, …, newMax } \} },强度值在 (newMin, newMax) 范围内。灰度数字图像的线性归一化根据以下公式进行:
I N = ( I − Min ⁡ )  newMax  −  newMin  Max ⁡ −  Min  +  newMin  I_N=(I-\operatorname{Min}) \frac{\text { newMax }- \text { newMin }}{\operatorname{Max}-\text { Min }}+\text { newMin } IN=(IMin)Max Min  newMax  newMin + newMin 
例如,如果图像的强度范围是 50 到 180,而所需范围是 0 到 255,则该过程需要从每个像素强度中减去 50,使得范围为 0 到 130。然后将每个像素强度乘以 255/130,使得范围为 0 到 255。

归一化也可能是非线性的,当 I I I I N I_N IN 之间不存在线性关系时就会发生这种情况。非线性归一化的一个例子是当归一化遵循 sigmoid 函数时,在这种情况下,归一化图像根据以下公式计算:
I N = (  newMax  −  newMin  ) 1 1 + e − I − β α +  newMin  I_N=(\text { newMax }- \text { newMin }) \frac{1}{1+e^{-\frac{I-\beta}{\alpha}}}+\text { newMin } IN=( newMax  newMin )1+eαIβ1+ newMin 
其中 α \alpha α 定义输入强度范围的宽度, β \beta β 定义输入强度范围的中心强度。

图像处理软件中的自动归一化通常归一化为图像文件格式中指定的数字系统的完整动态范围。

PyTorch归一化图像

PyTorch 中的标准化是使用 torchvision.transforms.Normalize() 完成的。这用平均值和标准差对张量图像进行归一化


from PIL import Image
import matplotlib.pyplot as plt
import numpy as np

img_path = 'Koa.jpg'
img = Image.open(img_path)

img_np = np.array(img)

plt.hist(img_np.ravel(), bins=50, density=True)
plt.xlabel("pixel values")
plt.ylabel("relative frequency")
plt.title("distribution of pixels")

从结果,我们发现RGB图像的像素值范围是0到255。

使用 ToTensor() 将 PIL 图像转换为 PyTorch 张量,并绘制该张量图像的像素值。我们定义变换函数将 PIL 图像转换为 PyTorch 张量图像。

import torchvision.transforms as transforms
import matplotlib.pyplot as plt

transform = transforms.Compose([
	transforms.ToTensor()
])

img_tr = transform(img)
img_np = np.array(img_tr)

plt.hist(img_np.ravel(), bins=50, density=True)
plt.xlabel("pixel values")
plt.ylabel("relative frequency")
plt.title("distribution of pixels")

从结果,我们发现张量图像的像素值范围从0.0到1.0。我们注意到 RBG 和张量图像的像素分布看起来相同,但像素值范围不同。

我们计算图像的平均值和标准差:

img_tr = transform(img)
mean, std = img_tr.mean([1,2]), img_tr.std([1,2])

print("mean and std before normalize:")
print("Mean of the image:", mean)
print("Std of the image:", std)

标准化前的均值和标准差: 
图像的均值:tensor([0.4916, 0.4498, 0.4000]) 
图像的标准差:tensor([0.2474, 0.2362, 0.2322])

为了归一化图像,我们在这里使用上面计算的图像平均值和标准差。如果图像与 ImageNet 图像相似,我们也可以使用 ImageNet 数据集的平均值和标准差。ImageNet 的平均值和标准差为:平均值 = [0.485, 0.456, 0.406] 和标准差 = [0.229, 0.224, 0.225]。如果图像与 ImageNet 不相似(如医学图像),则始终建议计算数据集的平均值和标准差并使用它们来归一体图像。

from torchvision import transforms

transform_norm = transforms.Compose([
	transforms.ToTensor(),
	transforms.Normalize(mean, std)
])

img_normalized = transform_norm(img)
img_np = np.array(img_normalized)

plt.hist(img_np.ravel(), bins=50, density=True)
plt.xlabel("pixel values")
plt.ylabel("relative frequency")
plt.title("distribution of pixels")

我们已经用计算出的平均值和标准差对图像进行了归一化。上面的输出显示了归一化图像的像素值分布。我们可以注意到张量图像(归一化之前)和归一化图像的像素分布之间的差异。

现在可视化标准化图像:

img_normalized = transform_norm(img)
img_normalized = np.array(img_normalized)
img_normalized = img_normalized.transpose(1, 2, 0)

plt.imshow(img_normalized)
plt.xticks([])
plt.yticks([])

👉参阅、更新:计算思维 | 亚图跨际

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

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

相关文章

基于单片机的水产养殖饲料自动投喂系统

文章目录 前言资料获取设计介绍功能介绍设计清单具体实现截图系统框架图设计获取 前言 💗博主介绍:✌全网粉丝10W,CSDN特邀作者、博客专家、CSDN新星计划导师,一名热衷于单片机技术探索与分享的博主、专注于 精通51/STM32/MSP430/AVR等单片机…

基于SpringBoot的扶贫助农管理系统

作者:计算机学姐 开发技术:SpringBoot、SSM、Vue、MySQL、JSP、ElementUI等,“文末源码”。 专栏推荐:前后端分离项目源码、SpringBoot项目源码、SSM项目源码 系统展示 基于JavaSpringBootVueMySQL的扶贫助农管理系统【附源码文档…

《JavaEE进阶》----14.<SpringMVC配置文件实践之【验证码项目】>

本篇博客介绍的是Google的开源项目Kaptcha来实现的验证码。 这种是最简单的验证码。 也是很常见的一种验证码。可以去看项目结果展示。就可以明白这个项目了。 前言: 随着安全性的要求越来越高、很多项目都使用了验证码。如今验证码的形式也是有许许多多、更复杂的图…

(有直接运行的源码)SpringBoot-MyBatis in查询 <foreach> 范围查询, 以及其他方法查询

文章目录 第一种: names入参为数组,foreach标签遍历第二种: name入参是纯字符串, 代替数据中的内容需要的工具类第三种: name入参为纯字符串, foreach标签分割遍历第四种: mysql中使用 (待补充)第五种: 存储过程和函数,以及其他数据库(待补充)源代码地址第一种: names入参为数组…

王道考研操作系统笔记(一)

虚拟内存的定义和特征: 基于局部性的原理, 在程序装入时,可以将程序中很快用到的部分装入内存,暂时用不到的数据装入外存,就可以让程序开始执行,在程序执行过程中,当所访问的信息不在内存的时…

以租代卖+押金原路退回系统 阿雪技术观市场策划

一、市场目标 1. 提高软件的知名度和用户使用率。 2. 促进物品的流转和租赁,增加用户参与度和满意度。 3. 探索新的商业模式,为公司带来更多的收益和发展机会。 二、实施步骤 1. 软件平台搭建 - 开发和优化市场策划软件,确保其功能完善…

数据库系列之GaussDB数据库中逻辑对象关系简析

初次接触openGauss或GaussDB数据库的逻辑对象,被其中的表空间、数据库、schema和用户之间的关系,以及授权管理困惑住了,与熟悉的MySQL数据库的逻辑对象又有明显的不同。本文旨在简要梳理下GaussDB数据库逻辑对象之间的关系,以加深…

【高级编程】Java IO流(上)字节流 InputStream OutputStream

文章目录 文件操作流输入流InputStream输出流OutputStream 文件操作 文件是指相关记录或放在一起的数据的集合。是一种用于存储数据的基本单位,它可以包含各种类型的信息,例如文本、图像、音频或视频。文件在计算机中通常存储在磁盘或其他存储介质上&am…

华为OD机试 - 最大相连男生数 - 矩阵(Java 2024 E卷 200分)

华为OD机试 2024E卷题库疯狂收录中,刷题点这里 专栏导读 本专栏收录于《华为OD机试(JAVA)真题(E卷D卷A卷B卷C卷)》。 刷的越多,抽中的概率越大,私信哪吒,备注华为OD,加…

SprinBoot+Vue健康健身追踪系统的设计与实现

目录 1 项目介绍2 项目截图3 核心代码3.1 Controller3.2 Service3.3 Dao3.4 application.yml3.5 SpringbootApplication3.5 Vue 4 数据库表设计5 文档参考6 计算机毕设选题推荐7 源码获取 1 项目介绍 博主个人介绍:CSDN认证博客专家,CSDN平台Java领域优质…

基于SpringBoot+Vue社区家政服务预约平台【提供源码+答辩PPT+参考文档+项目部署】

项目技术架构: 该家政服务预约平台基于B/S架构,后端采用SpringBoot架构,前端采用VueElementUI技术,并结合MySQL数据库。该系统分为系统管理员、家政人员、雇主用户三种角色。 系统核心功能: 用户管理、家政人员管理…

Zenmap

Zenmap Zenmap 本身作为一款端口扫描的工具无法直接判断一台服务器是否存在具体的漏洞,但它可以帮助你识别服务器上的开放端口、运行的服务以及操作系统信息,通过这些信息,需要你自己凭借经验推测可能存在的漏洞,或者结合其他工具…

微服务CI/CD实践(六)Jenkins Docker 自动化构建部署Java微服务

微服务CI/CD实践系列: 微服务CI/CD实践(一)环境准备及虚拟机创建 微服务CI/CD实践(二)服务器先决准备 微服务CI/CD实践(三)gitlab部署及nexus3部署 微服务CI/CD实践(四&#xff09…

log4j日志封装说明—slf4j对于log4j的日志封装-正确获取调用堆栈

日志是项目中必用的东西,日志产品里最普及应该就是log4j了。(logback这里暂不讨论。) 先看一下常用的log4j的用法,一般来说log4j都会配合slf4j或者common-logging使用,这里已slf4j为例。添加gradle依赖: dependencies { compile(l…

When is a Liability not a Liability

When is a Liability not a Liability? Textual Analysis, Dictionaries, and 10-Ks(2011) 论文阅读 文章目录 When is a Liability not a Liability? Textual Analysis, Dictionaries, and 10-Ks(2011) 论文阅读 Abstract哈佛H4N List 与 Fin-Neg如何创建词汇列表&#xff1…

声音克隆工具CosyVoice

阿里的免费声音克隆工具CosyVoice CosyVoice 是阿里通义实验室在七月初开源的一款专注于自然语音合成的语音大模型,它具备多语言、多音色和细腻的情感控制能力。这个系统支持中文、英文、日文、粤语和韩语五种语言的语音生成,并且在语音合成的效果上远超…

【python】SemEval 2014数据集 xml文件格式转换为csv+txt

笔记为自我总结整理的学习笔记,若有错误欢迎指出哟~ 【python】SemEval 2014数据集 xml文件格式转换为csvtxt SemEval2014简介4个子任务数据格式xml转csvxml转txt SemEval2014简介 SemEval2014,ABSA( Aspect Based Sentiment Analysis&#…

Dance with Compiler - EP3 ARM64 汇编传参约定以及 restrict 汇编分析

在 ARM64 架构(也称为 AArch64)中,函数调用约定定义了寄存器如何用于传递参数和返回值。这些约定有助于实现高效的函数调用和返回。在 ARM64 的汇编中,寄存器传参遵循以下约定: 参数传递寄存器 x0 - x7: 这 8 个寄存器…

Golang | Leetcode Golang题解之第396题旋转函数

题目: 题解: func maxRotateFunction(nums []int) int {numSum : 0for _, v : range nums {numSum v}f : 0for i, num : range nums {f i * num}ans : ffor i : len(nums) - 1; i > 0; i-- {f numSum - len(nums)*nums[i]ans max(ans, f)}return…

超声波气象站

超声波气象站的设计优势包括: 非接触式测量:超声波气象站利用超声波技术进行测量,可以在不接触被测物体的情况下进行精确的测量。这意味着测量过程不会受到外界干扰或影响,提高了测量的准确性和可靠性。 多参数测量:超…