geometry_utils.h File Reference

a few functions useful in geometry calculations. More...

#include <math.h>
#include <stdlib.h>
#include <math/box2.h>
#include <geometry/eda_angle.h>

Go to the source code of this file.

Classes
class	DISABLE_ARC_RADIUS_CORRECTION
	When creating polygons to create a clearance polygonal area, the polygon must be same or bigger than the original shape. More...

Enumerations
enum	ERROR_LOC { ERROR_OUTSIDE , ERROR_INSIDE }
	When approximating an arc or circle, should the error be placed on the outside or inside of the curve? (Generally speaking filled shape errors go on the inside and knockout errors go on the outside. More...

Functions
int	GetArcToSegmentCount (int aRadius, int aErrorMax, const EDA_ANGLE &aArcAngle)
int	CircleToEndSegmentDeltaRadius (int aInnerCircleRadius, int aSegCount)
int	GetCircleToPolyCorrection (int aMaxError)
template<typename T >
VECTOR2< T >	GetVectorSnapped45 (const VECTOR2< T > &aVec, bool only45=false)
	Snap a vector onto the nearest 0, 45 or 90 degree line. More...
template<typename in_type , typename ret_type = in_type, typename pad_type = unsigned int, typename = typename std::enable_if<std::is_unsigned<pad_type>::value>::type>
VECTOR2< ret_type >	GetClampedCoords (const VECTOR2< in_type > &aCoords, pad_type aPadding=1u)
	Clamps a vector to values that can be negated, respecting numeric limits of coordinates data type with specified padding. More...
bool	ClipLine (const BOX2I *aClipBox, int &x1, int &y1, int &x2, int &y2)
	Test if any part of a line falls within the bounds of a rectangle. More...

Detailed Description

a few functions useful in geometry calculations.

Definition in file geometry_utils.h.

Enumeration Type Documentation

ERROR_LOC

enum ERROR_LOC

When approximating an arc or circle, should the error be placed on the outside or inside of the curve? (Generally speaking filled shape errors go on the inside and knockout errors go on the outside.

This preserves minimum clearances.)

Enumerator
ERROR_OUTSIDE
ERROR_INSIDE

Definition at line 43 of file geometry_utils.h.

{ ERROR_OUTSIDE, ERROR_INSIDE };

Function Documentation

CircleToEndSegmentDeltaRadius()

int CircleToEndSegmentDeltaRadius	(	int	aInnerCircleRadius,
		int	aSegCount
	)

Returns

the radius diffence of the circle defined by segments inside the circle and the radius of the circle tangent to the middle of segments (defined by segments outside this circle)

Parameters

aInnerCircleRadius	is the radius of the circle tangent to the middle of segments
aSegCount	is the seg count to approximate the circle

Definition at line 67 of file geometry_utils.cpp.

{
 // The minimal seg count is 3, otherwise we cannot calculate the result
 // in practice, the min count is clamped to 8 in kicad
 if( aSegCount <= 2 )
 aSegCount = 3;
 
 // The angle between the center of the segment and one end of the segment
 // when the circle is approximated by aSegCount segments
 double alpha = M_PI / aSegCount;
 
 // aRadius is the radius of the circle tangent to the middle of each segment
 // and aRadius/cos(aplha) is the radius of the circle defined by seg ends
 int delta = KiROUND( std::abs( aRadius * ( 1 - 1/cos( alpha ) ) ) );
 
 return delta;
}

References std::abs(), delta, and KiROUND().

Referenced by ConvertArcToPolyline(), SHAPE_ARC::ConvertToPolyline(), CornerListToPolygon(), SCH_SHAPE::Plot(), LIB_SHAPE::Plot(), TransformCircleToPolygon(), TransformOvalToPolygon(), and PAD::TransformShapeToPolygon().

ClipLine()

bool ClipLine	(	const BOX2I *	aClipBox,
		int &	x1,
		int &	y1,
		int &	x2,
		int &	y2
	)

Test if any part of a line falls within the bounds of a rectangle.

Please note that this is only accurate for lines that are one pixel wide.

Parameters

aClipBox	- The rectangle to test.
x1	- X coordinate of one end of a line.
y1	- Y coordinate of one end of a line.
x2	- X coordinate of the other end of a line.
y2	- Y coordinate of the other end of a line.

Returns

- False if any part of the line lies within the rectangle.

Definition at line 137 of file geometry_utils.cpp.

{
 // Stock Cohen-Sutherland algorithm; check *any* CG book for details
 int outcode1 = clipOutCode( aClipBox, x1, y1 );
 int outcode2 = clipOutCode( aClipBox, x2, y2 );
 
 while( outcode1 || outcode2 )
 {
 // Fast reject
 if( outcode1 & outcode2 )
 return true;
 
 // Choose a side to clip
 int thisoutcode, x, y;
 
 if( outcode1 )
 thisoutcode = outcode1;
 else
 thisoutcode = outcode2;
 
 /* One clip round
 * Since we use the full range of 32 bit ints, the proportion
 * computation has to be done in 64 bits to avoid horrible
 * results */
 if( thisoutcode & 1 ) // Clip the bottom
 {
 y = aClipBox->GetBottom();
 x = x1 + (x2 - x1) * std::int64_t(y - y1) / (y2 - y1);
 }
 else if( thisoutcode & 2 ) // Clip the top
 {
 y = aClipBox->GetY();
 x = x1 + ( x2 - x1 ) * std::int64_t( y - y1 ) / ( y2 - y1 );
 }
 else if( thisoutcode & 8 ) // Clip the right
 {
 x = aClipBox->GetRight();
 y = y1 + ( y2 - y1 ) * std::int64_t( x - x1 ) / ( x2 - x1 );
 }
 else // if( thisoutcode & 4), obviously, clip the left
 {
 x = aClipBox->GetX();
 y = y1 + ( y2 - y1 ) * std::int64_t( x - x1 ) / ( x2 - x1 );
 }
 
 // Put the result back and update the boundary code
 // No ambiguity, otherwise it would have been a fast reject
 if( thisoutcode == outcode1 )
 {
 x1 = x;
 y1 = y;
 outcode1 = clipOutCode( aClipBox, x1, y1 );
 }
 else
 {
 x2 = x;
 y2 = y;
 outcode2 = clipOutCode( aClipBox, x2, y2 );
 }
 }
 
 return false;
}

References clipOutCode(), BOX2< Vec >::GetBottom(), BOX2< Vec >::GetRight(), BOX2< Vec >::GetX(), and BOX2< Vec >::GetY().

Referenced by ZONE_FILLER::buildThermalSpokes(), KIGFX::PCB_PAINTER::draw(), STROKE_PARAMS::Stroke(), KIGFX::ORIGIN_VIEWITEM::ViewDraw(), and PCB_TRACK::ViewGetLOD().

GetArcToSegmentCount()

int GetArcToSegmentCount	(	int	aRadius,
		int	aErrorMax,
		const EDA_ANGLE &	aArcAngle
	)

Returns

the number of segments to approximate a arc by segments with a given max error (this number is >= 1)

Parameters

aRadius	is the radius od the circle or arc
aErrorMax	is the max error This is the max distance between the middle of a segment and the circle.
aArcAngleDegree	is the arc angle

Definition at line 42 of file geometry_utils.cpp.

{
 // calculate the number of segments to approximate a circle by segments
 // given the max distance between the middle of a segment and the circle
 
 // avoid divide-by-zero
 aRadius = std::max( 1, aRadius );
 aErrorMax = std::max( 1, aErrorMax );
 
 // error relative to the radius value:
 double rel_error = (double)aErrorMax / aRadius;
 // minimal arc increment in degrees:
 double arc_increment = 180 / M_PI * acos( 1.0 - rel_error ) * 2;
 
 // Ensure a minimal arc increment reasonable value for a circle
 // (360.0 degrees). For very small radius values, this is mandatory.
 arc_increment = std::min( 360.0/MIN_SEGCOUNT_FOR_CIRCLE, arc_increment );
 
 int segCount = KiROUND( fabs( aArcAngle.AsDegrees() ) / arc_increment );
 
 // Ensure at least two segments are used for algorithmic safety
 return std::max( segCount, 2 );
}

References EDA_ANGLE::AsDegrees(), KiROUND(), and MIN_SEGCOUNT_FOR_CIRCLE.

Referenced by PDF_PLOTTER::Arc(), DRC_TEST_PROVIDER_SOLDER_MASK::buildRTrees(), SHAPE_POLY_SET::chamferFilletPolygon(), ConvertArcToPolyline(), ConvertOutlineToPolygon(), SHAPE_ARC::ConvertToPolyline(), CornerListToPolygon(), BOARD_ADAPTER::createTrack(), KIGFX::PCB_PAINTER::draw(), KIGFX::OPENGL_GAL::DrawArcSegment(), EXPORTER_PCB_VRML::ExportVrmlPadHole(), EXPORTER_PCB_VRML::ExportVrmlViaHoles(), ZONE_FILLER::fillCopperZone(), ZONE_FILLER::fillNonCopperZone(), gen_arc(), BOARD_ADAPTER::GetCircleSegmentCount(), EAGLE_PLUGIN::loadPlain(), EAGLE_PLUGIN::loadPolygon(), EAGLE_PLUGIN::loadSignals(), EAGLE_PLUGIN::packageCircle(), EAGLE_PLUGIN::packagePolygon(), PlotSolderMaskLayer(), PlotStandardLayer(), TEARDROP_MANAGER::SetTeardrops(), TransformCircleToPolygon(), TransformOvalToPolygon(), ZONE::TransformShapeToPolygon(), PAD::TransformShapeToPolygon(), and ZONE::TransformSmoothedOutlineToPolygon().

GetCircleToPolyCorrection()

int GetCircleToPolyCorrection

(

int

aMaxError

)

Returns

the radius correction to approximate a circle.

Parameters

aMaxError

is the same error value used to calculate the number of segments.

When creating a polygon from a circle, the polygon is inside the circle. Only corners are on the circle. This is incorrect when building clearance areas of circles, that need to build the equivalent polygon outside the circle.

Definition at line 106 of file geometry_utils.cpp.

{
 // Push all the error to the outside by increasing the radius
 return s_disable_arc_correction ? 0 : aMaxError;
}

References s_disable_arc_correction.

Referenced by CornerListToPolygon(), TransformCircleToPolygon(), TransformOvalToPolygon(), and PAD::TransformShapeToPolygon().

GetClampedCoords()

template<typename in_type , typename ret_type = in_type, typename pad_type = unsigned int, typename = typename std::enable_if<std::is_unsigned<pad_type>::value>::type>

VECTOR2< ret_type > GetClampedCoords	(	const VECTOR2< in_type > &	aCoords,
		pad_type	aPadding = 1u
	)

Clamps a vector to values that can be negated, respecting numeric limits of coordinates data type with specified padding.

Numeric limits are (-2^31 + 1) to (2^31 - 1).

Takes care of rounding in case of floating point to integer conversion.

Parameters

aCoord	- vector to clamp.
aPadding	- padding from the limits. Must not be negative.

Returns

clamped vector.

Definition at line 151 of file geometry_utils.h.

{
 typedef std::numeric_limits<int32_t> coord_limits;
 
 long long max = static_cast<long long>( coord_limits::max() ) - aPadding;
 long long min = -max;
 
 in_type x = aCoords.x;
 in_type y = aCoords.y;
 
 if( x < min )
 x = in_type( min );
 else if( x > max )
 x = in_type( max );
 
 if( y < min )
 y = in_type( min );
 else if( y > max )
 y = in_type( max );
 
 if( !std::is_integral<in_type>() && std::is_integral<ret_type>() )
 return VECTOR2<ret_type>( KiROUND( x ), KiROUND( y ) );
 
 return VECTOR2<ret_type>( x, y );
}

References KiROUND(), VECTOR2< T >::x, and VECTOR2< T >::y.

Referenced by DRAWING_TOOL::drawArc(), DRAWING_TOOL::DrawDimension(), DRAWING_TOOL::drawShape(), DRAWING_TOOL::DrawZone(), KIGFX::WX_VIEW_CONTROLS::ForceCursorPosition(), KIGFX::WX_VIEW_CONTROLS::GetCursorPosition(), KIGFX::WX_VIEW_CONTROLS::GetMousePosition(), EDIT_TOOL::getSafeMovement(), KIGFX::WX_VIEW_CONTROLS::onMotion(), DRAWING_TOOL::PlaceImage(), DRAWING_TOOL::PlaceImportedGraphics(), DRAWING_TOOL::PlaceText(), KIGFX::WX_VIEW_CONTROLS::refreshMouse(), KIGFX::WX_VIEW_CONTROLS::SetCrossHairCursorPosition(), KIGFX::WX_VIEW_CONTROLS::SetCursorPosition(), and KIGFX::WX_VIEW_CONTROLS::WarpMouseCursor().

GetVectorSnapped45()

template<typename T >

VECTOR2< T > GetVectorSnapped45	(	const VECTOR2< T > &	aVec,
		bool	only45 = false
	)

Snap a vector onto the nearest 0, 45 or 90 degree line.

The magnitude of the vector is NOT kept, instead the coordinates are set equal (and/or opposite) or to zero as needed. The effect of this is that if the starting vector is on a square grid, the resulting snapped vector will still be on the same grid.

Parameters

a	vector to be snapped

Returns

the snapped vector

Definition at line 105 of file geometry_utils.h.

{
 using ext_type = typename VECTOR2<T>::extended_type;
 
 auto newVec = aVec;
 const VECTOR2<T> absVec{ std::abs( aVec.x ), std::abs( aVec.y ) };
 
 if( !only45 && absVec.x > ext_type( absVec.y ) * 2 )
 {
 // snap along x-axis
 newVec.y = 0;
 }
 else if( !only45 && absVec.y > ext_type( absVec.x ) * 2 )
 {
 // snap onto y-axis
 newVec.x = 0;
 }
 else if( absVec.x > absVec.y )
 {
 // snap away from x-axis towards 45
 newVec.y = std::copysign( aVec.x, aVec.y );
 }
 else
 {
 // snap away from y-axis towards 45
 newVec.x = std::copysign( aVec.y, aVec.x );
 }
 
 return newVec;
}

References std::abs(), VECTOR2< T >::x, and VECTOR2< T >::y.

Referenced by EC_45DEGREE::Apply(), build45DegLeader(), DRAWING_TOOL::constrainDimension(), EDIT_TOOL::doMoveSelection(), DRAWING_TOOL::drawShape(), and KIGFX::PREVIEW::TWO_POINT_GEOMETRY_MANAGER::SetEnd().