概述
clear,clc
% 读取变量
chla = ncread('A20021722002263.L3m_SNSU_CHL_chlor_a_4km.nc','chlor_a');% lon,lat,列*行,8640*4320
lat = ncread('A20021722002263.L3m_SNSU_CHL_chlor_a_4km.nc','lat');% 4320*1
lon = ncread('A20021722002263.L3m_SNSU_CHL_chlor_a_4km.nc','lon');% 8640*1
% lat_min = find(lat==7.9792);%最小是最大,坐标,50°N,1
% lat_max = find(lat==24.0208);%10,321
% lon_min = find(lon==104.9792);%最小是最小,坐标,105,121
% lon_max = find(lon==122.0208);%135,361
lon_position = find(lon>105 & lon<122);% 序号,整数,从小到大,408*1
lat_position = find(lat>8 & lat<24);% 序号,整数,从小到大,1*384
lon = lon(lon_position);% 经度,小数,从小到大,408*1
lat = lat(lat_position);% 纬度,小数,从大到小,1*384
% lat = flip(lat);% 纬度,小数,从小到大,384*1
chla = chla(lon_position,lat_position);
% lon_len = length(lon);
% lon_len = lon_len';
% lat_len = length(lat);
% lat_len = lat_len';
% chla = chla(lon_len,lat_len);
[X,Y] = meshgrid(lon,lat);% lat,lon 384*408
figure;
m_proj('miller','lon',[105 122],'lat',[8 24]);
m_coast('patch',[.7 .7 .7],'edgecolor','b');
m_grid('linest','none','box','fancy','tickdir','in');
xlabel('longitude');
ylabel('latitude');
title('Seasonal MODIS Chlorophyll Concentration of Southern China Sea 2002.2');
colorbar;
hold on
m_pcolor(X,Y,chla');
meshgrid处必须是经纬度作为lon, lat
绝对不能是他的位置lon_position,lat_p
否则就没有数据
最后
以上就是刻苦凉面为你收集整理的matlab画图一片空白的问题的全部内容,希望文章能够帮你解决matlab画图一片空白的问题所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
发表评论 取消回复