二蛤莫挨老子-CSDN博客

原创风生近惯性振荡，用surface drifter数据滤波，得到近惯性流速 matlab

用surface drifter数据滤波，得到近惯性流速

2022-10-07 21:36:38 498 2

原创基于Landsat8数据利用ENVI5.1反演海表温度SST——以杭州湾为例

1.下载数据（下载网址：http://ids.ceode.ac.cn/query.html）（1）设置参数（2）选择高清无云的图像下载2. 辐射校正，ENVI5.1(1)Toolbox—Radiometric Correction—Radiometric Calibration，选择_MTL 文件（2）设置参数，导出文件3.大气校正（1）导入辐射校正所得的图像，在弹出的对话框选择第二个选项（2）设置参数，Flight time可以从“*_MTL.txt” 中找到，Ground El

2021-03-03 23:52:06 5036 9

原创 matlab 准平稳时间序列分析预报方法

clear all;close all;clc%% 读取数据lat=ncread('NCEP.sst.mean.nc','lat');lon=ncread('NCEP.sst.mean.nc','lon');sst=ncread('NCEP.sst.mean.nc','sst');time=ncread('NCEP.sst.mean.nc','time');sst(sst == -9.969209968386869e+36) = NaN;%% 插值[Lat Lon]=meshgri.

2020-06-11 15:05:13 819

原创 matlab 降水量海表温度的滑动平均

滑动平均：降水量比如 1951 1952 1953 1954 1955 1956 1957 1958 1959 1960 1961 1962 1963 1964 1965选取k=5则序列一是求 1951 1952 1953 1954 1955的平均降水量，平均降水量对应着第三个年份，即1953年。同理 1952 1953 1954 1955 1956的平均降水量对应着第三个年份，即1954年下面是滑动平均（moving smooth）函数%% moving smoothfunction

2020-05-20 19:25:47 2077 4

原创 matlab 风向玫瑰图

clear all;close all;clcwind=xlsread('路径名／文件名.xlsx');dir=wind(:,3);x=0:22.5:360;dir=dir*pi/180;x=x*pi/180;figureset(gcf,'color','w')rose(dir,x);view(90,-90);

2020-05-15 18:26:52 6467 2

原创 matlab滑动t检验和全球海表平均温度突变检验

滑动t检验程序一：clear all;close all;clc%% 读取数据data = xlsread('Tair.xlsx'); % 读取数据time = data(:,1);%时间序列数据temperature = data(:,2);%最高气温数据%% 画时间序列figure(1)set(gcf,'color','w')plot(time,temperature)xlabel('year')ylabel('temperature(\circC)')title('19

2020-05-15 09:39:20 3305 1

原创 Matlab 用corrcoef函数计算蒸发和降水的相关性，风场散度和降水的相关性

clear all;close all;clc%%lon = ncread('hawaii_evap.nc','longitude');lat = ncread('hawaii_evap.nc','latitude');evap= ncread('hawaii_evap.nc','param03');%蒸发prec= ncread('hawaii_evap.nc','param02');%降水corr = zeros(length(lon), length(lat) );%for i =

2020-05-09 20:31:31 911 1

原创 Matlab 画全球风场图并计算散度风和旋度风

clear all;close all;clc%%fname = 'wind_uv300.nc';lat = ncread(fname, 'lat');lon = ncread(fname, 'lon');uwind = ncread(fname, 'U');vwind = ncread(fname, 'V');uwind = uwind(:,:,2);vwind = vwind(:,:,2);%[lat2d, lon2d] = meshgrid(lat,lon);% plo

2020-05-09 20:30:18 6721 3

原创 matlab polyfit计算盐度的线性趋势并画图

clear all;close all;clc%%lon = ncread('hawaii_soest_salt.nc','longitude');lat = ncread('hawaii_soest_salt.nc','latitude');time= ncread('hawaii_soest_salt.nc','time');sss_tmp1 = ncread('hawaii_soest_salt.nc','saltydsl');%设置起始位置和最终位置sss_tmp2 = squeeze

2020-05-09 20:29:03 1799

原创 matlab nc数据绘制位密平面图和剖面图

注：这是我老师给的程序，自己稍微改动了下clear all;close all;clc%%% change the directory according to your datayour_dir = '/Users/boopstore3/Desktop/matlab/shiyanchengxu/seven/';fname = [ your_dir 'hawaii_Hydrobas...

2020-05-06 20:49:56 3751 3

原创 matlab 批量导入CTD数据+温盐剖面图

clear all;close all;clc%% 读取数据filelist=dir('/Users/boopstore3/Desktop/matlab/shiyanchengxu/six/CTD_data'); filename1 = filelist(3:6,:);for i=1:4 filename=filename1(i).name; [lat,lon,gtime,...

2020-05-02 22:16:32 5700 6

原创 matlab GOCI水色数据画TSS变化过程+批量读取.he5文件

下载地址：http://kosc.kiost.ac.kr/eng/clear all;close all;clc%% 读取经纬度file='COMS_GOCI_L2P_GA_20110524031644.LON_1.he5';LON=h5read(file,'/HDFEOS/GRIDS/Image Data/Data Fields/Longitude...

2020-05-01 12:41:44 2676 3

原创 matlab AVISO高度计数据画sla+u动图

下载地址https://marine.copernicus.euclear all; close all; clc%% 读取数据file='northwest_pacific.nc';sla=ncread(file,'sla'); u=ncread(file,'ugos'); v=ncread(file,'vgos');lon=ncread(file,'longitude'); lat...

2020-05-01 12:18:57 4199 8

原创 matlab HYCOM数据画温度断面图

下载地址：http://www.hycom.org/dataserver/gofs-3pt1/analysis%%clear,clc%读取温度clear,clctemperature=[];datadir='路径名'; %指定批量数据所在的文件夹filelist=dir([datadir,'*.nc']); %指定批量数据的类型k=length(filelist);for s=...

2020-05-01 12:07:44 4946 2

原创 matlab HYCOM数据画盐度分布图和流场图

下载的网址http://www.hycom.org/dataserver/gofs-3pt1/analysis%% 画盐度分布图clear,clc%读取经纬度lat = ncread('文件名.nc','lat');lon = ncread('文件名.nc','lon');[Lat,Lon]=meshgrid(lat,lon);%读取盐度s1=ncread('文件名.nc'...

2020-05-01 11:55:14 5356 5

原创 Matlab 用MODIS数据画全球海表温度分布图

MODIS数据可以从下面的网址下载https://podaac.jpl.nasa.gov/datasetlist?ids=Collections&values=MODIS_L3_SST&view=listclear,clcfile='文件名.nc';sst=ncread(file,'sst');lon=ncread(file,'lon');lat=ncread(file...

2020-05-01 11:44:40 9160 8

原创 matlab 用m_map画地形水深图

clear,clc%投影m_proj('albers equal-area','lat',[5 43],'long',[105 125],'rect','on');%加载水深文件[CS,CH]=m_etopo2('contourf',[-7000:500:3000],'edgecolor','none');%绘制高分辨率陆地和海岸线m_gshhs_h('patch',[0.1 0.9 ...

2020-05-01 11:25:14 8398 5