24QString QgsPixelCentroidsFromPolygonsAlgorithm::name()
const
26 return QStringLiteral(
"generatepointspixelcentroidsinsidepolygons" );
29QString QgsPixelCentroidsFromPolygonsAlgorithm::displayName()
const
31 return QObject::tr(
"Generate points (pixel centroids) inside polygons" );
34QStringList QgsPixelCentroidsFromPolygonsAlgorithm::tags()
const
36 return QObject::tr(
"raster,polygon,centroid,pixel,create" ).split(
',' );
39QString QgsPixelCentroidsFromPolygonsAlgorithm::group()
const
41 return QObject::tr(
"Vector creation" );
44QString QgsPixelCentroidsFromPolygonsAlgorithm::groupId()
const
46 return QStringLiteral(
"vectorcreation" );
49QString QgsPixelCentroidsFromPolygonsAlgorithm::shortHelpString()
const
51 return QObject::tr(
"Generates pixel centroids for the raster area falling inside polygons. Used to generate points "
52 "for further raster sampling." );
60QgsPixelCentroidsFromPolygonsAlgorithm *QgsPixelCentroidsFromPolygonsAlgorithm::createInstance()
const
62 return new QgsPixelCentroidsFromPolygonsAlgorithm();
65void QgsPixelCentroidsFromPolygonsAlgorithm::initAlgorithm(
const QVariantMap & )
75 QgsRasterLayer *rasterLayer = parameterAsRasterLayer( parameters, QStringLiteral(
"INPUT_RASTER" ), context );
80 std::unique_ptr< QgsProcessingFeatureSource > source( parameterAsSource( parameters, QStringLiteral(
"INPUT_VECTOR" ), context ) );
85 fields.
append(
QgsField( QStringLiteral(
"id" ), QMetaType::Type::LongLong ) );
86 fields.
append(
QgsField( QStringLiteral(
"poly_id" ), QMetaType::Type::Int ) );
87 fields.
append(
QgsField( QStringLiteral(
"point_id" ), QMetaType::Type::Int ) );
94 const double step = source->featureCount() ? 100.0 / source->featureCount() : 1;
122 const double xMin = bbox.
xMinimum();
123 const double xMax = bbox.
xMaximum();
124 const double yMin = bbox.
yMinimum();
125 const double yMax = bbox.
yMaximum();
128 int startRow, startColumn;
129 int endRow, endColumn;
130 QgsRasterAnalysisUtils::mapToPixel( xMin, yMax, extent, xPixel, yPixel, startRow, startColumn );
131 QgsRasterAnalysisUtils::mapToPixel( xMax, yMin, extent, xPixel, yPixel, endRow, endColumn );
133 std::unique_ptr< QgsGeos > engine = std::make_unique< QgsGeos >( f.
geometry().
constGet() );
134 engine->prepareGeometry();
136 for (
int row = startRow; row <= endRow; row++ )
138 for (
int col = startColumn; col <= endColumn; col++ )
145 QgsRasterAnalysisUtils::pixelToMap( row, col, extent, xPixel, yPixel, x, y );
146 if ( engine->contains( x, y ) )
148 feature.
setGeometry( std::make_unique< QgsPoint >( x, y ) );
166 outputs.insert( QStringLiteral(
"OUTPUT" ), dest );
@ VectorPoint
Vector point layers.
@ VectorPolygon
Vector polygon layers.
@ RegeneratesPrimaryKey
Algorithm always drops any existing primary keys or FID values and regenerates them in outputs.
QFlags< ProcessingAlgorithmDocumentationFlag > ProcessingAlgorithmDocumentationFlags
Flags describing algorithm behavior for documentation purposes.
Wrapper for iterator of features from vector data provider or vector layer.
bool nextFeature(QgsFeature &f)
Fetch next feature and stores in f, returns true on success.
This class wraps a request for features to a vector layer (or directly its vector data provider).
QgsFeatureRequest & setSubsetOfAttributes(const QgsAttributeList &attrs)
Set a subset of attributes that will be fetched.
@ FastInsert
Use faster inserts, at the cost of updating the passed features to reflect changes made at the provid...
@ RegeneratePrimaryKey
This flag indicates, that a primary key field cannot be guaranteed to be unique and the sink should i...
The feature class encapsulates a single feature including its unique ID, geometry and a list of field...
void setAttributes(const QgsAttributes &attrs)
Sets the feature's attributes.
void setFields(const QgsFields &fields, bool initAttributes=false)
Assigns a field map with the feature to allow attribute access by attribute name.
bool hasGeometry() const
Returns true if the feature has an associated geometry.
void setGeometry(const QgsGeometry &geometry)
Set the feature's geometry.
bool isCanceled() const
Tells whether the operation has been canceled already.
void setProgress(double progress)
Sets the current progress for the feedback object.
Encapsulate a field in an attribute table or data source.
Container of fields for a vector layer.
bool append(const QgsField &field, Qgis::FieldOrigin origin=Qgis::FieldOrigin::Provider, int originIndex=-1)
Appends a field.
const QgsAbstractGeometry * constGet() const
Returns a non-modifiable (const) reference to the underlying abstract geometry primitive.
QgsRectangle boundingBox() const
Returns the bounding box of the geometry.
virtual QgsRectangle extent() const
Returns the extent of the layer.
QgsCoordinateReferenceSystem crs
Contains information about the context in which a processing algorithm is executed.
QgsCoordinateTransformContext transformContext() const
Returns the coordinate transform context.
Custom exception class for processing related exceptions.
Base class for providing feedback from a processing algorithm.
A feature sink output for processing algorithms.
An input feature source (such as vector layers) parameter for processing algorithms.
A raster layer parameter for processing algorithms.
Represents a raster layer.
double rasterUnitsPerPixelX() const
Returns the number of raster units per each raster pixel in X axis.
double rasterUnitsPerPixelY() const
Returns the number of raster units per each raster pixel in Y axis.
A rectangle specified with double values.
double xMinimum() const
Returns the x minimum value (left side of rectangle).
double yMinimum() const
Returns the y minimum value (bottom side of rectangle).
double xMaximum() const
Returns the x maximum value (right side of rectangle).
double yMaximum() const
Returns the y maximum value (top side of rectangle).