答案 0 :(得分:1)
regionprops
非常适合这类任务。您可以将分段二进制图像和原始图像传递给它,以检索每个区域中的像素列表(假设每个区域未连接,如示例图像所示)。这些将以每个区域的n x 1
向量的形式出现,作为结构数组返回。
stats = regionprops(BW, I, 'PixelValues');
(您可能希望检索regionprops
返回的其他值,例如BoundingBox
或Centroid
,以帮助更轻松地识别哪个像素集属于哪个区域。请参阅文档看看有哪些选择。)。
然后,您可以定义一些统计函数来显示每个区域内的变化,例如,计算每个区域的方差和标准偏差:
for n = 1:length(stats)
stats(n).var = var(stats(n).PixelValues);
stats(n).std = std(stats(n).PixelValues);
end
如果您对"强度变化有其他具体定义"记住,那么你需要开发一些计算它的函数,然后调用它而不是内置的var
或std
。