Hi Jim, re: time frame to transition to use of newer imagery? > I'm really not sure about how long that I would need to have the > old AREA formats or the amount of time that I will have to figure > things out. Fair enough! re: > If I could just convert the full disk files, it would be all that I > would need. OK. There is a bit of a wrinkle on this that is probably not obvious. I create the old GOES-13 images in the Unidata-Wisconsin datastream as remaps using the full resolution GRB images as the source. This insures that the final resolution/clarity of the old products is as high as is possible. Remapping the new full disk imagery that is being sent in the Unidata-Wisconsin channel would result in images that are not as clear/sharp as using the full resolution ones as the source. re: > I do like the full disk coverage (space and time). Can I assume that you are referring to the new (GOES-16/17) full disk images in the Unidata-Wisconsin datastream? Here is something that will hopefully whet your appetite for using the new imagery. In order to demonstrate, you should startup an interactive McIDAS-X session in which the image frames are reasonably large and numerous. For instance, I run my McIDAS-X sessions (with MCGUI) with 100 frames each of which is 950x1500 (linesXelements) in size. In order to run a session with this many larger frames, you might need to increase the amount of shared memory that is made available by your machine. This is described in Chapter 1 of the Unidata McIDAS-X User's Guide. After you have a "nice" session setup, do the following: 1) point at one of the ADDE servers we operate on behalf or the user community DATALOC ADD NPGOESR atm.ucar.edu 2) define a McIDAS string as follows: TE TRUCOLFD "SF !1;ERASE F;ABITRUCOL NPGOESR/FD.!2 !1 X 'LAT=0 75 MAG=-2' 3) use the McIDAS REPEAT command to repeatedly run the commands defined in this string. Here is one example: REPEAT TRUCOLFD 1 TO 50 BY 1 -49 TO 0 This will run the ABITRUCOL McBASI script to create true color images from channels 1, 2 and 3. Here is the relevant snippit from the HELP ABITRUCOL invocation: Remarks: ABITRUCOL is based on the CIMSS Natural True Color method: http://cimss.ssec.wisc.edu/goes/OCLOFactSheetPDFs/ABIQuickGuide_CIMSSRGB_v2.pdf ABITRUCOL is a McBASI script that executes McIDAS commands to copy bands 1, 2, and 3 from an ABI image, manipulate them into the same resolution, create an approximated green band with IMGOPER, and then use the RGBDISP command to display the ABI natural true-color RGB image. While creating the true-color RGB image, the script runs IMGCOPY, IMGREMAP and IMGOPER commands that overwrite Area files 4501-4540 (AREA4501, AREA4502, ..., AREA4540). To change the Area file numbers or make any other changes (e.g., change the default dataset or station), copy the file ABITRUCOL.MCB from ~mcidas/data/ to your $HOME/mcidas/data/ directory and edit the file. Explanation of the REPEAT invocation: Repeat works like a DO loop. It executes the command(s) defined in a string substituting the values passed on the command line. For instance: REPEAT TRUCOLFD 1 TO 50 BY 1 -49 TO 0 BY 1 - 1 TO 50 BY 1 This will cause the script to be executed 50 times with the first loop index (indicated by !1) varying from 1 to 50 in steps of 1. I defined the McIDAS string to use this value as the frame counter, so this will step from frame 1 to frame 50 in increments of a single frame. - -49 TO 0 BY 1 In each iteration, the second index (indicated by !2) will be varied from -49 to 0 in steps of 1 Here, the ADDE dataset position will vary from -49 to 0 in increments of 1. The values are passed to the ABITRUCOL McBASI script which will then do the work of bringing down all three channels for each dataset position that are needed to create the (pseudo) true color image. It then displays each composite image in the appropriate frame using the RGBDISP command. To get going, try the following single frame execution: REPEAT TRUCOLFD 1 TO 1 BY 1 -30 TO -30 BY 1 The choice of '-30' will depend on the time of day that you run thins invocation as it says get the 30th most recent images for each channel from the dataset. I just ran this single frame invocation at 20:35Z, and it produced the image that I have attached to this note. Give the above a try and let me know if/when you have any questions. Cheers, Tom -- **************************************************************************** Unidata User Support UCAR Unidata Program (303) 497-8642 P.O. Box 3000 address@hidden Boulder, CO 80307 ---------------------------------------------------------------------------- Unidata HomePage http://www.unidata.ucar.edu **************************************************************************** Ticket Details =================== Ticket ID: SET-184160 Department: Support IDD Priority: Normal Status: Closed =================== NOTE: All email exchanges with Unidata User Support are recorded in the Unidata inquiry tracking system and then made publicly available through the web. If you do not want to have your interactions made available in this way, you must let us know in each email you send to us.
Attachment:
GOES16_TRUCOL_20200327_1520.jpg
Description: JPEG image
