YXS
/
nmWTAI-Platform
8
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

实现插值结果实时预览

Browse Source
feature/Interpolation-20260803
lvjunjie 2 weeks ago
parent fde6316a37
commit d3bb18b24f
4 changed files with 283 additions and 54 deletions
Show Stats Download Patch File Download Diff File

8
Include/nmNum/nmSubWxs/nmWxPropertyInterpolationDlg.h
Unescape Escape View File

@ -18,8 +18,10 @@ class QPointF;
class QPushButton; class QPushButton;
class QTableWidget; class QTableWidget;
class QTableWidgetItem; class QTableWidgetItem;
class QTimer;
class nmGuiPlot; class nmGuiPlot;
class nmObjPoint; class nmObjPoint;
class nmWxPropertyInterpolationPreviewDlg;
class NM_SUB_WXS_EXPORT nmWxPropertyInterpolationDlg : public iDlgBase class NM_SUB_WXS_EXPORT nmWxPropertyInterpolationDlg : public iDlgBase
{ {
@ -49,6 +51,7 @@ private slots:
void onShowPointsToggled(bool checked); void onShowPointsToggled(bool checked);
void onShowLabelsToggled(bool checked); void onShowLabelsToggled(bool checked);
void onPreview(); void onPreview();
void refreshPreview();
private: private:
// 初始化界面及持久化数据组. // 初始化界面及持久化数据组.
@ -80,6 +83,9 @@ private:
void clearTableSelection(); void clearTableSelection();
void updatePointButtons(); void updatePointButtons();
void updateCurrentMarkerVisibility(); void updateCurrentMarkerVisibility();
void schedulePreviewRefresh();
void updatePreviewWindow(bool activateWindow);
void showPreviewError(const QString& message, bool activateWindow);
private: private:
QPointer<nmGuiPlot> m_pPlot; QPointer<nmGuiPlot> m_pPlot;
@ -119,6 +125,8 @@ private:
// 结果预览及应用. // 结果预览及应用.
QPushButton* m_pPreviewButton; QPushButton* m_pPreviewButton;
QPushButton* m_pApplyButton; QPushButton* m_pApplyButton;
QPointer<nmWxPropertyInterpolationPreviewDlg> m_pPreviewDialog;
QTimer* m_pPreviewRefreshTimer;
}; };
#endif // NMWXPROPERTYINTERPOLATIONDLG_H #endif // NMWXPROPERTYINTERPOLATIONDLG_H

13
Include/nmNum/nmSubWxs/nmWxPropertyInterpolationPreviewDlg.h
Unescape Escape View File

@ -29,6 +29,19 @@ public:
bool showMeasurementLabels, bool showMeasurementLabels,
QWidget* parent = 0); QWidget* parent = 0);
void updatePreview(const QString& dataSetName,
const QString& scalarTitle,
const QRectF& bounds,
int columnCount,
int rowCount,
const QVector<double>& interpolationValues,
const QVector<QPointF>& measurementPoints,
const QVector<double>& measurementValues,
const QVector<QPointF>& outlinePoints,
bool showMeasurementPoints,
bool showMeasurementLabels);
void showMessage(const QString& dataSetName, const QString& message);
private: private:
void initializeRenderer(const QString& scalarTitle, void initializeRenderer(const QString& scalarTitle,
const QRectF& bounds, const QRectF& bounds,

137
Src/nmNum/nmSubWxs/nmWxPropertyInterpolationDlg.cpp
Unescape Escape View File

@ -22,6 +22,7 @@
#include <QSize> #include <QSize>
#include <QTableWidget> #include <QTableWidget>
#include <QTableWidgetItem> #include <QTableWidgetItem>
#include <QTimer>
#include <QTransform> #include <QTransform>
#include <QVBoxLayout> #include <QVBoxLayout>
#include <QtCore/qmath.h> #include <QtCore/qmath.h>
@ -185,13 +186,22 @@ nmWxPropertyInterpolationDlg::nmWxPropertyInterpolationDlg(
m_pRangeSpin(NULL), m_pRangeSpin(NULL),
m_pModelCombo(NULL), m_pModelCombo(NULL),
m_pPreviewButton(NULL), m_pPreviewButton(NULL),
m_pApplyButton(NULL) m_pApplyButton(NULL),
m_pPreviewDialog(NULL),
m_pPreviewRefreshTimer(NULL)
{ {
setWindowTitle(interpolationText("Property Interpolation")); setWindowTitle(interpolationText("Property Interpolation"));
setWindowModality(Qt::NonModal); setWindowModality(Qt::NonModal);
setModal(false); setModal(false);
setMinimumSize(520, 620); setMinimumSize(520, 620);
resize(560, 680); resize(560, 680);
m_pPreviewRefreshTimer = new QTimer(this);
m_pPreviewRefreshTimer->setSingleShot(true);
m_pPreviewRefreshTimer->setInterval(200);
connect(m_pPreviewRefreshTimer, SIGNAL(timeout()),
this, SLOT(refreshPreview()));
initInterpolationUI(); initInterpolationUI();
initializeDataSets(); initializeDataSets();
@ -554,6 +564,7 @@ void nmWxPropertyInterpolationDlg::onDataSetChanged(int index)
syncDataSetsToManager(false); syncDataSetsToManager(false);
m_nCurrentDataSetIndex = index; m_nCurrentDataSetIndex = index;
loadDataSet(index); loadDataSet(index);
schedulePreviewRefresh();
} }
void nmWxPropertyInterpolationDlg::onAddDataSet() void nmWxPropertyInterpolationDlg::onAddDataSet()
@ -572,6 +583,7 @@ void nmWxPropertyInterpolationDlg::onAddDataSet()
loadDataSet(newIndex); loadDataSet(newIndex);
updateDataSetButtons(); updateDataSetButtons();
syncDataSetsToManager(true); syncDataSetsToManager(true);
schedulePreviewRefresh();
} }
void nmWxPropertyInterpolationDlg::onDeleteDataSet() void nmWxPropertyInterpolationDlg::onDeleteDataSet()
@ -597,6 +609,7 @@ void nmWxPropertyInterpolationDlg::onDeleteDataSet()
loadDataSet(newIndex); loadDataSet(newIndex);
updateDataSetButtons(); updateDataSetButtons();
syncDataSetsToManager(true); syncDataSetsToManager(true);
schedulePreviewRefresh();
} }
void nmWxPropertyInterpolationDlg::onDataSetNameChanged(const QString& name) void nmWxPropertyInterpolationDlg::onDataSetNameChanged(const QString& name)
@ -609,6 +622,7 @@ void nmWxPropertyInterpolationDlg::onDataSetNameChanged(const QString& name)
m_vecDataSets[m_nCurrentDataSetIndex].name = name; m_vecDataSets[m_nCurrentDataSetIndex].name = name;
m_pDataSetCombo->setItemText(m_nCurrentDataSetIndex, name); m_pDataSetCombo->setItemText(m_nCurrentDataSetIndex, name);
syncDataSetsToManager(true); syncDataSetsToManager(true);
schedulePreviewRefresh();
} }
void nmWxPropertyInterpolationDlg::onCurrentDataSetEdited() void nmWxPropertyInterpolationDlg::onCurrentDataSetEdited()
@ -619,6 +633,7 @@ void nmWxPropertyInterpolationDlg::onCurrentDataSetEdited()
saveCurrentDataSet(); saveCurrentDataSet();
syncDataSetsToManager(true); syncDataSetsToManager(true);
schedulePreviewRefresh();
} }
bool nmWxPropertyInterpolationDlg::eventFilter(QObject* watched, QEvent* event) bool nmWxPropertyInterpolationDlg::eventFilter(QObject* watched, QEvent* event)
@ -706,6 +721,7 @@ void nmWxPropertyInterpolationDlg::appendPoint(const QPointF& valuePoint)
m_pPointTable->scrollToBottom(); m_pPointTable->scrollToBottom();
updatePointButtons(); updatePointButtons();
syncDataSetsToManager(true); syncDataSetsToManager(true);
schedulePreviewRefresh();
} }
void nmWxPropertyInterpolationDlg::appendPointRow( void nmWxPropertyInterpolationDlg::appendPointRow(
@ -871,6 +887,7 @@ void nmWxPropertyInterpolationDlg::onPointItemChanged(QTableWidgetItem* item)
marker->update(); marker->update();
} }
syncDataSetsToManager(true); syncDataSetsToManager(true);
schedulePreviewRefresh();
} }
void nmWxPropertyInterpolationDlg::onDeleteSelectedPoint() void nmWxPropertyInterpolationDlg::onDeleteSelectedPoint()
@ -892,6 +909,7 @@ void nmWxPropertyInterpolationDlg::onDeleteSelectedPoint()
clearTableSelection(); clearTableSelection();
updatePointButtons(); updatePointButtons();
syncDataSetsToManager(true); syncDataSetsToManager(true);
schedulePreviewRefresh();
} }
void nmWxPropertyInterpolationDlg::onClearPoints() void nmWxPropertyInterpolationDlg::onClearPoints()
@ -907,6 +925,7 @@ void nmWxPropertyInterpolationDlg::onClearPoints()
clearTableSelection(); clearTableSelection();
updatePointButtons(); updatePointButtons();
syncDataSetsToManager(true); syncDataSetsToManager(true);
schedulePreviewRefresh();
} }
void nmWxPropertyInterpolationDlg::onShowPointsToggled(bool checked) void nmWxPropertyInterpolationDlg::onShowPointsToggled(bool checked)
@ -919,6 +938,7 @@ void nmWxPropertyInterpolationDlg::onShowPointsToggled(bool checked)
m_vecDataSets[m_nCurrentDataSetIndex].showPoints = checked; m_vecDataSets[m_nCurrentDataSetIndex].showPoints = checked;
updateCurrentMarkerVisibility(); updateCurrentMarkerVisibility();
syncDataSetsToManager(true); syncDataSetsToManager(true);
schedulePreviewRefresh();
} }
void nmWxPropertyInterpolationDlg::onShowLabelsToggled(bool checked) void nmWxPropertyInterpolationDlg::onShowLabelsToggled(bool checked)
@ -931,6 +951,7 @@ void nmWxPropertyInterpolationDlg::onShowLabelsToggled(bool checked)
m_vecDataSets[m_nCurrentDataSetIndex].showLabels = checked; m_vecDataSets[m_nCurrentDataSetIndex].showLabels = checked;
updateCurrentMarkerVisibility(); updateCurrentMarkerVisibility();
syncDataSetsToManager(true); syncDataSetsToManager(true);
schedulePreviewRefresh();
} }
void nmWxPropertyInterpolationDlg::clearTableSelection() void nmWxPropertyInterpolationDlg::clearTableSelection()
@ -950,10 +971,64 @@ void nmWxPropertyInterpolationDlg::updatePointButtons()
} }
void nmWxPropertyInterpolationDlg::onPreview() void nmWxPropertyInterpolationDlg::onPreview()
{
if(m_pPreviewRefreshTimer != NULL) {
m_pPreviewRefreshTimer->stop();
}
updatePreviewWindow(true);
}
void nmWxPropertyInterpolationDlg::schedulePreviewRefresh()
{
if(m_pPreviewDialog.isNull() || !m_pPreviewDialog->isVisible() ||
m_pPreviewRefreshTimer == NULL) {
return;
}
// 连续编辑参数时重置计时，仅对最终状态执行一次插值计算.
m_pPreviewRefreshTimer->start();
}
void nmWxPropertyInterpolationDlg::refreshPreview()
{
if(m_pPreviewDialog.isNull() || !m_pPreviewDialog->isVisible()) {
return;
}
updatePreviewWindow(false);
}
void nmWxPropertyInterpolationDlg::showPreviewError(
const QString& message,
bool activateWindow)
{
const QString errorText = message.isEmpty() ?
interpolationText("Preview Failed") : message;
QString dataSetName;
if(m_nCurrentDataSetIndex >= 0 &&
m_nCurrentDataSetIndex < m_vecDataSets.size()) {
dataSetName = m_vecDataSets[m_nCurrentDataSetIndex].name;
}
if(!m_pPreviewDialog.isNull()) {
m_pPreviewDialog->showMessage(dataSetName, errorText);
if(activateWindow) {
m_pPreviewDialog->show();
m_pPreviewDialog->raise();
m_pPreviewDialog->activateWindow();
}
}
else if(activateWindow) {
QMessageBox::warning(this, interpolationText("Preview Failed"), errorText);
}
}
void nmWxPropertyInterpolationDlg::updatePreviewWindow(bool activateWindow)
{ {
if(m_nCurrentDataSetIndex < 0 || if(m_nCurrentDataSetIndex < 0 ||
m_nCurrentDataSetIndex >= m_vecDataSets.size() || m_nCurrentDataSetIndex >= m_vecDataSets.size() ||
m_pDataManager.isNull()) { m_pDataManager.isNull()) {
showPreviewError(interpolationText("Preview Failed"), activateWindow);
return; return;
} }
@ -966,9 +1041,9 @@ void nmWxPropertyInterpolationDlg::onPreview()
QRectF previewBounds; QRectF previewBounds;
if(!createPreviewOutline(m_pDataManager->getOutlineData(), if(!createPreviewOutline(m_pDataManager->getOutlineData(),
outlinePoints, previewBounds)) { outlinePoints, previewBounds)) {
QMessageBox::warning(this, showPreviewError(
interpolationText("Preview Failed"), interpolationText("Please create a valid reservoir boundary before previewing."),
interpolationText("Please create a valid reservoir boundary before previewing.")); activateWindow);
return; return;
} }
@ -1016,27 +1091,47 @@ void nmWxPropertyInterpolationDlg::onPreview()
QApplication::restoreOverrideCursor(); QApplication::restoreOverrideCursor();
if(!calculationSucceeded) { if(!calculationSucceeded) {
QMessageBox::warning(this, showPreviewError(errorMessage, activateWindow);
interpolationText("Preview Failed"),
errorMessage);
return; return;
} }
// 预览窗口只负责裁剪和渲染，不修改数据组及正式求解参数. // 预览窗口只负责裁剪和渲染，不修改数据组及正式求解参数.
nmWxPropertyInterpolationPreviewDlg previewDialog( if(m_pPreviewDialog.isNull()) {
dataSet.name, m_pPreviewDialog = new nmWxPropertyInterpolationPreviewDlg(
m_pPropertyCombo->currentText(), dataSet.name,
previewBounds, m_pPropertyCombo->currentText(),
columnCount, previewBounds,
rowCount, columnCount,
interpolationValues, rowCount,
measurementPoints, interpolationValues,
measurementValues, measurementPoints,
outlinePoints, measurementValues,
dataSet.showPoints, outlinePoints,
dataSet.showLabels, dataSet.showPoints,
this); dataSet.showLabels,
previewDialog.exec(); this);
m_pPreviewDialog->setAttribute(Qt::WA_DeleteOnClose, true);
}
else {
m_pPreviewDialog->updatePreview(
dataSet.name,
m_pPropertyCombo->currentText(),
previewBounds,
columnCount,
rowCount,
interpolationValues,
measurementPoints,
measurementValues,
outlinePoints,
dataSet.showPoints,
dataSet.showLabels);
}
if(activateWindow) {
m_pPreviewDialog->show();
m_pPreviewDialog->raise();
m_pPreviewDialog->activateWindow();
}
} }
void nmWxPropertyInterpolationDlg::updateCurrentMarkerVisibility() void nmWxPropertyInterpolationDlg::updateCurrentMarkerVisibility()

179
Src/nmNum/nmSubWxs/nmWxPropertyInterpolationPreviewDlg.cpp
Unescape Escape View File

@ -12,10 +12,9 @@
#include <vtkCamera.h> #include <vtkCamera.h>
#include <vtkCellArray.h> #include <vtkCellArray.h>
#include <vtkDataArray.h> #include <vtkDataArray.h>
#include <vtkDelaunay2D.h>
#include <vtkDoubleArray.h> #include <vtkDoubleArray.h>
#include <vtkExtractPolyDataGeometry.h> #include <vtkExtractPolyDataGeometry.h>
#include <vtkImageData.h>
#include <vtkImageDataGeometryFilter.h>
#include <vtkImplicitSelectionLoop.h> #include <vtkImplicitSelectionLoop.h>
#include <vtkInteractorStyleImage.h> #include <vtkInteractorStyleImage.h>
#include <vtkLabeledDataMapper.h> #include <vtkLabeledDataMapper.h>
@ -55,7 +54,9 @@ vtkSmartPointer<vtkLookupTable> createPreviewLookupTable(double minimum,
return lookupTable; return lookupTable;
} }
void showPreviewFailure(vtkRenderer* renderer, vtkRenderWindow* renderWindow) void showPreviewFailure(vtkRenderer* renderer,
vtkRenderWindow* renderWindow,
const QString& message)
{ {
if(renderer == NULL || renderWindow == NULL) { if(renderer == NULL || renderWindow == NULL) {
return; return;
@ -63,7 +64,9 @@ void showPreviewFailure(vtkRenderer* renderer, vtkRenderWindow* renderWindow)
vtkSmartPointer<vtkTextActor> messageActor = vtkSmartPointer<vtkTextActor> messageActor =
vtkSmartPointer<vtkTextActor>::New(); vtkSmartPointer<vtkTextActor>::New();
messageActor->SetInput(previewText("Preview Failed").toUtf8().constData()); const QString displayMessage = message.isEmpty() ?
previewText("Preview Failed") : message;
messageActor->SetInput(displayMessage.toUtf8().constData());
messageActor->SetDisplayPosition(24, 24); messageActor->SetDisplayPosition(24, 24);
messageActor->GetTextProperty()->SetColor(0.75, 0.15, 0.12); messageActor->GetTextProperty()->SetColor(0.75, 0.15, 0.12);
messageActor->GetTextProperty()->SetFontSize(16); messageActor->GetTextProperty()->SetFontSize(16);
@ -88,8 +91,8 @@ nmWxPropertyInterpolationPreviewDlg::nmWxPropertyInterpolationPreviewDlg(
: iDlgBase(parent), : iDlgBase(parent),
m_pVtkWidget(NULL) m_pVtkWidget(NULL)
{ {
setWindowTitle(previewText("Interpolation Preview") + " - " + dataSetName); setWindowModality(Qt::NonModal);
setWindowModality(Qt::WindowModal); setModal(false);
setMinimumSize(900, 700); setMinimumSize(900, 700);
resize(1200, 850); resize(1200, 850);
@ -108,12 +111,57 @@ nmWxPropertyInterpolationPreviewDlg::nmWxPropertyInterpolationPreviewDlg(
mainLayout->addLayout(buttonLayout); mainLayout->addLayout(buttonLayout);
connect(closeButton, SIGNAL(clicked()), this, SLOT(close())); connect(closeButton, SIGNAL(clicked()), this, SLOT(close()));
updatePreview(dataSetName, scalarTitle, bounds, columnCount, rowCount,
interpolationValues, measurementPoints,
measurementValues, outlinePoints,
showMeasurementPoints, showMeasurementLabels);
}
void nmWxPropertyInterpolationPreviewDlg::updatePreview(
const QString& dataSetName,
const QString& scalarTitle,
const QRectF& bounds,
int columnCount,
int rowCount,
const QVector<double>& interpolationValues,
const QVector<QPointF>& measurementPoints,
const QVector<double>& measurementValues,
const QVector<QPointF>& outlinePoints,
bool showMeasurementPoints,
bool showMeasurementLabels)
{
setWindowTitle(previewText("Interpolation Preview") + " - " + dataSetName);
initializeRenderer(scalarTitle, bounds, columnCount, rowCount, initializeRenderer(scalarTitle, bounds, columnCount, rowCount,
interpolationValues, measurementPoints, interpolationValues, measurementPoints,
measurementValues, outlinePoints, measurementValues, outlinePoints,
showMeasurementPoints, showMeasurementLabels); showMeasurementPoints, showMeasurementLabels);
} }
void nmWxPropertyInterpolationPreviewDlg::showMessage(
const QString& dataSetName,
const QString& message)
{
if(m_pVtkWidget == NULL) {
return;
}
QString windowTitle = previewText("Interpolation Preview");
if(!dataSetName.isEmpty()) {
windowTitle += " - " + dataSetName;
}
setWindowTitle(windowTitle);
vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();
renderer->SetBackground(0.88, 0.88, 0.88);
vtkSmartPointer<vtkRenderWindow> renderWindow =
vtkSmartPointer<vtkRenderWindow>::New();
renderWindow->SetMultiSamples(0);
renderWindow->AddRenderer(renderer);
m_pVtkWidget->SetRenderWindow(renderWindow);
showPreviewFailure(renderer, renderWindow, message);
}
void nmWxPropertyInterpolationPreviewDlg::initializeRenderer( void nmWxPropertyInterpolationPreviewDlg::initializeRenderer(
const QString& scalarTitle, const QString& scalarTitle,
const QRectF& bounds, const QRectF& bounds,
@ -143,35 +191,73 @@ void nmWxPropertyInterpolationPreviewDlg::initializeRenderer(
if(columnCount < 2 || rowCount < 2 || if(columnCount < 2 || rowCount < 2 ||
bounds.width() <= 0.0 || bounds.height() <= 0.0 || bounds.width() <= 0.0 || bounds.height() <= 0.0 ||
interpolationValues.size() != columnCount * rowCount) { interpolationValues.size() != columnCount * rowCount) {
showPreviewFailure(renderer, renderWindow); showPreviewFailure(renderer, renderWindow, previewText("Preview Failed"));
return; return;
} }
// 标量按行优先顺序保存，与对话框生成规则待插值点的顺序一致. // 规则采样点使用Kriging结果，测点作为额外顶点直接使用已知值.
vtkSmartPointer<vtkPoints> surfacePoints =
vtkSmartPointer<vtkPoints>::New();
vtkSmartPointer<vtkDoubleArray> scalarArray = vtkSmartPointer<vtkDoubleArray> scalarArray =
vtkSmartPointer<vtkDoubleArray>::New(); vtkSmartPointer<vtkDoubleArray>::New();
scalarArray->SetName("InterpolatedValue"); scalarArray->SetName("InterpolatedValue");
scalarArray->SetNumberOfComponents(1); scalarArray->SetNumberOfComponents(1);
scalarArray->SetNumberOfTuples(interpolationValues.size());
for(int i = 0; i < interpolationValues.size(); ++i) { const double spacingX = bounds.width() / (columnCount - 1);
scalarArray->SetValue(i, interpolationValues[i]); const double spacingY = bounds.height() / (rowCount - 1);
for(int row = 0; row < rowCount; ++row) {
const double y = bounds.top() + spacingY * row;
for(int column = 0; column < columnCount; ++column) {
const int valueIndex = row * columnCount + column;
const double x = bounds.left() + spacingX * column;
surfacePoints->InsertNextPoint(x, y, 0.0);
scalarArray->InsertNextValue(interpolationValues[valueIndex]);
}
} }
vtkSmartPointer<vtkImageData> imageData = if(measurementPoints.size() == measurementValues.size()) {
vtkSmartPointer<vtkImageData>::New(); const double coordinateTolerance = qMax(
imageData->SetDimensions(columnCount, rowCount, 1); qMax(bounds.width(), bounds.height()) * 1.0e-10,
imageData->SetOrigin(bounds.left(), bounds.top(), 0.0); 1.0e-12);
imageData->SetSpacing(bounds.width() / (columnCount - 1), const int regularPointCount = surfacePoints->GetNumberOfPoints();
bounds.height() / (rowCount - 1), 1.0); for(int i = 0; i < measurementPoints.size(); ++i) {
imageData->GetPointData()->SetScalars(scalarArray); vtkIdType matchingPointId = -1;
for(int pointIndex = 0; pointIndex < regularPointCount; ++pointIndex) {
vtkSmartPointer<vtkImageDataGeometryFilter> geometryFilter = double regularPoint[3] = { 0.0, 0.0, 0.0 };
vtkSmartPointer<vtkImageDataGeometryFilter>::New(); surfacePoints->GetPoint(pointIndex, regularPoint);
geometryFilter->SetInputData(imageData); if(qAbs(regularPoint[0] - measurementPoints[i].x()) <= coordinateTolerance &&
geometryFilter->Update(); qAbs(regularPoint[1] - measurementPoints[i].y()) <= coordinateTolerance) {
matchingPointId = pointIndex;
// 默认使用完整规则网格；仅在裁剪结果有效且保留标量时切换. break;
vtkPolyData* visibleSurface = geometryFilter->GetOutput(); }
}
if(matchingPointId >= 0) {
scalarArray->SetValue(matchingPointId, measurementValues[i]);
}
else {
surfacePoints->InsertNextPoint(measurementPoints[i].x(),
measurementPoints[i].y(), 0.0);
scalarArray->InsertNextValue(measurementValues[i]);
}
}
}
vtkSmartPointer<vtkPolyData> surfaceData =
vtkSmartPointer<vtkPolyData>::New();
surfaceData->SetPoints(surfacePoints);
surfaceData->GetPointData()->SetScalars(scalarArray);
vtkSmartPointer<vtkDelaunay2D> triangulationFilter =
vtkSmartPointer<vtkDelaunay2D>::New();
triangulationFilter->SetInputData(surfaceData);
triangulationFilter->SetProjectionPlaneMode(VTK_DELAUNAY_XY_PLANE);
triangulationFilter->SetTolerance(0.0);
triangulationFilter->BoundingTriangulationOff();
triangulationFilter->Update();
// 默认使用完整三角网格；仅在裁剪结果有效且保留标量时切换.
vtkPolyData* visibleSurface = triangulationFilter->GetOutput();
vtkSmartPointer<vtkExtractPolyDataGeometry> extractFilter; vtkSmartPointer<vtkExtractPolyDataGeometry> extractFilter;
vtkSmartPointer<vtkImplicitSelectionLoop> selectionLoop; vtkSmartPointer<vtkImplicitSelectionLoop> selectionLoop;
vtkSmartPointer<vtkPoints> selectionPoints; vtkSmartPointer<vtkPoints> selectionPoints;
@ -188,7 +274,7 @@ void nmWxPropertyInterpolationPreviewDlg::initializeRenderer(
selectionLoop->SetNormal(0.0, 0.0, 1.0); selectionLoop->SetNormal(0.0, 0.0, 1.0);
extractFilter = vtkSmartPointer<vtkExtractPolyDataGeometry>::New(); extractFilter = vtkSmartPointer<vtkExtractPolyDataGeometry>::New();
extractFilter->SetInputConnection(geometryFilter->GetOutputPort()); extractFilter->SetInputConnection(triangulationFilter->GetOutputPort());
extractFilter->SetImplicitFunction(selectionLoop); extractFilter->SetImplicitFunction(selectionLoop);
extractFilter->ExtractInsideOn(); extractFilter->ExtractInsideOn();
extractFilter->ExtractBoundaryCellsOn(); extractFilter->ExtractBoundaryCellsOn();
@ -205,7 +291,7 @@ void nmWxPropertyInterpolationPreviewDlg::initializeRenderer(
if(visibleSurface == NULL || visibleSurface->GetNumberOfPoints() == 0 || if(visibleSurface == NULL || visibleSurface->GetNumberOfPoints() == 0 ||
visibleSurface->GetPointData() == NULL || visibleSurface->GetPointData() == NULL ||
visibleSurface->GetPointData()->GetScalars() == NULL) { visibleSurface->GetPointData()->GetScalars() == NULL) {
showPreviewFailure(renderer, renderWindow); showPreviewFailure(renderer, renderWindow, previewText("Preview Failed"));
return; return;
} }
@ -213,10 +299,23 @@ void nmWxPropertyInterpolationPreviewDlg::initializeRenderer(
visibleSurface->GetPointData()->GetScalars(); visibleSurface->GetPointData()->GetScalars();
visibleSurface->GetPointData()->SetScalars(interpolatedScalars); visibleSurface->GetPointData()->SetScalars(interpolatedScalars);
double displayRange[2] = { 0.0, 0.0 }; // 颜色条使用测点值范围，避免规则采样未命中测点时最大、最小值向内收缩.
visibleSurface->GetPointData()->GetScalars()->GetRange(displayRange); double displayMinimum = 0.0;
double displayMinimum = displayRange[0]; double displayMaximum = 0.0;
double displayMaximum = displayRange[1]; if(!measurementValues.isEmpty()) {
displayMinimum = measurementValues[0];
displayMaximum = measurementValues[0];
for(int i = 1; i < measurementValues.size(); ++i) {
displayMinimum = qMin(displayMinimum, measurementValues[i]);
displayMaximum = qMax(displayMaximum, measurementValues[i]);
}
}
else {
double displayRange[2] = { 0.0, 0.0 };
visibleSurface->GetPointData()->GetScalars()->GetRange(displayRange);
displayMinimum = displayRange[0];
displayMaximum = displayRange[1];
}
if(qAbs(displayMaximum - displayMinimum) < 1.0e-12) { if(qAbs(displayMaximum - displayMinimum) < 1.0e-12) {
const double padding = qMax(qAbs(displayMinimum) * 0.01, 1.0e-6); const double padding = qMax(qAbs(displayMinimum) * 0.01, 1.0e-6);
displayMinimum -= padding; displayMinimum -= padding;
@ -311,6 +410,19 @@ void nmWxPropertyInterpolationPreviewDlg::initializeRenderer(
vtkSmartPointer<vtkVertexGlyphFilter>::New(); vtkSmartPointer<vtkVertexGlyphFilter>::New();
glyphFilter->SetInputData(measurementData); glyphFilter->SetInputData(measurementData);
// 深色外层形成清晰描边，避免测点颜色与附近插值色带融为一体.
vtkSmartPointer<vtkPolyDataMapper> borderMapper =
vtkSmartPointer<vtkPolyDataMapper>::New();
borderMapper->SetInputConnection(glyphFilter->GetOutputPort());
borderMapper->ScalarVisibilityOff();
vtkSmartPointer<vtkActor> borderActor =
vtkSmartPointer<vtkActor>::New();
borderActor->SetMapper(borderMapper);
borderActor->GetProperty()->SetColor(0.10, 0.10, 0.10);
borderActor->GetProperty()->SetPointSize(9.0);
renderer->AddActor(borderActor);
vtkSmartPointer<vtkPolyDataMapper> measurementMapper = vtkSmartPointer<vtkPolyDataMapper> measurementMapper =
vtkSmartPointer<vtkPolyDataMapper>::New(); vtkSmartPointer<vtkPolyDataMapper>::New();
measurementMapper->SetInputConnection(glyphFilter->GetOutputPort()); measurementMapper->SetInputConnection(glyphFilter->GetOutputPort());
@ -321,7 +433,8 @@ void nmWxPropertyInterpolationPreviewDlg::initializeRenderer(
vtkSmartPointer<vtkActor> measurementActor = vtkSmartPointer<vtkActor> measurementActor =
vtkSmartPointer<vtkActor>::New(); vtkSmartPointer<vtkActor>::New();
measurementActor->SetMapper(measurementMapper); measurementActor->SetMapper(measurementMapper);
measurementActor->GetProperty()->SetPointSize(8.0); measurementActor->GetProperty()->SetPointSize(5.0);
measurementActor->SetPosition(0.0, 0.0, 0.001);
renderer->AddActor(measurementActor); renderer->AddActor(measurementActor);
} }

Write Preview
Loading…
Cancel
Save