comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: MPFIT2DFUN- use for a time series of images?
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: MPFIT2DFUN- use for a time series of images? [message #75500 is a reply to message #75496] Fri, 11 March 2011 04:49 Go to previous messageGo to previous message
Allard de Wit is currently offline  Allard de Wit
Messages: 7
Registered: December 2002
Junior Member
On Mar 10, 9:01 pm, Katie <kmwhite21...@gmail.com> wrote:
> I want to fit pixel values to a logistic model that each have X, Y
> coordinates and a day of the year value associated with them.  I see
> from the MPFIT2DFUN documentation how to define X, Y for a each Z
> dependent value and the dependent error value, but can a third
> variable (day of the year in this case) be defined in order to fit a
> curve in MPFIT2DFUN?  I currently have each each date as a separate
> band in a stacked image file that I exported as an ASCII file.  I want
> to determine the day of the year that the pixel values (vegetation
> indices) reach a value on the the curve.
>
> Any suggestions, or perhaps I am not understanding the MPFIT2DFUN
> correctly (I am very new to IDL)?

Katie,

You will need to iterate over the rows and columns of your image and
fit
the logistic model to the time-series at each row/column value.
Indeed, MPFITFUN is one of the tools to do that.

A couple of years ago I implemented this approach for use with the
GIMMS
dataset. In fact it uses a double logistic model to model the NDVI
time-series. My implementation does not include a variable day-of-year
(just regular time-steps) but this is a trivial extension. It also
uses IDLs native Curve_Fit routine rather then MPFITFUN, although the
latter is better.

If you're interested I can send you a copy of that code.

Allard
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: How to check the variables by myself ?
Next Topic: Problem with cgcmdwindow__define

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Thu Oct 09 22:31:36 PDT 2025

Total time taken to generate the page: 0.48264 seconds