Hi Randall Bonnell,
As per my the understanding of the question there are two problems in the question, the first is the error encountered and the second to convert the band into an imaginary component:
- The error "Unable to read band 2. Band values must be less than or equal to the total number of bands, 1", suggests that the GeoTIFF file may have only one band available and here more than one band is being accessed which is a possible reason for the error.
- To verify the number of bands in the geotiff file, you may use the "geotiffinfo" function from the Mapping Toolbox as suggested below.
% Get information about the geotiff file
info = geotiffinfo(filename);
- To create a complex array for the phase change, read the real and imaginary components separately and then combine them using the "complex" function as suggested below:
% Read in the geotiff first (real) band
[Real, R_r] = readgeoraster(pathName, "OutputType", "double", "Bands", 1);
% Read in the geotiff second (imaginary) band
[Imaginary, R_i] = readgeoraster(pathName, "OutputType", "double", "Bands", 2);
% Combine real and imaginary components into a complex array
Complex = complex(Real, Imaginary);
Please refer to the following MATLAB Documentation for better understanding of the functions:
geotiffinfo
complex function
I hope this helps!
Thank You
Maneet Bagga
