Convert pixels in lat-lon to kilometer square

2 次查看(过去 30 天)
How can I change the pixels (in lat longs) of some monthly .nc files to kilometere square for every month. Spatial resolution of the data is 0.5 x 0.625.
  2 个评论
Joydeb Saha
Joydeb Saha 2022-10-26
ncFiles = dir('E:\COT Work\New\COTdata\cot2\*.nc');
N = length(ncFiles) ;
Q = zeros(N,1) ;
for i = 1:N
ncFile = fullfile(ncFiles(i).folder,ncFiles(i).name);
data_Wwind=netcdf(ncFile);
cot=data_Wwind{'M2TMNXRAD_5_12_4_TAUHGH'}(:);
th3 = cot > 3;
B = double(th3);
Q(i)=sum(B(:) == 1);
end
I writing the code here. This counts the COT values which are above 3. Then it is counts the total number of pixels which are above 3 for for every data file. Basically, I want the cot values which are above pixel 3, and convert those pixel (sum up) to area (in km square). Therefore Q shall give the monthly area values of COT which are above threshold 3. May be a little modification in the present code is needed.

请先登录,再进行评论。

回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Agriculture 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by