Cloudy SolarSoftware

From RHESSI Wiki

Revision as of 20:11, 22 June 2010

NOTE

This nugget is unfinished and still undergoing changes.

Introduction

Figure 1. This is an overview of the complete workflow from searching and finding data with the VSO, pre-processing the data when necessary using the PrepServer, to visualizing the data with PLOTMAN. Points 1-2-3 in the figure correspond to steps 1-2-3 in the text.

In our project "Extending the Virtual Solar Observatory ([1])," we have combined features available in Solar Software (SSW) to produce an integrated environment, supporting data location, retrieval, preparation, and analysis. This workflow in shown Figure 1 with examples given below. Our goal is an integrated analysis experience in IDL, easy-to-use but flexible enough to allow more sophisticated procedures such as multi-instrument analysis. To that end, we have made the transition from a locally oriented setting where all the analysis is done on the user’s computer, to an extended analysis environment where IDL has access to services available on on other computers through the Internet. We have implemented a form of Cloud Computing that uses the VSO search and a new data retrieval and pre-processing server (PrepServer) that provides remote execution of instrument-specific data preparation. The raw and pre-processed data can be displayed with our GUI plotting suite, PLOTMAN, which can display different data types and perform basic data operations. Figure 1 shows the entire framework with numbers corresponding to steps 1-2-3 explained below in the text.

Our environment supports data from a growing number of solar instruments that currently includes RHESSI, SOHO/EIT, TRACE, SECCHI/EUVI, HINODE/XRT, and HINODE/EIS.

Minimum Requirements

• IDL 6.4
• Sun Java 1.5
• SolarSoftware (SSW) with GEN package (standard)

Step 1: Finding the Data

SHOW_SYNOP shown in Figure 2 is an IDL graphical user interface (GUI) to search for and retrieve instrument data within a specified time interval using the VSO or other data finding facilities.

Figure 2. This is a screen shot of the SHOW_SYNOP GUI used to find and retrieve specific data sets within a given time interval. The three colored boxes indicate: (red) the VSO search form (start/end time and instrument), (green) the search result list, and (blue) the downloaded level-0 and pre-processed data. The yellow button is to display the configuration menu.

Search results can be directly downloaded into the active SSW IDL environment. To start searching with Show_SYNOP, first open the GUI by typing SHOW_SYNOP in your SSW IDL environment command line. The red box highlights the VSO search form with "Start Time" and "End Time" specifying the search interval and "remote sites ->" defining the instrument (TRACE in this example). Click on the "Search" button to query the VSO for data files that will be displayed in the list below the search form shown in the green box.

Alternatively, the search can be initiated directly from the IDL command line using the procedure VSO_FILES. An example script found in the documentation (link) can be used to find, retrieve, and plot near-coincident images from SOHO/EIT and Hinode/XRT. The resulting images are shown in Figure ?.

Step 2: Pre-processing

Figure 3. This figure shows PLOTMAN displaying an XRT image (January 25 2007 06:56:56) that has been pre-processed with default settings by the PrepServer.

Figure 4. When pre-processing a TRACE data file that contains multiple TRACE images, a selection widget will open allowing picking a set of images.

Figure 5. This figure shows PLOTMAN displaying a TRACE image (171 A, January 25 2007, 06:15:09) that has been pre-processed with default settings by the PrepServer. If multiple images were selected in the selection window (see Figure 4), they can be accessed by clicking on "Window_Control" in PLOTMAN's menu panel at the top.

Figure 6. This figure shows PLOTMAN displaying an contour overlay of a reconstructed RHESSI clean 6-12keV image on an EIT image (171 A, January 25 2007, 07:00:13). The reconstruction of the RHESSI image as well as the pre-processing of the EIT images was done with the PrepServer.

Typically, instrument data found with the VSO are unprocessed level-0 data. The PrepServer offers remote pre-processing of those data from within SHOW_SYNOP or from the IDL command line before or after downloading. In SHOW_SYNOP, search results are displayed in the list underneath the VSO search form (see Figure 2, green box). There are two different pre-processing behaviors configurable:

• By default a click on "Download" retrieves the level-0 data file, stores it on the local hard drive and adds it to a list (Figure 2, blue box). The level-0 data file is pre-processed by selecting it in the blue box and clicking on "Display". This uploads the level-0 data file to the PrepServer where it is pre-processed, and downloads the pre-processed product to the local hard drive.
• In auto-preprocessing mode, which can be activated in the configuration menu (Figure 2, yellow box), SHOW_SYNOP sends the URL of the level-0 data file to the PrepServer and only downloads the pre-processed data.

Please note: SHOW_SYNOP and the PrepServer are dependent on data providers for downloading speed. Be aware that downloading data may take a few minutes, during which SHOW_SYNOP or your SSW IDL session may be unresponsive.

VSO_PREP allows for data pre-processing from an SSW IDL command line:

Command
vso_prep, infile, outfile=outfile, instrument=instrument, oprep=oprep, _extra=extra

Input
infile: A local or remote file to be pre-processed
instrument: (optional) The instrument name. Default: Determined from infile
_extra: (optional) This parameter will collect all additional input parameters

Output
outfile: (optional) A file name for the pre-processed file. Default: prep_infile
oprep: (optional) An object that will be created after the pre-processing. It will allow for direct access to the data (i.e. for display)

; This searches the VSO for an EIT file and sends the returned link to the PrepServer for pre-processing
urls = vso_files('25-Jan-2007 06:55', inst='eit')
vso_prep, urls, out='eit.fits'
% VSO_PREP_MAIN: Prepping completed successfully.
% VSO_PREP_MAIN: Wrote prepped data to - /home/user/eit.fits

; This command will reconstruct a RHESSI clean 6-12keV image on the PrepServer
vso_prep, instrument='rhessi', im_time_interval=['25-Jan-2007 06:53:44', '25-Jan-2007 06:57:40'], image_alg='clean', im_energy_binning=[6,12]
% VSO_PREP_MAIN: Prepping completed successfully.
% VSO_PREP_MAIN: Wrote prepped data to - /home/user/prepped_temp.fits

Step 3: Visualizing

The data are visualized with PLOTMAN, which handles different data types such as light curves, images, spectra, and spectrograms. It provides basic display operations such as zooming, image overlays, solar rotation, etc. To display data with PLOTMAN from within SHOW_SYNOP, simply select a data file (Figure 2, blue box) and click on "Display". If the data file has already been pre-processed, it will be displayed immediately in a PLOTMAN window, otherwise it is first sent to the PrepServer for processing. If multiple file have been selected, all are displayed in the same PLOTMAN window.

PLOTMAN can also be called from the command line:

; This command pre-processes an XRT image, writes it into an object and displays the XRT image in PLOTMAN (see Figure 3)
vso_prep, 'http://kurasuta.cfa.harvard.edu/DATA/2007/01/25/H0600/XRT20070125_065656.3.fits', oprep=xrt
xrt->plotman, /color

; The first command searches the VSO for the TRACE data file (contains multiple images) closest to January 25 06:30. The second command detects the instrument type TRACE and opens a window for image selection (see Figure 4), then pre-processes the specified images. The last command displays all pre-processed images in PLOTMAN (see Figure 5)
urls = vso_files('25-Jan-2007 06:30', instr='trace')
vso_prep, urls, oprep=trace
trace->plotman, /color

; This last example first reconstructs a RHESSI clean 6-12keV image on the PrepServer; then searches the VSO for the EIT data file closest to January 25 2007 6:55 and sends the URL to the PrepServer for pre-processing; and finally opens both images in the same PLOTMAN window (parameter plotman=p) where the RHESSI contours are manually overlaid with the EIT image to create the image in Figure 6.
vso_prep, instrument='rhessi', im_time_interval=['25-Jan-2007 06:53:44', '25-Jan-2007 06:57:40'], image_alg='clean', im_energy_binning=[6,12], oprep=rhessi
urls = vso_files('25-Jan-2007 06:55', inst='eit')
vso_prep, urls, oprep=eit
rhessi->plotman, plotman=p
eit->plotman, plotman=p


Documentation

• SolarSoftware
• Virtual Solar Observatory
• PrepServer
• PLOTMAN
• VSO_FILES
• VSO_PREP
• SOCK_COPY
• Complete Code Example

Contacts

• SHOW_SYNOP: Dominic Zarro (dominic dot zarro at nasa dot gov)
• PrepServer: Laszlo I. Etesi (laszlo dot etesi at nasa dot gov)
• PLOTMAN: Kim Tolbert (kim dot tolbert at nasa dto gov)