-
Notifications
You must be signed in to change notification settings - Fork 2
/
UWBRADEstimationMain.m
44 lines (34 loc) · 1.09 KB
/
UWBRADEstimationMain.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
% Script to integrate the work flow of UWBRAD temperature estimation
% Yuna Duan March 24th based on Mike's seprated scripts
%% 1. Get UWBRAD Tb
% Average MEMLS Tb by antenna pattern file using function GetModelOutput
% Save result in dat/TbUWBRAD-v2
% Reorganize the data in the way shown in the numbers sheet
% Save the reorganized results in TbUWBRAD
clear
RecalculateTb=true;
if RecalculateTb
GetModelOutput;
ReOrganize;
end
%% 2. Generate virtual observation
% Choose random realizations of all parameters, and generate a
% set of observations along the flight lines
RenerateObs=true;
if RenerateObs
GenerateObs;
end
%% 3. Estimate Temperature
% Calculate MCMC estimate of temperature parameters (~4 minute
% run-time per location). Saves dat/Expi.mat, where i=0,1,2...
% Change ExpNum and interation# N to adjust the estimation
x=15; %location 1-47 to estimate
ExpNum=7;
N=1E4; %MCMC iterations
O=load('dat/TbObs.mat');
D=load('dat/TbUWBRAD.mat');
EstimateT;
%% 4.PostProcess
% Explores the MCMC output via figures, including the final
% estimated temperature profiles
PostProcess;