Login or Register to Ask a Question and Join Our Community

Sponsored Content
Full Discussion: Shell or awk script to compute average of all the points within a circle
Top Forums Shell Programming and Scripting Shell or awk script to compute average of all the points within a circle Post 302960899 by Indra2011 on Saturday 21st of November 2015 10:32:30 AM
Old 11-21-2015
Shell or awk script to compute average of all the points within a circle
HI Help,

I have a file which looks like below ---

Input file --->

Code:
1970113.00000 3460.00000 1.09516 
1970116.00000 3791.00000 1.06350 
1970120.00000 4120.00000 1.07588 
1970115.00000 4450.00000 1.09591 
1970116.00000 4780.00000 1.09965 
1970120.00000 5109.00000 1.06733 
1970122.00000 5440.00000 1.03760 
1970124.00000 5770.00000 1.02025 
1970123.00000 6100.00000 1.00998 
1970120.00000 6430.00000 0.96426

What I want to do ? ----
Code:
For each line (#NR= 1 to 10) I would like to search X($1),Y($2) within the I/P file  which lies within a radius of 50 and then average out all the $3 within that circle of radius 50.

Problem ??

My script below does the job but it takes ages to go through the input file which is TOO large.

Code:
 
set first = 1    # NR=1 ; first line fo the input file
set last = 10   # NR = 10; Last line of the input file
 
set num = ${first}
while (${num} <= ${last})
 
set X0 = `cat S2 | awk -v line=${num} 'NR==line' | awk '{print $1}'` # $1,X, for each Record
set Y0 = `cat S2 | awk -v line=${num} 'NR==line' | awk '{print $2}'` # $2,Y, for each Record
set XS = `cat S2 | awk -v line=${num} 'NR==line' | awk '{print $3}'` # $3,Value that will be averaged out, for each Record
set R0 = `echo "50"` # Search radius
 
set AVG = `cat S2 | awk -v X=${X0} -v Y=${Y0} -v R=${R0} '{print $1,$2,$3,sqrt(($1-X)*($1-X) + ($2-Y)*($2-Y)) <=R}'  | awk '{if($4==1){print $0}}' | awk '{ sum+=$3} END {print sum/NR}'`# Search all the points within the input file which lies within R #Average all the $3,values#
 
echo "${X0} ${Y0} ${XS} ${AVG}" >> TMP #For each record print a new line with existing $1,$2,$3 and $AVG
 
@ num++ # end of loop and it goes to next NR
end

This could be very easy for you experts.

Thanks,
Last edited by Indra2011; 11-22-2015 at 11:29 AM..
 

10 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

Shell Script to see the mount points.

Hi all, First of all I dont even know the ABC of scripting .. But now I want a Script to see the mount points of the file systems Can any body help plsssssssss :o (1 Reply)
Discussion started by: priky
1 Replies

2. Shell Programming and Scripting

Compute the median of a set of numbers with AWK?

Is there a way in awk to compute the median of a set of numbers in a file in the following format. 34 67 78 100 23 45 67 (3 Replies)
Discussion started by: Lucky Ali
3 Replies

3. Shell Programming and Scripting

shell script for finding average runtime of other script

so I've made a shell script that downloads 6 files in succession from a given url, then deletes them. Now I want to time the script, and the average time it uses by running it ~100 times. My problem is tho, how do I store the time it takes for each run through of the script? I know time writes to... (3 Replies)
Discussion started by: navlelo
3 Replies

4. Homework & Coursework Questions

Shell Script average runtime

Use and complete the template provided. The entire template must be completed. If you don't, your post may be deleted! 1. The problem statement, all variables and given/known data: Make a bash script that calculates average runtime for the first two scripts you made. The average should be... (17 Replies)
Discussion started by: navlelo
17 Replies

5. Shell Programming and Scripting

awk based script to find the average of all the columns in a data file

Hi All, I need the modification for the below mentioned code (found in one more post https://www.unix.com/shell-programming-scripting/27161-script-generate-average-values.html) to find the average values for all the columns(but for a specific rows) and print the averages side by side. I have... (4 Replies)
Discussion started by: ks_reddy
4 Replies

6. Shell Programming and Scripting

AWK script to split data and find average

Input: 2.58359023380340e+02 1.43758864405595e+02 -7.65700666212508e+00 1.06460208083228e+02 1.26185441783936e+02 -3.41389169427027e+01 -1.40393299309592e+02 -3.07758776849508e+01 1.45067703495838e+02 1.79405834959073e+02 5.06666234594205e+01 OUT 2.0105894389e+02 (average of... (8 Replies)
Discussion started by: chrisjorg
8 Replies

7. Shell Programming and Scripting

Compute average ignoring outliers of different segments within a dat file using awk

I have data files that look like this, say data.txt 0.00833 6.34 0.00833 6.95 0.00833 7.08 0.00833 8.07 0.00833 8.12 0.00833 8.26 0.00833 8.70 0.00833 9.36 0.01667 20.53 0.01667 6.35 0.01667 6.94 0.01667 7.07 0.01667 8.06 0.01667 8.10 0.01667 8.25 0.01667 8.71 0.01667 9.31... (7 Replies)
Discussion started by: malandisa
7 Replies

8. Shell Programming and Scripting

Linux Mount Points usage check with shell script

Hi Everyone, Just in need of your help again, I have managed to get a script to check the linux disk usage stuff. But wanted to tweak it little more to get the desired output. Requirement: Now its querying only one mount point, As its not saving in array instead calling it as variables. So... (1 Reply)
Discussion started by: thiyagoo
1 Replies

9. Shell Programming and Scripting

Compute average based on field values

Im looking for a way to average the values in field 14 (when field 2 is equal to 2016) and fields 3 and 4 (when field 2 is equal to 2017). Any help is appreciated. 001001 2016 33.22 38.19 48.07 51.75 59.77 67.68 70.86 72.21 66.92 53.67 42.31 40.15 001001 2017 ... (10 Replies)
Discussion started by: ncwxpanther
10 Replies

10. UNIX for Beginners Questions & Answers

awk or sed script to count number of occurrences and creating an average

Hi Friends , I am having one problem as stated file . Having an input CSV file as shown in the code U_TOP_LOGIC/U_HPB2/U_HBRIDGE2/i_core/i_paddr_reg_2_/Q,1,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,0,0... (4 Replies)
Discussion started by: kshitij
4 Replies

LEARN ABOUT DEBIAN

measurefs.reiser4

measurefs.reiser4(8)						reiser4progs manual					      measurefs.reiser4(8)

NAME

       measurefs.reiser4 - the program for measuring reiser4 filesystem parameters (fragmentation, node packing, etc.).

SYNOPSIS

       measurefs.reiser4 [ options ] FILE

DESCRIPTION

       measurefs.reiser4  is  reiser4  filesystem  measure program. You can estimate reiser4 filesystem fragmentation, packingm etc. structures by
       using it.

COMMON OPTIONS

       -V, --version
	      prints program version.

       -?, -h, --help
	      prints program help.

       -y, --yes
	      assumes an answer 'yes' to all questions.

       -f, --force
	      forces measurefs to use whole disk, not block device or mounted partition.

       -c, --cache N
	      sets tree cache node number to passed value. This affects very much behavior of libreiser4. It affects speed, tree allocation, etc.

MEASUREMENT OPTIONS

       -S, --tree-stat
	      shows different tree statistics (node packing, internal nodes, leaves, etc)

       -T, --tree-frag
	      measures total tree fragmentation. The result is fragmentation factor - value from 0.00000 (minimal fragmentation) to 1.00000 (maxi-
	      mal one). Most probably, this factor may affect sequential read performance.

       -D, --data-frag
	      measures	average  files	fragmentation.	This  means, that fragmentation of each file in filesystem will be measured separately and
	      results will be averaged. The result is fragmentation factor - value from 0.00000 (minimal fragmentation) to 1.00000 (maximal one).

	      Note, that for the fresh filesystem (created not very long time ago) and even fully filled by data, this value will be pretty small.

       -F, --file-frag FILE
	      measures fragmentation of the specified file. The result is fragmentation factor - value from  0.00000  (minimal	fragmentation)	to
	      1.00000  (maximal one). Note, that fragmentation of a small file (depends of used tail policy), which consists of tail items, is not
	      very reliable value. That is because, they is always afoot due to balancing.

	      Examples:

	      measurefs.reiser4 -F /usr/bin /dev/hda2
	      measurefs.reiser4 -F /bin/bash /dev/hda2

       -E, --show-file
	      show file fragmentation for each file if --data-frag is specified.

PLUGIN OPTIONS

       -p, --print-profile
	      prints the plugin profile. This is the set of default plugins used for all parts of a filesystem -- format, nodes,  files,  directo-
	      ries, hashes, etc. If --override is specified, then prints modified plugins.

       -l, --print-plugins
	      prints all plugins libreiser4 know about.

       -o, --override TYPE=PLUGIN, ...
	      overrides the default plugin of the type "TYPE" by the plugin "PLUGIN" in the plugin profile.

	      Examples:

	      measurefs.reiser4 -o nodeptr=nodeptr41,hash=rupasov_hash /dev/hda2

REPORTING BUGS

       Report bugs to <reiserfs-devel@vger.kernel.org>

SEE ALSO

       debugfs.reiser4(8), mkfs.reiser4(8), fsck.reiser4(8)

AUTHOR

       This manual page was written by Yury Umanets <umka@namesys.com>

reiser4progs							   28 Apr, 2003 					      measurefs.reiser4(8)
All times are GMT -4. The time now is 12:06 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy