<jupyter_start><jupyter_text>Giant Magellan Telescope Images<jupyter_code>import numpy as np
import ceo
%pylab inline<jupyter_output>Populating the interactive namespace from numpy and matplotlib<jupyter_text>Diffraction limited image The point source:<jupyter_code>n = 512
src = ceo.Source("K",resolution=(n+1,n+1),magnitude=14)<jupyter_output><empty_output><jupyter_text>The telescope:<jupyter_code>tel = ceo.GMT(n+1,26.0)
src.masked(tel)<jupyter_output><empty_output><jupyter_text>The imager sampling the GMT image at 4 times the Shannon criteria:<jupyter_code>osf = 8
imgr = ceo.Imaging(1,n,DFT_osf=osf,N_PX_IMAGE=128)
imgr.propagate(src)
frame = imgr.frame.host()
fig,(ax1,ax2) = plt.subplots(ncols=2)
fig.set_size_inches(14,6)
h1 = ax1.imshow(frame,interpolation='None') 
colorbar(h1,ax=ax1)
h2 = ax2.imshow(frame,interpolation='None',vmax=np.max(frame)*0.125) 
colorbar(h2,ax=ax2)
imshow(tel.f.host(shape=((n+1),(n+1))))
P = tel.piston_mask
print P[0][0].shape
imshow(np.reshape(P[0][6],((n+1),(n+1))))
colorbar()
src.reset()
tel.set_gmt_piston(src.wavefront.phase,np.array(np.random.randn(8)*200e-9,dtype=np.single))
imshow(src.phase.host())
colorbar()
print np.arange(1,8,dtype=np.single)
print tel.get_gmt_piston(src.wavefront.phase)*1e9
src.masked(tel)
imgr.reset()
imgr.propagate(src)
frame = imgr.frame.host()
fig,(ax1,ax2) = plt.subplots(ncols=2)
fig.set_size_inches(14,6)
h1 = ax1.imshow(frame,interpolation='None') 
colorbar(h1,ax=ax1)
h2 = ax2.imshow(frame,interpolation='None',vmax=np.max(frame)*0.125) 
colorbar(h2,ax=ax2)<jupyter_output><empty_output><jupyter_text>Turbulence images<jupyter_code>gs = ceo.Source("K",resolution=(n+1,n+1),magnitude=14)
gs.masked(tel)
atm = ceo.Atmosphere(15e-2,30,wind_speed=10)
imgr1 = ceo.Imaging(1,n,N_PX_IMAGE=128)
p = 25.5/n
print p*1e2<jupyter_output>4.98046875<jupyter_text>Snapshot image at T=1s<jupyter_code>imgr1.reset()
atm.get_phase_screen(gs,p,n+1,p,n+1,1)
imgr1.propagate(gs)
ax = imshow(imgr1.frame.host(),interpolation='None') 
colorbar(ax)<jupyter_output><empty_output><jupyter_text>Snapshot image at T=2s<jupyter_code>imgr1.reset()
atm.get_phase_screen(gs,p,n+1,p,n+1,2)
imgr1.propagate(gs)
ax = imshow(imgr1.frame.host(),interpolation='None') 
colorbar(ax)<jupyter_output><empty_output><jupyter_text>Snapshot image at T=10s<jupyter_code>imgr1.reset()
atm.get_phase_screen(gs,p,n+1,p,n+1,10)
imgr1.propagate(gs)
ax = imshow(imgr1.frame.host(),interpolation='None') 
colorbar(ax)<jupyter_output><empty_output><jupyter_text>20s exposure image sampled every 0.1s<jupyter_code>imgr1.reset()
tau = 1e-1
T = 0.0
tid = ceo.StopWatch()
tid.tic()
while T<=20:
    T += tau
    atm.get_phase_screen(gs,p,n+1,p,n+1,T)
    imgr1.propagate(gs)
tid.toc()
ax = imshow(imgr1.frame.host(),interpolation='None') 
colorbar(ax)<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.17206982543640897, 'fraction_numerical': 0.04916280726754542, 'mean_word_length': 9.51685393258427, 'pattern_counts': {'":': 0, '<': 36, '<?xml version=': 0, '>': 35, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '9', 'dir': 'jupyter-structured-clean-dedup', 'id': '380927', 'license': 'Zlib', 'n_tokens_mistral': 1343, 'n_tokens_neox': 1281, 'n_words': 141, 'path': 'imaging/GmtPsf.ipynb', 'repo_name': 'bsitarski/CEO', 'size': '2807'}