matlab学习笔记
我叫杨傲天
做一个谦卑的绅士
展开
-
Grace deg1 deg2 GIA
链接:https://pan.baidu.com/s/1571FIkkWuR24Kd-ciBrmJw提取码:sodj原创 2021-08-10 09:48:38 · 418 阅读 · 7 评论 -
GRACE_Matlab_Toolbox工具箱使用说明
1.软件下载地址:https://github.com/fengweiigg/GRACE_Matlab_Toolbox2.由于目前使用的数据大多都是RL06产品,而测地所冯伟老师的工具箱处理的是RL05产品,所以需要进行一些修改:(1)gmt_replace_C20函数中RL05替换成RL06(2)gmt_replace_C21_S21_C22_S22函数中的get_days改为gmt_get_days, RL05替换成RL06(3)Degree 1和Degree 2 都选用最新的3.控制文..原创 2021-06-24 22:24:51 · 7059 阅读 · 64 评论 -
MATLAB画区域专题地图(省边界)
refvec=[1 60 70];worldmap([20 60],[70 140])%描绘出中国地图的轮廓,其中的China是描绘中国地图,要是换做其他地图,就可以换成‘World’,或者'North Pole', 'South Pole', 'Pacific'%陆地水格网数据显示geoshow(gca,data/24,refvec,'DisplayType','texturemap');%data是陆地水格网数据%海岸线边界显示S = shaperead('landareas','U...原创 2021-06-08 11:46:21 · 3186 阅读 · 2 评论 -
MATLAB根据经纬度格网数据画地图
worldmap([20 60],[70 140])%描绘出中国地图的轮廓%陆地水格网数据显示geoshow(gca,data,refvec,'DisplayType','texturemap');%data是陆地水格网数据%海岸线边界显示S = shaperead('landareas','UseGeoCoords',true);geoshow([S.Lat], [S.Lon],'Color','black');ylabel(colorbar,'陆地水/cm','fontsize',10)..原创 2021-06-05 16:23:23 · 5819 阅读 · 1 评论 -
读取GLDAS-1版本数据集Noah水文模型
clc;clear;% 批读取NC文件的准备工作datadir = 'D:\MATLAB\toolbox\fengweiigg-GRACE_Matlab_Toolbox-6e91b05\GRACE_data\gldas_noah\'; %指定批量数据所在的文件夹filelist = dir([datadir,'*.nc4']); %列出所有满足指定类型的文件% a = filelist(1).name; %查看要读取的文件的编号% b = fi...原创 2021-05-30 16:39:24 · 2166 阅读 · 7 评论 -
matlab 同时画多条折线
clc;clear;load('D:\MATLAB\toolbox\fengweiigg-GRACE_Matlab_Toolbox-6e91b05\GRACE_results\result-ys\grid_data.mat')water1=textread('C:\Users\admin\Desktop\西部.txt','%f');water2=textread('C:\Users\admin\Desktop\东北部.txt','%f');water3=textread('C:\Users\ad.原创 2021-05-29 20:24:39 · 1706 阅读 · 0 评论 -
某区域格网点陆地水储量变化之和
clc;%读取研究区域内网格点的经纬度[longitude,latitude]=textread('C:\Users\admin\Desktop\四个区域网格点经纬度.txt','%f%f');%纬度-90°~90°,经度-180°~180°,(90°,-180°)对应第一行,第一列。a、b为矩阵行列号a=longitude+181;b=91-latitude;%grid_data(纬度编号,经度编号,月份),对应内容为陆地水月变化for i=1:158 c=0; f...原创 2021-05-29 19:46:46 · 273 阅读 · 0 评论