-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathplot_image.py
executable file
·110 lines (90 loc) · 27 KB
/
plot_image.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
#! /home/shulevski/miniconda2/envs/py27base/bin/python
def plot_image():
import matplotlib.pyplot as plt
from pylab import *
import aplpy
import numpy as np
from astropy.coordinates import SkyCoord
from astropy.coordinates import FK5
from astropy import units as u
from astropy.io import fits
#path = '/home/shulevski/Desktop/'
#path = '/home/shulevski/Documents/Research/Projects/A12_LBA/03Mar19/8min_images/'
#path = '/home/shulevski/Documents/Research/Projects/A12_source_counts/images_cat_41MHz/'
path = '/home/shulevski/Documents/Research/Projects/Image_combiner/'
#lofar_im = 'A12_4ch_8SBs_600s_high-MFS-image.fits'
#lofar_im = 'A12_03Mar19_MFS_I_masked.fits'
#lofar_im = 'A12_03Mar19_PBcorr_scaled_8min-MFS-I-image-pb_ZEA.fits'
#lofar_im = 'A12_18Feb19_PBcorr_MASKED_scaled-MFS-image-pb.fits'
#lofar_im = 'A12_rms.fits'
#lofar_im = 'A12_gresid_UVcut.fits'
#lofar_im='A12_18Feb19_MASKED_UVcut-MFS-image-pb.fits'
#lofar_im='A12_100_20_box_rms_UVcut.fits'
#lofar_im='A12_PBcorr_41MHz_MASKED_ZEA_UVcut-MFS-rms.fits'
#lofar_im='A12_PBcorr_41MHz_MASKED_UVcut-MFS-image-pb.fits'
#lofar_im='A12_PBcorr_41MHz_MASKED_UVcut-MFS-rms.fits'
#lofar_im = 'A12_18Feb19_MASKED-MFS_full-image-pb.fits'
#lofar_im = 'A12_18Feb19_MASKED-MFS-image-pb.fits'
#lofar_im = 'A12_18Feb19_PBcorr_scaled-MFS-image.fits'
lofar_im = 'rd_reprojected_crop.fits'
fig = plt.figure(figsize=(11, 9))
im = aplpy.FITSFigure(path+lofar_im, figure=fig, dimensions=[0, 1], slices=[0, 0], subplot=[0.15, 0.1, 0.8, 0.8])
#im.recenter(178.76, 52.84, radius=32.)
#im.show_grayscale(invert=True, stretch='linear', vmin=0.2, vmax=6.)
im.show_colorscale(cmap='inferno', stretch='linear', pmin=45., pmax=99.75)
#im.add_grid()
#im.grid.set_color('green')
'''
im.add_beam()
im.beam.set_edgecolor('white')
im.beam.set_facecolor('black')
im.beam.set_hatch('/')
#im.set_theme('publication')
im.tick_labels.set_xformat('hh:mm:ss')
im.tick_labels.set_yformat('dd:mm')
im.set_tick_color('k')
im.set_nan_color('white')
im.axis_labels.set_font(size=18)
im.tick_labels.set_font(size=16)
'''
'''
im.add_colorbar()
im.colorbar.show()
im.colorbar.set_location('right')
im.colorbar.set_width(0.1)
im.colorbar.set_pad(0.03)
#im.colorbar.set_axis_label_text('r.m.s. image noise[Jy/PSF]')
im.colorbar.set_axis_label_text('Flux [Jy/PSF]')
im.colorbar.set_axis_label_font(size=20)
im.colorbar.set_axis_label_pad(20)
im.colorbar.set_font(size=16)
'''
'''
factor = 1.
sigma = 1.
levno = 10
levbase = 2.
#levels = np.power(np.sqrt(levbase), range(levno))
#levels = np.insert(levels, 0, -levels[0], axis=0)
levels = np.linspace(0.2, 6., 10)
print "Levels: ", factor * sigma * levels
im.show_contour(path+lofar_im, levels=factor * sigma * levels, colors=['black'], linewidths=1.)
'''
'''
with fits.open(path+lofar_im) as fimage:
zra, zdec = fimage[0].header['CRVAL1'], fimage[0].header['CRVAL2']
if zra < 0.:
zra = zra + 360.
print "zra, zdec ", zra, zdec
'''
#im.show_circles(np.array([zra, zra, zra, zra]), np.array([zdec, zdec, zdec, zdec]), np.array([40., 50., 60., 70.]), edgecolor='red', linewidth=2)
#im.show_lines([np.array([[359.1028872436099, 354.542316779902, 350.0024768612484, 345.5032639675567, 341.0628825389106, 336.6971964812759, 332.41928100961366, 328.2391889876082, 324.16391814485286, 320.1975452313687, 316.3414823553622, 312.59480884115317, 308.95463665340543, 305.416475915222, 301.97457675457395, 298.62223283461714, 295.35203945369096, 292.15610469399223, 289.0262158424218, 285.95396551687514, 282.93084299303223, 279.9482965078879, 276.9977721216405, 274.0707342797265, 271.1586726942086, 268.2530996658346, 265.34554156170327, 262.427527887253, 259.49058126504264, 256.5262116633498, 253.52591840133792, 250.48120377953094, 247.3836026143626, 244.2247324407258, 240.996369600644, 237.69055672725432, 234.29974706981528, 230.81699042770404, 227.23616384416596, 223.5522472897494, 219.76164000303464, 215.86250675875363, 211.8551352258102, 207.74227641389118, 203.52943135185242, 199.2250407105575, 194.84053270567557, 190.39019080707587, 185.8908179745569, 181.36119765524774, 176.82138022556842, 172.29185118717197, 167.792657594308, 163.34257657251138, 158.9584023496056, 154.65440801094235, 150.44201054806683, 146.3296393317048, 142.32278465205962, 138.42418781176917, 134.63412809911534, 130.950763371205, 127.37048742172428, 123.88827616821509, 120.49800385564926, 117.19271857481314, 113.96487278429039, 110.80650908017465, 107.70940437259537, 104.665177240884, 101.66536391382195, 98.70146838287779, 95.76499186408176, 92.84744636936668, 89.9403566627715, 87.03525444678395, 84.12366830209423, 81.19711272029011, 78.24707953841157, 75.26503521044415, 72.24242762694345, 69.17070660030976, 66.04136263122287, 62.84598909473626, 59.576373425295266, 56.22462307632334, 52.78333175172989, 49.245790349321375, 45.606244852470354, 41.86019966978323, 38.00475933830903, 34.03899398424741, 29.964304819776103, 25.784756247021928, 21.507332647861453, 17.142073205229927, 12.702039974201629, 8.203085191210835, 3.6634040896590228, 359.1028872436099], [37.184887886874016, 37.098091830271905, 36.838142355602386, 36.4078277268027, 35.811653518455266, 35.05564163565997, 34.14707938445124, 33.09424230034684, 31.906113600063833, 30.592119441565078, 29.161893892875376, 27.62508185298765, 25.99118308381889, 24.269436543498994, 22.468741527782157, 20.59761063116952, 18.664148992487743, 16.67605440919953, 14.640633431252617, 12.5648292729179, 10.455258164295943, 8.318251511091187, 6.1599018913773955, 3.9861114695179474, 1.802641845962587, -0.38483530746494216, -2.5706872450590823, -4.749269053588698, -6.914872684753068, -9.061676352293098, -11.183694634493623, -13.274729810198423, -15.32832523621545, -17.337721954378292, -19.29582020398952, -21.19514810987165, -23.027840504488104, -24.785631588672597, -26.45986586949419, -28.041532419624634, -29.521327808942356, -30.88975283944784, -32.13724720527855, -33.25436414622262, -34.231983900668666, -35.06156032630963, -35.73538979019655, -36.24688606023385, -36.59084052478812, -36.76364483795121, -36.76345405108761, -36.590272873313836, -36.245955462521685, -35.73411876847958, -35.05997906534358, -34.23012905872976, -33.252277519507125, -32.13497434453578, -30.887341703594096, -29.518827519339418, -28.038992155613965, -26.45733392453941, -24.783154589399814, -23.02546278286457, -21.19291121148204, -19.293762511997535, -17.335878408232308, -15.326727125454653, -13.27340462601598, -11.18266596441421, -9.060963806035, -6.914491841848386, -4.7492314255277615, -2.571000246588252, -0.3855022332340298, 1.8016218418791103, 3.9847433925342375, 6.1581949239649205, 8.316219028128376, 10.452917744133696, 12.562202700725106, 14.637746689570623, 16.672937648666327, 18.660836473339074, 20.59414062337665, 22.46515615459307, 24.265781549865356, 25.98750750801939, 27.621437577231916, 29.15833504761088, 30.5887016397541, 31.902892983220692, 33.09127438897389, 34.14441774255762, 35.05333638210156, 35.80974974506525, 36.4063639109073, 36.83714887700858, 37.09758971392485, 37.184887886874016]])], edgecolor='black', linestyle='dashed', linewidth=2) # 5 degree elevation
#im.show_lines([np.array([[359.1039777318443, 353.8608692405084, 348.6532604679974, 343.5146866371425, 338.47484789756845, 333.55831062024424, 328.7837932299857, 324.16402374541275, 319.70606901699546, 315.4119903122779, 311.2796747794273, 307.30371509869644, 303.4762453162519, 299.78767739657445, 296.2273133800252, 292.78382964176024, 289.4456429997903, 286.20117507528886, 283.0390334297348, 279.948127388421, 276.91773444377617, 273.9375305970833, 270.99759546545164, 268.08840072193584, 265.2007885600375, 262.32594539653553, 259.4553749179386, 256.5808737842897, 253.69451277173638, 250.78862580466136, 247.8558091434963, 244.88893290037223, 241.88116699380976, 238.8260235622421, 235.71741766196934, 232.549747695582, 229.3179963607658, 226.01785188343746, 222.64584782171187, 219.19951775069865, 215.67755868052276, 212.07999424246185, 208.40832576103475, 204.6656567332211, 200.85677452529393, 196.98817291666853, 193.06800105989595, 189.10592884305183, 185.11292546563695, 181.1009566440298, 177.08261506596833, 173.0707069665631, 169.0778234342709, 165.11592710583997, 161.1959828428887, 157.32765523438107, 153.51908750936585, 149.776767240387, 146.1054756180306, 142.5083102607749, 138.9867671170597, 135.54086508595995, 132.16929717113558, 128.8695936981226, 125.63828572395016, 122.47105968543916, 119.36289714902784, 116.30819598081558, 113.30087123000848, 110.33443549443749, 107.40205956105677, 104.49661476904048, 101.61069892153036, 98.7366477654123, 95.86653414922483, 92.99215702965755, 90.10502259051064, 87.19631992203695, 84.25689403847421, 81.27721954297117, 78.24737903953248, 75.15705149886595, 71.99551726070997, 68.75168823020448, 65.4141740833883, 61.97139782796054, 58.41177660134796, 54.72398560833616, 50.89732372470843, 46.92219718817942, 42.79073116800708, 38.497505787809686, 34.04039160722179, 29.421429254714777, 24.647661378332458, 19.73178936847242, 14.692504436018455, 9.554347586671865, 4.346998175169605, 359.1039777318443], [46.98429765041916, 46.8845649500192, 46.58618704157542, 46.0935932374687, 45.41384844326593, 44.55622016043804, 43.53166143014613, 42.352272179712386, 41.030792096808675, 39.58016196351833, 38.013172715429945, 36.34220633689888, 34.57906200509963, 32.734854846024994, 30.819972416441093, 28.844074362609767, 26.816122492958606, 24.744430878501554, 22.636728009778857, 20.500225194241285, 18.34168715142469, 16.16750213908655, 13.983749960203664, 11.796266917896181, 9.610707264659316, 7.4326009883663, 5.267407934481895, 3.1205683151817314, 0.9975496255131002, -1.0961101093363308, -3.1547629824339842, -5.172616174264682, -7.143701907640238, -9.061853327854898, -10.920687546354477, -12.71359739255274, -14.433753706569165, -16.074120237903365, -17.627483340162307, -19.086498605435473, -20.443756289944208, -21.691866771043568, -22.82356628547372, -23.83184180749093, -24.710072172331472, -25.45218055873026, -26.05279143228587, -26.507383321406508, -26.812427694185256, -26.965504045489094, -26.965382293690205, -26.812065741487004, -26.506790963752447, -26.051984630843467, -25.451180920091435, -24.70890627332738, -23.830540408541008, -22.822163478688715, -21.690399106849622, -20.442261919933177, -19.085016470933972, -17.626052427079692, -16.07277890991875, -14.432539121424295, -12.712545000507825, -10.919830665804087, -9.061222790537366, -7.143325764538973, -5.172519452887219, -3.154967486376039, -1.0966342516484753, 0.9966909628262337, 3.119363902319118, 5.265850301583978, 7.430686524946501, 9.608436315100311, 11.793643876687696, 13.980783369912324, 16.164204793303455, 18.338076204041375, 20.49632227523857, 22.63255934611521, 24.74002741645871, 26.811520015193, 28.83931358889283, 30.815099071353494, 32.72991966947332, 34.57412067238945, 36.337319244553186, 38.00840457797537, 39.575581157085004, 41.026469680963736, 42.34828053425941, 43.52807246438456, 44.55310309372503, 45.41126721949503, 46.09160371735698, 46.58483420102738, 46.88388037699779, 46.98429765041916]])], edgecolor='black', linestyle='dashed', linewidth=2) # 10 degree elevation
#im.show_lines([np.array([[359.10585995271214, 352.84667334705415, 346.6568643904986, 340.60000696485247, 334.72892145041953, 329.08287149072646, 323.6868891695708, 318.552842759423, 313.68159290980714, 309.065582551552, 304.69136732738264, 300.54179912062483, 296.59774933360666, 292.83937360710087, 289.2469795904991, 285.8015802929423, 282.4852139568162, 279.28109948260453, 276.1736816083599, 273.1486060876515, 270.19265352340767, 267.2936515847315, 264.4403787868106, 261.62246839229755, 258.8303178321522, 256.05500694747127, 253.28822700787546, 250.52222163079375, 247.74974023658473, 244.96400440118578, 242.1586873187434, 239.32790649200885, 236.46622967353383, 233.56869394062943, 230.63083756293474, 227.6487439804581, 224.6190967267842, 221.53924349297492, 218.40726673604715, 215.22205732005074, 211.9833866979664, 208.6919721964474, 205.34952918786294, 201.95880348945198, 198.52357739218624, 195.04864344621197, 191.53974160990396, 188.00345760012365, 184.44708312691856, 180.8784418846545, 177.30568832039216, 173.7370888744206, 170.18079718361753, 166.64463536966727, 163.13589289794484, 159.6611526930453, 156.22615152016897, 152.83567849177516, 149.4935123708355, 146.2023954983491, 142.96403994499997, 139.7791600097175, 136.6475244664706, 133.5680219001773, 130.5387329189282, 127.55700380829343, 124.61951714008555, 121.72235582770824, 118.86105803573142, 116.03066114244591, 113.22573359272657, 110.44039396110134, 107.6683168853546, 104.90272575459927, 102.1363721753505, 99.36150233340703, 96.56981046347772, 93.75237978705171, 90.89961155132626, 88.00114329007074, 85.04575825658306, 82.02128932202116, 78.91452272481217, 75.71111021156372, 72.39550272400743, 68.95092532362033, 65.3594219642091, 61.60201029711459, 57.65900064445487, 53.51054810401117, 49.13751868868088, 44.522752063280336, 39.65278265893, 34.52002124527297, 29.12528427749866, 23.48038449480399, 17.610290498030846, 11.554201049253754, 5.364879405801936, 359.10585995271214], [56.98371847517952, 56.86470564684296, 56.50946970313317, 55.92593957490469, 55.12646127305312, 54.1266868349855, 52.9443524362216, 51.59813636103765, 50.10672076491434, 48.48810727128833, 46.759179052673694, 44.93546918295423, 43.03108354056316, 41.05872882199599, 39.029805085601005, 36.95453274783918, 34.84209346781601, 32.70077188189072, 30.538090603169326, 28.360934595726967, 26.175663391804882, 23.98821103145081, 21.80417437792697, 19.62889082493233, 17.46750651927553, 15.325036177998815, 13.206415447454932, 11.116546573500592, 9.060337950412825, 7.0427379049143255, 5.0687628586532965, 3.1435198028525537, 1.2722228174167434, -0.5397968204207079, -2.287087557351495, -3.964082264622628, -5.56511154748562, -7.084426503600414, -8.516231690147489, -9.854728865689605, -11.094171746587161, -12.228931547719979, -13.253572465895637, -14.162935533892501, -14.952228470336424, -15.617118348258058, -16.153823198296617, -16.55919815850256, -16.830811584956983, -16.967006725012965, -16.966945162388992, -16.830629246164765, -16.55890202699035, -16.15342470408374, -15.616633066462692, -14.951675722014102, -14.162337909190908, -13.25295529165288, -12.228322321865068, -11.093599560287963, -9.854223831537784, -8.515824391244356, -7.084147474667098, -5.564990795449597, -3.964148832807471, -2.2873691345824625, -0.5403193941231439, 1.2714352667169473, 3.1424455693787543, 5.067382746075712, 7.0410354328211655, 9.058299537688745, 11.114161703820368, 13.203676822112445, 15.321939861295276, 17.46405208172912, 19.62508148966292, 21.800017175724403, 23.983716970231843, 26.170847646360354, 28.355816723765717, 30.53269479352563, 32.69512724012989, 34.83623434294868, 36.94849910084731, 39.023642895102924, 41.05249051340915, 43.02482841095174, 44.929263772498636, 46.753097383613735, 48.48223083561384, 50.101138095270095, 51.59294197337802, 52.939644905125974, 54.12256578207924, 55.1230231522353, 55.92327235414679, 56.507646864377406, 56.863780278896286, 56.98371847517952]])], edgecolor='black', linestyle='dashed', linewidth=2) # 20 degree elevation
#im.show_lines([np.array([[359.10950380766013, 351.08507571257894, 343.23178209963254, 335.69559174333074, 328.58005447246194, 321.94190696646183, 315.79685362329616, 310.13040569355184, 304.9092124950545, 300.0903431599443, 295.62779165035414, 291.47649572815135, 287.5945101864064, 283.94396253684727, 280.49127684295115, 277.2069970616245, 274.0654180076329, 271.0441462253511, 268.12365793318395, 265.28688777665275, 262.5188628131309, 259.80638555298873, 257.13776445450696, 254.5025878028497, 251.89153606513668, 249.29622782371047, 246.70909480176368, 244.12328206101054, 241.5325700358894, 238.93131560343832, 236.3144098386153, 233.67725045649345, 231.01572719021408, 228.3262184966204, 225.60559802301725, 222.85124921434684, 220.06108629989072, 217.23357968743554, 214.3677835327623, 211.4633629738908, 208.52061826228228, 205.54050283469024, 202.5246323021057, 199.47528143830064, 196.3953665738574, 193.28841136895554, 190.1584947499713, 187.0101808180852, 183.848431703725, 180.6785055478501, 177.50584291766492, 174.33594588396855, 171.17425458946417, 168.0260263464473, 164.89622209159464, 161.78940442206587, 158.70965051649628, 155.6604821173291, 152.6448135437255, 149.6649175390642, 146.72240773464773, 143.8182357004137, 140.95269998707442, 138.1254642417149, 135.3355813736485, 132.58152081525333, 129.86119611148996, 127.17199032933927, 124.51077705698545, 121.87393502270317, 119.25735457458377, 116.65643440238924, 114.0660669369711, 111.48061082118839, 108.89384870329845, 106.29892835646908, 103.6882847765792, 101.05354046003647, 98.38538052883118, 95.67339878521031, 92.90591021241565, 90.06972502463846, 87.1498793565915, 84.12931851600196, 80.98853117944167, 77.70513832159348, 74.25345123580205, 70.60403226629528, 66.72332522173521, 62.57347743974141, 58.11256114167831, 53.295524792794325, 48.07635942453447, 42.41210762906885, 36.26935607942988, 29.633505150370397, 22.520094914956413, 14.985655595006138, 7.133523828859802, 359.10950380766013], [66.9831800735002, 66.83061996028037, 66.37785640390457, 65.6424258788821, 64.65020851853669, 63.431692111447205, 62.01855601499868, 60.441185689920346, 58.72722628655674, 56.900980120090736, 54.98336294485519, 52.992168874251064, 50.94246914113227, 48.84704005399416, 46.716765887551546, 44.56099365914462, 42.387834019466894, 40.20441099536962, 38.017066826890705, 35.8315290199227, 33.65304636856802, 31.486499832795005, 29.33649316384373, 27.20742722504361, 25.103561126295844, 23.029062591340903, 20.988049399704103, 18.984623273024592, 17.022897190144057, 15.10701680062231, 13.241176349216936, 11.429629314823051, 9.676693800371027, 7.986752582361631, 6.364247640019985, 4.813668936739903, 3.339537224529799, 1.9463806903886594, 0.6387053665269131, -0.579040612898663, -1.7025086075317453, -2.727496347155053, -3.6499992046728997, -4.466264647262222, -5.172848316323301, -5.766669891884776, -6.245066673794978, -6.6058426953707565, -6.847311204250871, -6.968328517749934, -6.968317587932243, -6.8472800785225765, -6.605796327392689, -6.245013196221577, -5.766620436652885, -5.172816773507576, -4.466267371921996, -3.6500546824917035, -2.7276248302966684, -1.702731732042091, -0.579381008273162, 0.6382244661721718, 1.9457358255972417, 3.3387050739411333, 4.812626660055808, 6.362973197462884, 7.985225029018428, 9.67489355505843, 11.427538404150907, 13.238778628379485, 15.10429815314102, 17.019845709359963, 18.981229428637537, 20.984306193131378, 23.024965704169468, 25.09910906668691, 27.20262147767194, 29.33133834988469, 31.4810038876776, 33.64722074677233, 35.82538893810047, 38.01063155910786, 40.19770423736537, 42.38088434247898, 44.55383507788236, 46.70943857077627, 48.83959120868496, 50.93495409990915, 52.984652405919405, 54.975920781278866, 56.89370063527041, 58.72021215515743, 60.43455516718125, 62.01244319991607, 63.42624520314153, 64.64558477931449, 65.63878208045907, 66.37533440646271, 66.82932903292038, 66.9831800735002]])], edgecolor='black', linestyle='dashed', linewidth=2) # 30 degree elevation
#im.show_lines([np.array([[359.37816974724853, 292.38959416697, 279.15361244103434, 273.45644112199915, 269.8386366150942, 267.071316929, 264.73056806942645, 262.63151827166166, 260.68049791192436, 258.82475321810017, 257.03210503479994, 255.28163873701382, 253.55902813718902, 251.85401008242744, 250.15893877007187, 248.46791776575273, 246.77625790140976, 245.08012741902138, 243.37631999317787, 241.66209753035537, 239.93508186888948, 238.19317935994778, 236.43452813823544, 234.65746143852814, 232.86048252788586, 231.04224823464475, 229.2015589705974, 227.3373537438898, 225.44870905745063, 223.53484084984868, 221.59510880597827, 219.6290224723953, 217.63624867566026, 215.6166197751305, 213.57014229426136, 211.49700547442254, 209.39758928911925, 207.27247145009318, 205.12243293538566, 202.94846157784963, 200.75175327503266, 198.53371042123874, 196.29593722224078, 194.0402316334945, 191.76857376306674, 189.48311069825527, 187.18613784556666, 184.88007701116214, 182.56745158552263, 180.2508593236053, 177.9329433217728, 175.6163618776139, 173.30375797231432, 170.99772913348667, 168.70079841797607, 166.4153872004925, 164.1437903689714, 161.8881544174738, 159.65045879986798, 157.43250077087623, 155.23588380364563, 153.06200954235413, 150.91207313066735, 148.78706165658218, 146.68775537390417, 144.61473130103366, 142.56836875794482, 140.54885637988767, 138.55620013801482, 136.5902318986037, 134.65061805901368, 132.7368678047071, 130.84834053177354, 128.98425196686588, 127.14367848346768, 125.32555905002424, 123.52869413825385, 121.7517407497598, 119.9932024574515, 118.25141296170433, 116.524511061182, 114.81040402467612, 113.10671494188338, 111.41070742191042, 109.71917746874111, 108.02829654801931, 106.3333800275272, 104.62853799288698, 102.9061342632628, 101.15592034384521, 99.36359325095886, 97.50827684619435, 95.55785959714575, 93.4597155301259, 91.1204406254232, 88.35579139090297, 84.74364610264738, 79.06200510052918, 65.89683120062128, 359.37816974724853], [88.98220604699732, 87.56899674436153, 85.46526456923402, 83.2947099615027, 81.10961362675444, 78.92250665187824, 76.73861904031878, 74.56104670110227, 72.39208841588513, 70.23370035028601, 68.08768126773256, 65.95575884081724, 63.83963446389857, 61.74100878490614, 59.661597508909416, 57.60314193047059, 55.5674164112081, 53.55623395864173, 51.571450528002146, 49.614968383842665, 47.68873869832365, 45.794763468515306, 43.935096776327526, 42.111845376329626, 40.32716857058806, 38.583277311302695, 36.88243245912004, 35.22694211638291, 33.61915794978296, 32.061470415896366, 30.556302806143865, 29.106104035228253, 27.713340109530172, 26.380484229736226, 25.110005505493355, 23.904356289305323, 22.765958172133708, 21.69718672380756, 20.700355106542997, 19.777696738315395, 18.9313472327088, 18.163325890921538, 17.475517067160947, 16.869651767780056, 16.347289874177974, 15.909803396843564, 15.558361170529263, 15.293915386650962, 15.117190327801902, 15.028673621040047, 15.028610262841678, 15.116999591969591, 15.29359530064093, 15.557908483720139, 15.909213635012877, 16.346557413913196, 16.868769929068364, 17.474478222037426, 18.1621215857773, 18.929968321677155, 19.776133525064093, 20.698597491799493, 21.695224356021992, 22.763780598828678, 23.90195310685755, 25.107366505326453, 26.377599539159363, 27.7102003257237, 29.10270035114167, 30.552627127325405, 32.057515467510726, 33.61491737328875, 35.22241055553385, 36.87760563512005, 38.5781520870891, 40.32174300433768, 42.106118764383844, 43.92906968581793, 45.788437760044864, 47.682117539122416, 49.60805625188249, 51.564253207623715, 53.548758528336315, 55.559671224318535, 57.595136589766305, 59.65334283629632, 61.732516788053665, 63.83091830249019, 65.9468327968393, 68.07856072752574, 70.22440181285526, 72.38262955673567, 74.55144655569511, 76.72889846158624, 78.91268946093419, 81.09972978637636, 83.28480601970297, 85.45545033550405, 87.55981984500916, 88.98220604699732]])], edgecolor='black', linestyle='dashed', linewidth=2) # 52 degree elevation
#im.show_lines([np.array([[359.37816974724853, 292.38959416697, 279.15361244103434, 273.45644112199915, 269.8386366150942, 267.071316929, 264.73056806942645, 262.63151827166166, 260.68049791192436, 258.82475321810017, 257.03210503479994, 255.28163873701382, 253.55902813718902, 251.85401008242744, 250.15893877007187, 248.46791776575273, 246.77625790140976, 245.08012741902138, 243.37631999317787, 241.66209753035537, 239.93508186888948, 238.19317935994778, 236.43452813823544, 234.65746143852814, 232.86048252788586, 231.04224823464475, 229.2015589705974, 227.3373537438898, 225.44870905745063, 223.53484084984868, 221.59510880597827, 219.6290224723953, 217.63624867566026, 215.6166197751305, 213.57014229426136, 211.49700547442254, 209.39758928911925, 207.27247145009318, 205.12243293538566, 202.94846157784963, 200.75175327503266, 198.53371042123874, 196.29593722224078, 194.0402316334945, 191.76857376306674, 189.48311069825527, 187.18613784556666, 184.88007701116214, 182.56745158552263, 180.2508593236053, 177.9329433217728, 175.6163618776139, 173.30375797231432, 170.99772913348667, 168.70079841797607, 166.4153872004925, 164.1437903689714, 161.8881544174738, 159.65045879986798, 157.43250077087623, 155.23588380364563, 153.06200954235413, 150.91207313066735, 148.78706165658218, 146.68775537390417, 144.61473130103366, 142.56836875794482, 140.54885637988767, 138.55620013801482, 136.5902318986037, 134.65061805901368, 132.7368678047071, 130.84834053177354, 128.98425196686588, 127.14367848346768, 125.32555905002424, 123.52869413825385, 121.7517407497598, 119.9932024574515, 118.25141296170433, 116.524511061182, 114.81040402467612, 113.10671494188338, 111.41070742191042, 109.71917746874111, 108.02829654801931, 106.3333800275272, 104.62853799288698, 102.9061342632628, 101.15592034384521, 99.36359325095886, 97.50827684619435, 95.55785959714575, 93.4597155301259, 91.1204406254232, 88.35579139090297, 84.74364610264738, 79.06200510052918, 65.89683120062128, 359.37816974724853], [88.98220604699732, 87.56899674436153, 85.46526456923402, 83.2947099615027, 81.10961362675444, 78.92250665187824, 76.73861904031878, 74.56104670110227, 72.39208841588513, 70.23370035028601, 68.08768126773256, 65.95575884081724, 63.83963446389857, 61.74100878490614, 59.661597508909416, 57.60314193047059, 55.5674164112081, 53.55623395864173, 51.571450528002146, 49.614968383842665, 47.68873869832365, 45.794763468515306, 43.935096776327526, 42.111845376329626, 40.32716857058806, 38.583277311302695, 36.88243245912004, 35.22694211638291, 33.61915794978296, 32.061470415896366, 30.556302806143865, 29.106104035228253, 27.713340109530172, 26.380484229736226, 25.110005505493355, 23.904356289305323, 22.765958172133708, 21.69718672380756, 20.700355106542997, 19.777696738315395, 18.9313472327088, 18.163325890921538, 17.475517067160947, 16.869651767780056, 16.347289874177974, 15.909803396843564, 15.558361170529263, 15.293915386650962, 15.117190327801902, 15.028673621040047, 15.028610262841678, 15.116999591969591, 15.29359530064093, 15.557908483720139, 15.909213635012877, 16.346557413913196, 16.868769929068364, 17.474478222037426, 18.1621215857773, 18.929968321677155, 19.776133525064093, 20.698597491799493, 21.695224356021992, 22.763780598828678, 23.90195310685755, 25.107366505326453, 26.377599539159363, 27.7102003257237, 29.10270035114167, 30.552627127325405, 32.057515467510726, 33.61491737328875, 35.22241055553385, 36.87760563512005, 38.5781520870891, 40.32174300433768, 42.106118764383844, 43.92906968581793, 45.788437760044864, 47.682117539122416, 49.60805625188249, 51.564253207623715, 53.548758528336315, 55.559671224318535, 57.595136589766305, 59.65334283629632, 61.732516788053665, 63.83091830249019, 65.9468327968393, 68.07856072752574, 70.22440181285526, 72.38262955673567, 74.55144655569511, 76.72889846158624, 78.91268946093419, 81.09972978637636, 83.28480601970297, 85.45545033550405, 87.55981984500916, 88.98220604699732]])], edgecolor='red', linestyle='dashed', linewidth=2) # 55 degree elevation
#im.add_scalebar(10.)
plt.show()
if __name__=='__main__':
plot_image()