Branch data     Line data    Source code

       1                 :            : /***************************************************************************
       2                 :            :                              qgsbearingutils.cpp
       3                 :            :                              -------------------
       4                 :            :     begin                : October 2016
       5                 :            :     copyright            : (C) 2016 by Nyall Dawson
       6                 :            :     email                : nyall dot dawson at gmail dot com
       7                 :            :  ***************************************************************************/
       8                 :            : 
       9                 :            : /***************************************************************************
      10                 :            :  *                                                                         *
      11                 :            :  *   This program is free software; you can redistribute it and/or modify  *
      12                 :            :  *   it under the terms of the GNU General Public License as published by  *
      13                 :            :  *   the Free Software Foundation; either version 2 of the License, or     *
      14                 :            :  *   (at your option) any later version.                                   *
      15                 :            :  *                                                                         *
      16                 :            :  ***************************************************************************/
      17                 :            : 
      18                 :            : #include "qgsbearingutils.h"
      19                 :            : #include "qgscoordinatereferencesystem.h"
      20                 :            : #include "qgscoordinatetransformcontext.h"
      21                 :            : #include "qgspointxy.h"
      22                 :            : #include "qgscoordinatetransform.h"
      23                 :            : #include "qgsexception.h"
      24                 :            : 
      25                 :          0 : double QgsBearingUtils::bearingTrueNorth( const QgsCoordinateReferenceSystem &crs, const QgsCoordinateTransformContext  &transformContext, const QgsPointXY &point )
      26                 :            : {
      27                 :            :   // step 1 - transform point into WGS84 geographic crs
      28                 :          0 :   QgsCoordinateTransform transform( crs, QgsCoordinateReferenceSystem::fromEpsgId( 4326 ), transformContext );
      29                 :            : 
      30                 :          0 :   if ( !transform.isValid() )
      31                 :            :   {
      32                 :            :     //raise
      33                 :          0 :     throw QgsException( QObject::tr( "Could not create transform to calculate true north" ) );
      34                 :            :   }
      35                 :            : 
      36                 :          0 :   if ( transform.isShortCircuited() )
      37                 :          0 :     return 0.0;
      38                 :            : 
      39                 :          0 :   QgsPointXY p1 = transform.transform( point );
      40                 :            : 
      41                 :            :   // shift point a tiny bit north
      42                 :          0 :   QgsPointXY p2 = p1;
      43                 :          0 :   p2.setY( p2.y() + 0.000001 );
      44                 :            : 
      45                 :            :   //transform back
      46                 :          0 :   QgsPointXY p3 = transform.transform( p2, QgsCoordinateTransform::ReverseTransform );
      47                 :            : 
      48                 :            :   // find bearing from point to p3
      49                 :          0 :   return point.azimuth( p3 );
      50                 :          0 : }