@@ -534,78 +534,78 @@ def trimarray(array):
534
534
self .signalShowMessage .emit (f"Fixing image header for { expression } " , 3000 )
535
535
536
536
# get NAXIS diff
537
- _max = template .fits_header .get ('NAXIS' )
538
- _diff = numpy .ndim (result )
537
+ _naxis = template .fits_header .get ('NAXIS' )
538
+ _ndims = numpy .ndim (result )
539
539
540
540
# sub() WCS to new NAXIS
541
541
_wcs = WCS (template .fits_header )
542
- _new_wcs = _wcs .wcs .sub (_diff )
542
+ _new_wcs = _wcs .wcs .sub (_ndims )
543
543
544
544
# create new header
545
545
_new_header = _new_wcs .to_header ()
546
546
_header = Header .fromstring (_new_header )
547
547
548
548
# create keyword ignore list
549
549
ignore_list = []
550
- for i in range (_diff ):
551
- ignore_list .append (f'NAXIS{ i + 1 + _diff } ' )
552
- ignore_list .append (f'CTYPE{ i + 1 + _diff } ' )
553
- ignore_list .append (f'CRVAL{ i + 1 + _diff } ' )
554
- ignore_list .append (f'CDELT{ i + 1 + _diff } ' )
555
- ignore_list .append (f'CUNIT{ i + 1 + _diff } ' )
556
- ignore_list .append (f'CRPIX{ i + 1 + _diff } ' )
550
+ for i in range (_ndims ):
551
+ ignore_list .append (f'NAXIS{ i + 1 + _ndims } ' )
552
+ ignore_list .append (f'CTYPE{ i + 1 + _ndims } ' )
553
+ ignore_list .append (f'CRVAL{ i + 1 + _ndims } ' )
554
+ ignore_list .append (f'CDELT{ i + 1 + _ndims } ' )
555
+ ignore_list .append (f'CUNIT{ i + 1 + _ndims } ' )
556
+ ignore_list .append (f'CRPIX{ i + 1 + _ndims } ' )
557
557
ignore_list .append ('HISTORY' )
558
558
ignore_list .append ('COMMENT' )
559
559
560
560
# diff the two headers
561
- d_header = pyfits .HeaderDiff (template .fits_header ,
562
- _header ,
563
- ignore_keywords = ignore_list ,
564
- ignore_comments = '*' )
561
+ _diff_header = pyfits .HeaderDiff (template .fits_header ,
562
+ _header ,
563
+ ignore_keywords = ignore_list ,
564
+ ignore_comments = '*' )
565
565
566
566
# update new header with missing cards
567
- for key in d_header .diff_keywords [0 ]:
567
+ for key in _diff_header .diff_keywords [0 ]:
568
568
_header [key ] = template .fits_header [key ]
569
569
570
570
# check for PC and remove
571
571
# a brute force approach for now
572
572
if template .fits_header .get ('PC01_01' ):
573
- for i in range (_max ):
574
- for n in range (_max ):
573
+ for i in range (_naxis ):
574
+ for n in range (_naxis ):
575
575
try :
576
- _header .remove (f'PC0{ i + 1 + _diff } _0{ n + 1 } ' )
576
+ _header .remove (f'PC0{ i + 1 + _ndims } _0{ n + 1 } ' )
577
577
except :
578
578
continue
579
- for i in range (_max ):
580
- for n in range (_max ):
579
+ for i in range (_naxis ):
580
+ for n in range (_naxis ):
581
581
try :
582
- _header .remove (f'PC0{ i + 1 } _0{ n + 1 + _diff } ' )
582
+ _header .remove (f'PC0{ i + 1 } _0{ n + 1 + _ndims } ' )
583
583
except :
584
584
continue
585
585
elif template .fits_header .get ('PC1_1' ):
586
- for i in range (_max ):
587
- for n in range (_max ):
586
+ for i in range (_naxis ):
587
+ for n in range (_naxis ):
588
588
try :
589
- _header .remove (f'PC{ i + 1 + _diff } _{ n + 1 } ' )
589
+ _header .remove (f'PC{ i + 1 + _ndims } _{ n + 1 } ' )
590
590
except :
591
591
continue
592
- for i in range (_max ):
593
- for n in range (_max ):
592
+ for i in range (_naxis ):
593
+ for n in range (_naxis ):
594
594
try :
595
- _header .remove (f'PC{ i + 1 } _{ n + 1 + _diff } ' )
595
+ _header .remove (f'PC{ i + 1 } _{ n + 1 + _ndims } ' )
596
596
except :
597
597
continue
598
598
elif template .fits_header .get ('PC001001' ):
599
- for i in range (_max ):
600
- for n in range (_max ):
599
+ for i in range (_naxis ):
600
+ for n in range (_naxis ):
601
601
try :
602
- _header .remove (f'PC00{ i + 1 + _diff } 00{ n + 1 } ' )
602
+ _header .remove (f'PC00{ i + 1 + _ndims } 00{ n + 1 } ' )
603
603
except :
604
604
continue
605
- for i in range (_max ):
606
- for n in range (_max ):
605
+ for i in range (_naxis ):
606
+ for n in range (_naxis ):
607
607
try :
608
- _header .remove (f'PC00{ i + 1 } 00{ n + 1 + _diff } ' )
608
+ _header .remove (f'PC00{ i + 1 } 00{ n + 1 + _ndims } ' )
609
609
except :
610
610
continue
611
611
@@ -624,9 +624,9 @@ def trimarray(array):
624
624
# restore card comments
625
625
for card in template .fits_header .keys ():
626
626
if card in _header :
627
- com = template .fits_header .comments [card ]
628
- if com :
629
- _header .set (keyword = card , comment = com )
627
+ _comm = template .fits_header .comments [card ]
628
+ if _comm :
629
+ _header .set (keyword = card , comment = _comm )
630
630
631
631
# set new header
632
632
template .fits_header = _header
0 commit comments