Everybody's heard about the blurred

Code: (Select All)

'OldMoses' first attempt at Gaussian Blur algorithm

SCREEN _NEWIMAGE(1024, 512, 32)

b& = _LOADIMAGE("bird sample.jpg", 32)

r& = _NEWIMAGE(_WIDTH(b&), _HEIGHT(b&), 32)



'seeding matrix  (7:1) x (1:7)

DIM a%(-3 TO 3)

a%(-3) = 1: a%(-2) = 2: a%(-1) = 4: a%(0) = 8: a%(1) = 4: a%(2) = 2: a%(3) = 1



'alpha multiplier matrix kernel (7:7)

DIM kern%(LBOUND(a%) TO UBOUND(a%), LBOUND(a%) TO UBOUND(a%))

FOR x% = LBOUND(kern%, 1) TO UBOUND(kern%, 1) '                 iterate kernel matrix columns

    FOR y% = LBOUND(kern%, 2) TO UBOUND(kern%, 2) '             iterate kernel matrix rows

        kern%(x%, y%) = a%(x%) * a%(y%) '                       seed the gaussian kernel position

NEXT y%, x%

mult! = 255 / kern%(0, 0) '                                     set a multiplier to prevent alpha bleed through



'display clear image, create and display blurred image

_PUTIMAGE (1, 1), b& '                                          place clear bird image on left of mainscreen

FOR ox% = 0 TO _WIDTH(b&) - 1 '

    FOR oy% = 0 TO _HEIGHT(b&) - 1 '

        _SOURCE b& '                                            source: clear bird image

        c~& = POINT(ox%, oy%) '                                 get source pixel color at (ox%, oy%)

        _DEST r& '                                              destination: blurred receiving image

        FOR x% = LBOUND(kern%, 1) TO UBOUND(kern%, 1) '         apply the alpha matrix around original point

            FOR y% = LBOUND(kern%, 2) TO UBOUND(kern%, 2)

                PSET (ox% + x%, oy% + y%), _RGBA32(_RED32(c~&), _GREEN32(c~&), _BLUE32(c~&), mult! * kern%(x%, y%))

        NEXT y%, x%

NEXT oy%, ox%

_DEST 0 '                                                       destination back to mainscreen

_PUTIMAGE (_WIDTH(0) - _WIDTH(r&) - 2, 1), r& '                 place blurred receiving image on right of mainscreen

END

' ______________________________________________________________________________________________________________________________________________ '/ \ FUNCTION ApplyGauss& (SourceHandle AS LONG) ' ApplyGauss& | ' __________________________________________________________________________________________________________________________________________|____ '/ \ '| Applies a Gaussian blur to the image passed in. | '| | '| BlurredImage = ApplyGauss&(OriginalImage) | '| | '| SourceHandle - the image to be blurred | '| | '| NOTE: This function is a modified version of RhoSigma's Image Processing Library's ApplyFilter& function found in imageprocess.bm. | '| The function has been modified to only support the "gauss8" method of blurring with no optional parameters available. | '| RhoSigma's unedited library can be obtained here: https://qb64phoenix.com/forum/showthread.php?tid=1033 | '| Thanks to RhoSigma for offering this library. | '\_______________________________________________________________________________________________________________________________________________/ DIM AS INTEGER SourceWidth, SourceHeight, FilterY, FilterX, y, x, NewAlpha, FilterWeight, NewRed, NewGreen, NewBlue, Size, Add, Div DIM AS LONG NewHandle, SumRed, SumGreen, SumBlue DIM AS _UNSIGNED LONG OriginalRGB, NewRGB DIM AS _OFFSET PixelOffset, FilterOffset DIM SourceBuffer AS _MEM DIM NewBuffer AS _MEM STATIC Weight(0 TO 6, 0 TO 6) AS INTEGER '+---------------------+ '| Apply filter values | '+---------------------+ Size = 3 ' set filter values for "gauss8" Add = 0 Div = 16 IF Weight(2, 2) = 0 THEN Weight(2, 2) = 1: Weight(2, 3) = 2: Weight(2, 4) = 1 Weight(3, 2) = 2: Weight(3, 3) = 4: Weight(3, 4) = 2 Weight(4, 2) = 1: Weight(4, 3) = 2: Weight(4, 4) = 1 END IF ApplyGauss& = -1 ' assume handle is invalid IF SourceHandle < -1 OR SourceHandle = 0 THEN ' valid source image handle? IF _PIXELSIZE(SourceHandle) = 4 THEN ' yes, 32bit image? '+-----------------------+ '| Copy the source image | '+-----------------------+ SourceWidth = _WIDTH(SourceHandle) ' yes, get source image dimensions SourceHeight = _HEIGHT(SourceHandle) NewHandle = _NEWIMAGE(SourceWidth, SourceHeight, 32) ' create new image canvas of same size _PUTIMAGE , SourceHandle, NewHandle ' copy source image '+----------------------+ '| Process copied image | '+----------------------+ $CHECKING:OFF IF NewHandle < -1 THEN ' is handle valid? Size = Size \ 2 ' yes, adjust filter size '+------------------------------------+ '| Use direct memory access for speed | '+------------------------------------+ SourceBuffer = _MEMIMAGE(SourceHandle) ' get image memory locations NewBuffer = _MEMIMAGE(NewHandle) '+-------------------------------------+ '| Iterate through source image pixels | '+-------------------------------------+ y = -1 ' set y location DO ' iterate vertically through source image y = y + 1 ' increment y location PixelOffset = (y * SourceWidth * 4) ' calculate pixel offset x = -1 ' set x location DO ' iterate horizontally through source image x = x + 1 ' increment x location _MEMGET SourceBuffer, SourceBuffer.OFFSET + PixelOffset, OriginalRGB ' get source image pixel NewAlpha = _ALPHA32(OriginalRGB) ' record pixel's alpha value SumRed = 0 ' clear previous summed pixel weight values SumGreen = 0 SumBlue = 0 '+-------------------------------------------------------+ '| Iterate through neigboring pixels using filter matrix | '+-------------------------------------------------------+ FilterY = y - Size - 1 ' calculate filter vertical start point DO ' iterate vertically through filter matrix FilterY = FilterY + 1 ' increment y location FilterOffset = (FilterY * SourceWidth * 4) + ((x - Size) * 4) ' calculate filter offset FilterX = x - Size - 1 ' calculate filter horizontal start point DO ' iterate horizontally through filter matrix FilterX = FilterX + 1 ' increment x location IF FilterY >= 0 AND FilterY < SourceHeight AND FilterX >= 0 AND FilterX < SourceWidth THEN ' is position outside image? _MEMGET SourceBuffer, SourceBuffer.OFFSET + FilterOffset, OriginalRGB ' no, get source image pixel ELSE _MEMGET SourceBuffer, SourceBuffer.OFFSET + PixelOffset, OriginalRGB ' yes, get center source image pixel END IF '+----------------------+ '| Sum up pixel weights | '+----------------------+ FilterWeight = Weight(FilterY - y + 3, FilterX - x + 3) ' get weight value from filter array SumRed = SumRed + (_RED32(OriginalRGB) * FilterWeight) ' apply weight value to RGB colors SumGreen = SumGreen + (_GREEN32(OriginalRGB) * FilterWeight) SumBlue = SumBlue + (_BLUE32(OriginalRGB) * FilterWeight) FilterOffset = FilterOffset + 4 ' increment to next filter offset LOOP UNTIL FilterX = x + Size LOOP UNTIL FilterY = y + Size NewRed = CINT(SumRed / Div) + Add ' calculate new pixel channel colors NewGreen = CINT(SumGreen / Div) + Add NewBlue = CINT(SumBlue / Div) + Add NewRGB = _RGBA32(NewRed, NewGreen, NewBlue, NewAlpha) ' calculate new pixel color _MEMPUT NewBuffer, NewBuffer.OFFSET + PixelOffset, NewRGB ' place new pixel color onto new image PixelOffset = PixelOffset + 4 ' increment to next pixel offset LOOP UNTIL x = SourceWidth - 1 LOOP UNTIL y = SourceHeight - 1 '+-----------------------------+ '| Free RAM then return result | '+-----------------------------+ _MEMFREE NewBuffer ' remove image buffers from RAM _MEMFREE SourceBuffer ApplyGauss& = NewHandle ' return new image END IF $CHECKING:ON END IF END IF END FUNCTION

Screen _NewImage(800, 600, 32) For t = 1 To 1000 size = Rnd * 255 x = Rnd * _Width: y = Rnd * _Height Line (x, y)-(x + size, y + size), _RGB(Rnd * 255, Rnd * 255, Rnd * 255), BF Next For t = 1 To 50000 PSet (Rnd * _Width, Rnd * _Height), _RGB(Rnd * 255, Rnd * 255, Rnd * 255) Next GaussianBlur 8, 100, 100, _Width - 100, _Height - 100 Sub GaussianBlur (radius, startx, starty, endx, endy) tempimage& = _CopyImage(_Display) origsource& = _Source origdest& = _Dest Dim pix As Long sigma = radius / 2 cons = 1 / (2 * 3.14159 * sigma * sigma) _Source tempimage& 'apply Gaussian Blur For y = starty To endy For x = startx To endx r = 0: g = 0: b = 0 For i = -radius To radius For j = -radius To radius If (y + j >= 0) And (x + i >= 0) And (y + j < _Height) And (x + i < _Width) Then pix = Point(x + j, y + i) k = cons * Exp(-((j * j + i * i) / (2 * sigma * sigma))) r = r + k * _Red32(pix) g = g + k * _Green32(pix) b = b + k * _Blue32(pix) End If Next Next PSet (x, y), _RGB32(r, g, b) Next Next _FreeImage tempimage& _Source origsource& _Dest origdest& End Sub