1.0.6: Proj string output

Geo-Linux-Calculations · Jun 26, 2022 · c2a61a2 · c2a61a2 · zvezdochiot · Jun 26, 2022
1 parent d394723
commit c2a61a2
Show file tree

Hide file tree

Showing 3 changed files with 13 additions and 5 deletions.
diff --git a/src/helmeltrans.c b/src/helmeltrans.c
@@ -85,7 +85,7 @@ int main(int argc, char* argv[])
         exit(EXIT_FAILURE);
     }
     command = argv[1];
-    if (!strcmp(command,"xyz")) 
+    if (!strcmp(command,"xyz"))
     {
         fprintf(stdout,"B-L-H -> X-Y-Z...\n");
     }
@@ -124,7 +124,7 @@ int main(int argc, char* argv[])
     }
     else
         ofile = stdout;
-    
+
     if(ofile == NULL)
     {
         fprintf(stderr,"Error writing %s\n",ofilename);
@@ -144,7 +144,7 @@ int main(int argc, char* argv[])
     b2 = b * b;
     ro = M_PI / 180.0;
     roi = 180.0 / M_PI;
-    if (!strcmp(command,"xyz")) 
+    if (!strcmp(command,"xyz"))
     {
         while(fgets( ibuf, 256, ifile)!=NULL)
         {

diff --git a/src/helmert3d.h b/src/helmert3d.h
@@ -22,9 +22,10 @@
 #include <stdlib.h>
 #include <string.h>
 #include <math.h>
-#define VERS "1.0.5"
+#define VERS "1.0.6"
 #define C_AUTHOR "U. Niethammer"
 #define C_YEAR "2021"
 #define DEBUG 0
+#define RAD2ARC 206264.8062
 
 #endif // HELMERT3D_H
diff --git a/src/helmparms3d.c b/src/helmparms3d.c
@@ -288,6 +288,7 @@ int main(int argc, char* argv[])
     double **R_mat=NULL;
     double trace1=0.0, trace2=0.0;
     double scal=0.0, ppm=0.0;
+    double rx=0.0, ry=0.0, rz=0.0;
     FILE *outfile;
 
     fprintf(stdout,"\n*******************************\n");
@@ -514,7 +515,7 @@ int main(int argc, char* argv[])
         fprintf(stderr,"Error writing %s\n",out_param_name);
         exit(EXIT_FAILURE);
     }
-    
+
     matrix_init(n, m, src_mat_T);
     matrix_transpose(n, n, R_mat, src_mat_T);
     fprintf(stdout,"R =\n");
@@ -529,6 +530,12 @@ int main(int argc, char* argv[])
 
     fprintf(stdout,"s = %10.10lf (= %10.10lf ppm)\n\n",scal, ppm);
 
+    rx = 0.5 * (src_mat_T[2][1] - src_mat_T[1][2]) * RAD2ARC;
+    ry = 0.5 * (src_mat_T[0][2] - src_mat_T[2][0]) * RAD2ARC;
+    rz = 0.5 * (src_mat_T[1][0] - src_mat_T[0][1]) * RAD2ARC;
+    fprintf(stdout,"Proj string:\n");
+    fprintf(stdout,"+towgs84=%.3f,%.3f,%.3f,%.3f,%.3f,%.3f,%.3f\n\n", T_vec[0], T_vec[1], T_vec[2], rx, ry, rz, ppm);
+
     if(argc < 4)
     {
         fprintf(stdout,"Results matrix:\n");