How do i count the flashes of a lightning discharge within my time of interest and plot them in per square km ?? I want the plot to look like this(capture.jpg --attached below).
6 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
I have the data of lightning discharges around a location (say, South Africa). I have the latitude and longitude of the flashes occurring over this region. How do I plot them on map, to know their flash density ?? Anyone please help formulate a matlab code to know the no. of flashes per square kilometer.
0 Commenti
Risposta accettata
Chad Greene
il 17 Mar 2016
The first step is to convert lat,lon to some projected x,y values because lats and lons are not equally spaced. If you have the Mapping Toolbox you can use projfwd for the coordinate transformation.
When you have the coordinates in x,y space, then it's just a matter of creating a 2D histogram to get lightning strike density. You can search the File Exchange site for 2D histogram functions that others have already created. Or you could manually create a 2D histogram by creating a grid, then in a loop you could determine the sum of x,y values that are within the bounds of each grid cell. Here's an example:
x = 20*randn(300,1);
y = 10*randn(300,1)+x/3;
plot(x,y,'rx')
axis equal;
gridres = 10;
xgrid = -100:gridres:100;
ygrid = -50:gridres:50;
% Preallocate a histogram:
hist2 = NaN(length(ygrid),length(xgrid));
for row = 1:length(ygrid)-1
for col = 1:length(xgrid)-1
% Get indices of all x,y points inside this grid cell:
ind = x>=xgrid(col)-gridres/2 & x<xgrid(col+1)-gridres/2 & y>=ygrid(row)-gridres/2 & y<ygrid(row+1)-gridres/2;
%
% Log the sum of number of lighting strikes inside the grid cell:
hist2(row,col) = sum(ind);
end
end
%
hold on
h = imagesc(xgrid,ygrid,hist2/(gridres^2));
uistack(h,'bottom')
0 Commenti
Più risposte (1)
kaustubh ahirrao
il 28 Mar 2021
x = 20*randn(300,1);
y = 10*randn(300,1)+x/3;
plot(x,y,'rx')
axis equal;
gridres = 10;
xgrid = -100:gridres:100;
ygrid = -50:gridres:50;
% Preallocate a histogram:
hist2 = NaN(length(ygrid),length(xgrid));
for row = 1:length(ygrid)-1
for col = 1:length(xgrid)-1
% Get indices of all x,y points inside this grid cell:
ind = x>=xgrid(col)-gridres/2 & x<xgrid(col+1)-gridres/2 & y>=ygrid(row)-gridres/2 & y<ygrid(row+1)-gridres/2;
%
% Log the sum of number of lighting strikes inside the grid cell:
hist2(row,col) = sum(ind);
end
end
%
0 Commenti
Vedere anche
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!