I-Nodal节点仪切割数据时，只能先切割线，再由线合并成炮，而施工sps均为完整sps，即一炮由多条测线构成，切割节点仪数据时，需要对sps数据进行按线拆分，为便于施工，写了本次小脚本。

import pandas as pd
import sys
import time
time_start = time.time()
########################输入输出sps文件都是sps2.1格式，sps软件8.0版本#############################
###################读取原始sps信息，即全sps文件，输出sps2.1格式sps软件8.0版本######################
filename  = 'PHSPS_1122(1-26817)节点'  ####SPS文件名
rline_need = 12222  ####选定测线
filepalce=r'E:/prc/python/节点/'
print('读取 '+filename+' sps文件......')
data = []
for line in open(filepalce+filename + '.X', 'r'):data.append(line)
############计算头卡行数########
i = 0
while data[i][0] == 'H':i = i + 1
ii = i  ##计算头卡行数
###########读取关系文件获得文件号与桩号对应关系####
v = data[ii:]
xlenth = len(v)
del data
datalor=[]
for i in range(0, xlenth):if int(v[i][49:54])==rline_need:  ##筛选需要的检波线datalor.append([int(v[i][10:15]),int(float(v[i][17:24])),int(v[i][28:34]),\int(v[i][38:43]),int(v[i][43:48]),\int(v[i][49:54]),int(v[i][62:66]),int(v[i][72:76])])
print('分选 '+str(rline_need)+' sps文件......')
pddataor2=pd.DataFrame(datalor)
######################删除变量，减少内存开销#############################
del v,datalor
##################检查测线的正确性#######################################
if len(pddataor2)==0:print('无 '+str(rline_need)+' 测线号，请输入正确测线号。')sys.exit(1)
############################################################################
pddataor2.columns=['ffid','shotline','shotpoint','chanfrom','chanto','rline','rpointfrom','rpointto']
pddataor2['tt']=1
pddatal4=pddataor2.values.tolist()
del pddataor2
print('写入 '+str(rline_need)+' sps中...')
f=open(filepalce+str(rline_need)+'线汇总版.X','w+')       # 要更改文件名
f.write('H00 SPS format version num.     SPS2.1;                                         '+"\n")
f.write('H01 Description of survey area  cd2,,N/A,N/A;                                   '+"\n")
f.write('H02 Date of survey              ;                                               '+"\n")
f.write('H021Post-plot date of issue     ;                                               '+"\n")
f.write('H022Tape/disk identifier        DISK;                                           '+"\n")
f.write('H03 Client                      N/A;                                            '+"\n")
f.write('H04 Geophysical contractor      N/A;                                            '+"\n")
f.write('H05 Positioning contractor      N/A;                                            '+"\n")
f.write('H06 Pos. proc. contractor       N/A;                                            '+"\n")
f.write('H07 Field computer system(s)    428XL;                                          '+"\n")
f.write('H08 Coordinate location         Center of source and of receiver pattern;       '+"\n")
f.write('H09 Offset from coord. location 0.0M,0.0DEG;                                    '+"\n")
f.write('H10 Clock time w.r.t. GMT       ;                                               '+"\n")
f.write('H11 Spare                       ;                                               '+"\n")
f.write('H12 Geodetic datum,-spheroid    ;                                               '+"\n")
f.write('H13 Spare                       ;                                               '+"\n")
f.write('H14 Geodetic datum parameters   ;                                               '+"\n")
f.write('H15 Spare                       ;                                               '+"\n")
f.write('H16 Spare                       ;                                               '+"\n")
f.write('H17 Vertical datum description  ;                                               '+"\n")
f.write('H18 Projection type             UTM;                                            '+"\n")
f.write('H19 Projection zone             ;                                               '+"\n")
f.write('H20 Description of grid units   Metres;                                         '+"\n")
f.write('H201Factor to metre                 1.00000000;                                 '+"\n")
f.write('H220Long. of central meridian   ;                                               '+"\n")
f.write('H231Grid origin                 ;                                               '+"\n")
f.write('H232Grid coord. at origin       ;                                               '+"\n")
f.write('H241Scale factor                ;                                               '+"\n")
f.write('H242Lat., long. scale factor    ;                                               '+"\n")
f.write('H256Lat., long. initial line    ;                                               '+"\n")
f.write('H257Circular bearing of H256    ;                                               '+"\n")
f.write('H258Quadrant bearing of H256    ;                                               '+"\n")
f.write('H259Angle from skew             ;                                               '+"\n")
f.write('H300Project code and description;                                               '+"\n")
f.write('H310Line number format          ;                                               '+"\n")
f.write('H400Type,Model,Polarity         ;                                               '+"\n")
f.write('H401Crew name,Comment           ;                                               '+"\n")
f.write('H402Sample int.,Record Len.     1,0.000000,N/A;                                 '+"\n") 
f.write('H403Number of channels          1,5376;                                         '+"\n")
f.write('H404Tape type,format,density    ;                                               '+"\n")
f.write('H405Filter_alias Hz,dB pnt,slope;                                               '+"\n")
f.write('H406Filter_notch Hz,-3dB points ;                                               '+"\n")
f.write('H407Filter_low Hz,dB pnt,slope  ;                                               '+"\n")
f.write('H408Time delay FTB-SOD app Y/N  ;                                               '+"\n")
f.write('H409Multi component recording   ;                                               '+"\n")
f.write('H410Aux. channel 1 contents     ;                                               '+"\n")
f.write('H411Aux. channel 2 contents     ;                                               '+"\n")
f.write('H412Aux. channel 3 contents     ;                                               '+"\n")
f.write('H413Aux. channel 4 contents     ;                                               '+"\n")
f.write('H414Spare                       ;                                               '+"\n")
f.write('H415Spare                       ;                                               '+"\n")
f.write('H416Spare                       ;                                               '+"\n")
f.write('H417Spare                       ;                                               '+"\n")
f.write('H418Spare                       ;                                               '+"\n")
f.write('H419Spare                       ;                                               '+"\n")
f.write('H600Type,model,polarity         G1,geophone 1;                                  '+"\n")
f.write('H601Damp coeff,natural freq.    ;                                               '+"\n")
f.write('H602Nunits,len(X),width(Y)      ;                                               '+"\n")
f.write('H603Unit spacing X,Y            ;                                               '+"\n")
f.write('H604Spare                       ;                                               '+"\n")
f.write('H605Spare                       ;                                               '+"\n")
f.write('H606Spare                       ;                                               '+"\n")
f.write('H607Spare                       ;                                               '+"\n")
f.write('H608Spare                       ;                                               '+"\n")
f.write('H700Type,model,polarity         E1,explosive 1;                                 '+"\n")
f.write('H701Size,vert. stk fold         ;                                               '+"\n")
f.write('H702Nunits,len(X),width(Y)      ;                                               '+"\n") 
f.write('H703Unit spacing X,Y            ;                                               '+"\n")
f.write('H711Nom. shot depth,charge len. ;                                               '+"\n")
f.write('H712Nom. soil,drill method      ;                                               '+"\n")
f.write('H713Weathering thickness        ;                                               '+"\n")
f.write('H26      1         2         3         4         5         6         7         8'+"\n")
for i in range(len(pddatal4)):ffid = int(pddatal4[i][0])shotline = int(pddatal4[i][1])shotpoint = int(pddatal4[i][2])chanfrom = int(pddatal4[i][3])chanto = int(pddatal4[i][4])rline = int(pddatal4[i][5])rpointfrom = int(pddatal4[i][6])rpointto = int(pddatal4[i][7])tt=int(pddatal4[i][8])print('X1 ','     %6d%d%d%-6d       %-6d %d%5d%5d%d%5d        %-4d      %-4d   %d'\%(ffid,tt,tt,shotline,shotpoint,tt,chanfrom,chanto,tt,rline,rpointfrom,rpointto,tt),file=f)
f.close()
print("SPS写入成功！")
time_end=time.time()
timecost = str(round(time_end-time_start,4))
print('耗时 '+timecost+' S')