function?out=teqcplot（files）;

%?TEQCPLOT??Read?and?plot?TEQC?report?files
%
%???TEQC?is?the?Toolkit?for?GPS/GLONASS/SBAS?Data?used?to?solve?
%???many?pre-processing?problems?with?GPS?GLONASS?and?SBAS?data:
%???TEQC?stands?for?Translation?Editing?and?Quality?Check
%???More?info?here:?http://facility.unavco.org/software/teqc/teqc.html
%
%???TEQCPLOT（FILENAME）?returns?the?content?of?a?TEQC?reportfile?in
%???a?figure?（and?a?struct?variable）.
%???Input?argument?FILENAME?is?optional.?
%
%???Valid?TEQC?report?files?are:
%??
%???FORMAT??DEscriptION
%???*.sn1???Signal?to?noise?ratio?（S/N）?Carrier?L1
%???*.sn2???Signal?to?noise?ratio?（S/N）?Carrier?L2
%???*.iod???Derivative?of?ionospheric?delay?observable?（m/s）
%???*.ion???Ionospheric?delay?observable?（m）
%???*.mp1???Multipath?Carrier?L1
%???*.mp2???Multipath?Carrier?L2
%???*.azi???Satellite?azimuthal?data?（degrees）
%???*.ele???Satellite?elevation?data?（degrees）
%
%???Jim?Hedfors?Uppsala?University?2005
%???Comments?and?requests?for?more?TEQC?Matlab?tools?to:
%???jim.hedfors@gmail.com

%?FUNCTION:?READER?&?PLOTTER?++++++++++++++++++++++++++++++++++++++++++++++

if?nargin==0
????[filenpath]=uigetfile（‘*.sn1;*.sn2;*.iod;*.ion;*.mp1;*.mp2;*.azi;*.ele‘...
????????‘Pick?your?TEQC?report?file‘）;
else
????[pathfilenext]=fileparts（files）;
????path=[path?‘\‘];
????filen={[filen?ext]};
end
file=char（filen）;
%[A]=importdata（[path?file]‘\t‘）;
A=textread（[path?file]?‘%s‘‘whitespace‘‘\n‘）
SAT（1:length（A）1:32）=NaN;
n=0;i=5;
sats=str2num（A{i}）;
snyggfilen=strrep（filen‘_‘‘-‘）;
h=waitbar（0[‘Please?wait...?Importing?‘?char（snyggfilen）]）;
while?i????waitbar（i/length（A）h）;drawnow
????n=n+1;
????i=i+1;
????if?sats~=0;
????????SAT（nsats（2:end））=str2num（A{i}）;
????????i=i+1;
????end
????if?str2num（A{i}）~=-1;
????????sats=str2num（A{i}）;
????end
end
close（h）;
sat=SAT（1:n+1:）;
t_samp=char（A（3））;
mjl=char（A（4））;
T_SAMP=str2num（t_samp（max（find（t_samp==‘?‘））:end））;
MJL_START=str2num（mjl（max（find（mjl==‘?‘））:end））;
MJL（1）=MJL_START;
for?i=2:length（sat）;
????MJL（i）=MJL_START+i*T_SAMP/60/60/24;
end
JD=mjl2jd（MJL）;
[typemaxyminy]=get_filetype（file）;

figure;box?on;hold?on
pcolor（JD‘1:32sat‘）
set（gcf‘renderer‘‘zbuffer‘）;
shading?flat
set（gca‘xticklabel‘[JD]）;
set（gca‘xlim‘[JD（1）?JD（end）]）
dateaxis（‘X‘13）
cbar（‘v‘[miny?maxy]type）;
set（gca‘ylim‘[1?33]）
set（gca‘ytick‘[1.5:1:32.5]）
set（gca‘yticklabel‘[‘S01‘;‘S02‘;‘S03‘;‘S04‘;‘S05‘;‘S06‘;...
????????‘S07‘;‘S08‘;‘S09‘;‘S10‘;‘S11‘;‘S12‘;...
????????‘S13‘;‘S14‘;‘S15‘;‘S16‘;‘S17‘;‘S18‘;...
????????‘S19‘;‘S20‘;‘S21‘;‘S22‘;‘S23‘;‘S24‘;...
????????‘S25‘;‘S26‘;‘S27‘;‘S28‘;‘S29‘;‘S30‘;...
????????‘S31‘;‘S32‘]）;
set（gca‘fontsize‘7）;
colormap（flipud（jet））;
caxis（[miny?maxy]）;
xlabel（[datestr（JD（1））?‘???|--------??T?Samp:?‘?num2str（T_SAMP）?‘?s??--------|????‘?datestr（JD（end））]）
ylabel（‘SVS‘）
timestr=secs2hms（length（sat））;
T=title（[‘TEQC?Report?file:?‘?strrep（fi

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件????7527768??2018-01-22?14:22??teqcplot\result.ion

?????文件????7527768??2018-01-22?14:19??teqcplot\result.mp1

?????文件?????184937??2006-05-12?15:10??teqcplot\rover.sn2

?????文件???????6109??2017-09-04?15:17??teqcplot\teqcplot.m

?????目錄??????????0??2018-01-22?14:25??teqcplot

-----------?---------??----------?-----??----

?????????????15246582????????????????????5