Preamble (please read)

Contents

  1. Determining which data you want
    1. I need to browse through some data before I can decide what data I'd like to request.
    2. I'd like a more detailed look instead of just a "browse."
    3. I'm looking for a specific feature, event or phenomenon but I haven't determined the date/times of the images I need.
    4. I know the start time and end time of the event or phenomenon I'd like to study, but I don't know which wavelength(s).
  2. Accessing the data
    1. I want images and/or a movie, and I don't need the full science data files.
    2. I need FITS data files, and I think I know the time range and wavelengths of the images I want.
    3. I need FITS data files, I know the time range and wavelengths of the images I want, but the available data volume is overwhelming. Case 1: Using subfields/cutouts
      1. I need FITS data files, I know the time range and wavelengths of the images I want, but the available volume is overwhelming. Case 2: Decreasing the temporal or spatial resolution
  3. Analyzing the data
    1. Well, well, well. So what have we here?
    2. What analysis tools are available?
    3. Section III.c: What are the limitations of the data?
    4. Section III.d: I'm ready to publish. Who do I acknowledge?
  4. Appendix I: IDL and the SolarSoft software tree
    1. Appendix I.a: What is IDL?
    2. Appendix I.b: What is SolarSoft?
    3. Appendix I.c: I have IDL, but I don't have the SolarSoft software tree on my computer.
    4. Appendix I.d: I have the SolarSoft software tree on my computer, but I'm not sure if I have the appropriate software packages installed.
  5. Appendix II: Identifying and retrieving data using the Virtual Solar Observatory's IDL interface
    1. Appendix II.a: Searching the VSO catalog using vso_search.pro, with several examples
    2. Appendix II.b: Retrieving data with vso_get.pro
  6. Appendix III: vso_search.pro documentation

Determining which data you want

I need to browse through some data before I can decide what data I'd like to request.

I'd like a more detailed look instead of just a "browse."

I'm looking for a specific feature, event or phenomenon but I haven't determined the date/times of the images I need.

I know the start time and end time of the event or phenomenon I'd like to study, but I don't know which wavelength(s).

Accessing the data

I want images and/or a movie, and I don't need the full science data files.

I need FITS data files, and I think I know the time range and wavelengths of the images I want.

I need FITS data files, I know the time range and wavelengths of the images I want, but the available data volume is overwhelming. Case 1: Using subfields/cutouts

I need FITS data files, I know the time range and wavelengths of the images I want, but the available volume is overwhelming. Case 2: Decreasing the temporal or spatial resolution

  IDL> searchfile = vso_search('2010-jul-20 14:00','2010-jul-20 14:10', instr='aia', wave='304')
  Records Returned : JSOC : 50/50
Records Returned : JSOC : 7/7

  IDL> searchfile = vso_search('2010-jul-20 14:00','2010-jul-20 14:10', instr='aia', wave='304', pixel=4096)
  Records Returned : JSOC : 50/50
  IDL> searchfile = vso_search('2010-jul-20 14:00','2010-jul-20 14:10', instr='aia', wave='304', pixel=1024)
  Records Returned : JSOC : 7/7

  IDL> searchfile = vso_search('2010-jul-20 14:00','2010-jul-20 14:10', instr='aia', wave='304', pixel=4096, sample=180)
  Records Returned : JSOC : 4/4

  IDL> getfile = vso_get(searchfile[6],filelist=filelist)
  VSO_GET: This will download 1 file(s)
  1 : http://vso.tuc.noao.edu/cgi-bin/drms_test/drms_export.cgi?series=aia_synoptic2;record=304_11763447-11763447
  % RDWRT_BUFF: Please wait. Downloading...
  % File: /cgi-bin/drms_test/drms_export.cgi?series=aia_synoptic2;record=304_11763447-11763447
  % Size: 4210560 bytes
  % From: vso.tuc.noao.edu
  % To: /Users/barbara/data
  % HTTP::COPY: 4210560 bytes of 4210560 total bytes copied in     2.73 seconds
  % HTTP::COPY: Wrote 4210560 bytes to file /Users/barbara/data/aia_test.synoptic2.304A_2010-07-20T14_10_02.12Z.image_lev1p5.fits
  Downloading completed

  IDL> print, filelist[0]
  /Users/barbara/data/aia_test.synoptic2.304A_2010-07-20T14_10_02.12Z.image_lev1p5.fits

Analyzing the data

Well, well, well. So what have we here?

What analysis tools are available?

Section III.c: What are the limitations of the data?

Section III.d: I'm ready to publish. Who do I acknowledge?

Appendix I: IDL and the SolarSoft software tree

Appendix I.a: What is IDL?

Appendix I.b: What is SolarSoft?

Appendix I.c: I have IDL, but I don't have the SolarSoft software tree on my computer.

Appendix I.d: I have the SolarSoft software tree on my computer, but I'm not sure if I have the appropriate software packages installed.

SolarSoft library and add new packages. This guide assumes you have the "AIA", "VSO", and "ONTOLOGY" packages installed.

Important Note: It is very important to do a regular update of the Solarsoft tree. Modifications to the tree's software are happening continually, and unless you do an update you're not going to have the most recent versions. If possible, set your computer to automatically do updates on a regular (e.g. weekly) basis.

Appendix II: Identifying and retrieving data using the Virtual Solar Observatory's IDL interface

Appendix II.a: Searching the VSO catalog using vso_search.pro, with several examples

I DL input code is shown next to the IDL> prompt, and the text response to the input (if any) is shown below the command line. We'll start by defining a time range to work with:

  IDL> time1 = '2010-jul-17 00:00'
  IDL> time2 = '2010-jul-17 01:00'

Example 1. Find all of the AIA images available for a given date range:

  IDL> searchfile = vso_search(time1, time2, instr='aia')
  Records Returned : JSOC : 2306/2306
  Records Returned : JSOC : 300/300
  IDL> print_struct, searchfile[0]
     DETECTOR  INSTRUMENT SOURCE  PROVIDER   INFO    PHYSOBS     FILEID        SIZE  URL  GETINFO
     AIA  SDO  JSOC  AIA level 14096x4096  intensity  aia_lev1:304:1058400038    6.620E+04
  IDL> print_struct, searchfile[0].time
         START                 _END
         2010-07-17T00:00:02  2010-07-17T00:00:03

  IDL> searchfile = vso_search(instr='aia', wave='193', pixels = 4096,/latest)
  Records Returned : JSOC : 1/1

  resolution=0.25
  pixels=1024
  pscale=2.4

  IDL>  searchfile = vso_search(near=time1, instr='aia', pixel=1024, wave='171')
  Records Returned : JSOC : 1/1
  IDL>  print_struct, searchfile.time
          START                 _END
      2010-07-16T23:59:59  2010-07-16T24:00:00

  IDL>  searchfile = vso_search('2010-jul-16 03:00','2010-jul-16 07:00', instr='aia', wave='193', sample=600)
  Records Returned : JSOC : 24/24
  Records Returned : JSOC : 24/24

  IDL>  searchfile = vso_search(near=time1, instr='aia', wave='171 - 193', sample=600)
  Records Returned : JSOC : 24/24
  Records Returned : JSOC : 24/24

Appendix II.b: Retrieving data with vso_get.pro

directory)

Appendix III: vso_search.pro documentation

Email: <you AT SPAMFREE example DOT com>

...


JsocWiki: AccessingAIAData (last edited 2013-05-01 04:35:23 by localhost)