def map_size(nx, xinmap):
# Setting the size of image
# here refer from the uv_limits function in uvinvert.c in difmap
# Radio radian to millisecond
rtomas = 2.062648062470963551564733573307786131966597008796332528822e+8
nmask = 0
my_units = all_class.units()
# Return the max U and V axis in current map pixel
xinc = xinmap / rtomas
uinc = 1.0 / (xinc * nx)
u_limit = uinc * (nx / 4 - nmask)
# Calculate the parameters of U
my_units.nx = nx
my_units.xinc = xinc
my_units.uinc = uinc
my_units.u_limit = u_limit
# Calculate the parameters of V
my_units.ny = nx
my_units.yinc = xinc
my_units.vinc = uinc
my_units.v_limit = u_limit
# Setting the inmap
my_units.xinmap = xinmap
yinmap = xinmap
my_units.yinmap = yinmap
my_units.binwid = 2
return my_units
def map_size(nx,xinmap):
# Setting the size of image
# here refer from the uv_limits function in uvinvert.c in difmap
# Radio radian to millisecond
rtomas =2.062648062470963551564733573307786131966597008796332528822e+8
nmask = 0
my_units = all_class.units()
# Return the max U and V axis in current map pixel
xinc = xinmap/rtomas
uinc = 1.0/(xinc*nx)
u_limit = uinc*(nx/4-nmask)
# Calculate the parameters of U
my_units.nx = nx
my_units.xinc = xinc
my_units.uinc = uinc
my_units.u_limit = u_limit
# Calculate the parameters of V
my_units.ny = nx
my_units.yinc = xinc
my_units.vinc = uinc
my_units.v_limit = u_limit
# Setting the inmap
my_units.xinmap = xinmap
yinmap = xinmap
my_units.yinmap = yinmap
my_units.binwid=2
return my_units
#for j in range(len(uv_data_select[i].split(','))-1):
#uv_data_residual[i].append(uv_data_select[i].split(',')[j])
#uv_data_residual[i].insert(j,uv_data_select[i].split(',')[j])
uv_data_residual[i][2] = uv_data_amp[i]
uv_data_residual[i][3] = uv_data_phs[i]
if all_class.debug:
print '*'*40
print uv_data_residual[0]
# Setting the maplot parameters
nx = 1024
ny = nx
xinmap = 0.1
yinmap = xinmap
domap = 1
my_units = all_class.units()
my_units = map_size.map_size(nx,xinmap)
if all_class.debug:
print 'Show all the members of my_units'
print my_units.nx
print my_units.xinc
print my_units.uinc
print my_units.u_limit
print my_units.ny
print my_units.yinc
print my_units.vinc
print my_units.v_limit
print my_units.xinmap
print my_units.yinmap
print my_units.binwid
#for j in range(len(uv_data_select[i].split(','))-1):
#uv_data_residual[i].append(uv_data_select[i].split(',')[j])
#uv_data_residual[i].insert(j,uv_data_select[i].split(',')[j])
uv_data_residual[i][2] = uv_data_amp[i]
uv_data_residual[i][3] = uv_data_phs[i]
if all_class.debug:
print '*' * 40
print uv_data_residual[0]
# Setting the maplot parameters
nx = 1024
ny = nx
xinmap = 0.1
yinmap = xinmap
domap = 1
my_units = all_class.units()
my_units = map_size.map_size(nx, xinmap)
if all_class.debug:
print 'Show all the members of my_units'
print my_units.nx
print my_units.xinc
print my_units.uinc
print my_units.u_limit
print my_units.ny
print my_units.yinc
print my_units.vinc
print my_units.v_limit
print my_units.xinmap
print my_units.yinmap
print my_units.binwid
