python2.7——GSI-16格式水准数据平差

fd = open('E:\pythontext/17110701.GSI','r')file=fd.read()result=list(file.split())del result[0]start=[result[0],result[1]]del result[0]del result[0]l=len(result)/18start_station=result[0][-4:]start_high=int(start[1][7:])#final is the data about station and high(without adjusting)final_name=[]final_high=[]distance=[]final_data=range(0,l)sum_distance=range(0,l)final_name.append(start_station)final_name.append(result[3][-4:])final_high.append(str(start_high))for i in range(l):    if result[3+i*18][-4:] not in final_name:            final_name.append(result[3+i*18][-4:])    if result[0+i*18][-4:] not in final_name:        final_name.append(result[0 + i * 18][-4:])    final_high.append(result[17+i*18][6:])    distance.append(result[16 + i * 18][6:])for i in range(l):    final_high[i]=int(final_high[i])*1.0/100000    sum_distance[i]=int(distance[i])*1.0/100000#adjustmentbihecha=(int(final_high[l])-start_high)*1.0/100000sumdistance=sum_distance[l-1]print final_highfinal_high.pop(l)for i in range(l):    if i==0:        final_data[i]=round(final_high[i],4)    else:        final_data[i]=round(final_high[i]-bihecha*sum_distance[i]/sumdistance,4)# outputoutput=open('E:\pythontext/17110701.txt','w')for i in range(l):    output.write(final_name[i])    output.write(' ')    output.write(str(final_data[i]))    output.write('\n')output.close()fd.close()