blob: a07c97e367ec1c1b95d2b93f7950114b01219bcf (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
|
#include <u.h>
#include <libc.h>
#include "map.h"
int
Xpolyconic(struct place *place, double *x, double *y)
{
double r, alpha;
double lat2, lon2;
if(fabs(place->nlat.l) > .01) {
r = place->nlat.c / place->nlat.s;
alpha = place->wlon.l * place->nlat.s;
*y = place->nlat.l + r*(1 - cos(alpha));
*x = - r*sin(alpha);
} else {
lon2 = place->wlon.l * place->wlon.l;
lat2 = place->nlat.l * place->nlat.l;
*y = place->nlat.l * (1+(lon2/2)*(1-(8+lon2)*lat2/12));
*x = - place->wlon.l * (1-lat2*(3+lon2)/6);
}
return(1);
}
proj
polyconic(void)
{
return(Xpolyconic);
}
|