This file is indexed.

/usr/bin/pfsstat is in pfstools 2.1.0-3.

This file is owned by root:root, with mode 0o755.

The actual contents of the file can be viewed below.

 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
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
#!/usr/bin/octave -q
#  
# This file is a part of PFSTOOLS package.
#  ---------------------------------------------------------------------- 
#  Copyright (C) 2003,2004 Rafal Mantiuk and Grzegorz Krawczyk
#  
#   This program is free software; you can redistribute it and/or modify
#   it under the terms of the GNU General Public License as published by
#   the Free Software Foundation; either version 2 of the License, or
#   (at your option) any later version.
# 
#   This program is distributed in the hope that it will be useful,
#   but WITHOUT ANY WARRANTY; without even the implied warranty of
#   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
#   GNU General Public License for more details.
# 
#   You should have received a copy of the GNU General Public License
#   along with this program; if not, write to the Free Software
#   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
#  ---------------------------------------------------------------------- 
# 
#  @author Rafal Mantiuk, <mantiuk@mpi-sb.mpg.de>
# 
#  $Id: pfsstat,v 1.2 2005/11/07 14:40:52 rafm Exp $
#
# See man page for more information

pin = pfsopen( "stdin" );

#fprintf( stderr, "l = %d\n", length( argv ) );

if( nargin != 0 )
  error( "Expecting no parameters" ); 
endif

if( 2 != exist( 'gausswin' ) )
  error( "%s requires the 'signal' package - please install octave-signal\n", \
      program_name () );
endif

firstFrame = true;
while( true )
  pin = pfsget( pin );

  if( pin.EOF == true ) # Are there any more frames
    break;
  endif

  ## Low pass filter to eliminate influence of high frequencies on
  ## the statistic
  kernel = gausswin( 30, 2.2 );
  ##*gausswin( 30, 2.2 )';
  kernel = kernel / sum( kernel(:) );
  ## Separate filtering for rows and cols
  fy = conv2( pin.channels.Y, kernel, 'same' );
  fy = conv2( fy, kernel', 'same' );

  ## Remove zero and negative values - set them to smalles possitive value
  min_positive = min( vec( fy + (fy<=0)*1000 ) );
  fy = fy .* (fy >= min_positive) + min_positive .* (fy < min_positive);

  ## Work in log domain
  lfy = log10(fy);

  ## Calculate statistics
  minimum = min( vec( lfy ) );
  maximum = max( vec( lfy ) );
  dynamic_range = maximum - minimum;
  s_average = mean( vec( lfy ) );
  s_median = median( vec( lfy ) );

  if( !firstFrame )
    fprintf( stderr, "======================\n" );
  endif
  if( isfield( pin.tags, "FILE_NAME" ) )
    fprintf( stderr, "File: %s\n", pin.tags.FILE_NAME );
  endif
  if( isfield( pin.tags, "LUMINANCE" ) )
    fprintf( stderr, "Luminance profile: %s\n", pin.tags.LUMINANCE );  
  endif     
  fprintf( stderr, "Width: %d Height: %d\n", pin.columns, pin.rows );
  fprintf( stderr, "Dynamic Range: \t%g [log_10]\t 1:%g [lin]\t %g [dB]\n", ...
                                dynamic_range, 10^dynamic_range, ...
                                round(20*dynamic_range)  );
  fprintf( stderr, "Minimum: \t%g [log_10]\t %g [lin]\n", ...
                                minimum, 10^minimum  );
  fprintf( stderr, "Maximum: \t%g [log_10]\t %g [lin]\n", ...
                                maximum, 10^maximum  );

  fprintf( stderr, "Log average: \t%g [log_10]\t %g [lin]\n", ...
                                s_average, 10^s_average  );

  fprintf( stderr, "Median: \t%g [log_10]\t %g [lin]\n", ...
                                s_median, 10^s_median  );

  firstFrame = false;
endwhile

pfsclose( pin );