Azzera filtri
Azzera filtri

How do I define HANTS matlab parameters?

7 visualizzazioni (ultimi 30 giorni)
gauri
gauri il 1 Mag 2024
Risposto: John il 10 Mag 2024
I am using HANTS MATLAB code for the NDVI time series. I have defined the parameters as follows;
function [yOut, amp, phi]=ApplyHants(y,nb,nf,fet,dod,HiLo,low,high,delta)
nb = 22 ; total number of images
nf = 3 frequency above zero frequency
fet = 5
dod = 5
HiLo = 'none'
low = -1.0 ; minimum value of NDVI
high = 1.0 ; maximum value of NDVI
delta = 0.1
but it is giving output as zero values.
Please suggest me how to define these above mentioned parameters for my NDVI time series of one year time period having 22 images.
I would appreciate your kind help.
Gauri
  5 Commenti
gauri
gauri il 2 Mag 2024
Modificato: gauri il 3 Mag 2024
Thank you Aditya for your support. It did not work. The header file generated by applyHants is as follows
ENVI
description = {Exported from MATLAB}
samples = 22
lines = 1315
bands = 1153
data type = 4
interleave = bsq
So here bands and samples are interchanged. I am puuting a part of applyHants matlab which reads the parameters as follows
function [yOut, amp, phi]=ApplyHants(y,nb,nf,fet,dod,HiLo,low,high,delta)
if (max(size(size(y)))~=3)
error('Input data must be three dimensional [time,lat,lon]')
end
[ni,ny,nx]=size(y);
yOut= zeros(ni,ny,nx,'single');
amp = zeros(nf+1,ny,nx,'single');
phi = zeros(nf+1,ny,nx,'single');
ts=1:ni;
I layer stacked the 22 files in image as follows
image = (22, 1315, 1153)
and passed it to applyHants as follows
[yOut, amp, phi] = ApplyHants(image, nb, nf, fet, dod, HiLo, low, high, delta);
I request you to please suggest me how to make them compatable to each other so that code works properly.
Once again thank you so much for your valuable comments.
Gauri
gauri
gauri il 3 Mag 2024
Can I once again request you to please have a look on it and suggest me some solution?
Thanks a lot for your valuable comments.
Gauri

Accedi per commentare.

Risposte (2)

Manikanta Aditya
Manikanta Aditya il 3 Mag 2024
Modificato: Manikanta Aditya il 3 Mag 2024
Please consider these as my suggestations from my research and exploration and try to see if you can resolve the issue. Might not completely fix the issue.
From your description, it seems like the dimensions of your input data might not be matching the expected dimensions in the 'ApplyHants' function. The function expects the input data to be three-dimensional in the form [time, lat, lon], but your data seems to be in the form [bands, lines, samples].
You can try transposing your data to match the expected dimensions.
Here’s how you can do it:
% Transpose the dimensions
image = permute(image, [3 2 1]);
% Now pass the transposed image to the function
[yOut, amp, phi] = ApplyHants(image, nb, nf, fet, dod, HiLo, low, high, delta);
This will rearrange your data dimensions to match the expected [time, lat, lon] format. Please try this and let me know if it helps.
By using permute(image, [1, 3, 2]), you are rearranging the dimensions of the image array to [time, samples, lines], which should match the expected order [time, lat, lon] for the ApplyHants function.
After making this change, the ApplyHants function should work as expected, and you should get the desired output instead of zero values.
  8 Commenti
Manikanta Aditya
Manikanta Aditya il 5 Mag 2024
Modificato: Manikanta Aditya il 5 Mag 2024
The error you're encountering is related to the way MATLAB handles logical operations with arrays. In MATLAB, the && and || operators expect scalar inputs or arrays of the same size. However, the condition any(~ready) returns a single logical value (true or false), which is incompatible with the array ready in the expression any(~ready) && (nloop < nloopmax).
To fix this issue, you can replace the logical AND operation (&&) with the element-wise AND operation (&), and wrap the condition with the any function, like this:
while any(~ready & (nloop < nloopmax))
% ... rest of the loop code ...
end
This way, the element-wise AND operation (&) is performed between the negated ready array (~ready) and the condition nloop < nloopmax, and then the any function checks if any element in the resulting logical array is true.
Alternatively, you can restructure the condition using && and any separately, like this:
while (any(~ready)) && (nloop < nloopmax)
% ... rest of the loop code ...
end
In this case, the any(~ready) expression returns a single logical value, which can be used with the && operator along with the scalar condition nloop < nloopmax.
gauri
gauri il 5 Mag 2024
Modificato: gauri il 5 Mag 2024
Thank you very much. It has overcome that error. However, a new error has occured from the following code;
% Initialize matrices
mat = zeros(min(2*nf+1, ni), ni);
p = ones(size(y));
bool_out = (y < low) | (y > high);
p(bool_out) = 0;
outliers(bool_out) = 1;
nout = sum(p == 0);
Arrays have incompatible sizes for this operation.
Error in HANTS (line 80)
za = mat .* (p .* y);
I request you to please suggest me how to overcome this.
Gauri

Accedi per commentare.


John
John il 10 Mag 2024
The error "Arrays have incompatible sizes for this operation" is likely due to a mismatch in the dimensions of the arrays involved in the element-wise multiplication operation
mat .* (p .* y)
In the HANTS code, the line
mat = zeros(min(2*nf+1, ni), ni);
initializes a matrix mat with dimensions
(min(2*nf+1, ni), ni)
. However, the dimensions of the arrays p and y are not explicitly specified in the code snippet you provided. To resolve this issue, you need to ensure that the dimensions of mat , p , and y are compatible for element-wise operations.
Here are a few steps you can try: Check the dimensions of your input data y by using the size function:
[ni, ny, nx] = size(y);
Ensure that the dimensions of mat match the expected dimensions based on the input data y .
The first dimension of mat should be
min(2*nf+1, ni)
, and the second and third dimensions should match the dimensions of y . You can modify the initialization of mat as follows:
mat = zeros(min(2*nf+1, ni), ny, nx);
Initialize the array p with the same dimensions as y :
p = ones(ni, ny, nx);
Modify the line that creates the logical array bool_out to ensure that the dimensions match:
bool_out = (y < low) | (y > high);
Update the line that assigns values to p based on bool_out :
p(bool_out) = 0;
Update the line that calculates nout :
nout = sum(p == 0, 'all');
After making these changes, the dimensions of mat , p , and y should be compatible, and the element-wise multiplication operation
mat .* (p .* y)
should work correctly. If you still encounter issues, provide more information about the dimensions of your input data y and the specific error message you receive.

Prodotti


Release

R2024a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by