Name: 通过开源数字图像相关 测量 肌腱中的局部组织应变 (视频) | JoVE Journal
Uploaded: 2023-01-27T14:30:00.000+00:00
Duration: 7 min 50 s
Description: 了解开源数字图像相关算法如何测量肌腱局部组织应变。该技术使用ALDIC方法分析Achilles肌腱样本，验证测量准确性并提供反馈，助力生物力学与组织工程研究。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

本工作由美国国立卫生研究院（R21 AR079095）和美国国家科学基金会（2142627）资助。

本研究已获得宾夕法尼亚州立大学机构动物护理和使用委员会的批准。
1. 组织制备
<ol>
	<li>本方案中，取2至4月龄雄性C57BL/6小鼠的跟腱。 
	注意：也可使用小鼠或其他小型动物的不同肌腱或韧带。
	<ol>
		<li>在跟腱表面的皮肤上做切口，暴露跖肌腱及其周围结缔组织，然后用手术刀将其去除。</li>
		<li>将暴露的比目鱼肌和腓肠肌从后肢分离，并用手术刀小心地从跟腱上刮除。</li>
		<li>使用旋转工具上的切割轮附件将跟骨与足部其余部分分离。</li>
	</ol>
	</li>
	<li>将组织置于1.5 mL含5 µg/mL的5-(4,6-二氯三嗪基)氨基荧光素（DTAF）和0.1 M碳酸氢钠缓冲液的溶液中，在室温下于旋转混合器上染色20分钟。该溶液可对组织中的蛋白质（如细胞外基质）进行染色。 
	注意：在上述20分钟染色期间，应完成步骤1.3。</li>
	<li>用磷酸盐缓冲液（PBS）配制1:1,000的DRAQ5溶液用于细胞核染色。使用涡旋混合器将溶液充分混匀。</li>
	<li>在步骤1.2的20分钟孵育结束后，将组织从DTAF溶液转移至DRAQ5溶液中，在室温下避光孵育10分钟。</li>
</ol>
2. 肌腱加载与图像采集
注意：本实验方案需要一种可安装在共聚焦显微镜上方的拉伸装置。本研究中使用了 Peterson 和 Szczesny 所描述的微型拉伸装置13。
<ol>
	<li>将肌腱放入拉伸加载装置的夹具中。在将夹具安装到加载装置之前，使用数字卡尺测量跟骨附着点与对侧夹具之间的距离。该距离即为肌腱的标距长度。
	<ol>
		<li>或者，可先将夹具安装到加载装置上，再插入肌腱，并推至接触以确定电机的零位位移位置。在插入肌腱后记录电机的位移，可能获得更准确的夹具间标距长度。</li>
	</ol>
	</li>
	<li>将夹具安装到含有 PBS 的加载装置中，以维持组织的水合状态。尽可能将肌腱与显微镜图像的 x 轴或 y 轴对齐，以便算法输出的 x 方向应变和 y 方向应变与肌腱的轴向相对应。 
	注意：本研究中，肌腱沿x轴方向排列。若无法将肌腱完全与图像坐标轴对齐，则可使用标准应变转换方程，将算法输出的x方向应变和y方向应变转换为与肌腱纵向/垂直方向对齐的应变值。14.</li>
	<li>预加载肌腱，施加1 g张力，如有需要，可施加循环载荷以对样本进行预处理。本方案中未使用预处理，因为研究目的是验证测量的局部组织应变，而非测定组织的材料特性。如果需要测量依赖于加载历史的宏观尺度材料特性，则建议进行预处理。完成预处理及恢复后，重新施加1 g预加载。</li>
	<li>如需要，对四条等间距的线进行光漂白 80 µm 组织中心区域分离（参见 Peterson 和 Szczesny）13 有关更多详细信息）。 
	注意：光漂白线条用于验证ALDIC算法的测量结果，并非执行ALDIC所必需。线条的数量和间距可根据需要进行调整，且应选择线条的位置以避免样品中可能降低线条清晰度的人工假象。</li>
	<li>在靠近夹具的组织左右两端重复进行光漂白操作。</li>
	<li>使用共聚焦显微镜采集三维图像（x, y: 1.25 µm/pixel，z： 2.5 µm/pixe在1 g预加载条件下，DTAF和DRAQ5荧光的强度。</li>
	<li>以0.5%/s的应变速率进行应变 ramp 至2%应变。注意，应变速率和增量应变幅度可进行调整。</li>
	<li>让组织应力松弛10分钟。 
	注意：应力松弛的持续时间应选择为使样品在图像采集期间承受近似准静态载荷。为判断应力松弛持续时间是否合适，需确定应力松弛最后1分钟内力-时间曲线的斜率补充图1），并将该斜率乘以总成像持续时间。在本研究中，最大应变增量处施加的力变化从未超过5%。</li>
	<li>在组织变形后，再获取一幅容积图像。</li>
	<li>重复步骤 2.7–2.9，直至达到所需的目标应变值。本文中选择的最终应变为 12%。</li>
</ol>
3. 图像处理
<ol>
	<li>使用 ImageJ 或 Fiji 对 DRAQ5（细胞核）通道的每个体视图像创建最大强度 z 投影。这些投影将作为 ALDIC 的二维斑点图像。</li>
	<li>将最大强度 z 投影保存为 .tiff 文件，并按照以下命名规范进行命名。
	<ol>
		<li>使用数字作为图像名称的第一个字符。</li>
		<li>该数字应与图像在菌株分析过程中被处理的顺序相对应。例如，第一幅图像应以数字“1”开头，第二幅图像应以数字“2”开头。可以选择不同的起始数字，但必须按顺序递增。示例命名规范如下：“0_Experiment1_MaxZProjection”。</li>
	</ol>
	</li>
	<li>将所有重命名后的最大强度 z 投影文件保存到一个文件夹中。</li>
</ol>
4. 光漂白条带分析代码的安装与应用
注意：仅当需要使用光漂白线验证 ALDIC 算法的准确性时，才需执行以下步骤。该代码通过计算每组光漂白线中各光漂白线之间距离的平均归一化变化，来确定局部组织应变。在本研究中，进一步将所有光漂白线组（即位于中心及左右两端的线组）的平均局部值进行平均，以获得每个样本的单一平均局部组织应变值。该值随后用于评估 ALDIC 算法的准确性。
<ol>
	<li>从 GitHub 下载“PBL Code”文件夹（https://github.com/Szczesnytendon/TendonStrainCalc），并将其中所有内容移至 MATLAB 的工作目录。</li>
	<li>打开“Micro_Mech_Template.m”MATLAB 脚本。
	<ol>
		<li>点击 运行，并选择包含体图像的一个图像文件。体图像可以是以下任意一种文件类型：.lsm、.tiff、.nd2。</li>
		<li>软件将自动加载文件夹中的所有图像，并显示参考体图像的投影图像。当提示时，单击鼠标左键 创建多点线，以追踪样品的左端和右端；单击鼠标右键 结束一条线。输入完成后，若边缘位置正确，则点击 确定 接受结果。</li>
		<li>当提示时，在样品上绘制一条随机的对角线作为参考线。</li>
		<li>输入所创建的光漂白线条数量，并使用多点线追踪这些光漂白线。</li>
		<li>如果结果可接受，则确认接受；若结果有误，则进行调整并重新处理。</li>
	</ol>
	</li>
	<li>对所有图像重复步骤 4.2，并将所有已标注线条的图像移至同一个文件夹中。</li>
	<li>打开脚本“Micro_Mech_Strain.m”。
	<ol>
		<li>点击 运行 执行代码，并选择一个已保存的、包含已追踪光漂白线的图像。</li>
		<li>选择图像后，点击 确定 确认所选的配套图像是否正确。</li>
	</ol>
	</li>
</ol>
5. 创建数字化转换图像
注意：仅当需要使用数字化转换图像来验证 ALDIC 算法的准确性时，才需要执行以下步骤。这些图像通过人工变换参考图像，模拟出已知大小的均匀二维应变场。
<ol>
	<li>从 GitHub 下载代码“Digital_strain.m”（https://github.com/Szczesnytendon/TendonStrainCalc）。</li>
	<li>打开并运行该代码。</li>
	<li>当系统提示时，输入所需的最大施加应变值、施加应变增量和泊松比，然后点击 确定。 
	注意：本实验中，最大施加应变为 0.1（10%），施加应变增量为 0.02（2%），泊松比采用 1，该值与肌腱拉伸实验数据一致15,16。该代码使用 MATLAB 内嵌函数 imwarp 及输入参数（如应变增量、泊松比）生成数字化变换后的图像。</li>
	<li>当系统提示时，选择未变形的参考图像。</li>
	<li>对于每个应变增量，将显示参考图像与变换后图像的叠加图。变换后的图像将保存至目录中，文件名为“DigitallyTransformedX%Strain”，其中 X 为应变增量。</li>
</ol>
6. 应变计算与验证代码的安装与应用
<ol>
	<li>从 GitHub 下载“Strain Calculation and Validation Code”文件夹（https://github.com/Szczesnytendon/TendonStrainCalc），并将所有内容移至 MATLAB 工作目录。</li>
	<li>根据 Yang 和 Bhattacharya12 的方法安装 mex C/C++ 编译器。步骤总结如下：
	<ol>
		<li>在 MATLAB 命令窗口中输入“mex -setup”并按 Enter，检查是否已安装 mex C/C++ 编译器。</li>
		<li>如果出现错误提示，表明编译器不支持或不存在，则继续执行步骤 6.3 和步骤 6.4。</li>
		<li>如果没有错误，则继续执行步骤 6.5。</li>
	</ol>
	</li>
	<li>要下载 mex C/C++ 编译器，请访问“https:/tdm-gcc.tdragon.net/”，并选择 TDM-gcc 编译器。</li>
	<li>将下载的编译器安装到已知位置。</li>
	<li>返回 MATLAB 命令窗口，输入：“setenv("MW_MINGW64_LOC","[在此输入您的安装路径]")”。例如，可输入“setenv("MW_MINGW64_LOC","C:\TDM-GCC-64")”。如果该命令成功执行，则表示 mex 编译器已正确安装。</li>
	<li>进入“main_aldic.m”函数脚本，并将第 22 行修改为与步骤 6.5 中执行的命令一致。</li>
	<li>打开脚本“Strain_calc_and_validate.m”。</li>
	<li>点击 运行 以开始图像分析。</li>
	<li>当系统提示时，根据需要调整 ALDIC 参数的值。 
	注意：窗口大小应为子集大小的 0.25 至 1 倍。有关参数选择的更多信息，请参阅在线用户手册：(https://www.researchgate.net/publication/344796296_Augmented_Lagrangian_Digital 
	_Image_Correlation_AL-DIC_Code_Manual)。
	<ol>
		<li>本研究中使用了以下参数值： 
		子集大小（像素）：20 
		窗口大小（像素）：10 
		ALDIC 求解方法：有限差分法（1） 
		未使用并行计算（1） 
		初始猜测值计算方法：基于图像金字塔的多重网格搜索（0）</li>
	</ol>
	</li>
	<li>当系统提示时，选中“是”复选框，使算法自动保存所需变量集合的平均值、标准差和二维图谱（例如，x 方向应变、y 方向应变、剪切应变、异常区域等）。选择需要保存的变量，然后点击 确定。</li>
	<li>当系统提示时，根据需要调整参数。
	<ol>
		<li>本实验中使用了以下参数值： 
		用于计算应变的邻近点数（numP）：12 
		异常区域识别的相关系数阈值（corr_threshold）：0.5 
		异常区域分析的子区域大小（像素）（Subsize）：32</li>
	</ol>
	</li>
	<li>当系统提示时，选择包含重命名后的最大强度 z 投影图像的文件夹。请注意，软件会自动执行增量式 ALDIC 以确定变形图像的应变场。即，每张变形图像都作为下一张变形图像的新“参考”图像。与累积式 ALDIC 相比，这种方法提高了结果的准确性（补充图 2），在累积式 ALDIC 中，每张变形图像均与原始（0% 应变）参考图像进行比较。若要执行累积分析，请加载图像，但仅选择原始参考图像和感兴趣的变形图像。 
	注意：正应变计算公式为 λ - 1，其中 λ 为组织拉伸率。组织拉伸率按 <img alt="特征值方程，λ=√(NTCN)，公式，数学分析。" src="/files/ftp_upload/64921/64921eq01.jpg" style="margin: auto;" data-altupdatedat="1788304762000" data-alt="Equation 1"> 计算，其中 N = [1 0]T 或 [0 1]T 分别对应 x 方向和 y 方向，且 C = FT F，其中 F 是通过 ALDIC 算法输出的每个数据点周围“numP”个点计算得到的变形梯度。剪切应变按 <img alt="角度相减公式：θ - π/2，三角函数公式图解，教学用途。" src="/files/ftp_upload/64921/64921eq02.jpg" style="margin: auto;" data-altupdatedat="1788304762000" data-alt="Equation 2"> 计算，其中 <img alt="静态平衡，θ=arccos(C₁₂/√C₁₁C₂₂)，力学应力分析公式。" src="/files/ftp_upload/64921/64921eq03.jpg" style="margin: auto;" data-altupdatedat="1788304762000" data-alt="Equation 3">。</li>
	<li>当系统提示时，单击鼠标左键创建一个四点的多边形，以定义测量应变的兴趣区域。从左上角的点开始，并按顺时针方向依次指定后续点。 
	注意：保存在 MATLAB 工作区中的变量“Storage”包含平均 x 方向应变、x 方向应变标准差、平均 y 方向应变、y 方向应变标准差、平均剪切应变、剪切应变标准差以及异常区域百分比的所有数值。异常区域根据步骤 6.13 中选定的兴趣区域内进行的相关系数分析来定义。“NuclearTrackingResults”文件夹（可通过修改第 555 和 556 行重命名）存储了步骤 6.10 中指定的所有图表。该文件夹还包含一个名为“Results”的电子表格文件，其中存储了步骤 6.10 中指定的所有均值和标准差。</li>
</ol>

<ol><li>Devkota, A. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Distributing+a+fixed+amount+of+cyclic+loading+to+tendon+explants+over+longer+periods+induces+greater+cellular+and+mechanical+responses%5BTitle%5D)+AND+%22Journal+of+Orthopaedic+Research%22%5BJournal%5D)">Distributing a fixed amount of cyclic loading to tendon explants over longer periods induces greater cellular and mechanical responses.</a> Journal of Orthopaedic Research. 11 (4), 1609-1612 (2007).</li>
<li>Sun, H. B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Cycle-dependent+matrix+remodeling+gene+expression+response+in+fatigue-loaded+rat+patellar+tendons%5BTitle%5D)+AND+%22Journal+of+Orthopaedic+Research%22%5BJournal%5D)">Cycle-dependent matrix remodeling gene expression response in fatigue-loaded rat patellar tendons.</a> Journal of Orthopaedic Research. 28 (10), 1380-1386 (2010).</li>
<li>Shepherd, J. H., Screen, H. R. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Fatigue+loading+of+tendon%5BTitle%5D)+AND+%22International+Journal+of+Experimental+Pathology%22%5BJournal%5D)">Fatigue loading of tendon.</a> International Journal of Experimental Pathology. 94 (4), 260-270 (2013).</li>
<li>Paschall, L., Pedaprolu, K., Carrozzi, S., Dhawan, A., Szczesny, S. Mechanical stimulation as both the cause and the cure of tendon and ligament injuries. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aPaschall%2C+L+%22Regenerative+Rehabilitation%3A+From+Basic+Science+to+the+Clinic%22">Regenerative Rehabilitation: From Basic Science to the Clinic</a>. , Springer. Cham, Switzerland. 359-386 (2022).</li>
<li>Andarawis-Puri, N., Ricchetti, E. T., Soslowsky, L. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Rotator+cuff+tendon+strain+correlates+with+tear+propagation%5BTitle%5D)+AND+%22Journal+of+Biomechanics%22%5BJournal%5D)">Rotator cuff tendon strain correlates with tear propagation.</a> Journal of Biomechanics. 42 (2), 158-163 (2009).</li>
<li>Cheng, V. W. T., Screen, H. R. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+micro-structural+strain+response+of+tendon%5BTitle%5D)+AND+%22Journal+of+Materials+Science%22%5BJournal%5D)">The micro-structural strain response of tendon.</a> Journal of Materials Science. 42 (21), 8957-8965 (2007).</li>
<li>Luyckx, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Digital+image+correlation+as+a+tool+for+three-dimensional+strain+analysis+in+human+tendon+tissue%5BTitle%5D)+AND+%22Journal+of+Experimental+Orthopaedics%22%5BJournal%5D)">Digital image correlation as a tool for three-dimensional strain analysis in human tendon tissue.</a> Journal of Experimental Orthopaedics. 1 (1), 7(2014).</li>
<li>Duncan, N. A., Bruehlmann, S. B., Hunter, C. J., Shao, X., Kelly, E. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((In+situ+cell-matrix+mechanics+in+tendon+fascicles+and+seeded+collagen+gels%3A+Implications+for+the+multiscale+design+of+biomaterials%5BTitle%5D)+AND+%22Computer+Methods+in+Biomechanics+and+Biomedical+Engineering%22%5BJournal%5D)">In situ cell-matrix mechanics in tendon fascicles and seeded collagen gels: Implications for the multiscale design of biomaterials.</a> Computer Methods in Biomechanics and Biomedical Engineering. 17 (1), 39-47 (2014).</li>
<li>Arnoczky, S. P., Lavagnino, M., Whallon, J. H., Hoonjan, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((In+situ+cell+nucleus+deformation+in+tendons+under+tensile+load%3B+A+morphological+analysis+using+confocal+laser+microscopy%5BTitle%5D)+AND+%22Journal+of+Orthopaedic+Research%22%5BJournal%5D)">In situ cell nucleus deformation in tendons under tensile load; A morphological analysis using confocal laser microscopy.</a> Journal of Orthopaedic Research. 20 (1), 29-35 (2002).</li>
<li>Screen, H. R. C., Bader, D. L., Lee, D. A., Shelton, J. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Local+strain+measurement+within+tendon%5BTitle%5D)+AND+%22Strain%22%5BJournal%5D)">Local strain measurement within tendon.</a> Strain. 40 (4), 157-163 (2004).</li>
<li>Fung, A. K., Paredes, J. J., Andarawis-Puri, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Novel+image+analysis+methods+for+quantification+of+in+situ+3-D+tendon+cell+and+matrix+strain%5BTitle%5D)+AND+%22Journal+of+Biomechanics%22%5BJournal%5D)">Novel image analysis methods for quantification of in situ 3-D tendon cell and matrix strain.</a> Journal of Biomechanics. 67, 184-189 (2018).</li>
<li>Yang, J., Bhattacharya, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Augmented+Lagrangian+digital+image+correlation%5BTitle%5D)+AND+%22Experimental+Mechanics%22%5BJournal%5D)">Augmented Lagrangian digital image correlation.</a> Experimental Mechanics. 59 (2), 187-205 (2019).</li>
<li>Peterson, B. E., Szczesny, S. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Dependence+of+tendon+multiscale+mechanics+on+sample+gauge+length+is+consistent+with+discontinuous+collagen+fibrils%5BTitle%5D)+AND+%22Acta+Biomaterialia%22%5BJournal%5D)">Dependence of tendon multiscale mechanics on sample gauge length is consistent with discontinuous collagen fibrils.</a> Acta Biomaterialia. 117, 302-309 (2020).</li>
<li>Humphrey, J. D., O'Rourke, S. L. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aHumphrey%2C+JD+%22An+Introduction+to+Biomechanics%22">An Introduction to Biomechanics</a>. , Springer. New York, NY. (2015).</li>
<li>Reese, S. P., Weiss, J. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Tendon+fascicles+exhibit+a+linear+correlation+between+Poisson%27s+ratio+and+force+during+uniaxial+stress+relaxation%5BTitle%5D)+AND+%22Journal+of+Biomechanical+Engineering%22%5BJournal%5D)">Tendon fascicles exhibit a linear correlation between Poisson's ratio and force during uniaxial stress relaxation.</a> Journal of Biomechanical Engineering. 135 (3), 34501(2013).</li>
<li>Ahmadzadeh, H., Freedman, B. R., Connizzo, B. K., Soslowsky, L. J., Shenoy, V. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Micromechanical+poroelastic+finite+element+and+shear-lag+models+of+tendon+predict+large+strain+dependent+Poisson%27s+ratios+and+fluid+expulsion+under+tensile+loading%5BTitle%5D)+AND+%22Acta+Biomaterialia%22%5BJournal%5D)">Micromechanical poroelastic finite element and shear-lag models of tendon predict large strain dependent Poisson's ratios and fluid expulsion under tensile loading.</a> Acta Biomaterialia. 22, 83-91 (2015).</li>
<li>Szczesny, S. E., Elliott, D. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Interfibrillar+shear+stress+is+the+loading+mechanism+of+collagen+fibrils+in+tendon%5BTitle%5D)+AND+%22Acta+Biomaterialia%22%5BJournal%5D)">Interfibrillar shear stress is the loading mechanism of collagen fibrils in tendon.</a> Acta Biomaterialia. 10 (6), 2582-2590 (2014).</li>
<li>Han, W. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Macro-+to+microscale+strain+transfer+in+fibrous+tissues+is+heterogeneous+and+tissue-specific%5BTitle%5D)+AND+%22Biophysical+Journal%22%5BJournal%5D)">Macro- to microscale strain transfer in fibrous tissues is heterogeneous and tissue-specific.</a> Biophysical Journal. 105 (3), 807-817 (2013).</li>
<li>Pedaprolu, K., Szczesny, S. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+novel%2C+open-source%2C+low-cost+bioreactor+for+load-controlled+cyclic+loading+of+tendon+explants%5BTitle%5D)+AND+%22Journal+of+Biomechanical+Engineering%22%5BJournal%5D)">A novel, open-source, low-cost bioreactor for load-controlled cyclic loading of tendon explants.</a> Journal of Biomechanical Engineering. 144 (8), 084505(2022).</li>
<li>Gatt, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Negative+Poisson%27s+ratios+in+tendons%3A+An+unexpected+mechanical+response%5BTitle%5D)+AND+%22Acta+Biomaterialia%22%5BJournal%5D)">Negative Poisson's ratios in tendons: An unexpected mechanical response.</a> Acta Biomaterialia. 24, 201-208 (2015).</li>
</ol>

所有作者均无利益冲突需要披露。

本文的目的是提供一种开源且经过验证的方法，用于测量肌腱在拉伸载荷下的二维应变场。该软件的基础是基于一种公开可用的ALDIC算法12。该算法被嵌入到一个更大的MATLAB程序中，并增加了增量（而非累积）应变分析的功能。随后，将这一改进后的算法应用于肌腱的拉伸测试，并通过两种不同技术（即数字变换图像和利用光漂白线进行应变测量）评估其准确性。此外，还增加了一项功能，可在无需已知真实应变值的情况下，评估ALDIC测量在任意样本上的准确性。
对数字化转换图像的分析表明，该算法能够以极小的误差精确测量高达10%的应变，而这一精度无法通过实际小鼠跟腱的拉伸测试进行评估，因为跟腱样本中的应变幅值较低。然而，将ALDIC算法计算的小鼠跟腱应变与利用光漂白标记线测量的应变进行比较后发现，ALDIC技术的误差处于光漂白标记线自身测量变异范围内。作为最终验证，通过从变形图像重建参考图像，并将重建结果与实际参考图像进行比较，评估了ALDIC算法计算的完整二维位移场的准确性。在数字化转换的图像中，随着施加应变的增大，异常区域数量和应变误差增加，尤其是在累积式AL...

肌腱是具有机械敏感性的组织，会因机械负荷而发生适应性改变或退变1,2,3,4。由于机械刺激在肌腱细胞生物学中起着重要作用，人们高度关注肌腱细胞在天然组织环境中受载时所经历的应变情况。目前已开发出多种实验和分析技术用于测量肌腱中的局部组织应变，包括利用散斑图案或光漂白线（PBLs）进行表面应变的2D/3D数字图像相关（DIC）分析5,6,7,8、通过测量组织内单个细胞核中心到中心距离的变化来评估应变9,10，以及一种考虑了面外运动和三维变形的新型全场3D DIC方法11。然而，这些技术的准确性与灵敏度仅在少数研究中被报道，且目前尚无任何技术实现公开共享，这限制了这些方法的广泛应用与推广。
本研究的目的是开发一种经过验证的分析工具，用于测量肌腱离体样本中的局部组织应变，该工具应易于获取且操作简便。所采用的方法基于一种公开可用的增广拉格朗日数字图像相关（ALDIC）算法，该算法由 Yang 和 Bhattacharya 用 MATLAB 编写12。该算法已针对肌腱样本的分析进行了改进，并通过应用于数字化变换图像进行验证，同时将实际肌腱样本中测得的应变与光漂白线方法获得的结果进行比较。此外，还在算法中增加了额外功能，以便在缺乏已知应变值或无二次测量技术的情况下，仍能验证计算所得位移场的准确性。因此，该算法是一种高度实用且灵活的工具，可用于精确测量肌腱中的局部二维组织应变。

<table><tbody><tr><td>5-DTAF（5-(4,6-二氯三嗪基)氨基荧光素），单一异构体</td><td>ThermoFisher</td><td>D16</td><td></td></tr><tr><td>卡尺</td><td>Mitutoyo</td><td>500-196-30</td><td></td></tr><tr><td>共聚焦显微镜</td><td>Nikon</td><td>A1R HD</td><td></td></tr><tr><td>Corning LSE 涡旋混合器</td><td>Corning</td><td>6775</td><td></td></tr><tr><td>DRAQ5 荧光探针溶液（5 mM）</td><td>ThermoFisher</td><td>62554</td><td></td></tr><tr><td>MATLAB</td><td>MathWorks</td><td>R2022b</td><td></td></tr><tr><td>拉伸加载装置</td><td>N/A</td><td>N/A</td><td>如 Peterson 等，2020 年所述的拉伸加载装置（参考文献 13）&amp;nbsp;</td></tr><tr><td>Tube Revolver 旋转混合仪</td><td>ThermoFisher</td><td>88881001</td><td></td></tr></tbody></table>

测量肌腱中的局部组织应变 通过 开源数字图像相关法

Measuring Local Tissue Strains in Tendons via Open-Source Digital Image Correlation

肌腱细胞所承受的应变引起了广泛的科学关注 原位 以及这些应变如何影响组织重塑。基于这一研究兴趣，已开发出多种分析技术，用于测量加载过程中肌腱离体样本内的局部组织应变。然而，在多个案例中，这些技术的准确性和灵敏度尚未被报道，且目前尚无任何算法公开可用。这使得在肌腱离体样本中广泛开展局部组织应变测量变得困难。因此，本研究的目标是开发一种经过验证、易于使用且可公开获取的分析工具，用于测量肌腱离体样本中的局部组织应变。具体而言，本研究对一种公开可用的增强拉格朗日数字图像相关（ALDIC）算法进行了改进，通过追踪小鼠跟腱在单轴拉伸过程中细胞核的位移，实现二维应变的测量。此外，通过分析数字变换后的图像，并将计算所得应变值与独立技术（即光漂白条纹法）测得的结果进行比较，验证了所计算应变的准确性。最后，本研究在算法中引入了一种技术，利用计算得到的位移场重建参考图像，可在缺乏已知应变值或第二种测量方法的情况下评估算法的准确性。该算法能够测量高达0.1的应变，精度达0.00015。通过比较重建参考图像与实际参考图像的技术，成功识别出含有错误数据的样本，并表明在数据质量良好的样本中，约85%的位移场是准确的。最后，小鼠跟腱中测得的应变结果与以往文献报道一致。因此，该算法是一种高度实用且可适应性强的工具，可用于准确测量肌腱中的局部组织应变。

在分析实际组织样本的应变场之前，首先使用小鼠跟腱内细胞核的数字化应变/变换图像对ALDIC方案进行验证。具体而言，通过对图像进行变换，以数字方式在x方向上产生2%、4%、6%、8%和10%的均匀应变，并模拟泊松比为115,16通过将ALDIC算法计算得到的平均应变值与已知的数字应变进行比较，评估该算法的准确性。此外，通过评估应变值的标准差以确定应变场的非均质性。ALDIC算法（采用增量分析方法）计算得到的应变与施加于数字变换图像的实际应变之间的差异如图所示。 图1. ALDIC 软件计算得到的 x 方向平均应变始终低于实际施加的应变（图1A），且误差的幅度随施加的应变增大而增加。然而，对于所有应变增量，误差幅度始终小于0.00015。在y方向上也存在轻微的应变低估现象（图1C）。在x应变和y应变的整个感兴趣区域内，计算应变的标准差随施加应变的增...

本文介绍了一种用于测量肌腱离体组织内局部二维组织应变的开源数字图像相关算法。该技术的准确性已通过多种方法得到验证，并可供公众使用。

肌腱细胞在体内所承受的应变及其对组织的影响，一直是科学界广泛关注的研究课题……

本方案能够实现对肌腱在受载过程中局部组织应变的测量。该方法有助于理解肌腱的力学行为及其对加载环境的适应机制。该技术的主要优势在于所使用的软件可在线公开获取，并且即使在真实组织应变未知的情况下，也能提供测量准确性的反馈信息。测量肌腱中局部组织的应变对于理解细胞在肌腱退变和修复过程中如何重塑组织具有重要意义。本实验操作将由本实验室的本科生 Stanton Godshall 和博士研究生 Krishna Pedaprolu 演示。将收获的跟腱在 DTAF 中染色后，从 DTAF 溶液转移至 DRAQ5 溶液中，并在室温下避光孵育 10 分钟。将肌腱放入拉伸加载装置的夹具中。在将夹具安装到加载装置之前，使用数字卡尺测量跟骨附着点与对侧夹具之间的距离。将夹具安装到含有 PBS 的加载装置中，以保持组织的水合状态。尽可能将肌腱与显微镜图像的 x 轴或 y 轴对齐，以使算法输出的 X 方向应变和 Y 方向应变与肌腱轴线相对应。用 1 克张力对肌腱进行预加载。如需要，可在组织中心区域光漂白四条间距为 80 微米的线，并在靠近夹具的左右两端重复此过程。各组织区域中线条间距离的变化是局部组织应变的间接测量指标，可用于验证数据。接下来，使用共聚焦显微镜获取DTAF和DRAQ5荧光在1克预载荷下的三维图像。以每秒0.5%的速率进行至2%应变的应变递增。请注意，应变率和增量应变幅度可以调整。在让组织应力松弛10分钟后，对变形后的组织再采集一次体积图像。根据需要，重复此过程进行多次应变增量。为生成数字化转换图像，请从 GitHub 下载代码 digital_strain.m。下载后，打开并运行该代码。当系统提示时，在按下“确定”按钮前，输入所需的最大施加应变、施加应变增量和泊松比值。随后当再次出现提示时，选择未变形的参考图像。对于每个应变增量，将显示参考图像与变换后图像的叠加图。经数字变换的图像将保存至名为“digitally transformed X percent strain”的目录中，其中X为应变增量。打开显示名称的脚本，然后点击“运行”以开始图像分析。当系统提示时，根据需要调整增强拉格朗日数字图像相关（ALDIC）参数的数值。在提示后，选中“是”复选框，以自动保存所需变量集合的平均值、标准差和二维图谱。当系统提示时，选择所需的变量，如 X 应变、Y 应变、剪切应变、异常区域等，然后点击“确定”。在下一个提示出现后，选择包含已重命名的最大强度 Z 投影图像的文件夹。软件将自动执行增量 ALDIC 分析，以确定变形图像的应变场。当系统提示时，单击鼠标左键创建一个四点闭合多边形，以定义用于测量应变的感兴趣区域。文件夹“nuclear tracking results”可通过对第555和556行代码进行调整来重命名，该文件夹存储了之前指定的所有图表。此文件夹还包含一个名为“results”的电子表格，用于存储之前指定的所有均值和标准差。在使用数字施加应变的图像进行验证时，ALDIC算法始终低估了平均X方向应变，且误差幅度随施加应变的增加而增大。Y 应变也大多被低估了。然而，在所有情况下，应变误差的幅度都非常小。尽管计算得到的 X 应变和 Y 应变的标准差数值较低，但随着施加应变的增加而增大。不良区域分析显示，在使用累积法分析的数字变换图像中，当施加应变超过6%后，出现无效应变计算的不良区域数量持续增加，而增量始终保持在1%。在四个样本中的一个样本（被视为异常值）中，在最大应变增量时，近一半的图像区域被识别为不良区域。由ALDIC计算得到的X方向应变大于通过PBLs测定的结果，两者差异在0.005以内，与所有样本PBL数据的平均标准差相当。通过测定腱组织在拉伸载荷下的局部应变大小及其空间分布发现，在所有样本中，X方向应变始终低于施加的应变。在所有应变增量中，Y 方向的平均应变近似为零，但标准差较高。平均剪切应变在整个应变增量过程中持续稳定增加。最重要的一点是，在施加更大应变之前，务必先保存图像。如果未保存，图像将丢失且无法恢复。尽管该技术已专门验证用于测量肌腱内的组织应变，但它也可为许多其他动物和人体组织的力学生物学及力学特性研究提供见解。

measuring-local-tissue-strains-tendons-via-open-source-digital-image

Research

科研

JoVE Journal

Engineering

工程

通过开源数字图像相关 测量 肌腱中的局部组织应变

我们描述了使用数字图像相关技术来表征经受单轴拉伸试验的血管组织样本表面局部应变场的方法。这些测量有助于精确量化样本的力学响应，并建立本构应力-应变关系。...

利用数字图像相关技术表征血管组织样本上的局部应变

在本实验方案中，采用数字图像相关技术估算股骨表面在骨折测试过程中的应变。该方法的创新之处在于股骨表面应用了高对比度的随机散斑图案，并结合了精确设定的照明条件、高速视频采集以及数字图像相关分析来进行应变计算。...

一种利用数字图像相关法在骨折测试期间估算尸体股骨皮质应变的方法

采用一种结合裂纹扩展速率测量与应变场分析的新方法，研究微结构尺度下疲劳短裂纹的扩展行为，以揭示晶粒尺度以下的累积变形场。...

使用数字图像相关法对微结构小疲劳裂纹扩展进行全场应变测量

我们开发了一种力学成像分析流程，用于研究动脉粥样硬化斑块在结构和力学特性上的异质性。该分析流程能够关联纤维斑块组织中胶原纤维取向的局部主导角度与分散程度、破裂行为以及应变特征图谱。...

研究局部胶原结构与动脉粥样硬化斑块纤维组织力学特性相关性的方法

我们介绍一种定制的实验平台和组织培养方案，该方案可重现由跟腱止点撞击引起的纤维软骨变性，同时在小鼠后肢离体组织中维持细胞长期存活，为探索腱组织撞击的力学生物学机制提供了一个合适的模型。...

用于研究跟腱撞击力学生物学的小鼠后肢离体模型

本方案采用基于直接线性变换校准的立体成像相机系统，用于采集三维数据 原位 拉伸状态下周围神经的位移。通过捕捉这些位移，可确定在不同拉伸程度下所诱发的应变，从而明确与拉伸相关的神经损伤阈值，推动基于拉伸机制的神经修复科学的发展。...

直接线性变换在测量中的应用 原位 拉伸过程中的周围神经张力

本方案描述了一种标准化、便携式的基于超声的检测方法，用于在精英运动员中，于多个踝关节角度下量化跟腱的功能性刚度谱，从而实现对不同负荷条件下肌腱力学行为的可靠且可重复的评估。...

一种基于振动的弹性成像方案在评估精英运动员跟腱刚度随关节角度变化时的可靠性

本文描述了一种能够获取光学不透明流体中或通过光学不透明几何结构的二维速度场的回声粒子图像测速（EPIV）系统，并报道了在管道流动中的验证性测量结果。...

回声粒子图像测速法

本文介绍了一种利用高分辨率、高速粒子图像测速技术（PIV）研究边界附近瞬态流动的方法。PIV是一种非侵入式测量技术，通过优化图像与记录特性、激光片光源特性以及分析算法等多个参数约束，可应用于任何光学可访问的流动研究。...

近表面高速粒子图像测速技术

边缘场静电微机电系统执行器的坚固器件设计在使用传统阶跃偏置进行开关操作时，会导致固有的低挤压膜阻尼条件和较长的稳定时间。采用直流动态波形实现实时开关时间优化，可缩短边缘场微机电系统执行器在从上至下和从下至上状态切换过程中的稳定时间。...

严重欠阻尼边缘场静电MEMS执行器中用于改善开关时间的实时直流动态偏置方法

了解开源数字图像相关算法如何测量肌腱局部组织应变。该技术使用ALDIC方法分析Achilles肌腱样本，验证测量准确性并提供反馈，助力生物力学与组织工程研究。

测量肌腱中的局部组织应变通过开源数字图像相关法

探索更多视频

测量肌腱中的局部组织应变 通过 开源数字图像相关法

探索更多视频

测量肌腱中的局部组织应变通过开源数字图像相关法