Image functions: Convert Cartesian to Polar coordinates

[loudar]

I wrote this function as part of Vedit, the image editor made in QB64, but I wanted to share it because I can't recall seeing one like it here. I used the Wiki example for modifying image data in memory and adapted it to read the pixel data, then write them back into a new image with a bit of math for the new x and y coordinates.

Enjoy!

Code: QB64: [Select]

1. SUB RectangleToPolar (Image AS LONG)
2.     IF R < 0 OR R > 1 OR G < 0 OR G > 1 OR B < 0 OR B > 1 OR _PIXELSIZE(Image) <> 4 THEN EXIT SUB
3.     DIM Buffer AS _MEM: Buffer = _MEMIMAGE(Image) 'Get a memory reference to our image
4.  
5.     maxx = _WIDTH(Image)
6.     maxy = _HEIGHT(Image)
7.     mx = maxx / 2
8.     my = maxy / 2
9.  
10.     DIM O AS _OFFSET, O_Last AS _OFFSET
11.     O = Buffer.OFFSET 'We start at this offset
12.     O_Last = Buffer.OFFSET + maxx * maxy * 4 'We stop when we get to this offset
13.     DIM imgPoints(maxy, maxx, 0 TO 3) AS _UNSIGNED _BYTE
14.     DIM p, maxp AS _UNSIGNED _INTEGER64
15.     maxp = maxx * maxy
16.     'use on error free code ONLY!
17.     $CHECKING:OFF
18.     DO
19.         p = p + 1
20.         y = FIX((p / maxp) * maxy)
21.         x = p - (FIX((p / maxp) * maxy) * maxx)
22.         imgPoints(y, x, 0) = _MEMGET(Buffer, O, _UNSIGNED _BYTE)
23.         imgPoints(y, x, 1) = _MEMGET(Buffer, O + 1, _UNSIGNED _BYTE)
24.         imgPoints(y, x, 2) = _MEMGET(Buffer, O + 2, _UNSIGNED _BYTE)
25.         imgPoints(y, x, 3) = _MEMGET(Buffer, O + 3, _UNSIGNED _BYTE)
26.         O = O + 4
27.     LOOP UNTIL O = O_Last
28.     'create new image
29.     DIM newImg AS LONG
30.     newImg = _NEWIMAGE(maxx, maxy, 32)
31.     prevDest& = _DEST
32.     _DEST newImg
33.     COLOR _RGBA(0, 0, 0, 255), _RGBA(0, 0, 0, 0)
34.     CLS
35.     scaley = (maxx / 2 / _PI)
36.     IF scaley > maxy THEN scaley = maxy / 2
37.     y = -1: DO: y = y + 1
38.         yfactor = (1 - (y / (maxy))) * scaley
39.         x = -1: DO: x = x + 1
40.             PSET (mx - (yfactor * SIN(2 * _PI * (1 - (x / maxx)))), my - (yfactor * COS(2 * _PI * (1 - (x / maxx))))), _RGBA(imgPoints(y, x, 2), imgPoints(y, x, 1), imgPoints(y, x, 0), imgPoints(y, x, 3))
41.         LOOP UNTIL x = maxx
42.     LOOP UNTIL y = maxy
43.     _DEST Image
44.     COLOR _RGBA(0, 0, 0, 255), _RGBA(0, 0, 0, 0)
45.     CLS
46.     _PUTIMAGE (0, 0)-(maxx, maxy), newImg, Image
47.     _FREEIMAGE newImg
48.     ERASE imgPoints
49.     DIM imgPoints(0, 0, 0) AS _UNSIGNED _BYTE
50.     _DEST prevDest&
51.     'turn checking back on when done!
52.     $CHECKING:ON
53.     _MEMFREE Buffer
54. END SUB