boofcv.alg.feature.detect.edge.GradientToEdgeFeatures Maven / Gradle / Ivy
Go to download
Show more of this group Show more artifacts with this name
Show all versions of boofcv-feature Show documentation
Show all versions of boofcv-feature Show documentation
BoofCV is an open source Java library for real-time computer vision and robotics applications.
/*
* Copyright (c) 2011-2019, Peter Abeles. All Rights Reserved.
*
* This file is part of BoofCV (http://boofcv.org).
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
package boofcv.alg.feature.detect.edge;
import boofcv.alg.InputSanityCheck;
import boofcv.alg.feature.detect.edge.impl.*;
import boofcv.concurrency.BoofConcurrency;
import boofcv.struct.image.GrayF32;
import boofcv.struct.image.GrayS16;
import boofcv.struct.image.GrayS32;
import boofcv.struct.image.GrayS8;
/**
*
* Give the image's gradient in the x and y direction compute the edge's intensity and orientation.
* Two ways are provided for computing the edge's intensity: euclidean norm and sum of absolute values (induced 1-norm).
* The former is the most accurate, while the later is much faster.
*
*
*
* norm: sqrt( gx2 + gy2)
* abs: |gx| + |gy|
* angle: atan( gy / gx )
*
*
*
* When computing the angle care is taken to avoid divided by zero errors.
*
*
* WARNING: Do not modify. Automatically generated by GenerateGradientToEdgeFeatures.
*
*
* @author Peter Abeles
*/
@SuppressWarnings("Duplicates")
public class GradientToEdgeFeatures {
/**
* Computes the edge intensity using a Euclidean norm.
*
* @param derivX Derivative along x-axis. Not modified.
* @param derivY Derivative along y-axis. Not modified.
* @param intensity Edge intensity.
*/
static public void intensityE(GrayF32 derivX , GrayF32 derivY , GrayF32 intensity )
{
InputSanityCheck.checkSameShape(derivX,derivY);
intensity.reshape(derivX.width,derivX.height);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplGradientToEdgeFeatures_MT.intensityE(derivX,derivY,intensity);
} else {
ImplGradientToEdgeFeatures.intensityE(derivX,derivY,intensity);
}
}
/**
* Computes the edge intensity using a Euclidean norm.
*
* @param derivX Derivative along x-axis. Not modified.
* @param derivY Derivative along y-axis. Not modified.
* @param intensity Edge intensity.
*/
static public void intensityAbs(GrayF32 derivX , GrayF32 derivY , GrayF32 intensity )
{
InputSanityCheck.checkSameShape(derivX,derivY);
intensity.reshape(derivX.width,derivX.height);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplGradientToEdgeFeatures_MT.intensityAbs(derivX,derivY,intensity);
} else {
ImplGradientToEdgeFeatures.intensityAbs(derivX,derivY,intensity);
}
}
/**
* Computes the edge orientation using the {@link Math#atan} function.
*
* @param derivX Derivative along x-axis. Not modified.
* @param derivY Derivative along y-axis. Not modified.
* @param angle Edge orientation in radians (-pi/2 to pi/2).
*/
static public void direction(GrayF32 derivX , GrayF32 derivY , GrayF32 angle )
{
InputSanityCheck.checkSameShape(derivX,derivY);
angle.reshape(derivX.width,derivX.height);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplGradientToEdgeFeatures_MT.direction(derivX,derivY,angle);
} else {
ImplGradientToEdgeFeatures.direction(derivX,derivY,angle);
}
}
/**
* Computes the edge orientation using the {@link Math#atan2} function.
*
* @param derivX Derivative along x-axis. Not modified.
* @param derivY Derivative along y-axis. Not modified.
* @param angle Edge orientation in radians (-pi to pi).
*/
static public void direction2(GrayF32 derivX , GrayF32 derivY , GrayF32 angle )
{
InputSanityCheck.checkSameShape(derivX,derivY);
angle.reshape(derivX.width,derivX.height);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplGradientToEdgeFeatures_MT.direction2(derivX,derivY,angle);
} else {
ImplGradientToEdgeFeatures.direction2(derivX,derivY,angle);
}
}
/**
*
* Sets edge intensities to zero if the pixel has an intensity which is less than any of
* the two adjacent pixels. Pixel adjacency is determined based upon the sign of the image gradient. Less precise
* than other methods, but faster.
*
*
* @param intensity Edge intensities. Not modified.
* @param derivX Image derivative along x-axis.
* @param derivY Image derivative along y-axis.
* @param output Filtered intensity. If null a new image will be declared and returned. Modified.
* @return Filtered edge intensity.
*/
static public GrayF32 nonMaxSuppressionCrude4(GrayF32 intensity , GrayF32 derivX , GrayF32 derivY, GrayF32 output )
{
InputSanityCheck.checkSameShape(intensity,derivX,derivY);
output = InputSanityCheck.checkDeclare(intensity,output);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplEdgeNonMaxSuppressionCrude_MT.inner4(intensity, derivX, derivY, output);
ImplEdgeNonMaxSuppressionCrude_MT.border4(intensity, derivX, derivY, output);
} else {
ImplEdgeNonMaxSuppressionCrude.inner4(intensity, derivX, derivY, output);
ImplEdgeNonMaxSuppressionCrude.border4(intensity, derivX, derivY, output);
}
return output;
}
/**
* Computes the edge intensity using a Euclidean norm.
*
* @param derivX Derivative along x-axis. Not modified.
* @param derivY Derivative along y-axis. Not modified.
* @param intensity Edge intensity.
*/
static public void intensityE(GrayS16 derivX , GrayS16 derivY , GrayF32 intensity )
{
InputSanityCheck.checkSameShape(derivX,derivY);
intensity.reshape(derivX.width,derivX.height);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplGradientToEdgeFeatures_MT.intensityE(derivX,derivY,intensity);
} else {
ImplGradientToEdgeFeatures.intensityE(derivX,derivY,intensity);
}
}
/**
* Computes the edge intensity using a Euclidean norm.
*
* @param derivX Derivative along x-axis. Not modified.
* @param derivY Derivative along y-axis. Not modified.
* @param intensity Edge intensity.
*/
static public void intensityAbs(GrayS16 derivX , GrayS16 derivY , GrayF32 intensity )
{
InputSanityCheck.checkSameShape(derivX,derivY);
intensity.reshape(derivX.width,derivX.height);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplGradientToEdgeFeatures_MT.intensityAbs(derivX,derivY,intensity);
} else {
ImplGradientToEdgeFeatures.intensityAbs(derivX,derivY,intensity);
}
}
/**
* Computes the edge orientation using the {@link Math#atan} function.
*
* @param derivX Derivative along x-axis. Not modified.
* @param derivY Derivative along y-axis. Not modified.
* @param angle Edge orientation in radians (-pi/2 to pi/2).
*/
static public void direction(GrayS16 derivX , GrayS16 derivY , GrayF32 angle )
{
InputSanityCheck.checkSameShape(derivX,derivY);
angle.reshape(derivX.width,derivX.height);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplGradientToEdgeFeatures_MT.direction(derivX,derivY,angle);
} else {
ImplGradientToEdgeFeatures.direction(derivX,derivY,angle);
}
}
/**
* Computes the edge orientation using the {@link Math#atan2} function.
*
* @param derivX Derivative along x-axis. Not modified.
* @param derivY Derivative along y-axis. Not modified.
* @param angle Edge orientation in radians (-pi to pi).
*/
static public void direction2(GrayS16 derivX , GrayS16 derivY , GrayF32 angle )
{
InputSanityCheck.checkSameShape(derivX,derivY);
angle.reshape(derivX.width,derivX.height);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplGradientToEdgeFeatures_MT.direction2(derivX,derivY,angle);
} else {
ImplGradientToEdgeFeatures.direction2(derivX,derivY,angle);
}
}
/**
*
* Sets edge intensities to zero if the pixel has an intensity which is less than any of
* the two adjacent pixels. Pixel adjacency is determined based upon the sign of the image gradient. Less precise
* than other methods, but faster.
*
*
* @param intensity Edge intensities. Not modified.
* @param derivX Image derivative along x-axis.
* @param derivY Image derivative along y-axis.
* @param output Filtered intensity. If null a new image will be declared and returned. Modified.
* @return Filtered edge intensity.
*/
static public GrayF32 nonMaxSuppressionCrude4(GrayF32 intensity , GrayS16 derivX , GrayS16 derivY, GrayF32 output )
{
InputSanityCheck.checkSameShape(intensity,derivX,derivY);
output = InputSanityCheck.checkDeclare(intensity,output);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplEdgeNonMaxSuppressionCrude_MT.inner4(intensity, derivX, derivY, output);
ImplEdgeNonMaxSuppressionCrude_MT.border4(intensity, derivX, derivY, output);
} else {
ImplEdgeNonMaxSuppressionCrude.inner4(intensity, derivX, derivY, output);
ImplEdgeNonMaxSuppressionCrude.border4(intensity, derivX, derivY, output);
}
return output;
}
/**
* Computes the edge intensity using a Euclidean norm.
*
* @param derivX Derivative along x-axis. Not modified.
* @param derivY Derivative along y-axis. Not modified.
* @param intensity Edge intensity.
*/
static public void intensityE(GrayS32 derivX , GrayS32 derivY , GrayF32 intensity )
{
InputSanityCheck.checkSameShape(derivX,derivY);
intensity.reshape(derivX.width,derivX.height);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplGradientToEdgeFeatures_MT.intensityE(derivX,derivY,intensity);
} else {
ImplGradientToEdgeFeatures.intensityE(derivX,derivY,intensity);
}
}
/**
* Computes the edge intensity using a Euclidean norm.
*
* @param derivX Derivative along x-axis. Not modified.
* @param derivY Derivative along y-axis. Not modified.
* @param intensity Edge intensity.
*/
static public void intensityAbs(GrayS32 derivX , GrayS32 derivY , GrayF32 intensity )
{
InputSanityCheck.checkSameShape(derivX,derivY);
intensity.reshape(derivX.width,derivX.height);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplGradientToEdgeFeatures_MT.intensityAbs(derivX,derivY,intensity);
} else {
ImplGradientToEdgeFeatures.intensityAbs(derivX,derivY,intensity);
}
}
/**
* Computes the edge orientation using the {@link Math#atan} function.
*
* @param derivX Derivative along x-axis. Not modified.
* @param derivY Derivative along y-axis. Not modified.
* @param angle Edge orientation in radians (-pi/2 to pi/2).
*/
static public void direction(GrayS32 derivX , GrayS32 derivY , GrayF32 angle )
{
InputSanityCheck.checkSameShape(derivX,derivY);
angle.reshape(derivX.width,derivX.height);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplGradientToEdgeFeatures_MT.direction(derivX,derivY,angle);
} else {
ImplGradientToEdgeFeatures.direction(derivX,derivY,angle);
}
}
/**
* Computes the edge orientation using the {@link Math#atan2} function.
*
* @param derivX Derivative along x-axis. Not modified.
* @param derivY Derivative along y-axis. Not modified.
* @param angle Edge orientation in radians (-pi to pi).
*/
static public void direction2(GrayS32 derivX , GrayS32 derivY , GrayF32 angle )
{
InputSanityCheck.checkSameShape(derivX,derivY);
angle.reshape(derivX.width,derivX.height);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplGradientToEdgeFeatures_MT.direction2(derivX,derivY,angle);
} else {
ImplGradientToEdgeFeatures.direction2(derivX,derivY,angle);
}
}
/**
*
* Sets edge intensities to zero if the pixel has an intensity which is less than any of
* the two adjacent pixels. Pixel adjacency is determined based upon the sign of the image gradient. Less precise
* than other methods, but faster.
*
*
* @param intensity Edge intensities. Not modified.
* @param derivX Image derivative along x-axis.
* @param derivY Image derivative along y-axis.
* @param output Filtered intensity. If null a new image will be declared and returned. Modified.
* @return Filtered edge intensity.
*/
static public GrayF32 nonMaxSuppressionCrude4(GrayF32 intensity , GrayS32 derivX , GrayS32 derivY, GrayF32 output )
{
InputSanityCheck.checkSameShape(intensity,derivX,derivY);
output = InputSanityCheck.checkDeclare(intensity,output);
if(BoofConcurrency.USE_CONCURRENT ) {
ImplEdgeNonMaxSuppressionCrude_MT.inner4(intensity, derivX, derivY, output);
ImplEdgeNonMaxSuppressionCrude_MT.border4(intensity, derivX, derivY, output);
} else {
ImplEdgeNonMaxSuppressionCrude.inner4(intensity, derivX, derivY, output);
ImplEdgeNonMaxSuppressionCrude.border4(intensity, derivX, derivY, output);
}
return output;
}
/**
*
* Converts an image containing edge angles (-pi/2 to pi/2) into a discrete set of angles.
* The conversion is done by rounding the angle to the nearest orientation in the set.
*
*
* Discrete value to angle (degrees): 0=0,1=45,2=90,-1=-45
*
*
* @param angle Input image containing edge orientations. Orientations are assumed to be
* from -pi/2 to pi/2. Not modified.
* @param discrete Output set of discretized angles. Values will be from -1 to 2, inclusive. If null a new
* image will be declared and returned. Modified.
* @return Discretized direction.
*/
static public GrayS8 discretizeDirection4(GrayF32 angle , GrayS8 discrete )
{
discrete = InputSanityCheck.checkDeclare(angle,discrete,GrayS8.class);
final float A = (float)(Math.PI/8.0);
final float B = (float)(Math.PI/4.0);
final int w = angle.width;
final int h = angle.height;
for( int y = 0; y < h; y++ ) {
int indexSrc = angle.startIndex + y*angle.stride;
int indexDst = discrete.startIndex + y*discrete.stride;
int end = indexSrc + w;
for( ; indexSrc < end; indexSrc++ , indexDst++ ) {
float a = angle.data[indexSrc];
int val;
if( a >= 0 ) {
val = (int)((a+A)/B);
} else {
val = (int)((a-A)/B);
}
discrete.data[indexDst] = (byte)(val == -2 ? 2 : val);
}
}
return discrete;
}
/**
*
* Converts an image containing edge angles (-pi to pi) into a discrete set of 8 angles.
* The conversion is done by rounding the angle to the nearest orientation in the set.
*
*
* Discrete value to angle (degrees): 0=0,1=45,2=90,3=135,4=180,-1=-45,-2=--90,-3=-135
*
*
* @param angle Input image containing edge orientations. Orientations are assumed to be
* from -pi to pi. Not modified.
* @param discrete Output set of discretized angles. Values will be from -3 to 4, inclusive. If null a new
* image will be declared and returned. Modified.
* @return Discretized direction.
*/
static public GrayS8 discretizeDirection8(GrayF32 angle , GrayS8 discrete )
{
discrete = InputSanityCheck.checkDeclare(angle,discrete,GrayS8.class);
final float A = (float)(Math.PI/8.0);
final float B = (float)(Math.PI/4.0);
final int w = angle.width;
final int h = angle.height;
for( int y = 0; y < h; y++ ) {
int indexSrc = angle.startIndex + y*angle.stride;
int indexDst = discrete.startIndex + y*discrete.stride;
int end = indexSrc + w;
for( ; indexSrc < end; indexSrc++ , indexDst++ ) {
float a = angle.data[indexSrc];
int val;
if( a >= 0 ) {
val = (int)((a+A)/B);
} else {
val = (int)((a-A)/B);
}
discrete.data[indexDst] = (byte)(val == -4 ? 4 : val);
}
}
return discrete;
}
/**
*
* Sets edge intensities to zero if the pixel has an intensity which is less than either of
* the two adjacent pixels. Pixel adjacency is determined by the gradients discretized direction.
*
*
* @param intensity Edge intensities. Not modified.
* @param direction 4-Discretized direction. See {@link #discretizeDirection4(GrayF32, GrayS8)}. Not modified.
* @param output Filtered intensity. If null a new image will be declared and returned. Modified.
* @return Filtered edge intensity.
*/
static public GrayF32 nonMaxSuppression4(GrayF32 intensity , GrayS8 direction , GrayF32 output )
{
InputSanityCheck.checkSameShape(intensity,direction);
output = InputSanityCheck.checkDeclare(intensity,output);
if( BoofConcurrency.USE_CONCURRENT ) {
ImplEdgeNonMaxSuppression_MT.inner4(intensity, direction, output);
ImplEdgeNonMaxSuppression_MT.border4(intensity, direction, output);
} else {
ImplEdgeNonMaxSuppression.inner4(intensity, direction, output);
ImplEdgeNonMaxSuppression.border4(intensity, direction, output);
}
return output;
}
/**
*
* Sets edge intensities to zero if the pixel has an intensity which is less than either of
* the two adjacent pixels. Pixel adjacency is determined by the gradients discretized direction.
*
*
* @param intensity Edge intensities. Not modified.
* @param direction 8-Discretized direction. See {@link #discretizeDirection8(GrayF32, GrayS8)}. Not modified.
* @param output Filtered intensity. If null a new image will be declared and returned. Modified.
* @return Filtered edge intensity.
*/
static public GrayF32 nonMaxSuppression8(GrayF32 intensity , GrayS8 direction , GrayF32 output )
{
InputSanityCheck.checkSameShape(intensity,direction);
output = InputSanityCheck.checkDeclare(intensity,output);
if( BoofConcurrency.USE_CONCURRENT ) {
ImplEdgeNonMaxSuppression_MT.inner8(intensity, direction, output);
ImplEdgeNonMaxSuppression_MT.border8(intensity, direction, output);
} else {
ImplEdgeNonMaxSuppression.inner8(intensity, direction, output);
ImplEdgeNonMaxSuppression.border8(intensity, direction, output);
}
return output;
}
}
© 2015 - 2025 Weber Informatics LLC | Privacy Policy