This shows you the differences between two versions of the page.
plotting [2011/10/18 17:54] m |
plotting [2013/11/19 13:03] (current) m |
||
---|---|---|---|
Line 445: | Line 445: | ||
unset dgrid3d | unset dgrid3d | ||
splot "/root/projects/sommercamp/120801scry.log.clean" with lines | splot "/root/projects/sommercamp/120801scry.log.clean" with lines | ||
+ | </file> | ||
+ | |||
+ | ====== entropy plot/location scripts ====== | ||
+ | |||
+ | <file python gpsrng2.py> | ||
+ | from pylab import * | ||
+ | import csv, os, Gnuplot, Gnuplot.funcutils | ||
+ | |||
+ | g = Gnuplot.Gnuplot(debug=1) | ||
+ | |||
+ | NMI = 1852.0 | ||
+ | D2R = pi/180.0 | ||
+ | |||
+ | #data = sys.argv[1] | ||
+ | #title = sys.argv[2] | ||
+ | |||
+ | def read_csv_file(filename): | ||
+ | data = [] | ||
+ | for row in csv.reader(open(filename)): | ||
+ | data.append(row) | ||
+ | return data | ||
+ | |||
+ | def process_gps_data(data): | ||
+ | latitude = [] | ||
+ | longitude = [] | ||
+ | intensity = [] | ||
+ | lf = [] | ||
+ | gsr = [] | ||
+ | temp = [] | ||
+ | for row in data: | ||
+ | latitude.append(float(row[0][0:2]) + \ | ||
+ | float(row[0][2:])/60.0) | ||
+ | longitude.append((float(row[1][0:3]) + \ | ||
+ | float(row[1][3:])/60.0)) | ||
+ | intensity.append(float(row[2])) | ||
+ | # lf.append(float(row[3])) | ||
+ | |||
+ | return (array(latitude), array(longitude), \ | ||
+ | array(intensity)#, array(lf)) | ||
+ | |||
+ | y=read_csv_file('/root/collect2012-3/monkx') | ||
+ | (lat, long, intensity) = process_gps_data(y) | ||
+ | |||
+ | # what is lat, long for min/max intensity? | ||
+ | |||
+ | minty=min(intensity) | ||
+ | maxxy=max(intensity) | ||
+ | |||
+ | for index, item in enumerate(intensity): | ||
+ | if item==minty: | ||
+ | lower=index | ||
+ | if item==maxxy: | ||
+ | max=index | ||
+ | |||
+ | print "low entropy sites -: ", | ||
+ | print lat[lower], long[lower] | ||
+ | print "\nlow entropy sites +: ", | ||
+ | print lat[max], long[max] | ||
+ | </file> | ||
+ | |||
+ | <file python gpsrngcum.py> | ||
+ | |||
+ | |||
+ | from pylab import * | ||
+ | import csv, os, Gnuplot, Gnuplot.funcutils | ||
+ | |||
+ | g = Gnuplot.Gnuplot(debug=1) | ||
+ | |||
+ | NMI = 1852.0 | ||
+ | D2R = pi/180.0 | ||
+ | |||
+ | #data = sys.argv[1] | ||
+ | #title = sys.argv[2] | ||
+ | |||
+ | def calc_limit_high_005(range): | ||
+ | var = (1.96 / 2) * math.sqrt(range) | ||
+ | limit_high = range/2 + var | ||
+ | return limit_high | ||
+ | |||
+ | def read_csv_file(filename): | ||
+ | data = [] | ||
+ | for row in csv.reader(open(filename)): | ||
+ | data.append(row) | ||
+ | return data | ||
+ | |||
+ | def process_gps_data(data): | ||
+ | latitude = [] | ||
+ | longitude = [] | ||
+ | RNG = [] | ||
+ | ppp = [] | ||
+ | cnt=1 | ||
+ | rngcum=0 | ||
+ | for row in data: | ||
+ | latitude.append(float(row[0][0:2]) + \ | ||
+ | float(row[0][2:])/60.0) | ||
+ | longitude.append((float(row[1][0:3]) + \ | ||
+ | float(row[1][3:])/60.0)) | ||
+ | rngcum+=(100-float(row[2])) | ||
+ | RNG.append(rngcum) | ||
+ | ppp.append((calc_limit_high_005(cnt*200))-(cnt*100)) | ||
+ | cnt+=1 | ||
+ | |||
+ | return (array(latitude), array(longitude), \ | ||
+ | array(RNG),array(ppp)) | ||
+ | |||
+ | y=read_csv_file('/root/collect2011/psych/studies/symptoms_newcastle/logs/dayoneRNG') | ||
+ | (lat, long, RNG, ppp) = process_gps_data(y) | ||
+ | # translate spherical coordinates to Cartesian | ||
+ | py = (lat-min(lat))*NMI*60.0 | ||
+ | px = (long-min(long))*NMI*60.0*cos(D2R*lat) | ||
+ | #newintensity=intensity-min(intensity) | ||
+ | |||
+ | #pack px, py, intensity and gsr into newy | ||
+ | newy=[] | ||
+ | cummd=0 | ||
+ | for x,yz,zz,pppp in zip(px,py,RNG,ppp): | ||
+ | newy.append((x,yz,zz,pppp)) | ||
+ | |||
+ | g('set parametric') | ||
+ | g('set style data line') | ||
+ | g('set surface') | ||
+ | g('unset key') | ||
+ | g('unset contour') | ||
+ | #g('set dgrid3d 80,80,30') | ||
+ | g('set dgrid3d 80,80,30') | ||
+ | g('set xlabel "metres WE"') | ||
+ | g('set ylabel "metres NS"') | ||
+ | #g('set label "signal intensity" at -100,0,100') | ||
+ | g('set view 60,20') | ||
+ | g.title("9th September 2011 Newcastle symptoms cumulative RNG") | ||
+ | #g('set term png size 14043,9933') # A0 | ||
+ | g('set term png size 1024,768') # example | ||
+ | g('set output "/root/collect2011/psych/studies/symptoms_newcastle/logimages/newcumRNG.png"') | ||
+ | g('set style lines 1') | ||
+ | |||
+ | g.splot(Gnuplot.Data(newy, using=(1,2,3)),Gnuplot.Data(newy, using=(1,2,4))) | ||
+ | |||
+ | #g.splot(Gnuplot.Data(newy, using=(1,2,3), with='lines')) | ||
</file> | </file> | ||