Login or Register to Ask a Question and Join Our Community


UNIX for Dummies Questions & Answers

How to filter out almost dupicate X Y (Easting Northing) coordinates?

 
Thread Tools Search this Thread
Top Forums UNIX for Dummies Questions & Answers How to filter out almost dupicate X Y (Easting Northing) coordinates?
# 8  
Old 09-22-2010
Gentlemen,

I thank you very much for your efforts in try to helping me with this code.

As I said, note that I am working in a Microsoft Windows XP command prompt window environment. Not sure if that will have any effect on the syntax.

I put your code into a command text file.
I can't seem to get the code to work without errors when listed vertically, so
I had to put it all on one line, and then type.....
mawk -v range=5 -f awk_commands.txt FileIn.txt

But......

All it does is print the complete contents of the input file (plus that extra "secret" column) to the screen, and does not create the duplicate and unique files.

The other strange behaviour is that no matter what value I use for range, a negative number, 0, or 1000000, I get the same output:

PointNumber_0000001 1000000.123456 1000000.123456 10000.123456 1000001.75
PointNumber_0000010 1000001.234567 1000002.234567 10345.234567 1000001.75
PointNumber_0000100 1000010.345678 1000020.456789 10030.987654 1000012.75
PointNumber_0001000 1000050.345678 1000050.456789 10030.987654 1000051.25
PointNumber_0010000 1000123.123456 1000456.123456 10789.123456 1000125.50
PointNumber_0100000 1000123.123456 1000456.123456 10789.123456 1000125.50
PointNumber_1000000 1000000.123456 1000000.123456 10000.123456 1000001.75
PointNumber_2000000 1000011.345678 1000021.456789 10030.987654 1000012.75
PointNumber_3000000 1000051.000678 1000049.999000 10030.987654 1000051.25[/SIZE][/SIZE]
 
Login or Register to Ask a Question

Previous Thread | Next Thread

10 More Discussions You Might Find Interesting

1. UNIX for Beginners Questions & Answers

Help with processing coordinates in a file.

I have a variation table (variation.txt) which is a very big file. The first column in the chromosome number and the second column is the position of the variation. I have a second file annotation.txt which has a list of 37,000 genes (1st column), their chromosome number(2nd column), their start... (1 Reply)
Discussion started by: Sanchari
1 Replies

2. UNIX for Dummies Questions & Answers

overlapped genomic coordinates

Hi, I would like to know how can I get the ID of a feature if its genomic coordinates overlap the coordinates of another file. Example: Get the 4th column (ID) of this file1: chr1 10 100 gene1 chr2 3000 5000 gene2 chr3 200 1500 gene3 if it overlaps with a feature in this file2: chr2... (1 Reply)
Discussion started by: fadista
1 Replies

3. UNIX for Dummies Questions & Answers

Length of a segment based on coordinates

Hi, I would like to have the length of a segment based on coordinates of its parts. Example input file: chr11 genes_good3.gtf aggregate_gene 1 100 gene1 chr11 genes_good3.gtf exonic_part 1 60 chr11 genes_good3.gtf exonic_part 70 100 chr11 genes_good3.gtf aggregate_gene 200 1000 gene2... (2 Replies)
Discussion started by: fadista
2 Replies

4. Shell Programming and Scripting

Differential substring removal using coordinates

Hello all, this might be better suited for a bioinformatics forum, but I thought I'd try my luck here as well. I have several tabular text files of DNA sequence reads that appear as such: File_1.txt >H01BA45XW GATTACAGATTCGACATCCAACTGAGGCATT >H02BG78WR CCTTACAGACTGGGCATGAATATTGCATACC... (3 Replies)
Discussion started by: vectorborne5
3 Replies

5. Shell Programming and Scripting

Determination n points between two coordinates

Hi guys. Can anyone tell me how to determine points between two coardinates. For example: Which type of command line gives me 50 points between (8, -5, 7) and (2, 6, 9) points Thanks (5 Replies)
Discussion started by: rpf
5 Replies

6. Shell Programming and Scripting

place cursor in specific coordinates

Hi, I have this problem on how to place the cursor in a text editor (for example: pico). I made this script that would attach comments to a script file then open the script file, I would like to know how to place the cursor in a specific place, for example at the end of the comments, ... (1 Reply)
Discussion started by: lechelle
1 Replies

7. Shell Programming and Scripting

Search for particular tag and arrange as coordinates

Hi I have a file whose sample contents are shown here, 1.2.3.4->2.4.2.4 a(10) b(20) c(30) 1.2.3.4->2.9.2.4 a(10) c(20) 2.3.4.3->3.6.3.2 b(40) d(50) c(20) 2.3.4.3->3.9.0.2 a(40) e(50) c(20) 1.2.3.4->3.4.2.4 a(10) c(30) 6.2.3.4->2.4.2.5 c(10) . . . . Here I need to search... (5 Replies)
Discussion started by: AKD
5 Replies

8. Shell Programming and Scripting

Calculating distance between two LAT long coordinates

hi, i have a pair of latitude and longitude and i want to calculate the distance between these two points. In vbscript i achieved in the following way...Now i want to implement this in unix shell scripting.... <% Dim lat1, lon1, lat2, lon2 const pi = 3.14159265358979323846 ... (8 Replies)
Discussion started by: aemunathan
8 Replies

9. Shell Programming and Scripting

Defining X and Y Coordinates Inside A Window

Hello, I am starting up an Xnest window and trying to place a program inside of it. I have the window inside of it now but it always spawns with the top left corner at (0, 0). I need to find a way to set the x and y coordinates to something other than (0, 0). I tried using the -geometry option... (1 Reply)
Discussion started by: lesnaubr
1 Replies

10. Shell Programming and Scripting

Removing dupicate lines in the file ..(they are not continuous)

I have duplicates records in a file, but they are not consecutive. I want to remove the duplicates , using script. Can some one help me in writing a ksh script to implement this task. Ex file is like below. 1234 5689 4556 1234 4444 (7 Replies)
Discussion started by: Srini75
7 Replies
Login or Register to Ask a Question

LEARN ABOUT DEBIAN

utm2ll

UTM2LL(1)						      General Commands Manual							 UTM2LL(1)

NAME

       utm2ll - Convert UTM coordinates to latitude/longitude geographical coordinates

SYNOPSIS

       utm2ll [-L] | [utm_x utm_y zone [nad27 | nad83 | wgs84]]

DESCRIPTION

       This  program  uses Redfearn's formulas to convert a given set of Universal Transverse Mercator (UTM) coordinates into the equivalent lati-
       tude and longitude geographical coordinates.  (This operation is often referred to as inverse projection, since it projects  a  previously-
       projected  flat	surface  back onto the curved surface from whence it originally came.)	The inputs are the UTM "x" (also known as easting)
       value, the UTM "y" (also known as northing) value, and the utm zone.  The "x" value includes the normal 500,000	false  easting.   The  "y"
       value  includes	the  normal 10,000,000 false northing, if the point is in the southern hemisphere.  For points in the southern hemisphere,
       make the zone number negative.

       Warning:  Not all of the possible triples of utm_x, utm_y, and zone values represent correct UTM coordinates.  The program  will  generally
       still  produce  latitude/longitude  coordinates,  even for incorrect inputs.  You can check that your original inputs were correct by using
       ll2utm to convert the latitude/longitude coordinates back into UTM coordinates.

       The output takes the form of a single line, containing the latitude and longitude, separated by white space.  The  values  are  in  decimal
       degrees; with latitudes south of the equator being negative, and longitudes west of the prime meridian being negative.

       If you provide just the "-L" option, the program will print some license information and exit.

       Projections,  and  inverse  projections, depend on defining an ellipsoid that approximates the shape of the earth (the reference ellipsoid)
       and defining reference coordinates (the datum) that allow measurements to be made.  Different choices of the ellipsoid and datum can  yield
       projections  that  differ  by  tens  of meters.	There are a wide variety of choices, due to both the historical progression of measurement
       technology, and the desire to maximize accuracy over a given region (such as North America, or one of the United States).

       This program defaults to the North American Datum of 1927 (NAD-27) with the Clarke Ellipsoid of 1866, since these appear to be  appropriate
       for  much  of  the  freely-available data.  The data are apparently in the process of being converted to the Geodetic Reference System 1980
       (GRS-80) ellipsoid and NAD-83.  If you come across such data, you can specify "nad83" on the command line.  The GTOPO30 data use the  World
       Geodetic System 1984 (WGS-84) ellipsoid, which can be invoked by specifying "wgs84" on the command line.

SEE ALSO

       The ll2utm(1) command provides the inverse conversion.

								   Jul 24, 2001 							 UTM2LL(1)