具体实施方式

以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本申请实施例。然而，本领域的技术人员应当清楚，在没有这些具体细节的其它实施例中也可以实现本申请。在其它情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本申请的描述。另外，在本申请说明书和所附权利要求书的描述中，术语“第一”、“第二”、“第三”等仅用于区分描述，而不能理解为指示或暗示相对重要性。

本申请提出一种旋转控制模块、方法与系统，目的是为了实现操控线虫等微小物体在指定位置旋转，以便对该物体进行多角度的观测。本申请主要利用超声波激励腔体容器内的结构体(可以是一个声子晶体板)发生共振，利用共振产生的周期性分布的强声场对处于该腔体容器内的液体中的微小物体进行捕获和排列，同时利用声场在物体周围诱发的非对称分布的声流产生的力矩旋转该物体，相关的具体实施方式请参见下文。

如图1所示，是本申请实施例提供的一种旋转控制系统的示意图，该系统包括三个功能模块，分别为旋转控制模块、成像模块和三维重建模块。其中，旋转控制模块用于操控待观测物体改变姿态，例如进行移动和旋转；成像模块用于拍摄待观测物体在不同姿态下的图像；三维重建模块用于根据拍摄到的图像重建出待观测物体的三维模型。

图1中旋转控制模块的示意图如图2所示，主要包括超声激励模块和腔体容器两部分，腔体容器可以是一个指定形状(例如矩形或方形等)的腔体结构的容器，其尺寸可以根据待观测物体的尺寸合理设置，其材料可以采用玻璃或者聚二甲基硅氧烷等。该腔体容器主要用于容纳含有待观测的目标物体的液体，在操作时，将包含目标物体的液体充入该腔体容器中。另外，往该腔体容器内放入一个结构体，该结构体的尺寸和形状可以根据实际需求设置，该结构体的材料可以采用不锈钢、黄铜或者硅等。示例性的，该结构体可以是一个声子晶体板，即弹性常数及密度周期分布的材料或结构。

超声激励模块用于向腔体容器发射超声波，激励处于该腔体容器内部的结构体产生共振。共振现象会使得该结构体的表面产生声场，通过该声场能够控制目标物体处于该结构体的表面旋转。在实际操作中，超声激励模块可以由信号发生器和功率放大器构成，其中信号发生器用于输出指定频率(基于该结构体的共振频率设置)的脉冲信号，功率放大器用于对该脉冲信号进行放大，得到放大的激励信号，然后可以通过压电片(或者其它超声换能器，如高斯束声源等)将该激励信号转换为超声波，通过产生的超声波激励该结构体产生共振。

在该旋转控制模块的一种实现方式中，腔体容器可以由石英基底、容器壁和可拆卸的顶盖构成，且在该石英基底上粘接一个压电片。其中，该容器壁可以与石英基底和顶盖键合。通过这样设置，一方面，该压电片产生的超声波能够直接作用到结构体上，产生足够强度的声场；另一方面，由于该腔体容器的顶盖是可拆卸的，因此可以方便地实现开放空间下的微注射等显微操作。

在该旋转控制模块的一种实现方式中，结构体为矩形的板状结构，且所述板状结构上设有多列等间距排列的栅格。如图3所示，为本申请实施例提供的一种板状结构体的示意图。图3示出的是该板状结构体的横截面，其中t为板厚，a为两个栅格之间的距离，w为栅格的宽度，h为栅格的高度。根据腔体容器和目标物体的尺寸，可以合理地设置该板状结构的长度、宽度以及图3中的各个尺寸参数。例如，若目标物体为线虫，通常成年线虫的身体直径约为50μm，体长约1mm，则图3中该板状结构体的各个尺寸参数可以设为t＝60μm，h＝40μm，w＝70μm以及a＝300μm。采用图3这种结构可以激发结构体中的非泄露兰姆波，以产生局域化(亚波长)的强透射声场，从而产生强声辐射力和声流。需说明的是，图3只是该板状结构体的一种表现形式，按照需求也可以将该板状结构体设置为弯曲的，且其栅格也可以是曲线形式的或者非等间距的，本申请不对该结构体的具体形状结构进行限定。

如图4所示，为图3所示结构体(由不锈钢材料制作)的透射谱。通过该透射谱可以获知，该结构体的共振频率约为4.6MHz，也即采用频率约为4.6MHz的超声波激励该结构体能够使其产生共振现象，从而生成强度足够大的声场，如此，才能获得较强的声辐射力和声流来操控目标物体运动。

图1中的成像模块可以由明场光源、光学透镜和相机等部件组成，主要用于拍摄目标物体在结构体的表面旋转的视频。另外，针对线虫体内组织器官或者细胞三维形貌重建的应用场合，还可以为该成像模块添加荧光光源。由于目标物体在声场的作用下自行旋转，故在视频拍摄过程中无需移动相机。

图1中的三维重建模块可以使用计算机等设备，通过采用相机标定算法和三维重建算法，在该计算机上运行图像处理软件以及纹理映射和渲染软件，基于成像模块采集到的目标物体的旋转视频，可以提取目标物体在多个不同旋转角度下的图像，进而根据这些图像重建出该目标物体对应的三维模型。

关于图1所示系统的工作原理和具体操作步骤，可以参见下文所述的方法实施例。

如图5所示，是本申请实施例提供的一种旋转控制方法的流程图，该方法应用于前文所述的旋转控制模块或系统。该方法包括：

501、使用超声激励模块向腔体容器发射超声波；

针对前文所述的旋转控制模块，首先，往腔体容器中充入包含目标物体的液体，具体可以使用移液枪等设备将混有微小物体(一般为柱形或者近似柱形的颗粒，例如线虫，作为待观测的目标物体)的液体(例如水)移入该腔体容器内。该腔体容器内还放置一个结构体(例如可以采用图3所示的板状结构体)，可以使用移液枪将该液体直接注射到该结构体上，在后续过程中液体可以作为声传播介质。

在操作时，使用超声激励模块向腔体容器发射超声波，通过超声波可以激励腔体容器内的结构体产生共振，从而在该结构体的表面形成声场。在本申请的一种实现方式中，腔体容器的底部设有压电片，超声激励模块包括信号发生器和功率放大器，使用超声激励模块向腔体容器发射超声波，可以包括：

(1)使用信号发生器输出指定频率的脉冲信号，所述指定频率根据结构体的共振频率设定；

(2)使用功率放大器对脉冲信号进行放大，得到激励信号，以使压电片在激励信号的作用下产生超声波。

超声激励模块可以由信号发生器和功率放大器构成，其中信号发生器用于产生指定频率(基于结构体的共振频率以及压电片的中心频率设置)的脉冲信号；例如，若该结构体的共振频率为4.6MHz，而该压电片的中心频率为4.5MHz，则可以采用信号发生器输出频率为4.577MHz的脉冲信号。该脉冲信号通过功率放大器放大，得到激励信号，用于激励压电片(作为超声换能器)产生超声波。

502、通过超声波激励结构体发生共振，以使结构体的表面产生声场，并通过声场控制目标物体于结构体的表面旋转。

通过超声波激励结构体产生共振，会在该结构体的表面形成足够强度的周期性的声场，进而产生周期性的声辐射力和声流。在声辐射力和声流的作用下，能够控制该目标物体于该结构体的表面旋转，以便于显微拍摄该目标物体在多个不同旋转角度下的图像。通过声场控制目标物体于该结构体的表面旋转的具体原理，请参见下文所述。

在本申请的一种实现方式中，通过所述声场控制所述目标物体于所述结构体的表面旋转，可以包括：

(1)通过所述声场的声辐射力将漂浮于所述液体中的所述目标物体吸附至所述结构体的表面；

(2)通过所述声场的声流控制所述目标物体旋转。

在声辐射力的作用下，能够将漂浮于液体中的目标物体吸附至该结构体的表面。而在声流的作用下，能够控制该目标物体旋转。以图3所示的板状结构体为例，使用多物理场仿真软件(例如COMSOL)可以计算模拟该板状结构表面产生的声场，包括声压场和速度场等，还可以进一步计算目标物体在声场中受到的声辐射力，根据软件模拟的结果可以解释目标物体移动和旋转的原理。

考虑到结构体的周期性和对称性，使用二维计算模型，计算区域包括液体(一般是水)和结构体，其中液体由压力声学(频域)模拟，结构体由固体力学模型模拟。在结构体的底部区域设置平面波入射(模拟压电片产生的超声波)，在结构体的顶部区域设置平面波辐射，假设超声波无反射，通过多物理场仿真软件模拟得到的该结构体表面的声压分布示意图如图6(a)所示。在图6(a)中，示出了结构体两个栅格表面附近的声压分布，其中两个白色小圈标示出目标物体(例如线虫)的捕获位置，可以视作目标物体的横截面，在计算中近似视为圆柱形。另外，通过多物理场仿真软件可以模拟得到图6(a)中该结构体表面的不同位置受到的声辐射力的示意图，如图6(b)所示。在图6(b)中，Fx为在x方向上受到的二维声辐射力分量，Fy为在y方向上受到的二维声辐射力分量，x方向和y方向可以参照图6(a)所标示的方向。可见，图6(b)中Fx为0的位置，是目标物体在x方向的平衡位置(图中的灰色圆圈)，由于在该平衡位置两侧的声辐射力均指向该平衡位置，因此会迫使该目标物体向该平衡位置移动。另外，由于在Fx为0的位置，对应的Fy为负值，即y方向上的声辐射力指向该结构体的表面，因此，漂浮的目标物体会被吸附至该结构体的表面。从结构体的上方朝下看，可以观察到目标物体以一定的间距排列分布。

接下来说明通过声场控制目标物体旋转的原理。声场中的粒子所受的辐射应力可以用以下公式表示：

其中，σ表示辐射应力，〈*〉表示时间平均算子，I表示单位张量，为并矢符号；ρ₀和c₀分别为声传播介质的密度和声速，室温下是介质的固定参数，可通过查询资料获得；p和v分别为一阶声压和速度场。

根据波动方程和动量守恒方程，可以获得以下2个公式：

其中，t表示在时间维度下求偏导。基于上述3个公式，采用有限元法和有限差分法等数值求解方法可以求解得到辐射应力σ。

在计算出声场中的粒子所受的辐射应力之后，可以采用以下公式分别计算目标物体在声场中受到的辐射力和辐射力矩：

F_rad＝-∫∫_SσdS

T_rad＝-∫∫_Sr×σdS

其中，F_rad表示辐射力，T_rad表示辐射力矩，r表示目标物体的质心(在二维计算模型下，质心可以为目标物体横截面的中心)到目标物体表面一点的方向向量，dS为该表面一点的法向向量和面积微元的乘积。

经过模拟的数值计算，发现目标物体受到的辐射力矩接近于0，而在目前的理论体系中，认为声场中的力矩来源于辐射力和声流，因此若辐射力矩接近于0，则可以判定目标物体的旋转是由声流引起的，以下计算声场在结构体表面产生的声流分布。

不可压缩流体的粘性应力(黏性应力)可以表示为：

τ_ij＝μ(u_i,j+u_j,i)

其中，τ_ij表示粘性应力，μ为流体动力粘度，在室温下是介质的固定参数，可通过查阅资料获得；u表示声流速度；下标i和下标j的取值可以为1、2或者3，分别表示变量在空间中沿着三个坐标轴x，y和z的分量，i和j之间的逗号表示变量对坐标求偏导。例如，x方向和y方向之间的粘性应力可以表示为：

τ₁₂＝μ(u_1,2+u_2,1)

其中，u_1,2表示u_x在y方向上的偏导，u_2,1表示u_y在x方向上的偏导，u_x表示声流速度在x方向上的分量，u_y表示声流速度在y方向上的分量。

具体的，声流速度u可以采用以下公式计算得到：

其中，p₂为二阶压力场，〈*〉为时间平均算子，为并矢符号，ρ₀为声传播介质的密度，μ为流体动力粘度，v为速度场，此公式可使用软件COMSOL中的蠕动流模块求解，在求解出声流速度u后即可计算得到粘性应力，根据粘性应力可得目标物体所受的声粘性扭矩为：

T_vis＝∫∫_Sr×τdS

由于捕获的目标物体两侧的声流是非对称分布的，即一侧强一侧弱，这导致两侧声流在目标物体表面产生的力矩虽然方向相反，但是大小不同，从而使得声粘性扭矩不为零，因此会使目标物体发生旋转。

图7是目标物体在声流作用下进行旋转的示意图，其中图7左方的目标物体1(圆形表示目标物体的横截面)沿着自身中心顺时针旋转，图7右方的目标物体2沿着自身中心逆时针旋转。图7中还用箭头标示出目标物体附近的声流和目标物体所受的力矩，可以看出对于左方的目标物体1，其左侧的声流较强，而右侧的声流较弱，因此会产生顺时针方向的力矩，使得目标物体顺时针旋转，对于右方的目标物体则结论相反。

在本申请的一种实现方式中，在通过所述声场控制所述目标物体于所述结构体的表面旋转之后，还可以包括：

(1)拍摄所述目标物体在所述结构体的表面旋转的视频；

(2)从所述视频中提取所述目标物体在多个不同旋转角度下的图像，并根据所述目标物体在多个不同旋转角度下的图像重建得到所述目标物体的三维模型。

为了重建该目标物体的三维模型，可以使用相机拍摄该目标物体在结构体的表面旋转的视频，然后使用计算机从该视频中提取出目标物体在多个不同旋转角度(例如0°、90°、180°、270°和360°)的图像，这个过程也称为拆帧，之后可以对图像进行裁剪、增强、滤波或者边缘识别等预处理过程；接着，使用三维重建算法(例如Marching Cube算法、BallPivoting算法和Screened Poisson算法等)重建得到目标物体的三维形貌，最后再进行纹理映射和渲染处理，得到对应的三维模型，相应的模型重建流程示意图如图8所示。

另一方面，在重建目标物体的三维模型之前，需要对相机参数进行标定和求解，即需要获得相机矩阵P，以描述三维空间中一点在二维图像中的像素点坐标(x，y)和实际空间坐标(X，Y，Z)的对应关系。相机矩阵的具体形式可以表示为P＝K[R|t]，其中K为相机的内参矩阵，可用现有的相机标定方法获得，[R|t]为外参矩阵，表示世界坐标系和相机坐标系的变换矩阵，其中R为旋转矩阵，可根据目标物体的旋转角度求得，t为平移矩阵，可根据世界坐标系原点和相机坐标系原点的位置求得。在获得相机矩阵后，结合图像处理获得的不同旋转角度的目标物体的图像坐标，可以计算出不同旋转角度的目标物体在世界坐标系下的坐标，之后再执行三维模型的重建流程。

本申请实施例采用超声激励结构体发生共振，利用共振产生的空间分布的局域强声场，会在该结构体的表面形成声辐射力和声流，在声辐射力和声流的作用下，可以控制该目标物体于该结构体的表面旋转，以便对该目标物体进行多角度的观测。而且，可以采用具有可拆卸顶盖的腔体容器，从而实现开放空间下的微注射等显微操作。另外，在控制目标物体旋转后，还可以通过显微镜相机拍摄目标物体不同旋转角度的图像，采用三维重建算法对这些图像进行重建，获得对应的三维模型。

应理解，上述各个实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。

为便于理解本申请提出的技术方案，以下列举一个实际的应用场景。首先，采用C304不锈钢，基于标准化学刻蚀工艺，制作出如图3所示的带有周期性栅格的结构体，具体尺寸参数为t＝60μm，h＝40μm，w＝70μm，a＝300μm。将该结构体放入一个腔体容器内，该腔体容器由石英玻璃基底、聚二甲基硅氧烷材料的容器壁和玻璃顶盖构成，容器壁可以与基底和顶盖键合。超声换能器采用中心频率为4.5MHz的PZT4型压电陶瓷片，其通过环氧树脂与腔体容器的玻璃基底粘接在一起。采用信号发生器产生频率为4.577MHz的脉冲信号，经由功率放大器放大后，激励压电陶瓷片产生超声波。结构体在超声波的作用下产生共振，使得该结构体的表面形成较强的声场。

在操作时，首先打开腔体容器的玻璃顶盖，使用移液枪将混有线虫和玻璃圆棒的水注射到该结构体的表面，然后盖上玻璃顶盖。线虫和玻璃圆棒在声场的作用下，会移动至该结构体的表面旋转，具体原理可参照前文所述。其中，玻璃圆棒用于辅助验证计算结论，证明本申请提出的方案可以用于观测任意圆柱形或者近似圆柱形的微小物体。将整个腔体容器放置于显微镜下，可以观察到线虫和玻璃圆棒都沿着结构体的表面竖直排列，如图9所示。

接下来，用显微镜观测其中一条线虫，拍摄该线虫进行旋转的视频，对视频拆帧后可以得到该线虫在多个不同旋转角度下的图像，之后即可按照图8所示的流程重建得到该线虫的三维模型。获得的线虫在多个不同旋转角度下的图像以及重建得到的三维模型的示意图如图10所示。

以上所述实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围，均应包含在本申请的保护范围之内。