大家好，又见面了，我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。
Jetbrains全家桶1年46，售后保障稳定
clc
clear
path1='D:\...\GF1_PMS2_E116.7_N40.2_20130501_L1A0000152355-MSS2.tiff';
image=imread(path1);
R=image(:,:,1);
G=image(:,:,2);
B=image(:,:,3);
Nir=image(:,:,4);
[row,column,band]=size(image);
N=zeros(row,column);
for i=1:row %进行一一遍历
for j=1:column
N(i,j)=double((G(i,j)-Nir(i,j)))/double((G(i,j)+Nir(i,j)));
end
end
%N=double(G-Nir)/double(G+Nir);% 也可以采用这种语法，不过貌似有些问题，还没找到
imwrite(N,'ndwi.tif','compression','none','resolution',[1,1]);版权声明：本文内容由互联网用户自发贡献，该文观点仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请联系我们举报，一经查实，本站将立刻删除。
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/234390.html原文链接：https://javaforall.net
