Created
August 29, 2012 11:35
-
-
Save dpshelio/3511226 to your computer and use it in GitHub Desktop.
Access to HFC using IDL
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
;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