20#include <QApplication>
24#include "moc_qgspointcloudlayerexporter.cpp"
31#include <pdal/StageFactory.hpp>
32#include <pdal/io/BufferReader.hpp>
33#include <pdal/Dimension.hpp>
41 return QStringLiteral(
"GPKG" );
43 return QStringLiteral(
"DXF" );
45 return QStringLiteral(
"ESRI Shapefile" );
47 return QStringLiteral(
"CSV" );
56 : mLayerAttributeCollection( layer->attributes() )
57 , mIndex( layer->dataProvider()->index()->clone().release() )
64 mPointRecordFormat = 3;
88 mFilterGeometryEngine.reset(
new QgsGeos( geometry ) );
89 mFilterGeometryEngine->prepareGeometry();
98 QVector< QgsGeometry > allGeometries;
101 if ( selectedFeaturesOnly )
113 allGeometries.append( f.
geometry() );
123 QgsDebugError( QStringLiteral(
"Error transforming union of filter layer: %1" ).arg( cse.
what() ) );
124 QgsDebugError( QStringLiteral(
"FilterGeometry will be ignored." ) );
132 mRequestedAttributes.clear();
134 const QVector<QgsPointCloudAttribute> allAttributes = mLayerAttributeCollection.
attributes();
138 if ( attribute.name().compare( QLatin1String(
"X" ), Qt::CaseInsensitive ) &&
139 attribute.name().compare( QLatin1String(
"Y" ), Qt::CaseInsensitive ) &&
140 attribute.name().compare( QLatin1String(
"Z" ), Qt::CaseInsensitive ) &&
141 attributeList.contains( attribute.name() ) &&
142 ! mRequestedAttributes.contains( attribute.name() ) )
144 mRequestedAttributes.append( attribute.name() );
151 QStringList allAttributeNames;
152 const QVector<QgsPointCloudAttribute> allAttributes = mLayerAttributeCollection.
attributes();
155 allAttributeNames.append( attribute.name() );
162 const QVector<QgsPointCloudAttribute> allAttributes = mLayerAttributeCollection.
attributes();
167 if ( attribute.name().compare( QLatin1String(
"X" ), Qt::CaseInsensitive ) ||
168 attribute.name().compare( QLatin1String(
"Y" ), Qt::CaseInsensitive ) ||
169 attribute.name().compare( QLatin1String(
"Z" ), Qt::CaseInsensitive ) ||
170 mRequestedAttributes.contains( attribute.name(), Qt::CaseInsensitive ) )
172 requestAttributes.
push_back( attribute );
175 return requestAttributes;
178QgsFields QgsPointCloudLayerExporter::outputFields()
185 if ( mRequestedAttributes.contains( attribute.name(), Qt::CaseInsensitive ) )
186 fields.
append(
QgsField( attribute.name(), attribute.variantType(), attribute.displayType() ) );
194 mMemoryLayer =
nullptr;
198 if ( QApplication::instance()->thread() != QThread::currentThread() )
199 QgsDebugMsgLevel( QStringLiteral(
"prepareExport() should better be called from the main thread!" ), 2 );
216 QgsDebugError( QStringLiteral(
"Error transforming extent: %1" ).arg( cse.
what() ) );
220 QStringList layerCreationOptions;
229 ExporterMemory exp(
this );
241 ExporterPdal exp(
this );
244 catch ( std::runtime_error &e )
246 setLastError( QString::fromLatin1( e.what() ) );
247 QgsDebugError( QStringLiteral(
"PDAL has thrown an exception: {}" ).arg( e.what() ) );
254 layerCreationOptions << QStringLiteral(
"GEOMETRY=AS_XYZ" )
255 << QStringLiteral(
"SEPARATOR=COMMA" );
272 ExporterVector exp(
this );
286 mMemoryLayer =
nullptr;
292 const QFileInfo fileInfo( mFilename );
293 return new QgsPointCloudLayer( mFilename, fileInfo.completeBaseName(), QStringLiteral(
"pdal" ) );
298 QString uri( mFilename );
299 uri +=
"|layername=" + mName;
307 const QFileInfo fileInfo( mFilename );
308 return new QgsVectorLayer( mFilename, fileInfo.completeBaseName(), QStringLiteral(
"ogr" ) );
318void QgsPointCloudLayerExporter::ExporterBase::run()
320 QgsRectangle geometryFilterRectangle( -std::numeric_limits<double>::infinity(),
321 -std::numeric_limits<double>::infinity(),
322 std::numeric_limits<double>::infinity(),
323 std::numeric_limits<double>::infinity(),
325 if ( mParent->mFilterGeometryEngine )
332 QVector<IndexedPointCloudNode> nodes;
333 qint64 pointCount = 0;
334 QQueue<IndexedPointCloudNode> queue;
335 queue.push_back( mParent->mIndex->root() );
336 while ( !queue.empty() )
340 const QgsRectangle nodeExtent = mParent->mIndex->nodeMapExtent( node );
341 if ( mParent->mExtent.intersects( nodeExtent ) &&
342 mParent->mZRange.overlaps( mParent->mIndex->nodeZRange( node ) ) &&
343 geometryFilterRectangle.intersects( nodeExtent ) )
345 pointCount += mParent->mIndex->nodePointCount( node );
346 nodes.push_back( node );
350 queue.push_back( child );
354 const qint64 pointsToExport = mParent->mPointsLimit > 0 ? std::min( mParent->mPointsLimit, pointCount ) : pointCount;
356 request.
setAttributes( mParent->requestedAttributeCollection() );
357 std::unique_ptr<QgsPointCloudBlock> block =
nullptr;
358 qint64 pointsExported = 0;
361 block = mParent->mIndex->nodeData( node, request );
363 const char *ptr = block->data();
364 int count = block->pointCount();
368 int xOffset = 0, yOffset = 0, zOffset = 0;
372 for (
int i = 0; i < count; ++i )
375 if ( mParent->mFeedback &&
378 mParent->mFeedback->setProgress( 100 *
static_cast< float >( pointsExported ) / pointsToExport );
379 if ( mParent->mFeedback->isCanceled() )
381 mParent->setLastError( QObject::tr(
"Canceled by user" ) );
386 if ( pointsExported >= pointsToExport )
396 if ( ! mParent->mZRange.contains( z ) ||
397 ! mParent->mExtent.contains( x, y ) ||
398 ( mParent->mFilterGeometryEngine && ! mParent->mFilterGeometryEngine->contains( x, y ) ) )
405 mParent->mTransform->transformInPlace( x, y, z );
407 handlePoint( x, y, z, attributeMap, pointsExported );
412 QgsDebugError( QStringLiteral(
"Error transforming point: %1" ).arg( cse.
what() ) );
429QgsPointCloudLayerExporter::ExporterMemory::~ExporterMemory()
431 mParent->mMemoryLayer->moveToThread( QApplication::instance()->thread() );
434void QgsPointCloudLayerExporter::ExporterMemory::handlePoint(
double x,
double y,
double z,
const QVariantMap &map,
const qint64 pointNumber )
436 Q_UNUSED( pointNumber )
441 for (
const QString &attribute : std::as_const( mParent->mRequestedAttributes ) )
443 const double val = map[ attribute ].toDouble();
444 featureAttributes.append( val );
447 mFeatures.append( feature );
450void QgsPointCloudLayerExporter::ExporterMemory::handleNode()
452 QgsVectorLayer *vl = qgis::down_cast<QgsVectorLayer *>( mParent->mMemoryLayer );
463void QgsPointCloudLayerExporter::ExporterMemory::handleAll()
477QgsPointCloudLayerExporter::ExporterVector::~ExporterVector()
479 delete mParent->mVectorSink;
480 mParent->mVectorSink =
nullptr;
483void QgsPointCloudLayerExporter::ExporterVector::handlePoint(
double x,
double y,
double z,
const QVariantMap &map,
const qint64 pointNumber )
485 Q_UNUSED( pointNumber )
490 for (
const QString &attribute : std::as_const( mParent->mRequestedAttributes ) )
492 const double val = map[ attribute ].toDouble();
493 featureAttributes.append( val );
496 mFeatures.append( feature );
499void QgsPointCloudLayerExporter::ExporterVector::handleNode()
501 if ( ! mParent->mVectorSink->addFeatures( mFeatures ) )
503 mParent->setLastError( mParent->mVectorSink->lastError() );
508void QgsPointCloudLayerExporter::ExporterVector::handleAll()
520 : mPointFormat( exp->mPointRecordFormat )
524 mOptions.add(
"filename", mParent->mFilename.toStdString() );
525 mOptions.add(
"a_srs", mParent->mTargetCrs.toWkt().toStdString() );
526 mOptions.add(
"minor_version", QStringLiteral(
"4" ).toStdString() );
527 mOptions.add(
"format", QString::number( mPointFormat ).toStdString() );
528 if ( mParent->mTransform->isShortCircuited() )
530 mOptions.add(
"offset_x", QString::number( mParent->mIndex->offset().x() ).toStdString() );
531 mOptions.add(
"offset_y", QString::number( mParent->mIndex->offset().y() ).toStdString() );
532 mOptions.add(
"offset_z", QString::number( mParent->mIndex->offset().z() ).toStdString() );
533 mOptions.add(
"scale_x", QString::number( mParent->mIndex->scale().x() ).toStdString() );
534 mOptions.add(
"scale_y", QString::number( mParent->mIndex->scale().y() ).toStdString() );
535 mOptions.add(
"scale_z", QString::number( mParent->mIndex->scale().z() ).toStdString() );
538 mTable.layout()->registerDim( pdal::Dimension::Id::X );
539 mTable.layout()->registerDim( pdal::Dimension::Id::Y );
540 mTable.layout()->registerDim( pdal::Dimension::Id::Z );
542 mTable.layout()->registerDim( pdal::Dimension::Id::Classification );
543 mTable.layout()->registerDim( pdal::Dimension::Id::Intensity );
544 mTable.layout()->registerDim( pdal::Dimension::Id::ReturnNumber );
545 mTable.layout()->registerDim( pdal::Dimension::Id::NumberOfReturns );
546 mTable.layout()->registerDim( pdal::Dimension::Id::ScanDirectionFlag );
547 mTable.layout()->registerDim( pdal::Dimension::Id::EdgeOfFlightLine );
548 mTable.layout()->registerDim( pdal::Dimension::Id::ScanAngleRank );
549 mTable.layout()->registerDim( pdal::Dimension::Id::UserData );
550 mTable.layout()->registerDim( pdal::Dimension::Id::PointSourceId );
552 if ( mPointFormat == 6 || mPointFormat == 7 || mPointFormat == 8 || mPointFormat == 9 || mPointFormat == 10 )
554 mTable.layout()->registerDim( pdal::Dimension::Id::ScanChannel );
555 mTable.layout()->registerDim( pdal::Dimension::Id::ClassFlags );
558 if ( mPointFormat != 0 && mPointFormat != 2 )
560 mTable.layout()->registerDim( pdal::Dimension::Id::GpsTime );
563 if ( mPointFormat == 2 || mPointFormat == 3 || mPointFormat == 5 || mPointFormat == 7 || mPointFormat == 8 || mPointFormat == 10 )
565 mTable.layout()->registerDim( pdal::Dimension::Id::Red );
566 mTable.layout()->registerDim( pdal::Dimension::Id::Green );
567 mTable.layout()->registerDim( pdal::Dimension::Id::Blue );
570 if ( mPointFormat == 8 || mPointFormat == 10 )
572 mTable.layout()->registerDim( pdal::Dimension::Id::Infrared );
575 mView.reset(
new pdal::PointView( mTable ) );
578void QgsPointCloudLayerExporter::ExporterPdal::handlePoint(
double x,
double y,
double z,
const QVariantMap &map,
const qint64 pointNumber )
580 mView->setField( pdal::Dimension::Id::X, pointNumber, x );
581 mView->setField( pdal::Dimension::Id::Y, pointNumber, y );
582 mView->setField( pdal::Dimension::Id::Z, pointNumber, z );
585 mView->setField( pdal::Dimension::Id::Classification, pointNumber, map[ QStringLiteral(
"Classification" ) ].toInt() );
586 mView->setField( pdal::Dimension::Id::Intensity, pointNumber, map[ QStringLiteral(
"Intensity" ) ].toInt() );
587 mView->setField( pdal::Dimension::Id::ReturnNumber, pointNumber, map[ QStringLiteral(
"ReturnNumber" ) ].toInt() );
588 mView->setField( pdal::Dimension::Id::NumberOfReturns, pointNumber, map[ QStringLiteral(
"NumberOfReturns" ) ].toInt() );
589 mView->setField( pdal::Dimension::Id::ScanDirectionFlag, pointNumber, map[ QStringLiteral(
"ScanDirectionFlag" ) ].toInt() );
590 mView->setField( pdal::Dimension::Id::EdgeOfFlightLine, pointNumber, map[ QStringLiteral(
"EdgeOfFlightLine" ) ].toInt() );
591 mView->setField( pdal::Dimension::Id::ScanAngleRank, pointNumber, map[ QStringLiteral(
"ScanAngleRank" ) ].toFloat() );
592 mView->setField( pdal::Dimension::Id::UserData, pointNumber, map[ QStringLiteral(
"UserData" ) ].toInt() );
593 mView->setField( pdal::Dimension::Id::PointSourceId, pointNumber, map[ QStringLiteral(
"PointSourceId" ) ].toInt() );
595 if ( mPointFormat == 6 || mPointFormat == 7 || mPointFormat == 8 || mPointFormat == 9 || mPointFormat == 10 )
597 mView->setField( pdal::Dimension::Id::ScanChannel, pointNumber, map[ QStringLiteral(
"ScannerChannel" ) ].toInt() );
598 const int classificationFlags = ( map[ QStringLiteral(
"Synthetic" ) ].toInt() & 0x01 ) << 0 |
599 ( map[ QStringLiteral(
"KeyPoint" ) ].toInt() & 0x01 ) << 1 |
600 ( map[ QStringLiteral(
"Withheld" ) ].toInt() & 0x01 ) << 2 |
601 ( map[ QStringLiteral(
"Overlap" ) ].toInt() & 0x01 ) << 3;
602 mView->setField( pdal::Dimension::Id::ClassFlags, pointNumber, classificationFlags );
605 if ( mPointFormat != 0 && mPointFormat != 2 )
607 mView->setField( pdal::Dimension::Id::GpsTime, pointNumber, map[ QStringLiteral(
"GpsTime" ) ].toDouble() );
610 if ( mPointFormat == 2 || mPointFormat == 3 || mPointFormat == 5 || mPointFormat == 7 || mPointFormat == 8 || mPointFormat == 10 )
612 mView->setField( pdal::Dimension::Id::Red, pointNumber, map[ QStringLiteral(
"Red" ) ].toInt() );
613 mView->setField( pdal::Dimension::Id::Green, pointNumber, map[ QStringLiteral(
"Green" ) ].toInt() );
614 mView->setField( pdal::Dimension::Id::Blue, pointNumber, map[ QStringLiteral(
"Blue" ) ].toInt() );
617 if ( mPointFormat == 8 || mPointFormat == 10 )
619 mView->setField( pdal::Dimension::Id::Infrared, pointNumber, map[ QStringLiteral(
"Infrared" ) ].toInt() );
623void QgsPointCloudLayerExporter::ExporterPdal::handleNode()
628void QgsPointCloudLayerExporter::ExporterPdal::handleAll()
630 pdal::BufferReader reader;
631 reader.addView( mView );
633 pdal::StageFactory factory;
635 pdal::Stage *writer = factory.createStage(
"writers.las" );
637 writer->setInput( reader );
638 writer->setOptions( mOptions );
639 writer->prepare( mTable );
640 writer->execute( mTable );
657 mOwnedFeedback->cancel();
Represents a indexed point cloud node in octree.
@ NoSymbology
Export only data.
@ Reverse
Reverse/inverse transform (from destination to source)
Abstract base class for all geometries.
virtual QgsRectangle boundingBox() const
Returns the minimal bounding box for the geometry.
This class represents a coordinate reference system (CRS).
Contains information about the context in which a coordinate transform is executed.
Custom exception class for Coordinate Reference System related exceptions.
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 & setNoAttributes()
Set that no attributes will be fetched.
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.
bool hasGeometry() const
Returns true if the feature has an associated geometry.
void setGeometry(const QgsGeometry &geometry)
Set the feature's geometry.
Base class for feedback objects to be used for cancellation of something running in a worker thread.
void progressChanged(double progress)
Emitted when the feedback object reports a progress change.
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.
A geometry is the spatial representation of a feature.
Qgis::GeometryOperationResult transform(const QgsCoordinateTransform &ct, Qgis::TransformDirection direction=Qgis::TransformDirection::Forward, bool transformZ=false)
Transforms this geometry as described by the coordinate transform ct.
const QgsAbstractGeometry * constGet() const
Returns a non-modifiable (const) reference to the underlying abstract geometry primitive.
static QgsGeometry unaryUnion(const QVector< QgsGeometry > &geometries, const QgsGeometryParameters ¶meters=QgsGeometryParameters())
Compute the unary union on a list of geometries.
Does vector analysis using the geos library and handles import, export, exception handling*.
Base class for all map layer types.
QgsCoordinateReferenceSystem crs
static QgsVectorLayer * createMemoryLayer(const QString &name, const QgsFields &fields, Qgis::WkbType geometryType=Qgis::WkbType::NoGeometry, const QgsCoordinateReferenceSystem &crs=QgsCoordinateReferenceSystem(), bool loadDefaultStyle=true) SIP_FACTORY
Creates a new memory layer using the specified parameters.
Collection of point cloud attributes.
void push_back(const QgsPointCloudAttribute &attribute)
Adds extra attribute.
int pointRecordSize() const
Returns total size of record.
const QgsPointCloudAttribute * find(const QString &attributeName, int &offset) const
Finds the attribute with the name.
QVector< QgsPointCloudAttribute > attributes() const
Returns all attributes.
Attribute for point cloud data pair of name and size in bytes.
DataType
Systems of unit measurement.
static void getPointXYZ(const char *ptr, int i, std::size_t pointRecordSize, int xOffset, QgsPointCloudAttribute::DataType xType, int yOffset, QgsPointCloudAttribute::DataType yType, int zOffset, QgsPointCloudAttribute::DataType zType, const QgsVector3D &indexScale, const QgsVector3D &indexOffset, double &x, double &y, double &z)
Retrieves the x, y, z values for the point at index i.
static QVariantMap getAttributeMap(const char *data, std::size_t recordOffset, const QgsPointCloudAttributeCollection &attributeCollection)
Retrieves all the attributes of a point.
DataType type() const
Returns the data type.
virtual QVariantMap originalMetadata() const
Returns a representation of the original metadata included in a point cloud dataset.
void cancel() override
Notifies the task that it should terminate.
QgsPointCloudLayerExporterTask(QgsPointCloudLayerExporter *exporter)
Constructor for QgsPointCloudLayerExporterTask.
void exportComplete()
Emitted when exporting the layer is successfully completed.
void finished(bool result) override
If the task is managed by a QgsTaskManager, this will be called after the task has finished (whether ...
bool run() override
Performs the task's operation.
Handles exporting point cloud layers to memory layers, OGR supported files and PDAL supported files.
void setFeedback(QgsFeedback *feedback)
Sets a QgsFeedback object to allow cancellation / progress reporting.
QgsMapLayer * takeExportedLayer()
Gets a pointer to the exported layer.
ExportFormat format() const
Returns the format for the exported file or layer.
void setAttributes(const QStringList &attributes)
Sets the list of point cloud attributes that will be exported.
ExportFormat
Supported export formats for point clouds.
@ Csv
Comma separated values.
@ Las
LAS/LAZ point cloud.
~QgsPointCloudLayerExporter()
void setAllAttributes()
Sets that all attributes will be exported.
bool setFormat(const ExportFormat format)
Sets the format for the exported file.
static QString getOgrDriverName(ExportFormat format)
Gets the OGR driver name for the specified format.
QStringList attributes() const
Gets the list of point cloud attributes that will be exported.
QgsPointCloudLayerExporter(QgsPointCloudLayer *layer)
Constructor for QgsPointCloudLayerExporter, associated with the specified layer.
void prepareExport()
Creates the QgsVectorLayer for exporting to a memory layer, if necessary.
static QList< ExportFormat > supportedFormats()
Gets a list of the supported export formats.
void setFilterGeometry(const QgsAbstractGeometry *geometry)
Sets a spatial filter for points to be exported based on geom in the point cloud's CRS.
void doExport()
Performs the actual exporting operation.
Represents a map layer supporting display of point clouds.
QgsPointCloudDataProvider * dataProvider() override
Returns the layer's data provider, it may be nullptr.
Point cloud data request.
void setAttributes(const QgsPointCloudAttributeCollection &attributes)
Set attributes filter in the request.
Point geometry type, with support for z-dimension and m-values.
A rectangle specified with double values.
bool isFinite() const
Returns true if the rectangle has finite boundaries.
Abstract base class for long running background tasks.
virtual void cancel()
Notifies the task that it should terminate.
void setProgress(double progress)
Sets the task's current progress.
Class for storage of 3D vectors similar to QVector3D, with the difference that it uses double precisi...
QString lastError() const override
Returns the most recent error encountered by the sink, e.g.
bool addFeatures(QgsFeatureList &flist, QgsFeatureSink::Flags flags=QgsFeatureSink::Flags()) override
Adds a list of features to the sink.
Options to pass to writeAsVectorFormat()
QString driverName
OGR driver to use.
QString layerName
Layer name. If let empty, it will be derived from the filename.
QStringList layerOptions
List of OGR layer creation options.
Qgis::FeatureSymbologyExport symbologyExport
Symbology to export.
QgsVectorFileWriter::ActionOnExistingFile actionOnExistingFile
Action on existing file.
QStringList datasourceOptions
List of OGR data source creation options.
QgsFeedback * feedback
Optional feedback object allowing cancellation of layer save.
static QStringList defaultLayerOptions(const QString &driverName)
Returns a list of the default layer options for a specified driver.
static QgsVectorFileWriter * create(const QString &fileName, const QgsFields &fields, Qgis::WkbType geometryType, const QgsCoordinateReferenceSystem &srs, const QgsCoordinateTransformContext &transformContext, const QgsVectorFileWriter::SaveVectorOptions &options, QgsFeatureSink::SinkFlags sinkFlags=QgsFeatureSink::SinkFlags(), QString *newFilename=nullptr, QString *newLayer=nullptr)
Create a new vector file writer.
static QStringList defaultDatasetOptions(const QString &driverName)
Returns a list of the default dataset options for a specified driver.
Represents a vector layer which manages a vector based data sets.
QgsFeatureIterator getFeatures(const QgsFeatureRequest &request=QgsFeatureRequest()) const FINAL
Queries the layer for features specified in request.
QgsFeatureIterator getSelectedFeatures(QgsFeatureRequest request=QgsFeatureRequest()) const
Returns an iterator of the selected features.
QgsVectorDataProvider * dataProvider() FINAL
Returns the layer's data provider, it may be nullptr.
#define BUILTIN_UNREACHABLE
#define QgsDebugMsgLevel(str, level)
#define QgsDebugError(str)
const QgsCoordinateReferenceSystem & crs