diff --git a/.gitignore b/.gitignore index f9091dc..9dc7d33 100644 --- a/.gitignore +++ b/.gitignore @@ -1,12 +1,20 @@ -# Ignore everything -* -# But not .py files -!*.py -# Include .gitignore itself -!.gitignore -# Recursively allow .py files in subdirectories -!*/ +# # Ignore everything +# * +# # But not .py files +# !*.py +# # Include .gitignore itself +# !.gitignore +# # Recursively allow .py files in subdirectories +# !*/ -dist +# dist +# build +# !.gitkeep +# end old ignore + +__pycache__ build -!.gitkeep \ No newline at end of file +dist +*.spec +notebooks +passcode \ No newline at end of file diff --git a/.vscode/launch.json b/.vscode/launch.json new file mode 100644 index 0000000..c5b047e --- /dev/null +++ b/.vscode/launch.json @@ -0,0 +1,19 @@ +{ + // Use IntelliSense to learn about possible attributes. + // Hover to view descriptions of existing attributes. + // For more information, visit: https://go.microsoft.com/fwlink/?linkid=830387 + "version": "0.2.0", + "configurations": [ + { + "name": "Python Debugger: Backend", + "type": "debugpy", + "request": "launch", + "program": "./backend.py", + "console": "integratedTerminal", + // "justMyCode": false, + "args": [ + "${command:pickArgs}" + ] + } + ] +} \ No newline at end of file diff --git a/.vscode/settings.json b/.vscode/settings.json new file mode 100644 index 0000000..f993589 --- /dev/null +++ b/.vscode/settings.json @@ -0,0 +1,3 @@ +{ + "editor.formatOnSave": true +} diff --git a/SimHei.ttf b/SimHei.ttf new file mode 100644 index 0000000..c5030ae Binary files /dev/null and b/SimHei.ttf differ diff --git a/cspell.config.yaml b/cspell.config.yaml new file mode 100644 index 0000000..b56761f --- /dev/null +++ b/cspell.config.yaml @@ -0,0 +1,7 @@ +version: "0.2" +ignorePaths: ["/**"] +dictionaryDefinitions: [] +dictionaries: [] +words: [] +ignoreWords: [] +import: [] diff --git a/modules/tidi/tidi.m b/modules/tidi/tidi.m new file mode 100644 index 0000000..05ba320 --- /dev/null +++ b/modules/tidi/tidi.m @@ -0,0 +1,266 @@ +%% TIDI +clear +clc +SourcePath='.\\radar\\data\\'; +Files=dir(fullfile(SourcePath,'*.')); +LengthFiles=length(Files); +for i=4:4 + Path=strcat(SourcePath,Files(i).name,'\'); + Sub_Files=dir(fullfile(Path,'*.VEC')); + Sub_LengthFiles=length(Sub_Files); + for j=1:Sub_LengthFiles + Date(j)=str2num(Sub_Files(j).name(13:15)); + FilestypeD(j)=str2num(Sub_Files(j).name(30:32)); + FilestypeR(j)=str2num(Sub_Files(j).name(35:36)); + end + m=0; + for t=1:366 + [m]=find(Date==t); + if isempty(m)==0 + [maxValPosD]=find(FilestypeD(m)==max(FilestypeD(m))); + n=m(maxValPosD); + [maxValPosR]=find(FilestypeR(n)==max(FilestypeR(n))); + filename=strcat(Path,Sub_Files(n(maxValPosR)).name); + Vzonal=ncread(filename,'u'); + VMerdional=ncread(filename,'v'); + UTime=ncread(filename,'ut_time'); + Lat=ncread(filename,'lat'); + Lon=ncread(filename,'lon'); + Height=ncread(filename,'alt_retrieved'); + UTime=double(UTime)/(1000*3600); + [R,C]=size(Vzonal); + for ii=1:R + for jj=1:C + if Vzonal(ii,jj)==-9999||Vzonal(R,C)==0 + Vzonal(ii,jj)=nan; + end + if VMerdional(ii,jj)==-9999||VMerdional(R,C)==0 + VMerdional(ii,jj)=nan; + end + end + end + + save(strcat(Path,filename(38:40),'_Lat.mat'),'Lat'); + save(strcat(Path,filename(38:40),'_Lon.mat'),'Lon'); + save(strcat(Path,filename(38:40),'_Height.mat'),'Height'); + save(strcat(Path,filename(38:40),'_UTime.mat'),'UTime'); + save(strcat(Path,filename(38:40),'_Vzonal.mat'),'Vzonal'); + save(strcat(Path,filename(38:40),'_VMerdional.mat'),'VMerdional'); + % [row col] =size(Vzonal); + % jj =1; + % VzonalWuhan(1:row,1,ii) =nan; + % Var_VzonalWuhan(1:row,1,ii) =nan; + % UTWuhan(1,1,ii) =0; + % latWuhan(1,1,ii) =nan; + % lonWuhan(1,1,ii) =nan; + % for rr =1:col + % if lat(rr,1)>9.52 && lat(rr,1)<29.52 %漠河（122.37°E， 53.50°N，磁倾角71°）、北京（116.25°E， 40.25°N， 磁倾角 59°） + % if lon(rr,1)>99.13 && lon(rr,1)<119.13 %武汉（114.61°E, 30.53°N, 磁倾角46°）、三亚（109.13°E, 19.52°N, 磁倾角27°） + % VzonalWuhan(:,jj,ii) =Vzonal(:,rr); + % Var_VzonalWuhan(:,jj,ii) =Var_Vzonal(:,rr); + % UTdateWuhan(:,jj,ii) =UTdate(:,rr); + % UTWuhan(:,jj,ii) =UT(rr,1); + % latWuhan(:,jj,ii) =lat(rr,1); + % lonWuhan(:,jj,ii) =lon(rr,1); + % jj =jj+1; + % end + % end + % end + % else + % NoData(tt) =ii; + % tt =tt+1; + end + % VzonalWorld(1:row,1:col,ii) =Vzonal; + % UTWorld(1:col,1,ii) =UT; + % latWorld(1:col,1,ii) =lat; + % lonWorld(1:col,1,ii) =lon; + end + clearvars -except SourcePath Files LengthFiles +end +% altWuhan =alt; +% altWorld =alt; +% clear alt +% clear row +% clear col +% clear Vzonal +% clear UTdate +% clear UT +% clear lat +% clear lon +% [row col day] =size(VzonalWuhan); +% Nullpoint =zeros(col,day); +% NullpointUT =zeros(col,day); +% Nullpointheight =NaN(row,col,day); +% for tt =1:day +% for jj =1:col +% for ii =1:row +% if VzonalWuhan(ii,jj,tt) ==-9999 || VzonalWuhan(ii,jj,tt) ==0 +% VzonalWuhan(ii,jj,tt) =nan; +% end +% if Var_VzonalWuhan(ii,jj,tt) ==-9000000 || Var_VzonalWuhan(ii,jj,tt) ==0 +% Var_VzonalWuhan(ii,jj,tt) =nan; +% end +% end +% for ii=1:row-1 +% if (VzonalWuhan(ii,jj,tt)*VzonalWuhan(ii+1,jj,tt))<0 +% Nullpoint (jj,tt) =Nullpoint(jj,tt)+1; +% Y =altWuhan(ii+1,1)-altWuhan(ii,1); +% X =VzonalWuhan(ii+1,jj,tt)-VzonalWuhan(ii,jj,tt); +% Nullpointheight (ii,jj,tt) =altWuhan(ii,1)-(Y/X)*VzonalWuhan(ii,jj,tt); +% NullpointUT (jj,tt) =UTWuhan(1,jj,tt); +% end +% end +% end +% end +% NullpointTotal =zeros(1,24); +% for ii =1:col +% for jj=1:day +% if NullpointUT(ii,jj) ~=0 +% m =floor(NullpointUT(ii,jj)); +% NullpointTotal(1,m+1) =NullpointTotal(1,m+1)+1; +% end +% end +% end +% m =1; +% n =1; +% t =1; +% for tt =1:day +% for jj =1:col +% for ii =1:row +% if isnan(Nullpointheight(ii,jj,tt)) ==0 +% NullHeight(m) =Nullpointheight(ii,jj,tt); +% NullLT(m) =NullpointUT(jj,tt)+8; +% if NullLT(m) >=24 +% NullLT(m) =NullLT(m)-24; +% end +% m=m+1; +% end +% if isnan(VzonalWuhan(ii,jj,tt)) ==0 +% WuhanHeight(n) =altWuhan(ii,1); +% WuhanLT(n) =UTWuhan(1,jj,tt)+8; +% if WuhanLT(n) >=24 +% WuhanLT(n) =WuhanLT(n)-24; +% end +% n =n+1; +% end +% end +% end +% end +% [row col day] =size(VzonalWorld); +% m =1; +% for tt =1:day +% for jj =1:col +% for ii =1:row +% if VzonalWorld(ii,jj,tt) ==-9999 || VzonalWorld(ii,jj,tt) ==0 +% VzonalWorld(ii,jj,tt) =nan; +% end +% if isnan(VzonalWorld(ii,jj,tt)) ==0 +% WorldHeight(m) =altWorld(ii,1); +% WorldLT(m) =UTWorld(jj,1,tt)+lonWorld(jj,1,tt)/15; +% if WorldLT(m) >=24 +% WorldLT(m) =WorldLT(m)-24; +% end +% m =m+1; +% end +% end +% end +% end + + +%% +% ax(1)=gca; +% set(ax(1),'Position',[0.12 0.12 0.75 0.7]) +% nx=0:1:16; +% set(ax(1),'Xtick',nx) +% set(ax(1),'Xlim',[0,32]) +% ax(2)=axes('Position',get(ax(1),'Position'),... +% 'XAxisLocation','top',... +% 'YAxisLocation','right',... +% 'Color','none',... +% 'XColor','k','YColor','k'); +% set(ax,'box','off') +% nx=0:100:3200; +% set(ax(2),'Xtick',nx) +% set(ax(2),'Xlim',[0,3200]) +% set(ax(2),'Ylim',[70,120]) +% set(ax(1),'Ylim',[70,120]) +% X =zeros(21,1); +% E =sqrt(Var_VzonalWuhan(:,1,10)); +% hold on +% fill([100,300,300,100],[70,70,120,120],[0.8,0.8,0.8]); +% hold on +% errorbare('h',VzonalWuhan(:,1,10)+200,altWuhan,E); +% hold on +% X(1:21,1) =200; +% plot(X,altWuhan,'--k'); +% hold on +% plot(Zonal0(:,1)+200,height0(:,1),'r'); +% hold on +% E =sqrt(Var_VzonalWuhan(:,2,10)); +% hold on +% fill([500,700,700,500],[70,70,120,120],[0.8,0.8,0.8]); +% hold on +% errorbare('h',VzonalWuhan(:,2,10)+600,altWuhan,E); +% hold on +% X(1:21,1) =600; +% plot(X,altWuhan,'--k'); +% E =sqrt(Var_VzonalWuhan(:,4,10)); +% hold on +% fill([900,1100,1100,900],[70,70,120,120],[0.8,0.8,0.8]); +% hold on +% errorbare('h',VzonalWuhan(:,4,10)+1000,altWuhan,E); +% hold on +% X(1:21,1) =1000; +% plot(X,altWuhan,'--k'); +% E =sqrt(Var_VzonalWuhan(:,5,10)); +% hold on +% fill([1300,1500,1500,1300],[70,70,120,120],[0.8,0.8,0.8]); +% hold on +% errorbare('h',VzonalWuhan(:,5,10)+1400,altWuhan,E); +% hold on +% X(1:21,1) =1400; +% plot(X,altWuhan,'--k'); +% E =sqrt(Var_VzonalWuhan(:,6,10)); +% hold on +% fill([1700,1900,1900,1700],[70,70,120,120],[0.8,0.8,0.8]); +% hold on +% errorbare('h',VzonalWuhan(:,6,10)+1800,altWuhan,E); +% hold on +% X(1:21,1) =1800; +% plot(X,altWuhan,'--k'); +% E =sqrt(Var_VzonalWuhan(:,7,10)); +% hold on +% fill([2100,2300,2300,2100],[70,70,120,120],[0.8,0.8,0.8]); +% hold on +% errorbare('h',VzonalWuhan(:,7,10)+2200,altWuhan,E); +% hold on +% X(1:21,1) =2200; +% plot(X,altWuhan,'--k'); +% E =sqrt(Var_VzonalWuhan(:,8,10)); +% hold on +% fill([2500,2700,2700,2500],[70,70,120,120],[0.8,0.8,0.8]); +% hold on +% errorbare('h',VzonalWuhan(:,8,10)+2600,altWuhan,E); +% hold on +% X(1:21,1) =2600; +% plot(X,altWuhan,'--k'); +% E =sqrt(Var_VzonalWuhan(:,3,10)); +% hold on +% fill([2900,3100,3100,2900],[70,70,120,120],[0.8,0.8,0.8]); +% hold on +% errorbare('h',VzonalWuhan(:,3,10)+3000,altWuhan,E); +% hold on +% X(1:21,1) =3000; +% plot(X,altWuhan,'--k'); +% set(ax,'YMinorTick','on','XMinorTick','off','Xticklabel',{''},... +% 'TickDir','in',... +% 'TickLength',[0.025 0.025],... +% 'Layer','top'); +% set(ax(1),'xticklabel',{'','','8.82','','','','8.85','','','','10.47','','','','10.50'}); +% set(ax(2),'xticklabel',{'','-100','0','100','','-100','0','100','','-100','0','100','','-100','0','100'}); +% set(get(ax(1),'xlabel'),'string','UT (hour)'); +% set(get(ax(2),'xlabel'),'string','Zonal Wind (m/s)'); +% set(get(ax(1),'ylabel'),'string','Height (km)'); +% set(get(ax(2),'ylabel'),'string','Height (km)'); + +