Recent Changes - Search:

Top

Marine gravitymeter, S-99

Land gravitymeter, G-936

Data format, conversion

Air/Sea operation

Base readings

Data processing

Projects

Other information

edit SideBar

September 22, 2022, at 07:44 AM
Version: pmwiki-2.3.22

SOFTWARE FOR FORMAT CONVERSION


Input file


Input file has this format:

Time,Position: Filtered vessel position: East,North,Lat,Long,Data: Depth {input 3: ch1},Data: Gravimeter {input 2: ch1},Data: Magnetometer {input 1: ch1},
19:59:18 09.08.2008,  448842.7430,8119926.7155,  073°10.095816',  013°25.020571',     0.0000,  10080.4100,  53616.4760
19:59:28 09.08.2008,  448822.5048,8119909.8114,  073°10.086439',  013°24.983845',  1328.1900,  10080.3300,  53616.8390
19:59:38 09.08.2008,  448802.1602,8119892.8544,  073°10.077033',  013°24.946924',  1328.1500,  10080.2500,  53616.3550

Output file


Output file has this format:

09.08.2008 19:59:18 073 10.095816' 013 25.020571' 0.0000 10080.4100 53616.4760
09.08.2008 19:59:28 073 10.086439' 013 24.983845' 1328.1900 10080.3300 53616.8390
09.08.2008 19:59:38 073 10.077033' 013 24.946924' 1328.1500 10080.2500 53616.3550

Software


This Python script does the re-formatting job:

Attach:preprocess-gravitydata-27Oct2009.py

#!/usr/bin/python
#
'''
************************************************************************* 
PRE-PROCESSING OF GRAVITY DATA


Program version:
  ------------------------------------------
  Ver. Date          By   Description
  ------------------------------------------
  0.1  27 Oct 2009  O.M. Initial version  


                      Dept. of Earth Science
                       University of Bergen
                              Norway
************************************************************************* 
'''


from string import split, replace
from sys import argv, exit

if len(argv) <> 2:
	print "Need name of input file as argument"
	exit(-1)

FileName = argv[1]

#---------------------------------------------------
#---  S t a r t   o f   p r o g r a m
#---------------------------------------------------

f = open(FileName)							# Open input file.

for line in f:								# Loop over all lines in input file.


	if line[0].isdigit():					# First get rid of lines with text or comments
		#print line,
		line = line.replace(',','')			# Get rid of all ','
		s = line.split()
		date = s[1]
		time = s[0]
		lat = s[3]
		lon = s[4]
		depth = s[5]
		gravity = s[6]
		magnetometer = s[7]
		#print '------------'
		#print line,
		#---- We need to split lat/long on the degree symbol, which has ASCII code 0xb0
		print date, time, \
		      lat.split(chr(0xb0))[0], lat.split(chr(0xb0))[1], \
		      lon.split(chr(0xb0))[0], lon.split(chr(0xb0))[1], \
			  depth, gravity, magnetometer

f.close()
Edit - History - Print - Search
Page last modified on October 27, 2009, at 12:39 PM
Electronics workshop
Department of Earth Science - University of Bergen
N O R W A Y