How to define variable names for resampling?

Question

Devendra 2024-4-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2109496-how-to-define-variable-names-for-resampling

评论： Devendra 2024-4-20

I am using following code to resample the 10m green,red and nir bands to 20m resolution;

% Resample the 10m bands to 20m 
if ismember(files(i + j).name, {'*_B03_10m.jp2', '*_B04_10m.jp2', '*_B08_10m.jp2'})
    img = imresize(img, 0.5, 'Method', 'bilinear');
end

But this part of code is being not executed.

The names of band files are as follows

T43RGN_20210119T053141_B03_10m

T43RGN_20210119T053141_B04_10m

T43RGN_20210119T053141_B05_20m

T43RGN_20210119T053141_B06_20m

T43RGN_20210119T053141_B07_20m

T43RGN_20210119T053141_B08_10m

I request you to kindly suggest me how to execute above mentioned part of code also.

Deva

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Stephen23 2024-4-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2109496-how-to-define-variable-names-for-resampling#answer_1444936

编辑：Stephen23 2024-4-20

在 MATLAB Online 中打开

You should:

replace the ISMEMBER with ENDSWITH
get rid of those asterisks

For example:

endsWith(files(i+j).name, {'_B03_10m.jp2', '_B04_10m.jp2', '_B08_10m.jp2'})

or even simply:

endsWith(files(i+j).name, '_10m.jp2')

You might also like to use the IgnoreCase option.

NOTE: your code does not work because ISMEMBER matches text literally, so what you incorrectly think are wildcard characters that match anything are in fact literal asterisk characters that match... the asterisk character. The ISMEMBER docunentation does not mention wildcard or special characters anywhere: inventing features which do not exist is unlikely to work.

6 个评论
显示 4更早的评论隐藏 4更早的评论

Devendra 2024-4-20

编辑：Walter Roberson 2024-4-20

在 MATLAB Online 中打开

I am very sorry for my mistake. It worked well. Thanks a lot for your kind help. However, when I am writing this layer stacked data into envi file using enviwrite, it is giving following error. If you could please suggest me how to use it properly, it will be great help to me. the part of code is as follows;

% Stack the image
stacked_image(:,:,j+1) = img;
end
fn = {files(i + j).name};
disp(string(fn));
end
disp(size(stacked_image));
img_names = regexp(fn,'\d+','match','once');
img_names = cellstr(reshape(img_names,1,[]));
fname = string(img_names(i)) + '_' + '6bands' + '.dat'
disp(fname);
enviwrite(stacked_image,fname);

it is giving following information

'C:\working_bpt\bands\'

T43RGN_20210119T053141_B03_10m.jp2

T43RGN_20210119T053141_B04_10m.jp2

T43RGN_20210119T053141_B05_20m.jp2

T43RGN_20210119T053141_B06_20m.jp2

T43RGN_20210119T053141_B07_20m.jp2

T43RGN_20210119T053141_B08_10m.jp2

5490 5490 6

fname =

"43_6bands.dat"

43_6bands.dat

Incorrect number or types of inputs or outputs for function enviwrite.

Error in ds_layer_stack (line 35)

enviwrite(stacked_image,fname);

Looking forward for your kind suggestions to write it in envi format.

Devendra

Walter Roberson 2024-4-20

@Devendra

As discussed in https://www.mathworks.com/matlabcentral/answers/2109356-how-to-use-enviwrite-in-matlab enviwrite() expects a hypercube object. A version of enviwrite() that expects a cuboid of data is at https://github.com/akayasse/ENVI-read-and-write-for-MATLAB

Devendra 2024-4-20

Thanks for your help. I have written envi file using enviwrite function. However, map information (latitude and longitude) with projection are missing. Please suggest me how to incorporate map information into envi header file. I will appreciate your kind help.

Devendra

请先登录，再进行评论。

How to define variable names for resampling?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

6 个评论
显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

How to define variable names for resampling?

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

6 个评论 显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

6 个评论
显示 4更早的评论隐藏 4更早的评论