?[aR]=geotiffread（‘G:\tif\2001001.tif‘）;%先導入緯度數據
info=geotiffinfo（‘G:\tif\2001001.tif‘）;
e=dir（‘*.tif‘）;
years=17;data=zeros（346*607years）;
????for?i=1:years
????????bz=importdata（e（i）.name）;
????????bz=reshape（bz261*6311）;
????????data（:i）=bz;
????end
????qushi=zeros（346607）;pz=zeros（346607）;
for?i=1:length（data）
????????bz=data（i:）;
????if?max（bz）>0
????????bz=bz‘;
????????X=[ones（size（bz））?bz];
????????X（:2）=[1:years]‘;
????????[bbintrrintstats]?=?regress（bzX）;
????????pz（i）=stats（3）;
????????qushi（i）=b（2）;
????end
end
????sy2=find（pz>0.1）;
????qushi（sy2）=0;
????name1=‘置信區間p.tif‘;
????name2=‘通過0.1置信區間的斜率.tif‘;
????geotiffwrite（name2qushiR‘GeoKeyDirectoryTag‘info.GeoTIFFTags.GeoKeyDirectoryTag）;
????geotiffwrite（name1pzR‘GeoKeyDirectoryTag‘info.GeoTIFFTags.GeoKeyDirectoryTag）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????222??2018-03-12?17:36??2001-2018西北地區1000m的MODSI?NDVI\2001001.tfw

?????文件?????843421??2018-03-12?17:36??2001-2018西北地區1000m的MODSI?NDVI\2001001.tif

?????文件????????222??2018-03-12?17:36??2001-2018西北地區1000m的MODSI?NDVI\2001017.tfw

?????文件?????843421??2018-03-12?17:36??2001-2018西北地區1000m的MODSI?NDVI\2001017.tif

?????文件????????222??2018-03-12?17:41??2001-2018西北地區1000m的MODSI?NDVI\2001033.tfw

?????文件?????843421??2018-03-12?17:41??2001-2018西北地區1000m的MODSI?NDVI\2001033.tif

?????文件????????222??2018-03-12?17:41??2001-2018西北地區1000m的MODSI?NDVI\2001049.tfw

?????文件?????843421??2018-03-12?17:41??2001-2018西北地區1000m的MODSI?NDVI\2001049.tif

?????文件????????222??2018-03-12?17:40??2001-2018西北地區1000m的MODSI?NDVI\2001065.tfw

?????文件?????843421??2018-03-12?17:40??2001-2018西北地區1000m的MODSI?NDVI\2001065.tif

?????文件????????222??2018-03-12?17:40??2001-2018西北地區1000m的MODSI?NDVI\2001081.tfw

?????文件?????843421??2018-03-12?17:40??2001-2018西北地區1000m的MODSI?NDVI\2001081.tif

?????文件????????222??2018-03-12?17:39??2001-2018西北地區1000m的MODSI?NDVI\2001097.tfw

?????文件?????843421??2018-03-12?17:39??2001-2018西北地區1000m的MODSI?NDVI\2001097.tif

?????文件????????222??2018-03-12?18:58??2001-2018西北地區1000m的MODSI?NDVI\2001113.tfw

?????文件?????843421??2018-03-12?18:58??2001-2018西北地區1000m的MODSI?NDVI\2001113.tif

?????文件????????222??2018-03-12?18:57??2001-2018西北地區1000m的MODSI?NDVI\2001129.tfw

?????文件?????843421??2018-03-12?18:57??2001-2018西北地區1000m的MODSI?NDVI\2001129.tif

?????文件????????222??2018-03-12?18:57??2001-2018西北地區1000m的MODSI?NDVI\2001145.tfw

?????文件?????843421??2018-03-12?18:57??2001-2018西北地區1000m的MODSI?NDVI\2001145.tif

?????文件????????222??2018-03-12?18:56??2001-2018西北地區1000m的MODSI?NDVI\2001161.tfw

?????文件?????843421??2018-03-12?18:56??2001-2018西北地區1000m的MODSI?NDVI\2001161.tif

?????文件????????222??2018-03-12?18:55??2001-2018西北地區1000m的MODSI?NDVI\2001177.tfw

?????文件?????843421??2018-03-12?18:55??2001-2018西北地區1000m的MODSI?NDVI\2001177.tif

?????文件????????222??2018-03-12?19:14??2001-2018西北地區1000m的MODSI?NDVI\2001193.tfw

?????文件?????843421??2018-03-12?19:14??2001-2018西北地區1000m的MODSI?NDVI\2001193.tif

?????文件????????222??2018-03-12?19:13??2001-2018西北地區1000m的MODSI?NDVI\2001209.tfw

?????文件?????843421??2018-03-12?19:13??2001-2018西北地區1000m的MODSI?NDVI\2001209.tif

?????文件????????222??2018-03-12?19:13??2001-2018西北地區1000m的MODSI?NDVI\2001225.tfw

?????文件?????843421??2018-03-12?19:13??2001-2018西北地區1000m的MODSI?NDVI\2001225.tif

............此處省略764個文件信息