qgscoordinateutils.cpp

Go to the documentation of this file.

/***************************************************************************
                             qgscoordinateutils.cpp
                             ----------------------
    begin                : February 2016
    copyright            : (C) 2016 by Nyall Dawson
    email                : nyall dot dawson at gmail dot com
 ***************************************************************************/
 
/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/
 
#include "qgscoordinateutils.h"
 
#include "qgis.h"
#include "qgscoordinateformatter.h"
#include "qgscoordinatenumericformat.h"
#include "qgscoordinatereferencesystem.h"
#include "qgscoordinatereferencesystemutils.h"
#include "qgscoordinatetransform.h"
#include "qgsexception.h"
#include "qgsproject.h"
#include "qgsprojectdisplaysettings.h"
#include "qgsrectangle.h"
 
#include <QLocale>
#include <QRegularExpression>
#include <QString>
 
#include "moc_qgscoordinateutils.cpp"
 
using namespace Qt::StringLiterals;

 
int QgsCoordinateUtils::calculateCoordinatePrecision( double mapUnitsPerPixel, const QgsCoordinateReferenceSystem &mapCrs, QgsProject *project )
{
  if ( !project )
    project = QgsProject::instance(); // skip-keyword-check
  // Get the display precision from the project settings
  const bool automatic = project->readBoolEntry( u"PositionPrecision"_s, u"/Automatic"_s );
  int dp = 0;
 
  if ( automatic )
  {
    const bool formatGeographic = project->displaySettings()->coordinateType() == Qgis::CoordinateDisplayType::MapGeographic ||
                                  ( project->displaySettings()->coordinateType() == Qgis::CoordinateDisplayType::CustomCrs &&
                                    project->displaySettings()->coordinateCustomCrs().isGeographic() );
 
    // we can only calculate an automatic precision if one of these is true:
    // - both map CRS and format are geographic
    // - both map CRS and format are not geographic
    // - map CRS is geographic but format is not geographic (i.e. map units)
    if ( mapCrs.isGeographic() || !formatGeographic )
    {
      // Work out a suitable number of decimal places for the coordinates with the aim of always
      // having enough decimal places to show the difference in position between adjacent pixels.
      // Also avoid taking the log of 0.
      if ( !qgsDoubleNear( mapUnitsPerPixel, 0.0 ) )
        dp = static_cast<int>( std::ceil( -1.0 * std::log10( mapUnitsPerPixel ) ) );
    }
    else
    {
      switch ( project->displaySettings()->geographicCoordinateFormat()->angleFormat() )
      {
        case QgsGeographicCoordinateNumericFormat::AngleFormat::DegreesMinutesSeconds:
        case QgsGeographicCoordinateNumericFormat::AngleFormat::DegreesMinutes:
          dp = 2;
          break;
        case QgsGeographicCoordinateNumericFormat::AngleFormat::DecimalDegrees:
          dp = 4;
          break;
      }
    }
  }
  else
    dp = project->readNumEntry( u"PositionPrecision"_s, u"/DecimalPlaces"_s );
 
  // Keep dp sensible
  if ( dp < 0 )
    dp = 0;
 
  return dp;
}
 
int QgsCoordinateUtils::calculateCoordinatePrecisionForCrs( const QgsCoordinateReferenceSystem &crs, QgsProject *project )
{
  QgsProject *prj = project;
  if ( !prj )
  {
    prj = QgsProject::instance(); // skip-keyword-check
  }
 
  const bool automatic = prj->readBoolEntry( u"PositionPrecision"_s, u"/Automatic"_s );
  if ( !automatic )
  {
    return prj->readNumEntry( u"PositionPrecision"_s, u"/DecimalPlaces"_s, 6 );
  }
 
  return calculateCoordinatePrecision( crs );
}
 
int QgsCoordinateUtils::calculateCoordinatePrecision( const QgsCoordinateReferenceSystem &crs )
{
  const Qgis::DistanceUnit unit = crs.mapUnits();
  if ( unit == Qgis::DistanceUnit::Degrees )
  {
    return 8;
  }
  else
  {
    return 3;
  }
}
 
QString QgsCoordinateUtils::formatCoordinateForProject( QgsProject *project, const QgsPointXY &point, const QgsCoordinateReferenceSystem &destCrs, int precision )
{
  if ( !project )
    return QString();
 
  QString formattedX;
  QString formattedY;
  formatCoordinatePartsForProject( project, point, destCrs, precision, formattedX, formattedY );
 
  if ( formattedX.isEmpty() || formattedY.isEmpty() )
    return QString();
 
  const Qgis::CoordinateOrder axisOrder = project->displaySettings()->coordinateAxisOrder();
 
  QgsCoordinateReferenceSystem crs = project->displaySettings()->coordinateCrs();
  if ( !crs.isValid() && !destCrs.isValid() )
  {
    return u"%1%2 %3"_s.arg( formattedX, QgsCoordinateFormatter::separator(), formattedY );
  }
  else if ( !crs.isValid() )
  {
    crs = destCrs;
  }
 
  const Qgis::CoordinateOrder order = axisOrder == Qgis::CoordinateOrder::Default ? QgsCoordinateReferenceSystemUtils::defaultCoordinateOrderForCrs( crs ) : axisOrder;
  switch ( order )
  {
    case Qgis::CoordinateOrder::Default:
    case Qgis::CoordinateOrder::XY:
      return u"%1%2 %3"_s.arg( formattedX, QgsCoordinateFormatter::separator(), formattedY );
 
    case Qgis::CoordinateOrder::YX:
      return u"%1%2 %3"_s.arg( formattedY, QgsCoordinateFormatter::separator(), formattedX );
  }
  BUILTIN_UNREACHABLE
}
 
void QgsCoordinateUtils::formatCoordinatePartsForProject( QgsProject *project, const QgsPointXY &point, const QgsCoordinateReferenceSystem &destCrs, int precision, QString &x, QString &y )
{
  x.clear();
  y.clear();
  if ( !project )
    return;
 
  QgsCoordinateReferenceSystem crs = project->displaySettings()->coordinateCrs();
  if ( !crs.isValid() && !destCrs.isValid() )
  {
    x = QgsCoordinateFormatter::formatAsPair( point.x(), precision );
    y = QgsCoordinateFormatter::formatAsPair( point.y(), precision );
    return;
  }
  else if ( !crs.isValid() )
  {
    crs = destCrs;
  }
 
  QgsPointXY p = point;
  const bool isGeographic = crs.isGeographic();
  if ( destCrs  != crs )
  {
    const QgsCoordinateTransform ct( destCrs, crs, project );
    try
    {
      p = ct.transform( point );
    }
    catch ( QgsCsException & )
    {
      return;
    }
  }
 
  if ( isGeographic )
  {
    std::unique_ptr< QgsGeographicCoordinateNumericFormat > format( project->displaySettings()->geographicCoordinateFormat()->clone() );
    format->setNumberDecimalPlaces( precision );
 
    QgsNumericFormatContext context;
    context.setInterpretation( QgsNumericFormatContext::Interpretation::Longitude );
    x = format->formatDouble( p.x(), context );
    context.setInterpretation( QgsNumericFormatContext::Interpretation::Latitude );
    y = format->formatDouble( p.y(), context );
  }
  else
  {
    // coordinates in map units
    x = QgsCoordinateFormatter::formatAsPair( p.x(), precision );
    y = QgsCoordinateFormatter::formatAsPair( p.y(), precision );
  }
}
 
QString QgsCoordinateUtils::formatExtentForProject( QgsProject *project, const QgsRectangle &extent, const QgsCoordinateReferenceSystem &destCrs, int precision )
{
  const QgsPointXY p1( extent.xMinimum(), extent.yMinimum() );
  const QgsPointXY p2( extent.xMaximum(), extent.yMaximum() );
  return u"%1 : %2"_s.arg( QgsCoordinateUtils::formatCoordinateForProject( project, p1, destCrs, precision ),
                           QgsCoordinateUtils::formatCoordinateForProject( project, p2, destCrs, precision ) );
}
 
double QgsCoordinateUtils::degreeToDecimal( const QString &string, bool *ok, bool *isEasting )
{
  const QString negative( u"swSW"_s );
  const QString easting( u"eEwW"_s );
  double value = 0.0;
  bool okValue = false;
 
  if ( ok )
  {
    *ok = false;
  }
  else
  {
    ok = &okValue;
  }
 
  const QLocale locale;
  QRegularExpression degreeWithSuffix( u"^\\s*([-]?\\d{1,3}(?:[\\.\\%1]\\d+)?)\\s*([NSEWnsew])\\s*$"_s
                                       .arg( locale.decimalPoint() ) );
  QRegularExpressionMatch match = degreeWithSuffix.match( string );
  if ( match.hasMatch() )
  {
    const QString suffix = match.captured( 2 );
    value = std::abs( match.captured( 1 ).toDouble( ok ) );
    if ( *ok == false )
    {
      value = std::abs( locale.toDouble( match.captured( 1 ), ok ) );
    }
    if ( *ok )
    {
      value *= ( negative.contains( suffix ) ? -1 : 1 );
      if ( isEasting )
      {
        *isEasting = easting.contains( suffix );
      }
    }
  }
  return value;
}
 
double QgsCoordinateUtils::dmsToDecimal( const QString &string, bool *ok, bool *isEasting )
{
  const QString negative( u"swSW-"_s );
  const QString easting( u"eEwW"_s );
  double value = 0.0;
  bool okValue = false;
 
  if ( ok )
  {
    *ok = false;
  }
  else
  {
    ok = &okValue;
  }
 
  const QLocale locale;
  const QRegularExpression dms( u"^\\s*(?:([-+nsew])\\s*)?(\\d{1,3})(?:[^0-9.]+([0-5]?\\d))?[^0-9.]+([0-5]?\\d(?:[\\.\\%1]\\d+)?)[^0-9.,]*?([-+nsew])?\\s*$"_s
                                .arg( locale.decimalPoint() ), QRegularExpression::CaseInsensitiveOption );
  const QRegularExpressionMatch match = dms.match( string.trimmed() );
  if ( match.hasMatch() )
  {
    const QString dms1 = match.captured( 2 );
    const QString dms2 = match.captured( 3 );
    const QString dms3 = match.captured( 4 );
 
    double v = dms3.toDouble( ok );
    if ( *ok == false )
    {
      v = locale.toDouble( dms3, ok );
      if ( *ok == false )
        return value;
    }
    // Allow for Degrees/minutes format as well as DMS
    if ( !dms2.isEmpty() )
    {
      v = dms2.toInt( ok ) + v / 60.0;
      if ( *ok == false )
        return value;
    }
    v = dms1.toInt( ok ) + v / 60.0;
    if ( *ok == false )
      return value;
 
    const QString sign1 = match.captured( 1 );
    const QString sign2 = match.captured( 5 );
 
    if ( sign1.isEmpty() )
    {
      value = !sign2.isEmpty() && negative.contains( sign2 ) ? -v : v;
      if ( isEasting )
      {
        *isEasting = easting.contains( sign2 );
      }
    }
    else if ( sign2.isEmpty() )
    {
      value = !sign1.isEmpty() && negative.contains( sign1 ) ? -v : v;
      if ( isEasting )
      {
        *isEasting = easting.contains( sign1 );
      }
    }
    else
    {
      *ok = false;
    }
  }
  return value;
}