Pixel connectivity

Relevant topics on
Graph connectivity
  • Connectivity
  • Algebraic connectivity
  • Cycle rank
  • Rank (graph theory)
  • SPQR tree
  • St-connectivity
  • K-connectivity certificate
  • Pixel connectivity
  • Vertex separator
  • Strongly connected component
  • Biconnected graph
  • Bridge
  • v
  • t
  • e

In image processing, pixel connectivity is the way in which pixels in 2-dimensional (or hypervoxels in n-dimensional) images relate to their neighbors.

Formulation

The 9 possible connectivities in a 5x5x5 neighborhood

In order to specify a set of connectivities, the dimension N {\displaystyle N} and the width of the neighborhood n {\displaystyle n} , must be specified. The dimension of a neighborhood is valid for any dimension n 1 {\displaystyle n\geq 1} . A common width is 3, which means along each dimension, the central cell will be adjacent to 1 cell on either side for all dimensions.

Let M N n {\displaystyle M_{N}^{n}} represent a N-dimensional hypercubic neighborhood with size on each dimension of n = 2 k + 1 , k Z {\displaystyle n=2k+1,k\in \mathbb {Z} }

Let q {\displaystyle {\vec {q}}} represent a discrete vector in the first orthant from the center structuring element to a point on the boundary of M N n {\displaystyle M_{N}^{n}} . This implies that each element q i { 0 , 1 , . . . , k } , i { 1 , 2 , . . . , N } {\displaystyle q_{i}\in \{0,1,...,k\},\forall i\in \{1,2,...,N\}} and that at least one component q i = k {\displaystyle q_{i}=k}

Let S N d {\displaystyle S_{N}^{d}} represent a N-dimensional hypersphere with radius of d = q {\displaystyle d=\left\Vert {\vec {q}}\right\Vert } .

Define the amount of elements on the hypersphere S N d {\displaystyle S_{N}^{d}} within the neighborhood M N n {\displaystyle M_{N}^{n}} as E {\displaystyle E} . For a given q {\displaystyle {\vec {q}}} , E {\displaystyle E} will be equal to the amount of permutations of q {\displaystyle {\vec {q}}} multiplied by the number of orthants.

Let n j {\displaystyle n_{j}} represent the amount of elements in vector q {\displaystyle {\vec {q}}} which take the value j {\displaystyle j} . n j = i = 1 N ( q i = j ) {\displaystyle n_{j}=\sum _{i=1}^{N}(q_{i}=j)}

The total number of permutation of q {\displaystyle {\vec {q}}} can be represented by a multinomial as N ! j = 0 k n j ! {\displaystyle {\frac {N!}{\prod _{j=0}^{k}n_{j}!}}}

If any q i = 0 {\displaystyle q_{i}=0} , then the vector q {\displaystyle {\vec {q}}} is shared in common between orthants. Because of this, the multiplying factor on the permutation must be adjusted from 2 N {\displaystyle 2^{N}} to be 2 N n 0 {\displaystyle 2^{N-n_{0}}}

Multiplying the number of amount of permutations by the adjusted amount of orthants yields,

E = N ! j = 0 k n j ! 2 N n 0 {\displaystyle E={\frac {N!}{\prod _{j=0}^{k}n_{j}!}}2^{N-n_{0}}}

Let V {\displaystyle V} represent the number of elements inside of the hypersphere S N d {\displaystyle S_{N}^{d}} within the neighborhood M N n {\displaystyle M_{N}^{n}} . V {\displaystyle V} will be equal to the number of elements on the hypersphere plus all of the elements on the inner shells. The shells must be ordered by increasing order of q = r {\displaystyle \left\Vert {\vec {q}}\right\Vert =r} . Assume the ordered vectors q {\displaystyle {\vec {q}}} are assigned a coefficient p {\displaystyle p} representing its place in order. Then an ordered vector q p , p { 1 , 2 , . . . , x = 1 k ( x + 1 ) } {\displaystyle {\vec {q_{p}}},p\in \{1,2,...,\sum _{x=1}^{k}(x+1)\}} if all r {\displaystyle r} are unique. Therefore V {\displaystyle V} can be defined iteratively as

V q p = V q p 1 + E q p , V q 0 = 0 {\displaystyle V_{\vec {q_{p}}}=V_{\vec {q_{p-1}}}+E_{\vec {q_{p}}},V_{\vec {q_{0}}}=0} ,

or

V q p = x = 1 p E q x {\displaystyle V_{\vec {q_{p}}}=\sum _{x=1}^{p}E_{\vec {q_{x}}}}

If some q x = q y {\displaystyle \left\Vert {\vec {q_{x}}}\right\Vert =\left\Vert {\vec {q_{y}}}\right\Vert } , then both vectors should be considered as the same p {\displaystyle p} such that

V q p = V q p 1 + E q p , 1 + E q p , 2 , V q 0 = 0 {\displaystyle V_{\vec {q_{p}}}=V_{\vec {q_{p-1}}}+E_{\vec {q_{p,1}}}+E_{\vec {q_{p,2}}},V_{\vec {q_{0}}}=0}

Note that each neighborhood will need to have the values from the next smallest neighborhood added. Ex. V q = ( 0 , 2 ) = V q = ( 1 , 1 ) + E q = ( 0 , 2 ) {\displaystyle V_{{\vec {q}}=(0,2)}=V_{{\vec {q}}=(1,1)}+E_{{\vec {q}}=(0,2)}}

V {\displaystyle V} includes the center hypervoxel, which is not included in the connectivity. Subtracting 1 yields the neighborhood connectivity, G {\displaystyle G}

G = V 1 {\displaystyle G=V-1} [1]

Table of Selected Connectivities

Neighborhood Size:

M N n {\displaystyle M_{N}^{n}}

Connectivity Type Typical Vector:

q {\displaystyle {\vec {q}}}

Sphere Radius

d {\displaystyle d}

Elements on Sphere:

E {\displaystyle E}

Elements in Sphere:

V {\displaystyle V}

Neighborhood Connectivity:

G {\displaystyle G}

1 edge (0) 0 1*1=1 1 0
3 point (1) 1 1*2=2 3 2
5 point-point (2) 4 1*2=2 5 4
...
1x1 face (0,0) 0 1*1=1 1 0
3x3 edge (0,1) 1 2*2=4 5 4
point (1,1) 2 1*4=4 9 8
5x5 edge-edge (0,2) 4 2*2=4 13 12
point-edge (1,2) 5 2*4=8 21 20
point-point (2,2) 8 1*4=4 25 24
...
1x1x1 volume (0,0,0) 0 1*1=1 1 0
3x3x3 face (0,0,1) 1 3*2=6 7 6
edge (0,1,1) 2 3*4=12 19 18
point (1,1,1) 3 1*8=8 27 26
5x5x5 face-face (0,0,2) 4 3*2=6 33 32
edge-face (0,1,2) 5 6*4=24 57 56
point-face (1,1,2) 6 3*8=24 81 80
edge-edge (0,2,2) 8 3*4=12 93 92
point-edge (1,2,2) 9 3*8=24 117 116
point-point (2,2,2) 12 1*8=8 125 124
...
1x1x1x1 hypervolume (0,0,0,0) 0 1*1=1 1 0
3x3x3x3 volume (0,0,0,1) 1 4*2=8 9 8
face (0,0,1,1) 2 6*4=24 33 32
edge (0,1,1,1) 3 4*8=32 65 64
point (1,1,1,1) 4 1*16=16 81 80
...

Example

Consider solving for G | q = ( 0 , 1 , 1 ) {\displaystyle G|{\vec {q}}=(0,1,1)}

In this scenario, N = 3 {\displaystyle N=3} since the vector is 3-dimensional. n 0 = 1 {\displaystyle n_{0}=1} since there is one q i = 0 {\displaystyle q_{i}=0} . Likewise, n 1 = 2 {\displaystyle n_{1}=2} . k = 1 , n = 3 {\displaystyle k=1,n=3} since max q i = 1 {\displaystyle \max q_{i}=1} . d = 0 2 + 1 2 + 1 2 = 2 {\displaystyle d={\sqrt {0^{2}+1^{2}+1^{2}}}={\sqrt {2}}} . The neighborhood is M 3 3 {\displaystyle M_{3}^{3}} and the hypersphere is S 3 2 {\displaystyle S_{3}^{\sqrt {2}}}

E = 3 ! 1 ! 2 ! 0 ! 2 3 1 = 6 2 4 = 12 {\displaystyle E={\frac {3!}{1!*2!*0!}}2^{3-1}={\frac {6}{2}}4=12}

The basic q {\displaystyle {\vec {q}}} in the neighborhood N 3 3 {\displaystyle N_{3}^{3}} , q 1 = ( 0 , 0 , 0 ) {\displaystyle {\vec {q_{1}}}=(0,0,0)} . The Manhattan Distance between our vector and the basic vector is q q 0 1 = 2 {\displaystyle \left\Vert {\vec {q}}-{\vec {q_{0}}}\right\Vert _{1}=2} , so q = q 3 {\displaystyle {\vec {q}}={\vec {q_{3}}}} . Therefore,

G q 3 = V q 3 1 = E q 1 + E q 2 + E q 3 1 = E q = ( 0 , 0 , 0 ) + E q = ( 0 , 0 , 1 ) + E q = ( 0 , 1 , 1 ) {\displaystyle G_{\vec {q_{3}}}=V_{\vec {q_{3}}}-1=E_{\vec {q_{1}}}+E_{\vec {q_{2}}}+E_{\vec {q_{3}}}-1=E_{{\vec {q}}=(0,0,0)}+E_{{\vec {q}}=(0,0,1)}+E_{{\vec {q}}=(0,1,1)}}
E q = ( 0 , 0 , 0 ) = 3 ! 3 ! 0 ! 0 ! 2 3 3 = 6 6 1 = 1 {\displaystyle E_{{\vec {q}}=(0,0,0)}={\frac {3!}{3!*0!*0!}}2^{3-3}={\frac {6}{6}}1=1}
E q = ( 0 , 0 , 1 ) = 3 ! 2 ! 1 ! 2 3 2 = 6 2 2 = 6 {\displaystyle E_{{\vec {q}}=(0,0,1)}={\frac {3!}{2!*1!}}2^{3-2}={\frac {6}{2}}2=6}
G = 1 + 6 + 12 1 = 18 {\displaystyle G=1+6+12-1=18}

Which matches the supplied table

Higher values of k & N

The assumption that all q p = r {\displaystyle \left\Vert {\vec {q_{p}}}\right\Vert =r} are unique does not hold for higher values of k & N. Consider N = 2 , k = 5 {\displaystyle N=2,k=5} , and the vectors q A = ( 0 , 5 ) , q B = ( 3 , 4 ) {\displaystyle {\vec {q_{A}}}=(0,5),{\vec {q_{B}}}=(3,4)} . Although q A {\displaystyle {\vec {q_{A}}}} is located in M 2 5 {\displaystyle M_{2}^{5}} , the value for r = 25 {\displaystyle r=25} , whereas q B {\displaystyle {\vec {q_{B}}}} is in the smaller space M 2 4 {\displaystyle M_{2}^{4}} but has an equivalent value r = 25 {\displaystyle r=25} . q C = ( 4 , 4 ) M 2 4 {\displaystyle {\vec {q_{C}}}=(4,4)\in M_{2}^{4}} but has a higher value of r = 32 {\displaystyle r=32} than the minimum vector in M 2 5 {\displaystyle M_{2}^{5}} .

For this assumption to hold, { N = 2 , k 4 N = 3 , k 2 N = 4 , k 1 {\displaystyle {\begin{cases}N=2,k\leq 4\\N=3,k\leq 2\\N=4,k\leq 1\end{cases}}}

At higher values of k {\displaystyle k} & N {\displaystyle N} , Values of d {\displaystyle d} will become ambiguous. This means that specification of a given d {\displaystyle d} could refer to multiple q p M n N {\displaystyle {\vec {q_{p}}}\in M_{n}^{N}} .

Types of connectivity

2-dimensional

Example of neighborhood of pixels - association of eight and four pixels

4-connected

4-connected pixels are neighbors to every pixel that touches one of their edges. These pixels are connected horizontally and vertically. In terms of pixel coordinates, every pixel that has the coordinates

( x ± 1 , y ) {\displaystyle \textstyle (x\pm 1,y)} or ( x , y ± 1 ) {\displaystyle \textstyle (x,y\pm 1)}

is connected to the pixel at ( x , y ) {\displaystyle \textstyle (x,y)} .

6-connected

6-connected pixels are neighbors to every pixel that touches one of their corners (which includes pixels that touch one of their edges) in a hexagonal grid or stretcher bond rectangular grid.

There are several ways to map hexagonal tiles to integer pixel coordinates. With one method, in addition to the 4-connected pixels, the two pixels at coordinates ( x + 1 , y + 1 ) {\displaystyle \textstyle (x+1,y+1)} and ( x 1 , y 1 ) {\displaystyle \textstyle (x-1,y-1)} are connected to the pixel at ( x , y ) {\displaystyle \textstyle (x,y)} .

8-connected

8-connected pixels are neighbors to every pixel that touches one of their edges or corners. These pixels are connected horizontally, vertically, and diagonally. In addition to 4-connected pixels, each pixel with coordinates ( x ± 1 , y ± 1 ) {\displaystyle \textstyle (x\pm 1,y\pm 1)} is connected to the pixel at ( x , y ) {\displaystyle \textstyle (x,y)} .

3-dimensional

6-connected

6-connected pixels are neighbors to every pixel that touches one of their faces. These pixels are connected along one of the primary axes. Each pixel with coordinates ( x ± 1 , y , z ) {\displaystyle \textstyle (x\pm 1,y,z)} , ( x , y ± 1 , z ) {\displaystyle \textstyle (x,y\pm 1,z)} , or ( x , y , z ± 1 ) {\displaystyle \textstyle (x,y,z\pm 1)} is connected to the pixel at ( x , y , z ) {\displaystyle \textstyle (x,y,z)} .

18-connected

18-connected pixels are neighbors to every pixel that touches one of their faces or edges. These pixels are connected along either one or two of the primary axes. In addition to 6-connected pixels, each pixel with coordinates ( x ± 1 , y ± 1 , z ) {\displaystyle \textstyle (x\pm 1,y\pm 1,z)} , ( x ± 1 , y 1 , z ) {\displaystyle \textstyle (x\pm 1,y\mp 1,z)} , ( x ± 1 , y , z ± 1 ) {\displaystyle \textstyle (x\pm 1,y,z\pm 1)} , ( x ± 1 , y , z 1 ) {\displaystyle \textstyle (x\pm 1,y,z\mp 1)} , ( x , y ± 1 , z ± 1 ) {\displaystyle \textstyle (x,y\pm 1,z\pm 1)} , or ( x , y ± 1 , z 1 ) {\displaystyle \textstyle (x,y\pm 1,z\mp 1)} is connected to the pixel at ( x , y , z ) {\displaystyle \textstyle (x,y,z)} .

26-connected

26-connected pixels are neighbors to every pixel that touches one of their faces, edges, or corners. These pixels are connected along either one, two, or all three of the primary axes. In addition to 18-connected pixels, each pixel with coordinates ( x ± 1 , y ± 1 , z ± 1 ) {\displaystyle \textstyle (x\pm 1,y\pm 1,z\pm 1)} , ( x ± 1 , y ± 1 , z 1 ) {\displaystyle \textstyle (x\pm 1,y\pm 1,z\mp 1)} , ( x ± 1 , y 1 , z ± 1 ) {\displaystyle \textstyle (x\pm 1,y\mp 1,z\pm 1)} , or ( x 1 , y ± 1 , z ± 1 ) {\displaystyle \textstyle (x\mp 1,y\pm 1,z\pm 1)} is connected to the pixel at ( x , y , z ) {\displaystyle \textstyle (x,y,z)} .

See also

References

  1. ^ Jonker, Pieter (1992). Morphological Image Processing: Architecture and VLSI design. Kluwer Technische Boeken B.V. pp. 92–96. ISBN 978-1-4615-2804-3.
  • A. Rosenfeld, A. C. Kak (1982), Digital Picture Processing, Academic Press, Inc., ISBN 0-12-597302-0
  • Cheng, CC; Peng, GJ; Hwang, WL (2009), "Subband Weighting With Pixel Connectivity for 3-D Wavelet Coding", IEEE Transactions on Image Processing, 18 (1): 52–62, doi:10.1109/TIP.2008.2007067, PMID 19095518, retrieved 2009-02-16
  • Cheng, CC; Peng, GJ; Hwang, WL (2009), "Subband Weighting With Pixel Connectivity for 3-D Wavelet Coding", IEEE Transactions on Image Processing, 18 (1): 52–62, doi:10.1109/TIP.2008.2007067, PMID 19095518, retrieved 2009-02-16