Matlab读取shape文件并统计均值

这段Matlab代码演示了如何读取shape文件和栅格数据,通过创建mask文件来统计shape文件区域内栅格数据的均值。首先,它读取栅格文件并获取其范围和分辨率。然后,它读取shape文件并生成mask文件。最后,使用mask文件对栅格数据进行处理并计算区域内数据的平均值。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

%% 本代码实现功能是根据输入的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

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值