Graphic with line erros
Show older comments
Im running my code and Im trying to make the x axis of the figure go from -180 to 180 i change the code and it aligns but at the same time it has mistakes, being those mistakes line throught the the entire y axis
I think the problem is the fact that the lon values (the x axis is longitude) go from 0.5 to 359.5 and the lat (y axis) go from -89.5 to 89.5 all in increments of 1 (from 0.5 to 1.5 and -89.5 to -88.5, respectively)
can somebody help me
10 Comments
Vasco
on 4 Sep 2023
dpb
on 4 Sep 2023
You can't plot what you don't have--as descirbed, you have a gap at the 0-degree location. You would have to fill that in to fill in the gap in the image.
However, there are a couple of solid bands that don't seem to fit -- what are they?
Attach the data would likely let somebody diagnose; otherwise, it's pretty tough to supply explicit answers.
Vasco
on 4 Sep 2023
Star Strider
on 4 Sep 2023
The code would help as well, since (I hope) it demonstrates how to create the plot.
Vasco
on 4 Sep 2023
Vasco
on 4 Sep 2023
Star Strider
on 4 Sep 2023
I cannot run the poste c ode with the provided file.
A better option would be to use the zip function to zip the original files and then upload the .zip file here.
Vasco
on 4 Sep 2023
I was hoping for the actual .nc files in the .zip file.
F1 = fileread('GRD-3_2002094-...0_LND_v04.txt')
NC1 = ncread('GRD-3_2002094-...0_LND_v04.txt')
The ncread function throws errors with this one.
I do not know what the problem is, since I do not usually work with .nc files.
.
Vasco
on 4 Sep 2023
Accepted Answer
% unzip the zip-file
unzip('NC_file.zip')
close all;
clc;
format long g;
% input = 'C:\Users\vasco\OneDrive\Ambiente de Trabalho\ESTAGIO_Vasco\JPL TELLUS GRACE Level-3\';
% output = 'C:\Users\vasco\OneDrive\Ambiente de Trabalho\ESTAGIO_Vasco\outputs_teste\';
% I don't have those directories, so I use the current directory ('.') for
% input and output:
input = '.';
output = '.';
lista_ficheiros = dir(fullfile(input, '*.nc'));
lon_1 = double(ncread(fullfile(input, lista_ficheiros(1).name), 'lon'));
lat_1 = double(ncread(fullfile(input, lista_ficheiros(1).name), 'lat'));
nlon = length(lon_1);
nlat = length(lat_1);
ntime = length(lista_ficheiros);
lwe_data = zeros(nlon, nlat, ntime);
for i = 1:length(lista_ficheiros)
filename = fullfile(input, lista_ficheiros(i).name);
nc = netcdf.open(filename);
lwe = ncread(filename, 'lwe_thickness');
lwe_data(:, :, i) = lwe;
netcdf.close(nc);
lat=ncread(filename, 'lat');
lon=ncread(filename, 'lon');
% lon is [0->360]. make it [0->180, -180->0]
idx = lon>180;
lon(idx) = lon(idx) - 360;
% circshift lon so that the first element that was >180 becomes the
% first element in lon, making lon [-180->180]
shift_amount = 1-find(idx,1);
lon = circshift(lon,shift_amount);
% circshift lwe the same amount along its 1st dimension (lon dimension)
lwe = circshift(lwe,shift_amount,1);
time = double(ncread(filename, 'time'));
tm = string(time);
nome_ficheiro = strcat('GRACE-', tm);
figure(1);
axis equal;
clev = -2.0:0.1:0.5;
[x,y] = meshgrid(lon,lat);
contourf(x, y, lwe', clev, 'LineStyle', 'none', 'Fill', 'on');
clim([min(clev), max(clev)]);
colormap(winter(length(clev)-1));
colorbar('eastoutside');
title(strcat('GRACE-', tm));
%estatisticas
nval = nnz(~isnan(lwe));
s = nansum(lwe);
media = nansum(s) / nval;
%Guardar figuras
figura = fullfile(output, strcat(nome_ficheiro, '.png'));
saveas(gcf, figura);
estatisticas = fullfile(output, strcat(nome_ficheiro, '_estatisticas.txt'));
fid = fopen(estatisticas, 'w');
fprintf(fid, 'Estatisticas GRACE-%s:\n', tm);
fprintf(fid, 'Soma lwe_thickness: %f\n', s);
fprintf(fid, 'Media lwe_thickness: %f\n', media);
fclose(fid);
% don't close the figure (so you can see it)
% close(gcf);
end
lwe_data_file = fullfile(output, 'lwe_data.mat');
save(lwe_data_file, 'lwe_data')
lwe_transpor = permute(lwe_data, [3, 1, 2]);
pixel_media = mean(lwe_transpor, 1);
pixel_outputs = fullfile(output, 'pixel_media.mat');
save(pixel_outputs, 'pixel_media');
nome_final='Media total';
figure(2);
clevv = -3.0:0.1:0.5;
data = squeeze(pixel_media(1, :, :));
rotated = imrotate(data, 270);
corrigido = fliplr(rotated);
centered_longitude = linspace(-180, 180, 360);
final = circshift(corrigido, [0,180]);
contourf(centered_longitude, lat, final, clevv, 'LineStyle', 'none', 'Fill', 'on');
clim([min(clevv), max(clevv)]);
colorbar('eastoutside');
colormap(winter(length(clevv) - 1));
hAx = gca;
hAx.YDir = 'normal';
figura_1 = fullfile(output, strcat(nome_final, '.png'));
saveas(gcf, figura_1);
6 Comments
Vasco
on 5 Sep 2023
Voss
on 5 Sep 2023
You're welcome!
@Vasco: See below for how to limit the lat/lon range:
% unzip the zip-file
unzip('NC_file.zip')
close all;
clc;
format long g;
% input = 'C:\Users\vasco\OneDrive\Ambiente de Trabalho\ESTAGIO_Vasco\JPL TELLUS GRACE Level-3\';
% output = 'C:\Users\vasco\OneDrive\Ambiente de Trabalho\ESTAGIO_Vasco\outputs_teste\';
% I don't have those directories, so I use the current directory ('.') for
% input and output:
input = '.';
output = '.';
lista_ficheiros = dir(fullfile(input, '*.nc'));
lon_1 = double(ncread(fullfile(input, lista_ficheiros(1).name), 'lon'));
lat_1 = double(ncread(fullfile(input, lista_ficheiros(1).name), 'lat'));
nlon = length(lon_1);
nlat = length(lat_1);
ntime = length(lista_ficheiros);
lwe_data = zeros(nlon, nlat, ntime);
for i = 1:length(lista_ficheiros)
filename = fullfile(input, lista_ficheiros(i).name);
nc = netcdf.open(filename);
lwe = ncread(filename, 'lwe_thickness');
lwe_data(:, :, i) = lwe;
netcdf.close(nc);
lat=ncread(filename, 'lat');
lon=ncread(filename, 'lon');
% lon is [0->360]. make it [0->180, -180->0]
idx = lon>180;
lon(idx) = lon(idx) - 360;
% circshift lon so that the first element that was >180 becomes the
% first element in lon, making lon [-180->180]
shift_amount = 1-find(idx,1);
lon = circshift(lon,shift_amount);
% circshift lwe the same amount along its 1st dimension (lon dimension)
lwe = circshift(lwe,shift_amount,1);
% logical index of which lon values are within the range you want
lon_idx = lon>=-100 & lon<=-30;
% logical index of which lat values are within the range you want
lat_idx = lat>=-66 & lat<=20;
% keep only those values from lon, lat, lwe:
lon = lon(lon_idx);
lat = lat(lat_idx);
lwe = lwe(lon_idx,lat_idx);
time = double(ncread(filename, 'time'));
tm = string(time);
nome_ficheiro = strcat('GRACE-', tm);
figure(1);
axis equal;
clev = -2.0:0.1:0.5;
[x,y] = meshgrid(lon,lat);
contourf(x, y, lwe', clev, 'LineStyle', 'none', 'Fill', 'on');
clim([min(clev), max(clev)]);
colormap(winter(length(clev)-1));
colorbar('eastoutside');
title(strcat('GRACE-', tm));
%estatisticas
nval = nnz(~isnan(lwe));
s = nansum(lwe);
media = nansum(s) / nval;
%Guardar figuras
figura = fullfile(output, strcat(nome_ficheiro, '.png'));
saveas(gcf, figura);
estatisticas = fullfile(output, strcat(nome_ficheiro, '_estatisticas.txt'));
fid = fopen(estatisticas, 'w');
fprintf(fid, 'Estatisticas GRACE-%s:\n', tm);
fprintf(fid, 'Soma lwe_thickness: %f\n', s);
fprintf(fid, 'Media lwe_thickness: %f\n', media);
fclose(fid);
% don't close the figure (so you can see it)
% close(gcf);
end
lwe_data_file = fullfile(output, 'lwe_data.mat');
save(lwe_data_file, 'lwe_data')
lwe_transpor = permute(lwe_data, [3, 1, 2]);
pixel_media = mean(lwe_transpor, 1);
pixel_outputs = fullfile(output, 'pixel_media.mat');
save(pixel_outputs, 'pixel_media');
nome_final='Media total';
figure(2);
clevv = -3.0:0.1:0.5;
data = squeeze(pixel_media(1, :, :));
% rotated = imrotate(data, 270);
% corrigido = fliplr(rotated);
% centered_longitude = linspace(-180, 180, 360);
% final = circshift(corrigido, [0,180]);
final = circshift(data,shift_amount,1);
final = final(lon_idx,lat_idx);
contourf(lon, lat, final.', clevv, 'LineStyle', 'none', 'Fill', 'on');
clim([min(clevv), max(clevv)]);
colorbar('eastoutside');
colormap(winter(length(clevv) - 1));
hAx = gca;
hAx.YDir = 'normal';
figura_1 = fullfile(output, strcat(nome_final, '.png'));
saveas(gcf, figura_1);
Vasco
on 6 Sep 2023
Voss
on 6 Sep 2023
You're welcome!
Vasco
on 6 Sep 2023
More Answers (0)
Categories
Find more on Green in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
