test_module_toolbox_num.py File Reference

Namespaces
	test_module_toolbox_num

Variables
tuple	test_module_toolbox_num.folder_test = os.path.join('dir.test','TEST')
tuple	test_module_toolbox_num.test_mode = module_prompt.prompt_inlist(['All','fit_gaussian'])
string	test_module_toolbox_num.title = 'Test: fit_gaussian_longitude and fit_gaussian_latitude routines.'
tuple	test_module_toolbox_num.fig_size = (12,6)
int	test_module_toolbox_num.fig_resolution = 100
string	test_module_toolbox_num.axis_name_y = 'Latitude'
string	test_module_toolbox_num.axis_name_x = 'Carrington Longitude'
int	test_module_toolbox_num.axis_size = 10
string	test_module_toolbox_num.axis_color = 'k'
int	test_module_toolbox_num.dphi = 20
list	test_module_toolbox_num.xtick_loc = [ i*dphi for i in range(0,int(360./dphi)+1) ]
list	test_module_toolbox_num.xtick_label = [ str(int(xtick_loc[i])) for i in range(0,len(xtick_loc)) ]
list	test_module_toolbox_num.ytick_loc = [ -90,-80,-60,-40,-20,0,20,40,60,80,90 ]
list	test_module_toolbox_num.ytick_label = [ str(int(ytick_loc[i])) for i in range(0,len(ytick_loc)) ]
int	test_module_toolbox_num.n_in = 1000
int	test_module_toolbox_num.sigma_in = 5
int	test_module_toolbox_num.bin_delta = 1
tuple	test_module_toolbox_num.seed_point = module_timespace.Coordinate()
list	test_module_toolbox_num.points = [ module_timespace.copy_coordinate(seed_point) for ipoint in range(0,n_in) ]
tuple	test_module_toolbox_num.longitudes = (np.random.normal(seed_point.lon,sigma_in,n_in)+360.)
tuple	test_module_toolbox_num.latitudes = np.random.normal(seed_point.lat,sigma_in,n_in)
tuple	test_module_toolbox_num.angles = np.arange(0,2*np.pi,2*np.pi/100)
tuple	test_module_toolbox_num.ellipse_x = (mean_fit_lon + 3*sigma_fit_lon*np.cos(angles) + 360.)
int	test_module_toolbox_num.ellispe_y = mean_fit_lat+3
tuple	test_module_toolbox_num.fig = plt.figure(num=1,figsize=fig_size,dpi=fig_resolution)
tuple	test_module_toolbox_num.sub = plt.subplot(111)
tuple	test_module_toolbox_num.data = sub.scatter([point.lon for point in points],[point.lat for point in points],color='r',s=10)
tuple	test_module_toolbox_num.ellipse = sub.scatter(ellipse_x,ellispe_y,color='b',s=1)
tuple	test_module_toolbox_num.center = sub.scatter(mean_fit_lon,mean_fit_lat,color='b',marker='X',s=10)
tuple	test_module_toolbox_num.legend = sub.legend(handles=[data,ellipse,center],labels=['Data','3*sigma','Mean fitted value'])
string	test_module_toolbox_num.title1 = 'Lon0-Lon_fit = %6.2f'
string	test_module_toolbox_num.title2 = 'Lat0-Lat_fit = %6.2f'
string	test_module_toolbox_num.title3 = 'Sigma0-Sigma_fit = (%6.1e,%6.1e)'
string	test_module_toolbox_num.title4 = 'err(X0) = (%6.1e,%6.1e)'
string	test_module_toolbox_num.title5 = 'err(sigma) = (%6.1e,%6.1e)'
string	test_module_toolbox_num.command = 'eog %s/test_fit_gaussian_*.png'
	plt.show() plt.close() More...