I want to create a scatter plot of rainfall for specific days. How would I do that from the code below?

1 次查看（过去 30 天）

Thomas 2014-3-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/123098-i-want-to-create-a-scatter-plot-of-rainfall-for-specific-days-how-would-i-do-that-from-the-code-bel

评论： Star Strider 2014-3-26

% [jday, year, pressure, temp, wspeed, wdirn, dewtemp, wgust, vis, month, day, hours, mins] = ReadRain(fname) % % Create Rainfall Plots for the Specified Period Based on Different Wind % Speed Categories % % INPUTS: % fname : name of file to read % % OUTPUTS: % jday : time (in Julian days since 1st Jan) % year : year % temp : temperature (degC) % wspeed : wind speed (knots) % wdirn : wind direction (degrees) % dewtemp : dewpoint temperature (degC) % precip : rainfall (mm) % % % % Author : Thomas Dent-Jones <mailto:ee10tdj@leeds.ac.uk ee10tdj@leeds.ac.uk> % Created : 19/02/2014 % Last modified : 19/02/2014 %

    [jday, year, precip, month, day, hours, mins, id, cnt] = ReadRain('Rain_daily_Cornwall_2013_v2.csv');
    [jdayw, yearw, pressure, temp, wspeed, wdirn, dewtemp, wgust, vis, monthw, dayw, hoursw, minsw, station] = WindReader3808('Winds_hourly_Cornwall_2013_v2.csv');
    stnlist = unique(id);
    n = length(stnlist);
    totprecip = zeros(n,1);
    ndays = zeros(n,1);
    day1 = 91;
    day2 = 244;
    days = day1:day2;
    precip2 = NaN*ones(length(days),n);
   for i=1:n
        jday1 = jday(id==stnlist(i));
        precip1 = precip(id==stnlist(i));
        for j=1:length(precip1)
            precip2(days==jday1(j),i)=precip1(j);
        end
        totprecip(i) = sum(precip1);
        ndays(i) = length(precip1);
   end
 %Creating the scatter plot.  
   [height, lon, lat, sid] = Station('Cornwall_station_details2_with_height.csv');
   x = zeros(n,1);
   y = zeros(n,1);
   h = zeros(n,1);
   for i=1:n
      x(i) = lon(sid == stnlist(i));
      y(i) = lat(sid == stnlist(i));
      h(i) = height(sid == stnlist(i));
   end
   axis([-5.6 -4.2 50 50.9]) 
   img = imread('Cornwall.jpg');             %# Load a sample image
   hold on;                                 %# Add to the plot
   image([-5.6 -4.2], [50.9 50],img);       %# Plot the image
   scatter(y(jday=214&jday=215),x(jday=214&jday=215),100,totprecip(jday=214&jday=215),'filled','MarkerEdgeColor','black');
   colorbar;
   caxis([120 500])
   xlabel('longitude');
   ylabel('latitude');
   %days=214&days=215

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Star Strider 2014-3-26

The important issue here is not your code but what your data look like and what you want to plot. Obviously one independent variable is time (date and possibly hour) and one dependent variable is rainfall. That may be all you need, but you can plot rainfall against one other independent variable if you want to.

You would do that from the output of the code you posted, if I’m reading your question correctly.

请先登录，再进行评论。

请先登录，再回答此问题。