Differences

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