From 33a9e3dc2cec6443944ab0ccc37b21313b14a693 Mon Sep 17 00:00:00 2001 From: L_DelOff <51275636+LDelOff@users.noreply.github.com> Date: Sat, 12 Jun 2021 23:31:37 +0300 Subject: [PATCH] =?UTF-8?q?=D0=90=D0=B4=D0=B0=D0=BF=D1=82=D0=B8=D1=80?= =?UTF-8?q?=D0=BE=D0=B2=D0=B0=D0=BD=D1=8B=20=D0=B4=D0=BB=D1=8F=20=D1=80?= =?UTF-8?q?=D0=B0=D0=B1=D0=BE=D1=82=D1=8B=20=D1=81=20=D0=B0=D0=B2=D1=82?= =?UTF-8?q?=D0=BE=D0=BC=D0=B0=D1=82=D0=B8=D1=87=D0=B5=D1=81=D0=BA=D0=B8?= =?UTF-8?q?=D0=BC=20=D0=B0=D0=BB=D0=B3=D0=BE=D1=80=D0=B8=D1=82=D0=BC=D0=BE?= =?UTF-8?q?=D0=BC?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- .gitignore | 97 ++++++++++++++++++++++++++++++++++++ analysis.m | 132 +++++++++++++++++++++++++++++++------------------ loaddata_edr.m | 83 +++---------------------------- 3 files changed, 188 insertions(+), 124 deletions(-) diff --git a/.gitignore b/.gitignore index b6dd41f..7d385ff 100644 --- a/.gitignore +++ b/.gitignore @@ -916,3 +916,100 @@ /DATA/f18/ssies/2020/06/PS.CKGWC_SC.U_DI.A_GP.SIES3-F18-R99990-B9999090-APGA_AR.GLOBAL_DD.20200630_TP.000001-235959_DF.EDR /DATA/f18/ssies/2020/06/PS.CKGWC_SC.U_DI.A_GP.SIES3-F18-R99990-B9999090-APGA_AR.GLOBAL_DD.20200630_TP.000001-235959_DF.EDR.gz /DATA/f18/ssies/2020/06/PS.CKGWC_SC.U_DI.A_GP.SIES3-F18-R99990-B9999090-APGA_AR.GLOBAL_DD.20200630_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140901_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140901_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140901_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140902_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140902_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140902_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140903_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140903_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140903_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140904_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140904_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140904_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140905_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140905_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140905_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140906_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140906_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140906_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140907_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140907_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140907_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140908_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140908_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140908_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140909_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140909_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140909_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140910_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140910_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140910_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140911_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140911_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140911_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140912_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140912_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140912_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140913_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140913_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140913_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140914_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140914_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140914_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140915_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140915_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140915_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140916_TP.000001-225240_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140916_TP.000001-225240_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140916_TP.000001-225240_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140917_TP.002818-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140917_TP.002818-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140917_TP.002818-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140918_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140918_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140918_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140919_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140919_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140919_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140920_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140920_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140920_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140921_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140921_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140921_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140922_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140922_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140922_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140923_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140923_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140923_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140924_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140924_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140924_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140925_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140925_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140925_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140926_TP.000002-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140926_TP.000002-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140926_TP.000002-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140927_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140927_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140927_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140928_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140928_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140928_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140929_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140929_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140929_TP.000001-235959_DF.mat +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140930_TP.000001-235959_DF.EDR +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140930_TP.000001-235959_DF.EDR.gz +/DATA/f15/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES2-F15-R99990-B9999090-APGA_AR.GLOBAL_DD.20140930_TP.000001-235959_DF.mat +/loaddata_edr.asv +/Program.asv +/DATA/f16/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES3-F16-R99990-B9999090-APGA_AR.GLOBAL_DD.20140901_TP.000001-235959_DF.EDR.gz +/DATA/f16/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES3-F16-R99990-B9999090-APGA_AR.GLOBAL_DD.20140902_TP.000002-235959_DF.EDR.gz +/DATA/f16/ssies/2014/09/PS.CKGWC_SC.U_DI.A_GP.SIES3-F16-R99990-B9999090-APGA_AR.GLOBAL_DD.20140903_TP.000001-235959_DF.EDR.gz.html +/analysis.asv +/mymean.asv diff --git a/analysis.m b/analysis.m index 9411168..27bfe2a 100644 --- a/analysis.m +++ b/analysis.m @@ -3,73 +3,102 @@ function param=analysis(param) % собираю данные по землетрясению %% Поехали type_sensor=param.type_sensor; -sattelite_range=param.sattelite_range; -year_range=param.year_range; -month_range=param.month_range; -day_range=param.day_range; DATA.DATA=[]; %% Объединяю данные в диапазоне if 1 - for sattelite=sattelite_range %6:20 - for year=year_range % года - for month=month_range % месяца - for day=day_range % дни - DATA_temp=load_file(sattelite, year, month, day,... - type_sensor); %загрузка данных - DATA.DATA=[DATA.DATA DATA_temp.DATA]; - end - end - end + for i=1:length(param.dates) + year=param.dates(i).year; + month=param.dates(i).month; + day=param.dates(i).day; + DATA_temp=load_file(param.sattelite_range, year, month, day,... + type_sensor); %загрузка данных + DATA.DATA=[DATA.DATA DATA_temp.DATA]; end end -%load matlab %% Сортировка данных R_quake=param.R_quake; % радиус действия относительно магнитуды,км (еще тыщу ставил) % расположение эпицентра fiA=param.fiA; % широта LA=param.LA; % долгота -report=[]; % тут будут все данные для графика -global zzz ccc -zzz=[]; -ccc=[]; -for i=day_range% дни, по которыс строится среднесуточные значения - report(end+1).date=[myformat(i,2) '.03.2011']; + +fname_rep=param.filename_report; + +if exist(fname_rep,'file') + load(fname_rep,'report'); % тут будут все данные для графика +else + %report=[]; + %report.data=[]; + report.data.N_sat=[]; + report.data.R=[]; +end + + +for i=1:length(param.dates) HHeden=[]; O2den=[]; Eden=[]; Etemp=[]; + for j=1:length(DATA.DATA) - if DATA.DATA(j).year==year_range % сортировка по году - if DATA.DATA(j).month==month_range % по месяцу - if (DATA.DATA(j).day==i)% по числам - for k=1:3 % координаты спутника только 3 раза в минуту(кадр) записываются - fiB=DATA.DATA(j).GeoLat(k); % 1. Географическая широта (градусы, север)(Geographic latitude (degrees, north)) - LB=DATA.DATA(j).GeoLong(k); % 2. Географическая долгота (градусы, восток)(Geographic longitude (degrees, east)) - R=mydistance(fiA,LA,fiB,LB); - ccc(1,end+1)=fiB; - ccc(2,end)=LB; - - zzz(end+1)=R-R_quake; - if R<=R_quake - HHeden(end+1)= mymean(DATA.DATA(j).RPA_HHeden((k-1)*5+1:(k-1)*5+5),1e10); - O2den(end+1) = mymean(DATA.DATA(j).RPA_O2den((k-1)*5+1:(k-1)*5+5),1e10); - Eden(end+1) = mymean(DATA.DATA(j).EP_Eden((k-1)*5+1:(k-1)*5+5),1e10); - Etemp(end+1) = mymean(DATA.DATA(j).EP_Etemp((k-1)*5+1:(k-1)*5+5),1e10); - end - end + date_true=1; + if DATA.DATA(j).year~=param.dates(i).year % сортировка по году + date_true=0; + end + if DATA.DATA(j).month~=param.dates(i).month % по месяцу + date_true=0; + end + if (DATA.DATA(j).day~=param.dates(i).day)% по числам + date_true=0; + end + if date_true + for k=1:3 % координаты спутника только 3 раза в минуту(кадр) записываются + fiB=DATA.DATA(j).GeoLat(k); % 1. Географическая широта (градусы, север)(Geographic latitude (degrees, north)) + LB=DATA.DATA(j).GeoLong(k); % 2. Географическая долгота (градусы, восток)(Geographic longitude (degrees, east)) + R=mydistance(fiA,LA,fiB,LB); + if R<=R_quake + HHeden(end+1)= mymean(DATA.DATA(j).RPA_HHeden((k-1)*5+1:(k-1)*5+5),1e10); + O2den(end+1) = mymean(DATA.DATA(j).RPA_O2den((k-1)*5+1:(k-1)*5+5),1e10); + Eden(end+1) = mymean(DATA.DATA(j).EP_Eden((k-1)*5+1:(k-1)*5+5),1e10); + Etemp(end+1) = mymean(DATA.DATA(j).EP_Etemp((k-1)*5+1:(k-1)*5+5),1e10); end end end end - report(end).RPA_HHeden=mymean(HHeden,1e10); - report(end).RPA_O2den=mymean(O2den,1e10); - report(end).EP_Eden=mymean(Eden,1e10); - report(end).EP_Etemp=mymean(Etemp,1e10); + num=give_number(report,param.sattelite_range,R_quake); + report.data(num).N_sat=param.sattelite_range; + report.data(num).R=R_quake; + report.data(num).date{i}=[myformat(param.dates(i).day,2),'.',... + myformat(param.dates(i).month,2),'.',... + myformat(param.dates(i).year,4)]; + report.data(num).RPA_HHeden(i)=mymean(HHeden,1e10); + report.data(num).RPA_O2den(i)=mymean(O2den,1e10); + report.data(num).EP_Eden(i)=mymean(Eden,1e10); + report.data(num).EP_Etemp(i)=mymean(Etemp,1e10); end + param.report=report; +save(fname_rep,'report'); +end + +function num=give_number(report,f,R) +flag=0; % Проверяет, есть ли в отчете данные со спутника f +[~,b]=size(report.data); + +for i=1:b + if report.data(i).N_sat==f + if report.data(i).R==R + flag=1; + num=i; + end + end +end +if flag==0 + num=b+1; +end end function DATA=load_file(sattelite, year, month, day, type_sensor) +success=0; foldername=['DATA/f' num2str(sattelite) '/' type_sensor '/' num2str_new(year,4) '/' num2str_new(month,2)]; target=ls(foldername); [str row]=size(target); @@ -82,7 +111,7 @@ for i=1:str counter=counter+1; end end - filename=filename(1:end-counter); + filename=filename(1:end-counter); if length(filename)>7 keyword=[num2str_new(year,4) num2str_new(month,2) num2str_new(day,2)]; flag=0; @@ -90,11 +119,18 @@ for i=1:str if filename(j:j+length(keyword)-1)==keyword flag=1; end - - end - if (prod(filename(end-3:end)=='.mat')==1)&&(flag==1) - DATA = load([foldername,'/',filename], 'DATA'); % открыть для чтения + end + if (prod(filename(end-3:end)=='.mat')==1)&&(flag==1) + disp(['Найден подходящий файл: ',filename]) + DATA = load([foldername,'/',filename], 'DATA'); % открыть для чтения + success=1; + end end end +if ~success + disp('Нужных файлов не найдено, перехожу к загрузке') + loaddata_edr(sattelite, year, month, type_sensor); + DATA=load_file(sattelite, year, month, day, type_sensor); +end end diff --git a/loaddata_edr.m b/loaddata_edr.m index 5fc835e..6fc35f8 100644 --- a/loaddata_edr.m +++ b/loaddata_edr.m @@ -1,49 +1,19 @@ -function loaddata_edr %DATA=loaddata_edr(filename) +function loaddata_edr(sattelite, year, month, type_sensor) %DATA=loaddata_edr(filename) %% By L_DelOff % Загрузка данных из файла *.EDR % Описание файла приведено в: % !docs/!AFRL ASCII and Binary File Format Descriptions.pdf %% Поехали warning off -sattelite_range=15:18; % Выбор с какого спутника требуются данные -year_range=2015; % -month_range=9; % %% Скачать данные в формате архива -if 1 - for sattelite=sattelite_range %6:20 - for year=year_range % года - for month=month_range % месяца - download_fun(sattelite, year, month); %скачивание и обработка нужных данных - end - end - end -end +download_fun(sattelite, year, month, type_sensor); %скачивание и обработка нужных данных %% Распаковать архив -if 1 - for sattelite=sattelite_range %6:20 - for year=year_range % года - for month=month_range % месяца - extract_data_fun(sattelite, year, month); %скачивание и обработка нужных данных - end - end - end -end +extract_data_fun(sattelite, year, month, type_sensor); %скачивание и обработка нужных данных %% Преобраховать EDR в mat -if 1 - for sattelite=sattelite_range %6:20 - for year=year_range % года - for month=month_range % месяца - convert_data_fun(sattelite, year, month); %скачивание и обработка нужных данных - end - end - end +convert_data_fun(sattelite, year, month, type_sensor); %скачивание и обработка нужных данных end - -end - -function convert_data_fun(sattelite, year, month) -type_sensor='ssies'; +function convert_data_fun(sattelite, year, month, type_sensor) foldername=['DATA/f' num2str(sattelite) '/' type_sensor '/' num2str_new(year,4) '/' num2str_new(month,2)]; target=ls(foldername); [str row]=size(target); @@ -91,8 +61,7 @@ for i=1:str end end -function extract_data_fun(sattelite, year, month) -type_sensor='ssies'; +function extract_data_fun(sattelite, year, month, type_sensor) foldername=['DATA/f' num2str(sattelite) '/' type_sensor '/' num2str_new(year,4) '/' num2str_new(month,2)]; target=ls(foldername); [str row]=size(target); @@ -107,10 +76,9 @@ for i=1:str end end -function download_fun(sattelite, year, month) +function download_fun(sattelite, year, month, type_sensor) % https://satdat.ngdc.noaa.gov/dmsp/data/f18/ssies/2020/05/ % https://satdat.ngdc.noaa.gov/dmsp/data/f18/ssies/2020/05/PS.CKGWC_SC.U_DI.A_GP.SIES3-F18-R99990-B9999090-APGA_AR.GLOBAL_DD.20200501_TP.000001-235959_DF.EDR.gz -type_sensor='ssies'; message=['Ищем данные спутника F' num2str(sattelite) ' по дате:' num2str_new(year,4) '/' num2str_new(month,2) '...']; fprintf(message) url_site=['https://satdat.ngdc.noaa.gov/dmsp/data/f' num2str_new(sattelite,2) '/' type_sensor '/' num2str_new(year,4) '/' num2str_new(month,2) '/']; @@ -157,43 +125,6 @@ while i