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

input

INPUT(5)							File Formats Manual							  INPUT(5)

NAME

       input - input data file format

DESCRIPTION

       This  manual  page  describes  the  output format of the file created by the Yagi-Uda project input or first programs. The file is an ASCII
       file, so can be edited manually if required. It is read by yagi or optimise.  The program optimise also creates the same file format.

       NOTES
       FREQUENCY
       MIN_FREQUENCY
       MAX_FREQUENCY
       STEP_FREQUENCY
       ELEMENTS
       DRIVEN
       PARASITIC
       ANGULAR_STEP
       DATA_DRIVEN
       DATA_PARASITIC

       The
       #
       is used as a comment.

       A typical file, for an antenna designed for the 2m amateur band (144-146
       MHz) might be as shown below. Note:
       1) Reflector is at position 0.0
       2) All elements lengths/positions are in m.

       NOTES This line can be up to 400 characters.
       FREQUENCY 145.000000
       MIN_FREQUENCY 144.000000
       MAX_FREQUENCY 146.000000
       STEP_FREQUENCY 0.100000
       ELEMENTS 10
       DRIVEN 1
       PARASITIC 9
       ANGULAR_STEP   180.000000
       #DATA_DRIVEN x y length	diameter voltage phase
       DATA_DRIVEN  0.30993 0.000000 1.002492 0.0063000000000000 1.000000 0.000000
       # DATA_PARASITIC x y length  diameter
       DATA_PARASITIC
	  0.00000    0.00000	1.01902    0.00630
	  0.49859    0.00000	0.95551    0.00630
	  1.07338    0.00000	0.93595    0.00630
	  1.81020    0.00000	0.92050    0.00630
	  2.58152    0.00000	0.91322    0.00630
	  3.30390    0.00000	0.91243    0.00630
	  4.03692    0.00000	0.90423    0.00630
	  4.64108    0.00000	0.89081    0.00630
	  5.37418    0.00000	0.83912    0.00630

FILES

	    filename				ASCII description of antenna.	    filename.out			binary file

SEE ALSO

       yagi(1), output(1), input(1), optimise(1), first(1) and yagi(5).

AUTHORS

       Dr. David Kirkby G8WRB (david.kirkby@onetel.net), with help with converting to DOS from Dr. Joe Mack NA3T (mack@fcrfv2.ncifcrf.gov).

Yagi version 1.16						 24th October 2000							  INPUT(5)
All times are GMT -4. The time now is 01:24 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy