Построение выходных данных модели CMAQ из файла NetCDF - PullRequest
Новая
       15

Построение выходных данных модели CMAQ из файла NetCDF

0 голосов
/ 04 февраля 2019

Я впервые работаю с форматами файлов NetCDF, поэтому извините за незнание, если я не очень хорошо спрашиваю.

В настоящее время я работаю в MATLAB с данными NetCDF.Вот мой код для его просмотра: 

clear all;
close all;
folder_name = uigetdir;

% Display the file contents
filename = 'CCTM_US1k.CONC.1km_HGB.2013244';
ncdisp(filename)
pathname2 = uigetdir('/projectnb/atmchem/ahsouri/MCIP')
filename2 = fullfile(pathname2, 'GRIDBDY2D_1km_HGB_2013242')
%% Read the Info
info = ncinfo(filename);
info2 = ncinfo(filename2);

%% Display the variables names
nVar = size(info.Variables, 2);
disp('Variable names are:')
for ii = 1:nVar
    disp(info.Variables(ii).Name);
end

nVar2 = size(info2.Variables, 2);
disp('Variable names are:')
for ii = 1:nVar2
    disp(info2.Variables(ii).Name);
end

%% Read Variables
Data = ncread(filename2,'LAT');
Data2 = ncread(filename2, 'LON');
Data3 = ncread(filename, 'NO2');

%% Mapping Texas
figure 
usamap('Texas')
Texashi = shaperead('usastatehi','UseGeoCoords', true, ...
    'Selector', {@(name) strcmpi(name,'Texas'), 'Name'});
geoshow(Texashi, 'FaceColor', [0.3 1.0, 0.675])

textm(Texashi.LabelLat, Texashi.LabelLon, Texashi.Name, ...
    'HorizontalAlignment', 'center')
geoshow(Data,Data2)

Данные, которые я просматриваю, имеют следующие атрибуты: Их имена в файлах: 

Data = LAT (Obviously a latitude)
Data2 = LON (Longitude)
Data3 = NO2 (Obviously data regarding nitrogen dioxide)

NO2 происходит из этого файла: 

CCTM_US1k.CONC.1km_HGB.2013244

LAT и LON происходят из этого файла: 

GRIDBDY2D_1km_HGB_2013242

Ниже приведен код, который я использую для анализа вывода модели 

size(Data)
ans = 
        1204

whos Data
Name      Size      Bytes    Class    Attributes
Data    1204x1       4816    Single

size(Data2)
ans = 
        1204
whos(Data2)
Name    Size     Bytes     Class     Attributes
Data2    1204x1   4816     single

size(Data3)
ans =
    336 264 27 25
whos(Data3)
Name     Size    Bytes    Class    Attributes
Data3    4-D   239500800  single

Я просто хочу узнать, как можно использовать инструмент картографирования MATLAB для правильного отображения данных NO2 в сетке на карте Техаса, которую я создал.В настоящее время я в состоянии построить окно «Широта и долгота» просто отлично, я просто не знаю, как обращаться с данными NO2.

Любая помощь приветствуется, и я прошу прощения за невежество, я пытался предоставитькаждый бит информации в моем распоряжении.

Я также был бы более чем открыт, чтобы увидеть ваше решение аналогичной проблемы в альтернативном языке программирования!

Best, Тейлор

...