讲解手册:Compressive Unlimited Sampling for MRI and Hyperspectral Imaging
配套 slides:cus_literature_review.pdf(2026-10-29 版)。
怎么用这份文档
- 每一节对应一页 slide。节标题里的"第 N 页"就是 slide 右下角
N/31中的 N;章节标题页(MRI、Hyperspectral imaging、Comparison and ideas、Backup)没有编号,这里不单列。 - 每节结构固定:这页回答什么 → 页面内容逐项解释 → 怎么讲(建议口径) → 可能被问到的问题 → 生词表。
- 生词表格式:英文术语 | 中文 | 注释。全文所有术语在最后的"总术语表"里按字母顺序汇总,方便查。
- 文末有"数字出处表":slide 上每个数字来自哪篇论文的哪里。
- 讲话时间按 25 分钟讲 + 5 分钟讨论设计。
- 讲之前先读现场应对:时间不够时跳哪些页、过渡句、被提问和被质疑时怎么答、讨论怎么引导、意外情况怎么处理。
- 英文版本:所有需要你开口说的内容都配了英文,可以直接照读:
- 每页"怎么讲"后面有 English script(按建议时间写,约每分钟 130 词);
- 每个问题后面有 EN Q(对方可能怎么用英文问,帮你听懂)和 EN A(可以直接说的英文回答);
- 过渡句、跳过时的替代说法、讨论问题等单句后面有
EN:; - 现场应对 · 3 末尾有一张常用英文句型表(请对方重复、争取思考时间、承认错误、会后跟进等)。
- 动态范围一律用 bit 表示,括号里附 dB,换算见第 4 页的"dB 与 bit 的关系"。注意:slides 上 PSNR 的 dB(图像质量)是另一回事,不要和动态范围的 dB 混淆。
目录与时间分配
现场应对
讲之前先读这一节。它覆盖时间、提问、质疑、讨论、意外情况;每页的具体问题在各页小节里。
1. 时间控制
页面分三级(按"去掉后结论是否还完整"):
| 级别 | 页码 | 说明 |
|---|---|---|
| 必讲 | 1, 2, 4, 6, 8, 11, 16, 17, 18, 19 | 结论、CUS 条件、动态范围证据、两篇 MRI 方法、Xu 的结论、对比表、想法 |
| 可压缩成一句 | 3, 7, 9, 10, 12, 15 | 原理细节和结果表;只读标题和结论框 |
| 可跳过 | 5, 13, 14 | Liu(对 CUS 无直接动机)、MST、Bian;跳过时说一句结论 |
- 落后 3 分钟:第 5、13、14 页各只读标题(约省 2.5 分钟)。
- 落后 6 分钟以上:再把第 3、7、9、10、12、15 页压成"标题 + 结论框"(约再省 4 分钟),保住 17–19 页和讨论时间。
- 只剩 10 分钟(15 分钟版):讲 1 → 2 → 4 → 6 → 8 → 11 → 16 → 17 → 18 → 19。
- 提前讲完:讨论想法 3 时翻 Backup 的 Kilcullen(PDF 第 33 页)████████████ ███████;讨论想法 4 时翻 Monroy(PDF 第 35 页)。
跳过时可以说的一句话:
- 第 5 页:"Liu 2021 是低成本低场 MRI,瓶颈是 SNR 不是 ADC 量程,所以对我们没有直接动机,细节在 Backup。" "Liu 2021 is a low-cost, low-field MRI. Its bottleneck is SNR, not the ADC range. So it gives us no direct motivation. The details are in the Backup."
- 第 13 页:"MST 是目前快照式高光谱重建的代表网络,比之前最好的方法高 2.55 dB,计算量只有 4.4%;对我们的意义是学到的光谱先验可以替代固定字典。" "MST is the representative network for snapshot hyperspectral reconstruction today. It beats the previous best method by 2.55 dB, with only 4.4% of the computation. For us, it shows that a learned spectral prior can replace a fixed dictionary."
- 第 14 页:"Bian 2024 把宽带滤光片做在芯片上,透光 75%,更容易饱和,但它是二维逐像素的,不在 CUS 模型里。" "Bian 2024 puts broadband filters on the chip. They pass 75% of the light, so the sensor saturates more easily. But it is 2D and works pixel by pixel, so it is not in the CUS model."
2. 过渡句(节与节之间)
- 第 2 → 3 页:"先看 MRI:它的测量方式本身就很像 CUS。" "First, MRI. The way it measures is already very much like CUS."
- 第 9 → 10 页:"MRI 讲完了。高光谱有两种完全不同的压缩方式,只有一种适合 CUS。" "That is all for MRI. Hyperspectral imaging has two very different ways to compress. Only one of them fits CUS."
- 第 16 → 17 页:"把七篇放在一起,用 CUS 的三个条件对比一下。" "Now let's put all seven papers together and compare them using the three CUS conditions."
- 第 17 → 18 页:"基于这张表,我想提四个想法,请大家一起讨论哪个值得做。" "Based on this table, I want to propose four ideas. Let's discuss together which one is worth doing."
- 第 19 页之后(进入讨论):"这四个里面,大家觉得哪个最值得先做?或者有没有我漏掉的方向?" "Of these four, which one do you think we should do first? Or is there a direction I missed?"
3. 被提问时的通用做法
- 先确认问题:没听清就复述一遍再答("您是问 X 吗?")。 "Do you mean X?"
- 先给结论,再给依据:一句话答案 → 出处(哪篇论文、哪页)→ 如有需要翻 Backup。
- 分清三种来源,说的时候明确是哪一种:
- 论文结果:"Gabr 2009 实测……" "Gabr 2009 measured that ..."
- 我们的推断:"这是我们的推断,论文没说……" "This is our own inference. The paper does not say this ..."
- 我们自己的初步模拟(不在 slides 上):"我们在体模上粗略试过……只是初步结果。" "We tried this roughly on a phantom ... It is only a preliminary result."
- 不知道时,不要编。可以说:
- "这一点我调研里没覆盖到,会后查一下再回复。"(英文:"I didn't cover that in this survey; I'll check and follow up.")
- "论文没报告这个数,我不确定。"("The paper doesn't report that number, so I'm not sure.")
- 问题太长或偏题:先简短回应,再说"这个我们放到最后讨论时间细聊",记下来。 "Let's talk about this in detail in the discussion time at the end."
- 需要细节时翻 Backup(Backup 页没有页码编号,下面是 PDF 页码):
| 被问到 | 翻到 PDF 第几页 |
|---|---|
| Liu 2021 硬件、协议、成本 | 25–26 |
| Yuan 2021 的 CASSI 模型、GAP/ADMM 公式、算法对比 | 27–28 |
| MST 网络结构、完整对比表、消融 | 29–30 |
| Bian 2024 材料、制作、应用 | 31–32 |
| 单像素:MHz 循环移位 pattern(Kilcullen) | 33 |
| 单像素:中红外单光子(Wang) | 34 |
| 二维 modulo 展开(Monroy) | 35 |
| 参考文献 | 36 |
常用英文句型(现场用)
| 场合 | 英文 |
|---|---|
| 请对方重复 | Sorry, could you repeat the question? |
| 确认理解 | Just to check: do you mean X? |
| 争取思考时间 | Good question. Let me think for a moment. |
| 先给结论 | The short answer is yes (or no). Let me explain why. |
| 引用论文 | According to Gabr 2009, ... / The paper reports that ... |
| 说明是我们的推断 | This is our own inference. The paper does not say this. |
| 说明是初步模拟 | We tried this in a quick simulation. It is only a preliminary result, and it is not on the slides. |
| 指向 Backup | I have a backup slide on this. Let me show it. |
| 推到讨论时间 | Can we come back to this in the discussion at the end? |
| 承认错误 | You are right, that is my mistake. I will check the paper and correct it. |
| 接受建议 | That is a good suggestion. Let me write it down. |
| 会后跟进 | I'm not sure. I will check and get back to you after the meeting. |
| 时间不够要往下走 | I am short of time, so I will skip the details. They are in the Backup. |
| 邀请讨论 | I would like to hear your opinions. Which idea do you think is most promising? |
4. 针对整体工作的质疑(跨页问题)
-
Q:这个调研全面吗?为什么是这 7 篇? A:不是穷尽调研。按"近期、引用较多、和 CUS 相关"挑选:MRI 3 篇、高光谱 4 篇,每篇都核对过题目、作者、年份和期刊。经典论文去掉了,只保留 Behin 2005 和 Gabr 2009 作为动态范围的证据。对比表底部写明了这一点。
- EN Q: "Is this survey complete? Why did you choose these 7 papers?"
- EN A: "It is not a complete survey. We chose papers that are recent, well cited, and related to CUS: 3 on MRI and 4 on hyperspectral imaging. We checked the title, authors, year and journal of each one. We left out the classic papers. We only kept Behin 2005 and Gabr 2009 as evidence for the dynamic range problem. The bottom of the comparison table says this."
-
Q:CUS 和已有的 unlimited sampling / modulo 压缩感知有什么区别? A:已有工作要么对全部样本做 modulo,要么在随机压缩测量上做 modulo,但随机测量的相邻值差别大、折叠跳变多。我们的两点是:(1) 设计平滑编码的测量矩阵,让跳变稀疏;(2) 两阶段 OMP,把跳变取整到格点 2λℤ 上再重估信号。论文中在 DCT 字典下,10⁻³ 相对误差阈值的恢复率是 71.5%,最强的 OMP 基线是 49.6%。
- EN Q: "How is CUS different from existing unlimited sampling or modulo compressed sensing?"
- EN A: "Existing work applies modulo either to all samples or to random compressive measurements. But random measurements change a lot from one to the next, so there are many folding jumps. We do two things. First, we design a measurement matrix with smooth codes, so the jumps are sparse. Second, we use a two-stage OMP: it rounds the jumps to the grid of 2 lambda times an integer, and then re-estimates the signal. In the paper, with a DCT dictionary and a relative error threshold of 10 to the minus 3, the recovery rate is 71.5%. The strongest OMP baseline gets 49.6%."
-
Q:modulo ADC 真的有硬件吗? A:有原型。unlimited sampling 方向已有 modulo ADC 的硬件原型,例如 Bhandari 等 TSP 2021(Unlimited Sampling From Theory to Practice,带原型 ADC)和 Ordentlich 等 JSTSP 2018(A Modulo-Based Architecture for Analog-to-Digital Conversion);成像方面有 Zhao 2015 的 modulo 相机像素。但还没有看到用在 MRI 接收机或高光谱探测器上的。
- EN Q: "Does modulo ADC hardware really exist?"
- EN A: "Yes, there are prototypes. In unlimited sampling there are modulo ADC hardware prototypes, for example Bhandari and others, TSP 2021, 'Unlimited Sampling From Theory to Practice', with a prototype ADC. Another is Ordentlich and others, JSTSP 2018, 'A Modulo-Based Architecture for Analog-to-Digital Conversion'. In imaging, Zhao 2015 has a modulo camera pixel. But we have not seen one used in an MRI receiver or a hyperspectral detector."
-
Q:既然有增益切换和更高位数的 ADC,为什么还要 modulo? A:增益切换需要事先知道信号在哪里大(或者多扫几次),而且要控制电路;更高位数的 ADC 功耗和成本更高。在线圈上数字化这类功耗受限的场景,目前只有约 13.5 bit,比临床需要的约 15–16 bit 少。modulo ADC 不需要事先知道量程,也不削顶,这是它的潜在优势。(第 4 页专题有详细数据。)
- EN Q: "We already have gain switching and higher-bit ADCs. Why do we need modulo?"
- EN A: "Gain switching needs to know in advance where the signal is large, or it needs extra scans. It also needs control circuits. A higher-bit ADC uses more power and costs more. When we digitize on the coil, power is limited, and today we only get about 13.5 bits. Clinical use needs about 15 to 16 bits. A modulo ADC does not need to know the range in advance, and it does not clip. This is its potential advantage. Slide 4 has the detailed numbers."
-
Q:有噪声时 CUS 还能用吗? A:如实说:论文里的实验都是无噪声的,噪声(包括量化噪声)列为未来工作。已知的是,pattern 太平滑时差分会放大噪声,这是选 σ_A 的折中之一。
- EN Q: "Does CUS still work with noise?"
- EN A: "To be honest, all experiments in the paper are noise-free. Noise, including quantization noise, is listed as future work. What we know is this: if the pattern is too smooth, taking differences amplifies the noise. This is one of the trade-offs when we choose sigma A."
-
Q:MRI 的 k-space 采样真的够平滑吗? A:调研的论文都没研究这个问题。我们在 Shepp–Logan 体模上做过粗略检查:把 ADC 量程设为峰值的 1% 时,中心读出线上约 31% 的相邻采样折叠圈数不同,8 倍读出过采样后约 12%。原因是边缘锐利,k-space 振荡快。这是我们自己的初步模拟、无噪声,不在 slides 上;所以想法 1 才把"██████████████████"列为待讨论问题。
- EN Q: "Is k-space sampling in MRI really smooth enough?"
- EN A: "None of the papers in this survey study this question. We did a rough check on the Shepp-Logan phantom. With the ADC range set to 1% of the peak, about 31% of neighbouring samples on the center readout line have different fold counts. With 8 times readout oversampling, it is about 12%. The reason is sharp edges, which make k-space oscillate fast. This is our own preliminary simulation, without noise, and it is not on the slides. That is why idea 1 lists ███████ ████████████ as an open question."
-
Q:为什么没讲 CT、雷达、频谱感知? A:早期我们调研过 A–F 六个方向(单像素、MRI、CT、雷达、频谱感知、生物信号),这次聚焦 MRI 和高光谱,是因为两者"按顺序测量 + 动态范围大"最明确。其余方向在另一份 slides(application survey)里。
- EN Q: "Why didn't you cover CT, radar, or spectrum sensing?"
- EN A: "Earlier we surveyed six directions, A to F: single-pixel, MRI, CT, radar, spectrum sensing, and biosignals. This time we focus on MRI and hyperspectral imaging, because these two most clearly have sequential measurements and a large dynamic range. The other directions are in a separate deck, the application survey."
-
Q:这些想法里哪个最可行? A:给出自己的倾向,但留给讨论:想法 3(██████ ███████ ██████)████████████████████,████████████████████████,██████████████████ ███ ████;想法 1、█ ██████ ███ ████████████████;想法 4 ██████████████████。这是我们的判断,不在 slides 上。
- EN Q: "Which of these ideas is the most feasible?"
- EN A: "I have my own preference, but I want to leave this to the discussion. Idea 3, ██████████ ███ ██████-█████ ████████, ██ ███ ███████ ██ ████ █████ ██ ██████████. ██ ████ ███████ ███ ███ █████ ██ ███ ███████████ ██████, ██ ██ ███ █████ ███ ████████ ███ ████. █████ █ ███ █ ████ ███ ██████████ ██ ██████████████ ████. Idea 4 ██ █ ███ ██ ████████ ███████. This is our own judgment, and it is not on the slides."
5. 讨论环节
- 目标:让导师和组员决定哪个想法值得做、先做什么,而不是逐条答辩。
- 可以抛出的问题:
- "四个想法里,哪个最值得先在仿真里试?" "Of the four ideas, which one should we try first in simulation?"
- "MRI 方向需要拿到原始 k-space 数据(例如 fastMRI 的原始数据),组里有没有现成的?" "For MRI, we need raw k-space data, for example the raw data from fastMRI. Does our group already have some?"
- "单像素方向,组里或合作方有没有 DMD 硬件可以做实验?" "For single-pixel imaging, does our group or a partner have DMD hardware for experiments?"
- 记录:把被选中的想法、被指出的问题、需要会后查的点记下来,会后整理。
6. 意外情况
- 投影或电脑出问题:手机或平板里存一份 PDF;讲稿以这份手册为准。
- 有人指出 slide 上的错误:先承认,记下来,会后核对原文再更正;不要当场争论。 "Thank you for pointing that out. Let me write it down. I will check the original paper after the meeting and correct it."
- 被问到 slides 上的一个词是什么意思:直接用手册里的生词表解释。会后考虑把这个词换成白话(skill 规则:被问到的词通常该从 slide 上改掉)。
- 有人质疑某个数字:翻文末的"数字出处表",说出论文和页码;如果是我们的推断或模拟,明确说出来。 "This number comes from [paper], page [N]. / This one is our own inference (or simulation), not from the paper."
第 1 页 · Summary
这页回答什么:整场汇报的结论和要讨论的想法,一页看完。
页面内容逐项解释
- In MRI the k-space center sets the ADC range. MRI 中 k-space 中心的信号最强,ADC 的量程必须照顾它。临床 3 T 最多需要约 15 bit,扫描仪的 ADC 加过采样能覆盖(Gabr 2009);小鼠全身 7 T 扫描需要 20 bit(Behin 2005);贴在线圈上的数字化器目前只有约 13.5 bit(Sporrer 2017)。
- MRI reconstruction tools could accommodate modulo data. 现有 MRI 重建方法在结构上能接入 modulo 数据:Chung & Ye 的扩散模型把"数据一致性"做成独立的线性一步;BJORK 学习轨迹时已有硬件惩罚项,████████████。注意是 "could",论文本身没有处理 modulo 数据。
- Every compressive spectral measurement sums many bands. 压缩光谱成像每次测量都要累加很多波段:Yuan 综述说 8-bit 相机重建后每个波段只剩约 4 bit;Bian 的宽带片上滤光片比窄带滤光片透光多。
- Single-pixel hyperspectral video is the closest match to CUS. Xu 2024 是单探测器、一串顺序测量,最接近 CUS;但它的 pattern 是伪随机的,相邻测量不平滑。
- Four ideas for discussion. 四个想法,最后两页展开。
底部小字:CUS、modulo ADC、k-space 的一句话定义;单像素成像论文放在 Backup。
怎么讲:只读每条粗体句,最后说"单像素成像的三篇放在 Backup,时间关系不展开"。
English script(可直接照读):
This page is the summary of the whole talk. First, in MRI, the k-space center sets the ADC range. Clinical 3 tesla needs at most about 15 bits, and the scanner ADC with oversampling can cover it. Mouse whole-body scans at 7 tesla need 20 bits. On-coil digitizers today have only about 13.5 bits. Second, MRI reconstruction tools could accommodate modulo data. I say 'could', because the papers do not handle modulo data. Third, every compressive spectral measurement sums many bands. Fourth, single-pixel hyperspectral video is the closest match to CUS, but its patterns are pseudo-random, not smooth. Finally, I have four ideas for discussion, on the last two pages. The three single-pixel imaging papers are in the Backup. Because of time, I will not go through them.
可能被问到的问题
- Q:为什么只选这 7 篇? A:按"近期 + 引用较多 + 和 CUS 相关"挑选,MRI 3 篇、高光谱 4 篇;经典论文(如 Lustig 2007 Sparse MRI、Wagadarikar 2008 CASSI)已经去掉;唯一保留引用的老论文是 Behin 2005,因为它是 MRI 动态范围问题的直接证据。不是穷尽调研(结论页对比表底部有写)。
- EN Q: "Why did you choose these 7 papers?"
- EN A: "I picked papers that are recent, well cited, and related to CUS. There are 3 MRI papers and 4 hyperspectral papers. I removed classic papers, like Lustig 2007 Sparse MRI and Wagadarikar 2008 CASSI. The only old paper I kept is Behin 2005, because it is direct evidence of the MRI dynamic range problem. This is not a complete survey; the table on the conclusion page says so."
- Q:单像素成像为什么放 Backup? A:这次重点是 MRI 和高光谱;单像素成像的思路在高光谱的 Xu 2024 里已经体现。
- EN Q: "Why is single-pixel imaging only in the Backup?"
- EN A: "This talk focuses on MRI and hyperspectral imaging. The single-pixel idea already shows up in the hyperspectral paper by Xu 2024."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| k-space | k 空间 | 图像的二维傅里叶变换所在的空间;MRI 扫描仪直接采集的就是 k-space 数据 |
| ADC (analog-to-digital converter) | 模数转换器 | 把模拟电压变成数字;位数(bit depth)决定能分辨多少级 |
| modulo ADC | 取模 ADC | 输入超过 [−λ, λ) 时不削顶,而是"折回"范围内(类似时钟 12 点后回到 1 点) |
| bit depth | 位深 | ADC 每个采样用多少位表示;16 bit ≈ 96 dB 动态范围,20 bit ≈ 120 dB |
| gain stepping | 增益步进 | 扫描过程中按需要切换接收机放大倍数,避免大信号削顶、小信号淹没 |
| accommodate | 容纳、兼容 | 这里意思是"结构上能接入",比 "absorb(直接吸收)" 弱 |
| pseudo-random | 伪随机 | 由确定算法生成、看起来随机的序列 |
第 2 页 · CUS 原理
英文标题:CUS works when measurement patterns change smoothly in sequence
这页回答什么:CUS 是什么;一个应用要具备什么条件才适合 CUS。后面每篇论文都用这一页的三个条件来衡量。
页面内容逐项解释
- 公式 y = M_λ(AΦα):信号 x = Φα 在字典 Φ(如 DCT)下稀疏;A 是测量矩阵(每一行是一个测量 pattern);M_λ 是 modulo 运算。
- 四步:
- Measure:做压缩测量 Ax(测量次数少于信号长度)。
- Fold:modulo ADC 把每个测量值折进 [−λ, λ),永不饱和。
- Difference:对相邻测量做差分,折叠只表现为若干个大小为 2λk 的跳变(k 是整数)。
- Recover:两阶段 OMP,先联合估计信号和跳变,把跳变取整到格点 2λℤ 上,再重新估计信号。
- Key idea:测量 pattern(A 的行)一个接一个平滑变化,相邻测量值就接近,跳变少(稀疏),恢复容易。
- What an application must offer(三个条件):
- ADC 覆盖不了的动态范围;
- 测量按顺序进行,每次测量是很多项的和;
- 能控制 pattern 的顺序或平滑程度。
怎么讲:一句话版本——"CUS 用 modulo ADC 避免饱和,靠相邻测量平滑变化让折叠跳变稀疏,然后用两阶段 OMP 恢复"。最后强调三个条件,"后面每篇论文我都按这三点评价"。
English script(可直接照读):
CUS stands for Compressive Unlimited Sampling. The model is y equals M lambda of A Phi alpha. The signal is sparse in a dictionary Phi, like the DCT. A is the measurement matrix, and M lambda is the modulo operation. There are four steps. First, we take compressive measurements. Second, the modulo ADC folds each value into minus lambda to lambda, so it never saturates. Third, we take differences of neighboring measurements. Then folding only shows up as a few jumps of 2 lambda times an integer. Fourth, a two-stage OMP recovers the signal. In one sentence: CUS uses a modulo ADC to avoid saturation, smooth patterns make the jumps sparse, and two-stage OMP recovers the signal. An application needs three things: dynamic range the ADC cannot cover, sequential measurements that each sum many terms, and control over the order or smoothness of the patterns. I will judge every paper by these three points.
可能被问到的问题
- Q:为什么差分后折叠变成稀疏跳变? A:y = u − 2λk,差分后 Δy = Δu − 2λΔk。如果 k 大部分时间不变(相邻测量都在同一"圈"里),Δk 大多为 0,只在跨圈时非零,所以是稀疏的。
- EN Q: "Why do the folds become sparse jumps after differencing?"
- EN A: "We have y equals u minus 2 lambda k. After differencing, delta y equals delta u minus 2 lambda delta k. If k stays the same most of the time, delta k is mostly 0. It is non-zero only when we cross into a new fold. So the jumps are sparse."
- Q:pattern 太平滑有没有坏处? A:有。太平滑时差分后的测量矩阵各列太像(相干性 μ_aug 变大),也会放大噪声。论文里用随机游走步长 σ_A 在"跳变少"和"可区分性"之间折中。
- EN Q: "Is there a downside if the patterns are too smooth?"
- EN A: "Yes. If the patterns are too smooth, the columns of the differenced matrix look too similar. The coherence, mu aug, gets larger, and noise is also amplified. The paper uses the random-walk step size, sigma A, to balance fewer jumps against being able to tell columns apart."
- Q:为什么要取整到格点? A:真实跳变一定是 2λ 的整数倍;第一阶段的估计有误差,取整后若误差小于 λ 就能精确恢复跳变,第二阶段只需估计信号。
- EN Q: "Why do you round the jumps to the lattice?"
- EN A: "A true jump is always an integer multiple of 2 lambda. The first-stage estimate has some error. If that error is smaller than lambda, rounding gives the exact jumps. Then the second stage only needs to estimate the signal."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| compressive measurement | 压缩测量 | 测量次数少于信号维数,靠稀疏性恢复 |
| dictionary (Φ) | 字典 | 一组基向量,信号在其下系数稀疏;本文用 DCT |
| DCT (discrete cosine transform) | 离散余弦变换 | 常用的稀疏化变换,JPEG 也用它 |
| sparse | 稀疏 | 大部分元素为零 |
| fold / folding | 折叠 | modulo 把超出范围的值折回范围内 |
| jump | 跳变 | 相邻两个测量的折叠圈数不同,差分后出现 2λ 的整数倍 |
| lattice | 格点 | 这里指 2λℤ,即所有 2λ 的整数倍 |
| OMP (orthogonal matching pursuit) | 正交匹配追踪 | 贪婪稀疏恢复算法,每步选一个与残差最相关的原子 |
| two-stage OMP | 两阶段 OMP | 本文方法:联合 OMP → 跳变取整 → 只对信号再做一次 OMP |
| pattern | 测量图案 / 测量模式 | 测量矩阵的一行;在光学里就是投影的二值图案 |
第 3 页 · MRI 成像原理
英文标题:An MRI scanner measures the Fourier transform of the image along gradient-driven paths
这页回答什么:MRI 采集的数据是什么;接收链长什么样;有哪些采样轨迹。
页面内容逐项解释
- Signal model:
- 氢原子核(自旋)在主磁场 B₀ 中以与磁场强度成正比的频率旋转(进动,Larmor 频率)。
- 梯度线圈让磁场随位置线性变化(G(t)·r),不同位置转速不同,于是每个位置积累不同的相位——位置被编码进了相位。
- 线圈收到的是所有位置的总和:s(t) = ∫ρ(r) e^{−i2π k(t)·r} dr,正好是图像 ρ 的傅里叶变换在 k(t) 处的值。
- k(t) ∝ ∫₀ᵗ G(τ)dτ:k 是梯度对时间的累积("梯度面积")。控制梯度波形就能控制在 k-space 里沿什么路径采样。
- Each ADC sample is one Fourier coefficient:每个 ADC 采样就是图像的一个傅里叶系数(忽略线圈灵敏度)。
- Receiver chain(接收链):接收线圈 → 前置放大器 → 解调到基带 → I、Q 两路各一个 ADC。MRI 信号是复数,实部(I)和虚部(Q)分别数字化,所以 modulo ADC 会分别折叠两路。
- Sampling paths:
- Cartesian(笛卡尔):每次激发采一条直线,图像 = 逆 FFT。
- Non-Cartesian(radial 径向、spiral 螺旋):激发次数少,需要非均匀 FFT(NUFFT)重建。
- Acceleration(加速):少采一些线或 shot,再用先验(稀疏、网络)补全,这就是压缩感知 MRI。
- Conclusion:每个采样对整幅图求和、按连续路径顺序到达、欠采样时就是压缩测量——结构上符合 CUS;前提是相邻采样变化小。
怎么讲:重点讲公式下面那句"每个 ADC 采样是一个傅里叶系数,梯度控制采样路径",再指一下 I/Q 两个 ADC。轨迹图一带而过。
English script(可直接照读):
An MRI scanner measures the Fourier transform of the image. The main field makes the spins precess. Gradient coils make the field change with position, so position is coded in the phase. The coil adds up the signal from all positions. So each ADC sample is one Fourier coefficient of the image, and the gradients control the path in k-space. Here is the receiver chain: the coil, the preamplifier, demodulation to baseband, and then two ADCs, one for I and one for Q. A modulo ADC would fold each of them separately. The sampling paths can be Cartesian, radial, or spiral. So each sample sums the whole image, samples arrive in order along a path, and with undersampling they are compressive measurements. This fits CUS, if neighboring samples change slowly.
可能被问到的问题
- Q:为什么是傅里叶变换? A:梯度使相位随位置线性变化,相位因子 e^{−i2πk·r} 正是傅里叶核;线圈对所有位置积分,所以得到傅里叶系数。
- EN Q: "Why is it a Fourier transform?"
- EN A: "The gradient makes the phase change linearly with position. The phase factor, e to the minus i 2 pi k dot r, is exactly the Fourier kernel. The coil integrates over all positions. So we get Fourier coefficients."
- Q:多线圈怎么办? A:每个线圈有自己的接收链和 ADC,测到的是"图像 × 该线圈灵敏度"的傅里叶变换;并行成像(parallel imaging)利用多线圈冗余加速。
- EN Q: "What about multiple coils?"
- EN A: "Each coil has its own receiver chain and ADC. It measures the Fourier transform of the image times that coil's sensitivity. Parallel imaging uses this redundancy across coils to speed up the scan."
- Q:一次激发采一条线,"顺序"指什么? A:一条读出线(readout)内部的采样点是按时间顺序连续采的,CUS 的差分就沿这个顺序做。
- EN Q: "If one excitation gives one line, what does 'in sequence' mean here?"
- EN A: "Inside one readout line, the samples are taken one after another in time. CUS takes the differences along this order."
- Q:全采样的 MRI 也是压缩测量吗? A:不是。只有欠采样(加速)时测量数少于未知数,才是压缩的。
- EN Q: "Is fully sampled MRI also a compressive measurement?"
- EN A: "No. Only with undersampling, or acceleration, are there fewer measurements than unknowns. Only then is it compressive."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| spin | 自旋 | 这里指氢原子核的磁矩,可以想成小磁针 |
| precess / precession | 进动 | 磁矩绕磁场方向旋转,像陀螺摆头 |
| main field B₀ | 主磁场 | MRI 的静磁场,如 1.5 T、3 T;本页低场例子 0.055 T |
| Larmor frequency | 拉莫尔频率 | 进动频率,与磁场强度成正比(氢核约 42.58 MHz/T) |
| gradient coil | 梯度线圈 | 产生随位置线性变化的附加磁场,用来做空间编码 |
| phase | 相位 | 旋转转过的角度;不同位置相位不同,就能区分位置 |
| receive coil | 接收线圈 | 感应自旋信号的天线 |
| preamplifier | 前置放大器 | 第一级低噪声放大 |
| demodulate to baseband | 解调到基带 | 把 MHz 级的射频信号搬到零频附近,便于采样 |
| I / Q (in-phase / quadrature) | 同相 / 正交分量 | 复数信号的实部和虚部,各用一个 ADC |
| trajectory | 采样轨迹 | k(t) 在 k-space 中走过的路径 |
| readout (line) | 读出(线) | 一次数据采集窗口内采到的一串点,Cartesian 中是一条直线 |
| excitation | 激发 | 用射频脉冲把磁化翻倒,之后才能读出信号;一次激发后读一条线或一条 shot |
| shot | 一次激发对应的轨迹段 | 非笛卡尔采样里一条螺旋臂或一条辐条 |
| Cartesian / radial / spiral | 笛卡尔 / 径向 / 螺旋 | 三种常见轨迹 |
| FFT / NUFFT | 快速傅里叶变换 / 非均匀快速傅里叶变换 | 后者用于非网格点上的采样 |
| prior | 先验 | 对图像的假设(稀疏、平滑、网络学到的分布),用来补全欠采样 |
第 4 页 · MRI 的动态范围问题
英文标题:The k-space center sets the ADC range; clinical scanners cope, on-coil digitizers fall short
这页回答什么:为什么 MRI 需要很多 bit 的 ADC;临床、临床前、线圈上数字化三种情况各需要多少、现在够不够。
页面内容逐项解释
- Why the range is large(左栏):
- k = 0 处的采样 s(0) = ∫ρ(r)dr 是整幅图像的和,通常远大于其他采样。
- 幅值随 |k| 快速下降;外圈采样对应边缘和细节,接近热噪声水平。
- 固定增益时,中心不能削顶 → ADC 满量程由峰值决定 → 最低有效位(LSB)随之确定;同时量化噪声必须低于边缘处的热噪声,否则细节被量化噪声淹没。两个要求一起决定需要多少 bit。
- 示意图:横轴 k_x,纵轴 log|s|;红色虚线是 ADC 满量程,底部阴影是热噪声。注明了是示意图,不是数据。
- How many bits are needed(右栏,标题旁注明 1 bit ≈ 6 dB):
- 临床 3 T(人体,表面线圈和阵列,含 3D):信号低于约 15 bit(88 dB)。扫描仪用 14-bit ADC 加过采样,有效约 17 bit,并在 k-space 中心附近额外衰减(Gabr 2009)。→ 够用。
- 临床前:小鼠全身 7 T 扫描需要 20 bit;在 16-bit 系统上,扫描中切换接收增益让图像 SNR 提高 8 倍(Behin 2005)。→ 不够,靠增益切换。
- 线圈上数字化(通道越来越多,用光纤代替同轴电缆):目前约 13.5 bit(81 dB)(Sporrer 2017、Fan 2021),低于临床需要。→ 不够。
- 后面三篇近期 MRI 论文都没讨论 ADC 位深。
- Conclusion:临床 3 T 的控制台 ADC 靠过采样和增益控制能应付;短板在高场临床前扫描和低功耗的线圈上数字化。相邻 k-space 采样变化多快(决定跳变数量)在所调研的论文里没有研究。
怎么讲:"k-space 中心是整幅图的和,最大;边缘是细节,接近噪声,所以 ADC 要很多 bit。临床 3 T 实测约 15 bit,现在的扫描仪有效 17 bit,够用;小鼠全身 7 T 要 20 bit,不够;现在的趋势是把 ADC 放到线圈上,这类低功耗方案只有约 13.5 bit,也不够。后者是 modulo ADC 最可能有用的地方。"
English script(可直接照读):
The k-space center is the sum of the whole image, so it is the largest sample. The edge of k-space holds the details, and it is close to the noise level. So the ADC needs many bits. This plot is a schematic, not real data. For clinical 3 tesla, the measured signal is within about 15 bits, or 88 dB. Today's scanners have about 17 bits effective, so that is enough. A mouse whole-body scan at 7 tesla needs 20 bits, so that is not enough. The trend now is to put the ADC on the coil. These low-power designs have only about 13.5 bits, or 81 dB. That is also not enough. This last case is where a modulo ADC is most likely to help.
dB 与 bit 的关系
动态范围 = 最大信号 / 最小可分辨信号(通常是噪声)。两种写法:
- dB:20·log₁₀(最大/最小)。射频、信号测量领域习惯用 dB(如 Gabr 的实测)。
- bit:n 位 ADC 能分辨 2ⁿ 级,能覆盖的比值是 2ⁿ。
- 换算:20·log₁₀2 ≈ 6.02,所以 1 bit ≈ 6 dB;bit = dB / 6.02。(严格地,理想 ADC 对满量程正弦的量化信噪比是 6.02n + 1.76 dB,口头换算用 6 dB/bit 即可。) "Dynamic range is the largest signal over the smallest signal we can resolve. Each extra bit doubles the range, and that is about 6 dB. So 13.5 bits is about 81 dB, 15 bits is about 88 to 90 dB, 16 bits is about 96 dB, and 20 bits is about 120 dB."
| 动态范围 | dB | bit | 对应什么 |
|---|---|---|---|
| 线圈上 Δ-Σ 接收机 | 81 | ≈ 13.5 | Sporrer 2017、Fan 2021 |
| 临床 3 T 实测上限(人体) | 88 | ≈ 14.6,记作"约 15" | Gabr 2009 |
| 理想 16-bit ADC | 96(含 1.76 dB 为 98) | 16 | |
| 3D 体模最高值 | 96.2 | ≈ 16 | Gabr 2009 |
| 14-bit ADC + 过采样(500 kHz 带宽) | 103 | ≈ 17 | Gabr 2009 的扫描仪 |
| 小鼠全身 7 T | ≈ 120 | 20 | Behin 2005 |
可能被问到的问题
- Q:Behin 是 2005 年的,现在的 ADC 是不是已经够了? A:临床 3 T 基本够用;但有三类情况仍然紧张。 Gabr 等 2009 年在 3 T 人体上实测,表面线圈和阵列的信号低于约 15 bit(88 dB,含 3D),扫描仪用 14-bit、80 MHz 的 ADC 加过采样,有效约 17 bit,再加上 k-space 中心附近的额外衰减,就能覆盖。仍然紧张的是:(1) 高场临床前、整体线圈(Behin 的 20 bit);(2) 放进磁体、贴在线圈上的低功耗数字化,目前约 13.5 bit;(3) 高场直接射频采样时,时钟抖动会吃掉有效位数。详见下面的专题:现在的 MRI ADC 够不够用。
- EN Q: "Behin is from 2005. Aren't today's ADCs already good enough?"
- EN A: "For clinical 3 tesla, mostly yes. Gabr and colleagues measured in 2009 at 3 tesla on people. Surface coil and array signals stayed below about 15 bits, or 88 dB, including 3D. The scanner used a 14-bit, 80 megahertz ADC with oversampling, so about 17 bits effective, plus extra attenuation near the k-space center. Three cases are still tight. First, high-field preclinical whole-body coils, where Behin needs 20 bits. Second, low-power digitizers on the coil, inside the magnet, with about 13.5 bits today. Third, direct RF sampling at high field, where clock jitter eats effective bits."
- Q:20 bit 和 16 bit 差多少? A:差 4 bit,即幅度比差 2⁴ = 16 倍(约 24 dB)。
- EN Q: "How big is the difference between 20 bits and 16 bits?"
- EN A: "It is 4 bits. That is a factor of 2 to the 4, which is 16 times in amplitude, or about 24 dB."
- Q:gain stepping 怎么做? A:k-space 中心附近用低增益避免削顶,外圈用高增益提升弱信号,重建时再按增益比例校正。
- EN Q: "How does gain stepping work?"
- EN A: "Near the k-space center, we use low gain so the signal does not clip. In the outer part, we use high gain to boost the weak signal. In reconstruction, we correct each part by its gain ratio."
- Q:为什么说"折叠跳变多少没被研究"? A:动态范围讨论的是幅度跨度;CUS 关心的是相邻采样之间变化多快,这一点所调研的论文都没涉及。(如有人追问,可口头补充:我们在 Shepp–Logan 体模上做过粗略检查,中心读出线上相邻采样折叠圈数不同的比例约 31%,8 倍读出过采样后约 12%。这是我们自己的初步模拟,不在 slides 上。)
- EN Q: "Why do you say the number of folding jumps has not been studied?"
- EN A: "Dynamic range is about the span of amplitudes. CUS cares about how fast neighboring samples change. None of the papers I reviewed looks at this. If needed: we did a rough check on a Shepp-Logan phantom. On the center readout line, about 31 percent of neighboring samples are in different folds. With 8 times readout oversampling, it is about 12 percent. This is our own preliminary simulation. It is not on the slides."
专题:现在的 MRI ADC 够不够用
(对应第 4 页的追问;调研于 2026-10-29。每条都注明出处,我们的推断单独标出。动态范围用 bit 表示、括号附 dB,换算见上面的 dB 与 bit 的关系。)
一句话结论:临床 3 T 系统在控制台用"14–16 bit 高速 ADC + 过采样抽取 + 按 k-space 位置调增益",有效约 17 bit,能覆盖实测信号(约 15 bit 以内);位深仍然不够的,是高场/临床前整体线圈(约 20 bit)、贴在线圈上的低功耗数字化(约 13.5 bit)、以及高频直接采样时的时钟抖动。第二类正是 modulo ADC 可能有用的地方。
1. 实测:临床 3 T 信号需要多少 bit
- Gabr et al., J. Magn. Reson. 2009(Johns Hopkins,Philips Achieva 3 T,人体 + 体模)[1]:
- 典型 3 T 人体表面线圈信号:低于约 15 bit(88 dB),包括 3D 采集;心脏、脊柱阵列 低于约 12.5 bit(75 dB);头部阵列低于约 15 bit(88 dB)。
- 最高值:心脏阵列做 3D SSFP 体模,带宽 499 kHz,约 16 bit(96.2 dB)。
- 人体 2D 图像都低于约 13.6 bit(82 dB);人体 3D 最高约 14.6 bit(88 dB)。
- 早期文献估计 21–32 bit(125、150、193 dB),被这篇实测否定(偏高)。
- 1.5 T 比 3 T 约少 1 bit(6 dB;SNR 与 B₀ 成正比)。
- 建议设计时再留约 1.7 bit 余量(10 dB)。
- 他们的扫描仪:14-bit ADC、80 MHz 采样;在 500 kHz 带宽下过采样相当于多 3.16 bit,有效约 17 bit(103 dB)。接收链里还有"按 k-space 位置调整的放大"(中心附近额外衰减)。
- 所以:需要约 15 bit + 1.7 bit 余量 ≈ 16.3 bit,扫描仪有效 17 bit,够用但余量不大。
- 对比 Behin 2005(slides 第 4 页):小鼠全身 7 T 扫描需要 20 bit(约 120 dB),16-bit 谱仪上用增益切换让 SNR 提高 8 倍 [2]。这是临床前、高场、整体覆盖的极端情况,比 Gabr 的临床表面线圈多约 5 bit。
- 我们的推断:按 Gabr 的"动态范围随 SNR 即 B₀ 增长",7 T 比 3 T 约多 1.2 bit(20·log₁₀(7/3) ≈ 7 dB),即临床 7 T 阵列可能到约 16 bit、加余量后约 17.5 bit,略超过上面 17 bit 的能力。这是按论文的比例关系外推,没有找到 7 T 实测数据。
2. 现在的接收机怎么做
- 控制台 ADC + 过采样:ISMRM 2016 教育讲座(Greig Scott)[3]:
- 位数决定可支持的动态范围;理想 16-bit ADC 的量化信噪比 SNQR = 6.02n + 1.76 ≈ 98 dB。
- 但实际会因采样时钟抖动和采样保持电容的 kT/C 噪声下降:例如 100 MHz 输入、0.3 ps rms 抖动只剩 74.5 dB,相当于约 12 bit 有效。
- 现在常用过采样 + 抽取、在"交替 Nyquist 带"上采样来提升有效位数(每 4 倍过采样约多 1 bit);流水线(pipeline)ADC 在位数、带宽、抖动之间最平衡。
- 我们的推断:高场直接射频采样时(7 T 约 300 MHz)抖动的影响更大,这是"直接采样"路线的位深压力来源之一。
- 增益控制仍是标准做法:接收增益/衰减设置不当,k-space 中心会削顶 [4];Gabr 的扫描仪也在中心附近额外衰减 [1];Behin 的 gain stepping [2];还有自适应动态范围接收机的专利 [5]。Sensors 2021 的一篇用三种接收增益分别扫描、再在 k-space 拼接,SNR 提高 5–13%(3 次扫描,时间成本高)[6]。
- 更高位数的 ADC:超导 ADC(20 GHz 过采样,通信中可达 24 bit)在 4 T 小动物 MRI 上替换内置 16-bit ADC,显示出更大动态范围,但进一步提升被室温线圈的热噪声限制(IEEE Trans. Appl. Supercond. 2014)[7]。
3. 新趋势:在磁体里、贴着线圈数字化
线圈通道越来越多(16、32、64、128 通道),每通道一根同轴电缆带来干扰、发热和安全问题,所以趋势是在磁体内或线圈上直接数字化、用光纤传输:
| 方案 | 有效位数(动态范围) | 出处 |
|---|---|---|
| 片上 CMOS 双通道接收机(贴在线圈上,Δ-Σ 基带数字化),1.5–10.5 T | 最高约 13.6 bit(81.9 dB),每通道功耗 < 240 mW | Sporrer et al., IEEE TBioCAS 2017 [8] |
| Δ-Σ 数字光纤链路(现成器件原型) | 约 13.5 bit(81 dB),200 kHz 带宽;作者认为"典型 MRI 信号 > 80 dB(约 13.3 bit)",所以够用 | Fan et al., Rev. Sci. Instrum. 2021 [9] |
| 16 通道磁体内接收机(数字化 + 光纤) | 与磁体外的先进数字化器"基本相当";同平台 2016 年摘要用 14-bit、125 MSps ADC | Reber et al., IEEE TMI 2019 [10];ISMRM 2016 [11] |
| 模拟光纤链路 | 约 10–11 bit(61–66 dB,100–500 kHz),单独使用不够 | Gabr et al. 2009 [1] |
- 要点:这些低功耗、贴近线圈的方案约 13.5 bit,比 Gabr 测到的 3D 头部阵列约 15 bit 少约 1 bit,比加了余量的约 16.3 bit 少约 3 bit。Gabr 的结论就是"模拟链路单独不够,需要信号压缩或多个不同增益的并行通道"。
- 注意两种说法的差别:Fan 2021 说 81 dB "够用",依据是"典型 MRI 信号 > 80 dB";Gabr 2009 实测 3D 头部阵列到 88 dB 并建议再留 10 dB。所以"够不够"取决于对比的是典型信号还是最坏情况。
- 我们的推断(对 CUS 的意义):在线圈上用功耗受限的 ADC 时,位数不够、又不方便做复杂增益控制——这正是 modulo ADC 的潜在用武之地(折叠代替削顶,或代替多路增益并行)。目前没有找到在 MRI 接收机里用 modulo ADC 的工作。
4. 讲的时候怎么回答(30 秒版本)
"临床 3 T 实测信号约 15 bit 以内,现在的扫描仪 ADC 加过采样有效约 17 bit,再加增益控制基本够。但两种情况仍然紧张:高场临床前整体线圈,Behin 算出要 20 bit;还有现在的趋势——把 ADC 放到线圈上、用光纤传,低功耗方案只有约 13.5 bit,比需要的少 1 到 3 bit。后者是 modulo ADC 最可能有用的地方。"
English 30-second version:"For clinical 3 tesla, the measured signal is within about 15 bits, or 88 dB. Today's scanner ADC with oversampling gives about 17 bits effective, or 103 dB. With gain control, that is mostly enough. But two cases are still tight. First, high-field preclinical whole-body coils: Behin found they need 20 bits, or about 120 dB. Second, the current trend is to put the ADC on the coil and send data over fiber. These low-power designs have only about 13.5 bits, or 81 dB. That is 1 to 3 bits less than needed. This second case is where a modulo ADC is most likely to help."
参考来源
- R. E. Gabr, M. Schär, A. D. Edelstein, D. L. Kraitchman, P. A. Bottomley, W. A. Edelstein, "MRI dynamic range and its compatibility with signal transmission media," J. Magn. Reson. 198(2):137–145, 2009. PMC2873084
- R. Behin, J. Bishop, R. M. Henkelman, "Dynamic range requirements for MRI," Concepts Magn. Reson. B 26B:28–35, 2005. doi:10.1002/cmr.b.20042
- G. Scott, "RF Receivers: Signal Detection Chain, Digitization, System Noise Figures — from MRI Signal to Bits," ISMRM 2016 educational session. ISMRM 2016 #7108
- "Receiver gain," Questions and Answers in MRI. mriquestions.com
- "Adaptive dynamic range receiver for MRI," US Patent 7403010. Google Patents
- "MRI-Based Image Signal-to-Noise Ratio Enhancement with Different Receiving Gains in K-Space," Sensors 21(16):5296, 2021. doi:10.3390/s21165296
- "Superconductor Analog-to-Digital Converter for High-Resolution Magnetic Resonance Imaging," IEEE Trans. Appl. Supercond., 2014. doi:10.1109/TASC.2014.2361132
- B. Sporrer et al., "A Fully Integrated Dual-Channel On-Coil CMOS Receiver for Array Coils in 1.5–10.5 T MRI," IEEE TBioCAS, 2017. doi:10.1109/TBCAS.2017.2764443(81.9 dB 数值见 BUAA 论文页)
- M. Fan et al., "A high-dynamic-range digital RF-over-fiber link for MRI receive coils using delta–sigma modulation," Rev. Sci. Instrum., 2021. doi:10.1063/5.0047041 / arXiv:2105.13305
- J. Reber et al., "An In-Bore Receiver for Magnetic Resonance Imaging," IEEE TMI 39(4):997–1007, 2019. doi:10.1109/TMI.2019.2939090
- "Scalable, In-Bore Array Receiver Platform for MRI," ISMRM 2016 #2170. ISMRM
说明:Sensors 2021 那篇,网上有摘要说"固定增益需 19 bit、可变增益需 15 bit",但我们拿到的摘要里没有这两个数,所以没有采用。
专题生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| oversampling and decimation | 过采样与抽取 | 以远高于信号带宽的速率采样,再数字滤波降采样,等效增加位数(每 4 倍过采样约多 1 bit) |
| SNQR (signal-to-quantization-noise ratio) | 信号量化噪声比 | 理想 ADC:6.02n + 1.76 dB |
| clock jitter | 时钟抖动 | 采样时刻的随机误差,输入频率越高影响越大 |
| kT/C noise | kT/C 噪声 | 采样保持电容上的热噪声 |
| pipeline ADC | 流水线 ADC | 高速、中高位数的 ADC 结构 |
| direct RF sampling | 射频直接采样 | 不先下变频,直接对 MHz 级射频信号采样 |
| in-bore / on-coil receiver | 磁体内 / 线圈上接收机 | 在磁体里或线圈上直接数字化 |
| delta-sigma (Δ-Σ) modulator | Δ-Σ 调制器 | 高过采样、低位数量化 + 噪声整形的 ADC 结构 |
| RF-over-fiber | 光纤射频传输 | 用光纤传 MR 信号代替同轴电缆 |
| SFDR (spurious-free dynamic range) | 无杂散动态范围 | 链路线性度指标 |
| SSFP (steady-state free precession) | 稳态自由进动 | 一种快速成像序列,信号较强 |
| superconductor ADC | 超导 ADC | 低温工作、极高过采样率的 ADC |
第 4 页生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| dynamic range | 动态范围 | 最大可测信号与最小可分辨信号之比;本文用 bit 表示,1 bit ≈ 6 dB(见 dB 与 bit 的关系) |
| dB (decibel) | 分贝 | 比值的对数单位,20·log₁₀(幅度比);注意 PSNR 的 dB 是图像质量,不是动态范围 |
| noise floor | 噪声底 | 系统固有噪声水平,低于它的信号无法分辨 |
| thermal noise | 热噪声 | 电子热运动产生的噪声,MRI 中主要来自人体和线圈 |
| clip / clipping | 削顶 | 信号超过 ADC 满量程被截断 |
| full scale | 满量程 | ADC 能表示的最大输入 |
| least significant bit (LSB) | 最低有效位 | ADC 的最小量化台阶 |
| quantization noise | 量化噪声 | 连续值取整到台阶带来的误差 |
| SNR (signal-to-noise ratio) | 信噪比 | |
| eightfold | 八倍 | |
| schematic | 示意图 | 非真实数据 |
第 5 页 · M1 Liu 2021:低场无屏蔽 MRI
英文标题:A shielding-free 0.055 T scanner cancels interference with sensing coils and a network 论文:Liu et al., A low-cost and shielding-free ultra-low-field brain MRI scanner, Nature Communications 2021。
这页回答什么:低成本 MRI 是否会成为 modulo ADC 的应用场景。
页面内容逐项解释
- Problem:全球约每百万人 7 台 MRI;高场 MRI 需要液氦和屏蔽室,成本高。
- System:0.055 T 钐钴永磁体(约 750 kg),普通插座供电、功耗 < 1200 W,不需要磁屏蔽也不需要射频屏蔽。
- Interference cancellation(干扰消除,本文核心):
- 扫描仪周围放 10 个小感应线圈,专门记录电磁干扰(EMI)。
- 每个重复周期(TR)除了采 MR 信号,还加一个"只有干扰"的时间窗。
- 每次扫描单独训练一个 5 层 CNN:输入 10 个感应线圈的信号,输出接收线圈上的干扰。训练数据就是本次扫描的"只有干扰"窗口。
- 把预测的干扰从每条读出线里减掉,再做 FFT 重建。
- Results:噪声与同一扫描仪放在射频屏蔽室里相比只差 5%;4 种临床对比度共 30 分钟;25 名患者与 3 T 图像对照阅片;材料成本 < 2 万美元。
- 左图:论文 Fig. 2——感应线圈位置、每个 TR 的两个时间窗、CNN 训练与预测、逐线消除。
- Conclusion for CUS:论文没报告 ADC 位深和饱和;0.055 T 时瓶颈是 SNR 低,不是 ADC 量程。它说明了低成本接收机的需求,但对 modulo ADC 没有直接支撑。
怎么讲:这篇和 CUS 关系较弱,讲快一点:"它代表低成本 MRI 方向,核心是用感应线圈 + 网络代替屏蔽室;但低场的问题是信号太弱,不是太强,所以对 modulo 不是直接动机。"
English script(可直接照读):
This paper is less related to CUS, so I will be quick. It represents the low-cost MRI direction. The scanner runs at 0.055 tesla, with a permanent magnet and a normal wall socket. It needs no shielded room. Instead, 10 sensing coils record the interference, and a small CNN, trained again for each scan, predicts the interference and subtracts it. The noise is only 5% worse than in a shielded room. But the problem at low field is that the signal is too weak, not too strong. The paper does not report ADC bit depth or saturation. So it shows a need for cheap receivers, but it is not a direct motivation for modulo ADCs.
可能被问到的问题
- Q:既然不相关,为什么选它? A:低成本、便携 MRI 是近年热点(引用 266 次),我们原以为便宜接收机会用低位 ADC,读完发现瓶颈在 SNR,这本身是一个结论。
- EN Q: "If it's not really related, why did you include this paper?"
- EN A: "Low-cost and portable MRI is a hot topic in recent years. This paper has 266 citations. We first thought cheap receivers would use low-bit ADCs. After reading it, we found the bottleneck is SNR. That itself is a useful conclusion."
- Q:强干扰会不会让 ADC 饱和? A:论文没报告。理论上未屏蔽的强干扰可能抬高峰值输入,但没有数据,不能下结论。
- EN Q: "Could strong interference saturate the ADC?"
- EN A: "The paper does not report this. In theory, strong unshielded interference could raise the peak input. But there is no data, so we cannot draw a conclusion."
- Q:为什么每次扫描都重新训练? A:干扰随环境、时间、病人身体(人体像天线)变化,每次训练能适应当下的干扰。
- EN Q: "Why do they retrain the network for every scan?"
- EN A: "The interference changes with the environment, with time, and with the patient's body. The body acts like an antenna. Training for each scan lets the network adapt to the current interference."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| ultra-low-field (ULF) | 超低场 | 一般指 < 0.1 T |
| shielding / RF cage | 屏蔽 / 射频屏蔽室 | 阻挡外界电磁干扰的金属房间 |
| permanent magnet | 永磁体 | 不需要电流和液氦 |
| samarium–cobalt (SmCo) | 钐钴 | 一种稀土永磁材料,温度稳定性好 |
| EMI (electromagnetic interference) | 电磁干扰 | |
| sensing coil | 感应线圈 | 只用来测干扰,不测 MR 信号 |
| repetition (TR) | 重复周期 | 两次激发之间的时间 |
| CNN (convolutional neural network) | 卷积神经网络 | |
| contrast (T1, T2, FLAIR, diffusion) | 对比度(加权方式) | 不同成像序列突出不同组织特性 |
第 6 页 · M2 Chung & Ye 2022:扩散模型重建方法
英文标题:One prior learned on magnitude images reconstructs MRI from any sampling pattern 论文:Chung & Ye, Score-based diffusion models for accelerated MRI, Medical Image Analysis 2022。
这页回答什么:扩散模型(score-based model)怎么做欠采样 MRI 重建,为什么它不依赖采样模式。
页面内容逐项解释
- 顶部图:论文 Fig. 1,从纯噪声出发,逐步去噪得到图像;每一步插入测量数据 y 做数据一致性。
- Training (no k-space needed):
- 学习 score:s_θ(x, t) ≈ ∇ₓ log p_t(x),即"加了噪声水平 t 的图像分布"的对数密度的梯度,直观上是"往更像真实图像的方向走"。
- 训练方式:denoising score matching(给图像加噪声,训练网络预测怎么去噪)。噪声水平 σ(t) 从 0.01 到 378(VE-SDE)。
- 数据:约 2.5 万张 fastMRI 膝关节幅值图像,320×320。只用图像,不需要原始 k-space。
- Measurement model:y = P_Ω F S x,其中 S 是线圈灵敏度、F 是傅里叶变换、P_Ω 是采样 mask。
- Sampling: repeat 2000 times from noise(重建过程):
- Predictor:用 s_θ 做一步反向扩散(去噪一步)。
- Data consistency:x ← x + λA*(y − Ax),把估计拉回和测量一致。
- Corrector:一步 Langevin 动力学 x ← x + ε s_θ + √(2ε) z(沿 score 走一小步再加一点随机噪声)。
- 再做一次数据一致性。
- Key innovation:先验训练时从未见过采样 mask,所以同一个网络适用于任意 mask 和加速倍数;物理模型只通过第 2 步进入,而且对任何满足 ‖A*A‖ ≤ 1 的线性 A 都成立(论文 Prop. 1:此时这一步是非扩张的)。
怎么讲:"训练时只学图像长什么样,不管怎么采样;重建时每走一步去噪,就用测量数据校正一次。好处是换采样方式不用重新训练。对我们重要的是:数据一致性是独立的一步,██████████████████████████。"
English script(可直接照读):
During training, the model only learns what images look like. It does not care how we sample. It learns the score, which points toward more realistic images. Training uses about 25,000 fastMRI knee magnitude images, with no raw k-space. During reconstruction, we start from pure noise and repeat 2000 times. Each time, we take one denoising step, and then we correct the estimate with the measured data. The benefit is that we do not need to retrain when the sampling pattern changes. The same network works for any mask and any acceleration. What matters for us is this: data consistency is a separate step. █████, ██ █████ ███████ ██ ████ █ ███████ ████ ███████ ███████.
可能被问到的问题
- Q:score 是什么? A:概率密度对数的梯度 ∇log p(x)。沿 score 方向走,x 会变得更"像"训练数据。
- EN Q: "What exactly is the score?"
- EN A: "It is the gradient of the log probability density, gradient of log p of x. If we move x along the score, x looks more like the training data."
- Q:为什么要 2000 步?很慢吧? A:是,每张切片 10–20 分钟。论文 Fig. 10 显示 4 倍加速时约 50 步就够,约 500 步后不再提升;也可以先用 U-Net 结果热启动只跑 40 步(CCDF)。
- EN Q: "Why 2000 steps? Isn't that very slow?"
- EN A: "Yes, it takes 10 to 20 minutes per slice. But Fig. 10 in the paper shows that at 4 times acceleration, about 50 steps are enough. After about 500 steps, there is no more gain. We can also warm-start from a U-Net result and run only 40 steps. This is called CCDF."
- Q:只用幅值图像训练,怎么重建复数数据? A:把实部和虚部分别当作图像,用同一个先验处理(下一页讲)。
- EN Q: "If it's trained only on magnitude images, how does it reconstruct complex data?"
- EN A: "We treat the real part and the imaginary part as two separate images. We process both with the same prior. I will show this on the next slide."
- Q:‖A*A‖ ≤ 1 这个条件 MRI 满足吗? A:线圈灵敏度归一化(S*S = I)后满足,论文 Eq. 9。
- EN Q: "Does MRI satisfy the condition that the norm of A-star A is at most 1?"
- EN A: "Yes, if the coil sensitivities are normalized, so S-star S equals the identity. This is Equation 9 in the paper."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| score-based diffusion model | 基于分数的扩散模型 | 生成模型的一种,通过逐步去噪从噪声生成样本 |
| score | 分数 | ∇ₓ log p(x),对数概率密度的梯度 |
| prior | 先验 | 这里是网络学到的"MR 图像分布" |
| magnitude image | 幅值图像 | 复数图像取模,DICOM 中常见 |
| denoising score matching | 去噪分数匹配 | 训练 score 网络的方法:学习从加噪图像预测噪声方向 |
| VE-SDE (variance-exploding SDE) | 方差爆炸型随机微分方程 | 扩散过程的一种形式,噪声方差随时间增长 |
| reverse diffusion | 反向扩散 | 从噪声逐步去噪回到数据 |
| predictor / corrector | 预测器 / 校正器 | 数值求解反向 SDE 的两类步骤 |
| Langevin dynamics | 朗之万动力学 | 沿 score 方向走一步再加随机扰动的采样方法 |
| data consistency | 数据一致性 | 把估计投影/拉回到与测量相符 |
| sampling mask (P_Ω) | 采样掩模 | 指定采哪些 k-space 点 |
| coil sensitivity map (S) | 线圈灵敏度图 | 每个线圈对不同位置的响应强度 |
| non-expansive | 非扩张 | 映射不会把两点距离拉大,保证迭代稳定 |
| fastMRI | — | NYU 与 Meta 的公开 MRI 数据集 |
第 7 页 · M2 Chung & Ye 2022:扩展与结果
英文标题:Complex and multi-coil data reuse the same prior, with large gains on unseen masks
这页回答什么:方法能不能处理真实的复数、多线圈数据;效果和监督学习方法比怎么样。
页面内容逐项解释
-
左图:论文 Fig. 2,多线圈的两种做法。
-
Extensions:
- Complex data:predictor 和 corrector 分别作用在实部和虚部上,用同一个先验。
- Multi-coil:(a) 每个线圈单独重建,再用平方和开根号(root-sum-of-squares, SSOS)合并;(b) 每隔 m 步加一次耦合各线圈的数据一致性(SENSE 型)。论文对 1D mask 用 (b),2D mask 用 (a)。
-
表格(PSNR,dB):
数据 mask 与加速倍数 本文 监督基线 单线圈 变密度 Poisson,15× 30.46 22.84 单线圈 Gaussian 1D,4×(训练用的 mask) 33.32 33.01 多线圈 Gaussian 2D,15× 29.17 16.66 多线圈 Gaussian 1D,8×(训练用的 mask) 32.43 33.15 监督基线:单线圈 DuDoRNet,多线圈 E2E-VarNet,都在 Gaussian 1D mask 上训练。规律:换了 mask,本文大幅领先;用训练时的 mask,两者接近,多线圈甚至略输。
-
Conclusion for CUS:数据一致性是独立的线性步骤,████ ███ ████████████████████████;██████"████"██████████████:██████ ██ ██████████。代价:每张切片 10–20 分钟。
怎么讲:指着表格说"监督方法在没见过的 mask 上掉到 17–23 dB,本文保持在 29–30 dB;在训练用的 mask 上两者差不多"。然后读结论框。
English script(可直接照读):
Please look at the table. On masks they have never seen, the supervised methods drop to 17 to 23 dB. This method stays at 29 to 30 dB. On the mask used for training, the two are close. In the multi-coil 8 times case, the supervised method is even slightly better. For complex data, the same prior is applied to the real and imaginary parts. For CUS, the key point is that data consistency is a separate linear step. ██ ████████ ███ ████████████ ███ ████ ██ ████████. ██ ██ ███ █████████ █████████ ████ ███ █████████ █████, ███ ███ ███ ███████ █ █ ██ █████ ███ ███████ ████████. The cost is 10 to 20 minutes per slice.
可能被问到的问题
- Q:为什么监督方法在新 mask 上掉这么多? A:监督网络学的是"特定 mask 的伪影 → 干净图像"的映射,换 mask 后伪影模式变了。
- EN Q: "Why do the supervised methods drop so much on a new mask?"
- EN A: "A supervised network learns a mapping from the artifacts of one specific mask to a clean image. When the mask changes, the artifact pattern changes."
- Q:多线圈 8× 为什么输了? A:在训练分布内,专门训练的 E2E-VarNet 有优势;本文的 SSOS 合并也会带来轻微模糊(论文自述)。
- EN Q: "Why does this method lose in the multi-coil 8 times case?"
- EN A: "Inside the training distribution, E2E-VarNet has an advantage because it was trained for that case. Also, the SSOS combination in this method causes slight blurring. The paper says this itself."
- Q:不确定性怎么得到? A:并行采 8 个后验样本,算逐像素标准差(Backup 没放,论文 Fig. 8)。
- EN Q: "How do they get uncertainty estimates?"
- EN A: "They draw 8 posterior samples in parallel and compute the standard deviation at each pixel. This is Fig. 8 in the paper. It is not in our backup slides."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| multi-coil / parallel imaging | 多线圈 / 并行成像 | 用多个接收线圈同时采集 |
| root-sum-of-squares (SSOS) | 平方和开根号 | 多线圈图像合并的常用方法 |
| SENSE | 灵敏度编码 | 利用线圈灵敏度图联合重建的并行成像方法 |
| PSNR (peak signal-to-noise ratio) | 峰值信噪比 | 图像质量指标,越高越好 |
| variable-density Poisson | 变密度泊松采样 | 中心密、外围疏的随机 2D 采样 |
| Gaussian 1D / 2D mask | 高斯 1D / 2D 掩模 | 按高斯密度随机选线(1D)或选点(2D) |
| supervised baseline | 监督学习基线 | 用成对数据(欠采样→全采样)训练的网络 |
| DuDoRNet / E2E-VarNet | — | 两个代表性监督重建网络 |
| posterior sample | 后验样本 | 从 p(x |
第 8 页 · M3 BJORK 2022:联合优化方法
英文标题:BJORK learns the k-space trajectory and an unrolled reconstruction together 论文:Wang et al., B-spline parameterized joint optimization of reconstruction and k-space trajectories (BJORK) for accelerated 2D MRI, IEEE TMI 2022。
这页回答什么:怎样同时学"怎么采"(轨迹)和"怎么重建"(网络),并满足硬件限制。
页面内容逐项解释
- 顶部图:论文 Fig. 1,训练流程:参考图像 → 按当前轨迹 ω 做 NUFFT 前向采样加噪声 → 重建网络 → 与参考图像算损失 → 同时更新 θ(网络)和 ω(轨迹)。
- Joint objective:
min_{ω,θ} E ℓ(f_θ(ω; A(ω)x + ε), x) + μ₁ φ_g(|D₁ω|) + μ₂ φ_s(|D₂ω|)
- 第一项:重建误差(ℓ1 + ℓ2²)。
- φ:hinge 惩罚,只在超过限值时罚。D₁ω(一阶差分)对应梯度幅值,限值 26.7 mT/m;D₂ω(二阶差分)对应 slew rate,限值 150 T/m/s。
- B-splines:ω = Bc,用二次 B 样条参数化轨迹;变量少很多,且不会出现快速振荡。
- Reconstruction:展开 6 次迭代,每次是"共轭梯度(CG)数据一致性 + 共享权重的 CNN 去噪器"(MoDL 风格)。
- Gradients w.r.t. ω:对轨迹求梯度时,用精确的 NUFFT Jacobian 近似,而不是自动微分(论文 App. C:SSIM 0.957 vs 0.930)。
- Key innovation:采样密度从数据中学出来(不像早期的 SPARKLING 需要预先假设密度),在明确的硬件惩罚下,与将使用这些数据的网络一起优化。
怎么讲:"目标函数 = 重建误差 + 两个硬件惩罚。轨迹用 B 样条表示,天然平滑;重建是展开网络。对我们最重要的是这个'惩罚项框架'——██████████████████████████。"
English script(可直接照读):
BJORK learns the k-space trajectory and the reconstruction network together. The objective is the reconstruction error plus two hardware penalties. One penalty limits the gradient amplitude, and the other limits the slew rate. They only penalize values above the limit. The trajectory is represented with B-splines, so it is smooth by design and has far fewer variables. The reconstruction is an unrolled network with 6 iterations. The sampling density is learned from data, not assumed in advance. For us, the most important part is this penalty framework. █████, ██ █████ ███ ███ ████ ███████ ███ ███████ █████.
可能被问到的问题
- Q:slew rate 是什么? A:梯度变化率(T/m/s),受梯度放大器和人体神经刺激限制;对应轨迹的二阶导数。
- EN Q: "What is slew rate?"
- EN A: "It is how fast the gradient changes, in tesla per meter per second. It is limited by the gradient amplifier and by nerve stimulation in the body. It corresponds to the second derivative of the trajectory."
- Q:为什么用 B 样条? A:直接优化每个采样点,参数多、容易震荡、易卡在局部最优;B 样条把轨迹表示为少量控制点的平滑组合。
- EN Q: "Why use B-splines?"
- EN A: "If we optimize every sample point directly, there are many parameters. The trajectory can oscillate and get stuck in a local optimum. B-splines write the trajectory as a smooth combination of a few control points."
- Q:展开网络是什么? A:把迭代重建算法(这里是 CG + 去噪交替)固定迭代次数后展开成网络层,端到端训练。
- EN Q: "What is an unrolled network?"
- EN A: "We take an iterative reconstruction algorithm, here CG and denoising in turn. We fix the number of iterations and turn each iteration into network layers. Then we train it end to end."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| joint optimization | 联合优化 | 同时优化采样和重建 |
| B-spline | B 样条 | 分段多项式基函数,天然平滑 |
| hinge penalty | 铰链惩罚 | max( |
| gradient amplitude | 梯度幅值 | 梯度场强度,mT/m |
| slew rate | 梯度切换率 | 梯度变化速度,T/m/s |
| unrolled network | 展开网络 | 把迭代算法展开为固定层数的网络 |
| conjugate gradient (CG) | 共轭梯度法 | 解线性方程组的迭代法,这里用于数据一致性 |
| denoiser | 去噪器 | 这里是 CNN |
| Jacobian | 雅可比矩阵 | 输出对输入的偏导数矩阵 |
| automatic differentiation | 自动微分 | 深度学习框架自动求导 |
| SPARKLING | — | 一种早期的非笛卡尔轨迹设计方法,需预设采样密度 |
| MoDL | — | 一种展开式模型驱动深度学习重建框架 |
第 9 页 · M3 BJORK 2022:优化过程与结果
英文标题:Coarse-to-fine B-splines move trajectories far from radial and raise image quality
这页回答什么:学到的轨迹是什么样;比传统和已有学习方法好多少。
页面内容逐项解释
-
顶部图:论文 Fig. 5,从 32 条径向辐条出发,Decim 32 → 16 → 8 → 4 逐级细化,SSIM 0.957 → 0.963 → 0.965 → 0.967 → 0.968;最右是不用 B 样条、直接优化采样点的结果(0.959),几乎没离开径向。
-
Schedule:每个 B 样条系数先对应 32 个采样点,收敛后减半,共 4 级 × 3 个 epoch。
-
Prospective scan:3 T、32 通道线圈,32 个 shot 11 秒 vs 320 条辐条 1 分 40 秒(10 倍加速);网络重建和 CS 重建下伪影都更少。
-
表格(PSNR,dB):
设置 Standard SPARKLING BJORK 16 辐条,网络重建 32.7 33.9 34.3 32 辐条,网络重建 35.0 36.0 36.9 32 辐条,CS(ℓ1 小波) 33.9 35.7 36.3 8 臂螺旋,网络重建 40.9 41.7 41.9 fastMRI 脑部,1520 张测试切片,多线圈。
-
Conclusion for CUS:████████████"████████████"██████。注意:学到的波形在每个 shot 内部会抖动(见顶部图),可能比直辐条产生更多跳变。
怎么讲:"由粗到细的 B 样条让轨迹能离开初始形状;每种设置 BJORK 都最好,换成传统 CS 重建也有提升。提醒一点:学出来的轨迹有小抖动,对 CUS 未必有利。"
English script(可直接照读):
The B-splines go from coarse to fine. This lets the trajectory move away from its starting radial shape. Without B-splines, the trajectory almost stays radial. In every setting in the table, BJORK is the best. It also improves when we switch to traditional CS reconstruction. In a real scan at 3 tesla, they got 10 times acceleration with fewer artifacts. One warning: the learned trajectory has small wiggles inside each shot. This may cause more folding jumps than straight spokes, so it may not help CUS.
可能被问到的问题
- Q:Decim 是什么? A:decimation,每个 B 样条控制点覆盖多少个采样点;越大越粗、越平滑。
- EN Q: "What does Decim mean?"
- EN A: "It stands for decimation. It is how many sample points each B-spline control point covers. A larger value is coarser and smoother."
- Q:prospective 是什么意思? A:前瞻性实验,即在真实扫描仪上按学到的轨迹实际采集,而不是用已有数据回顾性模拟。
- EN Q: "What does prospective mean here?"
- EN A: "It means a prospective experiment. They actually scan on a real scanner with the learned trajectory. It is not a simulation on existing data."
- Q:比 SPARKLING 只高 0.2–0.9 dB,算大吗? A:在 MRI 里算明显(比 Standard 高 1.0–2.4 dB);且 CS 重建也提升,说明轨迹本身更好,不只是和网络配合。
- EN Q: "The gain over SPARKLING is only 0.2 to 0.9 dB. Is that really large?"
- EN A: "In MRI, this is a clear gain. Over the standard trajectory, the gain is 1.0 to 2.4 dB. Also, CS reconstruction improves too. This shows the trajectory itself is better, not just a good match with the network."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| coarse-to-fine | 由粗到细 | 多级优化策略 |
| decimation (Decim) | 抽取率 | 每个控制点对应的采样点数 |
| spoke | 辐条 | 径向轨迹中的一条直线 |
| SSIM (structural similarity) | 结构相似性 | 图像质量指标,最大 1 |
| prospective | 前瞻性 | 实际扫描验证 |
| in vivo | 活体 | 在活体(人)上的实验 |
| CS (compressed sensing) | 压缩感知 | 这里指 ℓ1 小波正则的传统重建 |
| ℓ1 wavelet | ℓ1 小波正则 | 以小波系数 ℓ1 范数为正则项 |
第 10 页 · 高光谱成像原理
英文标题:A hyperspectral cube must be squeezed into one 2D frame or one detector
这页回答什么:高光谱数据是什么;压缩高光谱成像有哪两类做法;modulo 在哪里起作用。
页面内容逐项解释
-
Data cube X(x, y, λ):每个像素有一条光谱,例如 256×256 像素 × 28 波段(450–650 nm),或 128×128 × 64 波段。左边是立方体示意图。
-
Scanning:逐行或逐波段扫描需要很多次曝光,运动场景下太慢。
-
CASSI / DMD 定义写在左栏:CASSI = coded aperture snapshot spectral imager(编码孔径快照光谱成像仪);DMD = digital micromirror device(数字微镜器件,一组可开关的小镜子)。
-
表格对比:
快照式(编码孔径、滤光阵列) 单像素式 探测器 二维阵列,一帧 一个光电二极管 编码 掩模 + 棱镜(CASSI),或每个像素不同滤光片 DMD 在空间和光谱上的 pattern 一次测量累加 一个像素上的所有波段 所有像素和波段 测量方式 并行(一帧) 顺序(一串) modulo 作用在 每个像素,二维 沿测量序列 两类都是压缩的:测量数少于体素数,靠先验(稀疏、全变分、网络)恢复。
-
Conclusion:单像素式符合 CUS 的测量模型(一串顺序测量),██████████ ███████ ██████████;██████████████████████ ██████,████████████████。
怎么讲:重点讲表格最后一行"modulo 作用在哪":"快照式是一帧图像里每个像素各自折叠,要做二维展开;单像素式是一串测量依次折叠,正好是 CUS 的模型。"
English script(可直接照读):
This page shows what a hyperspectral cube is. Each pixel has a full spectrum, for example 256 by 256 pixels and 28 bands. Scanning band by band needs many exposures, so it is too slow for moving scenes. Compressive systems come in two types. Snapshot systems, like CASSI or filter arrays, use a 2D detector and take one frame. Single-pixel systems use a DMD and one photodiode, and take a sequence of measurements. Both are compressive. They take fewer measurements than voxels and use priors to recover the cube. The key row is the last one: where does modulo act? In a snapshot system, each pixel in one frame folds on its own, so we need 2D unwrapping. In a single-pixel system, a sequence of measurements folds one after another. That is exactly the CUS model, ██ ████████████ ████████ ███ ██ ████ ███████.
可能被问到的问题
- Q:快照式为什么不能直接用 CUS? A:CUS 靠相邻测量差分;快照式一帧里的像素是空间上的邻居,不是时间序列,折叠后需要在二维空间上展开(类似相位展开),恢复问题不同。
- EN Q: "Why can't you apply CUS directly to snapshot systems?"
- EN A: "CUS uses differences between neighbouring measurements. In a snapshot frame, the pixels are neighbours in space, not a time sequence. So after folding, we must unwrap in 2D space, like phase unwrapping. That is a different recovery problem."
- Q:单像素为什么一次测量累加所有像素和波段? A:所有经 DMD 反射的光汇聚到一个探测器上,探测器只输出一个总强度。
- EN Q: "Why does one single-pixel measurement add up all pixels and all bands?"
- EN A: "All the light reflected by the DMD goes to one detector. The detector only gives one total intensity."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| hyperspectral image | 高光谱图像 | 每像素几十到上百个波段 |
| data cube | 数据立方体 | (x, y, λ) 三维数据 |
| band | 波段 | |
| voxel | 体素 | 三维数据中的一个元素 |
| snapshot | 快照式 | 一次曝光得到全部数据 |
| coded aperture | 编码孔径 | 带随机透光图案的掩模 |
| prism / disperser | 棱镜 / 色散元件 | 让不同波长偏折不同 |
| filter array | 滤光片阵列 | 每个像素上一个滤光片 |
| photodiode | 光电二极管 | 单个探测器 |
| DMD (digital micromirror device) | 数字微镜器件 | 每个小镜子可翻转开/关,用来投影二值 pattern |
| total variation (TV) | 全变分 | 惩罚相邻像素差的正则项,鼓励分段平滑 |
| unwrapping | 展开 | 从折叠值恢复真实值,类似相位解缠 |
第 11 页 · 高光谱的动态范围问题
英文标题:Every compressive spectral measurement sums many voxels, so the detector range runs out
这页回答什么:为什么压缩光谱成像会遇到探测器动态范围不够;论文里有什么证据。
页面内容逐项解释
- Snapshot model:Y = Σₙ Xₙ ⊙ Mₙ + E,每个像素累加 N_λ 个被编码的波段;向量形式 y = Φx + e,Φ = [D₁, …, D_{N_λ}],Dₙ = diag(vec Mₙ)。
- Single-pixel model:yᵢ = φ_{s,i}ᵀ X φ_{w,i},每个值累加整个立方体;φ_s 是空间 pattern,φ_w 是光谱 pattern,每次测量一对。
- What the papers report:
- 快照式(Yuan 2021):8-bit 相机重建后,"每个像素可能只有 4 bit 动态范围";8 → 12 bit 改善了结果。
- 单像素(Xu 2024):一半镜子打开时,每个测量约为场景总光强的一半,是"一个大常数 + 小变化"。
- 片上滤光片(Bian 2024):透光约 75%,马赛克滤光片 < 10%。
- Conclusion:累加让每次测量信号变大,而感兴趣的细节仍然很小——正是 CUS 设计针对的情形。四篇论文都没有给出饱和的解决办法,也没有用 modulo ADC。
怎么讲:"不管快照还是单像素,一次测量都是很多体素的和,信号大、细节小。Yuan 的综述直接说 8-bit 相机最后每个波段只剩 4 bit。"
English script(可直接照读):
Snapshot or single-pixel, one measurement is always a sum of many voxels. In the snapshot model, each pixel adds up all the coded bands. In the single-pixel model, each value adds up the whole cube. So the signal is large, but the details we care about are small. The papers show this. Yuan's review says that with an 8-bit camera, each pixel may have only 4 bits of dynamic range after reconstruction, and going from 8 to 12 bits improved the results. In Xu 2024, half the mirrors are on, so each measurement is about half of the total scene light: a large constant plus a small change. Bian 2024 has about 75% light throughput, while mosaic filters have less than 10%. This is exactly the case CUS is designed for. None of the four papers solves saturation, and none uses a modulo ADC.
可能被问到的问题
- Q:"4 bit" 是怎么来的? A:综述的定性说法:每个像素累加多个编码波段,8-bit 的量化台阶被分摊到多个波段上,重建后每个波段的有效位数只剩约 4 bit。不是严格推导。
- EN Q: "Where does the 4-bit number come from?"
- EN A: "It is a qualitative statement in the review. Each pixel adds up many coded bands, so the 8-bit steps are shared across many bands. After reconstruction, each band keeps only about 4 effective bits. It is not a strict derivation."
- Q:单像素那条"大常数 + 小变化"为什么是动态范围问题? A:ADC 要容纳大常数(避免饱和),小变化只占很少几个量化台阶,分辨率被浪费。modulo ADC 正好能"折掉"大常数。
- EN Q: "Why is a large constant plus a small change a dynamic range problem?"
- EN A: "The ADC must hold the large constant to avoid saturation. Then the small change uses only a few quantization steps, so resolution is wasted. A modulo ADC can fold away the large constant."
- Q:75% 透光率和动态范围什么关系? A:透光越多,同样场景下像素收到的光越多,越早饱和。这是我们的推断,论文正文没讨论。
- EN Q: "What does 75% light throughput have to do with dynamic range?"
- EN A: "More light gets through, so for the same scene each pixel gets more light and saturates earlier. This is our own inference. The paper does not discuss it in the main text."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| ⊙ (Hadamard product) | 逐元素乘积 | |
| diag(vec M) | 对角矩阵 | 把掩模拉成向量放在对角线上 |
| mosaic filter | 马赛克滤光片 | 类似 Bayer 阵列,每个像素一个窄带滤光片 |
| light throughput | 光通量 / 透光率 | 有多少入射光到达探测器 |
| saturation | 饱和 | 信号超过探测器或 ADC 上限 |
第 12 页 · H1 Yuan 2021:快照压缩成像综述
英文标题:Snapshot imaging sums coded bands in one frame; the review flags detector bit depth 论文:Yuan, Brady & Katsaggelos, Snapshot Compressive Imaging: Theory, Algorithms, and Applications, IEEE Signal Processing Magazine 2021。
这页回答什么:快照压缩成像(SCI)的统一模型、算法格局和开放问题。
页面内容逐项解释
-
左图:论文 Fig. 4,CASSI 光路(物镜、编码孔径、中继、棱镜、探测器)和采样原理(立方体 → 掩模 → 色散剪切 → 求和)。
-
Coded aperture (CASSI):一块掩模编码场景;棱镜把第 n 个波段平移 d(λₙ − λ_c) 个像素,于是一块掩模对不同波段相当于 N_λ 块不同掩模;探测器把各波段相加。
-
公式:y = Φx + e,Φ = [D₁, …, D_{N_λ}]。
-
ΦΦᵀ 是对角的:所以"投影到满足 y = Φx 的集合"只需要每个像素做一次除法,非常便宜。这是 GAP、ADMM 类算法快的原因。
-
表格(10 个模拟场景,28 波段):
方法(类别) PSNR (dB) 时间 GAP-TV(优化) 23.73 35 s DeSCI(优化) 25.86 3 h λ-net(端到端网络) 29.25 0.07 s ADMM-net(展开网络) 32.13 0.08 s 优化方法在 CPU 上计时,网络在 GPU 上。
-
Plug-and-play:迭代求解器 + 任意去噪器,换掩模不需重训;学习型网络需要重训。
-
Open problem named by the review:每个像素累加很多编码波段,8-bit 相机重建后每个像素可能只有 4 bit 动态范围;8 → 12 bit 改善了结果。
-
Conclusion for CUS:综述给出了 modulo 快照系统的前向模型 y = M_λ(Φx),并把探测器动态范围列为开放问题;████-███-████ ██████████████ ███ ██ ███ ████。
怎么讲:"这是快照压缩成像的标准综述。一个模型 y = Φx 统一了光谱、视频等各种快照系统;算法上深度方法领先 3 dB 以上,但换掩模要重训。最重要的是:综述自己把探测器位深列为开放问题,就是上一页那个 4 bit。"
English script(可直接照读):
This is the standard review of snapshot compressive imaging. One model, y equals Phi x, covers spectral, video and other snapshot systems. Because Phi Phi transpose is diagonal, the projection step is just one division per pixel, so GAP and ADMM methods are fast. On 10 simulated scenes with 28 bands, deep methods lead by more than 3 dB. ADMM-net gets 32.13 dB in 0.08 seconds, while DeSCI gets 25.86 dB in 3 hours. But learned networks must be retrained when the mask changes. Plug-and-play methods do not need retraining. Most important for us: the review itself lists detector bit depth as an open problem. That is the 4 bits from the last page. For CUS, the review gives the modulo snapshot model, y equals the modulo of Phi x, ███ ████-███-████ █████████ █████ ███████ ███ ███ ███ ███ ██████████.
可能被问到的问题
- Q:GAP 是什么? A:Generalized Alternating Projection,交替做两步:投影到满足测量的集合,再用去噪器(如 TV)去噪。
- EN Q: "What is GAP?"
- EN A: "GAP means Generalized Alternating Projection. It alternates two steps. First, it projects onto the set that fits the measurements. Then it denoises with a denoiser, such as TV."
- Q:plug-and-play 和展开网络区别? A:PnP 用现成的预训练去噪器,算法固定、灵活;展开网络把迭代展开成层端到端训练,快且准,但绑定特定系统。
- EN Q: "What is the difference between plug-and-play and unfolded networks?"
- EN A: "Plug-and-play uses a ready, pre-trained denoiser. The algorithm is fixed and flexible. An unfolded network turns the iterations into layers and trains them end to end. It is fast and accurate, but tied to one specific system."
- Q:理论保证是什么? A:综述引用 Jalali & Yuan 2019:误差上界随一帧压缩的波段数 N_t 增大而增大;先验越好(压缩码越好),一帧能压的波段越多。
- EN Q: "What theoretical guarantee does the review give?"
- EN A: "The review cites Jalali and Yuan 2019. The error bound grows with N_t, the number of bands compressed into one frame. A better prior, meaning a better compression code, lets one frame hold more bands."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| SCI (snapshot compressive imaging) | 快照压缩成像 | |
| shear / shift | 剪切 / 平移 | 色散使不同波段在探测器上错开 |
| projection | 投影 | 把估计拉到满足约束的集合上 |
| GAP-TV | 广义交替投影 + 全变分 | 经典迭代方法 |
| DeSCI | — | 基于非局部低秩先验的优化方法,慢但较准 |
| end-to-end network | 端到端网络 | 直接从测量映射到重建 |
| unfolded network | 展开网络 | 同 unrolled |
| ADMM (alternating direction method of multipliers) | 交替方向乘子法 | 常用优化算法 |
| plug-and-play (PnP) | 即插即用 | 在迭代算法里插入任意去噪器 |
第 13 页 · H2 MST 2022
英文标题:MST beats DGSMP by 2.55 dB with 4.4% of its computation 论文:Cai et al., Mask-guided Spectral-wise Transformer for Efficient Hyperspectral Image Reconstruction, CVPR 2022。
这页回答什么:快照式高光谱重建的代表性网络是怎么设计的。
页面内容逐项解释
-
顶部图:论文 Fig. 2,(a) 高光谱图像空间稀疏、光谱相关;(b) CASSI 流程;(c) 谱向注意力 + 掩模引导机制。
-
Spectral-wise attention:特征重排为 HW×C,每个波段是一个 token(而不是每个像素/窗口): Aⱼ = softmax(σⱼ Kⱼᵀ Qⱼ),headⱼ = (Mⱼ ⊙ Vⱼ) Aⱼ。 计算量随图像尺寸线性增长(2HWC²/N),同时有全局感受野。
-
Mask guidance:把掩模按色散平移 → 卷积 + sigmoid 门控 → 再平移回来,得到 M,用来按"每个像素编码得是否可靠"给 V 加权。
-
表格:
方法 PSNR (dB) 参数 (M) GFLOPS MST-L 35.18 2.03 28.15 MST-S 34.26 0.93 12.96 DGSMP 32.63 3.76 646.65 DGSMP 是论文中最强的基线。消融:谱向注意力 +0.89 dB,掩模引导再 +1.08 dB。
-
Conclusion for CUS:学到的光谱先验远强于固定字典;███ ██████████,██████████████████;但它假设掩模已知且固定,每个系统要重训。
怎么讲:"MST 把每个波段当 token 做注意力,抓住了'光谱高度相关'这个特点,还用掩模告诉网络哪里编码可靠;比 DGSMP 高 2.55 dB,计算量只有它的 4.4%。"
English script(可直接照读):
MST is a representative network for snapshot hyperspectral reconstruction. It treats each band as a token for attention, not each pixel or window. This uses the fact that the spectrum is highly correlated. The cost grows only linearly with image size, and it still sees the whole image. It also uses the mask to tell the network where the coding is reliable. MST-L reaches 35.18 dB. That is 2.55 dB above DGSMP, the strongest baseline, with only 4.4% of its computation. In the ablation, spectral-wise attention adds 0.89 dB, and mask guidance adds 1.08 dB more. For CUS: learned spectral priors are much stronger than a fixed dictionary. █████ ███ █████ ███ █████, ████ █ ███████ █████ ██████ ███ ██████. But it assumes a known, fixed mask, and must be retrained for each system.
可能被问到的问题
- Q:为什么按波段做注意力? A:高光谱图像空间上稀疏、光谱上相关;按波段做注意力计算量只和图像大小线性相关,比空间注意力省得多(论文 Tab. 2b:谱向 33.18 dB vs 空间 32.67–32.86 dB)。
- EN Q: "Why do attention over bands instead of over space?"
- EN A: "Hyperspectral images are sparse in space and correlated across the spectrum. Attention over bands costs only linear time in image size, so it is much cheaper than spatial attention. In Table 2b of the paper, spectral-wise gets 33.18 dB, and spatial gets 32.67 to 32.86 dB."
- Q:GFLOPS 是什么? A:每次前向计算的十亿次浮点运算数,衡量计算量。
- EN Q: "What does GFLOPS mean?"
- EN A: "It is the number of billions of floating-point operations in one forward pass. It measures computation."
- Q:对 CUS 有什么直接用? A:CUS 现在用固定 DCT 字典 + OMP;██████████████,████████ ███ ██████、████████,██████████████████。
- EN Q: "How is this directly useful for CUS?"
- EN A: "CUS now uses a fixed DCT dictionary and OMP. ██ ███ ██████ ██ █ ████████, ██ ███ █████ ███ ███ ██ ████ ███ █████ ███ ██████ ███ ████████████. ████ ██ ███ ███ ████ ████ ██ ███████ ███ ██████████████."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| transformer | Transformer | 基于注意力机制的网络 |
| self-attention | 自注意力 | 用 Q、K、V 计算元素之间的关联 |
| token | 词元 / 令牌 | 注意力机制处理的基本单位 |
| Q / K / V (query / key / value) | 查询 / 键 / 值 | 注意力的三组向量 |
| softmax | — | 归一化为概率的函数 |
| receptive field | 感受野 | 输出能"看到"的输入范围 |
| sigmoid gate | sigmoid 门控 | 用 0–1 的权重控制信息通过 |
| ablation | 消融实验 | 逐个去掉组件看贡献 |
| GFLOPS | 十亿次浮点运算 | 计算量 |
| DGSMP | — | 一种基于深度高斯尺度混合先验的展开方法 |
第 14 页 · H3 Bian 2024:片上宽带高光谱传感器
英文标题:Broadband filters on the chip give 96 channels at 124 fps with 75% light throughput 论文:Bian et al., A broadband hyperspectral image sensor with high spatio-temporal resolution, Nature 2024(我们读的是 arXiv 版本)。
这页回答什么:不需要额外光学元件的快照式高光谱方案。
页面内容逐项解释
-
左图:我们画的示意图(论文 arXiv 版没有图片):4×4 共 16 种材料的单元,以及三条示意性的宽带、弱相关透过率曲线。
-
Measurement(我们的记号,论文没给公式):像素 p 记录 y_p ≈ r_pᵀ s_p,一个值对应最多 96 个未知数;4×4 单元里的邻居提供多样性。
-
Making it:16 种材料从 35 种候选中用进化算法挑选(按重建精度);每种材料一次光刻,单元 10 µm,固化在裸单色芯片上;用单色仪标定每个像素的响应;用 SRNet(谱向注意力网络)解码。
-
表格:
CMOS 版 InGaAs 版 波长范围 400–1000 nm 400–1700 nm 通道数 61 96 像素、帧率 2048², 47 fps 1024², 124 fps 透光率 71.8% 74.8% 其他方案透光率:马赛克滤光片 < 10%,编码孔径 < 50%,RGB < 30%。
-
Demonstrations:低光(290 lux,1 ms 曝光)、叶绿素指数(误差 3.73%)、100 Hz 血氧、纺织品分类 98.15%、月球成像。
-
Conclusion for CUS:透光多 → 亮场景更早饱和(我们的推断;动态范围只在补充材料里)。██████ ██████████████████████████,████████ ███ ██████。
怎么讲:"芯片上直接做 16 种宽带滤光材料,每个像素对光谱编码不同,网络解出 96 个通道。透光率 75% 很高,对弱光好;但对我们来说意味着亮场景更容易饱和——这是我们的推断。"
English script(可直接照读):
This is a snapshot hyperspectral sensor with no extra optics. They put 16 broadband filter materials directly on the chip, in 4 by 4 cells. Each pixel codes the spectrum differently, and a network called SRNet decodes up to 96 channels. The InGaAs version gives 96 channels at 124 frames per second, from 400 to 1700 nanometers. Light throughput is about 75%, which is very high. Mosaic filters have less than 10%, and coded apertures less than 50%. This is good for low light. But for us, it means bright scenes saturate earlier. This is our own inference; the paper does not say this. ████, █ ██████ █████ █████ ████ ████ █████ ██ ██, █████ ██ ███████ ███ ███████ ███ █████. Note that the figure on the left is our own sketch, because the arXiv version has no figures.
可能被问到的问题
- Q:为什么是宽带而不是窄带滤光片? A:窄带只透过一小段光谱,光损失大;宽带透过大部分光,靠不同滤光片曲线的差异编码光谱,再由网络解码。
- EN Q: "Why use broadband filters instead of narrowband ones?"
- EN A: "A narrowband filter passes only a small part of the spectrum, so it loses a lot of light. A broadband filter passes most of the light. The spectrum is coded by the differences between the filter curves, and a network decodes it."
- Q:图是论文的吗? A:不是,arXiv 版本没有图片,这是我们画的示意,页面上注明了。正式版 Nature 论文 Fig. 1 有结构图。
- EN Q: "Is that figure from the paper?"
- EN A: "No. The arXiv version has no figures, so this is our own sketch. The slide says so. The published Nature paper has a structure figure in Fig. 1."
- Q:动态范围论文怎么说? A:正文只说"见补充材料 6.4",我们手上的 arXiv 版本没有补充材料。
- EN Q: "What does the paper say about dynamic range?"
- EN A: "The main text only says, see Supplementary Section 6.4. The arXiv version we have does not include the supplementary material."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| on-chip | 片上 | 直接做在图像传感器上 |
| broadband filter | 宽带滤光片 | 透过较宽波长范围 |
| photolithography | 光刻 | 半导体图形化工艺 |
| evolutionary search | 进化搜索 | 用遗传算法类方法选材料组合 |
| monochromator | 单色仪 | 输出单一波长的光,用于标定 |
| calibration | 标定 | 测量每个像素的光谱响应 |
| CMOS / InGaAs | 互补金属氧化物半导体 / 铟镓砷 | 两种传感器;后者覆盖近红外 |
| fps (frames per second) | 帧每秒 | |
| SRNet | — | 论文的重建网络,U 形结构 + 谱向注意力 |
第 15 页 · H4 Xu 2024:系统与测量
英文标题:One DMD codes space and spectrum together for a single detector 论文:Xu, Lu, Saragadam & Kelly, A compressive hyperspectral video imaging system using a single-pixel detector, Nature Communications 2024。
这页回答什么:单像素高光谱视频系统怎么搭,怎么编码。高光谱部分和 CUS 最接近的一篇。
页面内容逐项解释
- 左图:论文 Fig. 1,(a) 光路:DMD 左半边、狭缝、光栅、45° 反射镜、DMD 右半边、单探测器、ADC;(b) 矩阵模型。
- Optics:场景成像到 DMD 左半边(加空间 pattern);光经过狭缝和光栅,光谱落到同一块 DMD 右半边(加光谱 pattern);剩下的光进入一个探测器。
- Measurement:y = φ_sᵀ X_t φ_w;一帧内 b_t = R_t Φ vec(X_t),R_t 选出实际测量的那些 (空间, 光谱) pattern 对。
- Patterns:
- 空间:多分辨率结构化随机 pattern(STOne);光谱:打乱顺序的 Walsh–Hadamard。两者每次测量都换。
- DMD 只能开/关(0/1),通过加互补 pattern 得到 ±1 pattern。
- pattern 速率 5 kHz;每帧用 2048 个测量,窗口步长 1024(50% 重叠)。
- Key innovation:空间、光谱、时间在一块 DMD 上联合编码,每次测量约对应 900 个体素;只压缩光谱的方案最多约 20。
怎么讲:"这是唯一一个'单探测器 + 顺序测量'的高光谱系统,和 CUS 的测量模型完全一致。关键是空间和光谱 pattern 每次都换,所以压缩比能到 900:1。"
English script(可直接照读):
For hyperspectral, this paper is the closest to CUS. It is the only system with one detector and sequential measurements, which matches the CUS measurement model exactly. One DMD does both jobs. The left half codes the image. Then a slit and a grating spread the spectrum onto the right half, which codes the wavelengths. All remaining light goes to one detector. The DMD is only 0 or 1, so they add complementary patterns to get plus or minus 1 patterns. The key point is that the spatial and spectral patterns change at every measurement. So each measurement covers about 900 voxels, while spectrum-only compression reaches about 20. That is why the compression ratio can reach 900 to 1.
可能被问到的问题
- Q:为什么一块 DMD 能同时做空间和光谱编码? A:左右两半分别用:左半边编码图像,经色散后光谱展开在右半边,右半边按波长开关镜子。
- EN Q: "How can one DMD do both spatial and spectral coding?"
- EN A: "It uses the two halves separately. The left half codes the image. After dispersion, the spectrum spreads over the right half, and the right half switches mirrors by wavelength."
- Q:±1 pattern 怎么来? A:DMD 只有 0/1;测一次 pattern 再测一次它的互补(0/1 翻转),两者相减相当于 ±1 pattern。
- EN Q: "How do you get plus-or-minus-1 patterns?"
- EN A: "The DMD only has 0 and 1. We measure one pattern, then its complement, with 0 and 1 flipped. Subtracting the two gives a plus-or-minus-1 pattern."
- Q:900:1 怎么算? A:重建的体素数除以测量数(论文定义)。
- EN Q: "How is the 900 to 1 ratio computed?"
- EN A: "It is the number of reconstructed voxels divided by the number of measurements. This is the paper's definition."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| single-pixel detector | 单像素探测器 | 只有一个感光单元 |
| slit | 狭缝 | 限制光束宽度 |
| grating | 光栅 | 色散元件,把不同波长分开 |
| complementary pattern | 互补图案 | 0/1 取反 |
| STOne | — | 一种多分辨率结构化随机 pattern,任意连续 4k² 个能构成 2k×2k 完整集合 |
| Walsh–Hadamard | 沃尔什–哈达玛 | ±1 正交矩阵,常用于单像素成像 |
| permuted | 打乱顺序的 | |
| compression ratio (CR) | 压缩比 | |
| sliding window / stride | 滑动窗口 / 步长 |
第 16 页 · H4 Xu 2024:重建与结果
英文标题:4D total variation with optical flow recovers 128×128×64-band video at 4.3 fps
这页回答什么:怎么从一串测量重建高光谱视频;效果如何;对 CUS 意味着什么。
页面内容逐项解释
- 左图:论文 Fig. 2a,测量流 → 滑动窗口 → 灰度视频 → 光流 → 4DTV + 光流的高光谱重建。
- Two stages:
- 用互补 pattern 的测量先重建灰度视频,估计帧间光流(运动)。
- 重建高光谱视频:min |∇₄u| + |∂_{t,flow} u| s.t. ‖b − RSu‖ < ε。x、y、t 用一阶差分,λ 用二阶差分;用原始–对偶方法求解。
- Fast version:CNN + LSTM(看 5 帧),GPU 上每 16 个 patch 0.04–0.06 s,但压缩比较低(最多 100:1);优化方法每帧 390 s。
- Results:361–827 nm,每波段 7.3 nm;18.4 万个测量重建 157 帧;光谱与光纤光谱仪一致;不用光流时颜色丢失。
- Conclusion for CUS:和 CUS 最接近:单探测器、一串顺序测量、每个值约为场景总光的一半。但相邻 pattern 是伪随机的,差别大;███ ████████████ ███████ ██████████████████。
怎么讲:重点读结论框:"测量模型完全对上,动态范围问题也有(每个值是总光的一半),唯一不对的是 pattern 顺序——████ ███ ██████████。████想法 3 ████████。"
English script(可直接照读):
Reconstruction has two stages. First, they rebuild a grayscale video and estimate optical flow. Then they rebuild the hyperspectral video with 4D total variation plus optical flow. The result is 128 by 128 pixels, 64 bands, at 4.3 frames per second. Now the conclusion box. The measurement model matches CUS completely. The dynamic range problem is also there, because each value is about half of the total light. The only thing that does not match is the pattern order. Neighbouring patterns are pseudo-random and very different. ██ ███ ███, ██ ████ ███████ ████ ██ ████ ████████████ ████████ ██████ ██ ████ █ ███ ███████. ████ ██ ███ ████████ █████ ██ Idea 3.
可能被问到的问题
- Q:光流是干什么的? A:估计物体在帧间怎么运动,让重建时可以沿运动方向约束平滑,避免运动模糊、保留颜色。
- EN Q: "What is the optical flow for?"
- EN A: "It estimates how objects move between frames. Then reconstruction can enforce smoothness along the motion. This avoids motion blur and keeps the colours."
- Q:为什么光谱方向用二阶差分? A:光谱曲线大多分段线性,二阶差分惩罚弯曲而不惩罚斜坡。
- EN Q: "Why use a second-order difference along the spectrum?"
- EN A: "Most spectral curves are piecewise linear. A second-order difference penalizes bending, but not slopes."
- Q:快速版和 900:1 能同时达到吗? A:不能。网络版受 GPU 显存限制,只在分块、压缩比 ≤ 100:1 下测试;900:1 是优化方法的结果。slides 上已注明。
- EN Q: "Can the fast version also reach 900 to 1?"
- EN A: "No. The network version is limited by GPU memory. It was only tested on patches, with compression ratio up to 100 to 1. The 900 to 1 result is from the optimization method. The slides note this."
- Q:论文有没有提 ADC 或饱和? A:没有。只提到分辨率提高会降低每次测量的 SNR。
- EN Q: "Does the paper mention the ADC or saturation?"
- EN A: "No. It only says that higher resolution lowers the SNR of each measurement."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| 4D total variation (4DTV) | 四维全变分 | 在 x, y, λ, t 上的平滑正则 |
| optical flow | 光流 | 帧间像素运动场 |
| primal–dual method | 原始–对偶方法 | 一类凸优化求解器 |
| LSTM (long short-term memory) | 长短期记忆网络 | 处理序列的循环网络 |
| patch | 图像块 | |
| fiber spectrometer | 光纤光谱仪 | 用来测参考光谱 |
第 17 页 · 横向对比表
英文标题:No paper uses a modulo ADC; MRI and single-pixel video measure in a controllable order
这页回答什么:把 7 篇论文按 CUS 的三个条件(加一列"是否讨论 ADC 位深")放在一起比,看 CUS 最适合用在哪里。
页面内容逐项解释
| 论文 | 探测器量程被证明是瓶颈? | 测量按顺序? | pattern 顺序可控? | 讨论 ADC 位深? |
|---|---|---|---|---|
| M1 低场扫描仪 | 否(SNR 低) | 是(k-space 线) | 是(序列) | 否 |
| M2 扩散重建 | — | 是 | 任意 mask 都行 | 否 |
| M3 BJORK 轨迹 | — | 是 | 是,学习得到 | 否 |
| H1 快照综述 | 是(每波段 4 bit) | 否(一帧) | 掩模设计 | 是,开放问题 |
| H2 MST | — | 否 | 否,固定掩模 | 否 |
| H3 片上传感器 | 只在补充材料 | 否 | 否,固定滤光片 | 否 |
| H4 单像素视频 | 隐含(每值约为一半光) | 是(一串) | 是,可编程 | 否 |
| Behin 2005(背景) | 是(20 bit,7 T 小鼠) | 是 | — | 是 |
| Gabr 2009(背景) | 临床 3 T 否(≤ 约 15 bit) | 是 | — | 是 |
"—"表示不适用或未讨论。底部注明:7 篇近期论文按引用和相关性选取,结果引自论文,未复现。
- Conclusion:位深不够有文献记录的,是临床前 MRI、线圈上数字化和快照光谱成像,近期的成像论文都没有解决;顺序且可控的测量存在于 MRI 和单像素视频中——这是 CUS 适合的地方。
怎么讲:按列读,一列一句:"第一列,证明量程不够的是 Behin(临床前 20 bit)和 Yuan(快照式每波段只剩 4 bit),而 Gabr 说明临床 3 T 约 15 bit、现有 ADC 够用;第二、三列,能按顺序、可控测量的是 MRI 和单像素;最后一列,没有一篇用 modulo ADC。所以 CUS 适合 MRI 轨迹和单像素视频。"
English script(可直接照读):
This table compares the 7 papers, column by column. First column: who shows that the detector range is not enough? Behin shows it for preclinical MRI: 20 bits. Yuan shows it for snapshot spectral imaging: maybe only 4 bits per band. Gabr shows that clinical 3 tesla needs about 15 bits, so current ADCs are enough there. Second and third columns: which systems measure in order, with an order we can control? That is MRI and single-pixel video. Last column: no paper uses a modulo ADC. So the range problem is documented in preclinical MRI, on-coil digitization and snapshot spectral imaging, and recent imaging papers do not solve it. Ordered, controllable measurements exist in MRI and single-pixel video. That is where CUS fits: MRI trajectories and single-pixel video.
可能被问到的问题
- Q:快照式既然量程问题最明确,为什么不是首选? A:它不满足"按顺序测量",CUS 的差分思路用不上,████████████████(想法 4)。
- EN Q: "Snapshot imaging has the clearest range problem. Why is it not your first choice?"
- EN A: "Snapshot imaging does not measure in sequence. So the difference idea in CUS does not work there. ██ █████ ████ █ ███ ██ ██████. ████ ██ Idea 4."
- Q:M2 为什么"量程"一栏是"—"? A:它是重建方法,不涉及硬件采集。
- EN Q: "Why is the range column a dash for M2?"
- EN A: "M2 is a reconstruction method. It does not deal with hardware acquisition."
- Q:Gabr 那一行为什么是"否"? A:它实测临床 3 T 信号约 15 bit 以内,现有扫描仪 ADC 加过采样有效约 17 bit,够用。它放在表里,是为了说明 MRI 的位深压力集中在临床前和线圈上数字化,而不是临床控制台。
- EN Q: "Why does the Gabr row say no?"
- EN A: "Gabr measured clinical 3 tesla signals. They need about 15 bits or less. Current scanner ADCs with oversampling give about 17 effective bits, so that is enough. We keep Gabr in the table to show where the pressure is. For MRI, the bit-depth problem is in preclinical systems and on-coil digitization, not in the clinical console."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| controllable | 可控的 | |
| implied | 隐含的 | 论文没直接说,但可推出 |
| supplement | 补充材料 |
第 18 页 · MRI 方向的想法
英文标题:Ideas for MRI: █████ ███████ ████ ██████████ ██████ ███ ████ █████████ ██████████████
这页回答什么:基于调研,MRI 方向有哪些可以讨论的想法。每个想法按"想法 / 依据 / 待讨论"三行。
Idea 1: ███████-█████ ██████████ ████████(██████████████████)
- 想法:████████ ████ ████ ████████,██████████████████████████(██████ █-█████ ██████████████████████)。
- 依据:█████ ██████ █████ ██████████████ █ ████████;█-█████ ██████████████ ███ ████(█████)。
- 待讨论:████████████ ████ ████████;██████████████████████████?██████████████████,████████████████?
Idea 2: █████████ ██████████████ ████ ██████ ████(████ ██████ ██████████████)
- 想法:████ █████ ████,██████████████████████████████,██████████ ██ ████████████、████、██████。
- 依据:██████████████████████████,██████ ‖█*█‖ ≤ █ ██ █ ████;██████████████(█████ & ██)。
- 待讨论:██████ ████████,██████████████████████████;██████████████████ ██–██ ████。
怎么讲:每个想法读"想法"一行,再挑"待讨论"里一个问题抛给听众。
English script(可直接照读):
████ ███ █ █████ ███ ███. Idea 1 ██ ███████-█████ ██████████ ████████. █████ ███████ ██████ ██████ █-██████ ████████████ ████ █ █████ ███████. ██ █████ ███ ███ ████ ███████ ████ ██ ███ ████████ ███ ████ ████ █████████. ██ █████████ ███ ██████ ██ ████ █████, ██ ███ ██████ ███████ ████████ ███████ ████ ███ █-█████ ██████. ███ █-█████ ██████ ████ ████ ███ ███ █████, ██ █████ █████. ████████ ███ ███: ███ ████ █████ ███████ █████ ████ ███████ ███████ ████? Idea 2 ██ █████████ ██████████████ ████ ██████ ████. ██ ████ ███ █████ █████. ██ ██████ ████ ███ ████ ███████████ ████ ██ █ ███████-█████ ███. ███ ███████, ██ ███ ███ ███████ ██ ██ ████████ ███ ████ ████████, ██████, ███ ███████ █████. █████ ███ ██ ████ ███ █████ ████ ███ ██████ ██ ████████. ████████ ███ ███: ████ ██████ ████ ██ ███-██████, ██ ███ █████'█ █████ ██ ██████ █████. ██ ████ █ ███████?
可能被问到的问题
- Q:████████████████?████████████? A:████████████████ ███████ ████(███ ████████ █ ████████),██████████ █-█████ ██████,██████████,████████████。(████████:████████████████████,██████████████████ █× ██ ██%,█× ██ ██%。████ ██████ ██,██████████████████。)
- EN Q: "████ ██ ███████ ████████████, ███ ███ ████ ██ ████?"
- EN A: "██ ██████ ████ ███████ ████ ███████ █████ ███ ███████ █████████. ███ ███████ ████ █ █████ ████████████ ██ ███████. ████ ████████ ███████ ███ ██████ ██ █-█████. █████ ██████ ██████ ████, ██ █████ ███ █████ ████ █████. ██ ███ ███████ ██████████, ██ ███ ██████ ███████ ████, ███ ████ ████ ███ █████ ██ ███████ ██ █ █████ ███ █████ ██ ███████ ██ █ █████. ████ ██ █ ███████████ ██████ ████ ███ ███ ██████████. ██ ██ ███ ██ ███ ██████."
- Q:想法 2 ██████"██ ███ ████████████████"██████████? A:██████████████,████████████████;████████████████████████████████████████████████。
- EN Q: "███ ██ Idea 2 █████████ ████ █████ ██████████ ████ ███ ███ ████ ███████ █████████ ██████████████?"
- EN A: "████ ██ █ █████ ██ █ █████. ██ ███ ██████████ ██ █████, ███ ██████ ████ ██████ ███ ██. ██ ██ █████████ ███ █ █████, ███ █████ ████████ ████ ██████. ████ ██ ████████ ███ ████████ ██ ███ ████ ████████."
- Q:██████████? A:██████████████████████;████████████████████(████ ███ ████████████)。
- EN Q: "████ ███████ ███████ ██ ██ ███-██████?"
- EN A: "██ ███ ████████ ██ █ █████ ███████ ████████. ██ ████ █ ████ ██████████████ ██ █ ██████████. ███ ███████, ███ ██████ █████ ████ ██ ███."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| ███████-█████ | ██████████ | |
| ███████ ████████████ | ██████████ | ████████████████ ███████ ████ |
| ███████ ████ | ██████████████ | ██████████████████████ |
| ███-██████ | ████ | ██████████████████████ |
| ████ ████████ | ████████ | ██ █,████████████████ |
第 19 页 · 光谱成像方向的想法
英文标题:Ideas for spectral imaging: █████ ██████-█████ ████████ ████████; ████ ████████ ██████ ██ ██
Idea 3: ████████ ███████ ██████-█████ ████████(████████████████ ███████)
- 想法:██ ████████ ████████████ ███████ ████,██████ ███████ ████████████;██████ ██████ ██████ + ███ ████。
- 依据:██████████████████████████████,███████ ██████████(██ ████);████████ ███████ ██████ ███ ████████(█████████ ████,██ ██████ ██)。
- 待讨论:██████████████████████████(████████████████████████),████████████████████?██████████████████████████–████████████████。
Idea 4: ██ ██████ ██████ ███ ████████ ████████ ███████(██████████████████ ██████ ████)
- 想法:████████████████████████████ ██████ ████,████████████████████████。
- 依据:██████████████████████████████████(████ ████);██████████████ ██%(████ ████);██████████████ ██████ ████████████(██████ ████,██ ██████ ██)。
- 待讨论:██████████████████,████████,████████████ ███ ██████████。
怎么讲:读完后直接转入讨论,例如问:"这四个想法里,哪个最值得先做?"
English script(可直接照读):
████ ███ █ █████ ███ ████████ ███████. Idea 3 ██ ████████ ███████ ██████-█████ ████████. ██ █████ ████████ ██ ██████████ ██████ ████████, ██ ████ ████████ ████████ ██████ ██ ████ █ ███ ███████. ████ ██ ███ █ ██████ ████████ ███ ███████ ████ ███. ██ ██████-█████ ███████, ████ █████ ██ █████ ████ ██ ███ █████ █████, ███ ███ ████████ ███ ██████-██████, ██ ██ ████ █████. ██████-█████ ████████ ███████ ███ ██ ███ █████; █████████ ████ ██ ██ ██████ ██. ███ ████ ████████: ████ ██████ ████████ ████ ███████████, █████ ██ ███ █████-███ ██ ███ █████, ██ ████ ███ █████-██████████ ███████? Idea 4 ██ ██ ██████ ██████ ███ ████████ ████████ ███████. ████████ ███████ █████ ██ ██ ████ ███████ █████, ██ ████ ████ ████. ███ ███████ ██ ███ █████ ██ ██, ███ █████ █ ████████. ██ ██ ██ █ █████████ ████████ ███████ ████ ███. ███ █ ████████ ███ ████████: ██ █████ █ █████, █████ ███ ██████ ██ ███ █████?
可能被问到的问题
- Q:████████████"██████████████"██████? A:████████ ████ ██████████(██████████████),██ █████████ ██████████;████████ ██████████████████(██████、█████ ██、███████ █████ ████),████████ ███ ████████。
- EN Q: "██ █████ ██ ████████ ████████ █████ █████████ ██████ ██ ████ █ ███ ████?"
- EN A: "██ ███ ███ ███ ████ ██ ████ ████ ████████, █████ ████████ █████ ██████ ██ ████ █ ███. ██ ███ ██████ █████ ████ █████████. ████████ ████ ████ ████ ███████ █████ ██████, ███ ███████ ███████ █████, █████ █████ ███ ███████ █████ █████. █████ ███ ████ ███ █████ █████ █████."
- Q:██████████████████████? A:████ ███████ ████,██████████████████████████(█_███)████,████████████;████████████████ █_█ ██████。
- EN Q: "████ ██ ███ ████ ██ ██████ ███████████?"
- EN A: "██ ████████ ████████ ███ ███ ███████, ███ ███████ ██ ███ ██████████ ███████████ ██████ ██████ ████ ███████. ████ ██, ██ ███ ████ ██. ████ ██████ ████████ ██ ██████. ████ ██ ███████ ███ █████ █ █████-███ ██ ███ █████."
- Q:想法 4 ██ ██████ ████████████████? A:██████ ████████(██████)██████████████;████████████████████████████████████,████████████████████,██████████████████████。
- EN Q: "███ ██ Idea 4 █████████ ████ ██████'█ ██████?"
- EN A: "██████ ████ ██ ██████████ ██ ██████, ███-███████████ ██████. ██ █ ████████ ██████, ████ █████ ██ █ ███ ██ ███████ █████ █████. ██ █████ ██████████, ██ █████ ████ ███████████ ██████████████. ███ ██ ███████ █████ █ █████ ██ █ ███ ███████."
生词表
| 英文 | 中文 | 注释 |
|---|---|---|
| ███████████ | ████████ | ██████████████████████,██████████████ |
| █████-██████████ ███████ | ████████████ | ██ ██ █████ ███████ ████████████████████ |
| ████ ████ | ██████ | ████████████████ |
| ██████ █████ | ████████ |
Backup 页
正文不讲,有人追问时翻到对应页。
| Backup 内容 | 何时翻 |
|---|---|
| Liu 2021 详细版两页(硬件、干扰消除) | 有人问 0.055 T 系统细节、协议、成本 |
| Yuan 2021 详细版两页(CASSI、算法家族与公式) | 有人问 GAP/ADMM 公式、理论保证 |
| MST 详细版两页(架构图、完整表格与消融) | 有人问网络结构、TSA-Net 对比 |
| Bian 2024 详细版两页(材料与制作、性能与应用) | 有人问制作工艺、应用细节 |
| S1 Kilcullen 2022:MHz 级单像素成像(循环移位 pattern) | ████想法 3 ██ |
| S2 Wang 2023:中红外单光子单像素成像 | 有人问单像素在其他波段的价值 |
| S3 Monroy 2025:高阶差分 modulo 成像 | ████想法 4 ██;也和我们的差分步骤、DCT 字典相似 |
| 参考文献 |
数字出处表
| slide 页 | 数字 | 出处 |
|---|---|---|
| 1, 4, 17 | 小鼠全身 7 T 需要 20 bit;16-bit 系统 gain stepping SNR 提高 8 倍 | Behin, Bishop & Henkelman, Concepts Magn. Reson. B, 2005(摘要) |
| 5 | 约每百万人 7 台;0.055 T;约 750 kg;< 1200 W;10 个感应线圈;5 层 CNN;噪声差 5%;30 min;25 名患者;< 2 万美元 | Liu 2021, pp. 2–12, Fig. 1–2 |
| 6 | σ(t) 0.01–378;约 2.5 万张 320×320;2000 步 | Chung & Ye 2022, Eq. 19, p. 6–7 |
| 7 | 30.46/22.84、33.32/33.01、29.17/16.66、32.43/33.15 dB;10–20 min | Chung & Ye 2022, Table 1–2, Sec. 4.2 |
| 8 | 26.7 mT/m、150 T/m/s;6 次展开 | Wang (BJORK) 2022, p. 5 |
| 9 | SSIM 0.957→0.968,非参数 0.959;PSNR 表;3 T、32 线圈、11 s vs 1:40;1520 张测试切片 | BJORK 2022, Fig. 5, Table I–II |
| 10 | 256×256×28(450–650 nm);128×128×64 | MST 2022 / Xu 2024 |
| 11, 12 | "每个像素可能只有 4 bit";8→12 bit 改善 | Yuan 2021, Sec. XI-A, pp. 37–38 |
| 12 | GAP-TV 23.73 dB/35 s;DeSCI 25.86/3 h;λ-net 29.25/0.07 s;ADMM-net 32.13/0.08 s | Yuan 2021, Table V, p. 35 |
| 13 | MST-L 35.18 dB/2.03 M/28.15 G;MST-S 34.26/0.93/12.96;DGSMP 32.63/3.76/646.65;+0.89、+1.08 dB;2.55 dB、4.4% | Cai 2022, Table 1, 2a, 2c, p. 6 |
| 11, 14 | 透光 71.8%/74.8%;马赛克 < 10%、编码孔径 < 50%、RGB < 30%;61/96 通道;2048²@47 fps、1024²@124 fps;16 种材料选自 35 种;3.73%、98.15%、290 lux | Bian 2024 (arXiv), pp. 3–8, 17–19 |
| 15, 16 | 5 kHz;2048 测量/帧,步长 1024;CR ≈ 900:1;128×128×64;4.3 fps;361–827 nm,7.3 nm;157 帧/18.4 万测量;390 s/帧;0.04–0.06 s/16 patch(CR ≤ 100:1) | Xu 2024, pp. 3–13, Table 2 |
| 1, 4, 17 | 临床 3 T 信号 < 约 15 bit(88 dB);心脏/脊柱 < 约 12.5 bit(75 dB);3D SSFP 体模约 16 bit(96.2 dB);14-bit 80 MHz ADC 过采样后约 17 bit(103 dB);1.5 T 少约 1 bit;留约 1.7 bit(10 dB)余量 | Gabr et al., J. Magn. Reson. 2009(PMC2873084) |
| 1, 4 | 线圈上 CMOS 接收机约 13.6 bit(81.9 dB);Δ-Σ 光纤链路约 13.5 bit(81 dB @ 200 kHz);16-bit 理想 98 dB,100 MHz/0.3 ps 抖动时约 12 bit(74.5 dB) | Sporrer 2017;Fan 2021;Scott, ISMRM 2016 |
| 口头补充(不在 slides) | 体模中心读出线跳变比例 31%(1×)、12%(8× 过采样);k-space 峰值高于中位数 55 dB | 我们的模拟:make_figs.py,结果在 fig/mri_kspace_folding_stats.csv |
总术语表(按字母顺序)
| 英文 | 中文 | 注释 | 首次出现页 |
|---|---|---|---|
| ablation | 消融实验 | 逐个去掉组件看贡献 | 13 |
| accommodate | 容纳、兼容 | 1 | |
| ADC (analog-to-digital converter) | 模数转换器 | 1 | |
| ADMM | 交替方向乘子法 | 12 | |
| automatic differentiation | 自动微分 | 8 | |
| B-spline | B 样条 | 分段多项式,天然平滑 | 8 |
| band | 波段 | 10 | |
| bit depth | 位深 | 1 | |
| broadband filter | 宽带滤光片 | 14 | |
| calibration | 标定 | 14 | |
| Cartesian / radial / spiral | 笛卡尔 / 径向 / 螺旋(轨迹) | 3 | |
| CASSI | 编码孔径快照光谱成像仪 | 10 | |
| clipping | 削顶 | 4 | |
| CNN | 卷积神经网络 | 5 | |
| coarse-to-fine | 由粗到细 | 9 | |
| coded aperture | 编码孔径 | 10 | |
| coil sensitivity map | 线圈灵敏度图 | 6 | |
| complementary pattern | 互补图案 | 15 | |
| compression ratio (CR) | 压缩比 | 15 | |
| compressive measurement | 压缩测量 | 2 | |
| conjugate gradient (CG) | 共轭梯度法 | 8 | |
| data consistency | 数据一致性 | 6 | |
| data cube | 数据立方体 | 10 | |
| DCT | 离散余弦变换 | 2 | |
| decimation | 抽取率 | 9 | |
| demodulate to baseband | 解调到基带 | 3 | |
| denoising score matching | 去噪分数匹配 | 6 | |
| dictionary | 字典 | 2 | |
| DMD | 数字微镜器件 | 10 | |
| dB (decibel) | 分贝 | 1 bit ≈ 6 dB | 4 |
| dynamic range | 动态范围 | 用 bit 表示 | 4 |
| EMI | 电磁干扰 | 5 | |
| end-to-end network | 端到端网络 | 12 | |
| excitation | 激发 | 3 | |
| FFT / NUFFT | 快速 / 非均匀快速傅里叶变换 | 3 | |
| fold / folding | 折叠 | 2 | |
| GAP-TV | 广义交替投影 + 全变分 | 12 | |
| gain stepping | 增益步进 | 1 | |
| GFLOPS | 十亿次浮点运算 | 13 | |
| gradient coil | 梯度线圈 | 3 | |
| grating | 光栅 | 15 | |
| ████ ████ | ██████ | 19 | |
| hinge penalty | 铰链惩罚 | 8 | |
| hyperspectral image | 高光谱图像 | 10 | |
| I / Q | 同相 / 正交分量 | 3 | |
| ███████████ | ████████ | 19 | |
| in vivo | 活体 | 9 | |
| Jacobian | 雅可比矩阵 | 8 | |
| k-space | k 空间 | 1 | |
| Langevin dynamics | 朗之万动力学 | 6 | |
| Larmor frequency | 拉莫尔频率 | 3 | |
| lattice | 格点 | 2 | |
| light throughput | 透光率 / 光通量 | 11 | |
| LSB | 最低有效位 | 4 | |
| LSTM | 长短期记忆网络 | 16 | |
| magnitude image | 幅值图像 | 6 | |
| modulo ADC | 取模 ADC | 1 | |
| monochromator | 单色仪 | 14 | |
| mosaic filter | 马赛克滤光片 | 11 | |
| multi-coil / parallel imaging | 多线圈 / 并行成像 | 7 | |
| noise floor | 噪声底 | 4 | |
| ███-██████ | ████ | 18 | |
| non-expansive | 非扩张 | 6 | |
| ███████ ████ | ██████████████ | 18 | |
| OMP | 正交匹配追踪 | 2 | |
| optical flow | 光流 | 16 | |
| pattern | 测量图案 | 2 | |
| photolithography | 光刻 | 14 | |
| plug-and-play (PnP) | 即插即用 | 12 | |
| posterior sample | 后验样本 | 7 | |
| precession | 进动 | 3 | |
| predictor / corrector | 预测器 / 校正器 | 6 | |
| primal–dual method | 原始–对偶方法 | 16 | |
| prior | 先验 | 3 | |
| prospective | 前瞻性(实验) | 9 | |
| pseudo-random | 伪随机 | 1 | |
| PSNR | 峰值信噪比 | 7 | |
| quantization noise | 量化噪声 | 4 | |
| readout | 读出(线) | 3 | |
| ███████ ████████████ | ██████████ | 18 | |
| root-sum-of-squares | 平方和开根号 | 7 | |
| sampling mask | 采样掩模 | 6 | |
| saturation | 饱和 | 11 | |
| score | 分数(对数密度梯度) | 6 | |
| score-based diffusion model | 基于分数的扩散模型 | 6 | |
| SCI | 快照压缩成像 | 12 | |
| self-attention | 自注意力 | 13 | |
| SENSE | 灵敏度编码 | 7 | |
| shot | 一次激发的轨迹段 | 3 | |
| slew rate | 梯度切换率 | 8 | |
| slit | 狭缝 | 15 | |
| SNR | 信噪比 | 4 | |
| snapshot | 快照式 | 10 | |
| sparse | 稀疏 | 2 | |
| spin | 自旋 | 3 | |
| spoke | 辐条 | 9 | |
| SSIM | 结构相似性 | 9 | |
| STOne | 多分辨率结构化随机图案 | 15 | |
| supervised baseline | 监督学习基线 | 7 | |
| thermal noise | 热噪声 | 4 | |
| token | 词元 | 13 | |
| total variation (TV) | 全变分 | 10 | |
| trajectory | 采样轨迹 | 3 | |
| ultra-low-field | 超低场 | 5 | |
| unrolled / unfolded network | 展开网络 | 8 | |
| unwrapping | 展开 / 解缠 | 10 | |
| VE-SDE | 方差爆炸型随机微分方程 | 6 | |
| voxel | 体素 | 10 | |
| Walsh–Hadamard | 沃尔什–哈达玛 | 15 |