fast spectral resampling [message #71362] |
Fri, 11 June 2010 03:59 |
lbusett@yahoo.it
Messages: 30 Registered: February 2006
|
Member |
|
|
Hi all,
I need to resample of a huge number (hundreds of thousands) of
spectra produced by a radiative transfer model (2100 bands, with a 1
nm spectral interval) to the spectral characteristics of a given
sensor, expressed as the center wavelength and gaussian FWHM of the
different bands. The number of output bands can be as high as a few
thousands.
I'm currently building a matrix (hcrf_in) containing some thousands of
the input spectra (one spectra per column) and performing the
resampling with the ENVI_RESAMPLE_SPECTRA routine, which gives in
output a new matrix (hcrf_out) containing the resampled data. The
command is like the following.
ENVI_RESAMPLE_SPECTRA, in_wl, hcrf_in, out_wl, hcrf_out,
OUT_DT=float , OUT_FWHM=out_fwhm , interleave = 2
The routine works well, but it is really slow, and I don't have access
to the code in order to try to optimize the computation.
Does anybody know if there are any other routines or methods that can
be implemented in order to speed-up the computation ? Keep in mind
that that spectral sampling and FWHM of the different output bands can
be variable.
Thanks in advance for the help,
Lorenzo
|
|
|