本文目录一览：
- 1、matlab中通过iradon函数转化后，怎么设置矩阵大小
- 2、iradon函数怎么修改旋转中心
- 3、医学图像处理：CT与CBCT
- 4、为啥用iradon函数会出现负值
- 5、matlab I=iradon(...) 得到的矩阵I的各个值是什么意思？是灰度还是吸收率？还是什么？
matlab中通过iradon函数转化后，怎么设置矩阵大小
当P在CD上时x+y=1,
当P在BC上时AP=(1-y)AB+yAC=2(1-y)AD+yAC,
∴x/2+y=1.
看图得x+y=1,x/2+y=1.
由x,y=0得0=x=2,0=y=1.
iradon函数怎么修改旋转中心
M文件致理解由系列语句组相独立运行体M脚本文件与M函数文件M脚本文件没参数传递功能M函数文件功能M函数文件格式严格规定必须function其格式：Function输变量＝函数名称（输入变量）语句；end；函数输入变量向量输入（输入）输入形式[x1,x2,x3]元向量形式；输同理M函数文件命名必须其函数名改变M脚本文件则完固定功能模块其运行产变量均全局变量区别于M函数局部变量并且没参数传递函数文件与命令文件主要区别于:1.函数文件般都要带参数,都要返结,命令文件没参数与返结2.函数文件变量局部变量,运行期间效,运行完毕自清除,命令文件变量全局变量,执行完毕仍保存内存3.函数文件要定义函数名,且保存该函数文件文件名必须函数名.mM函数文件变量自变量,变量用[]括起-iradon函数
医学图像处理：CT与CBCT
CT图像的本质是衰减系数成像。一束X射线（比如下图多束平行射线中其中一束）在穿透物体的过程中 ，其部分能量被物体吸收，导致射出的X射线产生衰减。根据比尔-朗博定律(Beer–Lambert law)，光束的衰减量A主要取决于三个因素：①物质对特定波长光的吸收率，即摩尔衰减系数；②光程；③衰减物质的浓度，即：-iradon函数
因此，该射线射出之后的能量强弱，实质上是该射线经过的物体，其当前截断面上，不同物质的分布的一个综合作用结果。假如我们能够通过某种方式计算得到物体某一个截断面上每条射线（假设密度足够大）上每个点的摩尔衰减系数 ，则该截断面的物质结构信息就可以被呈现出来。-iradon函数
实际上CT图像如果读取为grayscale，其数值（即CT值，单位为Hounsfield Unit）就是由物质的衰减系数 确定的(这个 不知道是不是就是上面的 )：
也就是说，在已知某个截面上入射X光线强度（由X射线源决定）以及射出强度（由接收器）的前提下，求解出射线经过的各个点的衰减系数，成为CT图像重建问题的关键。考虑到射线经过区域不同位置的衰减系数不同，X射线的衰减量应该是各个位置衰减系数的积分。对一组平行X射线束中的每一束上面的衰减量同时进行积分，就相当于对2D衰减系数分布图像 沿入射方向进行投影。换句话说，就是平行X射线束对某个截面的一次扫描，可以拿到该截面一个固定角度的一维投影结果。-iradon函数
对一幅图像进行多个角度的投影，将所有的投影向量concat起来，可以得到一个矩阵，如下图所示：（注意concat时对向量长度不足最长投影向量的部分进行补零处理）
沿不同角度投影之后得到的投影矩阵，可以通过反变换重建出原始图像。比如一个简单的例子：
上面说的投影与重建实际上就是Radon变换和iRadon的特殊情况。Radon变换的更一般的表述如下：
The Radon transform is the integral transform which takes a function f defined on the plane to a function Rf defined on the (two-dimensional) space of lines in the plane, whose value at a particular line is equal to the line integral of the function over that line. -iradon函数
简单地说就是Radon变换将一个（离散）2D平面上的函数 映射到线平面 上的函数 ，其中 。如果这样说不好理解的话，我们可以认为2D数字图像的像素值即函数 ，其在每个点 上的函数值即灰度值，或者是CT图像的CT值。实际上Matlab中对Radon变化的解释就是这样的：-iradon函数
我们可以看一个Radon变化的例子，比如对这张原图：
下面是这张图对应的Radon变化结果，其中纵坐标代表投影度数，横坐标其实就是对应位置的像素值，不过他用的单位是方块。假如X射线束从上往下为0°，旋转方向为逆时针，则45°时投影，中间点的像素值应该是0（对角线上两个黑块，像素值-1和两个白块，像素值+1）；同理135°时投影区域像素值不为-2（全黑）的长度最大。-iradon函数
中心切片定理：
CT重建过程示意：
参考文献：
成像是牙科患者临床评估的重要诊断辅助手段。 1960年代全景射线摄影技术的引入及其在1970年代和1980年代的广泛采用预示了牙科放射学的重大进步，其为临床医师提供了颌骨和颌面部结构单一，但非常详细的成像。然而，无论是单独使用还是组合使用该技术，口内和口外手术均受所有平面二维（2D）投影的共性问题所困扰：结构的放大，扭曲，叠加和歪曲。-iradon函数
另一方面，人们在三维（3D）射线摄影成像（例如，立体镜，可调孔径的CT）领域也进行了诸多尝试，虽然可以利用CT图像，但考虑到其成本、获取难度、以及较高的放射剂量，CT图像在牙科领域的应用受到了限制。-iradon函数
专门用于对颌面区域成像的锥形束CT（CBCT）的引入，标志着颌面区域成像无论是在数据采集还是在图像重建过程，从2D到3D的转变。CBCT引起了牙科各个领域的空前兴趣，因为其掀起了颌面成像的一场革命，促进了牙科诊断从2D图像过渡到3D图像，并真正实现了将成像的作用从诊断扩展到手术和外科手术的图像指导，这种图像指导仅仅需要一个第三方应用软件即可实现。-iradon函数
CBCT是一项最新技术。其成像过程的关键是一个两端分别固定在X射线源和检测器的旋转机架。该旋转机架将发散的金字塔形或圆锥形电离辐射源穿过感兴趣区域的中心，引导到另一侧的X射线检测器上。 X射线源和检测器围绕固定在感兴趣区域中心的旋转支点旋转。在旋转过程中，以完整或有时是局部的弧线获取了多个（从150到600多个）视场（FOV）的顺序平面投影图像。-iradon函数
此过程不同于传统的医学CT成像，使用呈螺旋形的扇形X射线束获取FOV的各个图像切片，然后堆叠这些切片以获得3D表示。CT图像的每个切片都需要单独的扫描和2D重建过程。相比之下，由于CBCT单次扫描时包含了整个FOV，因此仅需要机架的一个旋转序列即可获取用于图像重建的足够数据，见下图：-iradon函数
首先，CBCT最大的优势在于其快速成像：
但是其主要缺点（尤其是使用较大的FOV时）在于，由于检测板接收到大量散射后的放射结果，与噪声及 对比度分辨率(contrast resolution) 等相关的图像质量会受到较大限制。
为啥用iradon函数会出现负值
M文件大致可以理解为由一系列的语句组成的相对独立的一个运行体。分为M脚本文件与M函数文件。M脚本文件没有参数传递功能，但M函数文件有此功能。M函数文件的格式有严格规定，它必须以“function”开头，其格式如下：Function输出变量＝函数名称（输入变量）语句；end；当函数输入变量为向量输入（多输入）时，输入形式为[x1,x2,x3]的一元向量形式；输出同理。M函数文件的命名必须是其函数名，不可改变。M脚本文件则为完成一固定功能的模块，其运行时产生的变量均为全局变量，区别于M函数的局部变量，并且没有参数传递。函数文件与命令文件的主要区别在于:1.函数文件一般都要带参数,都要有返回结果,而命令文件没有参数与返回结果2.函数文件的变量是局部变量,运行期间有效,运行完毕就自动被清除,而命令文件的变量是全局变量,执行完毕后仍被保存在内存中3.函数文件要定义函数名,且保存该函数文件的文件名必须是函数名.m。M函数文件可以有多个因变量和多个自变量,当有多个因变量时用[]括起来。-iradon函数
matlab I=iradon(...) 得到的矩阵I的各个值是什么意思？是灰度还是吸收率？还是什么？
1.调用matlab的图像读取函数imread读图片，把图片放在work文件夹下I=imread('样品.jpg');2.变为灰度图G=rgb2gray(I);3.双击workspace中的G项，在弹出的arrayedit中就看到了灰度矩阵-iradon函数
