- 博客(45)
- 资源 (2)
- 收藏
- 关注
RSS订阅原创 python AttributeError: ‘Colorbar‘ object has no attribute ‘set_clim‘
运行程序显示错误：AttributeError: ‘Colorbar’ object has no attribute ‘set_clim’解决:是matplotlib版本过高所致matplotlib3.3.0版本取消了Colorbar的set_clim方法，因此需要降低版本3.2.0即可解决问题首先pip uninstall matplotlib之后pip install matplotlib==3.2.0...
2021-06-02 10:45:52
2055
原创 tensorflow安装相关说明
首先，安装tensorflow前需要建立tensorflow环境执行以下命令创建python版本为3.6、名字为tensorflow的虚拟环境。conda create -n tensorflow python=3.6激活tensorflow环境activate tensorflow安装CPU版本:pip install --ingore-installed --upgrade tensorflow安装GPU版本:pip install --ingore-installed --upg
2020-10-17 12:32:43
209
1
原创 有限角CT的正则化图像重建
CT的基本问题是由获取到的投影数据重建出被扫描目标的断层图像。FBP中要求采集得到的投影数据是完备的。比如在扇束CT中，要求扫描角度范围大于180°+扇角(如何理解?)实际应用中在有限的旋转角度范围内扫描，得到的投影数据是不完备的。有时为了节省扫描时间或者减少辐射剂量，也在有限的旋转角度范围内扫描。后果:有限角CT重建采用FBP得到的重建图像会存在明显的滑坡伪影。解决:基于最优化理论的正则化方法可以在重建过程中引入被扫描目标的先验知识，因此可以用来处理不完备数据重建问题，而且能在一定程度上抑制噪声和
2020-09-24 18:22:43
1412
原创 浅谈 Nyquist–Shannon(奈奎斯特-香农)采样定理
Nyquist–Shannon sampling theoremNyquist–Shannon(奈奎斯特-香农)采样定理是数字信号处理领域中的一个定理，它是连接连续时间信号和离散时间信号的基本桥梁。定理内容 :如果一个系统以超过信号最高频率至少两倍的速率对模拟信号进行均匀采样，那么原始模拟信号就能从采样产生的离散值中完全恢复。为了防止由于混叠引起的信号被破坏，需要以奈奎斯特速率或更高的速率进行采样。如果不遵守这个基本要求，就无法消除混叠(混叠永久与原始频谱混合，两者无法区分)。下面是解释采样定理的时
2020-09-16 14:33:44
30919
3
原创 overleaf latex文档中插入matlab代码(有中文注释、关键字高亮)
在\begin{document}前插入\usepackage{listings} \usepackage{xcolor}\lstset{ language=Matlab, %代码语言使用的是matlab frame=shadowbox, %把代码用带有阴影的框圈起来 rulesepcolor=\color{red!20!green!20!blue!20},%代码块边框为淡青色...
2020-03-13 15:34:21
10357
2
原创 LearnOpencv
OpenCV版本：4.2.0如何安装OpenCV及一次性配置：参考博文注意添加属性表后，属性表应随着编译文件变动。新建编译文件后，把配置好的属性表添加进去如何编译自带Samples:参考博文...
2020-02-27 16:52:07
130
原创 NVIDIA CUDA入门学习
CUDA 安装电脑系统：win10-x64参考博文设备管理器-显示适配器：为NVIDIA GeForce MX150是支持CUDA的，所以在这里下载了CUDA安装包一路安装即可，按照博文一路安装下去。...
2020-02-12 09:42:44
230
原创 CNN LeNet-5模型
实在无法理解S2层中引入参数个数为2∗6=122*6=122∗6=12个;大概是每个特征图中，取局部像素平均值时，权值w和偏置b都是可训练的，所以每个特征图引入参数为2个把~...
2019-12-08 18:35:45
126
原创 求解CT问题无解情况
CT成像问题中经常遇到Ax=pA\boldsymbol{x}=\boldsymbol{p}Ax=p无解的情况。这时候，就需要寻找广义解，即找到一个x\boldsymbol{x}x有Ax−p→0A\boldsymbol{x}-\boldsymbol{p}\to 0Ax−p→0。即∣∣Ax−p∣∣||A\boldsymbol{x}-\boldsymbol{p}||∣∣Ax−p∣∣或∣∣Ax−p∣∣2→...
2019-12-07 19:51:15
192
原创 牛顿法求解非线性方程
基础知识[步骤1] 准备：选定初始近似值x0x_0x0，计算f0=f(x0),f0′=f′(x0)f_0=f(x_0),f_0'=f'(x_0)f0=f(x0),f0′=f′(x0)[步骤2] 迭代：按公式：x1=x0−f0/f0′ x_1=x_0-f_0/f_0'x1=x0−f0/f0′迭代一次，得到新的近似值x1x_1x1，计算f1=f(x1),f1′...
2019-12-02 20:27:28
1382
原创 下降搜索法
基础知识求解Ax=bAx=bAx=b，若A对称正定，问题可化为求解f(x)=min(12xTAx−bTx)f(x)=min(\frac{1}{2}x^TAx-b^Tx)f(x)=min(21xTAx−bTx)对于f(x)f(x)f(x)，其梯度为Ax−bAx-bAx−bpython代码import numpy as npimport matplotlib.pyplot as plt...
2019-11-30 20:03:31
281
原创 jacobi迭代法、高斯赛德尔迭代法python实现
jacobi迭代法代码#程序目的：实现jacobi迭代法import numpy as npimport matplotlib.pyplot as pltdef draw_fuction(): x=np.arange(0,8,0.01) y=np.arange(0,8,0.01) x,y=np.meshgrid(x,y) z=4*x+y-24 w=2*x+5*y-30 pl...
2019-11-28 13:50:40
1188
原创 CBP(卷积反投影)实现
CBP公式如下：f(x1,x2)=∫0πf^(r,ϕ)∗H(r)∣r=(x1,x2)⋅ϕdϕf(x_1,x_2)=\int_0^{\pi}\hat{f}(r,\phi)*H(r)|_{r=(x_1,x_2)\cdot{\phi}}d\phif(x1,x2)=∫0πf^(r,ϕ)∗H(r)∣r=(x1,x2)⋅ϕdϕ求f(x1,x2)f(x_1,x_2)f(x1,x2)就简单地...
2019-11-20 13:26:00
2660
原创 几种滤波器频域时域表示
因为做卷积反投影，所以需要推导滤波函数在时域的表达。R-L 滤波器HRL(w)={∣w∣∣w∣≤wm0∣w∣>wmH_{RL}(w)=\left\{\begin{array}{c}|w| && |w|\leq w_m\\0 && |w|>w_m\end{array}\right.HRL(w)={∣w∣0∣w∣≤wm∣w∣>w...
2019-11-14 14:59:32
5103
原创 python实现龙贝格算法
基础部分代码如下程序完全照书上描述写得：#龙贝格算法实现：import numpyimport math#外推法def Trap(f,a,b,Iold,k): #返回用k阶复化梯形公式估计a、b范围内函数的积分 #定义Iold的必要性可以看书上例子p108例子3 if k == 1: Inew=(f(a)+f(b))*(b-a)/2 ...
2019-11-14 09:55:09
3668
1
原创 使用python构建学生-老师-课程选课系统
# 实现学生选课系统import shelveimport random#读取数据库（excel表）内信息# 学生类# 属性：姓名、学号、电话、所选课程列表# 方法：查看：显示该学生所有课程信息；添加课程：将选好的课程添加到课程列表中class Student: def __init__(self,name,number=0,phone=0,classname=None):#设置N...
2019-11-11 14:08:56
2173
3
原创 数据结构迷宫问题(链式栈)
链式栈结构的实现#include"fstream"#include<string>#include<iostream>using namespace std;#define OK 1#define ERROR 0typedef int SElemType;typedef int Status;typedef struct SNode {//链式存储栈...
2019-11-08 09:40:00
1285
1
原创 计算扇形与直线相交长，模拟扇形正投影
交线长计算公式：代码如下：import mathimport structimport numpy as npimport matplotlib.pyplot as plt #plt用于显示图片import sympy as sy#solve方法只能求解线性方程组from scipy.optimize import root,fsolve#扇形相关程序（扇形参数：角度phi,半径...
2019-10-29 16:21:03
633
原创 基函数做最佳平方逼近
使用基函数1,x,x2,…,xn1,x,x^2,\dots,x^n1,x,x2,…,xn做最佳平方逼近原函数f(x)=1+x2f(x)=\sqrt{1+x^2}f(x)=1+x2代码如下：import numpy as npimport mathimport matplotlib.pyplot as plt #plt用于显示图片import sympyfrom scipy impo...
2019-10-28 14:37:13
2286
1
原创 数据结构第二次作业
题目试编写在带头结点的单链表中删除一个最大值结点的高效算法。代码#include <stdio.h>#include<malloc.h>#include<stdlib.h>#include<conio.h>#define elemType int#define status int#define ERROR 0#define OK...
2019-10-14 15:44:27
505
原创 Homework of AdvancedProgramming
文章目录第一次作业题目：读取一个二进制图像文件，并尝试进行图像处理（更改其像素数值）再写回文件代码结果第一次作业题目：读取一个二进制图像文件，并尝试进行图像处理（更改其像素数值）再写回文件代码代码如下：其中image.raw是已存在的二进制图像文件。其中有四种灰度值（像素数值）变换函数。import numpy as npimport structimport matplotlib....
2019-10-13 20:07:47
200
翻译 Michigan Image Reconstruction Toolbox (MIRT)工具包使用说明
文章目录MIRT 官方文档解读下载安装许可证MIRT 官方文档解读原网页翻译如下：密歇根图像重建工具箱(MIRT)是一组用Mathwork的Matlab语言编写的图像重建(和相关的图像问题)的开源算法。这个软件是由杰夫·费斯勒和他的团队在密歇根大学开发的。工具箱包括以下内容。迭代和非迭代算法的层析成像(PET, SPECT, x射线CT)。磁共振(MR)图像重建的方法，包括补偿非共振效...
2019-10-12 23:25:00
1603
2
原创 数据结构第一次作业
文章目录题目：data代码代码结果取行结果：取列结果：思路题目：设计一函数，该函数对K个M*N的灰度图像（矩阵）处理。功能是取每个图像的第x行（或列），生成一个新的图像。要求使用指针实现上述程序。data白蚁的CT扫描图像(RAW)代码#include <fstream>#include <iostream>using namespace std;in...
2019-10-12 22:36:52
687
2
原创 CT正投影
文章目录所用知识编写程序，实现特殊几何形状。如圆、椭圆、三角形、正方形等的正投影编写程序，实现离散化图像的正投影所用知识编写程序，实现特殊几何形状。如圆、椭圆、三角形、正方形等的正投影编写程序，实现离散化图像的正投影...
2019-10-09 08:37:16
2129
2
原创 生成灰度图像并colormap
编程实现二进制数据的读写自己生成一幅图像，显示（通过第三方软件显示，验证对错）改变“传输函数”，将图像变为“彩色”，并进一步了解其工作原理...
2019-09-29 13:51:02
1967
2
原创 Python模拟Lagrange插值结果并分析误差
模拟Lagrange插值结果题目：exp(x)在[0,4]的n阶多项式插值及误差思路主要思路分为两部分：如何计算n阶多项式插值？误差如何计算？首先第一个部分：对于lagrange插值函数的计算，首先我们观察公式：Ln(x)=∑k=0nyklk(x)L_n(x)=\sum\limits_{k=0}^n y_k l_k(x)Ln(x)=k=0∑nyklk(x)其中，lk(...
2019-09-29 09:20:17
2579
原创 数字图像学习相关
数字图像入门什么是数字图像？数字图像，又称数码图像或数位图像，是以二维图像用有限数字数值像素的表示。由数组或矩阵表示，其光照位置和强度均为离散的。数字图像是由模拟图像数字化得到的、以像素为基本元素的、可以用数字计算机或数字电路存储和处理的图像。数字图像的显示和储存像素(pixel)像素是图形单元(pictureelement)的简称，它是位图最小的完整单位。像素具有两种属性：一种是相对于...
2019-09-28 09:06:53
840
原创 使用python进行数值分析
数值分析第一次作业题目：已知exe^xex泰勒展开公式为：ex=1+x+x22!+⋯+xnn!e^x=1+x+\frac{x^2}{2!}+\dots+\frac{x^n}{n!}ex=1+x+2!x2+⋯+n!xn;现设计一python函数，计算使用泰勒展开后e2e^2e2,e−2e^{-2}e−2以下两种方法的绝对误差结果。解题思路：要计算绝对误差结果，需要解析值-准确值；使用不同...
2019-09-22 00:01:54
2198
原创 微信小程序相关
文章目录写在前面的话微信小程序云开发部署微信小程序云开发写在前面的话微信小程序+云开发微信小程序云开发微信小程序云开发提供几大基础能力支持：云函数：在云端运行的代码，开发者只需编写自身业务逻辑代码数据库：一个即可在小程序前端操作，也能在云函数中读写的JSON数据库存储：在小程序前端直接上传/下载云端文件，在云开发控制台可视化云调用：基于云函数免鉴权使用小程序开放接口的能力，包括服...
2019-07-16 11:29:59
147
翻译 LFToolbox0.4官方使用文档说明
1 变化和未来计划见ReadMe文件2 安装Toolbox解压Toolbox至一个正确的路径，运行程序 LFMatlabPathSetup来安装Matlab路径。每次重新开启Matlab时都要设定这一步骤，所以你可以考虑在startup.m添加语句，例如：run('~/MyMatlabCodee/LFToolbox0.4/LFMatlabPathSetup.m')3 一个快速指南3....
2019-04-11 22:20:17
2640
7
原创 win下sublime卸载完成后默认使用sublime的文件还是sublime图标
c盘下搜索sublime，user文件夹下有两个sublime—text文件夹，删掉就好了，亲测有用。
2019-04-11 20:14:44
347
转载 网址存档
漫谈计算机摄影学：https://zhuanlan.zhihu.com/p/24982662光场子孔径图像提取（C++实现：https://blog.csdn.net/weixin_38285131/article/details/79483094
2019-04-08 20:14:30
142
转载 项目依赖安装参考
VS2017配置OpenMeshhttps://blog.csdn.net/codeNeverDie/article/details/80621966VS2017配置Eigen3https://blog.csdn.net/Kerwines/article/details/82807596VS2017配置CGALhttps://blog.csdn.net/lyh24k/article/de...
2019-03-11 12:00:07
210
翻译 Lytro In Matlab
本文翻译自：http://www.timzaman.com/2012/12/imaging-lytro-in-matlab/Lytro In Matlab这篇文章告诉你如何使用Matlab的Lytro相机拍摄的图片。Lytro图像被称为“光场图像”或者动态图像，但不管它们告诉你什么;这只是一幅由数字传感器前的“微透镜阵列”拍摄的普通图像。听起来容易吗?它是很容易的。你可以对光场图片做什么...
2019-02-25 14:29:33
458
原创 有关光场的网站
光场研究相关网站）lytro数据集：lytro数据集：网址：https://www.irisa.fr/temics/demos/lightField/index.html内容：lytro一代数据集，可以使用LFToolbox编码
2019-02-22 09:47:58
420
原创 常见latex编译问题
常见latex编译问题）数学公式下标正对latex换行数学公式下标正对\max \limits_{f}latex换行在需换行位置使用\
2019-02-20 15:42:07
703
原创 markdown插入数学公式
想法来源于:https://www.jianshu.com/p/054484d0892a借助MathJax引擎在markdown写作首部添加代码：&lt;script type="text/javascript" src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=default"&gt;&lt;/s
2019-02-20 10:35:09
644
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人