ham-ref/gpredict
1
0
Fork 0
mirror of https://github.com/csete/gpredict synced 2026-08-13 17:39:46 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions

#3 finish solar terminator: work with maps at longitudes other than 0

Browse source
This commit is contained in:
Joe 2020-10-25 12:53:59 +01:00 • committed by Alexandru Csete
parent 7ab700e50a
commit e8aefd2c2e
1 changed files with 14 additions and 10 deletions
Download patch file Download diff file Expand all files Collapse all files

24
src/gtk-sat-map.c
View file

@ -2691,9 +2691,9 @@ static void redraw_terminator(GtkSatMap * satmap)
for (longitude = -180; longitude <= 180; ++longitude)
{
lx = cos(de2ra * (longitude + sgn(sz) * 90));
ly = sin(de2ra * (longitude + sgn(sz) * 90));
int centered_longitude = longitude + (satmap->left_side_lon - 180.0);
lx = cos(de2ra * (centered_longitude + sgn(sz) * 90));
ly = sin(de2ra * (centered_longitude + sgn(sz) * 90));
/* lz = 0.0; */
rx = ly * sz /* -lz*sy */ ;
@ -2703,19 +2703,23 @@ static void redraw_terminator(GtkSatMap * satmap)
gdouble length = sqrt(rx * rx + ry * ry + rz * rz);
lonlat_to_xy(satmap,
longitude, asin(rz / length) * (1.0 / de2ra), &x, &y);
centered_longitude, asin(rz / length) * (1.0 / de2ra), &x, &y);
if( 180 == longitude )
{
// make sure the last point is on the right side
x = satmap->x0 + satmap->width;
}
line->coords[2 * (longitude + 181)] = x;
line->coords[2 * (longitude + 181) + 1] = y;
}
lonlat_to_xy(satmap, -180.0, sz < 0.0 ? 90.0 : -90.0, &x, &y);
line->coords[0] = x;
line->coords[1] = y;
line->coords[0] = satmap->x0;
line->coords[1] = sz < 0.0 ? satmap->y0 : ( satmap->y0 + satmap->height );
lonlat_to_xy(satmap, 180.0, sz < 0.0 ? 90.0 : -90.0, &x, &y);
line->coords[724] = x;
line->coords[725] = y;
line->coords[724] = satmap->x0 + satmap->width;
line->coords[725] = sz < 0.0 ? satmap->y0 : ( satmap->y0 + satmap->height );
g_object_set(satmap->terminator, "points", line, NULL);
goo_canvas_points_unref(line);