Shell or awk script to compute average of all the points within a circle Post: 302960960

Sponsored Content

Top Forums Shell Programming and Scripting Shell or awk script to compute average of all the points within a circle Post 302960960 by RudiC on Sunday 22nd of November 2015 02:05:57 PM

11-22-2015

Registered User

I'm not sure I interpreted your requirements correctly. Your sample file doesn't have any line's coordinates within radius 50 from any other, so any test is impossible for that set (it did some avaraging for R = 1500).
Try

Code:

awk -vR0=50 '
        {X[NR]=$1
         Y[NR]=$2
         V[NR]=$3
        }
END     {for (n=1; n<=NR; n++)  {X0 = X[n]
                                 Y0 = Y[n]
                                 SUM = CNT = 0
                                 for (i=1; i<=NR; i++)  {R = sqrt((X[i]-X0)*(X[i]-X0) + (Y[i]-Y0)*(Y[i]-Y0))
                                                         if (R<=R0)     {SUM += V[i]
                                                                         CNT++
                                                                        }
                                                        }
                                 print X0, Y0, V[n], SUM/CNT
                                }
        }
' file
1970113.00000 3460.00000 1.09516 1.09516
1970116.00000 3791.00000 1.06350 1.0635
1970120.00000 4120.00000 1.07588 1.07588
1970115.00000 4450.00000 1.09591 1.09591
1970116.00000 4780.00000 1.09965 1.09965
1970120.00000 5109.00000 1.06733 1.06733
1970122.00000 5440.00000 1.03760 1.0376
1970124.00000 5770.00000 1.02025 1.02025
1970123.00000 6100.00000 1.00998 1.00998
1970120.00000 6430.00000 0.96426 0.96426

Not sure how this approach would handle 80000 lines, though...

---------- Post updated at 20:05 ---------- Previous update was at 20:00 ----------

I have to correct myself - that won't be ALLINES * ALLINES result lines but ALLINES result lines, still ALLINES * ALLINES computations...

These 2 Users Gave Thanks to RudiC For This Post:

RudiC

View Public Profile for RudiC

Find all posts by RudiC

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

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. 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...

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...

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...

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...

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...

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...

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. 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...

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)