%% 本代码实现功能是根据输入的shape文件(wgs84),对栅格数据进行统计均值。
%思路:首先根据待统计的栅格数据的信息(范围,分辨率等),然后对应把shape文件生成
%与栅格范围大小一致的mask文件,然后就可以矩阵之间运算了
% 作者:李艳忠 地点:美国奥本大学图书馆 2019-4-23 21:07
clc;clear;
inputRasterFile='data\19980101.tif';
%% 1.读取栅格文件
[DataRaster,Ref]=geotiffread(inputRasterFile); %读取全国范围的栅格数据
lonSpan=Ref.LongitudeLimits; %栅格经度范围 73 135
latSpan=Ref.LatitudeLimits; %栅格纬度范围 18 54
resolution=Ref.CellExtentInLatitude; %栅格分辨率0.2500
xgv=[lonSpan(1)+resolution/2:resolution:lonSpan(2)-resolution/2];
ygv=[latSpan(1)+resolution/2:resolution:la
Matlab读取shape文件并统计均值
最新推荐文章于 2025-03-18 22:35:42 发布