QGIS » QGIS Application

        /*! Update proj.4 parameters in our database from proj.4

         * @returns number of updated CRS on success and

         *   negative number of failed updates in case of errors.

         * @note added in 1.7

         */

        static int syncDb();

usr/lib/qgis/crssync

#!/bin/sh

set -e

case "$1" in

	configure)

		/usr/lib/qgis/crssync

		;;

esac	

/***************************************************************************

                             crssync.cpp

                             sync srs.db with proj

                             -------------------

    begin                : 2011

    copyright            : (C) 2011 by Juergen E. Fischer, norBIT GmbH

    email                : jef at norbit dot de

***************************************************************************/

/***************************************************************************

 *                                                                         *

 *   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 "qgsapplication.h"

#include "qgscoordinatereferencesystem.h"

#include "qgsconfig.h"

#include <QRegExp>

#include <iostream>

#include <limits>

#include <cpl_error.h>

void CPL_STDCALL showError( CPLErr errClass, int errNo, const char *msg )

{

  QRegExp re( "EPSG PCS/GCS code \\d+ not found in EPSG support files.  Is this a valid\nEPSG coordinate system?" );

  if ( errNo != 6 && !re.exactMatch( msg ) )

  {

    std::cerr << msg;

  }

}

int main( int argc, char ** argv )

{

  QgsApplication a( argc, argv, false );

#if defined(Q_WS_MACX)

  // If we're on Mac, we have the resource library way above us...

  a.setPkgDataPath( QgsApplication::prefixPath() + "/../../../../" + QString( QGIS_DATA_SUBDIR ) );

#elif defined(Q_WS_WIN)

  a.setPkgDataPath( QgsApplication::prefixPath() + "/" QGIS_DATA_SUBDIR );

#else

  a.setPkgDataPath( QgsApplication::prefixPath() + "/../" QGIS_DATA_SUBDIR );

#endif

  std::cout << "Synchronizing CRS database with PROJ definitions." << std::endl;

  CPLPushErrorHandler( showError );

  int res = QgsCoordinateReferenceSystem::syncDb();

  CPLPopErrorHandler();

  if ( res == 0 )

  {

    std::cout << "No CRS updates were necessary." << std::endl;

  }

  else if ( res > 0 )

  {

    std::cout << res << " CRSs updated." << std::endl;

  }

  else if ( res == std::numeric_limits<int>::min() )

  {

    std::cout << "CRSs synchronization not possible." << std::endl;

  }

  else if ( res < 0 )

  {

    std::cout << -res << " CRSs could not be updated." << std::endl;

  }

  return 0;

}

#include <proj_api.h>

int QgsCoordinateReferenceSystem::syncDb()

{

  int updated = 0, errors = 0;

  sqlite3 *database;

  if ( sqlite3_open( QgsApplication::srsDbFilePath().toUtf8().constData(), &database ) != SQLITE_OK )

  {

    qCritical( "Can't open database: %s [%s]\n", QgsApplication::srsDbFilePath().toLocal8Bit().constData(), sqlite3_errmsg( database ) );

    return -1;

  }

  const char *tail;

  sqlite3_stmt *select;

  QString sql = "select auth_name,auth_id,parameters from tbl_srs WHERE auth_name IS NOT NULL AND auth_id IS NOT NULL";

  if ( sqlite3_prepare( database, sql.toAscii(), sql.size(), &select, &tail ) != SQLITE_OK )

  {

    qCritical( "Could not prepare: %s [%s]\n", sql.toAscii().constData(), sqlite3_errmsg( database ) );

    sqlite3_close( database );

    return -1;

  }

  OGRSpatialReferenceH crs = OSRNewSpatialReference( NULL );

  while ( sqlite3_step( select ) == SQLITE_ROW )

  {

    const char *auth_name = ( const char * ) sqlite3_column_text( select, 0 );

    const char *auth_id   = ( const char * ) sqlite3_column_text( select, 1 );

    const char *params    = ( const char * ) sqlite3_column_text( select, 2 );

    QString proj4;

    QString input = QString( "%1:%2" ).arg( QString( auth_name ).toLower() ).arg( auth_id );

    QgsDebugMsgLevel( QString( "fetching crs %1" ).arg( input ), 3 );

#if !defined(PJ_VERSION) || PJ_VERSION!=470

    // 4.7.0 has a bug that crashes after 16 consecutive pj_init_plus with different strings

    if ( proj4.isEmpty() )

    {

      input = QString( "+init=%1:%2" ).arg( QString( auth_name ).toLower() ).arg( auth_id );

      projPJ pj = pj_init_plus( input.toAscii() );

      if ( !pj )

      {

        input = QString( "+init=%1:%2" ).arg( QString( auth_name ).toUpper() ).arg( auth_id );

        pj = pj_init_plus( input.toAscii() );

      }

      if ( pj )

      {

        char *def = pj_get_def( pj, 0 );

        if ( def )

        {

          proj4 = def;

          pj_dalloc( def );

          input.prepend( ' ' ).append( ' ' );

          if ( proj4.startsWith( input ) )

          {

            proj4 = proj4.mid( input.size() );

          }

        }

        else

        {

          QgsDebugMsg( QString( "could not retrieve proj string for %1 from PROJ" ).arg( input ) );

        }

      }

      else

      {

        QgsDebugMsgLevel( QString( "could not retrieve crs for %1 from PROJ" ).arg( input ), 3 );

      }

      pj_free( pj );

    }

#else

    // otherwise check OGR

    if ( !proj4.isEmpty() )

    {

      OGRErr ogrErr = OSRSetFromUserInput( crs, input.toAscii() );

      if ( ogrErr != OGRERR_NONE )

      {

        input = QString( "%1:%2" ).arg( QString( auth_name ).toUpper() ).arg( auth_id );

        ogrErr = OSRSetFromUserInput( crs, input.toAscii() );

      }

      if ( ogrErr == OGRERR_NONE )

      {

        char *output = 0;

        if ( OSRExportToProj4( crs, &output ) == OGRERR_NONE )

        {

          QString proj4ogr = output;

#if 0

          QString proj4proj = proj4;

          proj4proj.replace( QRegExp( "\\.0+([, ])" ), "\\1" );

          proj4proj.replace( QRegExp( "(\\.\\d*[1-9])0+([, ])" ), "\\1\\2" );

          proj4proj.replace( ".0000000001 ", " " );

          proj4proj.replace( ".10000000000001 ", ".1 " );

          proj4proj.replace( ".50000000000001 ", ".5 " );

          proj4proj.replace( ".399999999999999 ", ".4 " );

          proj4proj.replace( ".16500000000001 ", ".165 " );

          proj4proj.replace( " +towgs84=0,0,0,0,0,0,0 ", " +datum=WGS84 " );

          if ( proj4proj.contains( " +datum=WGS84 " ) && proj4proj.endsWith( " +towgs84=0,0,0" ) )

          {

            proj4proj.replace( " +towgs84=0,0,0", "" );

          }

          proj4ogr.replace( " +towgs84=0,0,0,0,0,0,0 ", " +datum=WGS84 " );

          proj4ogr.replace( " +a=6378249.145 +b=6356514.96582849 ", " +ellps=clrk80 " );

          proj4ogr.replace( " +gamma=0 ", " " );

          proj4ogr.replace( " +alpha=0 ", " " );

          if ( proj4ogr != proj4proj )

          {

            qWarning( "%s", QString( "Mismatch for %1:%2 in OGR and PROJ found:\n  PROJ:'%3'\n  OGR :'%4'" ).arg( auth_name ).arg( auth_id ).arg( proj4proj ).arg( proj4ogr ).toAscii().constData() );

          }

#endif

          proj4 = proj4ogr.trimmed();

        }

        else

        {

          QgsDebugMsg( QString( "could not retrieve proj.4 string for %1 from OGR" ).arg( input ) );

        }

        if ( output )

          CPLFree( output );

      }

    }

#endif

    if ( proj4.isEmpty() )

    {

      continue;

    }

    if ( proj4 != params )

    {

      char *errMsg = NULL;

      sql = QString( "UPDATE tbl_srs SET parameters=%1 WHERE auth_name=%2 AND auth_id=%3" )

            .arg( quotedValue( proj4 ) )

            .arg( quotedValue( auth_name ) )

            .arg( quotedValue( auth_id ) );

      if ( sqlite3_exec( database, sql.toUtf8(), 0, 0, &errMsg ) != SQLITE_OK )

      {

        qCritical( "Could not execute: %s [%s/%s]\n",

                   sql.toLocal8Bit().constData(),

                   sqlite3_errmsg( database ),

                   errMsg ? errMsg : "(unknown error)" );

        errors++;

      }

      else

      {

        updated++;

        QgsDebugMsgLevel( "Updated: " + sql, 3 );

      }

      if ( errMsg )

        sqlite3_free( errMsg );

    }

  }

  OSRDestroySpatialReference( crs );

  sqlite3_finalize( select );

  sqlite3_close( database );

  if ( errors > 0 )

    return -errors;

  else

    return updated;

}

ADD_EXECUTABLE(crssync main.cpp)

INCLUDE_DIRECTORIES(

  ../core

  ${GDAL_INCLUDE_DIR}

  ${PROJ_INCLUDE_DIR}

)

TARGET_LINK_LIBRARIES(crssync

  qgis_core

  ${PROJ_LIBRARY}

  ${GDAL_LIBRARY}

)

INSTALL(CODE "MESSAGE(\"Installing crssync ...\")")

INSTALL(TARGETS crssync RUNTIME DESTINATION ${QGIS_LIBEXEC_DIR})

    /*! Update proj.4 parameters in our database from proj.4

     * @returns number of updated CRS on success and

     *   negative number of failed updates in case of errors.

     * @note added in 1.7

     */

    static int syncDb();

    static QString quotedValue( QString value );

SUBDIRS(astyle core analysis ui gui app providers plugins helpviewer crssync)

start "CRS synchronization" /wait %OSGEO4W_ROOT%\apps\@package@\crssync