import?arcgisscripting?os?sys?math?string
gp?=?arcgisscripting.create（）
gp.AddToolbox（“C:/Program?Files/ArcGIS/ArcToolbox/Toolboxes/Data?Management?Tools.tbx“）


p?=?range（200）
o?=?range（200）
PI?=?3.1415926535897931
x1=range（200）
x2=range（200）
yy1=range（100）
y2=range（200）
b1=range（200）
l1=range（200）
b2=range（200）
l2=range（200）
k=range（10）
rg=range（200）
aa3?=?0
bb2?=?0.0
aa1?=?1
#輸入坐標文本文件
inputFile?=?sys.argv[1]
#輸出坐標文本文件
outputFile?=?sys.argv[2]
#inputt?the?number?of?the?points?
n?=?int（sys.argv[4]）
#inputt?the?coordiate?54?or?80
d?=?int（sys.argv[3]）
#ensure?the?zhengsuanfansuan?or?huandan
k0?=?int（sys.argv[5]）
#when?fansuan?ensure?the?l0（起始中央子午線經度）
try:
????if?len（sys.argv）?>?9:
????????l0?=?int（sys.argv[8]）
????else:
????????l0=0
except?Exception?ErrorDesc:
????gp.AddError（str（ErrorDesc））
????

#?輸入結束中央子午線經度
ll1?=?int（sys.argv[6]）
#?選擇3度帶或者6度帶
d1?=?int（sys.argv[7]）
????
def?zs1（）:
????print
def?fs1（）:
????print
def?hd（）:
????print


def?k54（）:
????k[0]=6399698.90178271
????k[1]=6367558.49687
????k[2]=32005.7801
????k[3]=133.9213
????k[4]=0.7032
????k[5]=0.005051773902
????k[6]=0.000029838676
????k[7]=6.738525414683492e-3
????k[8]=0.000000241496

def?k80（）:
????k[0]=6399596.65198801
????k[1]=6367452.13279
????k[2]=32009.8575
????k[3]=133.9602
????k[4]=0.6976
????k[5]=0.005052505593
????k[6]=0.000029847335
????k[7]=6.739501819472925e-3
????k[8]=0.0000002416

def?jdzh1（jd）:
????i=float（jd）
????i?=?int（i）
????
????j=float（jd）*100-int（i）*100
????j?=?int（j）
????##gp.AddError（i）
????##gp.AddError（j）
????k1=float（jd）*10000-float（i）*10000-float（j）*100
????##gp.AddError（k1）
????
????if?k1>=60:
????????j=j+1
????????k1=0
????jdd=float（i）*PI/180+float（j）*PI/180/60+float（k1）*PI/180/3600
????##gp.AddError（jdd）
????return?jdd

def?jdzh2（jd）:
????
????
????i=float（jd）*180/PI
????i?=?int（i）
????
????j=（float（jd）*180/PI-int（i））*60
????j=int（j）
????
????k=（（float（jd）*180/PI-int（i））*60-int（j））*60
????k?=?float（k）
????
????jdd=float（i）+float（j）/100+float（k）/10000
???
????return?jdd
i?=?0
if?k0==1:
???f=open（inputFile）
???for?a?in?f.readlines（）:
???????fu?=?a.split（‘‘）
???????b1[i]=fu[0]
???????
???????l1[i]=fu[1]
???????
???????i?+=1
else:
?????f=open（inputFile）
?????for?a?in?f.readlines（）:
?????????fu?=?a.split（‘‘）
?????????yy1[i]=float（fu[0]）
???
?????????x1[i]=float（fu[1]）
?????????i?+=?1
def?zs（）:
????o=range（20）
????##gp.AddError（b2[i]）
????
????o[1]=math.tan（b2[i]）
????##gp.AddError（o[1]）
????
????o[2]=math.cos（b2[i]）
????##gp.AddError（o[2]）
????
????o[3]=k[7]*o[2]*o[2]
????##gp.AddError（o[3]）
????
????o[4]=o[1]*o[1]
????##gp.AddError（o[4]）
????
????o[5]=1+o[3]
????##gp.AddError（o[5]）
????
????o[6]=k[0]/math.sqrt（o[5]）
????##gp.AddError（o[6]）
????
????o[7]=（l2[i]-ll2）*o[2]
????o[8]=o[7]*o[7]
????o[9]=o[1]*o[2]
????o[10]=o[9]*o[9]
????o[11]=（k[2]+o[10]*（k[3]+o[10]*k[4]））
????###gp.AddError（o[11]）
???