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