new mask

2 files changed, 396 insertions, 390 deletions

diff --git a/tksao/frame/framergb.C b/tksao/frame/framergb.C
index 12583ef..b1c3c20 100644
--- a/tksao/frame/framergb.C
+++ b/tksao/frame/framergb.C
@@ -281,396 +281,6 @@ BBox FrameRGB::imageBBox(FrScale::SecMode mode)
   return rr;
 }
 
-void FrameRGB::loadRGBCube(MemType which, const char* fn, FitsImage* img)
-{
-  if (!img || !img->isValid() || !(img->isImage() || img->isPost()) || (img->depth() != 3))
-    goto error;
-
-  context[0].bfits_ = img;
-
-  if (img->isPost())
-    which = POST;
-
-  switch (which) {
-  case ALLOC:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImageFitsNextAlloc(&context[1], interp,
-				   fn, context[0].bfits_->fitsFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImageFitsNextAlloc(&context[2], interp,
-				   fn, context[1].bfits_->fitsFile(),3);
-    break;
-  case ALLOCGZ:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImageFitsNextAllocGZ(&context[1], interp,
-				     fn, context[0].bfits_->fitsFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImageFitsNextAllocGZ(&context[2], interp,
-				     fn, context[1].bfits_->fitsFile(),3);
-    break;
-  case CHANNEL:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImageFitsNextChannel(&context[1], interp,
-				     fn, context[0].bfits_->fitsFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImageFitsNextChannel(&context[2], interp,
-				     fn, context[1].bfits_->fitsFile(),3);
-    break;
-  case MMAP:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImageFitsNextMMap(&context[1], interp,
-				  fn, context[0].bfits_->fitsFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImageFitsNextMMap(&context[2], interp,
-				  fn, context[1].bfits_->fitsFile(),3);
-    break;
-  case SMMAP:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImageFitsNextSMMap(&context[1], interp,
-				   fn, context[0].bfits_->fitsFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImageFitsNextSMMap(&context[2], interp,
-				   fn, context[1].bfits_->fitsFile(),3);
-    break;
-  case MMAPINCR:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImageFitsNextMMapIncr(&context[1], interp,
-				      fn, context[0].bfits_->fitsFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImageFitsNextMMapIncr(&context[2], interp,
-				      fn, context[1].bfits_->fitsFile(),3);
-    break;
-  case SHARE:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImageFitsNextShare(&context[1], interp,
-				   fn, context[0].bfits_->fitsFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImageFitsNextShare(&context[2], interp,
-				   fn, context[1].bfits_->fitsFile(),3);
-    break;
-  case SSHARE:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImageFitsNextSShare(&context[1], interp,
-				    fn, context[0].bfits_->fitsFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImageFitsNextSShare(&context[2], interp,
-				    fn, context[1].bfits_->fitsFile(),3);
-    break;
-  case SOCKET:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImageFitsNextSocket(&context[1], interp,
-				    fn, context[0].bfits_->fitsFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImageFitsNextSocket(&context[2], interp,
-				    fn, context[1].bfits_->fitsFile(),3);
-    break;
-  case SOCKETGZ:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImageFitsNextSocketGZ(&context[1], interp,
-				      fn, context[0].bfits_->fitsFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImageFitsNextSocketGZ(&context[2], interp,
-				      fn, context[1].bfits_->fitsFile(),3);
-    break;
-  case VAR:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImageFitsNextVar(&context[1], interp,
-				 fn, context[0].bfits_->fitsFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImageFitsNextVar(&context[2], interp,
-				 fn, context[1].bfits_->fitsFile(),3);
-    break;
-  case POST:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImageFitsNextPost(&context[1], interp,
-				  img, context[0].bfits_->baseFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImageFitsNextPost(&context[2], interp,
-				  img, context[1].bfits_->baseFile(),3);
-    break;
-  case PHOTO:
-    if (context[0].bfits_ && context[0].bfits_->isValid())
-      context[1].bfits_ = 
-	new FitsImagePhotoCubeNext(&context[1], interp,
-				   fn, context[0].bfits_->baseFile(),2);
-    if (context[1].bfits_ && context[1].bfits_->isValid())
-      context[2].bfits_ = 
-	new FitsImagePhotoCubeNext(&context[2], interp,
-				   fn, context[1].bfits_->baseFile(),3);
-    break;
-  default:
-    // na
-    break;
-  }
-
-  // is everything ok?
-  if (context[0].bfits_ && context[0].bfits_->isValid() &&
-      (context[0].bfits_->isImage() || context[0].bfits_->isPost()) &&
-      context[1].bfits_ && context[1].bfits_->isValid() &&
-      (context[1].bfits_->isImage() || context[1].bfits_->isPost()) &&
-      context[2].bfits_ && context[2].bfits_->isValid() &&
-      (context[2].bfits_->isImage() || context[2].bfits_->isPost())) {
-
-    loadRGBFinish();
-    return;
-  }
-
- error:
-  context[0].unload();
-  context[1].unload();
-  context[2].unload();
-
-  reset();
-  updateColorScale();
-    
-  Tcl_AppendResult(interp, "Unable to load rgb cube file", NULL);
-  result = TCL_ERROR;
-  return;
-}
-
-void FrameRGB::loadRGBImage(MemType which, const char* fn, FitsImage* img)
-{
-  FitsImage* r = img;
-  FitsImage* g = NULL;
-  FitsImage* b = NULL;
-
-  if (!img || !img->isValid() || !(img->isImage() || img->isPost()))
-    goto error;
-
-  switch (which) {
-  case ALLOC:
-    if (r && r->isValid())
-      g = new FitsImageMosaicNextAlloc(&context[1], interp, 
-				       fn, r->fitsFile(), 
-				       FitsFile::NOFLUSH,1);
-    if (g && g->isValid())
-      b = new FitsImageMosaicNextAlloc(&context[2], interp, 
-				       fn, g->fitsFile(), 
-				       FitsFile::NOFLUSH,1);
-    break;
-  case ALLOCGZ:
-    if (r && r->isValid())
-      g = new FitsImageMosaicNextAllocGZ(&context[1], interp, 
-					 fn, r->fitsFile(), 
-					 FitsFile::NOFLUSH,1);
-    if (g && g->isValid())
-      b = new FitsImageMosaicNextAllocGZ(&context[2], interp, 
-					 fn, g->fitsFile(), 
-					 FitsFile::NOFLUSH,1);
-    break;
-  case CHANNEL:
-    if (r && r->isValid())
-      g = new FitsImageMosaicNextChannel(&context[1], interp, 
-					 fn, r->fitsFile(), 
-					 FitsFile::NOFLUSH,1);
-    if (g && g->isValid())
-      b = new FitsImageMosaicNextChannel(&context[2], interp, 
-					 fn, g->fitsFile(), 
-					 FitsFile::NOFLUSH,1);
-    break;
-  case MMAP:
-    if (r && r->isValid())
-      g = new FitsImageMosaicNextMMap(&context[1], interp, 
-				      fn, r->fitsFile(),1);
-    if (g && g->isValid())
-      b = new FitsImageMosaicNextMMap(&context[2], interp, 
-				      fn, g->fitsFile(),1);
-    break;
-  case MMAPINCR:
-    if (r && r->isValid())
-      g = new FitsImageMosaicNextMMapIncr(&context[1], interp, 
-					  fn, r->fitsFile(),1);
-    if (g && g->isValid())
-      b = new FitsImageMosaicNextMMapIncr(&context[2], interp,
-					  fn, g->fitsFile(),1);
-    break;
-  case SHARE:
-    if (r && r->isValid())
-      g = new FitsImageMosaicNextShare(&context[1], interp, 
-				       fn, r->fitsFile(),1);
-    if (g && g->isValid())
-      b = new FitsImageMosaicNextShare(&context[2], interp, 
-				       fn, g->fitsFile(),1);
-    break;
-  case SOCKET:
-    if (r && r->isValid())
-      g = new FitsImageMosaicNextSocket(&context[1], interp, 
-					fn, r->fitsFile(), 
-					FitsFile::FLUSH,1);
-    if (g && g->isValid())
-      b = new FitsImageMosaicNextSocket(&context[2], interp,
-					fn,g->fitsFile(), 
-					FitsFile::FLUSH,1);
-    break;
-  case SOCKETGZ:
-    if (r && r->isValid())
-      g = new FitsImageMosaicNextSocketGZ(&context[1], interp, 
-					  fn, r->fitsFile(), 
-					  FitsFile::FLUSH,1);
-    if (g && g->isValid())
-      b = new FitsImageMosaicNextSocketGZ(&context[2], interp,
-					  fn, g->fitsFile(), 
-					  FitsFile::FLUSH,1);
-    break;
-  case VAR:
-    if (r && r->isValid())
-      g = new FitsImageMosaicNextVar(&context[1], interp, 
-				     fn, r->fitsFile(),1);
-    if (g && g->isValid())
-      b = new FitsImageMosaicNextVar(&context[2], interp, 
-				     fn, g->fitsFile(),1);
-    break;
-  default:
-    // na
-    break;
-  }
-
-  // ok, figure out which is which channel
-  context[0].bfits_ = context[1].bfits_ = context[2].bfits_ = NULL;
-
-  {
-    const char* ext = r->fitsFile()->extname();
-    if (ext) {
-      if (!strncmp(ext,"RED",3))
-	context[0].bfits_ = r;
-      else if (!strncmp(ext,"GREEN",3)) {
-	context[1].bfits_ = r;
-	r->setContext(&context[1]);
-      }
-      else if (!strncmp(ext,"BLUE",3)) {
-	context[2].bfits_ = r;
-	r->setContext(&context[2]);
-      }
-      else
-	context[0].bfits_ = r;
-    }
-    else
-      context[0].bfits_ = r;
-  }
-
-  {
-    const char* ext = g->fitsFile()->extname();
-    if (ext) {
-      if (!strncmp(ext,"RED",3)) {
-	context[0].bfits_ = g;
-	g->setContext(&context[0]);
-      }
-      else if (!strncmp(ext,"GREEN",3))
-	context[1].bfits_ = g;
-      else if (!strncmp(ext,"BLUE",3)) {
-	context[2].bfits_ = g;
-	g->setContext(&context[3]);
-      }
-      else
-	context[1].bfits_ = g;
-    }
-    else
-      context[1].bfits_ = g;
-  }
-
-  {
-    const char* ext = b->fitsFile()->extname();
-    if (ext) {
-      if (!strncmp(ext,"RED",3)) {
-	context[0].bfits_ = b;
-	b->setContext(&context[0]);
-      }
-      else if (!strncmp(ext,"GREEN",3)) {
-	context[1].bfits_ = b;
-	b->setContext(&context[1]);
-      }
-      else if (!strncmp(ext,"BLUE",3))
-	context[2].bfits_ = b;
-      else
-	context[2].bfits_ = b;
-    }
-    else
-      context[2].bfits_ = b;
-  }
-
-  // is everything ok?
-  if (context[0].bfits_ && context[0].bfits_->isValid() && 
-      (context[0].bfits_->isImage() || context[0].bfits_->isPost()) &&
-      context[1].bfits_ && context[1].bfits_->isValid() && 
-      (context[1].bfits_->isImage() || context[1].bfits_->isPost()) &&
-      context[2].bfits_ && context[2].bfits_->isValid() && 
-      (context[2].bfits_->isImage() || context[2].bfits_->isPost())) {
-
-    loadRGBFinish();
-    return;
-  }
-
- error:
-  context[0].unload();
-  context[1].unload();
-  context[2].unload();
-
-  reset();
-  updateColorScale();
-    
-  Tcl_AppendResult(interp, "Unable to load rgb image file", NULL);
-  result = TCL_ERROR;
-  return;
-}
-
-void FrameRGB::loadRGBFinish()
-{
-  for (int ii=0; ii<3; ii++) {
-    context[ii].loadInit(1,NOMOSAIC,Coord::WCS);
-    context[ii].loadFinish();
-  }
-
-  channel = 0;
-  currentContext = &context[channel];
-  keyContext = &context[channel];
-  keyContextSet =1;
-
-  alignWCS();
-  if (!preservePan) {
-    centerImage();
-    // cursor is in REF, crosshair in REF
-    crosshair = cursor;
-  }
-  updateColorScale();
-  update(MATRIX);
-}
-
-void FrameRGB::loadDone(int rr)
-{
-  if (rr) {
-    if (!keyContextSet) {
-      keyContext = currentContext;
-      keyContextSet =1;
-    }
-
-    Base::loadDone(rr);
-  }
-}
-
 void FrameRGB::pushMatrices()
 {
   for (int ii=0; ii<3; ii++)
diff --git a/tksao/frame/framergbload.C b/tksao/frame/framergbload.C
new file mode 100644
index 0000000..98fbfae
--- /dev/null
+++ b/tksao/frame/framergbload.C
@@ -0,0 +1,396 @@
+// Copyright (C) 1999-2018
+// Smithsonian Astrophysical Observatory, Cambridge, MA, USA
+// For conditions of distribution and use, see copyright notice in "copyright"
+
+#include "framergb.h"
+
+void FrameRGB::loadRGBCube(MemType which, const char* fn, FitsImage* img)
+{
+  if (!img || !img->isValid() || !(img->isImage() || img->isPost()) || (img->depth() != 3))
+    goto error;
+
+  context[0].bfits_ = img;
+
+  if (img->isPost())
+    which = POST;
+
+  switch (which) {
+  case ALLOC:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImageFitsNextAlloc(&context[1], interp,
+				   fn, context[0].bfits_->fitsFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImageFitsNextAlloc(&context[2], interp,
+				   fn, context[1].bfits_->fitsFile(),3);
+    break;
+  case ALLOCGZ:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImageFitsNextAllocGZ(&context[1], interp,
+				     fn, context[0].bfits_->fitsFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImageFitsNextAllocGZ(&context[2], interp,
+				     fn, context[1].bfits_->fitsFile(),3);
+    break;
+  case CHANNEL:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImageFitsNextChannel(&context[1], interp,
+				     fn, context[0].bfits_->fitsFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImageFitsNextChannel(&context[2], interp,
+				     fn, context[1].bfits_->fitsFile(),3);
+    break;
+  case MMAP:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImageFitsNextMMap(&context[1], interp,
+				  fn, context[0].bfits_->fitsFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImageFitsNextMMap(&context[2], interp,
+				  fn, context[1].bfits_->fitsFile(),3);
+    break;
+  case SMMAP:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImageFitsNextSMMap(&context[1], interp,
+				   fn, context[0].bfits_->fitsFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImageFitsNextSMMap(&context[2], interp,
+				   fn, context[1].bfits_->fitsFile(),3);
+    break;
+  case MMAPINCR:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImageFitsNextMMapIncr(&context[1], interp,
+				      fn, context[0].bfits_->fitsFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImageFitsNextMMapIncr(&context[2], interp,
+				      fn, context[1].bfits_->fitsFile(),3);
+    break;
+  case SHARE:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImageFitsNextShare(&context[1], interp,
+				   fn, context[0].bfits_->fitsFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImageFitsNextShare(&context[2], interp,
+				   fn, context[1].bfits_->fitsFile(),3);
+    break;
+  case SSHARE:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImageFitsNextSShare(&context[1], interp,
+				    fn, context[0].bfits_->fitsFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImageFitsNextSShare(&context[2], interp,
+				    fn, context[1].bfits_->fitsFile(),3);
+    break;
+  case SOCKET:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImageFitsNextSocket(&context[1], interp,
+				    fn, context[0].bfits_->fitsFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImageFitsNextSocket(&context[2], interp,
+				    fn, context[1].bfits_->fitsFile(),3);
+    break;
+  case SOCKETGZ:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImageFitsNextSocketGZ(&context[1], interp,
+				      fn, context[0].bfits_->fitsFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImageFitsNextSocketGZ(&context[2], interp,
+				      fn, context[1].bfits_->fitsFile(),3);
+    break;
+  case VAR:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImageFitsNextVar(&context[1], interp,
+				 fn, context[0].bfits_->fitsFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImageFitsNextVar(&context[2], interp,
+				 fn, context[1].bfits_->fitsFile(),3);
+    break;
+  case POST:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImageFitsNextPost(&context[1], interp,
+				  img, context[0].bfits_->baseFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImageFitsNextPost(&context[2], interp,
+				  img, context[1].bfits_->baseFile(),3);
+    break;
+  case PHOTO:
+    if (context[0].bfits_ && context[0].bfits_->isValid())
+      context[1].bfits_ = 
+	new FitsImagePhotoCubeNext(&context[1], interp,
+				   fn, context[0].bfits_->baseFile(),2);
+    if (context[1].bfits_ && context[1].bfits_->isValid())
+      context[2].bfits_ = 
+	new FitsImagePhotoCubeNext(&context[2], interp,
+				   fn, context[1].bfits_->baseFile(),3);
+    break;
+  default:
+    // na
+    break;
+  }
+
+  // is everything ok?
+  if (context[0].bfits_ && context[0].bfits_->isValid() &&
+      (context[0].bfits_->isImage() || context[0].bfits_->isPost()) &&
+      context[1].bfits_ && context[1].bfits_->isValid() &&
+      (context[1].bfits_->isImage() || context[1].bfits_->isPost()) &&
+      context[2].bfits_ && context[2].bfits_->isValid() &&
+      (context[2].bfits_->isImage() || context[2].bfits_->isPost())) {
+
+    loadRGBFinish();
+    return;
+  }
+
+ error:
+  context[0].unload();
+  context[1].unload();
+  context[2].unload();
+
+  reset();
+  updateColorScale();
+    
+  Tcl_AppendResult(interp, "Unable to load rgb cube file", NULL);
+  result = TCL_ERROR;
+  return;
+}
+
+void FrameRGB::loadRGBImage(MemType which, const char* fn, FitsImage* img)
+{
+  FitsImage* r = img;
+  FitsImage* g = NULL;
+  FitsImage* b = NULL;
+
+  if (!img || !img->isValid() || !(img->isImage() || img->isPost()))
+    goto error;
+
+  switch (which) {
+  case ALLOC:
+    if (r && r->isValid())
+      g = new FitsImageMosaicNextAlloc(&context[1], interp, 
+				       fn, r->fitsFile(), 
+				       FitsFile::NOFLUSH,1);
+    if (g && g->isValid())
+      b = new FitsImageMosaicNextAlloc(&context[2], interp, 
+				       fn, g->fitsFile(), 
+				       FitsFile::NOFLUSH,1);
+    break;
+  case ALLOCGZ:
+    if (r && r->isValid())
+      g = new FitsImageMosaicNextAllocGZ(&context[1], interp, 
+					 fn, r->fitsFile(), 
+					 FitsFile::NOFLUSH,1);
+    if (g && g->isValid())
+      b = new FitsImageMosaicNextAllocGZ(&context[2], interp, 
+					 fn, g->fitsFile(), 
+					 FitsFile::NOFLUSH,1);
+    break;
+  case CHANNEL:
+    if (r && r->isValid())
+      g = new FitsImageMosaicNextChannel(&context[1], interp, 
+					 fn, r->fitsFile(), 
+					 FitsFile::NOFLUSH,1);
+    if (g && g->isValid())
+      b = new FitsImageMosaicNextChannel(&context[2], interp, 
+					 fn, g->fitsFile(), 
+					 FitsFile::NOFLUSH,1);
+    break;
+  case MMAP:
+    if (r && r->isValid())
+      g = new FitsImageMosaicNextMMap(&context[1], interp, 
+				      fn, r->fitsFile(),1);
+    if (g && g->isValid())
+      b = new FitsImageMosaicNextMMap(&context[2], interp, 
+				      fn, g->fitsFile(),1);
+    break;
+  case MMAPINCR:
+    if (r && r->isValid())
+      g = new FitsImageMosaicNextMMapIncr(&context[1], interp, 
+					  fn, r->fitsFile(),1);
+    if (g && g->isValid())
+      b = new FitsImageMosaicNextMMapIncr(&context[2], interp,
+					  fn, g->fitsFile(),1);
+    break;
+  case SHARE:
+    if (r && r->isValid())
+      g = new FitsImageMosaicNextShare(&context[1], interp, 
+				       fn, r->fitsFile(),1);
+    if (g && g->isValid())
+      b = new FitsImageMosaicNextShare(&context[2], interp, 
+				       fn, g->fitsFile(),1);
+    break;
+  case SOCKET:
+    if (r && r->isValid())
+      g = new FitsImageMosaicNextSocket(&context[1], interp, 
+					fn, r->fitsFile(), 
+					FitsFile::FLUSH,1);
+    if (g && g->isValid())
+      b = new FitsImageMosaicNextSocket(&context[2], interp,
+					fn,g->fitsFile(), 
+					FitsFile::FLUSH,1);
+    break;
+  case SOCKETGZ:
+    if (r && r->isValid())
+      g = new FitsImageMosaicNextSocketGZ(&context[1], interp, 
+					  fn, r->fitsFile(), 
+					  FitsFile::FLUSH,1);
+    if (g && g->isValid())
+      b = new FitsImageMosaicNextSocketGZ(&context[2], interp,
+					  fn, g->fitsFile(), 
+					  FitsFile::FLUSH,1);
+    break;
+  case VAR:
+    if (r && r->isValid())
+      g = new FitsImageMosaicNextVar(&context[1], interp, 
+				     fn, r->fitsFile(),1);
+    if (g && g->isValid())
+      b = new FitsImageMosaicNextVar(&context[2], interp, 
+				     fn, g->fitsFile(),1);
+    break;
+  default:
+    // na
+    break;
+  }
+
+  // ok, figure out which is which channel
+  context[0].bfits_ = context[1].bfits_ = context[2].bfits_ = NULL;
+
+  {
+    const char* ext = r->fitsFile()->extname();
+    if (ext) {
+      if (!strncmp(ext,"RED",3))
+	context[0].bfits_ = r;
+      else if (!strncmp(ext,"GREEN",3)) {
+	context[1].bfits_ = r;
+	r->setContext(&context[1]);
+      }
+      else if (!strncmp(ext,"BLUE",3)) {
+	context[2].bfits_ = r;
+	r->setContext(&context[2]);
+      }
+      else
+	context[0].bfits_ = r;
+    }
+    else
+      context[0].bfits_ = r;
+  }
+
+  {
+    const char* ext = g->fitsFile()->extname();
+    if (ext) {
+      if (!strncmp(ext,"RED",3)) {
+	context[0].bfits_ = g;
+	g->setContext(&context[0]);
+      }
+      else if (!strncmp(ext,"GREEN",3))
+	context[1].bfits_ = g;
+      else if (!strncmp(ext,"BLUE",3)) {
+	context[2].bfits_ = g;
+	g->setContext(&context[3]);
+      }
+      else
+	context[1].bfits_ = g;
+    }
+    else
+      context[1].bfits_ = g;
+  }
+
+  {
+    const char* ext = b->fitsFile()->extname();
+    if (ext) {
+      if (!strncmp(ext,"RED",3)) {
+	context[0].bfits_ = b;
+	b->setContext(&context[0]);
+      }
+      else if (!strncmp(ext,"GREEN",3)) {
+	context[1].bfits_ = b;
+	b->setContext(&context[1]);
+      }
+      else if (!strncmp(ext,"BLUE",3))
+	context[2].bfits_ = b;
+      else
+	context[2].bfits_ = b;
+    }
+    else
+      context[2].bfits_ = b;
+  }
+
+  // is everything ok?
+  if (context[0].bfits_ && context[0].bfits_->isValid() && 
+      (context[0].bfits_->isImage() || context[0].bfits_->isPost()) &&
+      context[1].bfits_ && context[1].bfits_->isValid() && 
+      (context[1].bfits_->isImage() || context[1].bfits_->isPost()) &&
+      context[2].bfits_ && context[2].bfits_->isValid() && 
+      (context[2].bfits_->isImage() || context[2].bfits_->isPost())) {
+
+    loadRGBFinish();
+    return;
+  }
+
+ error:
+  context[0].unload();
+  context[1].unload();
+  context[2].unload();
+
+  reset();
+  updateColorScale();
+    
+  Tcl_AppendResult(interp, "Unable to load rgb image file", NULL);
+  result = TCL_ERROR;
+  return;
+}
+
+void FrameRGB::loadRGBFinish()
+{
+  for (int ii=0; ii<3; ii++) {
+    context[ii].loadInit(1,NOMOSAIC,Coord::WCS);
+    context[ii].loadFinish();
+  }
+
+  channel = 0;
+  currentContext = &context[channel];
+  keyContext = &context[channel];
+  keyContextSet =1;
+
+  alignWCS();
+  if (!preservePan) {
+    centerImage();
+    // cursor is in REF, crosshair in REF
+    crosshair = cursor;
+  }
+  updateColorScale();
+  update(MATRIX);
+}
+
+void FrameRGB::loadDone(int rr)
+{
+  if (rr) {
+    if (!keyContextSet) {
+      keyContext = currentContext;
+      keyContextSet =1;
+    }
+
+    Base::loadDone(rr);
+  }
+}
+