============================================================== MEASURE_PDF: Measure the "goodness" of a group-wise PDF GENERAL score = measure_pdf(input_pdf, pdf_evaluation_method) INPUT/S -input pdf: -pdf_evaluation_method: OUTPUT/S -score: PENDING WORK - KNOWN BUG/S -None. COMMENT/S RELATED FUNCTION/S ABOUT -Created: February 2004 -Last update: February 2004 -Revision: 0.0.1 -Author: R. S. Schestowitz, University of Manchester ==============================================================
0001 function score = measure_pdf(input_pdf, pdf_evaluation_method, warped_image) 0002 % ============================================================== 0003 % MEASURE_PDF: Measure the "goodness" of a group-wise PDF 0004 % 0005 % GENERAL 0006 % 0007 % score = measure_pdf(input_pdf, pdf_evaluation_method) 0008 % 0009 % INPUT/S 0010 % 0011 % -input pdf: 0012 % 0013 % -pdf_evaluation_method: 0014 % 0015 % OUTPUT/S 0016 % 0017 % -score: 0018 % 0019 % PENDING WORK 0020 % 0021 % - 0022 % 0023 % KNOWN BUG/S 0024 % 0025 % -None. 0026 % 0027 % COMMENT/S 0028 % 0029 % 0030 % RELATED FUNCTION/S 0031 % 0032 % 0033 % 0034 % ABOUT 0035 % 0036 % -Created: February 2004 0037 % -Last update: February 2004 0038 % -Revision: 0.0.1 0039 % -Author: R. S. Schestowitz, University of Manchester 0040 % ============================================================== 0041 0042 accumulator = 0; 0043 0044 % initialise it to zero 0045 0046 if (strcmp(pdf_evaluation_method,'TFC Suggestion')), 0047 for i = 1:size(input_pdf,1), 0048 current_value = input_pdf(i,:); 0049 accumulator = accumulator - log(warped_image(i)); 0050 % still need to use the above current_value in some way 0051 end 0052 elseif(strcmp(pdf_evaluation_method,'Default')), 0053 for i = 1:size(input_pdf,1), 0054 current_value = input_pdf(i,:); 0055 accumulator = accumulator + log(mean(current_value)); 0056 end 0057 else 0058 error('Unknown PDF evaluation method.'); 0059 end 0060 0061 score = - accumulator;