-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
In this version, we - update the evaluation code - add annotations to map the steps in Figure 1 to different code sections in FindPolynyaMain.m
- Loading branch information
Showing
19 changed files
with
1,354 additions
and
854 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,144 @@ | ||
| function In = NameList | ||
| % Name list of main function | ||
| % you can run the function by FindPolyunyaMain.m | ||
| % Var Name | ||
| % TIME | ||
| % TimeTotal total time needs to detect polynyas | ||
| % StartTime time start to detect polynyas, default is | ||
| % the fisrt of TimeTotal NewYear | ||
| % READ SIC FILES | ||
| % SICFile a struct storing the information of SIC | ||
| % | files need to be read | ||
| % |-SICFile.Dir the path of the folder storing SIC files | ||
| % |-SICFile.Name1 name of SIC files before the timestamp | ||
| % |-SICFile.Name2 name of SIC files after the timestamp | ||
| % |-SICFile.LandMask land mask for SIC files | ||
| % SICLon longitude for SIC | ||
| % SICLat latitude for SIC | ||
| % DETECT OPEN SEA & ICE | ||
| % Lim the threshold for the openwater/ice | ||
| % LASTING OPEN WATER | ||
| % SeriesLength the length of time series to detect lasting | ||
| % openwater | ||
| % FrequencyThreshold the threshold for detecting lasting | ||
| % openwater | ||
| % MapRange remapping to get membership of openwater | ||
| % WARM SEASON | ||
| % HeatLossFlag is the HeatLoss Mode for judging warm | ||
| % season on | ||
| % RESTART | ||
| % RestartDir path of the folder to store restart files | ||
| % RestartStride the stride to store a restart file | ||
|
|
||
| %% | ||
| clear all; | ||
| %% Time | ||
| In.TimeTotal = ... The time you want to identify polynyas | ||
| [datetime('2003-02-01') : datetime('2003-10-15'), ... | ||
| datetime('2004-02-01') : datetime('2011-08-31'), ... | ||
| datetime('2012-08-01') : datetime('2023-01-31')]; | ||
| % In.TimeTotal = ... The time you want to identify polynyas | ||
| % datetime('2017-02-01') : datetime('2018-01-31'); | ||
| TimeTotalmmdd = str2double(string(datestr(In.TimeTotal, 'mmdd'))); | ||
| In.TimeTotal = In.TimeTotal(TimeTotalmmdd >= 0401 & TimeTotalmmdd <= 1031); | ||
| % In.TimeTotal = [In.TimeTotal(1 : 2813), datetime('2017-11-01') : datetime('2017-11-15'), In.TimeTotal(2814 : end)]; | ||
| % datetime('2014-02-01') : datetime('2015-01-31'); | ||
|
|
||
| In.TimeGap = [datetime('2003-12-01'), datetime('2012-07-01')]; | ||
| In.StartTime = In.TimeTotal(1); | ||
| In.NewYear = '02-01'; | ||
|
|
||
| %% Read Data | ||
| % we will read data as [SICDir \ SICFileName1 Timestr SICFileName2] | ||
| In.SICFile.Dir = ... | ||
| 'G:\Antaratica_ASI_SIC_6250'; | ||
| In.SICFile.Name1 = ...% {'asi-AMSR2-s3125-'}; | ||
| {'asi-s6250-', 'asi-s6250-', 'asi-AMSR2-s6250-'}; | ||
| In.SICFile.Name2 = ...% {'-v5.4.hdf'}; | ||
| {'-v5.4.hdf', '-v5.4.hdf', '-v5.4.hdf'}; | ||
| In.SICFile.VarName = ... | ||
| 'ASI Ice Concentration'; | ||
| % read the Lon & Lat | ||
| In.SICLon = hdfread(... | ||
| [In.SICFile.Dir, '\LongitudeLatitudeGrid-s6250-Antarctic.hdf'], 'Longitudes'); | ||
| In.SICLat = hdfread(... | ||
| [In.SICFile.Dir, '\LongitudeLatitudeGrid-s6250-Antarctic.hdf'], 'Latitudes'); | ||
| % In.SICLon = load('G:\Antaratica_ASI_SIC_6250\LonLat_25000.mat'); | ||
| % In.SICLat = In.SICLon.Lat; | ||
| % In.SICLon = In.SICLon.Lon; | ||
| % read land mask | ||
| In.SICFile.LandMask = ncread([In.SICFile.Dir, '\LandMaskAMSR6250.nc'], 'LandMask'); | ||
| % set the resolution | ||
| In.Resolution = 6.25; | ||
|
|
||
| %% Land Fast Ice | ||
| In.FastIceFlag = FastIceParameter; | ||
|
|
||
| %% Detect Opensea & Ice | ||
| In.Lim = ... The threshold for the openwater/ice | ||
| 30; | ||
|
|
||
| %% Lasting Open Water | ||
| In.SeriesLength = ... The length of identify time series | ||
| 14; | ||
| % the FrequencyThreshold is calculated as FrequencyThreshold ./ SeriesLength | ||
| In.FrequencyThres = ... the threshold for frequent openwater. it is a 1*2 vector | ||
| [0.7, 0.5]; | ||
|
|
||
| In.MapRange = ... Remapping for membership of openwater | ||
| [30, 30]; | ||
|
|
||
| %% Match | ||
| global IDCpacity | ||
| IDCpacity = ... It is a constant for OverlapDye.m. If it | ||
| 1000000; | ||
|
|
||
| %% Physical ID to Logical ID | ||
| global ReincarnationTol | ||
| ReincarnationTol = ... The tolerance for the Reincarnation | ||
| 4; | ||
| In.CombineMergeThres = 0.5; | ||
| In.MinPolynyaArea = 3; | ||
|
|
||
| %% Rebirth | ||
| In.RebirthOverlapThres = ... | ||
| 0.2; | ||
| In.RebirthOverlapThresYear = ... | ||
| 0.4; | ||
| In.TimeFilterAfter = ... | ||
| 20; | ||
|
|
||
| %% Warm Season | ||
| % In.HeatLossFlag = HeatLossParms; | ||
|
|
||
| %% Air Temperature Judge | ||
| In.TempeJudgeFlag = AirTemperatureParameter; | ||
|
|
||
| %% Cross Year Track | ||
| In.CrossYearOverlapThres = ... | ||
| 0.2; | ||
| In.SeriesLengthThres = ... | ||
| 0.3; | ||
| In.SeriesLengthThresYear = ... | ||
| 0.2; | ||
|
|
||
| %% Restart | ||
| In.RestartDir = ... | ||
| 'C:\Users\13098\Documents\冰间湖识别\Data\tempData'; | ||
| In.RestartStride = ... | ||
| 120; | ||
|
|
||
| %% CacheFile | ||
| In.Cache = ... | ||
| 'C:\Users\13098\Documents\冰间湖识别\Data\tempData'; | ||
|
|
||
| %% SaveResults | ||
| In.Save.Path = ... | ||
| 'G:\DEEP-AAShare\SIC30_6.25km_20d'; | ||
| In.Save.FileName1 = ... | ||
| 'DEEP_s6250_AMSR_SIC_'; | ||
| In.Save.FileName2 = ... | ||
| '_v1.0.nc'; | ||
|
|
||
| end | ||
| % run the function by FindPolyunyaMain.m |
Oops, something went wrong.