#include <sdet_edgel_regions.h>

Public Types
enum	RegionLabel { UNLABELED = 0, EDGE, LABEL }
Public Member Functions
	sdet_edgel_regions (int array_scale=2, bool verbose=false, bool debug=false)
	~sdet_edgel_regions ()
bool	compute_edgel_regions (vil1_image const &image, vcl_vector< vtol_edge_2d_sptr > &sgrp, vcl_vector< vtol_intensity_face_sptr > &faces)
	The key process loop.
bool	compute_edgel_regions (vil_image_resource_sptr const &image, vcl_vector< vtol_edge_2d_sptr > &sgrp, vcl_vector< vtol_intensity_face_sptr > &faces)
bool	compute_edgel_regions (gevd_bufferxy *buf, vcl_vector< vtol_edge_2d_sptr > &sgrp, vcl_vector< vtol_intensity_face_sptr > &faces)
void	SetVerbose ()
void	ClearVerbose ()
void	SetDebug ()
void	ClearDebug ()
unsigned int	BaseLabel (unsigned int label)
	Get the most basic label equivalent to a given label.
unsigned int	GetMaxRegionLabel () const
void	SetMaxRegionLabel (unsigned int label)
vbl_array_2d< unsigned int > const &	GetRegionArray () const
int	GetXSize () const
int	GetYSize () const
vil1_image	GetEdgeImage (vcl_vector< vtol_edge_2d_sptr > &edgels)
	Return label image (255/0) indicating edgels.
bool	InsertRegionEquivalence (unsigned int label_b, unsigned int label_a)
	assign equivalence of region label b to region label a.
void	GrowEquivalenceClasses ()
	Form equivalence classes by transitive closure on each label.
void	PropagateEquivalence ()
	Populate the label_map_ to reflect the equivalences between labels.
unsigned int	GetLabel (vtol_edge_2d_sptr e, unsigned int nr) const
	Get the a region label for an edge used to construct the boundaries.
void	print_region_array ()
	Print the region label array.
void	print_region_equivalence ()
	Print the contents of the forward equivalence index.
void	print_reverse_region_equivalence ()
	Print the contents of the reverse equivalence index.
void	print_base_equivalence ()
	Print the reduced equivalence relation.
void	print_intensity_data ()
	Print the fitted intensity data for all faces.
Protected Member Functions
bool	GroupContainsEdges (vcl_vector< vtol_edge_2d_sptr > &sg)
	Check if the vtol_edge list sg (spatial group) contains edge(s).
bool	InitRegionArray (vcl_vector< vtol_edge_2d_sptr > &sg)
	Initialize the region label array.
unsigned char	label_code (unsigned int label)
	Get code for a given label.
bool	add_to_forward (unsigned int key, unsigned int value)
	Add a new pair to the forward equivalence list.
bool	add_to_reverse (unsigned int key, unsigned int value)
	Add a new pair to the reverse equivalence list.
unsigned char	EncodeNeighborhood (unsigned int ul, unsigned int ur, unsigned int ll, unsigned int lr)
	Encode a 2x2 neighborhood with the state of the region array for a given location.
void	UpdateConnectedNeighborhood (unsigned int x, unsigned int y)
	Propagate connected components.
void	AssignEdgeLabels (unsigned int x, unsigned int y)
	After connected components have been generated pass over the array and assign region labels to the sdet_region_edge(s).
void	ApplyRegionEquivalence ()
	Scan the region_label_array_ and apply the region equivalence map.
bool	out_of_bounds (unsigned int x, unsigned int y)
	Bounds check on region_label_array_.
void	insert_adjacency (unsigned int region, vtol_edge_2d_sptr e)
	Insert an Edge into the adjacency list for a region.
void	CollectEdges ()
	Get the edges adjacent to each region.
void	CollectFaceEdges ()
	Trace through the topology network keeping regions on the left.
void	ConstructFaces ()
	Construct face(s) from edge(s) in the face_edge_index_ array.
void	AccumulateMeans ()
	Accumulate intensity statistics from each region and update the vtol_intensity_face parameters.
void	AccumulateRegionData ()
	Insert region pixels into the vtol_intensity_face arrays.
void	InsertFaceData ()
	Do both a scatter matrix update and insertion into the region pixel arrays of each intensity face.
unsigned int	X (unsigned int x)
	The region array can be offset from the corner of the ROI.
unsigned int	Y (unsigned int y)
	The region array can be offset from the corner of the ROI.
float	Xf (float x)
	Transforms the image x coordinate to the array coordinate with a scale factor.
float	Yf (float y)
bool	insert_edgel (float pre_x, float pre_y, float x, float y, sdet_region_edge_sptr const &e)
	A utility for inserting an edgel into the region_label_array_.
void	insert_equivalence (unsigned int ll, unsigned int ur, unsigned int &lr)
	This is the fundamental assignment of label equivalence.
bool	merge_equivalence (vcl_map< unsigned int, vcl_vector< unsigned int > * > &tab, unsigned int cur_label, unsigned int label)
	Find the set of labels equivalent to label from a given hash table and merge into equivalence_set_.
bool	get_next_label (vcl_vector< unsigned int > *labels, unsigned int &label)
	Find the next label not accounted for in the current equivalence set.
void	print_edge_colis (unsigned int x, unsigned int y, sdet_region_edge_sptr const &r1, sdet_region_edge_sptr const &r2)
	A collision is defined by the condition where a region is bounded by two different edges at adjacent pixels without crossing a vertex.
bool	remove_hairs (vcl_vector< vtol_edge_2d_sptr > &edges)
	Remove hairs from region boundary.
void	get_buffer_row (unsigned int row)
	get a row from a BufferXY.
void	get_image_row (unsigned int row)
	Get an image row.
unsigned short	get_intensity (unsigned int x)
	Get the intensity of a single pixel.
int	bytes_per_pix ()
bool	compute_edgel_regions (vcl_vector< vtol_edge_2d_sptr > &sgrp, vcl_vector< vtol_intensity_face_sptr > &faces)
Protected Attributes
bool	verbose_
bool	debug_
bool	image_source_
bool	buf_source_
int	s_
vil1_image	image_
gevd_bufferxy *	buf_
vbl_array_2d < sdet_region_edge_sptr >	edge_boundary_array_
vbl_array_2d< unsigned int >	region_label_array_
unsigned int	min_region_label_
unsigned int	max_region_label_
float	Xob_
float	Yob_
unsigned int	xo_
unsigned int	yo_
unsigned int	xend_
unsigned int	yend_
unsigned int	xs_
unsigned int	ys_
vcl_map< unsigned int, vcl_vector< unsigned int > * >	region_pairs_forward_
vcl_map< unsigned int, vcl_vector< unsigned int > * >	region_pairs_reverse_
vcl_map< unsigned int, vcl_vector< unsigned int > * >	equivalence_set_
vcl_map< unsigned int, unsigned int >	label_map_
vcl_map< int, sdet_region_edge_sptr >	region_edges_
vcl_map< unsigned int, vcl_vector< vtol_edge_2d_sptr > * >	region_edge_adjacency_
vcl_vector < vtol_intensity_face_sptr > *	faces_
vtol_intensity_face_sptr *	intensity_face_index_
vcl_vector< vtol_edge_2d_sptr > **	face_edge_index_
vcl_vector< vtol_edge_2d_sptr > *	failed_insertions_
unsigned char *	ubuf_
unsigned short *	sbuf_

Detailed Description

Definition at line 60 of file sdet_edgel_regions.h.

Member Enumeration Documentation

enum sdet_edgel_regions::RegionLabel

Enumerator:

UNLABELED
EDGE
LABEL

Definition at line 63 of file sdet_edgel_regions.h.

Constructor & Destructor Documentation

sdet_edgel_regions::sdet_edgel_regions	(	int	array_scale = `2`,
		bool	verbose = `false`,
		bool	debug = `false`
	)

Definition at line 183 of file sdet_edgel_regions.cxx.

sdet_edgel_regions::~sdet_edgel_regions ( )

Definition at line 213 of file sdet_edgel_regions.cxx.

Member Function Documentation

void sdet_edgel_regions::AccumulateMeans ( ) [protected]

Accumulate intensity statistics from each region and update the vtol_intensity_face parameters.

Definition at line 1742 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::AccumulateRegionData ( ) [protected]

Insert region pixels into the vtol_intensity_face arrays.

Definition at line 1788 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::add_to_forward	(	unsigned int	key,
		unsigned int	value
	)		`[protected]`

Add a new pair to the forward equivalence list.

That is, a ==> b. Note that there can be multiple equivalences which are stored in a vcl_vector

Definition at line 881 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::add_to_reverse	(	unsigned int	key,
		unsigned int	value
	)		`[protected]`

Add a new pair to the reverse equivalence list.

That is, b==>a. Note that there can be multiple equivalences which are stored in a vcl_vector

Definition at line 914 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::ApplyRegionEquivalence ( ) [protected]

Scan the region_label_array_ and apply the region equivalence map.

The result is that all equivalences are reconciled with the smallest labels.

Definition at line 1460 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::AssignEdgeLabels	(	unsigned int	x,
		unsigned int	y
	)		`[protected]`

After connected components have been generated pass over the array and assign region labels to the sdet_region_edge(s).

As in ::UpdateConnectedNeighborhood, the algorithm uses a 2x2 neighborhood E.G., ee But the purpose is to assign labels. No updating of the aa region_label_array_ is carried out.

Note that the 2x2 neighborhood is encoded as uchar = [ul|ur|ll|lr] 2 2 2 2 bits

Definition at line 1315 of file sdet_edgel_regions.cxx.

unsigned int sdet_edgel_regions::BaseLabel ( unsigned int label )

Get the most basic label equivalent to a given label.

If the label_map_ is not defined, this function uses the forward label hash table. Otherwise 0 is returned.

Definition at line 418 of file sdet_edgel_regions.cxx.

int sdet_edgel_regions::bytes_per_pix ( ) [protected]

Definition at line 716 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::ClearDebug ( ) [inline]

Definition at line 85 of file sdet_edgel_regions.h.

void sdet_edgel_regions::ClearVerbose ( ) [inline]

Definition at line 83 of file sdet_edgel_regions.h.

void sdet_edgel_regions::CollectEdges ( ) [protected]

Get the edges adjacent to each region.

Definition at line 1515 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::CollectFaceEdges ( ) [protected]

Trace through the topology network keeping regions on the left.

At this point, we have a list of edges for each region. We will trace out the list and form OneCycle(s). Then the set of OneCycle(s) will form the boundary of the face corresponding to the region.

Definition at line 1544 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::compute_edgel_regions	(	vil1_image const &	image,
		vcl_vector< vtol_edge_2d_sptr > &	sgrp,
		vcl_vector< vtol_intensity_face_sptr > &	faces
	)

The key process loop.

Carries out the steps:

Connected components
Edge-label assignment
Collect region boundaries
Construct vtol_intensity_faces
Calculate intensity fit

Definition at line 271 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::compute_edgel_regions	(	vil_image_resource_sptr const &	image,
		vcl_vector< vtol_edge_2d_sptr > &	sgrp,
		vcl_vector< vtol_intensity_face_sptr > &	faces
	)

Definition at line 295 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::compute_edgel_regions	(	gevd_bufferxy *	buf,
		vcl_vector< vtol_edge_2d_sptr > &	sgrp,
		vcl_vector< vtol_intensity_face_sptr > &	faces
	)

Definition at line 252 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::compute_edgel_regions	(	vcl_vector< vtol_edge_2d_sptr > &	sgrp,
		vcl_vector< vtol_intensity_face_sptr > &	faces
	)		`[protected]`

Definition at line 319 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::ConstructFaces ( ) [protected]

Construct face(s) from edge(s) in the face_edge_index_ array.

This method has been made virtual so that sub-classes of vtol_intensity_face can be constructed by sub-classes of sdet_edgel_regions.

Definition at line 1623 of file sdet_edgel_regions.cxx.

unsigned char sdet_edgel_regions::EncodeNeighborhood	(	unsigned int	ul,
		unsigned int	ur,
		unsigned int	ll,
		unsigned int	lr
	)		`[protected]`

Encode a 2x2 neighborhood with the state of the region array for a given location.

      The Neighborhood    The states are:
           ul ur          UNLABELED, EDGE, LABEL
           ll lr              0       1      2

The encoding maps the 2x2 pattern to an unsigned char. The layout of the uchar is: [ul|ur|ll|lr] with 2 bits for the state of each position.

Definition at line 953 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::get_buffer_row ( unsigned int row ) [protected]

get a row from a BufferXY.

Definition at line 1676 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::get_image_row ( unsigned int row ) [protected]

Get an image row.

Definition at line 1702 of file sdet_edgel_regions.cxx.

unsigned short sdet_edgel_regions::get_intensity ( unsigned int x ) [protected]

Get the intensity of a single pixel.

Definition at line 1723 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::get_next_label	(	vcl_vector< unsigned int > *	labels,
		unsigned int &	label
	)		`[protected]`

Find the next label not accounted for in the current equivalence set.

The set of labels is searched to find a label larger than label, but not in the set, labels.

Definition at line 535 of file sdet_edgel_regions.cxx.

vil1_image sdet_edgel_regions::GetEdgeImage ( vcl_vector< vtol_edge_2d_sptr > & sg )

Return label image (255/0) indicating edgels.

Paint the edgels into the region label array and then output an image where the value is 255 if the pixel is an edge, 0 otherwise

Definition at line 433 of file sdet_edgel_regions.cxx.

unsigned int sdet_edgel_regions::GetLabel	(	vtol_edge_2d_sptr	e,
		unsigned int	nr
	)		const

Get the a region label for an edge used to construct the boundaries.

A return corresponding to UNLABELED means the domain outside the ROI or nr is larger than the number of adjacent regions.

Definition at line 1487 of file sdet_edgel_regions.cxx.

unsigned int sdet_edgel_regions::GetMaxRegionLabel ( ) const [inline]

Definition at line 90 of file sdet_edgel_regions.h.

vbl_array_2d<unsigned int> const& sdet_edgel_regions::GetRegionArray ( ) const [inline]

Definition at line 92 of file sdet_edgel_regions.h.

int sdet_edgel_regions::GetXSize ( ) const [inline]

Definition at line 93 of file sdet_edgel_regions.h.

int sdet_edgel_regions::GetYSize ( ) const [inline]

Definition at line 94 of file sdet_edgel_regions.h.

bool sdet_edgel_regions::GroupContainsEdges ( vcl_vector< vtol_edge_2d_sptr > & sg ) [protected]

Check if the vtol_edge list sg (spatial group) contains edge(s).

Todo:: not yet implemented

Definition at line 657 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::GrowEquivalenceClasses ( )

Form equivalence classes by transitive closure on each label.

The resulting label equivalence is stored in the map, equivalence_set_.

The idea is to add labels to the current equivalence set, cur_set, from either forward or reverse equivalence classes until no new equivalences can be found. The current equivalence class, cur_set, is repeatedly scanned when new labels are added to pick up new equivalence groups. Old equivalence entries are removed from forward and reverse forward and reverse equivalence maps as they are used. When the cur_set is completely closed, i.e. no new labels can be added, then a new label not in cur_set is used to seed a new equivalence class.

Definition at line 576 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::InitRegionArray ( vcl_vector< vtol_edge_2d_sptr > & sg ) [protected]

Initialize the region label array.

There are three types of region label symbols:

UNLABELED - no label has yet been assigned,
EDGE, the existence of an edgel boundary pixel.
An unsigned integer which represents an existing region label.

Definition at line 736 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::insert_adjacency	(	unsigned int	region,
		vtol_edge_2d_sptr	e
	)		`[protected]`

Insert an Edge into the adjacency list for a region.

Definition at line 1497 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::insert_edgel	(	float	pre_x,
		float	pre_y,
		float	xedgel,
		float	yedgel,
		sdet_region_edge_sptr const &	re
	)		`[protected]`

A utility for inserting an edgel into the region_label_array_.

An edgel and a previous edgel in the chain are used to interpolate intermediate edgels to take account of pixel quantization

Definition at line 672 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::insert_equivalence	(	unsigned int	ll,
		unsigned int	ur,
		unsigned int &	lr
	)		`[protected]`

This is the fundamental assignment of label equivalence.

Definition at line 978 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::InsertFaceData ( ) [protected]

Do both a scatter matrix update and insertion into the region pixel arrays of each intensity face.

AccumulateScatterData is done first so that the number of pixels can be determined.

Definition at line 1827 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::InsertRegionEquivalence	(	unsigned int	label_b,
		unsigned int	label_a
	)

assign equivalence of region label b to region label a.

Definition at line 405 of file sdet_edgel_regions.cxx.

unsigned char sdet_edgel_regions::label_code ( unsigned int label ) [protected]

Get code for a given label.

Definition at line 867 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::merge_equivalence	(	vcl_map< unsigned int, vcl_vector< unsigned int > * > &	tab,
		unsigned int	cur_label,
		unsigned int	label
	)		`[protected]`

Find the set of labels equivalent to label from a given hash table and merge into equivalence_set_.

cur_label is the equivalence set being formed. label is the table key of equivalences to be merged. The labels equivalent to label are in the input map, tab.

Definition at line 475 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::out_of_bounds	(	unsigned int	x,
		unsigned int	y
	)		`[protected]`

Bounds check on region_label_array_.

Definition at line 1478 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::print_base_equivalence ( )

Print the reduced equivalence relation.

Definition at line 127 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::print_edge_colis	(	unsigned int	x,
		unsigned int	y,
		sdet_region_edge_sptr const &	r1,
		sdet_region_edge_sptr const &	r2
	)		`[protected]`

A collision is defined by the condition where a region is bounded by two different edges at adjacent pixels without crossing a vertex.

This can happen since boundary positions are sub-pixel and region definition is at pixel granularity. The edge collision causes a needed edge to be superseded.

Definition at line 1203 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::print_intensity_data ( )

Print the fitted intensity data for all faces.

Definition at line 137 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::print_region_array ( )

Print the region label array.

Definition at line 65 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::print_region_equivalence ( )

Print the contents of the forward equivalence index.

Definition at line 83 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::print_reverse_region_equivalence ( )

Print the contents of the reverse equivalence index.

Definition at line 105 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::PropagateEquivalence ( )

Populate the label_map_ to reflect the equivalences between labels.

Definition at line 451 of file sdet_edgel_regions.cxx.

bool sdet_edgel_regions::remove_hairs ( vcl_vector< vtol_edge_2d_sptr > & edges ) [protected]

Remove hairs from region boundary.

This condition can occur because the region labels are not sub-pixel. A "hair" is an extra edge joined at a vertex which is part of a continuous chain. Unattached vertices are detected by incrementing a count at each vertex for each attached edge in the input array, "edges". Such hairs are removed from the input array.

Definition at line 1255 of file sdet_edgel_regions.cxx.

void sdet_edgel_regions::SetDebug ( ) [inline]

Definition at line 84 of file sdet_edgel_regions.h.

void sdet_edgel_regions::SetMaxRegionLabel ( unsigned int label ) [inline]

Definition at line 91 of file sdet_edgel_regions.h.

void sdet_edgel_regions::SetVerbose ( ) [inline]

Definition at line 82 of file sdet_edgel_regions.h.