自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+
  • 博客(67)
  • 资源 (3)
  • 收藏
  • 关注

原创 深入理解PSNR(峰值信噪比)(附matlab代码)

深入理解PSNR作者:老李日期:2022-1-19本文引入MSE、SNR、变异系数(Coefficient of Variation),并希望从统计学的角度上解释这个变量这个变量。PSNR是“Peak Signal to Noise Ratio”的缩写,即峰值信噪比。MSE用 真实值-预测值 然后平方之后求和平均。在图像中的MSE定义说人话就是:两个图像中每一个相同位置的像素值相减,平方,求和,再求平均。表达的是两幅图在每一个位置上的像素值的差异的平均。数值越大,表示两张图片在相似

2022-01-19 00:49:20 42597 6

原创 深入理解SSIM(两图像结构相似度指标)(附matlab代码)

SSIM的深入理解作者:老李日期:2022-1-18SSIMSSIM用于评价两张图像的相似程度。对于SSIM这个指标的设计思路,如下图所示。图像的均值表示图像的亮度。图像的方差表示图像的对比度。注意:要先将图像归一化(把像素的阈值调整为[0,1]),才可以通过计算,得出数值进行比较。当α=β=γ=1,C_3= 0.5C_2(常用),则SSIM表达式为:“对比”是怎么体现的关于,这里我觉得是大部分博文所没有涉及到的地方。该问题等价于这也是以下三个等式如何产生的。也就是说以下三个

2022-01-18 04:49:19 23776 19

原创 图像质量的定义

ISP自学笔记

2022-01-16 18:13:56 1013

原创 matplotlib 修改字体

matplotlib 修改字体2021-9-9.ttc文件首先你得有一个字体的文件**(.ttc文件)**该类型文件的描述百度百科是这么写的:TTC全称是TrueType Collection,是一个字体格式。TTC字体是TrueType字体集成文件(. TTC文件),是在一单独文件结构中包含多种字体,以便更有效地共享轮廓数据,当多种字体共享同一笔画时,TTC技术可有效地减小字体文件的大小。TTC是几个TTF合成的字库,安装后字体列表中会看到两个以上的字体。两个字体中大部分字都一样时,可以将两

2021-09-09 11:00:01 565

原创 作业一:误差

第一次作业:泰勒展开的误差分析

2020-11-16 11:28:44 202

原创 所遇到的matlab中的路径操作的一点总结

matlab中的路径操作最近在看matconvnet,我发现里面有频繁的路径操作。而对于一个比较大的程序而言,里面必然有很多的关于路径的操作,于是想写一篇分享总结一下。我们在使用matconvnet的时候,无论是test还是train,运行的一开始,会有一个run(fullfile(fileparts(mfilename('fullpath')), ... '..', '..', 'matlab', 'vl_setupnn.m')) ;(mfilename(‘fullpath’):返回当先

2020-07-25 17:13:55 861 1

原创 基本的GMRES基于cpp实现

基本的GMRES基于cpp实现理论部分和matlab版本可以查看我的其他文章:基本理论完整实现以下为代码,头文件里面有cuda的东西,

2020-07-22 11:44:18 1198 3

原创 如何将matlab生成的矩阵传入C中

从matlab生成一个对称正定矩阵并传到C里面

2020-07-11 16:51:54 1026

原创 cell与struct

matconvnetm描述某一层的基本用法是struct+cellstruct可以理解为matlab版的键值对(dictionary)cell可以理解为c里面的结构体,他的操作很随意,可以任意加东西进去

2020-07-06 18:10:04 502

原创 如何使用matlab动态描述一个过程

1.程序中简单的演示2..gif文件的制作

2020-06-28 17:26:29 673

原创 pytorch拓展-cpp extension

pytorch的cpp拓展

2020-06-09 12:11:16 1343 1

原创 关于GMRES的总结

关于GMRES的一些总结包括理论和程序

2020-05-25 17:44:14 11025 2

原创 givens变换化上hessenborg矩阵为上三角矩阵

givens旋转在GMRES上的应用

2020-05-25 17:16:14 2876 1

原创 Basic GMRES

一般的GMRES方法

2020-05-22 22:03:39 4592

原创 两种生成对称正定矩阵的方法

Pascal三角方法和U'VU

2020-05-22 16:43:08 10148 1

原创 等式约束凸二次规划问题的解法

日期:2020-5-8

2020-05-10 21:29:17 3326 4

原创 pytorch中backward的几个要点

torch.backward()

2020-04-28 15:55:33 849

原创 运用递归的方法求解对称三对角矩阵的特征值——matlab实现

运用递归的方法求解对称三对角矩阵的特征值

2020-04-16 15:04:47 4316 2

原创 Neuman边界条件

Neuman边界条件对这样一个问题我的目标:我们有如下的效果:一阶近似t=0.06t = 0.1t = 0.9二阶近似1.t=0.62.t=0.13.t = 0.9误差的差别改变步长取t=0.91.M=20,dt = 0.001一阶二阶差别2.M=40,dt=0.0003,t=0.6第二道题t = 0.9...

2020-04-13 23:55:53 8320

原创 pytorch中关于复数的处理

pytorch中的关于复数的表达及其处理,并设计高通滤波

2020-04-08 10:10:58 12532

原创 matlab与python图像之间的数据搬运问题

matlab与python图像之间的数据搬运问题我们用matlab作图并保存图像为png,然后用python读取出来,会发现数据的大小发生了改变这时候1.我们可以在matlab中把转变为图像的矩阵保存为.mat2.用python读取.mat文件(函数存在于scipy中)这样问们就获得了与matlab一样的矩阵,显示的时候可以显示出相同的图片。最重要的是数据本身没有改变。代码如下:#i...

2020-04-07 17:28:43 414

原创 微分方程组求解(Euler法、RK法)

微分方程组求解(Euler法、RK法、线性多步法)

2020-04-04 21:40:59 5410

原创 中点方法、改进欧拉及Heun方法

解一道题格式的表达代码示例:%% 用三种方法求解公式h = 0.01; t = 0:h:1;N = length(t);y = zeros(1,N);y(1) = 1;df = @(t,y)(-y+t+1);%改进Euler法(迭代方法)w = y(1);for k = 1:N-1 wk = w + h*df(t(k), w); z1 = wk; z...

2020-03-30 11:27:00 6390

原创 二阶泰勒展开法求解微分方程

日期:3030-3-29二阶泰勒展开法求解微分方程1.

2020-03-29 20:55:29 3983

原创 显式欧拉格式——matlab实现

日期:2020-3-16显式欧拉格式格式为:我运用显式欧拉格式计算三道例题,然后展现其效果,最后将我的代码展示出来1.2.3.呈现出这样的曲线的原因是该函数由一开始的极其陡峭变为后来的非常平稳,其中斜率的变化率非常之大。而欧拉格式走的是直线,在一个步长内斜率是不发生变化的,也就是当这个折线走完一个步长,与真实的曲线之间在数值上已经有很大的差别了。后面随着曲线的陡峭程...

2020-03-16 15:13:56 8117

原创 用几种欧拉格式算一道题——基于matlab

作者:老李日期:2020-3-15用几种欧拉格式算一道题在这篇报告中,我先使用了显示欧拉格式来求解微分方程,隐式欧拉格式中,我使用了改进的欧拉格式和修正的欧拉格式来进行求解。这道题是这样的:我们要做的是,通过已有的信息,求出 y(3)。我使用matlab来实现这个事情首先设立参数和函数句柄代码如下%% function and parament% we are going t...

2020-03-15 15:13:56 3779

原创 ct图像中的金属伪影

金属伪影的基本认识

2020-03-02 16:24:00 11876 2

原创 如何用matlab查看自己的显卡设备

水文其实就一句:gpuDevice然后效果如下:matlab有自己的神经网络模块:matconvnet然后pytorch里面对10这个版本只有10.1。但是我装的是10.2,可能是因为这个原因,我gpu版本的pytorch始终没有安装成功...

2020-02-22 09:29:00 2148

原创 数据集的准备

2020-2-17数据集的准备载入数据集是一件很重要的事情首先,我将要训练的图片放在了一个文件夹里,然后我要做的事情是把这些文件导入到我的程序中一般会这样的文件的名称形如dataset.py在这样的文件里主要做两件事情:1.将(数据,标签)以列表的形式传入列表变量中2.把传递数据后的列表放在一个类里,这个类继承torch.utils.data.Dataset的方法我所做的事情都是为...

2020-02-17 22:33:39 495

原创 动手搭U-Net

动手搭建U-Net我终于开始搭建第一个网络了,我决定拿UNet下手论文中对该结构的描述是这样的:Network ArchitectureThe network architecture is illustrated in Figure 1. It consists of a contracting path (left side) and an expansive path (right...

2020-02-16 21:27:51 869

原创 正则化

日期:2020-2-15pytorch中用于归一化的函数为nn.BatchNorm2d文档如下:

2020-02-15 21:44:56 179

原创 手写CNN

日期:2020-2-14手写CNN这是我第一次学着写CNN,我将步骤和代码记录下,并写了许多注释1.调包import osimport torchimport torch.nn as nnimport torch.utils.data as Dataimport torchvisionimport matplotlib.pyplot as plt2.声明超参数# Hyper ...

2020-02-14 21:00:07 730

原创 pytorch中的minibatch

日期:2020-2-13pytorch中的torch.utils.data这个库可以非常好的对数据实现批处理主要用到2个函数import torchimport torch.utils.data as DataData.TensorDataset()#设置数据集,数据与标签相对应Data.DataLoader()#传入数据集,设置批处理大小,是否打乱数据,顺序,线程数举个例子#ba...

2020-02-13 16:43:45 2866 2

原创 分类和回归动态表示

日期:2020-2-12这是我学着第一次搭建的网络,并将学习的过程以动态的形式表现了出来我在代码中写了详尽的注释# regression asd classficiationimport torchfrom torch.autograd import Variableimport torch.nn.functional as F #引入神经网络模块import matplotlib....

2020-02-12 21:51:57 222

原创 pytorch中Variable的自动求导

日期:2020-2-11作者:老李pytorch中variable这个变量具有自动求导功能,只需要在参数列表中加入requires_grad = True举个例子我们如下定义:我们对矩阵进行简单的求导var=[1234] var = \begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} var=[13​24​]我们对该矩阵平方的均...

2020-02-11 21:43:14 462 4

原创 线性插值方法去伪影

日期:2020-2-9dudonet的结构如下第一步要通过线性插值实现去金属伪影代码如下:clear; clc; close all%% 建立原始图像并加入金属球Img0=phantom('Modified Shepp-Logan',401);vec = -200:200;[xx,yy] = meshgrid(vec,-vec);x1 = 30;y1 = 50;r1 = 80;...

2020-02-09 23:53:49 2149

原创 CUDA实现(A+B)/2

作者:老李日期:2020-2-8这是我的第一个CUDA程序,也是照着英伟达社区的讲座来写的,写这篇文章的目的是,我希望自己能够尽量把一个程序讲清楚。目标//target: (A+B)/2 = C我先把代码贴上:代码如下在这里插入代码片...

2020-02-08 21:15:43 609

原创 等弧度扇束的投影——matlab实现

时间:2020-1-31作者:老李等弧度扇束的投影的投影工作的步凑分为以下四个次序% 探测器形式为等弧度探测器1.找出扇束放射源的位置2.确定扇束的投射范围3.表达射线4.对射线加总起来代码如下%% generate fanbeam projection%% 2020-1-31% 探测器形式为等弧度探测器% 1.找出扇束放射源的位置% 2.确定扇束的投射范围% 3.表...

2020-01-31 22:34:35 1543 2

原创 matlab中自带的CT投影与反投影函数(扇束、平行束及其之间的转化)

2020-1-27作者:老李

2020-01-27 14:04:30 12218 8

原创 dudoNet中的反投影

2020-1-24dudoNet中的反投影表达式在dudo-Net这篇文章中,用到了对反投影的一个表达,这个表达令我眼前一亮。是这样的:其中,整个过程实际上就是滤波反投影(FBP)的过程。通过这个表达,我们可以得到CT的图像域与sino图域之间的梯度的一个关系:从此,我们得到了一个关于两个神经网络:sinogram domain enhancement net(SE-NET)...

2020-01-24 15:44:46 3573 1

单能谱CT重建-扇束直接重建.rar

写扇束FBP的比较少,我这一款您绝对满意

2021-08-01

matlab训练图像.rar

如何构造金属伪影,我做了一些处理,并批量构造了数据集。这里面有我的程序,读者可以通过修改路径自己进行构造。程序中的伪影效果是我对硬化效应的理解

2020-05-08

GMREScpp.rar

基于cpp的GMRES实现,里面头文件有cuda的部分,删去即可,数据文件为10*10的pascal矩阵,具有对称正定的性质,各位也可以自己生成数据,往后还可能会有cuda优化的版本

2019-09-22

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

提示
确定要删除当前文章?
取消 删除