summaryrefslogtreecommitdiff
path: root/src/panorama.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/panorama.c')
-rw-r--r--src/panorama.c386
1 files changed, 0 insertions, 386 deletions
diff --git a/src/panorama.c b/src/panorama.c
deleted file mode 100644
index 2c58c65..0000000
--- a/src/panorama.c
+++ /dev/null
@@ -1,386 +0,0 @@
-/*
- * "$Id: panorama.c,v 1.8 2005/04/12 19:56:42 hofmann Exp $"
- *
- * flpsed program.
- *
- * Copyright 2004,2005 by Johannes Hofmann
- *
- * 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.
- *
- * This program is distributed in the hope that it will be useful,
- * but WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
- * General Public License for more details.
- *
- * You should have received a copy of the GNU General Public
- * License along with this library; if not, write to the Free Software
- * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
- * USA.
- */
-
-
-#include <stdio.h>
-#include <unistd.h>
-#include <stdlib.h>
-#include <string.h>
-#include <math.h>
-
-double pi, deg2rad;
-
-double distance(double phi_a, double lam_a, double phi_b, double lam_b);
-double sin_alpha(double lam_a, double lam_b, double phi_b, double c);
-double cos_alpha(double phi_a, double phi_b, double c);
-double alpha(double phi_a, double lam_a, double phi_b, double lam_b);
-double center_angle(double alph_a, double alph_b, double d1, double d2);
-int read_file(double phi_a, double lam_a, char *name,
- int dist, double h_d_ratio, double mm_per_deg,
- double a_center);
-int
-main(int argc, char **argv) {
- char c;
- char *sep, *tmp, *data_file = NULL;
- int errflg = 0;
- int dist = 100000;
- double phi_a, lam_a, phi_b, lam_b;
- double h_d_ratio = 0.1, mm_per_deg = 200.0, dist_c_m1 = 0.0;
- double dist_c_m2 = 0.0, a_center = 0.0;
- char *pos = NULL, *m1 = NULL, *m2 = NULL, *m_c = NULL;
-
- pi = asin(1.0) * 2.0;
- deg2rad = pi / 180.0;
-
- while ((c = getopt(argc, argv, "a:b:f:d:r:p:m:n:s:t:c:")) != EOF) {
- switch (c) {
- case 'a':
- tmp = strdup(optarg);
- if ((sep = strchr(tmp, ',')) == NULL) {
- errflg++;
- break;
- }
- *sep = '\0';
- phi_a = atof(tmp) * deg2rad;
- lam_a = atof(sep+1) * deg2rad;
- free(tmp);
- break;
- case 'b':
- tmp = strdup(optarg);
- if ((sep = strchr(tmp, ',')) == NULL) {
- usage();
- }
- *sep = '\0';
- phi_b = atof(tmp) * deg2rad;
- lam_b = atof(sep+1) * deg2rad;
- free(tmp);
- break;
- case 'f':
- data_file = optarg;
- break;
- case 'd':
- dist = atoi(optarg);
- break;
- case 'r':
- h_d_ratio = atof(optarg);
- break;
- case 's':
- dist_c_m1 = -atof(optarg);
- break;
- case 't':
- dist_c_m2 = -atof(optarg);
- break;
- case 'p':
- pos = optarg;
- break;
- case 'm':
- m1 = optarg;
- break;
- case 'n':
- m2 = optarg;
- break;
- case 'c':
- m_c = optarg;
- break;
- }
-}
-
-#if 0
- {
- double a1 = 0.0, a2 = 90.0;
- double d1 = -20.0, d2 = 10.0;
-
- fprintf(stderr, "==> %f\n", center_angle(a1*deg2rad, a2*deg2rad, d1, d2) / deg2rad);
- exit(1);
- }
-#endif
-
- if (errflg) {
- usage();
- exit(1);
- }
-
- if (!data_file) {
- usage();
- exit(1);
- }
-
- if (pos){
- int n = get_pos(data_file, pos, &phi_a, &lam_a);
- if (n == 0) {
- fprintf(stderr, "No matching entry for %s in datafile.\n",
- pos);
- exit(1);
- } else if (n > 1) {
- fprintf(stderr, "More than one matching entry for %s in datafile.\n",
- pos);
- exit(1);
- }
- }
-
- if (m1 && m2) {
- double phi, lam;
- double a1, a2, a;
- int n;
-
- n = get_pos(data_file, m1, &phi, &lam);
- if (n == 0) {
- fprintf(stderr, "No matching entry for %s in datafile.\n",
- m1);
- exit(1);
- } else if (n > 1) {
- fprintf(stderr, "More than one matching entry for %s in datafile.\n",
- m1);
- exit(1);
- }
- a1 = alpha(phi_a, lam_a, phi, lam);
-
- n = get_pos(data_file, m2, &phi, &lam);
- if (n == 0) {
- fprintf(stderr, "No matching entry for %s in datafile.\n",
- m2);
- exit(1);
- } else if (n > 1) {
- fprintf(stderr, "More than one matching entry for %s in datafile.\n",
- m2);
- exit(1);
- }
- a2 = alpha(phi_a, lam_a, phi, lam);
-
- a_center = center_angle(a1, a2, dist_c_m1, dist_c_m2);
- mm_per_deg = (dist_c_m1 - dist_c_m2) / (tan(a1 - a_center) - tan(a2 - a_center));
- fprintf(stderr, "center_angle = %f\n", a_center / deg2rad);
-
- }
-
-
-
- read_file(phi_a, lam_a, data_file, dist, h_d_ratio, mm_per_deg, a_center);
-}
-
-
-#define PS_HEADER "%!\n \
-/mountain { % expects name, height, dist, deg on stack\n \
- newpath % Start a new path\n \
- gsave % Keep rotations temporary\n \
- rotate % Rotate by degrees on stack\n \
- dup \n \
- 0 moveto \n \
- add % End line at 20 + dist + height\n \
- dup % Duplicate end of line pos\n \
- 0 lineto\n \
- 10 add % Start text 10 from end of line \n \
- 0 moveto\n \
- show\n \
- stroke\n \
- grestore % Get back the unrotated state\n \
-} def \n \
-/Times-Roman findfont \n \
-8 scalefont \n \
-setfont \n \
-300 344 translate\n \
-5 0 moveto\n \
-0 0 5 0 360 arc\n \
-stroke
-%11 0 moveto\n"
-
-
-#define PS_HEADER_PHOTO "%!\n \
- /mountain { % expects name, height, dist, deg on stack\n \
- newpath\n \
- 0.352778 div % convert from mm to PS points\n \
- dup\n \
- dup\n \
- 200 exch moveto\n \
- 500 exch lineto\n \
- 500 exch moveto\n \
- pop \n \
- pop \n \
- show\n \
- stroke\n \
-} def\n \
-\n \
- /Times-Roman findfont \n \
- 6 scalefont \n \
- setfont \n
- 0 420 translate \n"
-
-
-int get_pos(char *filename, char *name, double *phi, double *lam) {
- FILE *fp;
- char buf[1024];
- char *vals[10];
- char **ap, *bp;
- double height;
- int found = 0;
-
- fp = fopen(filename, "r");
- if (!fp) {
- perror("fopen");
- return 1;
- }
-
- while (fgets(buf, sizeof(buf), fp)) {
- bp = buf;
- for (ap = vals; (*ap = strsep(&bp, ",")) != NULL;)
- if (++ap >= &vals[10])
- break;
-
- if (strstr(vals[1], name)) {
- *phi = atof(vals[3]) * deg2rad;
- *lam = atof(vals[4]) * deg2rad;
- height = atof(vals[5]);
- fprintf(stderr, "Found matching entry: %s (%fm)\n", vals[1], height);
- found++;
- }
- }
-
- fclose(fp);
- return found;
-}
-
-
-
-int read_file(double phi_a, double lam_a, char *name,
- int dist, double h_d_ratio, double mm_per_deg, double a_center) {
- FILE *fp;
- char buf[1024];
- char *vals[10];
- char **ap, *bp;
- double phi_b, lam_b, height, alph, dist_center;
-
- fp = fopen(name, "r");
- if (!fp) {
- perror("fopen");
- return 1;
- }
-
- printf("%s", PS_HEADER_PHOTO);
- printf("/mm_per_deg %f def\n", mm_per_deg);
- printf("/a_center %f def\n", a_center);
- while (fgets(buf, sizeof(buf), fp)) {
- bp = buf;
- for (ap = vals; (*ap = strsep(&bp, ",")) != NULL;)
- if (++ap >= &vals[10])
- break;
-
- phi_b = atof(vals[3]) * deg2rad;
- lam_b = atof(vals[4]) * deg2rad;
-
- if (phi_a == phi_b && lam_a == lam_b) {
- continue;
- }
-
- height = atof(vals[5]);
-
- if (distance(phi_a, lam_a, phi_b, lam_b) * 6368000.0 > dist) {
- continue;
- }
-
-
- if (height / (distance(phi_a, lam_a, phi_b, lam_b)* 6368000) < h_d_ratio) {
- continue;
- }
-
- alph = alpha(phi_a, lam_a, phi_b, lam_b);
- alph = alph - a_center;
- if (alph > pi / 2.0 || alph < - pi /2.0) {
- continue;
- }
-
- dist_center = tan(alph) * mm_per_deg;
-
- printf("(%s) %f %f %f mountain\n", vals[1],
- height / 20.0,
- distance(phi_a, lam_a, phi_b, lam_b) * 25000.0,
- dist_center);
-
-
-#if 0
- printf("%f %f %s %fm\n",
- alpha(phi_a, lam_a, phi_b, lam_b),
- distance(phi_a, lam_a, phi_b, lam_b) * 6368000,
- vals[1],
- height);
-#endif
- }
-
- printf("showpage\n");
-
- fclose(fp);
- return 0;
-}
-
-
-double
-distance(double phi_a, double lam_a, double phi_b, double lam_b) {
- return acos(sin(phi_a)*sin(phi_b) + cos(phi_a)*cos(phi_b)*cos(lam_a - lam_b));
-}
-
-double
-sin_alpha(double lam_a, double lam_b, double phi_b, double c) {
- return sin(lam_b - lam_a) * cos(phi_b) / sin(c);
-}
-
-double
-cos_alpha(double phi_a, double phi_b, double c) {
- return (sin(phi_b) - sin(phi_a) * cos(c)) / (cos(phi_a) * sin(c));
-}
-
-double alpha(double phi_a, double lam_a, double phi_b, double lam_b) {
- double dist, sin_alph, cos_alph, alph;
-
- dist = distance(phi_a, lam_a, phi_b, lam_b);
- sin_alph = sin_alpha(lam_a, lam_b, phi_b, dist);
- cos_alph = cos_alpha(phi_a, phi_b, dist);
-
- if (sin_alph > 0) {
- alph = acos(cos_alph);
- } else {
- alph = 2.0 * pi - acos(cos_alph);
- }
-
- return alph;
-}
-
-double
-center_angle(double alph_a, double alph_b, double d1, double d2) {
- double tan_a, tan_b;
-
- fprintf(stderr, "a=%f, b=%f d1=%f d2=%f\n",
- alph_a / deg2rad, alph_b / deg2rad, d1, d2);
-
- tan_a = tan(alph_a - alph_b);
- fprintf(stderr, "tan_a %f\n", tan_a);
- tan_b = (d2 - d1 + ((sqrt((d2*(d2 - (2.0*d1*(1.0 + (2.0 * tan_a * tan_a))))) + (d1*d1))))) / (2.0*d2*tan_a);
-
- fprintf(stderr, "tan_b=%f\n", tan_b);
- return alph_a + atan(tan_b);
-}
-
-void
-usage() {
- fprintf(stderr,
- "usage: panorama -a <phi>,<lambda>\n"
- " -b <phi>,<lambda>\n");
-}
generated by cgit v1.2.3 (git 2.25.1) at 2024-09-12 14:58:46 +0000