def Create_interaction(test=False, profile=None, **kwargs): 'Main program to start the moving lines and set up all the map and interactions' warnings.simplefilter(action="ignore", category=FutureWarning) goto_cwd() ui = Create_gui() ui.tb.set_message('Creating basemap') profile = Get_basemap_profile() m = mi.build_basemap(ax=ui.ax1, profile=profile) if profile: sla, slo = profile['Start_lat'], profile['Start_lon'] else: sla, slo = None, None line = init_plot(m, start_lat=sla, start_lon=slo, color='red') flabels = 'labels.txt' faero = 'aeronet_locations.txt' try: ui.tb.set_message('putting labels and aeronet') mi.plot_map_labels(m, flabels) mi.plot_map_labels(m, faero, marker='*', skip_lines=2, color='y') except: print 'Label files not found!' get_datestr(ui) ui.tb.set_message('making the Excel connection') wb = ex.dict_position(datestr=ui.datestr, color=line.get_color(), profile=profile, version=__version__, platform_file=platform_filename, **kwargs) ui.tb.set_message('Building the interactivity on the map') lines = mi.LineBuilder(line, m=m, ex=wb, tb=ui.tb, blit=True) ui.tb.set_message('Saving temporary excel file') savetmp(ui, wb) build_buttons(ui, lines) lines.get_bg(redraw=True) bind_move_window(ui, lines) ui.tb.set_message('Ready for interaction') def stopandquit(): 'simple function to handle the stop and quit' lines.ex.wb.close() ui.root.quit() ui.root.destroy() ui.root.protocol('WM_DELETE_WINDOW', stopandquit) if not test: ui.root.mainloop() return lines, ui
def Create_interaction(test=False,profile=None,**kwargs): 'Main program to start the moving lines and set up all the map and interactions' warnings.simplefilter(action = "ignore", category = FutureWarning) goto_cwd() ui = Create_gui() ui.tb.set_message('Creating basemap') profile = Get_basemap_profile() m = mi.build_basemap(ax=ui.ax1,profile=profile) if profile: sla,slo = profile['Start_lat'],profile['Start_lon'] else: sla,slo = None,None line = init_plot(m,start_lat=sla,start_lon=slo,color='red') flabels = 'labels.txt' faero = 'aeronet_locations.txt' try: ui.tb.set_message('putting labels and aeronet') mi.plot_map_labels(m,flabels) mi.plot_map_labels(m,faero,marker='*',skip_lines=2,color='y') except: print 'Label files not found!' get_datestr(ui) ui.tb.set_message('making the Excel connection') wb = ex.dict_position(datestr=ui.datestr,color=line.get_color(),profile=profile, version=__version__,platform_file=platform_filename,**kwargs) ui.tb.set_message('Building the interactivity on the map') lines = mi.LineBuilder(line,m=m,ex=wb,tb=ui.tb,blit=True) ui.tb.set_message('Saving temporary excel file') savetmp(ui,wb) build_buttons(ui,lines) lines.get_bg(redraw=True) bind_move_window(ui,lines) ui.tb.set_message('Ready for interaction') def stopandquit(): 'simple function to handle the stop and quit' lines.ex.wb.close() ui.root.quit() ui.root.destroy() ui.root.protocol('WM_DELETE_WINDOW',stopandquit) if not test: ui.root.mainloop() return lines,ui
def Create_interaction(**kwargs): #fig,ax = plt.subplots() m = mi.build_basemap() plt.gca().set_title('line segments') lat0, lon0 = mi.pll('22 58.783S'), mi.pll('14 38.717E') x0, y0 = m(lon0, lat0) line, = m.plot([x0], [y0], 'ro-') text = ('Press s to stop interaction\\n' 'Press i to restart interaction\\n') #plt.text(1.0,0.1,text) wb = ex.dict_position(**kwargs) lines = mi.LineBuilder(line, m=m, ex=wb) print 'after line builder' plt.show() print 'after show' plt.gcf().canvas._tkcanvas.master.geometry("900x1000") print 'after canvas' g = gui.gui(lines) g.make_gui() print 'after gui' return lines
def Create_interaction(**kwargs): #fig,ax = plt.subplots() m = mi.build_basemap() plt.gca().set_title('line segments') lat0,lon0 = mi.pll('22 58.783S'), mi.pll('14 38.717E') x0,y0 = m(lon0,lat0) line, = m.plot([x0],[y0],'ro-') text = ('Press s to stop interaction\\n' 'Press i to restart interaction\\n') #plt.text(1.0,0.1,text) wb = ex.dict_position(**kwargs) lines = mi.LineBuilder(line,m=m,ex=wb) print 'after line builder' plt.show() print 'after show' plt.gcf().canvas._tkcanvas.master.geometry("900x1000") print 'after canvas' g = gui.gui(lines) g.make_gui() print 'after gui' return lines
def Create_interaction(test=False,**kwargs): ui = Create_gui() m = mi.build_basemap(ax=ui.ax1) line = init_plot(m,color='red') flabels = 'labels.txt' faero = 'aeronet_locations.txt' try: mi.plot_map_labels(m,flabels) mi.plot_map_labels(m,faero,marker='*',skip_lines=2,color='y') except: print 'Label files not found!' get_datestr(ui) wb = ex.dict_position(datestr=ui.datestr,color=line.get_color(),**kwargs) lines = mi.LineBuilder(line,m=m,ex=wb,tb=ui.tb) savetmp(ui,wb) build_buttons(ui,lines) if not test: ui.root.mainloop() return lines,ui
def Create_interaction(test=False, **kwargs): ui = Create_gui() m = mi.build_basemap(ax=ui.ax1) line = init_plot(m, color='red') flabels = 'labels.txt' faero = 'aeronet_locations.txt' try: mi.plot_map_labels(m, flabels) mi.plot_map_labels(m, faero, marker='*', skip_lines=2, color='y') except: print 'Label files not found!' get_datestr(ui) wb = ex.dict_position(datestr=ui.datestr, color=line.get_color(), **kwargs) lines = mi.LineBuilder(line, m=m, ex=wb, tb=ui.tb) savetmp(ui, wb) build_buttons(ui, lines) if not test: ui.root.mainloop() return lines, ui
rts[i]['tau_fl'] # In[96]: for i,daystr in enumerate(dds): print rts[i]['lat'][rts[i]['fl']][:,0].shape,rts[i]['lon'][rts[i]['fl']][:,0].shape,rts[i]['tau_fl'].shape # In[97]: fig = plt.figure() ax = plt.subplot(111) m = build_basemap(lower_left=[-2,-25],upper_right=[15,-8],ax=ax,larger=False) sa = [] for i,daystr in enumerate(dds): x,y = m(rts[i]['lon'][rts[i]['fl']][:,0]+i*0.03,rts[i]['lat'][rts[i]['fl']][:,0]) sca = ax.scatter(x,y,c=rts[i]['tau_fl'], s=10,alpha=0.7,vmin=0.0,vmax=60.0,edgecolor='None') sa.append(sca) #pu.prelim() cb = plt.colorbar(sa[0]) cb.set_label('COD') plt.savefig(fp+'..//zen_ict/v3/{}_COD_map.png'.format(vr),transparent=True,dpi=600) # In[98]:
rts[i]['tau_fl'] # In[234]: for i,daystr in enumerate(dds): print rts[i]['lat'][rts[i]['fl']][:,0].shape,rts[i]['lon'][rts[i]['fl']][:,0].shape,rts[i]['tau_fl'].shape # In[354]: fig = plt.figure(figsize=(10,6)) ax = plt.subplot(111) m = build_basemap(lower_left=[123,33],upper_right=[131,39],ax=ax,larger=False) sa = [] for i,daystr in enumerate(dds): x,y = m(rts[i]['lon'][rts[i]['fl']][:,0]+i*0.03,rts[i]['lat'][rts[i]['fl']][:,0]) sca = ax.scatter(x,y,c=rts[i]['tau_fl'], s=30,marker='s',alpha=0.5,vmin=0.0,vmax=50.0,edgecolor='None') sa.append(sca) #pu.prelim() plt.title('Cirrus Optical Depth from 4STAR') cb = plt.colorbar(sa[0]) cb.set_label('COD') plt.savefig(fp+'plot/KORUS_COD_map.png',transparent=True,dpi=600) # In[355]:
# In[158]: taufl.x_edge.shape,taufl.y_edge.shape # In[167]: help(m.pcolor) # In[166]: fig = plt.figure() ax = plt.subplot(111) m = build_basemap(lower_left=[-55,40],upper_right=[-35,55],ax=ax,larger=False) x,y = m((taufl.x_edge[:-1]+taufl.x_edge[1:])/2.0,(taufl.y_edge[:-1]+taufl.y_edge[1:])/2.0) cc = m.contourf(x,y,taufl.statistic,vmin=0,vmax=26,cmap=plt.cm.gnuplot2_r) ax.set_title('Mean COD from NAAMES 2015 for boundary layer liquid clouds') cb = plt.colorbar(cc) cb.set_label('COD') #plt.savefig(fp+'../plot/NAAMES_map_COD_stasts.png',dpi=600,transparent=True) # In[132]: reffl.statistic # In[122]:
iams_c3,iams_c3h = lu.load_ict(fp+'dc8\\SEAC4RS-AMS_DC8_20130821_R1.ict',return_header=True) # # Now start plotting the files # In[17]: import map_interactive as mi # ## Check out the cld case #1 # In[27]: fig,ax = plt.subplots(1,1) m = mi.build_basemap(lower_left=[-100,15],upper_right=[-85,32],ax=ax) xt,yt = m(-95.3831,29.7628) ax.text(xt,yt,'+') ax.text(xt,yt,'Houston, TX',horizontalalignment='right',verticalalignment='top') cs = m.scatter(istar_cld['Longitude'],istar_cld['Latitude'],marker='o',latlon=True, c=istar_cld['Start_UTC'],s=20,cmap=plt.cm.rainbow,edgecolor='none') cb = plt.colorbar(cs) cb.set_label('UTC [h]') # In[30]: it = (istar_cld['Start_UTC']>19)&(istar_cld['Start_UTC']<22) # In[31]: