基于PYTHON的TEQC图形化展示概要1.docVIP

基于PYTHON的TEQC图形化展示（COMPACT3） 1、拷贝以下代码另存为teqcplot.py 2、需要在Python环境运行 说明：teqcplot version 1.4.2, 7 Sep. 2014. 程序作者：Stuart K. Wier import os import sys import string import datetime import numpy as np import matplotlib as mptl import matplotlib.cm as cm import matplotlib.pyplot as plt from numpy import * from datetime import timedelta from datetime import datetime from matplotlib.pyplot import grid, figure, plot, savefig azifile= elefile= parmfile= parmtype= debug=False dogps=True doglonass=True dogalileo=True dosbas=True dobeidou=True doqzss=True trackCountLimit=6 colorMax=1.234567 colorMin=1.234567 maxHour= -999.0 minHour= -999.0 global lineSize colorname= showLegend=False showLabel= True doParmPlot=False svs=[] svslist=[] parmsvslist=[] options=[] files=[] SVpositionList=[] SVparmList=[] SVidList=[] SVparmidList=[] doShowSVslist=[] doNotShowSVslist=[] global asciiStartTime global maxT global minT widthDistance = 7.3 # for width of plot; inches is units of distance in matplotlib; pixeldensity = 100 def read_input_files (): global asciiStartTime # open the azi file datapath1 = os.path.dirname (azifile) filename1 = os.path.basename(azifile) fileext = string.split(filename1, .) file1 = open (azifile) # open the elevation file datapath2 = os.path.dirname (elefile) filename2 = os.path.basename(elefile) fileext = string.split(filename2, .) file2 = open (elefile) # handle the parm file, if any if len(parmfile) 5: datapathp = os.path.dirname (parmfile) pfilename = os.path.basename(parmfile) if plottype != Time-parameter plot: print Color by value the +parmtype+ data from +pfilename pfile = open (parmfile) # count lines in input files allLines = file1.readlines() file1linecount = len(allLines) file1.seek(0) # check ele file allLines = file2.readlines() if len(allLines) != file1linecount: print Problem: the input