How to extract satellite measured chlorophyll-a (chl-a) in MATLAB

9 次查看（过去 30 天）
Karthik M
Karthik M 2021-7-28
评论： Karthik M 2021-8-3
Hi Folks,
I was using thie following script to extract chlor_a variable with lat[0S 30N] and lon[30W 100E] further need to capture the data values of chla for particular region esp for North Indian Ocean lat, long range
ncdisp = ('A20150012015008.L3m_8D_CHL_chlor_a_4km.nc');
lat = ncread('A20150012015008.L3m_8D_CHL_chlor_a_4km.nc','lat');
lon = ncread('A20150012015008.L3m_8D_CHL_chlor_a_4km.nc','lon');
chlor_a = ncread('A20150012015008.L3m_8D_CHL_chlor_a_4km.nc','chlor_a');
% Load longitude array:
lon = double(ncread('A20150012015008.L3m_8D_CHL_chlor_a_4km.nc','lon'));
ind1 = find(lon>=30 & lon<=100);
% Do the same for lat:
lat = double(ncread('A20150012015008.L3m_8D_CHL_chlor_a_4km.nc','lat'));
ind2 = find(lat>=0 & lat<=30);
% Clip lat and lon to their specified range:
lat = lat(ind2);
lon = lon(ind1);
% Make a grid:
[Lat,Lon] = meshgrid(lat,lon);
Looking forward to your valuable suggestions
regards

请先登录，再回答此问题。

采纳的回答

KSSV
KSSV 2021-7-28
Let lon,lat be your grid coordinates of the nc file.
[X,Y] = meshgrid(lon,lat) ;
xi = linspace(30,100) ;
yi = linspace(0,30) ;
[Xi,Yi] = meshgrid(xi,yi) ;
iwant = interp2(X,Y,chlor_a',Xi,Yi) ;
If the data is 3D run a loop.
  12 个评论
显示 10更早的评论
Karthik M
Karthik M 2021-8-2
Hi @KSSV
I tried using this script still error persist
nccreate('myfile.nc','myvar')
>> A = 99;
>> ncwrite('myfile.nc','myvar',A)
>> S = ncinfo('myfile.nc');
>> file_fmt = S.Format
file_fmt =
'netcdf4_classic'
>> S.Format = 'netcdf4';
>> ncwriteschema('newfile.nc',S)
kindly help
Karthik M
Karthik M 2021-8-3
Hi @KSSV
I tried using this script still error persist
nccreate('myfile.nc','myvar')
>> A = 99;
>> ncwrite('myfile.nc','myvar',A)
>> S = ncinfo('myfile.nc');
>> file_fmt = S.Format
file_fmt =
'netcdf4_classic'
>> S.Format = 'netcdf4';
>> ncwriteschema('newfile.nc',S)
kindly help

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 IEC Certification Kit (for ISO 26262 and IEC 61508) 的更多信息

标签

产品


版本

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by