In [1]:
import glob
import matplotlib.pyplot as plt
import matplotlib.colors as colors
%matplotlib inline
mapfiles = glob.glob('/Users/schriste/Desktop/sunpy_test_img/*.*')
cubefiles = glob.glob('/Users/schriste/Desktop/sunpy_test_img/cubes/*')
from sunpy.map import Map
import numpy as np
import sunpy
In [2]:
nplots = 0
for ind, f in enumerate(mapfiles):
print(ind, f)
try:
m = Map(f)
print(m.__repr__())
nplots += 1
except:
print('fail for f')
pass
(0, '/Users/schriste/Desktop/sunpy_test_img/20110919_000500_s4c1A.fts')
SunPy CORMap
---------
Observatory: STEREO A
Instrument: SECCHI
Detector: COR1
Measurement: white-light
Wavelength: nan nm
Obs Date: 2011-09-19 00:05:00
dt: 1.500520 s
Dimension: [ 512. 512.] pix
scale: [ 15.00860023 15.00860023] arcsec / pix
array([[ 676., 676., 674., ..., 700., 700., 700.],
[ 676., 676., 674., ..., 700., 700., 700.],
[ 676., 676., 676., ..., 700., 700., 700.],
...,
[ 678., 678., 676., ..., 708., 708., 710.],
[ 678., 678., 678., ..., 708., 708., 710.],
[ 676., 678., 678., ..., 706., 706., 708.]], dtype=float32)
(1, '/Users/schriste/Desktop/sunpy_test_img/20110919_000500_s4c1B.fts')
SunPy CORMap
---------
Observatory: STEREO B
Instrument: SECCHI
Detector: COR1
Measurement: white-light
Wavelength: nan nm
Obs Date: 2011-09-19 00:05:53
dt: 1.699660 s
Dimension: [ 512. 512.] pix
scale: [ 15.00352001 15.00352001] arcsec / pix
array([[ 722., 722., 726., ..., 708., 708., 710.],
[ 726., 726., 726., ..., 710., 708., 710.],
[ 730., 728., 726., ..., 708., 708., 708.],
...,
[ 716., 716., 716., ..., 708., 708., 708.],
[ 716., 716., 716., ..., 706., 708., 708.],
[ 716., 716., 716., ..., 706., 708., 708.]], dtype=float32)
(2, '/Users/schriste/Desktop/sunpy_test_img/20110919_000530_n4euA.fts')
SunPy EUVIMap
---------
Observatory: STEREO A
Instrument: SECCHI
Detector: EUVI
Measurement: 195.0 Angstrom
Wavelength: 195.0 Angstrom
Obs Date: 2011-09-19 00:05:31
dt: 8.003840 s
Dimension: [ 2048. 2048.] pix
scale: [ 1.58777404 1.58777404] arcsec / pix
array([[ 724., 724., 724., ..., 722., 722., 722.],
[ 724., 724., 724., ..., 722., 722., 722.],
[ 724., 724., 724., ..., 722., 722., 722.],
...,
[ 724., 724., 724., ..., 724., 724., 724.],
[ 724., 724., 724., ..., 724., 724., 724.],
[ 724., 724., 724., ..., 722., 722., 722.]], dtype=float32)
(3, '/Users/schriste/Desktop/sunpy_test_img/20110919_000530_n4euB.fts')
SunPy EUVIMap
---------
Observatory: STEREO B
Instrument: SECCHI
Detector: EUVI
Measurement: 195.0 Angstrom
Wavelength: 195.0 Angstrom
Obs Date: 2011-09-19 00:06:24
dt: 8.007300 s
Dimension: [ 2048. 2048.] pix
scale: [ 1.59000003 1.59000003] arcsec / pix
array([[ 686., 688., 688., ..., 694., 696., 698.],
[ 686., 688., 688., ..., 692., 694., 696.],
[ 686., 688., 688., ..., 690., 692., 694.],
...,
[ 686., 686., 686., ..., 686., 686., 686.],
[ 686., 686., 686., ..., 686., 686., 686.],
[ 686., 686., 686., ..., 684., 684., 684.]], dtype=float32)
(4, '/Users/schriste/Desktop/sunpy_test_img/20110919_000800_d4c2A.fts')
SunPy CORMap
---------
Observatory: STEREO A
Instrument: SECCHI
Detector: COR2
Measurement: white-light
Wavelength: nan nm
Obs Date: 2011-09-19 00:08:01
dt: 2.001520 s
Dimension: [ 256. 256.] pix
scale: [ 117.59999847 117.59999847] arcsec / pix
array([[ 684., 684., 684., ..., 684., 684., 684.],
[ 684., 684., 684., ..., 684., 684., 684.],
[ 684., 684., 684., ..., 684., 684., 684.],
...,
[ 684., 684., 684., ..., 684., 684., 684.],
[ 684., 684., 684., ..., 684., 684., 684.],
[ 684., 684., 684., ..., 684., 684., 684.]], dtype=float32)
(5, '/Users/schriste/Desktop/sunpy_test_img/20110919_000800_d4c2B.fts')
SunPy CORMap
---------
Observatory: STEREO B
Instrument: SECCHI
Detector: COR2
Measurement: white-light
Wavelength: nan nm
Obs Date: 2011-09-19 00:08:54
dt: 1.999910 s
Dimension: [ 256. 256.] pix
scale: [ 117.59999847 117.59999847] arcsec / pix
array([[ 692., 692., 692., ..., 692., 692., 692.],
[ 692., 692., 692., ..., 694., 692., 692.],
[ 692., 692., 692., ..., 694., 692., 692.],
...,
[ 692., 692., 692., ..., 696., 692., 692.],
[ 692., 692., 692., ..., 694., 692., 692.],
[ 692., 692., 692., ..., 692., 692., 692.]], dtype=float32)
(6, '/Users/schriste/Desktop/sunpy_test_img/20110919_000901_s4h1A.fts')
SunPy HIMap
---------
Observatory: STEREO A
Instrument: SECCHI
Detector: HI1
Measurement: white-light
Wavelength: nan nm
Obs Date: 2011-09-19 00:09:01
dt: 1199.940000 s
Dimension: [ 1024. 1024.] pix
scale: [ 0.01995507 0.01995507] deg / pix
array([[ 9510., 10658., 10495., ..., 41466., 42677., 39874.],
[ 9209., 9464., 9556., ..., 38254., 38567., 37338.],
[ 9660., 9785., 9840., ..., 37093., 37483., 36754.],
...,
[ 9682., 9904., 9969., ..., 81340., 83771., 86069.],
[ 9784., 9993., 10083., ..., 81822., 84202., 86140.],
[ 10056., 10088., 9939., ..., 4682., 4307., 4541.]])
(7, '/Users/schriste/Desktop/sunpy_test_img/20110919_000901_s4h1B.fts')
SunPy HIMap
---------
Observatory: STEREO B
Instrument: SECCHI
Detector: HI1
Measurement: white-light
Wavelength: nan nm
Obs Date: 2011-09-19 00:09:54
dt: 1199.940000 s
Dimension: [ 1024. 1024.] pix
scale: [ 0.01992513 0.01992513] deg / pix
array([[ 5208., 6936., 5613., ..., 9135., 9128., 8999.],
[ 65603., 66324., 66282., ..., 9122., 9122., 8956.],
[ 65686., 67354., 66767., ..., 9189., 9056., 9061.],
...,
[ 30398., 31022., 31212., ..., 15168., 11138., 7629.],
[ 30672., 31394., 31125., ..., 18889., 11349., 7353.],
[ 33041., 34608., 34173., ..., 9082., 8543., 7186.]])
(8, '/Users/schriste/Desktop/sunpy_test_img/20110919_000921_s4h2A.fts')
SunPy HIMap
---------
Observatory: STEREO A
Instrument: SECCHI
Detector: HI2
Measurement: white-light
Wavelength: nan nm
Obs Date: 2011-09-19 00:09:21
dt: 4949.840000 s
Dimension: [ 1024. 1024.] pix
scale: [ 0.07209535 0.07209535] deg / pix
array([[ 28065., 30512., 30736., ..., 25526., 25433., 23280.],
[ 26350., 27075., 27195., ..., 23113., 23027., 22375.],
[ 26528., 27306., 27357., ..., 23014., 22980., 22353.],
...,
[ 28839., 29437., 29692., ..., 49586., 49126., 47941.],
[ 28685., 29286., 29355., ..., 49265., 48960., 47754.],
[ 28706., 29148., 29340., ..., 5758., 5420., 5733.]])
(9, '/Users/schriste/Desktop/sunpy_test_img/20110919_000921_s4h2B.fts')
SunPy HIMap
---------
Observatory: STEREO B
Instrument: SECCHI
Detector: HI2
Measurement: white-light
Wavelength: nan nm
Obs Date: 2011-09-19 00:10:14
dt: 4949.840000 s
Dimension: [ 1024. 1024.] pix
scale: [ 0.07189592 0.07189592] deg / pix
array([[ 5501., 5657., 5331., ..., 27144., 27014., 26857.],
[ 43007., 43429., 43351., ..., 27128., 27230., 26930.],
[ 43086., 43416., 43392., ..., 27102., 27216., 26852.],
...,
[ 22290., 22750., 22621., ..., 24993., 24808., 24547.],
[ 22371., 22839., 22928., ..., 25049., 25094., 24311.],
[ 23144., 25524., 25262., ..., 28280., 27973., 25612.]])
WARNING: VerifyWarning: Verification reported errors: [astropy.io.fits.verify]
WARNING:astropy:VerifyWarning: Verification reported errors:
WARNING: VerifyWarning: HDU 0: [astropy.io.fits.verify]
WARNING:astropy:VerifyWarning: HDU 0:
WARNING: VerifyWarning: Card 77: [astropy.io.fits.verify]
WARNING:astropy:VerifyWarning: Card 77:
WARNING: VerifyWarning: Unfixable error: Unprintable string 'offset_bias.pro\t1.30, 02/19/10, 398.366'; commentary cards may only contain printable ASCII characters [astropy.io.fits.verify]
WARNING:astropy:VerifyWarning: Unfixable error: Unprintable string 'offset_bias.pro\t1.30, 02/19/10, 398.366'; commentary cards may only contain printable ASCII characters
WARNING: VerifyWarning: Note: PyFITS uses zero-based indexing.
[astropy.io.fits.verify]
WARNING:astropy:VerifyWarning: Note: PyFITS uses zero-based indexing.
(10, '/Users/schriste/Desktop/sunpy_test_img/32232663.fts')
SunPy LASCOMap
---------
Observatory: SOHO
Instrument: LASCO
Detector: C3
Measurement: white-light
Wavelength: nan nm
Obs Date: 2010-09-04 17:30:05
dt: 19.092800 s
Dimension: [ 1024. 1024.] pix
scale: [ 56. 56.] arcsec / pix
array([[402, 407, 408, ..., 408, 407, 406],
[401, 405, 404, ..., 403, 404, 404],
[402, 404, 403, ..., 405, 405, 405],
...,
[402, 404, 405, ..., 404, 404, 405],
[418, 424, 422, ..., 422, 421, 423],
[420, 428, 428, ..., 426, 427, 428]], dtype=int16)
(11, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.131A_2013-09-21T16_00_08.62Z.image_lev1.fits')
SunPy AIAMap
---------
Observatory: SDO
Instrument: AIA 1
Detector: AIA
Measurement: 131.0 Angstrom
Wavelength: 131.0 Angstrom
Obs Date: 2013-09-21 16:00:08
dt: 2.900794 s
Dimension: [ 4096. 4096.] pix
scale: [ 0.600698 0.600698] arcsec / pix
array([[-1, 2, 3, ..., 0, 1, 1],
[ 0, -1, 2, ..., 0, -1, -1],
[ 0, 0, 1, ..., -2, -2, 0],
...,
[ 0, 0, 2, ..., 0, -1, 2],
[-3, 0, 2, ..., 2, 0, -2],
[ 2, 1, -2, ..., 0, -1, 0]], dtype=int16)
(12, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.171A_2013-09-21T16_00_11.34Z.image_lev1.fits')
SunPy AIAMap
---------
Observatory: SDO
Instrument: AIA 3
Detector: AIA
Measurement: 171.0 Angstrom
Wavelength: 171.0 Angstrom
Obs Date: 2013-09-21 16:00:11
dt: 2.000186 s
Dimension: [ 4096. 4096.] pix
scale: [ 0.599489 0.599489] arcsec / pix
array([[ 1, 0, 1, ..., -1, 1, 0],
[-1, 1, -1, ..., -1, 4, 0],
[ 1, 0, 0, ..., 0, 0, 1],
...,
[ 2, 2, -1, ..., 0, 0, 0],
[ 0, -1, -2, ..., 0, 1, -1],
[-2, 1, -1, ..., 0, 2, 0]], dtype=int16)
(13, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.193A_2013-09-21T16_00_06.84Z.image_lev1.fits')
SunPy AIAMap
---------
Observatory: SDO
Instrument: AIA 2
Detector: AIA
Measurement: 193.0 Angstrom
Wavelength: 193.0 Angstrom
Obs Date: 2013-09-21 16:00:06
dt: 1.999592 s
Dimension: [ 4096. 4096.] pix
scale: [ 0.600714 0.600714] arcsec / pix
array([[ 0, 2, -1, ..., 0, 0, 0],
[ 1, -1, -1, ..., -1, 0, -1],
[ 0, -1, 2, ..., -1, -1, 0],
...,
[ 1, 1, 2, ..., 1, 1, 0],
[ 0, 0, 0, ..., 0, 3, -1],
[ 1, 2, -3, ..., -1, -1, 1]], dtype=int16)
(14, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.211A_2013-09-21T16_00_11.62Z.image_lev1.fits')
SunPy AIAMap
---------
Observatory: SDO
Instrument: AIA 2
Detector: AIA
Measurement: 211.0 Angstrom
Wavelength: 211.0 Angstrom
Obs Date: 2013-09-21 16:00:11
dt: 2.901205 s
Dimension: [ 4096. 4096.] pix
scale: [ 0.600758 0.600758] arcsec / pix
array([[ 1, 0, 1, ..., 0, -1, 0],
[-2, 2, -1, ..., 1, 2, 0],
[-1, 3, 0, ..., -2, 0, 0],
...,
[-1, -2, -1, ..., 1, -1, 1],
[-1, 0, 1, ..., 0, 0, 0],
[-2, -1, 2, ..., 0, 0, -1]], dtype=int16)
(15, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.304A_2013-09-21T16_00_07.14Z.image_lev1.fits')
SunPy AIAMap
---------
Observatory: SDO
Instrument: AIA 4
Detector: AIA
Measurement: 304.0 Angstrom
Wavelength: 304.0 Angstrom
Obs Date: 2013-09-21 16:00:07
dt: 2.902062 s
Dimension: [ 4096. 4096.] pix
scale: [ 0.600165 0.600165] arcsec / pix
array([[ 0, 0, -1, ..., 1, 1, 0],
[ 1, 0, 1, ..., 0, -1, 1],
[ 0, 0, 2, ..., -2, 2, 0],
...,
[-1, 1, 1, ..., -1, 0, 2],
[ 0, -1, 1, ..., -2, 1, -2],
[ 0, 2, 0, ..., 0, 1, -1]], dtype=int16)
(16, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.335A_2013-09-21T16_00_02.63Z.image_lev1.fits')
SunPy AIAMap
---------
Observatory: SDO
Instrument: AIA 1
Detector: AIA
Measurement: 335.0 Angstrom
Wavelength: 335.0 Angstrom
Obs Date: 2013-09-21 16:00:02
dt: 2.901319 s
Dimension: [ 4096. 4096.] pix
scale: [ 0.600737 0.600737] arcsec / pix
array([[ 1, -2, 1, ..., 0, -1, -1],
[-1, 2, -1, ..., 1, 1, -2],
[-2, -2, 2, ..., -1, 0, 0],
...,
[ 0, 0, 2, ..., -1, 0, 0],
[ 0, 0, 1, ..., -1, 0, -1],
[-2, 1, 0, ..., 0, 0, 0]], dtype=int16)
(17, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.94A_2013-09-21T16_00_01.12Z.image_lev1.fits')
SunPy AIAMap
---------
Observatory: SDO
Instrument: AIA 4
Detector: AIA
Measurement: 94.0 Angstrom
Wavelength: 94.0 Angstrom
Obs Date: 2013-09-21 16:00:01
dt: 2.900842 s
Dimension: [ 4096. 4096.] pix
scale: [ 0.600109 0.600109] arcsec / pix
array([[ 1, 0, -1, ..., -1, 0, 1],
[ 0, 0, -2, ..., -3, 0, 0],
[ 1, 1, -1, ..., -1, -1, -2],
...,
[ 2, 0, -1, ..., 0, 0, -1],
[ 1, -2, 1, ..., 0, 0, -1],
[ 1, -3, 0, ..., -1, -2, 0]], dtype=int16)
(18, '/Users/schriste/Desktop/sunpy_test_img/hmi.Ic_45s_nrt.20130514_020000_TAI.2.continuum.fits')
SunPy HMIMap
---------
Observatory: SDO
Instrument: HMI FRONT2
Detector: HMI
Measurement: continuum
Wavelength: 6173.0 Angstrom
Obs Date: 2013-05-14 01:59:07
dt: 0.000000 s
Dimension: [ 4096. 4096.] pix
scale: [ 0.504283 0.504283] arcsec / pix
array([[ nan, nan, nan, ..., nan, nan, nan],
[ nan, nan, nan, ..., nan, nan, nan],
[ nan, nan, nan, ..., nan, nan, nan],
...,
[ nan, nan, nan, ..., nan, nan, nan],
[ nan, nan, nan, ..., nan, nan, nan],
[ nan, nan, nan, ..., nan, nan, nan]], dtype=float32)
(19, '/Users/schriste/Desktop/sunpy_test_img/hsi_image_20101016_191218.fits')
SunPy RHESSIMap
---------
Observatory: RHESSI
Instrument: RHESSI
Detector: RHESSI
Measurement: [ 12. 25.] keV
Wavelength: [ 12. 25.] keV
Obs Date: 2010-10-16 19:12:18
dt: 0.000000 s
Dimension: [ 64. 64.] pix
scale: [ 4. 4.] arcsec / pix
array([[ 0.01848893, 0.01848893, 0.08024438, ..., 0.03925416,
0.03616409, 0.01848893],
[ 0.06725851, 0.09361217, 0.08673389, ..., 0.05927928,
0.04436807, 0.01848893],
[ 0.08278352, 0.10089278, 0.09368125, ..., 0.0821826 ,
0.06657066, 0.01848893],
...,
[ 0.21122217, 0.20945911, 0.20059948, ..., 0.01409415,
0.05177701, 0.01848893],
[ 0.21778092, 0.21185483, 0.19914874, ..., 0.04582634,
0.08060541, 0.01848893],
[ 0.01848893, 0.01848893, 0.20828068, ..., 0.06199721,
0.1085261 , 0.01848893]], dtype=float32)
(20, '/Users/schriste/Desktop/sunpy_test_img/mdi.fd_M_lev182.2003.10.29_22_00_30_TAI.data.fits')
SunPy MDIMap
---------
Observatory: SOHO
Instrument:
Detector: MDI
Measurement: magnetogram
Wavelength: nan nm
Obs Date: 2003-10-29 21:59:43
dt: 0.000000 s
Dimension: [ 1024. 1024.] pix
scale: [ 1.983938 1.983938] arcsec / pix
array([[ nan, nan, nan, ..., nan, nan, nan],
[ nan, nan, nan, ..., nan, nan, nan],
[ nan, nan, nan, ..., nan, nan, nan],
...,
[ nan, nan, nan, ..., nan, nan, nan],
[ nan, nan, nan, ..., nan, nan, nan],
[ nan, nan, nan, ..., nan, nan, nan]])
(21, '/Users/schriste/Desktop/sunpy_test_img/swap_lv1_20120101_001607.fits')
SunPy SWAPMap
---------
Observatory: PROBA2
Instrument: SWAP
Detector: SWAP
Measurement: 174.0 Angstrom
Wavelength: 174.0 Angstrom
Obs Date: 2012-01-01 00:16:07
dt: 10.000000 s
Dimension: [ 1024. 1024.] pix
scale: [ 3.16226784 3.16226784] arcsec / pix
array([[ 0., 0., 0., ..., 0., 0., 0.],
[ 0., 0., 0., ..., 0., 0., 0.],
[ 0., 0., 0., ..., 0., 0., 0.],
...,
[ 0., 0., 0., ..., 0., 0., 0.],
[ 0., 0., 0., ..., 0., 0., 0.],
[ 0., 0., 0., ..., 0., 0., 0.]], dtype=float32)
(22, '/Users/schriste/Desktop/sunpy_test_img/XRT20130916_160004.3.fits')
SunPy XRTMap
---------
Observatory: Hinode
Instrument: XRT
Detector: XRT
Measurement: Open-Ti poly
Wavelength: nan keV
Obs Date: 2013-09-16 16:00:04
dt: 0.008156 s
Dimension: [ 256. 256.] pix
scale: [ 8.2288 8.2288] arcsec / pix
masked_array(data =
[[756 756 756 ..., 749 751 753]
[754 754 754 ..., 756 757 758]
[753 753 753 ..., 750 749 749]
...,
[750 750 750 ..., 750 750 750]
[750 750 750 ..., 750 750 750]
[750 750 750 ..., 750 750 750]],
mask =
False,
fill_value = 999999)
(23, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184106.2.fits')
SunPy XRTMap
---------
Observatory: Hinode
Instrument: XRT
Detector: XRT
Measurement: Open-Gband
Wavelength: nan keV
Obs Date: 2014-12-11 18:41:06
dt: 0.044840 s
Dimension: [ 384. 384.] pix
scale: [ 1.0286 1.0286] arcsec / pix
masked_array(data =
[[ 958 1010 1055 ..., 916 919 892]
[ 940 1000 1035 ..., 907 893 905]
[ 965 989 1002 ..., 932 914 905]
...,
[ 317 318 317 ..., 109 105 108]
[ 325 313 306 ..., 107 103 103]
[ 318 314 285 ..., 105 102 104]],
mask =
False,
fill_value = 999999)
(24, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184129.7.fits')
SunPy XRTMap
---------
Observatory: Hinode
Instrument: XRT
Detector: XRT
Measurement: Open-Gband
Wavelength: nan keV
Obs Date: 2014-12-11 18:41:29
dt: 0.032692 s
Dimension: [ 384. 384.] pix
scale: [ 1.03206 1.03206] arcsec / pix
masked_array(data =
[[1486 1515 1532 ..., 1372 1376 1348]
[1461 1492 1537 ..., 1358 1340 1357]
[1493 1496 1519 ..., 1373 1359 1357]
...,
[ 451 452 459 ..., 103 97 102]
[ 436 433 428 ..., 102 97 100]
[ 410 405 390 ..., 102 98 99]],
mask =
False,
fill_value = 999999)
(25, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184138.4.fits')
SunPy XRTMap
---------
Observatory: Hinode
Instrument: XRT
Detector: XRT
Measurement: Open-Be thick
Wavelength: nan keV
Obs Date: 2014-12-11 18:41:38
dt: 32.768000 s
Dimension: [ 512. 512.] pix
scale: [ 4.1144 4.1144] arcsec / pix
masked_array(data =
[[351 348 351 ..., 352 349 350]
[350 352 349 ..., 350 352 350]
[345 351 348 ..., 352 348 350]
...,
[348 348 343 ..., 348 345 347]
[347 341 347 ..., 346 345 343]
[346 346 349 ..., 345 431 351]],
mask =
False,
fill_value = 999999)
(26, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184221.9.fits')
SunPy XRTMap
---------
Observatory: Hinode
Instrument: XRT
Detector: XRT
Measurement: Open-Be thick
Wavelength: nan keV
Obs Date: 2014-12-11 18:42:21
dt: 32.768000 s
Dimension: [ 512. 512.] pix
scale: [ 4.1144 4.1144] arcsec / pix
masked_array(data =
[[352 351 351 ..., 347 350 351]
[348 352 349 ..., 352 351 349]
[352 350 348 ..., 349 353 352]
...,
[347 344 346 ..., 346 343 347]
[345 346 347 ..., 347 342 345]
[346 348 340 ..., 347 346 344]],
mask =
[[False False False ..., False False False]
[False False False ..., False False False]
[False False False ..., False False False]
...,
[False False False ..., False False False]
[False False False ..., False False False]
[False False False ..., False False False]],
fill_value = 999999)
(27, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184304.6.fits')
SunPy XRTMap
---------
Observatory: Hinode
Instrument: XRT
Detector: XRT
Measurement: Open-Al thick
Wavelength: nan keV
Obs Date: 2014-12-11 18:43:04
dt: 16.384500 s
Dimension: [ 512. 512.] pix
scale: [ 4.1144 4.1144] arcsec / pix
masked_array(data =
[[350 346 347 ..., 349 350 348]
[348 350 350 ..., 348 348 349]
[348 349 349 ..., 350 350 346]
...,
[343 344 345 ..., 343 343 344]
[345 344 346 ..., 344 344 343]
[345 344 346 ..., 345 346 345]],
mask =
False,
fill_value = 999999)
(28, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184324.0.fits')
SunPy XRTMap
---------
Observatory: Hinode
Instrument: XRT
Detector: XRT
Measurement: Open-Al thick
Wavelength: nan keV
Obs Date: 2014-12-11 18:43:24
dt: 16.384800 s
Dimension: [ 512. 512.] pix
scale: [ 4.1144 4.1144] arcsec / pix
masked_array(data =
[[351 349 349 ..., 347 352 347]
[349 345 352 ..., 349 351 349]
[350 349 347 ..., 350 346 349]
...,
[346 344 345 ..., 347 345 345]
[342 345 346 ..., 347 342 342]
[344 343 343 ..., 343 344 346]],
mask =
False,
fill_value = 999999)
(29, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184512.3.fits')
SunPy XRTMap
---------
Observatory: Hinode
Instrument: XRT
Detector: XRT
Measurement: Al poly-Ti poly
Wavelength: nan keV
Obs Date: 2014-12-11 18:45:12
dt: 0.064196 s
Dimension: [ 512. 512.] pix
scale: [ 4.1144 4.1144] arcsec / pix
masked_array(data =
[[348 349 346 ..., 348 349 347]
[347 349 349 ..., 348 347 345]
[346 348 349 ..., 347 349 348]
...,
[344 343 343 ..., 345 343 344]
[343 343 344 ..., 341 344 345]
[341 345 345 ..., 347 342 346]],
mask =
False,
fill_value = 999999)
(30, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184519.2.fits')
SunPy XRTMap
---------
Observatory: Hinode
Instrument: XRT
Detector: XRT
Measurement: Al poly-Ti poly
Wavelength: nan keV
Obs Date: 2014-12-11 18:45:19
dt: 0.724640 s
Dimension: [ 512. 512.] pix
scale: [ 4.1144 4.1144] arcsec / pix
masked_array(data =
[[347 347 348 ..., 349 345 350]
[345 348 345 ..., 347 347 347]
[347 349 347 ..., 348 346 347]
...,
[345 345 343 ..., 345 339 346]
[345 343 341 ..., 345 342 343]
[344 347 345 ..., 346 345 343]],
mask =
[[False False False ..., False False False]
[False False False ..., False False False]
[False False False ..., False False False]
...,
[False False False ..., False False False]
[False False False ..., False False False]
[False False False ..., False False False]],
fill_value = 999999)
In [3]:
i = 1
print(nplots)
fig = plt.figure(figsize=(10, nplots*10))
for ind, f in enumerate(mapfiles):
print(ind, f)
m = Map(f)
fig.add_subplot(nplots, 1, i)
if type(m) == type([]):
m[0].plot(title=m[0].name + ' ' + str(ind))
else:
m.plot(title=m.name + ' ' + str(ind))
plt.colorbar()
i += 1
plt.show()
31
(0, '/Users/schriste/Desktop/sunpy_test_img/20110919_000500_s4c1A.fts')
(1, '/Users/schriste/Desktop/sunpy_test_img/20110919_000500_s4c1B.fts')
(2, '/Users/schriste/Desktop/sunpy_test_img/20110919_000530_n4euA.fts')
(3, '/Users/schriste/Desktop/sunpy_test_img/20110919_000530_n4euB.fts')
(4, '/Users/schriste/Desktop/sunpy_test_img/20110919_000800_d4c2A.fts')
(5, '/Users/schriste/Desktop/sunpy_test_img/20110919_000800_d4c2B.fts')
(6, '/Users/schriste/Desktop/sunpy_test_img/20110919_000901_s4h1A.fts')
(7, '/Users/schriste/Desktop/sunpy_test_img/20110919_000901_s4h1B.fts')
(8, '/Users/schriste/Desktop/sunpy_test_img/20110919_000921_s4h2A.fts')
(9, '/Users/schriste/Desktop/sunpy_test_img/20110919_000921_s4h2B.fts')
(10, '/Users/schriste/Desktop/sunpy_test_img/32232663.fts')
(11, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.131A_2013-09-21T16_00_08.62Z.image_lev1.fits')
(12, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.171A_2013-09-21T16_00_11.34Z.image_lev1.fits')
(13, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.193A_2013-09-21T16_00_06.84Z.image_lev1.fits')
(14, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.211A_2013-09-21T16_00_11.62Z.image_lev1.fits')
(15, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.304A_2013-09-21T16_00_07.14Z.image_lev1.fits')
(16, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.335A_2013-09-21T16_00_02.63Z.image_lev1.fits')
(17, '/Users/schriste/Desktop/sunpy_test_img/aia.lev1.94A_2013-09-21T16_00_01.12Z.image_lev1.fits')
(18, '/Users/schriste/Desktop/sunpy_test_img/hmi.Ic_45s_nrt.20130514_020000_TAI.2.continuum.fits')
(19, '/Users/schriste/Desktop/sunpy_test_img/hsi_image_20101016_191218.fits')
(20, '/Users/schriste/Desktop/sunpy_test_img/mdi.fd_M_lev182.2003.10.29_22_00_30_TAI.data.fits')
(21, '/Users/schriste/Desktop/sunpy_test_img/swap_lv1_20120101_001607.fits')
(22, '/Users/schriste/Desktop/sunpy_test_img/XRT20130916_160004.3.fits')
(23, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184106.2.fits')
(24, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184129.7.fits')
(25, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184138.4.fits')
(26, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184221.9.fits')
(27, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184304.6.fits')
(28, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184324.0.fits')
(29, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184512.3.fits')
(30, '/Users/schriste/Desktop/sunpy_test_img/XRT20141211_184519.2.fits')
/Users/schriste/Developer/repositories/sunpy/sunpy/map/mapbase.py:1257: Warning: This map is not properly oriented. Plot axes may be incorrect
Warning)
In [4]:
i = 1
for ind, f in enumerate(mapfiles):
m = Map(f)
norm = m.plot_settings['norm']
try:
print("{i} {instr} {norm} {mi} {ma} {gamma}".format(i=ind, instr=m.name[0:20], norm=norm, mi=norm.vmin, ma=norm.vmax, gamma=norm.gamma))
except:
print("{i} {instr} {norm} {mi} {ma}".format(i=ind, instr=m.name[0:20], norm=norm, mi=norm.vmin, ma=norm.vmax))
i += 1
0 STEREO A COR1 white- <matplotlib.colors.PowerNorm object at 0x12ed24690> None None 0.5
1 STEREO B COR1 white- <matplotlib.colors.PowerNorm object at 0x12dfa9810> None None 0.5
2 STEREO A EUVI 195.0 <matplotlib.colors.PowerNorm object at 0x12df93a50> 686.0 16398.0 0.25
3 STEREO B EUVI 195.0 <matplotlib.colors.PowerNorm object at 0x12dfa9810> 654.0 16404.0 0.25
4 STEREO A COR2 white- <matplotlib.colors.PowerNorm object at 0x12df93b90> None None 0.5
5 STEREO B COR2 white- <matplotlib.colors.PowerNorm object at 0x149e42110> None None 0.5
6 STEREO A HI1 white-l <matplotlib.colors.PowerNorm object at 0x13ab8ab90> None None 0.25
7 STEREO B HI1 white-l <matplotlib.colors.PowerNorm object at 0x12ff538d0> None None 0.25
8 STEREO A HI2 white-l <matplotlib.colors.PowerNorm object at 0x135851250> None None 0.25
9 STEREO B HI2 white-l <matplotlib.colors.PowerNorm object at 0x14ad35a90> None None 0.25
10 SOHO C3 white-light <matplotlib.colors.PowerNorm object at 0x12ec34dd0> None None 0.5
11 SDO AIA 131.0 Angstr <astropy.visualization.mpl_normalize.ImageNormalize object at 0x14a701250> None None
12 SDO AIA 171.0 Angstr <astropy.visualization.mpl_normalize.ImageNormalize object at 0x14b2eb250> None None
13 SDO AIA 193.0 Angstr <astropy.visualization.mpl_normalize.ImageNormalize object at 0x1473a71d0> None None
14 SDO AIA 211.0 Angstr <astropy.visualization.mpl_normalize.ImageNormalize object at 0x12dfa9a50> None None
15 SDO AIA 304.0 Angstr <astropy.visualization.mpl_normalize.ImageNormalize object at 0x12fe373d0> None None
16 SDO AIA 335.0 Angstr <astropy.visualization.mpl_normalize.ImageNormalize object at 0x135851150> None None
17 SDO AIA 94.0 Angstro <astropy.visualization.mpl_normalize.ImageNormalize object at 0x12eccab50> None None
18 SDO HMI continuum 20 <astropy.visualization.mpl_normalize.ImageNormalize object at 0x1473c32d0> None None
19 RHESSI RHESSI [ 12. <astropy.visualization.mpl_normalize.ImageNormalize object at 0x11ddb5750> None None
20 SOHO MDI magnetogram <matplotlib.colors.Normalize object at 0x14a0ddc50> 3512.11 -3512.11
21 PROBA2 SWAP 174.0 An <astropy.visualization.mpl_normalize.ImageNormalize object at 0x11576c790> None None
22 Hinode XRT Open-Ti p <astropy.visualization.mpl_normalize.ImageNormalize object at 0x12dba1c90> None None
23 Hinode XRT Open-Gban <astropy.visualization.mpl_normalize.ImageNormalize object at 0x146f70050> None None
24 Hinode XRT Open-Gban <astropy.visualization.mpl_normalize.ImageNormalize object at 0x12dbb7c10> None None
25 Hinode XRT Open-Be t <astropy.visualization.mpl_normalize.ImageNormalize object at 0x146f70050> None None
26 Hinode XRT Open-Be t <astropy.visualization.mpl_normalize.ImageNormalize object at 0x12dba1410> None None
27 Hinode XRT Open-Al t <astropy.visualization.mpl_normalize.ImageNormalize object at 0x146f70050> None None
28 Hinode XRT Open-Al t <astropy.visualization.mpl_normalize.ImageNormalize object at 0x12dba1410> None None
29 Hinode XRT Al poly-T <astropy.visualization.mpl_normalize.ImageNormalize object at 0x146f70050> None None
30 Hinode XRT Al poly-T <astropy.visualization.mpl_normalize.ImageNormalize object at 0x12df4a050> None None
In [5]:
from sunpy import cm
In [6]:
cm.get_cmap('stereocor{det!s}'.format(det=m1.detector[-1]))
---------------------------------------------------------------------------
NameError Traceback (most recent call last)
<ipython-input-6-bc16533e2b1e> in <module>()
----> 1 cm.get_cmap('stereocor{det!s}'.format(det=m1.detector[-1]))
NameError: name 'm1' is not defined
In [ ]:
In [ ]:
nplots = 2
fig = plt.figure(figsize=(10, nplots*10))
fig.add_subplot(2, 1, 1)
m1.plot()
plt.colorbar()
fig.add_subplot(2, 1, 2)
m1.plot(norm=colors.PowerNorm(0.1))
plt.colorbar()
In [ ]:
m1 = Map(mapfiles[27])
In [ ]:
plt.imshow(np.arange(0,100).reshape(10, 10), cmap=m1.plot_settings['cmap'])
In [ ]:
np.arange(0,100).reshape(10, 10)
In [ ]:
import numpy.ma
In [ ]:
m2.data = numpy.ma.masked_where(m2.data == np.nan, m2.data)
In [ ]:
m2.data
In [ ]:
plt.imshow(m2.data)
In [ ]:
m1.meta.get('content').find('Magneto')
In [ ]:
m1.plot_settings
In [ ]:
file = '/Users/schriste/Desktop/sunpy_test_img/iris_l2_20130801_074720_4040000014_SJI_1400_t000.fits'
In [ ]:
from sunpy.instr.iris import SJI_to_cube
In [ ]:
iris = SJI_to_cube(file)
In [ ]:
type(iris.maps[0])
In [ ]:
iris[0].data.shape
In [ ]:
Content source: ehsteve/ipython-notebooks
Similar notebooks: