Skip to content

Instantly share code, notes, and snippets.

@dpshelio
Created August 29, 2012 11:35
Show Gist options
  • Save dpshelio/3511226 to your computer and use it in GitHub Desktop.
Save dpshelio/3511226 to your computer and use it in GitHub Desktop.
Access to HFC using IDL
;set time range:
timerange = ['2001-04-15','2001-04-16']
cmd_url = hio_form_hfc_query(timerange=timerange,listname='VIEW_AR_HQI')
smart_ars = decode_votable(ssw_hio_query(cmd_url,/verb),/quiet)
; we should have the url on the database, but it seems it's not there (I'm taking note for that!!)
; but we can work around this way:
mdimags=mdi_time2file(timerange[0],timerange[1],/stanford)
lab = where(smart_ars[0].filename eq ssw_strsplit(mdimags,'/',/tail),num)
sock_copy,mdimags[lab],local_file=filename
fits2map,filename,mdimap
; Now we plot the mdi image from where these detections were made as it follows:
; Plot the image
plot_map,mdimap,/limb,grid=20
; plot the first detection
hfr_oplot_feature,smart_ars[0],ref_map=mdimap
; if we want to plot all the detections on that image, one way to do it is:
; 1) find the detections with the same filename:
ars_image = where(smart_ars.filename eq smart_ars[0].filename,num_ars_image)
; 2) plot them all!
hfr_oplot_feature,smart_ars[ars_image],ref_map=mdimap
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment