How to plot 2D density plot with very large files?
8 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Vikash Raj
le 19 Fév 2023
Réponse apportée : Vikash Raj
le 20 Fév 2023
Hi
I want to plot 2D density plots with a colourbar for very large files. I want my x ais to be longitude, y axis to be latitude and the z axis as absolute VTEC as colour bar. I was trying to run the codes given below but I getting a error. as
"Error using repmat
Requested 373471x373471 (1039.2GB) array exceeds maximum array size preference (15.9GB). This might cause MATLAB to become unresponsive.
Error in meshgrid (line 61)
xx = repmat(xrow,size(ycol));"
Any suggestions on how can I fix this so that I can obtain my dentsity contour plots ?
Im attaching my part of data. Is it possible for the graph to have blue colour as background colour as shown with the figure attached.
load('Data 1.mat')
whos
xi =Longitude;
yi = Latitude;
[numel(xi), numel(uniquetol(xi))]
[numel(yi), numel(uniquetol(yi))]
[X,Y] = meshgrid(xi,yi);
Z = Absolute_VTEC;
pcolor(X,Y,Z)
shading interp
colorbar
colormap(jet)
2 commentaires
dpb
le 19 Fév 2023
373471.^2*8/1024.^2
MB -- you don't have enough memory; decimate X,Y until you can hold the resulting array in memory. With only roughly 10,000 pixels on a monitor in each direction, you can't discriminate more than that many discrete locations, anyway.
Réponse acceptée
Walter Roberson
le 19 Fév 2023
Modifié(e) : Walter Roberson
le 19 Fév 2023
You do not have a grid of lat and long: you have scattered data. You can see from the ratio of the number of unique elements to the number of stored elements that it does not happen to be a grid of data that has been reshaped into a row.
A density plot reflects how dense points are, and takes as input vectors of x and vectors of y -- no third value. You cannot do a density plot of scattered data when you have three variables. You can do things like take means.
load('Data 1.mat')
scatter(Longitude, Latitude, 1, Absolute_VTEC)
mask = isfinite(Latitude) & isfinite(Longitude);
Latitude = Latitude(mask); Longitude = Longitude(mask);
N = 100;
[latbins, latedges] = discretize(Latitude, N);
[lonbins, lonedges] = discretize(Longitude, N+1);
denmat = accumarray([latbins(:), lonbins(:)], Absolute_VTEC(mask), [], @mean, NaN);
latmid = (latedges(1:end-1) + latedges(2:end))/2;
lonmid = (lonedges(1:end-1) + lonedges(2:end))/2;
h = pcolor(lonmid, latmid, denmat);
h.EdgeColor = 'none';
colorbar()
2 commentaires
Walter Roberson
le 20 Fév 2023
Yes. The code divides the range of lat and long values into intervals; the intervals applied in two different directions form a rectangular grid. The code figures out where on the grid that each value belongs, and for each grid point takes the mean of the vtec values that fall into that grid point. With this particular coding the grid locations with no elements are assigned nan, but that is easy to change.
Plus de réponses (1)
Voir également
Catégories
En savoir plus sur 2-D and 3-D Plots dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!