linarphy/naroo_reader
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Update intensity values

Browse source
This commit is contained in:
linarphy 2023-07-04 11:55:05 +02:00
parent 187e72848d
commit 891883e5a8
No known key found for this signature in database
GPG key ID: 3D4AAAC3AD16E79C
1 changed files with 4 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
spectrum.py
View file

@ -561,7 +561,6 @@ bias = {
} }
mean_bias = np.mean( [ bias[ 'top' ] , bias[ 'down' ] ] , axis = 0 ) mean_bias = np.mean( [ bias[ 'top' ] , bias[ 'down' ] ] , axis = 0 )
mean_bias = 0
if verbose: if verbose:
print( 'bias substraction finished' ) print( 'bias substraction finished' )
@ -569,7 +568,7 @@ if verbose:
mean_data = np.mean( data[ mean_data = np.mean( data[
spectrum[ 'y' ][ 'min' ] : spectrum[ 'y' ][ 'max' ], spectrum[ 'y' ][ 'min' ] : spectrum[ 'y' ][ 'max' ],
spectrum[ 'x' ][ 'min' ] : spectrum[ 'x' ][ 'max' ] spectrum[ 'x' ][ 'min' ] : spectrum[ 'x' ][ 'max' ]
] , axis = 0 ) ] , axis = 0 ) - mean_bias
if intensity != None: if intensity != None:
if verbose: if verbose:
@ -689,6 +688,7 @@ if intensity != None:
true_intensity_value = ( intensity_stairs.shape[0] - np.polyval( index_polyval , intensity_value ) ) * step true_intensity_value = ( intensity_stairs.shape[0] - np.polyval( index_polyval , intensity_value ) ) * step
final_intensity[index] = np.exp( true_intensity_value ) final_intensity[index] = np.exp( true_intensity_value )
if final_intensity[ index ] > np.exp( 2.2 ):
if verbose: if verbose:
print( 'intensity calibration finished' ) print( 'intensity calibration finished' )
@ -701,6 +701,8 @@ if output == None:
else: else:
if verbose: if verbose:
print( 'storing result in ' + output ) print( 'storing result in ' + output )
plt.plot( final_intensity[ 1 : - 1 ] )
plt.show()
main_hdu = fits.PrimaryHDU( final_intensity[1:-1] ) # remove -1 main_hdu = fits.PrimaryHDU( final_intensity[1:-1] ) # remove -1
main_hdu.header[ 'CRVAL1' ] = wavelengths[0] main_hdu.header[ 'CRVAL1' ] = wavelengths[0]
main_hdu.header[ 'CDELT1' ] = wavelengths[1] - wavelengths[0] main_hdu.header[ 'CDELT1' ] = wavelengths[1] - wavelengths[0]