Show Posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.


Topics - _vince

Pages: 1 [2]
16
Programs / Floyd–Steinberg dithering
« on: July 16, 2018, 04:30:08 am »
implementation of https://en.wikipedia.org/wiki/Floyd%E2%80%93Steinberg_dithering

Code: QB64: [Select]
  1. deflng a-z
  2.  
  3. img1 = _loadimage("img/zmago.jpg", 32)
  4.  
  5. w = _width(img1)
  6. h = _height(img1)
  7.  
  8. img2 = _newimage(w, h, 32)
  9. img3 = _newimage(w, h, 32)
  10. img4 = _newimage(w, h, 32)
  11.  
  12. screen _newimage(w*2, h*2, 32)
  13.  
  14. dither_bw img1, img2, 0.1
  15.  
  16. dither img1, img3, 2
  17. dither img1, img4, 4
  18.  
  20. _putimage (0, 0), img1
  21. _putimage (w, 0), img2
  22. _putimage (0, h), img3
  23. _putimage (w, h), img4
  24.  
  28.  
  29. 'colour dither
  30. 'source image, destination image, number of colours per channel
  31. sub dither(img1, img2, num)
  32.         w = _width(img1)
  33.         h = _height(img1)
  34.  
  35.         _dest img2
  36.         _source img2
  37.  
  38.         _putimage , img1
  39.  
  40.         for y=0 to h-1
  41.         for x=0 to w-1
  42.  
  43.                 z = point(x, y)
  44.  
  45.                 r = (_red(z)*num\255)*255\num
  46.                 g = (_green(z)*num\255)*255\num
  47.                 b = (_blue(z)*num\255)*255\num
  48.  
  49.                 pset (x, y), _rgb(r, g, b)
  50.  
  51.                 qr = _red(z) - r
  52.                 qg = _green(z) - g
  53.                 qb = _blue(z) - b
  54.  
  55.                 r = _red(point(x + 1, y)) +   qr*7/16
  56.                 g = _green(point(x + 1, y)) + qg*7/16
  57.                 b = _blue(point(x + 1, y)) +  qb*7/16
  58.                 pset (x + 1, y), _rgb(r, g, b)
  59.  
  60.                 r = _red(point(x - 1, y + 1)) +   qr*3/16
  61.                 g = _green(point(x - 1, y + 1)) + qg*3/16
  62.                 b = _blue(point(x - 1, y + 1)) +  qb*3/16
  63.                 pset (x - 1, y + 1), _rgb(r, g, b)
  64.  
  65.                 r = _red(point(x, y + 1)) +   qr*5/16
  66.                 g = _green(point(x, y + 1)) + qg*5/16
  67.                 b = _blue(point(x, y + 1)) +  qb*5/16
  68.                 pset (x, y + 1), _rgb(r, g, b)
  69.  
  70.                 r = _red(point(x + 1, y + 1)) +   qr*1/16
  71.                 g = _green(point(x + 1, y + 1)) + qg*1/16
  72.                 b = _blue(point(x + 1, y + 1)) +  qb*1/16
  73.                 pset (x + 1, y + 1), _rgb(r, g, b)
  74.         next
  75.         next
  77.  
  78. 'black and white dither
  79. 'source image, destination image, bw threshold percent
  80. sub dither_bw(img1, img2, t as double)
  81.         w = _width(img1)
  82.         h = _height(img1)
  83.  
  84.         _dest img2
  85.         _source img2
  86.  
  87.         _putimage , img1
  88.  
  89.         for y=0 to h-1
  90.         for x=0 to w-1
  91.  
  92.                 z = point(x, y)
  93.  
  94.                 c = -((_red(z)+_green(z)+_blue(z))/3 > 255*t)*255
  95.  
  96.                 pset (x, y), _rgb(c, c, c)
  97.  
  98.                 qr = _red(z) - c
  99.                 qg = _green(z) - c
  100.                 qb = _blue(z) - c
  101.  
  102.                 r = _red(point(x + 1, y)) +   qr*7/16
  103.                 g = _green(point(x + 1, y)) + qg*7/16
  104.                 b = _blue(point(x + 1, y)) +  qb*7/16
  105.                 pset (x + 1, y), _rgb(r, g, b)
  106.  
  107.                 r = _red(point(x - 1, y + 1)) +   qr*3/16
  108.                 g = _green(point(x - 1, y + 1)) + qg*3/16
  109.                 b = _blue(point(x - 1, y + 1)) +  qb*3/16
  110.                 pset (x - 1, y + 1), _rgb(r, g, b)
  111.  
  112.                 r = _red(point(x, y + 1)) +   qr*5/16
  113.                 g = _green(point(x, y + 1)) + qg*5/16
  114.                 b = _blue(point(x, y + 1)) +  qb*5/16
  115.                 pset (x, y + 1), _rgb(r, g, b)
  116.  
  117.                 r = _red(point(x + 1, y + 1)) +   qr*1/16
  118.                 g = _green(point(x + 1, y + 1)) + qg*1/16
  119.                 b = _blue(point(x + 1, y + 1)) +  qb*1/16
  120.                 pset (x + 1, y + 1), _rgb(r, g, b)
  121.         next
  122.         next
  124.  

17
Programs / Cosine Transform
« on: June 27, 2018, 07:03:36 pm »
Neat program I found:

Code: QB64: [Select]
  1. defint a-z
  2.  
  3. dim shared mx,my,mbl,mbr,mw
  5. pi = 3.141593
  6.  
  7. sw = 800
  8. sh = 600
  9.  
  10. dim xx(sw - 1), x(sw - 1)
  11.  
  12. screen _newimage(sw,sh,32)
  13.  
  15.         getmouse
  16.  
  17.         line(0, 0)-(sw, sh), _rgb(0,0,0), bf
  18.  
  19.         line(sw/2, 0)-(sw/2, sh), _rgb(255,0,0)
  20.         line(0, sh/4)-(sw, sh/4), _rgb(255,0,0)
  21.         line(0, 3*sh/4)-(sw, 3*sh/4), _rgb(255,0,0)
  22.  
  23.         ox = 0
  24.         oy = 3*sh/4 - x(0)
  25.         for i = 0 to ubound(xx)
  26.                 if xx(i) <> 0 then line(i, sh/4)-step(0,-xx(i))
  27.  
  28.                 line(ox, oy)-(i, 3*sh/4 - x(i))
  29.                 ox = i
  30.                 oy = 3*sh/4 - x(i)
  31.         next
  32.  
  33.         'if m = 0 then
  34.                 line(mx, sh/4)-(mx, my)
  35.         'end if
  36.  
  37.         if mbl then
  38.                 do
  39.                         getmouse
  40.                 loop until mbl = 0
  41.  
  42.                 xx(mx) =  sh/4 - my
  43.  
  44.                 for i = 0 to ubound(x)
  45.                         x(i) = x(i) + xx(mx) * cos(2*pi*i*(mx - sw/2)*(sw))
  46.                 next
  47.         end if
  48.  
  49.         _display
  52.  
  53.  
  54. sub getmouse ()
  55.         do
  56.                 mx = _mousex
  57.                 my = _mousey
  58.                 mbl = _mousebutton(1)
  59.                 mbr = _mousebutton(2)
  60.                 mw = mw + _mousewheel
  61.         loop while _mouseinput
  63.  
  64.  

18
Programs / Lemon explorer
« on: March 19, 2018, 07:01:44 am »
Implementation of http://paulbourke.net/fractals/lemon/

Code: QB64: [Select]
  1. defint a-z
  2.  
  3. const sw = 640
  4. const sh = 480
  5.  
  6. dim shared mx,my,mbl,mbr,mw
  7.  
  10. dim uuu as double, vvv as double
  14. z = 0.01
  15. zz = 0.1
  16.  
  18. p1 = _newimage(sw, sh, 32)
  19. screen _newimage(sw, sh, 32)
  20.  
  21. redraw = -1
  22. iter = 100
  23.  
  25.         mw = 0
  26.         getMouse
  27.  
  28.         if redraw then
  29.                 line(0,0)-(sw,sh),_rgb(0,0,0),bf
  30.                 for y = 0 to sh-1
  31.                 for x = 0 to sw-1
  32.  
  33.                         xx = (x - sw/2)*z + x0
  34.                         yy = (y - sh/2)*z + y0
  35.                         u = yy
  36.                         v = xx
  37.                         j=0
  38.                         for i = 0 to iter
  39.                                 cmul uuu,vvv,u,v,u,v
  40.                                 cmul uu,vv,uuu,vvv,yy,xx
  41.                                 cmul uu,vv,uu,vv,uuu+1,vvv
  42.                                 uuu=uuu-1
  43.                                 cmul uuu,vvv,uuu,vvv,uuu,vvv
  44.                                 cdiv uu,vv,uu,vv,uuu,vvv
  45.  
  46.                                 if uu*uu + vv*vv > 10^10 then
  47.                                         i = iter
  48.                                         exit for
  49.                                 end if
  50.                                 if abs(u-uu)<0.00001 and abs(v-vv)<0.00001 then
  51.                                         j = j + 1
  52.                                 else
  53.                                         j = 0
  54.                                 end if
  55.                                 if j > 4 then exit for
  56.  
  57.                                 u = uu
  58.                                 v = vv
  59.                         next
  60.                         if i < iter then
  61.                                 c = 255-10*i
  62.                                 pset(x, y), _rgb(c,c,c)
  63.                         end if
  64.                 next
  65.                 next
  66.  
  67.                 'locate 1,1
  68.                 'print "iter =";iter
  69.  
  70.                 _dest p1
  71.                 _putimage , 0
  72.                 _dest 0
  73.  
  74.                 _putimage , p1
  75.                 _autodisplay
  76.  
  77.                 redraw = 0
  78.         end if
  79.  
  80.         if mw < 0 then
  81.                 zz = zz + 0.01
  82.         elseif mw > 0 then
  83.                 if zz > 0.01 then zz = zz - 0.01
  84.         end if
  85.  
  86.         'draw box
  87.         if omx <> mx or omy <> my or mw <> 0 then
  88.                 _putimage , p1
  89.                 line (mx - (sw*zz/2), my - (sh*zz/2))-step(sw*zz,sh*zz),_rgb(255,255,255),b
  90.                 _autodisplay
  91.  
  92.                 omx = mx
  93.                 omy = my
  94.         end if
  95.  
  96.         if mbl then
  97.                 do
  98.                         getMouse
  99.                 loop while mbl
  100.  
  101.                 x0 = x0 + (mx - sw/2)*z
  102.                 y0 = y0 - (sh/2 - my)*z
  103.                 z = z*zz
  104.                 redraw = -1
  105.         elseif mbr then
  106.                 do
  107.                         getMouse
  108.                 loop while mbr
  109.  
  110.                 x0 = x0 + (mx - sw/2)*z
  111.                 y0 = y0 - (sh/2 - my)*z
  112.                 z = z/zz
  113.                 redraw = -1
  114.         end if
  115.  
  116.         k = _keyhit
  117.         if k = 43 then
  118.                 iter = iter + 50
  119.                 redraw = -1
  120.         elseif k = 45 then
  121.                 if iter > 50 then iter = iter - 50
  122.                 redraw = -1
  123.         end if
  124.  
  125. loop until k = 27
  127.  
  128. sub getMouse ()
  129.         do
  130.                 mx = _mousex
  131.                 my = _mousey
  132.                 mbl = _mousebutton(1)
  133.                 mbr = _mousebutton(2)
  134.                 mw = mw + _mousewheel
  135.         loop while _mouseinput
  137.  
  138. sub cmul (u as double, v as double, x as double, y as double, a as double, b as double)
  139.         uu# = x*a - y*b
  140.         vv# = x*b + y*a
  141.         u = uu#
  142.         v = vv#
  144.  
  145. sub cdiv (u as double, v as double, x as double, y as double, a as double, b as double)
  146.         d# = a*a + b*b
  147.         uu# = (x*a + y*b)/d#
  148.         vv# = (y*a - x*b)/d#
  149.         u = uu#
  150.         v = vv#
  152.  

19
Programs / Chaotic scattering
« on: February 15, 2018, 06:52:14 pm »
https://en.wikipedia.org/wiki/Chaotic_scattering

Demo of the Gaspard-Rice system as described above.  Left-click to change location.

Code: QB64: [Select]
  1. defint a-z
  2. sw = 640
  3. sh = 480
  8.  
  12.  
  13. pi = 3.141593
  14.  
  15. screen _newimage(sw,sh,12)
  16.  
  17. r = 150
  18. rr = 100
  19.  
  20. 'line(0,0)-(sw,sh),0,bf
  21. 'for b=0 to 2*pi step 2*pi/3
  22. '       circle (r*cos(b)+sw/2, r*sin(b)+sh/2),rr
  23. 'next
  24. 'end
  25.  
  26. xx = sw/2
  27. yy = sh/2
  28.  
  30.         do
  31.                 mx=_mousex
  32.                 my=_mousey
  33.                 mb=_mousebutton(1)
  34.         loop while _mouseinput
  35.  
  36.         line(0,0)-(sw,sh),0,bf
  37.         for b=0 to 2*pi step 2*pi/3
  38.                 circle (r*cos(b)+sw/2, r*sin(b)+sh/2),rr
  39.         next
  40.  
  41.  
  42.         if mb then
  43.                 f = -1
  44.                 do while mb
  45.                         do
  46.                                 mb=_mousebutton(1)
  47.                         loop while _mouseinput
  48.                 loop
  49.                 for b=0 to 2*pi step 2*pi/3
  50.                         x1=r*cos(b)+sw/2
  51.                         y1=r*sin(b)+sh/2
  52.                         if (mx-x1)^2+(my-y1)^2 < rr*rr then f = 0
  53.                 next
  54.                 if f then
  55.                         xx = mx
  56.                         yy = my
  57.                         f = -1
  58.                 end if
  59.         end if
  60.  
  61.         x0 = xx
  62.         y0 = yy
  63.  
  64.         a = _atan2(my-yy,mx-xx)
  65.  
  66.         t=0
  67.         do
  68.                 t=t+1
  69.                 x = t*cos(a)+x0
  70.                 y = t*sin(a)+y0
  71.                 if x<0 or x>sw or y<0 or y>sh then exit do
  72.                 for b=0 to 2*pi step 2*pi/3
  73.                         x1=r*cos(b)+sw/2
  74.                         y1=r*sin(b)+sh/2
  75.                         if (x-x1)^2+(y-y1)^2 < rr*rr then
  76.                                 a1 = _atan2(y-y1,x-x1)
  77.                                 a2 = 2*a1-a-pi
  78.  
  79.                                 line(x0, y0)-(x,y),14
  80.  
  81.                                 x0 = x
  82.                                 y0 = y
  83.                                 a = a2
  84.                                 t=0
  85.                                 exit for
  86.                         end if
  87.                 next
  88.         loop
  89.  
  90.         line (x0,y0)-(x,y),14
  91.  
  92.         _display
  93.         _limit 50
  96.  

20
QB64 Discussion / Burning ship explorer
« on: July 21, 2017, 07:41:07 pm »
left mouse, zoom in
right mouse, zoom out
mousewheel, zoom size
+, increase iterations
-, decrease iterations

(edit: added other variations commented out)
Code: QB64: [Select]
  1. defint a-z
  2.  
  3. const sw = 640
  4. const sh = 480
  5.  
  6. dim shared mx,my,mbl,mbr,mw
  7.  
  13. z = 0.01
  14. zz = 0.1
  15.  
  17. p1 = _newimage(sw, sh, 12)
  18. screen _newimage(sw, sh, 12)
  19.  
  20. redraw = -1
  21. iter = 100
  22.  
  24.         mw = 0
  25.         getMouse
  26.  
  27.         if redraw then
  28.                 for y = 0 to sh-1
  29.                 for x = 0 to sw-1
  30.                         u = 0
  31.                         v = 0
  32.  
  33.                         xx = (x - sw/2)*z + x0
  34.                         yy = (y - sh/2)*z + y0
  35.                         for i = 0 to iter
  36.                                 '''mandelbrot
  37.                                 'u = u
  38.                                 'v = v
  39.                                 'uu = u*u - v*v + xx
  40.                                 'vv = 2*u*v + yy
  41.                                 '''
  42.  
  43.                                 '''burning ship
  44.                                 u = abs(u)
  45.                                 v = abs(v)
  46.                                 uu = u*u - v*v + xx
  47.                                 vv = 2*u*v + yy
  48.                                 '''
  49.  
  50.                                 '''tricorn
  51.                                 'u = u
  52.                                 'v = -v
  53.                                 'uu = u*u - v*v + xx
  54.                                 'vv = 2*u*v + yy
  55.                                 '''
  56.  
  57.                                 '''tetration
  58.                                 'u = u
  59.                                 'v = v
  60.                                 'cexp uu, vv, u, v, u, v
  61.                                 'cexp uu, vv, uu, vv, xx, yy
  62.                                 '''
  63.  
  64.                                 '''other, the possibilities are infinite
  65.                                 'if u > v then
  66.                                 '       u = -abs(u)
  67.                                 '       v = v
  68.                                 'elseif -1 then
  69.                                 '       u = u
  70.                                 '       v = -abs(v)
  71.                                 'end if
  72.                                 'uu = u*u - v*v + xx
  73.                                 'vv = 2*u*v + yy
  74.                                 '''
  75.  
  76.                                 u = uu
  77.                                 v = vv
  78.  
  79.                                 if (u*u + v*v) > 4 then exit for
  80.                         next
  81.                         if i > iter then i = 0
  82.                         pset(x, y), i
  83.                 next
  84.                 next
  85.  
  86.                 locate 1,1
  87.                 print "iter =";iter
  88.  
  89.                 _dest p1
  90.                 _putimage , 0
  91.                 _dest 0
  92.  
  93.                 _putimage , p1
  94.                 _display
  95.  
  96.                 redraw = 0
  97.         end if
  98.  
  99.         if mw < 0 then
  100.                 zz = zz + 0.01
  101.         elseif mw > 0 then
  102.                 if zz > 0.01 then zz = zz - 0.01
  103.         end if
  104.  
  105.         'draw box
  106.         if omx <> mx or omy <> my or mw <> 0 then
  107.                 _putimage , p1
  108.                 line (mx - (sw*zz/2), my - (sh*zz/2))-step(sw*zz,sh*zz),14,b
  109.                 _display
  110.  
  111.                 omx = mx
  112.                 omy = my
  113.         end if
  114.  
  115.         if mbl then
  116.                 do
  117.                         getMouse
  118.                 loop while mbl
  119.  
  120.                 x0 = x0 + (mx - sw/2)*z
  121.                 y0 = y0 - (sh/2 - my)*z
  122.                 z = z*zz
  123.                 redraw = -1
  124.         elseif mbr then
  125.                 do
  126.                         getMouse
  127.                 loop while mbr
  128.  
  129.                 x0 = x0 + (mx - sw/2)*z
  130.                 y0 = y0 - (sh/2 - my)*z
  131.                 z = z/zz
  132.                 redraw = -1
  133.         end if
  134.  
  135.         k = _keyhit
  136.         if k = 43 then
  137.                 iter = iter + 50
  138.                 redraw = -1
  139.         elseif k = 45 then
  140.                 if iter > 50 then iter = iter - 50
  141.                 redraw = -1
  142.         end if
  143.  
  144. loop until k = 27
  146.  
  147. sub getMouse ()
  148.         do
  149.                 mx = _mousex
  150.                 my = _mousey
  151.                 mbl = _mousebutton(1)
  152.                 mbr = _mousebutton(2)
  153.                 mw = mw + _mousewheel
  154.         loop while _mouseinput
  156.  
  157. 'u + iv = (x + iy) ^ (a + ib)
  158. sub cexp (u as double, v as double, x as double, y as double, a as double, b as double)
  159.         dim mag as double, arg as double
  160.         dim lnz as double, argz as double
  161.  
  162.         lnz = 0.5*log((x*x + y*y)+0.00001)
  163.         argz = _atan2(y, x)
  164.  
  165.         mag = exp(a*lnz - b*argz)
  166.         arg = a*argz + b*lnz
  167.  
  168.         u = mag * cos(arg)
  169.         v = mag * sin(arg)
  171.  

21
QB64 Discussion / Vince's Windowing System
« on: July 20, 2017, 02:11:30 am »
draggable scrollable windows, twm style resizing, and linked list powered textboxes in each window

Pages: 1 [2]