含有经纬度tiff格式图像批量读取并画图.m

clear all
clc;
Info=imfinfo（‘悬浮物月分布图_2010010120100131_20150701.tiff‘）;
%  输入经纬度范围
user_lon_range = [107109.5];
user_lat_range = [1819.8];
xdata_resolution=0.016676251;
ydata_resolution=0.016687376;
%算出经纬度最小最大行列数
minlon=round（（min（user_lon_range）-98）/ydata_resolution）; 
maxlon=round（（max（user_lon_range）-98）/ydata_resolution）;
minlat=round（（25-min（user_lat_range））/xdata_resolution）;
maxlat=round（（25-max（user_lat_range））/xdata_resolution）;
row=minlat-maxlat;
col=maxlon-minlon;
%输入路径
data_dir = ‘D:\01-hlx\03-SSC\Matdata\dry\‘; 
files = dir（ [data_dir ‘悬浮物月分布图_*.tiff‘]）; 
data_cnt = length（files）;  
for i=1:data_cnt
    datafile_name = files（i1）.name;           %文件名
    dataflie_path = [data_dir datafile_name];  
    a=imread（dataflie_path）;
    data（::i）=a;
	SSC_da