comp.lang.idl-pvwave archive: archive » Convolution with non-constant Kernel?

Home » Public Forums » archive » Convolution with non-constant Kernel?

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: Convolution with non-constant Kernel? [message #73573 is a reply to message #73468]

Fri, 19 November 2010 02:07

rogass
Messages: 200
Registered: April 2008

Senior Member

Hi, try the following routines for 2D images and 2D mask or mask sized
[maskx,masky,n_elements(image)],whereas boundary pixels are not
considered, large images and small masks will cause heavy memory
usage:

function cr_discrete_convolve,mask,im
szm = ulong(size(mask,/dimensions))
szi = ulong(size(im,/dimensions))
return, reform(/
over,total(total(im[cr_window(szi[0],szi[1],szm[0],szm[1])]* $
n_elements(szm) eq 2? rebin(/
sample,mask,szm[0],szm[1],szi[0]*szi[1]):$
mask,2),1),szi[0],szi[1])
end

function cr_window,sx,sy,wx,wy,norot=norot
wx = ulong(wx)
wy = ulong(wy)
sx = ulong(sx)
sy = ulong(sy)
base = ((ll=ulindgen(wx,wy) mod wx)) + transpose(ll)*sx
return, rebin(/sample,(keyword_set(norot)? temporary(base) : $
rot(reform(/over,temporary(base),wx,wy),180)),wx,wy,sx*sy) + $
rebin(/sample,ulindgen(1,1,sx*sy),wx,wy,sx*sy)
end

Cheers

CR

Report message to a moderator

[Message index]

		Convolution with non-constant Kernel? By: SonicKenking on Thu, 11 November 2010 16:56
		Re: Convolution with non-constant Kernel? By: Juggernaut on Sat, 13 November 2010 08:23
		Re: Convolution with non-constant Kernel? By: rogass on Fri, 12 November 2010 14:03
		Re: Convolution with non-constant Kernel? By: rogass on Fri, 19 November 2010 02:07
		Re: Convolution with non-constant Kernel? By: Gray on Wed, 17 November 2010 12:52

Previous Topic:	plot of implicit function
Next Topic:	Re: IDLPong, a Pong game in IDL.

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Tue Dec 02 08:06:28 PST 2025

Total time taken to generate the page: 1.12621 seconds