EOVSA Data Products - Revision history

Bchen: /* SSWIDL Users */

2025-06-16T10:39:47Z

SSWIDL Users

← Older revision		Revision as of 10:39, 16 June 2025
Line 188:		Line 188:

	====SSWIDL Users====		====SSWIDL Users====
	For event-based data products, the only difference is in how the time array is read and converted. In these FITS files, the time array is provided in Julian Date (JD), not in SSWIDL anytim format. Use the following function to convert JD to SSWIDL “anytim” seconds, then proceed with the same spectrogram analysis and plotting steps as in the synoptic example above.		For event-based data products, the only difference is in how the time array is read and converted. In these FITS files, the time array is provided in Julian Date (JD), not in SSWIDL anytim format. Use the following function to convert JD to SSWIDL “anytim” seconds, then proceed with the same spectrogram analysis and plotting steps as in the synoptic example above. Note that for most event-based spectrograms, background-subtraction has already been performed.

	<syntaxhighlight lang="idl">		<syntaxhighlight lang="idl">
Line 207:		Line 207:
	timeut = jd2anytim(time.time) ; Convert JD to SSWIDL "anytim" seconds		timeut = jd2anytim(time.time) ; Convert JD to SSWIDL "anytim" seconds

	; The rest of your workflow (creating the spectrogram, plotting~~, background subtraction, plotman, etc.~~)		; The rest of your workflow (creating the spectrogram and plotting)
	; is identical to the synoptic example above.		; is identical to the synoptic example above.
	</syntaxhighlight>		</syntaxhighlight>

Sjyu: /* SSWIDL Users */

2025-06-16T02:30:14Z

SSWIDL Users

← Older revision		Revision as of 02:30, 16 June 2025
Line 122:		Line 122:
	====SSWIDL Users====		====SSWIDL Users====
	For IDL users, here is an example of using SSWIDL's "[https://hesperia.gsfc.nasa.gov/ssw/gen/idl/fits/mrdfits.pro mrdfits]" to read the total-power FITS file and the "[https://hesperia.gsfc.nasa.gov/ssw/gen/idl/objects/spectrogram.pro spectrogram]" function to visualize it. We use the 2021 Oct 28 X1 flare as example, downloaded from [https://ovsa.njit.edu/browser/?suntoday_date=2021-10-28 the browser].		For IDL users, here is an example of using SSWIDL's "[https://hesperia.gsfc.nasa.gov/ssw/gen/idl/fits/mrdfits.pro mrdfits]" to read the total-power FITS file and the "[https://hesperia.gsfc.nasa.gov/ssw/gen/idl/objects/spectrogram.pro spectrogram]" function to visualize it. We use the 2021 Oct 28 X1 flare as example, downloaded from [https://ovsa.njit.edu/browser/?suntoday_date=2021-10-28 the browser].
	<~~pre~~>		<syntaxhighlight lang="idl">
	; Read the FITS file using mrdfits		; Read the FITS file using mrdfits
	filename = 'EOVSA_TPall_20211028.fts'		filename = 'EOVSA_TPall_20211028.fts'
Line 146:		Line 146:
	loadct,5		loadct,5
	s.plot,/log,timerange=['2021-10-28T15:15:00','2021-10-28T16:00:00'],/xsty,/ysty,ytitle='Frequency [GHz]',charsize=1.5		s.plot,/log,timerange=['2021-10-28T15:15:00','2021-10-28T16:00:00'],/xsty,/ysty,ytitle='Frequency [GHz]',charsize=1.5
	</~~pre~~>		</syntaxhighlight>
	[[File:IDL_synop_tp_20211028.jpg\|center\|500px]]		[[File:IDL_synop_tp_20211028.jpg\|center\|500px]]

	<~~pre~~>		<syntaxhighlight lang="idl">
	; If you wish to do a pre-flare background subtraction and bring out more details of the flare		; If you wish to do a pre-flare background subtraction and bring out more details of the flare
	timebkg = ['2021-10-28T15:15', '2021-10-28T15:20'] ; selecting a background window		timebkg = ['2021-10-28T15:15', '2021-10-28T15:20'] ; selecting a background window
Line 162:		Line 162:
	loadct,5		loadct,5
	s.plot,/linear,timerange=['2021-10-28T15:15:00','2021-10-28T16:00:00'],/xsty,/ysty,ytitle='Frequency [GHz]',charsize=1.5		s.plot,/linear,timerange=['2021-10-28T15:15:00','2021-10-28T16:00:00'],/xsty,/ysty,ytitle='Frequency [GHz]',charsize=1.5
	</~~pre~~>		</syntaxhighlight>
	[[File:IDL_synop_tp_20211028_bkgsub.jpg\|center\|500px]]		[[File:IDL_synop_tp_20211028_bkgsub.jpg\|center\|500px]]

	<~~pre~~>		<syntaxhighlight lang="idl">
	; One can also load the object into plotman and visualize it interactively. See the snapshot below for an example.		; One can also load the object into plotman and visualize it interactively. See the snapshot below for an example.
	; The widget allows a convenient way to explore the light curves by selecting "plot_Control" -> "Image or Spectrogram Profiles" -> "Rows or Columns")		; The widget allows a convenient way to explore the light curves by selecting "plot_Control" -> "Image or Spectrogram Profiles" -> "Rows or Columns")
	s->set,timerange = ['2021-10-28T15:15:00','2021-10-28T16:00:00']		s->set,timerange = ['2021-10-28T15:15:00','2021-10-28T16:00:00']
	s->plotman		s->plotman
	</~~pre~~>		</syntaxhighlight>
	[[File:IDL_synop_tp_20211028_bkgsub_plotman.jpg\|center\|500px]]		[[File:IDL_synop_tp_20211028_bkgsub_plotman.jpg\|center\|500px]]


	For synoptic images:		For synoptic images:
	<~~pre~~>		<syntaxhighlight lang="idl">
	read_sdo,'eovsa_20191225.spw11-20.tb.disk.fits',header,data,/UNCOMP_DELETE		read_sdo,'eovsa_20191225.spw11-20.tb.disk.fits',header,data,/UNCOMP_DELETE
	index2map,header,data,eomap		index2map,header,data,eomap
	plot_map,eomap		plot_map,eomap
	</~~pre~~>		</syntaxhighlight>

	[[File:eovsa_20191225_image_sswidl.jpg\| center \|250px]]		[[File:eovsa_20191225_image_sswidl.jpg\| center \|250px]]


	===Event-Based Data Products===		===Event-Based Data Products===

Sjyu: /* SSWIDL Users */

2025-06-16T02:29:23Z

SSWIDL Users

← Older revision		Revision as of 02:29, 16 June 2025
Line 191:		Line 191:
	For event-based data products, the only difference is in how the time array is read and converted. In these FITS files, the time array is provided in Julian Date (JD), not in SSWIDL anytim format. Use the following function to convert JD to SSWIDL “anytim” seconds, then proceed with the same spectrogram analysis and plotting steps as in the synoptic example above.		For event-based data products, the only difference is in how the time array is read and converted. In these FITS files, the time array is provided in Julian Date (JD), not in SSWIDL anytim format. Use the following function to convert JD to SSWIDL “anytim” seconds, then proceed with the same spectrogram analysis and plotting steps as in the synoptic example above.

	<~~pre~~>		<syntaxhighlight lang="idl">
	; Function to convert Julian Date (JD) to SSWIDL "anytim" seconds (since 1979-01-01 00:00:00 UTC)		; Function to convert Julian Date (JD) to SSWIDL "anytim" seconds (since 1979-01-01 00:00:00 UTC)
	function jd2anytim, jd		function jd2anytim, jd
Line 210:		Line 210:
	; The rest of your workflow (creating the spectrogram, plotting, background subtraction, plotman, etc.)		; The rest of your workflow (creating the spectrogram, plotting, background subtraction, plotman, etc.)
	; is identical to the synoptic example above.		; is identical to the synoptic example above.
	</~~pre~~>		</syntaxhighlight>

	'''Tip:''' After <code>timeut = jd2anytim(time.time)</code>, you can continue with the code from the synoptic data section—creating the spectrogram, plotting, background subtraction, and interactive exploration with <code>plotman</code>—using the same commands.		'''Tip:''' After <code>timeut = jd2anytim(time.time)</code>, you can continue with the code from the synoptic data section—creating the spectrogram, plotting, background subtraction, and interactive exploration with <code>plotman</code>—using the same commands.

Sjyu: /* SSWIDL Users */

2025-06-16T02:19:03Z

SSWIDL Users

@@ Line 189: / Line 189: @@
 ====SSWIDL Users====
-We are working on an example with SSWIDL and will post it soon.
+For event-based data products, the only difference is in how the time array is read and converted. In these FITS files, the time array is provided in Julian Date (JD), not in SSWIDL anytim format. Use the following function to convert JD to SSWIDL “anytim” seconds, then proceed with the same spectrogram analysis and plotting steps as in the synoptic example above.
 =Requesting EOVSA Data or Analysis Assistance=

Sjyu: /* Reading and Using level 1 data */

2025-06-16T02:11:37Z

Reading and Using level 1 data

@@ Line 115: / Line 115: @@
 ===Introduction===
 All our level 1 data products are in FITS format. All the images have standard, WCS-compatible coordinates. Users can use their favorite method to read these files. In the following, we provide minimal examples to read them with Astropy and Sunpy.
 ===Synoptic Data Products===
@@ Line 189: / Line 182: @@
 [[File:eovsa_20191225_image_sswidl.jpg| center |250px]]
 =Requesting EOVSA Data or Analysis Assistance=

Bchen: /* SSWIDL users */

2025-06-15T13:27:11Z

SSWIDL users

← Older revision		Revision as of 13:27, 15 June 2025
Line 127:		Line 127:
	An example of how to read and plot the flare spectrograms and images in Python (with [https://www.astropy.org/ Astropy], [https://sunpy.org/ SunPy], and optionally, [https://docs.sunpy.org/projects/radiospectra/en/latest/ radiospectra]) can be accessed at [https://colab.research.google.com/drive/1bF_WjKRk51Hb3h10EzcGaBCw9QFjWwqW?usp=sharing this Google Colab Jupyter notebook].		An example of how to read and plot the flare spectrograms and images in Python (with [https://www.astropy.org/ Astropy], [https://sunpy.org/ SunPy], and optionally, [https://docs.sunpy.org/projects/radiospectra/en/latest/ radiospectra]) can be accessed at [https://colab.research.google.com/drive/1bF_WjKRk51Hb3h10EzcGaBCw9QFjWwqW?usp=sharing this Google Colab Jupyter notebook].

	====SSWIDL ~~users~~====		====SSWIDL Users====
	For IDL users, here is an example of using SSWIDL's "[https://hesperia.gsfc.nasa.gov/ssw/gen/idl/fits/mrdfits.pro mrdfits]" to read the total-power FITS file and the "[https://hesperia.gsfc.nasa.gov/ssw/gen/idl/objects/spectrogram.pro spectrogram]" function to visualize it. We use the 2021 Oct 28 X1 flare as example, downloaded from [https://ovsa.njit.edu/browser/?suntoday_date=2021-10-28 the browser].		For IDL users, here is an example of using SSWIDL's "[https://hesperia.gsfc.nasa.gov/ssw/gen/idl/fits/mrdfits.pro mrdfits]" to read the total-power FITS file and the "[https://hesperia.gsfc.nasa.gov/ssw/gen/idl/objects/spectrogram.pro spectrogram]" function to visualize it. We use the 2021 Oct 28 X1 flare as example, downloaded from [https://ovsa.njit.edu/browser/?suntoday_date=2021-10-28 the browser].
	<pre>		<pre>

Bchen: /* Event-Based Data Products */

2025-06-15T13:26:59Z

Event-Based Data Products

← Older revision		Revision as of 13:26, 15 June 2025
Line 117:		Line 117:

	===Event-Based Data Products===		===Event-Based Data Products===
	* An example of how to read and plot the flare spectrograms and images in Python (with [https://www.astropy.org/ Astropy] and [https://sunpy.org/ SunPy]) can be accessed at [https://colab.research.google.com/drive/1wMvuxuNip5cJJHMoTOLT6u_OdWobpBY0?usp=sharing this Google Colab Jupyter notebook].		====Python Users====
	* We are working on an example with SSWIDL and will post it soon.		An example of how to read and plot the flare spectrograms and images in Python (with [https://www.astropy.org/ Astropy] and [https://sunpy.org/ SunPy]) can be accessed at [https://colab.research.google.com/drive/1wMvuxuNip5cJJHMoTOLT6u_OdWobpBY0?usp=sharing this Google Colab Jupyter notebook].

			====SSWIDL Users====
			We are working on an example with SSWIDL and will post it soon.

	===Synoptic Data Products===		===Synoptic Data Products===

Bchen: /* SSWIDL users */

2025-06-15T13:25:20Z

SSWIDL users

← Older revision		Revision as of 13:25, 15 June 2025
Line 154:		Line 154:

	<pre>		<pre>
	; If you wish to do a pre-flare background subtraction		; If you wish to do a pre-flare background subtraction and bring out more details of the flare
	timebkg = ['2021-10-28T15:15', '2021-10-28T15:20']		timebkg = ['2021-10-28T15:15', '2021-10-28T15:20'] ; selecting a background window
	dt0 = min(abs(timeut - anytim(timebkg[0])), bkg_i0)		dt0 = min(abs(timeut - anytim(timebkg[0])), bkg_i0)
	dt1 = min(abs(timeut - anytim(timebkg[1])), bkg_i1)		dt1 = min(abs(timeut - anytim(timebkg[1])), bkg_i1)
	spec_bkg = mean(spec[bkg_i0:bkg_i1,*], dimension=1)		spec_bkg = mean(spec[bkg_i0:bkg_i1,*], dimension=1) ; this is the average background
	spec_sub = spec - rebin(reform(spec_bkg, 1, nfreq), ntime, nfreq)		spec_sub = spec - rebin(reform(spec_bkg, 1, nfreq), ntime, nfreq) ;do the subtraction
	; Recreate the spectrogram and do the plotting		; Recreate the spectrogram and do the plotting
	s = spectrogram(spec_sub, timeut, fghz)		s = spectrogram(spec_sub, timeut, fghz)

Bchen: /* SSWIDL users */

2025-06-15T13:23:57Z

SSWIDL users

← Older revision		Revision as of 13:23, 15 June 2025
Line 170:		Line 170:

	<pre>		<pre>
	; One can also load the object into plotman and visualize it interactively		; One can also load the object into plotman and visualize it interactively. See the snapshot below for an example.
			; The widget allows a convenient way to explore the light curves by selecting "plot_Control" -> "Image or Spectrogram Profiles" -> "Rows or Columns")
	s->set,timerange = ['2021-10-28T15:15:00','2021-10-28T16:00:00']		s->set,timerange = ['2021-10-28T15:15:00','2021-10-28T16:00:00']
	s->plotman		s->plotman

Bchen: /* SSWIDL users */

2025-06-15T13:14:28Z

SSWIDL users

← Older revision		Revision as of 13:14, 15 June 2025
Line 178:		Line 178:

	For synoptic images:		For synoptic images:
	<pre ~~style="background-color: #FCEBD9;"~~>		<pre>
	read_sdo,'eovsa_20191225.spw11-20.tb.disk.fits',header,data,/UNCOMP_DELETE		read_sdo,'eovsa_20191225.spw11-20.tb.disk.fits',header,data,/UNCOMP_DELETE
	index2map,header,data,eomap		index2map,header,data,eomap