Russian Qt Forum

QwtPlotCurve::closestPoint

double distance(const QPointF& p, const QPointF& a, const QPointF& b)
{
	const QPointF v = b - a;
	const QPointF w = p - a;

	const double c1 = v.x() * w.x() + v.y() * w.y();
	const double c2 = v.x() * v.x() + v.y() * v.y();

	const QPointF np = (c1 <= 0 ? a : c2 <= c1 ? b : a + c1 / c2 * v );

	return sqrt( (p.x() - np.x()) * (p.x() - np.x()) + (p.y() - np.y()) * (p.y() - np.y()) );
}

virtual QwtText trackerText(const QwtDoublePoint& point) const
{
	// Берём откуда-нибудь курву
	const MyPlotCurve* curve = ...;

	// Переходим к пиксельным координатам
	const QwtScaleMap xmap = curve->plot()->canvasMap(QwtPlot::xBottom);
	const QwtScaleMap ymap = curve->plot()->canvasMap(QwtPlot::yLeft);

	const QPointF p(xmap.transform(point.x()), ymap.transform(point.y()));

	// Запихиваем контрольные точки курвы в полигон
	QPolygonF polygon;

	const QwtData& data = curve->data();

	for(std::size_t i = 0; i < data.size(); ++i) {
		polygon.push_back(QPointF(xmap.transform(data.x(i)), ymap.transform(data.y(i))));
	}

	// Учитываем случай, когда курва интерполируется не линейно
	const bool      fitted = curve->testCurveAttribute(QwtPlotCurve::Fitted);
	QwtCurveFitter* fitter = curve->curveFitter();

	if(fitted && fitter) {
		polygon = fitter->fitCurve(polygon);
	}

	// Находим расстояние до курвы как минимум расстояний до всех отрезков курвы
	double mind = std::numeric_limits<double>::max();

	if(polygon.size() > 1) {
		for(int i = 0; i < polygon.size() - 1; ++i) {
			const QPointF& a = polygon[i];
			const QPointF& b = polygon[i + 1];

			const double d = distance(p, a, b);

			if(mind > d) {
				mind = d;
			}
		}
	}

	// Показываем расстояние во всплывающей подсказке
	return QString::number(static_cast<int>(mind));
}