From 707c217494eaecdf160a3b46efe0d86271f9528b Mon Sep 17 00:00:00 2001
From: fmederos <fmederos@gmail.com>
Date: Tue, 18 Apr 2017 02:44:27 -0300
Subject: [PATCH] Incorporate handy zoom funcions to plotfits.c similar to
 calibrate.c

---
 plotfits.c | 76 +++++++++++++++++++++++++++++++++++++++++++++++++-----
 1 file changed, 69 insertions(+), 7 deletions(-)

diff --git a/plotfits.c b/plotfits.c
index 748091b..7a9b576 100644
--- a/plotfits.c
+++ b/plotfits.c
@@ -283,6 +283,7 @@ void get_site(int site_id)
 int main(int argc,char *argv[])
 {
   int i;
+  float xmin,xmax,ymin,ymax,zmin,zmax,width;
   float tr[]={-0.5,1.0,0.0,-0.5,0.0,1.0};
   float heat_l[] = {0.0, 0.2, 0.4, 0.6, 1.0};
   float heat_r[] = {0.0, 0.5, 1.0, 1.0, 1.0};
@@ -333,7 +334,7 @@ int main(int argc,char *argv[])
     img.ra0=15.0*sex2dec(sra);
     img.de0=sex2dec(sde);
 
-    // Hour angle
+    // Hour angle and parallactic angle
     h=gmst(img.mjd)+m.lng-img.ra0;
     q=atan2(sin(h*D2R),(tan(m.lat*D2R)*cos(img.de0*D2R)-sin(img.de0*D2R)*cos(h*D2R)))*R2D;
     printf("Hour angle: %.3f deg, parallactic angle: %.3f deg\n",h,q);
@@ -365,15 +366,33 @@ int main(int argc,char *argv[])
   cat=read_pixel_catalog(filename);
   ast=read_astrometric_catalog(starfile,mag,sx,sy,-q);
 
-  // Plot image
+  // Open PGPlot server
   cpgopen("/xs");
-  cpgwnad(0.0,img.naxis1,0.0,img.naxis2);
-  cpgsfs(2);
-  cpgctab (heat_l,heat_r,heat_g,heat_b,5,1.0,0.5);
+//  cpgwnad(0.0,img.naxis1,0.0,img.naxis2);
+//  cpgsfs(2);
+//  cpgctab (heat_l,heat_r,heat_g,heat_b,5,1.0,0.5);
+  cpgask(0);
+  cpgsch(0.8);
+
+  // Default limits
+  width=(img.naxis1>img.naxis2) ? img.naxis1 : img.naxis2;
+  xmin=0.5*(img.naxis1-width);
+  xmax=0.5*(img.naxis1+width);
+  ymin=0.5*(img.naxis2-width);
+  ymax=0.5*(img.naxis2+width);
+  zmin=img.zmin;
+  zmax=img.zmax;
 
   // For ever loop
   for (;;) {      
     if (redraw==1) {
+      cpgeras();
+
+      cpgsvp(0.1,0.95,0.1,0.95);
+      cpgwnad(xmin,xmax,ymin,ymax);
+      cpglab("x (pix)","y (pix)"," ");
+      cpgsfs(2);
+      cpgctab (heat_l,heat_r,heat_g,heat_b,5,1.0,0.5);
       cpgimag(img.z,img.naxis1,img.naxis2,1,img.naxis1,1,img.naxis2,img.zmin,img.zmax,tr);
       cpgbox("BCTSNI",0.,0,"BCTSNI",0.,0);
     
@@ -415,7 +434,7 @@ int main(int argc,char *argv[])
     }
 
     // Reread
-    if (c=='r') {
+    if (c=='R') {
       ast=reread_astrometric_catalog(starfile,mag+1);
       redraw=1;
     }
@@ -437,7 +456,7 @@ int main(int argc,char *argv[])
       nselect++;
     }
     
-    // 
+    // Plot identified stars
     if (c=='p') {
       if (plotstars==1)
 	plotstars=0;
@@ -452,6 +471,49 @@ int main(int argc,char *argv[])
       redraw=1;
     }
 
+    // Center
+    if (c=='c') {
+      xmin=x-0.5*width;
+      xmax=x+0.5*width;
+      ymin=y-0.5*width;
+      ymax=y+0.5*width;
+      redraw=1;
+      continue;
+    }
+
+    // Reset
+    if (c=='r') {
+      width=(img.naxis1>img.naxis2) ? img.naxis1 : img.naxis2;
+      xmin=0.5*(img.naxis1-width);
+      xmax=0.5*(img.naxis1+width);
+      ymin=0.5*(img.naxis2-width);
+      ymax=0.5*(img.naxis2+width);
+      redraw=1;
+      continue;
+    }
+
+    // Zoom
+    if (c=='z' || c=='+' || c=='=') {
+      width/=1.25;
+      xmin=x-0.5*width;
+      xmax=x+0.5*width;
+      ymin=y-0.5*width;
+      ymax=y+0.5*width;
+      redraw=1;
+      continue;
+    }
+
+    // Unzoom
+    if (c=='x' || c=='-') {
+      width*=1.25;
+      xmin=x-0.5*width;
+      xmax=x+0.5*width;
+      ymin=y-0.5*width;
+      ymax=y+0.5*width;
+      redraw=1;
+      continue;
+    }
+
     // Help
     if (c=='h') {
       printf("Calibrates astrometry. Initially requires manual matching of at least three stars. Use 'a' to select star on the image, then 'b' to select star from the catalog, then 'f' to fit");