This page (revision-27) was last changed on 03-Aug-2023 16:53 by Peter Young

This page was created on 02-Apr-2013 16:51 by Peter Young

Only authorized users are allowed to rename pages.

Only authorized users are allowed to delete pages.

Page revision history

Version Date Modified Size Author Changes ... Change note
27 03-Aug-2023 16:53 5 KB Peter Young to previous
26 18-Apr-2021 18:50 5 KB RichardAdacy to previous | to last CPA
25 21-Nov-2016 13:13 5 KB RichardAdacy to previous | to last CPA
24 07-Oct-2016 17:32 5 KB Peter Young to previous | to last Added section about downloading the cutout data
23 23-Mar-2016 18:53 5 KB Peter Young to previous | to last
22 23-Mar-2016 16:16 5 KB Peter Young to previous | to last
21 23-Mar-2016 16:15 5 KB Peter Young to previous | to last

Page References

Incoming links Outgoing links

Version management

Difference between version and

At line 1 changed one line
!!!Making use of AIA data for EIS analysis
[{ALLOW edit EISMainUsers}]
[{ALLOW view Anonymous}]
!!!Procedure for obtaining SDO cutouts from the JSOC service
At line 3 changed one line
(This is work in progress...)
SDO cutouts (both AIA and HMI) can be obtained very quickly through the JSOC. In this article I illustrate how cutouts corresponding to a specific EIS raster can be obtained. Note that the method is generally useful and does not apply specifically to EIS: all you need is a field-of-view size and a pointing location.
At line 5 changed one line
!!Obtaining AIA cutout images
I consider an EIS active region observation that began at 02:30 UT on 15-Apr-2011. The raster took 62 mins to complete. We would like AIA 193 cutouts at 1 min cadence for the period 02:30 to 03:32, and HMI magnetograms at 45s cadence for the same period.
At line 7 changed one line
The request for AIA cutout images can be made directly from IDL as seen from the example below:
Note that Lockheed also provide a cutout service, but this is much slower than the JSOC and does not return HMI data.
At line 9 removed one line
IDL> ssw_cutout_service,'8-feb-2011 18:38','8-feb-2011 19:50',cgiquery,queryout,ref_helio='S30E36',ref_time='8-feb-2011 18:38',fovx=100,fovy=100,wave='193',email=email_address,max_frames=73
At line 11 changed one line
This requests 1 min cadence 193 A images for the time period 18:38-19:50 on 8-Feb-2011. The FOV is 100" x 100" and the center of the FOV is given by ref_helio in longitude-latitude coordinates. The reference time at which the pointing coordinates are specified is given by ref_time. It's important to give an e-mail address so that you receive notification when the cutouts are ready.
!!Get EIS pointing and FOV size
At line 13 changed one line
If your coordinates are in heliocentric units, then you can convert them to longitude-latitude by
The wiki page [Obtaining pointing information from EIS data|EISPointing] explains how to obtain the center and FOV of an EIS raster, so you can follow the instructions there. You should get values of (xcen,ycen)=(199,305), and (fovx,fovy)=(120,160).
At line 15 changed one line
IDL> lonlat=xy2lonlat((x,y),ref_time)
!!Request AIA 193 cutouts
At line 18 added one line
Go to the [JSOC exportdata webpage|http://jsoc.stanford.edu/ajax/exportdata.html]. In the 'RecordSet' box, type:
At line 20 added 109 lines
{{{
aia.lev1_euv_12s[2011.04.15_02:30/62m@1m][193]
}}}
Click on the 'Recount' button and make sure a non-zero number of Records are found.
For 'Method', choose 'url-tar' from the drop down list. You will see a new option called 'Processing' that appears. Tick the box.
Within the processing options, select 'im-patch - Extract sub-frame'. A new box appears with various input widgets, some of which will be automatically filled in.
Make sure 'Tracking' is ticked, and also tick the 'Register' option.
For the other fields, type:
LocUnits: 'arcsec from center'\\
T_Ref: '2011-04-15T03:01:02Z' (choose time from middle of raster)\\
X: 199\\
Y: 305\\
BoxUnits: pixels\\
Width: 300 (some contingency has been factored in here)\\
Height: 400 (some contingency has been factored in here)
For 'Protocol' leave it as 'FITS'. Type in your name and e-mail, and then click on 'Submit Export Request'. ''If you're doing this for the first time, then you will have to register.'' After a couple of minutes wait, click on 'Submit status request' and you should see a link to a tar file containing your data (it may take longer for large data-sets).
Note that cutouts will be rice-compressed images (there used to be an option to download uncompressed images, but this has been removed). See below for how to read compressed images quickly.
Users should be aware that the AIA de-spiking routine (which is automatically applied to all files) can remove real data, which particularly affects small-scale, dynamic phenomena. The spikes can be recovered, however, so please check the page [re-spiking AIA images|AIARespike].
!!Request HMI cutouts
Repeat the above procedure, except that in the RecordSet box you should do:
{{{
hmi.M_45s[2011.04.15_02:30/62m@45s]
}}}
For the field of view, you need to account for the difference in plate scale between AIA and HMI, so multiply by 1.2:
Width: 360\\
Height: 480
!!Downloading cutout images from IDL
After you have submitted your cutout request to the JSOC you will receive a request_id of the form 'JSOC_20161007_486'. You can check the status of your request from IDL by doing:
{{{
IDL> sdo_jsoc_check_status, request_id
}}}
which will check to see if the tar file is available for downloading. If it is, then you should see a message like:
{{{
The data are available at the webpage:
http://jsoc.stanford.edu/SUM9/D862749672/S00000
}}}
You can download and unpack the tar file to your current directory by doing:
{{{
IDL> sdo_jsoc_check_status, request_id, /download
}}}
Note that the tar file is automatically removed after unpacking.
!!Reading cutout images
A cutout image can be read with:
{{{
IDL> read_sdo, file, index, data, /use_shared_lib
}}}
The /use_shared_lib keyword is highly recommended as it greatly increases the speed at which rice-compressed images are read. It may not work on all computers though (Windows?).
!!Further information
Shea Hess Webber and Peter Young have written a [document|https://pyoung.org/gsfc/JSOC%20Export%20Process.pdf] giving more information about using the JSOC service.
!!Comments
''Why use pixels instead of arcsec for the FOV size?''\\
I found erratic results when using arcsec - the HMI FOV would end up with a different size to AIA.
''Why use the register option?''\\
This performs sub-pixel interpolations so that when you make a movie from the cutouts you do not see the occasional pixel jumps that otherwise would be seen. This isn't recommended if you plan to re-spike your cutouts.
''What if I want multiple AIA filters?''\\
Instead of giving {{{[193]}}} in the RecordSet box, give e.g. {{{[171,304,335]}}}.
''How do I get AIA UV images?''\\
For these you need to give 'aia.lev1_uv_24s' in the RecordSet box.
''How do I get HMI 45s intensitygrams, or HMI 720s magnetograms?''\\
In the RecordSet box, use 'hmi.Ic_45s' and 'hmi.M_720s'.
''I want to respike my cutouts. Can I get spike cutouts?''\\
Please check the page [re-spiking AIA images|AIARespike].
''What about HARPs? How do I extract information for these?''\\
A HARP is a HMI Active Region Patch, i.e., a cutout of HMI data that has been processed to yield, e.g., the vector magnetic field. See the separate page on [Extracting HARP data from the JSOC|JsocHarp].
''Can I get header information from JSOC without downloading data?''\\
Putting the following URL into your browser will give a text page with the specified keywords (T_OBS, WAVELNTH, EXPTIME and NSPIKES) given as columns.
{{{
http://jsoc.stanford.edu/cgi-bin/ajax/show_info?ds=aia.lev1_euv_12s[2011.04.15_02:30/62m@1m][193]&key=T_OBS,WAVELNTH,EXPTIME,NSPIKES
}}}