matlab两张栅格影像赋权重求和.txt

close all; 
clear all; 
clc
filepath1='F:\process\extract1.tif';

info = geotiffinfo(filepath1);
lat1=info.BoundingBox(1,2);
lat2=info.BoundingBox(2,2);
lon1=info.BoundingBox(1,1);
lon2=info.BoundingBox(2,1);
data1=imread(filepath1);
[m1,n1]=size(data1);
data1=data1(:);
data1(data1<0)=0;
%data1=double(data1);
%data1=(log(data1)-log(100))./(log(75000)-log(100));
data1(data1>1)=1;
data1(data1<0)=0;

filepath2='F:\process\extract2.tif';

data2=imread(filepath2);
[m2,n2]=size(data2);
data2=data2(:);
data2(data2<0)=0;
data_cal=data1.*0.5+data2.*0.5;
%data_cal=data1.*data2;
[X,ps]=mapminmax(data_cal',0,1);
%data_end=reshape(X,m,n);
data_end=reshape(X',m1,n1);