1 function [ log_sum ] = new_objective( w, red, g )
\r
2 %[ likelihood ] = new_objective( w, red, g )
\r
3 %returns data log-likelihood for s specified ignition point
\r
5 % w : processed wrf-out information
\r
6 % red : sturct coming from
\r
7 % g : struct coming from load_subset_detections
\r
9 % likelihood : likelihood of data given ignition pt
\r
12 radius = round(sig*sqrt(32*log(10)))+1;
\r
13 weight = gauss_weight(sig);
\r
16 % loop over appropriate subset of detections
\r
19 heat = make_heats(red.tign_g,g(k).time);
\r
20 mask = get_fire_pixels(g(k).fxdata);
\r
21 dp_matrix = make_detection_probabilities(heat);
\r
22 %skip if no detections in
\r
24 [mm, nn] = size(mask);
\r
25 % loop over detection pixels in mask
\r
26 for ii = radius:mm-radius
\r
27 for jj =radius:nn-radius
\r
29 pp = compute_pixel_probability(ii,jj,heat,sig,weight,dp_matrix);
\r
30 like = like+log(pp);
\r
32 % pp = 1 - compute_pixel_probability(ii,jj,heat,radius,weight);
\r
33 % like = like+log(pp);
\r
38 log_sum = log_sum + like;
\r