Change epsg datum transform 9607 to 9606 by switching the sign of the…

… rotation parameters
qgis · Nov 25, 2013 · fc115a4 · fc115a4
1 parent 40a8cb9
commit fc115a4
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 0 deletions.
diff --git a/resources/srs.db b/resources/srs.db
diff --git a/src/core/qgscoordinatereferencesystem.cpp b/src/core/qgscoordinatereferencesystem.cpp
@@ -1996,6 +1996,15 @@ bool QgsCoordinateReferenceSystem::syncDatumTransform( const QString& dbPath )
     p7 = csList[csSize - 2];
     p7 = p7.isEmpty() ? "NULL" : p7;
 
+    //switch sign of rotation parameters. See http://trac.osgeo.org/proj/wiki/GenParms#towgs84-DatumtransformationtoWGS84
+    if ( coord_op_method == "9607" )
+    {
+      p4 = qgsDoubleToString( -( p4.toDouble() ) );
+      p5 = qgsDoubleToString( -( p5.toDouble() ) );
+      p6 = qgsDoubleToString( -( p6.toDouble() ) );
+      coord_op_method = "9606";
+    }
+
     //entry already in db?
     sqlite3_stmt* stmt;
     QString cOpCode;