CN117523025B - 一种基于光强分布调校的数字仿体生成方法 - Google Patents
一种基于光强分布调校的数字仿体生成方法 Download PDFInfo
- Publication number
- CN117523025B CN117523025B CN202410016633.5A CN202410016633A CN117523025B CN 117523025 B CN117523025 B CN 117523025B CN 202410016633 A CN202410016633 A CN 202410016633A CN 117523025 B CN117523025 B CN 117523025B
- Authority
- CN
- China
- Prior art keywords
- image
- uniformity
- visual field
- transformation matrix
- imitation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000009826 distribution Methods 0.000 title claims abstract description 34
- 230000009466 transformation Effects 0.000 claims abstract description 86
- 239000011159 matrix material Substances 0.000 claims abstract description 82
- 230000000007 visual effect Effects 0.000 claims abstract description 75
- 238000005286 illumination Methods 0.000 claims abstract description 70
- 238000012937 correction Methods 0.000 claims abstract description 53
- 238000013507 mapping Methods 0.000 claims abstract description 44
- 238000003384 imaging method Methods 0.000 claims abstract description 16
- 230000005284 excitation Effects 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 20
- 238000007781 pre-processing Methods 0.000 claims description 9
- 238000003860 storage Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 4
- 230000001678 irradiating effect Effects 0.000 claims description 4
- 239000003086 colorant Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims 4
- 230000003287 optical effect Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000004075 alteration Effects 0.000 description 2
- 238000000799 fluorescence microscopy Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于光强分布调校的数字仿体生成方法，所述方法包括：对投影图像进行图像逆操作以实现相机标定；其中，图像逆操作包括畸变校正、视野均匀度校正和/或照明均匀度校正，分别得到畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵；根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待输入空间光调制器的图像进行校正，经成像投影后形成数字仿体。本发明通过图像逆操作校正，减小数字仿体的畸变，提高了数字仿体的均匀度和光强准确性。
Description
技术领域
本发明涉及生物荧光成像技术领域，尤其涉及一种基于光强分布调校的数字仿体生成方法。
背景技术
在医用光学设备领域，仿体常用于设备的质量检测，是标准化验证的重要方式，在设备的生产制造及临床使用等过程中发挥着重要作用。仿体的重要特性包括稳定性、准确性等，数字仿体以光学调控的方式将光场再现出来，具有较高的稳定性；而对于仿体的准确性，由于光学系统存在像差，再加上荧光分子的发射效率收到照明光强度的影响，将会造成数字仿体的空间、强度与真实情况有所偏差，因此需要使用适当的校正手段，来提高数字仿体的准确性。
畸变是由于光学系统视场放大不均匀所产生的，校正方法包括硬件和软件两个方法，硬件上可以通过镜头设计来实现，需要对整个光学系统进行仿真，较为复杂；软件上可以通过对图像进行变换，如数字图像处理方法包括几何校正、透视变换等，通常用于图像后处理。在生成数字仿体时，由于匀光光束发散角度受限，因此会产生画面强度差异。在荧光成像过程中，激发光会影响荧光发射强度，而普通的数字仿体由于直接发射荧光，无法将荧光激发过程考虑进去，造成数字仿体与实际情况不完成符合。
发明内容
本发明的目的在于针对现有技术的不足，提供了一种基于光强分布调校的数字仿体生成方法。
第一方面，本发明实施例提供了一种基于光强分布调校的数字仿体生成方法，所述方法包括：
对投影图像进行图像逆操作以实现相机标定；其中，图像逆操作包括畸变校正、视野均匀度校正和/或照明均匀度校正，分别得到畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵；
根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待输入空间光调制器的图像进行校正，经成像投影后形成数字仿体。
进一步地，对投影图像进行畸变校正得到畸变坐标映射关系的过程包括：
将标准黑白棋盘格作为原始图像输入至空间光调制器，将原始图像投影出来，形成棋盘格仿体；
通过相机记录棋盘格仿体，得到棋盘格仿体图像；
提取原始图像和棋盘格仿体图像的控制点坐标；
将棋盘格仿体图像的控制点坐标映射到原始图像上，得到畸变坐标映射关系。
进一步地，对投影图像进行视野均匀度校正得到视野均匀度变换矩阵的过程包括：
将一白色图片作为原始图像输入至空间光调制器，将原始图像投影出来，形成白光视野仿体；
通过相机记录下白光视野仿体，得到白光视野仿体图像；
移除白光视野仿体，在均匀环境光的条件下使用相机记录下当前的图像，得到环境光视野图像；
计算白光视野仿体图像与环境光视野图像的比值，得到视野均匀度变换矩阵。
进一步地，对投影图像进行照明均匀度校正得到照明均匀度变换矩阵的过程包括：
当采集荧光图像的时候，使用激发光照射荧光分子，通过相机记录下成像面的照明光能量分布，得到第一激发光能量图像；
再将激发光照射在白屏上，以模拟照明光照射荧光分子，通过相机记录下成像面的照明光能量分布，得到第二激发光能量图像；
计算第一激发光能量图像与第二激发光能量图像的比值，得到照明均匀度变换矩阵。
进一步地，根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待输入空间光调制器的图像进行校正，经成像投影后形成数字仿体包括：
待成像图像的像素、数字仿体的像素与空间光调制器的像素一一对应；
根据畸变坐标映射关系对待成像图像进行畸变校正；根据视野均匀度变换矩阵对待成像图像进行视野均匀度预处理；根据照明均匀度变换矩阵对待成像图像进行照明均匀度校正；
将校正后的待成像图像输入至空间光调制器，形成数字仿体。
进一步地，根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待输入空间光调制器的图像进行校正，经成像投影后形成数字仿体还包括：
依据变换矩阵设置空间光调制器对应的占空比，实现对视野的均匀度和照明的均匀度的修正。
第二方面，本发明实施例提供了一种数字仿体，通过上述的基于光强分布调校的数字仿体生成方法生成。
第三方面，本发明实施例提供了一种双通道数字仿体，由可见光数字仿体和近红外数字仿体进行空间配准，使得可见光数字仿体和近红外数字仿体在空间上重合形成；
其中，可见光数字仿体的形成过程包括：
对投影图像进行图像逆操作以实现相机标定；其中，图像逆操作包括畸变校正、视野均匀度校正和/或照明均匀度校正，分别得到畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵；
采用红、绿、蓝三色作为光源，根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待成像的可见光图像进行校正，将校正后的可见光图像输入至空间光调制器，形成可见光数字仿体；
近红外数字仿体的形成过程包括：
对投影图像进行图像逆操作以实现相机标定；其中，图像逆操作包括畸变校正、视野均匀度校正和/或照明均匀度校正，分别得到畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵；
采用多色近红外LED作为光源，根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待成像的近红外光图像进行校正，将校正后的可见光图像输入至空间光调制器，形成近红外数字仿体。
第四方面，本发明实施例提供了一种电子设备，包括存储器和处理器，所述存储器与所述处理器耦接；其中，所述存储器用于存储程序数据，所述处理器用于执行所述程序数据以实现上述的基于光强分布调校的数字仿体生成方法。
第五方面，本发明实施例提供了一种计算机可读存储介质，其上存储有计算机程序，所述程序被处理器执行时实现上述的基于光强分布调校的数字仿体生成方法。
与现有技术相比，本发明的有益效果为：
本发明提供了一种基于光强分布调校的数字仿体生成方法，通过对投影图像进行图像逆操作以实现相机标定，得到畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵，并根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待输入空间光调制器的图像进行校正，修改空间光调制器的像素分布，从而实现了对数字仿体的修正，提高了数字仿体的准确性，且视场角不影响荧光强度。本发明通过图像预处理以及空间光调制器的联合控制，简便有效，调节灵活。本发明方法可以实现多个数字仿体的空间配准，实现多通道数字仿体，解决投影遮挡问题。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种基于光强分布调校的数字仿体生成方法的流程图；
图2是本发明实施例提供的数字仿体光强分布调校的示意图；
图3是本发明实施例提供的双通道数字仿体匹配的示意图；
图4是本发明实施例提供的一种电子设备。
图中，1-空间光调制器；101-可见光空间光调制器；102-近红外光空间光调制器；2-数字仿体；3-相机。
具体实施方式
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本发明的一些方面相一致的装置和方法的例子。
在本发明使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本发明。在本发明和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解，尽管在本发明可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本发明范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
下面结合附图，对本发明进行详细说明。在不冲突的情况下，下述的实施例及实施方式中的特征可以相互组合。
如图1所示，本发明实施例提供了一种数字仿体光强分布调校方法，所述方法包括以下步骤：
步骤S1，将投影图像进行图像逆操作以实现相机3的标定，所述图像逆操作包括畸变校正、视野均匀度校正和/或照明均匀度校正，分别得到畸变坐标映射关系、视野均匀度变换矩阵、照明均匀度变换矩阵；
其中，如图2所示，畸变校正的目的为校正光学系统、投影角度带来的像差，本发明对投影图像进行畸变校正得到畸变坐标映射关系，具体过程包括：
将标准黑白棋盘格作为原始图像输入至空间光调制器，将原始图像投影出来，形成棋盘格仿体；
通过相机记录棋盘格仿体，得到棋盘格仿体图像；
提取原始图像和棋盘格仿体图像的控制点坐标；
将棋盘格仿体图像的控制点坐标映射到原始图像上，得到畸变坐标映射关系。
数字仿体通过镜头投影的方式产生荧光，镜头中央视野的光通量大于边缘视野，由此数字仿体的强度分布也受到影响。为消除不同视场角度的光强差异，本实例对投影图像进行视野均匀度校正得到视野均匀度变换矩阵，具体过程包括：
将一白色图片作为原始图像输入至空间光调制器，将原始图像投影出来，形成白光视野仿体；
通过相机记录下白光视野仿体，得到白光视野仿体图像；
移除白光视野仿体，在均匀环境光的条件下使用相机记录下当前的图像，得到环境光视野图像；
计算白光视野仿体图像与环境光视野图像的比值，得到视野均匀度变换矩阵。
明均匀度校正的主要目的是将激发光强度纳入荧光仿体发射强度的考虑范围，即激发光强度是荧光强的重要影响因素，需要依据激发光的强弱来调节数字仿体的强度。本实例中对投影图像进行照明均匀度校正得到照明均匀度变换矩阵的过程包括：
当采集荧光图像的时候，使用激发光照射荧光分子，通过相机记录下成像面的照明光能量分布，得到第一激发光能量图像；
再将激发光照射在白屏上，以模拟照明光照射荧光分子，通过相机记录下成像面的照明光能量分布，得到第二激发光能量图像；
计算第一激发光能量图像与第二激发光能量图像的比值，得到照明均匀度变换矩阵。
需要说明的是，本实例涉及投影图像进行预处理，即图像的逆操作。假设真实荧光分布图像为A，经过光学系统生成的仿体与实际荧光相比存在一定的误差，记为σ，此时仿体的表达式为Aσ。求解σ的逆变换σ-1，并对仿体进行预处理，得到仿体Aσ-1σ=A，可以看出，通过逆向操作的预处理，结合空间光调制器的阵列特征，仿体的逆操作被编码在空间光处理器上，使仿体的光强分布与实际荧光一致。具体来说，仿体图像预处理有多种形式，包括常见的三种：畸变校正，视野均匀度校正和照明均匀度校正。由于逆向计算要有原始仿体，因此将相机放置于可成像位置，校正相机本身的畸变和视野均匀度，得到已标定的相机。
步骤S2，根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待输入空间光调制器1的图像进行校正，经镜头成像投影后形成数字仿体2。
进一步地，待成像图像的像素、数字仿体的像素与空间光调制器的像素一一对应；在本实例中，根据畸变坐标映射关系对待成像图像进行畸变校正；根据视野均匀度变换矩阵对待成像图像进行视野均匀度预处理；根据照明均匀度变换矩阵对待成像图像进行照明均匀度校正；将校正后的待成像图像输入至空间光调制器，形成可见光数字仿体。
进一步地，所述步骤S2还包括：设置空间光调制器对应的占空比，实现对视野的均匀度和照明的均匀度的修正。
需要说明的是，根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待输入空间光调制器的图像进行校正预处理后，空间光调制器的输入图像发生了投影变换，当把校正后的可见光图像输入至空间光调制器，空间光调制器的图像也产生了相同的投影变换，即原来空间光调制器的像素翻转情况发生了若干像素的偏移，从而导致数字仿体的图案也发生了若干像素的偏移，用于补偿镜头本身、投射角度等因素造成的数字仿体畸形。此外，空间光调制器像素翻转时间将影响相应位置的光强度，因此可以通过修改占空比，对视野的均匀度和照明的均匀度作修正。
实施例1
本实施例提供了基于上述的基于光强分布调校的数字仿体生成方法生成了一种双通道数字仿体。
具体地，如图3所示，所述双通道数字仿体包括可见光数字仿体和近红外数字仿体。
其中，可见光数字仿体的形成过程包括：
对投影图像进行图像逆操作以实现相机标定；其中，图像逆操作包括畸变校正、视野均匀度校正和/或照明均匀度校正，分别得到畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵；
采用红、绿、蓝三色作为光源，根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待成像的可见光图像进行校正，将校正后的可见光图像输入至可见光空间光调制器101，形成可见光数字仿体；
近红外数字仿体的形成过程包括：
对投影图像进行图像逆操作以实现相机标定；其中，图像逆操作包括畸变校正、视野均匀度校正和/或照明均匀度校正，分别得到畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵；
采用多色近红外LED作为光源，根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待成像的近红外光图像进行校正，将校正后的可见光图像输入至近红外光空间光调制器102，形成近红外数字仿体。
最后，将可见光部分数字仿体和近红外数字仿体进行空间配准，即平移操作，将可见光部分数字仿体和近红外数字仿体在空间上重合，得到双通道数字仿体。
综上所述，本发明提供了一种基于光强分布调校的数字仿体生成方法，通过对投影图像进行图像逆操作以实现相机标定，得到畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵，并根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待输入空间光调制器的图像进行校正，修改空间光调制器的像素分布，从而实现了对数字仿体的修正，提高了数字仿体的准确性，且视场角不影响荧光强度。本发明通过图像预处理以及空间光调制器的联合控制，简便有效，调节灵活。本发明方法可以实现多个数字仿体的空间配准，实现多通道数字仿体，解决投影遮挡问题。
如图4所示，本申请实施例提供一种电子设备，其包括存储器，用于存储一个或多个程序；处理器。当一个或多个程序被处理器执行时，实现如上述第一方面中任一项的方法。
还包括通信接口，该存储器、处理器和通信接口相互之间直接或间接地电性连接，以实现数据的传输或交互。例如，这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。存储器可用于存储软件程序及模块，处理器通过执行存储在存储器内的软件程序及模块，从而执行各种功能应用以及数据处理。该通信接口可用于与其他节点设备进行信令或数据的通信。
其中，存储器可以是但不限于，随机存取存储器（Random Access Memory，RAM），只读存储器（Read Only Memory，ROM），可编程只读存储器（Programmable Read-OnlyMemory，PROM），可擦除只读存储器（Erasable Programmable Read-Only Memory，EPROM），电可擦除只读存储器（Electric Erasable Programmable Read-Only Memory，EEPROM）等。
处理器可以是一种集成电路芯片，具有信号处理能力。该处理器可以是通用处理器，包括中央处理器（Central Processing Unit，CPU）、网络处理器（Network Processor，NP）等；还可以是数字信号处理器（Digital Signal Processing，DSP）、专用集成电路（Application Specific Integrated Circuit，ASIC）、现场可编程门阵列（Field－Programmable Gate Array，FPGA）或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
在本申请所提供的实施例中，应该理解到，所揭露的方法及系统，也可以通过其它的方式实现。以上所描述的方法及系统实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本申请的多个实施例的方法及系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
另外，在本申请各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。
另一方面，本申请实施例提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上述第一方面中任一项的方法。所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器（ROM，Read-Only Memory）、随机存取存储器（RAM，Random AccessMemory）、磁碟或者光盘等各种可以存储程序代码的介质。
本领域技术人员在考虑说明书及实践这里公开的内容后，将容易想到本申请的其它实施方案。本申请旨在涵盖本申请的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本申请的一般性原理并包括本申请未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的。
应当理解的是，本申请并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。
Claims (7)
1.一种基于光强分布调校的数字仿体生成方法，其特征在于，所述方法包括：
对投影图像进行图像逆操作以实现相机标定；其中，图像逆操作包括畸变校正、视野均匀度校正和/或照明均匀度校正，分别得到畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵；
根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待输入空间光调制器的图像进行校正，经成像投影后形成数字仿体；
其中，对投影图像进行畸变校正得到畸变坐标映射关系的过程包括：
将标准黑白棋盘格作为原始图像输入至空间光调制器，将原始图像投影出来，形成棋盘格仿体；
通过相机记录棋盘格仿体，得到棋盘格仿体图像；
提取原始图像和棋盘格仿体图像的控制点坐标；
将棋盘格仿体图像的控制点坐标映射到原始图像上，得到畸变坐标映射关系；
其中，对投影图像进行视野均匀度校正得到视野均匀度变换矩阵的过程包括：
将一白色图片作为原始图像输入至空间光调制器，将原始图像投影出来，形成白光视野仿体；
通过相机记录下白光视野仿体，得到白光视野仿体图像；
移除白光视野仿体，在均匀环境光的条件下使用相机记录下当前的图像，得到环境光视野图像；
计算白光视野仿体图像与环境光视野图像的比值，得到视野均匀度变换矩阵；
其中，对投影图像进行照明均匀度校正得到照明均匀度变换矩阵的过程包括：
当采集荧光图像的时候，使用激发光照射荧光分子，通过相机记录下成像面的照明光能量分布，得到第一激发光能量图像；
再将激发光照射在白屏上，以模拟照明光照射荧光分子，通过相机记录下成像面的照明光能量分布，得到第二激发光能量图像；
计算第一激发光能量图像与第二激发光能量图像的比值，得到照明均匀度变换矩阵。
2.根据权利要求1所述的基于光强分布调校的数字仿体生成方法，其特征在于，根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待输入空间光调制器的图像进行校正，经成像投影后形成数字仿体包括：
待成像图像的像素、数字仿体的像素与空间光调制器的像素一一对应；
根据畸变坐标映射关系对待成像图像进行畸变校正；根据视野均匀度变换矩阵对待成像图像进行视野均匀度预处理；根据照明均匀度变换矩阵对待成像图像进行照明均匀度校正；
将校正后的待成像图像输入至空间光调制器，形成数字仿体。
3.根据权利要求1所述的基于光强分布调校的数字仿体生成方法，其特征在于，根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待输入空间光调制器的图像进行校正，经成像投影后形成数字仿体还包括：
依据变换矩阵设置空间光调制器对应的占空比，实现对视野的均匀度和照明的均匀度的修正。
4.一种数字仿体，其特征在于，通过权利要求1-3任一项所述的基于光强分布调校的数字仿体生成方法生成。
5.一种双通道数字仿体，其特征在于，由可见光数字仿体和近红外数字仿体进行空间配准，使得可见光数字仿体和近红外数字仿体在空间上重合形成；
其中，可见光数字仿体的形成过程包括：
对投影图像进行图像逆操作以实现相机标定；其中，图像逆操作包括畸变校正、视野均匀度校正和/或照明均匀度校正，分别得到畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵；
采用红、绿、蓝三色作为光源，根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待成像的可见光图像进行校正，将校正后的可见光图像输入至空间光调制器，形成可见光数字仿体；
近红外数字仿体的形成过程包括：
对投影图像进行图像逆操作以实现相机标定；其中，图像逆操作包括畸变校正、视野均匀度校正和/或照明均匀度校正，分别得到畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵；
采用多色近红外LED作为光源，根据畸变坐标映射关系、视野均匀度变换矩阵和/或照明均匀度变换矩阵对待成像的近红外光图像进行校正，将校正后的可见光图像输入至空间光调制器，形成近红外数字仿体；
其中，对投影图像进行畸变校正得到畸变坐标映射关系的过程包括：
将标准黑白棋盘格作为原始图像输入至空间光调制器，将原始图像投影出来，形成棋盘格仿体；
通过相机记录棋盘格仿体，得到棋盘格仿体图像；
提取原始图像和棋盘格仿体图像的控制点坐标；
将棋盘格仿体图像的控制点坐标映射到原始图像上，得到畸变坐标映射关系；
其中，对投影图像进行视野均匀度校正得到视野均匀度变换矩阵的过程包括：
将一白色图片作为原始图像输入至空间光调制器，将原始图像投影出来，形成白光视野仿体；
通过相机记录下白光视野仿体，得到白光视野仿体图像；
移除白光视野仿体，在均匀环境光的条件下使用相机记录下当前的图像，得到环境光视野图像；
计算白光视野仿体图像与环境光视野图像的比值，得到视野均匀度变换矩阵；
其中，对投影图像进行照明均匀度校正得到照明均匀度变换矩阵的过程包括：
当采集荧光图像的时候，使用激发光照射荧光分子，通过相机记录下成像面的照明光能量分布，得到第一激发光能量图像；
再将激发光照射在白屏上，以模拟照明光照射荧光分子，通过相机记录下成像面的照明光能量分布，得到第二激发光能量图像；
计算第一激发光能量图像与第二激发光能量图像的比值，得到照明均匀度变换矩阵。
6.一种电子设备，包括存储器和处理器，其特征在于，所述存储器与所述处理器耦接；其中，所述存储器用于存储程序数据，所述处理器用于执行所述程序数据以实现上述权利要求1-3任一项所述的基于光强分布调校的数字仿体生成方法。
7.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述程序被处理器执行时实现如权利要求1-3任一项所述的基于光强分布调校的数字仿体生成方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410016633.5A CN117523025B (zh) | 2024-01-05 | 2024-01-05 | 一种基于光强分布调校的数字仿体生成方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410016633.5A CN117523025B (zh) | 2024-01-05 | 2024-01-05 | 一种基于光强分布调校的数字仿体生成方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117523025A CN117523025A (zh) | 2024-02-06 |
| CN117523025B true CN117523025B (zh) | 2024-05-10 |
Family
ID=89766886
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410016633.5A Active CN117523025B (zh) | 2024-01-05 | 2024-01-05 | 一种基于光强分布调校的数字仿体生成方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117523025B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108765301A (zh) * | 2018-03-27 | 2018-11-06 | 长春理工大学 | 一种光学系统及其畸变校正方法及系统 |
| CN109461126A (zh) * | 2018-10-16 | 2019-03-12 | 重庆金山医疗器械有限公司 | 一种图像畸变校正方法及系统 |
| CN110149511A (zh) * | 2019-05-13 | 2019-08-20 | 北京理工大学 | 一种畸变校正方法、装置及显示系统 |
| CN114881887A (zh) * | 2022-06-06 | 2022-08-09 | 中国传媒大学 | 一种大尺度异形面投影几何校正方法 |
| CN116124752A (zh) * | 2023-04-04 | 2023-05-16 | 之江实验室 | 一种基于多光谱调控的组织仿生模体及其生成方法 |
| CN116543055A (zh) * | 2023-03-02 | 2023-08-04 | 西安全志科技有限公司 | 图像畸变标定方法、校正方法、计算机装置、计算机可读存储介质 |
| CN116866537A (zh) * | 2023-06-21 | 2023-10-10 | 中国海洋大学 | 一种基于角点检测的投影图像畸变校正方法及投影系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3822687A1 (en) * | 2019-11-15 | 2021-05-19 | Leica Microsystems CMS GmbH | Optical imaging device for a microscope |
-
2024
- 2024-01-05 CN CN202410016633.5A patent/CN117523025B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108765301A (zh) * | 2018-03-27 | 2018-11-06 | 长春理工大学 | 一种光学系统及其畸变校正方法及系统 |
| CN109461126A (zh) * | 2018-10-16 | 2019-03-12 | 重庆金山医疗器械有限公司 | 一种图像畸变校正方法及系统 |
| CN110149511A (zh) * | 2019-05-13 | 2019-08-20 | 北京理工大学 | 一种畸变校正方法、装置及显示系统 |
| CN114881887A (zh) * | 2022-06-06 | 2022-08-09 | 中国传媒大学 | 一种大尺度异形面投影几何校正方法 |
| CN116543055A (zh) * | 2023-03-02 | 2023-08-04 | 西安全志科技有限公司 | 图像畸变标定方法、校正方法、计算机装置、计算机可读存储介质 |
| CN116124752A (zh) * | 2023-04-04 | 2023-05-16 | 之江实验室 | 一种基于多光谱调控的组织仿生模体及其生成方法 |
| CN116866537A (zh) * | 2023-06-21 | 2023-10-10 | 中国海洋大学 | 一种基于角点检测的投影图像畸变校正方法及投影系统 |
Non-Patent Citations (1)
| Title |
|---|
| 基于迭代校正的鱼眼棋盘格图像角点检测;申夏晶;程梦娇;肖江剑;宋文祥;;激光与光电子学进展;20160523(第09期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117523025A (zh) | 2024-02-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11210839B2 (en) | Photometric image processing | |
| US11244436B2 (en) | Method for inspecting mounting state of component, printed circuit board inspection apparatus, and computer readable recording medium | |
| CN105814887A (zh) | 图像生成装置 | |
| JP6484706B2 (ja) | 画像を記録するための装置および方法 | |
| JP6783895B2 (ja) | 反射を低減したコントラスト画像を生成するための方法およびこれに関連する装置 | |
| US9105092B2 (en) | Registration method of images for surgery | |
| CN112848281B (zh) | 一种光固化3d打印机光补偿方法 | |
| CN105466889B (zh) | 一种空间频域成像中复杂组织体表面光照度的采集方法 | |
| US11714278B2 (en) | Method and system for calibrating a wearable heads-up display to produce aligned virtual images in an eye space | |
| US20210211620A1 (en) | Method and apparatus for correcting color convergence error, and device | |
| CN112469978A (zh) | 颜色评价装置、颜色评价方法和颜色评价方法中使用的显示物 | |
| CN117523025B (zh) | 一种基于光强分布调校的数字仿体生成方法 | |
| CN113920037B (zh) | 一种内窥镜画面矫正方法、装置、矫正系统及存储介质 | |
| JP2019117120A (ja) | 画像処理システム、画像処理方法、及びそのプログラム | |
| CN115100037A (zh) | 基于多线扫相机图像拼接的大幅面瓷砖成像方法及系统 | |
| CN113840125A (zh) | 投影仪的控制方法、装置、投影仪、存储介质 | |
| CN103905738B (zh) | 高动态范围图像生成系统及方法 | |
| Maier-Hein et al. | Accounting for anisotropic noise in fine registration of time-of-flight range data with high-resolution surface data | |
| CN117156289A (zh) | 色彩风格校正方法、系统、电子设备、存储介质及芯片 | |
| CN111105365A (zh) | 一种纹理影像的色彩校正方法、介质、终端和装置 | |
| CN108377383B (zh) | 一种多投影3d系统光场对比度调整方法及其系统 | |
| CN109269404B (zh) | 图像处理方法、装置及指纹识别设备 | |
| WO2021053852A1 (ja) | 外観検査装置、外観検査装置の較正方法及びプログラム | |
| CN109813533B (zh) | 一种批量测试doe衍射效率和均匀性的方法及装置 | |
| JP7190701B2 (ja) | 投写画像調整システム及び方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |