- 博客(67)
- 资源 (3)
- 收藏
- 关注
RSS订阅原创 深入理解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
45425
6
原创 深入理解SSIM(两图像结构相似度指标)(附matlab代码)
SSIM的深入理解作者：老李日期：2022-1-18SSIMSSIM用于评价两张图像的相似程度。对于SSIM这个指标的设计思路，如下图所示。图像的均值表示图像的亮度。图像的方差表示图像的对比度。注意：要先将图像归一化（把像素的阈值调整为[0,1]），才可以通过计算，得出数值进行比较。当α=β=γ=1，C_3= 0.5C_2（常用），则SSIM表达式为：“对比”是怎么体现的关于，这里我觉得是大部分博文所没有涉及到的地方。该问题等价于这也是以下三个等式如何产生的。也就是说以下三个
2022-01-18 04:49:19
26957
20
原创 matplotlib 修改字体
matplotlib 修改字体2021-9-9.ttc文件首先你得有一个字体的文件**（.ttc文件）**该类型文件的描述百度百科是这么写的：TTC全称是TrueType Collection，是一个字体格式。TTC字体是TrueType字体集成文件(. TTC文件)，是在一单独文件结构中包含多种字体,以便更有效地共享轮廓数据,当多种字体共享同一笔画时,TTC技术可有效地减小字体文件的大小。TTC是几个TTF合成的字库，安装后字体列表中会看到两个以上的字体。两个字体中大部分字都一样时，可以将两
2021-09-09 11:00:01
702
原创 所遇到的matlab中的路径操作的一点总结
matlab中的路径操作最近在看matconvnet，我发现里面有频繁的路径操作。而对于一个比较大的程序而言，里面必然有很多的关于路径的操作，于是想写一篇分享总结一下。我们在使用matconvnet的时候，无论是test还是train，运行的一开始，会有一个run(fullfile(fileparts(mfilename('fullpath')), ... '..', '..', 'matlab', 'vl_setupnn.m')) ;(mfilename(‘fullpath’)：返回当先
2020-07-25 17:13:55
1010
1
原创 基本的GMRES基于cpp实现
基本的GMRES基于cpp实现理论部分和matlab版本可以查看我的其他文章：基本理论完整实现以下为代码，头文件里面有cuda的东西，
2020-07-22 11:44:18
1418
3
原创 cell与struct
matconvnetm描述某一层的基本用法是struct+cellstruct可以理解为matlab版的键值对(dictionary)cell可以理解为c里面的结构体，他的操作很随意，可以任意加东西进去
2020-07-06 18:10:04
599
原创 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
9295
原创 matlab与python图像之间的数据搬运问题
matlab与python图像之间的数据搬运问题我们用matlab作图并保存图像为png，然后用python读取出来，会发现数据的大小发生了改变这时候1.我们可以在matlab中把转变为图像的矩阵保存为.mat2.用python读取.mat文件(函数存在于scipy中)这样问们就获得了与matlab一样的矩阵，显示的时候可以显示出相同的图片。最重要的是数据本身没有改变。代码如下：#i...
2020-04-07 17:28:43
467
原创 中点方法、改进欧拉及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
7592
1
原创 显式欧拉格式——matlab实现
日期：2020-3-16显式欧拉格式格式为：我运用显式欧拉格式计算三道例题，然后展现其效果，最后将我的代码展示出来1.2.3.呈现出这样的曲线的原因是该函数由一开始的极其陡峭变为后来的非常平稳，其中斜率的变化率非常之大。而欧拉格式走的是直线，在一个步长内斜率是不发生变化的，也就是当这个折线走完一个步长，与真实的曲线之间在数值上已经有很大的差别了。后面随着曲线的陡峭程...
2020-03-16 15:13:56
8634
原创 用几种欧拉格式算一道题——基于matlab
作者：老李日期：2020-3-15用几种欧拉格式算一道题在这篇报告中，我先使用了显示欧拉格式来求解微分方程，隐式欧拉格式中，我使用了改进的欧拉格式和修正的欧拉格式来进行求解。这道题是这样的：我们要做的是，通过已有的信息，求出 y(3)。我使用matlab来实现这个事情首先设立参数和函数句柄代码如下%% function and parament% we are going t...
2020-03-15 15:13:56
4198
原创 如何用matlab查看自己的显卡设备
水文其实就一句：gpuDevice然后效果如下：matlab有自己的神经网络模块：matconvnet然后pytorch里面对10这个版本只有10.1。但是我装的是10.2，可能是因为这个原因，我gpu版本的pytorch始终没有安装成功...
2020-02-22 09:29:00
2369
原创 数据集的准备
2020-2-17数据集的准备载入数据集是一件很重要的事情首先，我将要训练的图片放在了一个文件夹里，然后我要做的事情是把这些文件导入到我的程序中一般会这样的文件的名称形如dataset.py在这样的文件里主要做两件事情：1.将（数据，标签）以列表的形式传入列表变量中2.把传递数据后的列表放在一个类里，这个类继承torch.utils.data.Dataset的方法我所做的事情都是为...
2020-02-17 22:33:39
589
原创 动手搭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
952
原创 手写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
812
原创 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
2935
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
280
原创 pytorch中Variable的自动求导
日期：2020-2-11作者：老李pytorch中variable这个变量具有自动求导功能，只需要在参数列表中加入requires_grad = True举个例子我们如下定义：我们对矩阵进行简单的求导var=[1234] var = \begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} var=[1324]我们对该矩阵平方的均...
2020-02-11 21:43:14
496
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
2529
原创 CUDA实现(A+B)/2
作者：老李日期：2020-2-8这是我的第一个CUDA程序，也是照着英伟达社区的讲座来写的，写这篇文章的目的是，我希望自己能够尽量把一个程序讲清楚。目标//target: (A+B)/2 = C我先把代码贴上：代码如下在这里插入代码片...
2020-02-08 21:15:43
681
原创 等弧度扇束的投影——matlab实现
时间：2020-1-31作者：老李等弧度扇束的投影的投影工作的步凑分为以下四个次序% 探测器形式为等弧度探测器1.找出扇束放射源的位置2.确定扇束的投射范围3.表达射线4.对射线加总起来代码如下%% generate fanbeam projection%% 2020-1-31% 探测器形式为等弧度探测器% 1.找出扇束放射源的位置% 2.确定扇束的投射范围% 3.表...
2020-01-31 22:34:35
1827
2
原创 dudoNet中的反投影
2020-1-24dudoNet中的反投影表达式在dudo-Net这篇文章中，用到了对反投影的一个表达，这个表达令我眼前一亮。是这样的：其中，整个过程实际上就是滤波反投影（FBP）的过程。通过这个表达，我们可以得到CT的图像域与sino图域之间的梯度的一个关系：从此，我们得到了一个关于两个神经网络：sinogram domain enhancement net（SE-NET）...
2020-01-24 15:44:46
3678
1
GMREScpp.rar
2019-09-22
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人