jiaxinhuang
diff --git a/‎pycl/pycl30.py
Lines changed: 55 additions & 0 deletions b/‎pycl/pycl30.py
Lines changed: 55 additions & 0 deletions
diff --git a/‎pycl/pycl31.py
Lines changed: 51 additions & 0 deletions b/‎pycl/pycl31.py
Lines changed: 51 additions & 0 deletions
@@ -0,0 +1,55 @@
+#-*- coding:utf-8 -*-
+'''
+Created on 2017年5月9日
+
+@author: huangjiaxin
+'''
+import csv
+from PIL import Image
+import numpy as np
+
+def test1():
+    data = np.genfromtxt('yankeedoodle.csv', delimiter=',')
+    data = data[:, :-1]
+    data = data.reshape((-1,))
+    data = data[:-1]
+    data = data.reshape((139, 53)).transpose().reshape((-1,))
+    
+    im = Image.new('L', (139, 53))
+    im.putdata(data*255)
+    im.show()
+    
+    data = data.reshape((53, 139)).transpose().reshape((-1,))
+    data_str = ["%.5f" % i for i in data]
+    msg = []
+    for i in range(0, 7365, 3):
+        msg.append(chr(int(data_str[i][5] + data_str[i+1][5] + data_str[i+2][6])))
+    print ''.join(msg)
+
+
+def test2():
+    tempList = []
+    tempSList = []
+    with open('yankeedoodle.csv', 'rb') as f:
+        file = csv.reader(f)
+        for row in file:
+            if row[-1] == '':
+                tempList += [float(x) for x in row[:-1]]
+                tempSList += row[:-1]
+            else:
+                tempList += [float(x) for x in row[:]]
+                tempSList += row[:]
+    tempArray = np.array(tempList).reshape(139, 53)
+    im = Image.new('L', (53, 139))
+    im.putdata(tempArray)
+    im.show()
+    
+    print len(tempSList)
+    msg = []
+    for i in range(0, 7365, 3):
+        msg.append(chr(int(tempSList[i][-2] + tempSList[i+1][-2] + tempSList[i+2][-1])))
+    print ''.join(msg)
+    
+
+if __name__ == '__main__':
+    test2()
@@ -0,0 +1,51 @@
+#-*- coding:utf-8 -*-
+'''
+Created on 2017年5月10日
+
+@author: huangjiaxin
+'''
+from PIL import Image
+import numpy as np
+
+def mapIdx(c, maxIter=128):
+    z = 0
+    for k in range(maxIter):
+        z = z * z + c
+        if np.abs(z) > 2:
+            break
+    return k
+
+def solve():
+    left = 0.34
+    bottom = 0.57
+    width = 0.036
+    height = 0.027
+    
+    im = Image.open('mandelbrot.gif')
+    imdata = np.array(list(im.getdata()))
+    im_w, im_h = im.size
+    dw = width / im_w
+    dh = height / im_h
+    xx = np.linspace(left, left + width - dw, im_w)
+    yy = np.linspace(bottom, bottom + height - dh, im_h)
+    yy = yy[::-1]
+    xx.shape = (1, im_w)
+    yy.shape = (im_h, 1)
+    grids = xx + 1j * yy
+    
+    for i in range(im_h):
+        for j in range(im_w):
+            grids[i, j] = mapIdx(grids[i, j])
+    
+    im.putdata(grids.reshape((-1,)))
+#     im.show()
+    
+    imdata.shape = (480, 640)
+    diffs = imdata - grids
+    msg = diffs[np.where(np.abs(diffs) == 16)]
+    newIm = Image.new('1', [23, 73])
+    newIm.putdata([i<0 for i in msg])
+    newIm.save('level31_out.png')
+    
+if __name__ == "__main__":
+    solve()