5

我需要获得一组具有亚像素精度的二进制图像中组件的质心（质心）。

Mathematica 8 有一个很好的补充： 

i =  Import@"http://i.stack.imgur.com/2pxrN.png";  
m1 = ComponentMeasurements[MorphologicalComponents[i], "Centroid"] /. 
                                                                Rule[_, x_] -> x
(*
-> {{403.229, 453.551}, {660.404, 371.383}, {114.389, 434.646}, {295.5, 206.}}
*)

但是当这些结果显示与其他地方进行的其他计算不一致时，我遇到了一些麻烦。

所以我推出了自己的，也许不太好： 

i = Import@"http://i.stack.imgur.com/2pxrN.png";
f[i_] := N@{#[[2]], ImageDimensions[i][[2]] - #[[1]]} & /@
         ( Mean /@
             Function[x, Map[
                Position[x, #, 2] &,
                Complement[Union@Flatten[x], {0}]]]
             [MorphologicalComponents[i]]);
f[i]
Show[i, Graphics[{Red, Disk[#, 10] & /@ f[i]}]]
(*
-> {{403.729, 453.051}, {660.904, 370.883}, {114.889, 434.146}, {296., 205.5}}
*)

在此处输入图像描述

您可以看到这些结果之间存在 0.5 的偏移： 

Thread[Subtract[m1, f[i]]]
(*
-> {{-0.5, -0.5, -0.5, -0.5}, {0.5, 0.5, 0.5, 0.5}}
*)

起初我认为问题与图像尺寸是偶数还是奇数有关，但事实并非如此。

我更喜欢使用ComponentMeasurements[ ..,"Centroid"]和纠正结果，但我担心未来的 Mma 版本可能会修改这种行为并破坏结果。

我还可以使用已知图像运行先前的“校准”并计算偏移量，因此它会自动更正，但我想先了解发生了什么。

这是一个错误吗？
关于为什么会发生这种情况的任何想法？

4

1 回答 1

6

我觉得文档页面ComponentMeasurements包含解决方案：

位置、面积和长度测量在标准图像坐标系中进行，其中位置 {0,0} 对应于左下角，x 从 0 到宽度，y 从 0 到高度。

您正在计算整个像素并ComponentMeasurements测量像素位置。在此系统中，左下像素的中心位于 {1/2,1/2}。

于 2011-06-24T19:23:23.090 回答