DM 脚本计算具有给定 ROI 的曲线下面积答案

【问题标题】：DM script to calculate area under a curve with given ROIDM 脚本计算具有给定 ROI 的曲线下面积
【发布时间】：2021-10-27 17:38:27
【问题描述】：

我想计算一个给定ROI的曲线下面积，如附图，我想得到13到30的绿色曲线下面积，我们可以通过DM脚本来实现这个功能吗？

以前没有做过类似的事情，任何建议谢谢！

【问题讨论】：

标签： dm-script

【解决方案1】：

如果图像的标签为“A”，那么您可以简单地执行以下操作： sum( A[] )

或者，更详细一点：

image img := GetFrontImage()
number integral = sum( img[] )
result("\n The intensity of the selected region in image [")
result( img.imageGetLabel() + "]:" + img.ImageGetName() + " is: ")
result( integral)

请注意，这将为您提供原始数字中的积分（或只是通道值的总和）。如果您已校准强度，则必须将校准考虑在内。如果 - 在大多数情况下 - 你只有一个强度等级但没有偏移量，它是直截了当的：

image img := GetFrontImage()
number integral = sum( img[] )
number scale = img.ImageGetIntensityScale()
number integral_cal = integral * scale
result("\n The intensity of the selected region in image [")
result( img.imageGetLabel() + "]:" + img.ImageGetName() + " is: ")
result( integral_cal )

如果你也有一个起源，它变成：

image img := GetFrontImage()
number integral = sum( img[] )
number scale = img.ImageGetIntensityScale()
number origin = img.ImageGetIntensityOrigin()

number t,l,b,r
img.GetSelection(t,l,b,r)
number nChannels = r - l
number integral_cal = (integral - origin * nChannels) * scale

result("\n The intensity of the selected region in image [")
result( img.imageGetLabel() + "]:" + img.ImageGetName() + " is: ")
result( integral_cal )

【讨论】：