|
Re: Pixelwise temporal trend, problem with REFORM [message #75291 is a reply to message #75285] |
Wed, 23 February 2011 04:54  |
Gray
Messages: 253 Registered: February 2010
|
Senior Member |
|
|
On Feb 23, 7:04 am, Aga <u...@compgroups.net/> wrote:
> Hi All, I am trying to perform multi-temporal pixelwise linear regression on stack images using LINFIT(X,Y). For that both X and Y should have the same dimension. My stack image = Y and has 3 Dims
>
> my_stack=FINDGEN(n_elements(all_files),xsize, ysize)
>
> Whereas X refers to the time and is defined as
>
> x=INDGEN(n_elements(all_files))
> X INT=Array[4]
> IDL> print, x
> 0 1 2 3
> I am trying to use REFORM to rescale X to the size of Y but unsuccessfully
> Error message: REFORM: New subscripts must not change the number elements in X.
> I also tried to do a combination of reform and rebin but rebin allows to multiply by the factor of the X Dim.
>
> Any suggestion would be much appreciated
> Thanks, Aga
you need rebin(x,n_elements(all_files),xsize,ysize,/sample)
|
|
|