mirror of
https://github.com/qgis/QGIS.git
synced 2025-02-28 00:17:30 -05:00
e2e47d71c6
Reimplement translate in terms of transform(QTransform). Includes new unit test for rotate() and translate() methods. Includes sip bindings update.
480 lines
20 KiB
Plaintext
480 lines
20 KiB
Plaintext
|
|
/** polyline is represented as a vector of points */
|
|
typedef QVector<QgsPoint> QgsPolyline;
|
|
|
|
/** polygon: first item of the list is outer ring, inner rings (if any) start from second item */
|
|
typedef QVector< QVector<QgsPoint> > QgsPolygon;
|
|
|
|
/** a collection of QgsPoints that share a common collection of attributes */
|
|
typedef QVector<QgsPoint> QgsMultiPoint;
|
|
|
|
/** a collection of QgsPolylines that share a common collection of attributes */
|
|
typedef QVector< QVector<QgsPoint> > QgsMultiPolyline;
|
|
|
|
/** a collection of QgsPolygons that share a common collection of attributes */
|
|
typedef QVector< QVector< QVector<QgsPoint> > > QgsMultiPolygon;
|
|
|
|
%If (WS_MACX)
|
|
typedef unsigned long size_t;
|
|
%End
|
|
%If (WS_X11 || WS_WIN)
|
|
typedef unsigned int size_t;
|
|
%End
|
|
|
|
class QgsGeometry
|
|
{
|
|
%TypeHeaderCode
|
|
#include <qgsgeometry.h>
|
|
%End
|
|
|
|
public:
|
|
//! Constructor
|
|
QgsGeometry();
|
|
|
|
/** copy constructor will prompt a deep copy of the object */
|
|
QgsGeometry( const QgsGeometry & );
|
|
|
|
|
|
//! Destructor
|
|
~QgsGeometry();
|
|
|
|
/** static method that creates geometry from Wkt */
|
|
static QgsGeometry* fromWkt( QString wkt ) /Factory/;
|
|
|
|
/** construct geometry from a point */
|
|
static QgsGeometry* fromPoint( const QgsPoint& point ) /Factory/;
|
|
/** construct geometry from a multipoint */
|
|
static QgsGeometry* fromMultiPoint( const QgsMultiPoint& multipoint ) /Factory/;
|
|
/** construct geometry from a polyline */
|
|
static QgsGeometry* fromPolyline( const QgsPolyline& polyline ) /Factory/;
|
|
/** construct geometry from a multipolyline*/
|
|
static QgsGeometry* fromMultiPolyline( const QgsMultiPolyline& multiline ) /Factory/;
|
|
/** construct geometry from a polygon */
|
|
static QgsGeometry* fromPolygon( const QgsPolygon& polygon ) /Factory/;
|
|
/** construct geometry from a multipolygon */
|
|
static QgsGeometry* fromMultiPolygon( const QgsMultiPolygon& multipoly ) /Factory/;
|
|
/** construct geometry from a rectangle */
|
|
static QgsGeometry* fromRect( const QgsRectangle& rect ) /Factory/;
|
|
/**Construct geometry from a QPointF
|
|
* @param point source QPointF
|
|
* @note added in QGIS 2.7
|
|
*/
|
|
static QgsGeometry* fromQPointF( const QPointF& point ) /Factory/;
|
|
|
|
/**Construct geometry from a QPolygonF. If the polygon is closed than
|
|
* the resultant geometry will be a polygon, if it is open than the
|
|
* geometry will be a polyline.
|
|
* @param polygon source QPolygonF
|
|
* @note added in QGIS 2.7
|
|
*/
|
|
static QgsGeometry* fromQPolygonF( const QPolygonF& polygon ) /Factory/;
|
|
|
|
/**
|
|
Set the geometry, feeding in a geometry in GEOS format.
|
|
This class will take ownership of the buffer.
|
|
@note not available in python bindings
|
|
*/
|
|
// void fromGeos( GEOSGeometry* geos );
|
|
/**
|
|
Set the geometry, feeding in the buffer containing OGC Well-Known Binary and the buffer's length.
|
|
This class will take ownership of the buffer.
|
|
*/
|
|
void fromWkb( unsigned char * wkb /Array/, size_t length /ArraySize/ );
|
|
%MethodCode
|
|
// create copy of Python's string and pass it to fromWkb()
|
|
unsigned char * copy = new unsigned char[a1];
|
|
memcpy(copy, a0, a1);
|
|
sipCpp->fromWkb(copy, a1);
|
|
%End
|
|
|
|
/**
|
|
Returns the buffer containing this geometry in WKB format.
|
|
You may wish to use in conjunction with wkbSize().
|
|
*/
|
|
SIP_PYOBJECT asWkb();
|
|
%MethodCode
|
|
sipRes = PyString_FromStringAndSize((const char *)sipCpp->asWkb(), sipCpp->wkbSize());
|
|
%End
|
|
|
|
/**
|
|
* Returns the size of the WKB in asWkb().
|
|
*/
|
|
size_t wkbSize() const;
|
|
|
|
/**Returns a geos geometry. QgsGeometry keeps ownership, don't delete the returned object!
|
|
@note not available in python bindings
|
|
*/
|
|
// const GEOSGeometry* asGeos() const;
|
|
|
|
/** Returns type of wkb (point / linestring / polygon etc.) */
|
|
QGis::WkbType wkbType() const;
|
|
|
|
/** Returns type of the vector */
|
|
QGis::GeometryType type();
|
|
|
|
/** Returns true if wkb of the geometry is of WKBMulti* type */
|
|
bool isMultipart();
|
|
|
|
/** compare geometries using GEOS */
|
|
bool isGeosEqual( QgsGeometry & );
|
|
|
|
/** check validity using GEOS */
|
|
bool isGeosValid();
|
|
|
|
/** check if geometry is empty using GEOS */
|
|
bool isGeosEmpty();
|
|
|
|
/** get area of geometry using GEOS */
|
|
double area();
|
|
|
|
/** get length of geometry using GEOS */
|
|
double length();
|
|
|
|
double distance( QgsGeometry& geom );
|
|
|
|
/**
|
|
Returns the vertex closest to the given point, the corresponding vertex index, squared distance snap point / target point
|
|
and the indices of the vertices before/after. The vertices before/after are -1 if not present
|
|
*/
|
|
QgsPoint closestVertex( const QgsPoint& point, int& atVertex /Out/, int& beforeVertex /Out/, int& afterVertex /Out/, double& sqrDist /Out/ );
|
|
|
|
/**
|
|
Returns the indexes of the vertices before and after the given vertex index.
|
|
|
|
This function takes into account the following factors:
|
|
|
|
1. If the given vertex index is at the end of a linestring,
|
|
the adjacent index will be -1 (for "no adjacent vertex")
|
|
2. If the given vertex index is at the end of a linear ring
|
|
(such as in a polygon), the adjacent index will take into
|
|
account the first vertex is equal to the last vertex (and will
|
|
skip equal vertex positions).
|
|
*/
|
|
void adjacentVertices( int atVertex, int& beforeVertex /Out/, int& afterVertex /Out/ );
|
|
|
|
/** Insert a new vertex before the given vertex index,
|
|
* ring and item (first number is index 0)
|
|
* If the requested vertex number (beforeVertex.back()) is greater
|
|
* than the last actual vertex on the requested ring and item,
|
|
* it is assumed that the vertex is to be appended instead of inserted.
|
|
* Returns false if atVertex does not correspond to a valid vertex
|
|
* on this geometry (including if this geometry is a Point).
|
|
* It is up to the caller to distinguish between
|
|
* these error conditions. (Or maybe we add another method to this
|
|
* object to help make the distinction?)
|
|
*/
|
|
bool insertVertex( double x, double y, int beforeVertex );
|
|
|
|
/** Moves the vertex at the given position number
|
|
* and item (first number is index 0)
|
|
* to the given coordinates.
|
|
* Returns false if atVertex does not correspond to a valid vertex
|
|
* on this geometry
|
|
*/
|
|
bool moveVertex( double x, double y, int atVertex );
|
|
|
|
/** Deletes the vertex at the given position number and item
|
|
* (first number is index 0)
|
|
* Returns false if atVertex does not correspond to a valid vertex
|
|
* on this geometry (including if this geometry is a Point),
|
|
* or if the number of remaining verticies in the linestring
|
|
* would be less than two.
|
|
* It is up to the caller to distinguish between
|
|
* these error conditions. (Or maybe we add another method to this
|
|
* object to help make the distinction?)
|
|
*/
|
|
bool deleteVertex( int atVertex );
|
|
|
|
/**
|
|
* Returns coordinates of a vertex.
|
|
* @param atVertex index of the vertex
|
|
* @return Coordinates of the vertex or QgsPoint(0,0) on error
|
|
*/
|
|
QgsPoint vertexAt( int atVertex );
|
|
|
|
/**
|
|
* Returns the squared cartesian distance between the given point
|
|
* to the given vertex index (vertex at the given position number,
|
|
* ring and item (first number is index 0))
|
|
*/
|
|
double sqrDistToVertexAt( QgsPoint& point /In/, int atVertex );
|
|
|
|
/**
|
|
* Searches for the closest vertex in this geometry to the given point.
|
|
* @param point Specifiest the point for search
|
|
* @param atVertex Receives index of the closest vertex
|
|
* @return The squared cartesian distance is also returned in sqrDist, negative number on error
|
|
*/
|
|
double closestVertexWithContext( const QgsPoint& point, int& atVertex /Out/ );
|
|
|
|
/**
|
|
* Searches for the closest segment of geometry to the given point
|
|
* @param point Specifies the point for search
|
|
* @param minDistPoint Receives the nearest point on the segment
|
|
* @param afterVertex Receives index of the vertex after the closest segment. The vertex
|
|
* before the closest segment is always afterVertex - 1
|
|
* @param leftOf Out: Returns if the point lies on the left of right side of the segment ( < 0 means left, > 0 means right )
|
|
* @param epsilon epsilon for segment snapping
|
|
* @return The squared cartesian distance is also returned in sqrDist, negative number on error
|
|
*/
|
|
double closestSegmentWithContext( const QgsPoint& point, QgsPoint& minDistPoint /Out/, int& afterVertex /Out/ );
|
|
|
|
/**Adds a new ring to this geometry. This makes only sense for polygon and multipolygons.
|
|
@return 0 in case of success (ring added), 1 problem with geometry type, 2 ring not closed,
|
|
3 ring is not valid geometry, 4 ring not disjoint with existing rings, 5 no polygon found which contained the ring*/
|
|
int addRing( const QList<QgsPoint>& ring );
|
|
|
|
/**Adds a new island polygon to a multipolygon feature
|
|
@return 0 in case of success, 1 if not a multipolygon, 2 if ring is not a valid geometry, 3 if new polygon ring
|
|
not disjoint with existing polygons of the feature*/
|
|
int addPart( const QList<QgsPoint> &points, QGis::GeometryType geomType = QGis::UnknownGeometry );
|
|
|
|
/**Adds a new island polygon to a multipolygon feature
|
|
@return 0 in case of success, 1 if not a multipolygon, 2 if ring is not a valid geometry, 3 if new polygon ring
|
|
not disjoint with existing polygons of the feature
|
|
@note not available in python bindings
|
|
*/
|
|
// int addPart( GEOSGeometry *newPart );
|
|
|
|
/**Adds a new island polygon to a multipolygon feature
|
|
@return 0 in case of success, 1 if not a multipolygon, 2 if ring is not a valid geometry, 3 if new polygon ring
|
|
not disjoint with existing polygons of the feature
|
|
@note available in python bindings as addPartGeometry (added in 2.1)
|
|
*/
|
|
int addPart( QgsGeometry *newPart /Transfer/ ) /PyName=addPartGeometry/;
|
|
|
|
/**Translate this geometry by dx, dy
|
|
@return 0 in case of success*/
|
|
int translate( double dx, double dy );
|
|
|
|
/**Transform this geometry as described by CoordinateTranasform ct
|
|
@return 0 in case of success*/
|
|
int transform( const QgsCoordinateTransform& ct );
|
|
|
|
/**Transform this geometry as described by QTransform ct
|
|
@note added in 2.8
|
|
@return 0 in case of success*/
|
|
int transform( const QTransform& ct );
|
|
|
|
/**Rotate this geometry around the Z axis
|
|
@note added in 2.8
|
|
@param rotation clockwise rotation in degrees
|
|
@param center rotation center
|
|
@return 0 in case of success*/
|
|
int rotate( double rotation, const QgsPoint& center );
|
|
|
|
/**Splits this geometry according to a given line. Note that the geometry is only split once. If there are several intersections
|
|
between geometry and splitLine, only the first one is considered.
|
|
@param splitLine the line that splits the geometry
|
|
@param[out] newGeometries list of new geometries that have been created with the split
|
|
@param topological true if topological editing is enabled
|
|
@param[out] topologyTestPoints points that need to be tested for topological completeness in the dataset
|
|
@return 0 in case of success, 1 if geometry has not been split, error else*/
|
|
int splitGeometry( const QList<QgsPoint>& splitLine,
|
|
QList<QgsGeometry*>&newGeometries /Out/,
|
|
bool topological,
|
|
QList<QgsPoint> &topologyTestPoints /Out/);
|
|
|
|
/**Replaces a part of this geometry with another line
|
|
@return 0 in case of success */
|
|
int reshapeGeometry( const QList<QgsPoint>& reshapeWithLine );
|
|
|
|
/**Changes this geometry such that it does not intersect the other geometry
|
|
@param other geometry that should not be intersect
|
|
@return 0 in case of success*/
|
|
int makeDifference( QgsGeometry* other );
|
|
|
|
/**Returns the bounding box of this feature*/
|
|
QgsRectangle boundingBox();
|
|
|
|
/** Test for intersection with a rectangle (uses GEOS) */
|
|
bool intersects( const QgsRectangle& r ) const;
|
|
|
|
/** Test for intersection with a geometry (uses GEOS) */
|
|
bool intersects( const QgsGeometry* geometry ) const;
|
|
|
|
/** Test for containment of a point (uses GEOS) */
|
|
bool contains( const QgsPoint* p ) const;
|
|
|
|
/** Test for if geometry is contained in another (uses GEOS) */
|
|
bool contains( const QgsGeometry* geometry ) const;
|
|
|
|
/** Test for if geometry is disjoint of another (uses GEOS) */
|
|
bool disjoint( const QgsGeometry* geometry ) const;
|
|
|
|
/** Test for if geometry equals another (uses GEOS) */
|
|
bool equals( const QgsGeometry* geometry ) const;
|
|
|
|
/** Test for if geometry touch another (uses GEOS) */
|
|
bool touches( const QgsGeometry* geometry ) const;
|
|
|
|
/** Test for if geometry overlaps another (uses GEOS) */
|
|
bool overlaps( const QgsGeometry* geometry ) const;
|
|
|
|
/** Test for if geometry is within another (uses GEOS) */
|
|
bool within( const QgsGeometry* geometry ) const;
|
|
|
|
/** Test for if geometry crosses another (uses GEOS) */
|
|
bool crosses( const QgsGeometry* geometry ) const;
|
|
|
|
/** Returns a buffer region around this geometry having the given width and with a specified number
|
|
of segments used to approximate curves */
|
|
QgsGeometry* buffer( double distance, int segments ) /Factory/;
|
|
|
|
/** Returns a buffer region around the geometry, with additional style options.
|
|
* @param distance buffer distance
|
|
* @param segments For round joins, number of segments to approximate quarter-circle
|
|
* @param endCapStyle Round (1) / Flat (2) / Square (3) end cap style
|
|
* @param joinStyle Round (1) / Mitre (2) / Bevel (3) join style
|
|
* @param mitreLimit Limit on the mitre ratio used for very sharp corners
|
|
* @note added in 2.4
|
|
* @note needs GEOS >= 3.3 - otherwise always returns 0
|
|
*/
|
|
QgsGeometry* buffer( double distance, int segments, int endCapStyle, int joinStyle, double mitreLimit ) /Factory/;
|
|
|
|
/** Returns an offset line at a given distance and side from an input line.
|
|
* See buffer() method for details on parameters.
|
|
* @note added in 2.4
|
|
* @note needs GEOS >= 3.3 - otherwise always returns 0
|
|
*/
|
|
QgsGeometry* offsetCurve( double distance, int segments, int joinStyle, double mitreLimit ) /Factory/;
|
|
|
|
/** Returns a simplified version of this geometry using a specified tolerance value */
|
|
QgsGeometry* simplify( double tolerance ) /Factory/;
|
|
|
|
/** Returns the center of mass of a geometry
|
|
* @note for line based geometries, the center point of the line is returned,
|
|
* and for point based geometries, the point itself is returned */
|
|
QgsGeometry* centroid() /Factory/;
|
|
|
|
/** Returns a point within a geometry */
|
|
QgsGeometry* pointOnSurface() /Factory/;
|
|
|
|
/** Returns the smallest convex polygon that contains all the points in the geometry. */
|
|
QgsGeometry* convexHull() /Factory/;
|
|
|
|
/* Return interpolated point on line at distance */
|
|
QgsGeometry* interpolate( double distance ) /Factory/;
|
|
|
|
/** Returns a geometry representing the points shared by this geometry and other. */
|
|
QgsGeometry* intersection( QgsGeometry* geometry ) /Factory/;
|
|
|
|
/** Returns a geometry representing all the points in this geometry and other (a
|
|
* union geometry operation).
|
|
* @note this operation is not called union since its a reserved word in C++.*/
|
|
QgsGeometry* combine( QgsGeometry* geometry ) /Factory/;
|
|
|
|
/** Returns a geometry representing the points making up this geometry that do not make up other. */
|
|
QgsGeometry* difference( QgsGeometry* geometry ) /Factory/;
|
|
|
|
/** Returns a Geometry representing the points making up this Geometry that do not make up other. */
|
|
QgsGeometry* symDifference( QgsGeometry* geometry ) /Factory/;
|
|
|
|
/** Exports the geometry to WKT
|
|
* @note precision parameter added in 2.4
|
|
* @return true in case of success and false else
|
|
*/
|
|
QString exportToWkt( const int &precision = 17 ) const;
|
|
|
|
/** Exports the geometry to GeoJSON
|
|
* @return a QString representing the geometry as GeoJSON
|
|
* @note precision parameter added in 2.4
|
|
*/
|
|
QString exportToGeoJSON( const int &precision = 17 ) const;
|
|
|
|
/** try to convert the geometry to the requested type
|
|
* @param destType the geometry type to be converted to
|
|
* @param destMultipart determines if the output geometry will be multipart or not
|
|
* @return the converted geometry or NULL pointer if the conversion fails.
|
|
* @note added in 2.2
|
|
*/
|
|
QgsGeometry* convertToType( QGis::GeometryType destType, bool destMultipart = false ) /Factory/;
|
|
|
|
/* Accessor functions for getting geometry data */
|
|
|
|
/** return contents of the geometry as a point
|
|
if wkbType is WKBPoint, otherwise returns [0,0] */
|
|
QgsPoint asPoint() const;
|
|
|
|
/** return contents of the geometry as a polyline
|
|
if wkbType is WKBLineString, otherwise an empty list */
|
|
QgsPolyline asPolyline() const;
|
|
|
|
/** return contents of the geometry as a polygon
|
|
if wkbType is WKBPolygon, otherwise an empty list */
|
|
QgsPolygon asPolygon() const;
|
|
|
|
/** return contents of the geometry as a multi point
|
|
if wkbType is WKBMultiPoint, otherwise an empty list */
|
|
QgsMultiPoint asMultiPoint() const;
|
|
|
|
/** return contents of the geometry as a multi linestring
|
|
if wkbType is WKBMultiLineString, otherwise an empty list */
|
|
QgsMultiPolyline asMultiPolyline() const;
|
|
|
|
/** return contents of the geometry as a multi polygon
|
|
if wkbType is WKBMultiPolygon, otherwise an empty list */
|
|
QgsMultiPolygon asMultiPolygon() const;
|
|
|
|
/** return contents of the geometry as a list of geometries */
|
|
QList<QgsGeometry*> asGeometryCollection() const /Factory/;
|
|
|
|
/**Return contents of the geometry as a QPointF if wkbType is WKBPoint,
|
|
* otherwise returns a null QPointF.
|
|
* @note added in QGIS 2.7
|
|
*/
|
|
QPointF asQPointF() const;
|
|
|
|
/**Return contents of the geometry as a QPolygonF. If geometry is a linestring,
|
|
* then the result will be an open QPolygonF. If the geometry is a polygon,
|
|
* then the result will be a closed QPolygonF of the geometry's exterior ring.
|
|
* @note added in QGIS 2.7
|
|
*/
|
|
QPolygonF asQPolygonF() const;
|
|
|
|
/** delete a ring in polygon or multipolygon.
|
|
Ring 0 is outer ring and can't be deleted.
|
|
@return true on success */
|
|
bool deleteRing( int ringNum, int partNum = 0 );
|
|
|
|
/** delete part identified by the part number
|
|
@return true on success */
|
|
bool deletePart( int partNum );
|
|
|
|
/**Converts single type geometry into multitype geometry
|
|
e.g. a polygon into a multipolygon geometry with one polygon
|
|
@return true in case of success and false else*/
|
|
bool convertToMultiType();
|
|
|
|
/** Modifies geometry to avoid intersections with the layers specified in project properties
|
|
* @return 0 in case of success,
|
|
* 1 if geometry is not of polygon type,
|
|
* 2 if avoid intersection would change the geometry type,
|
|
* 3 other error during intersection removal
|
|
* @param ignoreFeatures possibility to give a list of features where intersections should be ignored (not available in python bindings)
|
|
*/
|
|
int avoidIntersections();
|
|
|
|
class Error
|
|
{
|
|
public:
|
|
Error();
|
|
Error( QString m );
|
|
Error( QString m, QgsPoint p );
|
|
|
|
QString what();
|
|
QgsPoint where();
|
|
bool hasWhere();
|
|
};
|
|
|
|
/** Validate geometry and produce a list of geometry errors */
|
|
void validateGeometry( QList<QgsGeometry::Error> &errors /Out/ );
|
|
|
|
/** compute the unary union on a list of geometries. May be faster than an iterative union on a set of geometries.
|
|
@param geometryList a list of QgsGeometry* as input
|
|
@returns the new computed QgsGeometry, or null
|
|
*/
|
|
static QgsGeometry *unaryUnion( const QList<QgsGeometry*>& geometryList ) /Factory/;
|
|
}; // class QgsGeometry
|
|
|