search

foxsi / foxsi-smex

Tools for scientists to investigate the capabilities of FOXSI
MIT License
2 stars 8 forks source link

Point Spread Function #20

Closed ehsteve closed 8 years ago

ehsteve commented 8 years ago

First draft of psf. Not yet completed.

ehsteve commented 8 years ago

@aringlis Here is my first draft of the PSF work in python. I think you can use this to get started porting this over to IDL. Thanks.

aringlis commented 8 years ago

Will do.

LinErinG commented 8 years ago

@aringlis @ehsteve Note that Sam left a convenient slot for the PSF in his framework, a function entitled "foxsi_get_psf_array." If the IDL version of the PSF code is slid into that function (which was kept quite general but which has a placeholder PSF there at the moment), it should plug-and-play with the rest of the routines. A clarification is that "away from solar center" means "off-axis angle;" that terminology could be changed.

aringlis commented 8 years ago

@LinErinG - thanks, I will slot the new IDL code in there. I may tweak the output slightly to return an SSW map of the PSF.

aringlis commented 8 years ago

@ehsteve are pitch and yaw astropy quantities, or just variables? If the latter, I'm wondering what the expected units are.

ehsteve commented 8 years ago

Mostly finished and needs review.

ehsteve commented 8 years ago

This should be finished. Here is some code to test

from pyfoxsi.psf import convolve
from sunpy.map import Map
import matplotlib.pyplot as plt

f = '/Users/schriste/Google Drive/Work/FOXSI SMEX/Data/hsi_image_20050513_164526to164626_pixon_3arcsecx64_25to50kev_d1to9.fits'
hsi = Map(f)

foxsi_map = convolve(hsi)

foxsi_map.plot()
plt.show()
ehsteve commented 8 years ago

figure_1

ayshih commented 8 years ago

After @ehsteve merges my PR to his branch, the code will produce this (with oversample_psf set to 3): figure_1