2.7: coordinate deformation

Geo-Linux-Calculations · Nov 23, 2021 · e7b4945 · e7b4945
1 parent 4f7d015
commit e7b4945
Show file tree

Hide file tree

Showing 13 changed files with 661 additions and 151 deletions.
diff --git a/Makefile b/Makefile
@@ -3,8 +3,9 @@ PROGNAME1     = geofindkey
 PROGNAME2     = geoszbtoyxh
 PROGNAME3     = geositer500
 PROGNAME4     = geofindkey2p
+PROGNAME5     = geodeform500
 SRCS          = src
-PROGS         = $(PROGNAME1) $(PROGNAME2) $(PROGNAME3) $(PROGNAME4)
+PROGS         = $(PROGNAME1) $(PROGNAME2) $(PROGNAME3) $(PROGNAME4) $(PROGNAME5)
 CPP           = gcc
 CFLAGS        = -Wall -s
 LIBS          = -lm
@@ -33,6 +34,9 @@ $(PROGNAME3): $(SRCS)/$(PROGNAME3).c
 $(PROGNAME4): $(SRCS)/$(PROGNAME4).c
 	$(CPP) $(CFLAGS) $^ -o $@ $(LIBS)
 
+$(PROGNAME5): $(SRCS)/$(PROGNAME5).c
+	$(CPP) $(CFLAGS) $^ -o $@ $(LIBS)
+
 install: $(PROGS)
 	$(INSTALL) -d $(PREFIX)/bin
 	$(INSTALL) -m 0755 $(PROGS) $(PREFIX)/bin/

diff --git a/README.md b/README.md
@@ -220,7 +220,7 @@ OKD-12    3.889  288.39138889  133.60805556
 Exec finding unknown distances:
 ```sh
 ./geositer500 doc/szbunknow.dat szb.dat
-````
+```
 
 Result sbz.dat:
 ```
@@ -240,6 +240,53 @@ OKD-11 3.8890 71.60250000 313.60777778
 OKD-12 3.8890 288.39138889 133.60805556
 ```
 
+### Deformation:
+
+Input Data doc/gps.dat:
+```
+A             2303.443 5622.903 149.307 2302.854 5623.223 148.856 1.0
+B             2474.753 5830.133 114.434 2474.540 5829.768 114.844 1.0
+C             2748.516 5313.559 109.125 2748.808 5314.058 109.855 1.0
+(2a+3b+5c)/10 2577.372 5530.400 118.754
+(a+8b+c)/10   2484.998 5757.752 117.390
+(6a+2b+2c)/10 2426.719 5602.480 134.296
+a             2303.443 5622.903 149.307
+b             2474.753 5830.133 114.434
+c             2748.516 5313.559 109.125
+(a+b)/2       2389.098 5726.518 131.870
+(a+c)/2       2525.980 5468.231 129.216
+(b+c)/2       2611.635 5571.846 111.779
+A*            2302.854 5623.223 148.856
+B*            2474.540 5829.768 114.844
+C*            2748.808 5314.058 109.855
+```
+
+Exec finding unknown distances:
+```sh
+./geodeform500 doc/gps.dat gpsout.dat
+```
+
+Result gpsout.dat:
+```
+Deformation coordinate base: 3
+(CONST)
+A 2303.4430 5622.9030 149.3070 2302.8540 5623.2230 148.8560
+B 2474.7530 5830.1330 114.4340 2474.5400 5829.7680 114.8440
+C 2748.5160 5313.5590 109.1250 2748.8080 5314.0580 109.8550
+(2a+3b+5c)/10 2577.3720 5530.4000 118.7540 2577.2210 5530.5915 118.9920
+(a+8b+c)/10 2484.9980 5757.7520 117.3900 2484.7597 5757.4651 117.7270
+(6a+2b+2c)/10 2426.7190 5602.4800 134.2960 2426.2635 5602.6663 134.0994
+a 2303.4430 5622.9030 149.3070 2302.8540 5623.2230 148.8560
+b 2474.7530 5830.1330 114.4340 2474.5400 5829.7680 114.8440
+c 2748.5160 5313.5590 109.1250 2748.8080 5314.0580 109.8550
+(a+b)/2 2389.0980 5726.5180 131.8700 2388.7176 5726.5110 131.8718
+(a+c)/2 2525.9800 5468.2310 129.2160 2525.8176 5468.4730 129.4140
+(b+c)/2 2611.6350 5571.8460 111.7790 2611.4850 5571.9893 112.0412
+A* 2302.8540 5623.2230 148.8560 2302.2650 5623.5430 148.4050
+B* 2474.5400 5829.7680 114.8440 2474.3270 5829.4030 115.2540
+C* 2748.8080 5314.0580 109.8550 2749.1000 5314.5570 110.5850
+```
+
 ---
 
 https://github.com/Geo-Linux-Calculations/geofindkey
diff --git a/doc/gps.dat b/doc/gps.dat
@@ -0,0 +1,15 @@
+A             2303.443 5622.903 149.307 2302.854 5623.223 148.856 1.0
+B             2474.753 5830.133 114.434 2474.540 5829.768 114.844 1.0
+C             2748.516 5313.559 109.125 2748.808 5314.058 109.855 1.0
+(2a+3b+5c)/10 2577.372 5530.400 118.754
+(a+8b+c)/10   2484.998 5757.752 117.390
+(6a+2b+2c)/10 2426.719 5602.480 134.296
+a             2303.443 5622.903 149.307
+b             2474.753 5830.133 114.434
+c             2748.516 5313.559 109.125
+(a+b)/2       2389.098 5726.518 131.870
+(a+c)/2       2525.980 5468.231 129.216
+(b+c)/2       2611.635 5571.846 111.779
+A*            2302.854 5623.223 148.856
+B*            2474.540 5829.768 114.844
+C*            2748.808 5314.058 109.855
diff --git a/man/man1/geodeform500.1 b/man/man1/geodeform500.1
@@ -0,0 +1,82 @@
+.TH "GeoDeform500" 1 2.7 "22 Nov 2021" "User Manual"
+
+.SH NAME
+geodeform500
+
+.SH DESCRIPTION
+Deformation of coordinates by the "unit spheres" method. 
+
+.SH SYNOPSIS
+geodeform500 [option] input-file report-file
+
+.SH OPTIONS
+.TP
+-d N
+number of decimals in the result coordinates, default=4
+.TP
+-e N.N
+eps, default=0.000001
+.TP
+-m str
+base mode {CONST,VAR}, default=CONST
+.TP
+-h
+help
+
+.SH FILES
+.TP
+input-file
+(sample):
+ A             2303.443 5622.903 149.307 2302.854 5623.223 148.856 1.05
+ B             2474.753 5830.133 114.434 2474.540 5829.768 114.844 1.05
+ C             2748.516 5313.559 109.125 2748.808 5314.058 109.855 1.05
+ (2a+3b+5c)/10 2577.372 5530.400 118.7545
+ (a+8b+c)/10   2484.998 5757.752 117.3905
+ (6a+2b+2c)/10 2426.719 5602.480 134.2965
+ a             2303.443 5622.903 149.3075
+ b             2474.753 5830.133 114.4345
+ c             2748.516 5313.559 109.1255
+ (a+b)/2       2389.098 5726.518 131.8705
+ (a+c)/2       2525.980 5468.231 129.2165
+ (b+c)/2       2611.635 5571.846 111.7795
+ A*            2302.854 5623.223 148.8565
+ B*            2474.540 5829.768 114.8445
+ C*            2748.808 5314.058 109.8555
+.TP
+report-file
+(sample):
+ Deformation coordinate base: 3
+ (CONST)
+ A 2303.4430 5622.9030 149.3070 2302.8540 5623.2230 148.8560
+ B 2474.7530 5830.1330 114.4340 2474.5400 5829.7680 114.8440
+ C 2748.5160 5313.5590 109.1250 2748.8080 5314.0580 109.8550
+ (2a+3b+5c)/10 2577.3720 5530.4000 118.7540 2577.2210 5530.5915 118.9920
+ (a+8b+c)/10 2484.9980 5757.7520 117.3900 2484.7597 5757.4651 117.7270
+ (6a+2b+2c)/10 2426.7190 5602.4800 134.2960 2426.2635 5602.6663 134.0994
+ a 2303.4430 5622.9030 149.3070 2302.8540 5623.2230 148.8560
+ b 2474.7530 5830.1330 114.4340 2474.5400 5829.7680 114.8440
+ c 2748.5160 5313.5590 109.1250 2748.8080 5314.0580 109.8550
+ (a+b)/2 2389.0980 5726.5180 131.8700 2388.7176 5726.5110 131.8718
+ (a+c)/2 2525.9800 5468.2310 129.2160 2525.8176 5468.4730 129.4140
+ (b+c)/2 2611.6350 5571.8460 111.7790 2611.4850 5571.9893 112.0412
+ A* 2302.8540 5623.2230 148.8560 2302.2650 5623.5430 148.4050
+ B* 2474.5400 5829.7680 114.8440 2474.3270 5829.4030 115.2540
+ C* 2748.8080 5314.0580 109.8550 2749.1000 5314.5570 110.5850
+
+.SH EXAMPLE
+geodeform500 doc/gps.dat gpsout.dat
+
+.SH COPYRIGHT
+Copyright 2018-2021 zvezdochiot.
+ All rights reserved.
+
+.SH SEE ALSO
+ geoeasy
+ geofindkey
+ geofindkey2p
+ geositer500
+ geoszbtoyxh
+
+.SH CONTACTS
+ Website: https://github.com/Geo-Linux-Calculations/geofindkey
+ Email: [email protected]
diff --git a/man/man1/geofindkey.1 b/man/man1/geofindkey.1
@@ -1,4 +1,4 @@
-.TH "GeoFindKey" 1 2.6 "09 Nov 2021" "User Manual"
+.TH "GeoFindKey" 1 2.7 "22 Nov 2021" "User Manual"
 
 .SH NAME
 geofindkey
@@ -74,6 +74,7 @@ Copyright 2018-2021 zvezdochiot.
 
 .SH SEE ALSO
  geoeasy
+ geodeform500
  geofindkey2p
  geoszbtoyxh
  geositer500

diff --git a/man/man1/geofindkey2p.1 b/man/man1/geofindkey2p.1
@@ -1,4 +1,4 @@
-.TH "GeoFindKey2Pow" 1 2.6 "09 Nov 2021" "User Manual"
+.TH "GeoFindKey2Pow" 1 2.7 "22 Nov 2021" "User Manual"
 
 .SH NAME
 geofindkey2p
@@ -71,6 +71,7 @@ Copyright 2018-2021 zvezdochiot.
 
 .SH SEE ALSO
  geoeasy
+ geodeform500
  geofindkey
  geoszbtoyxh
  geositer500

diff --git a/man/man1/geositer500.1 b/man/man1/geositer500.1
@@ -1,4 +1,4 @@
-.TH "GeoSIter500" 1 2.6 "09 Nov 2021" "User Manual"
+.TH "GeoSIter500" 1 2.7 "22 Nov 2021" "User Manual"
 
 .SH NAME
 geositer500
@@ -68,6 +68,7 @@ Copyright 2018-2021 zvezdochiot.
 
 .SH SEE ALSO
  geoeasy
+ geodeform500
  geofindkey
  geofindkey2p
  geoszbtoyxh

diff --git a/man/man1/geoszbtoyxh.1 b/man/man1/geoszbtoyxh.1
@@ -1,4 +1,4 @@
-.TH "GeoSZBtoYXH" 1 2.6 "09 Nov 2021" "User Manual"
+.TH "GeoSZBtoYXH" 1 2.7 "22 Nov 2021" "User Manual"
 
 .SH NAME
 geoszbtoyxh
@@ -68,6 +68,7 @@ Copyright 2018-2021 zvezdochiot.
 
 .SH SEE ALSO
  geoeasy
+ geodeform500
  geofindkey
  geofindkey2p
  geositer500