how to make monthly graph from daily file data

조회 수: 5 (최근 30일)

Soni huu 2012년 7월 6일

0
링크

이 질문에 대한 바로 가기 링크

https://kr.mathworks.com/matlabcentral/answers/42885-how-to-make-monthly-graph-from-daily-file-data

thank to per isakson

from this code (daily/1 file) can u make monthly graph (30file data)?? =====

 function    RainData = ReadManySoniData( folder_name, file_spec )
      sad = dir( fullfile( folder_name, file_spec ) );
      RainData = struct([]);
      for sa = transpose( sad )
          RainData = cat( 2, RainData, ReadOneSoniData( folder_name, sa.name ) );
      end
      [ dummy, ixs ] = sort( [ RainData(:).DayNumber ] );
      RainData = RainData( ixs );
end
function    rain_data = ReadOneSoniData( folder_name, file_name )
      fid = fopen( fullfile( folder_name, file_name ), 'r' );
      if not( fid >= 3 )
          error( 'ReadOneSoniData:NoFileFound'     ...
               , 'Cannot find file "%s"'           ...
               , fullfile( folder_name, file_name ) )
      end    
      cac = textscan( fid, '%s', 'Whitespace','', 'Delimiter','\n' );
      fclose( fid );
      cac = cac{:};
      tmp = cellfun( 'isempty'    ...
              ,   regexp( cac, '\d{2}:\d{2}:\d{2}\s+C\s*$' ) );
      isc = not( tmp );
      tmp = cellfun( 'isempty'    ...
              ,   regexp( cac, '\d{2}:\d{2}:\d{2}\s+\*\*\s+----' ) );
      iss = not( tmp );
      cac( isc | iss ) = [];
      str = transpose( char( cac ) );
      nl  = sprintf('\n');
      str = cat( 1, str, repmat( nl(:), [length(nl),size(str,2)] ) );
      cac = cell(1,9);
      [cac{:}] = strread( str, '%8c%2c%4f%7f%4c%4u%4u%4u+%2u'      ...
                      ,        'delimiter', ' ', 'whitespace', ''  );
      try
          date_vec            = nan(1,3);
          date_vec( [2,3,1] ) = sscanf( file_name, '%2u-%2u-%4u%*s' );
      catch me
          if strcmp( me.identifier, 'MATLAB:index_assign_element_count_mismatch' )
              warning( 'ReadOneSoniData:CannotParseFileName'      ...
                   , 'Cannot extract a date from file name: "%s"' ...
                   , file_name                                    )
              rain_data = struct([]);
              return
          else
              rethrow( me )
          end
      end
      str = transpose( char( cac{1} ) );
      vec = nan( size(str,2), 3 );
      [ vec(:,1), vec(:,2), vec(:,3) ]                                ...
      =   strread( str, '%2u:%2u:%2u', 'delimiter','','whitespace','' );
      rain_data.Created   = datestr( now, 'yyyy-mm-dd HH:MM:SS' );
      rain_data.DataFile  = fullfile( folder_name, file_name );
      rain_data.Datevec   = [ repmat( date_vec, [size(vec,1),1] ), vec ];
      rain_data.DayNumber = datenum( date_vec );     
      rain_data.Rain      = cac{3};
      rain_data.DailyRain = sum( rain_data.Rain );
      % and more as you see fit. 
end

댓글 수: 2
없음 표시없음 숨기기

Soni huu 2012년 7월 6일

편집: Soni huu 2012년 7월 6일

this is sample yearly data (362 file data) (2.11m)

http://www.4shared.com/zip/7SuYaN2y/2011_2.html?

per isakson 2012년 7월 6일

"from this code (daily/1 file) can u make monthly graph (30file data)??"

Put more effort in describing what you need!

댓글을 달려면 로그인하십시오.

이 질문에 답변하려면 로그인하십시오.

채택된 답변

per isakson 2012년 7월 6일

2
링크

이 답변에 대한 바로 가기 링크

https://kr.mathworks.com/matlabcentral/answers/42885-how-to-make-monthly-graph-from-daily-file-data#answer_52712

편집: per isakson 2012년 7월 6일

Here is a function that returns total monthly rain. Try

>> mr = MonthlyRain( RainData );
>> plot( mr(1).Rain, 'd' );
>> bar( mr.Rain );

The values of the monthly rain could they be correct?

function    monthly_rain = MonthlyRain( RainData )
      day_number  = [ RainData(:).DayNumber ];
      month_number= month( day_number ); 
      year_number = year( day_number );
      year_list   = unique( year_number );
      monthly_rain = struct( 'Year', num2cell( year_list ), 'Rain', nan(12,1) );  
      ix_yy = 0;
      for yy = year_list
          is_yy = ( yy == year_number );
          ix_yy = ix_yy + 1;
          for mm = 1 : 12
              is_mm = ( mm == month_number );
              is_ym = ( is_yy & is_mm );
              if any( is_ym )
                  monthly_rain(ix_yy).Rain(mm) = sum([RainData( is_ym ).DailyRain]);
              end
          end
      end
  end

댓글 수: 63
이전 댓글 61개 표시이전 댓글 61개 숨기기

per isakson 2012년 7월 13일

Do you have problem reading the file, '12-19-2010.dat'?

I have successfully red the 2011 data, which I downloaded some days ago. However, I will not download more data. I've run out of time.

Add these lines

    fid = fopen( 'c:\temp\ReadRainDataFailures.log', 'a' );
    if fid >= 3
        fprintf(fid, '%s | %s | %s\n',datestr(now),folder_name, sa.name );
        fclose( fid );
    end

after

% disp( le.stack(1) )

That will give you a list of the files, which cannot be red. You might want to change the name and folder of the file. Do you have a "c:\temp"?

If you run into specific problems try make a question at Answers.

Soni huu 2012년 7월 13일

the code work now......... yes i have c:\temp

댓글을 달려면 로그인하십시오.

추가 답변 (1개)

per isakson 2012년 7월 6일

3
링크

이 답변에 대한 바로 가기 링크

https://kr.mathworks.com/matlabcentral/answers/42885-how-to-make-monthly-graph-from-daily-file-data#answer_52717

Here is a function that collects total daily rain for one month at a time. Try

>> [ day_number, daily_rain ] = DailyRain( RainData, 2011, 11 );
>> plot( day_number, daily_rain )
>> figure, plot( day_number, daily_rain, '.' )
>> figure, plot( day_number, '.' )

As is and all that! You must check the the values. Missing data might cause surprises.

function [ day_number, daily_rain ] = DailyRain( RainData, year_number, month_number )
      day_number  = [ RainData(:).DayNumber ];
      is_yy = ( year( day_number )  == year_number );
      is_mm = ( month( day_number ) == month_number );
      is_ym = ( is_yy & is_mm );
      if any( is_ym )
          daily_rain  = [ RainData( is_ym ).DailyRain ];
          day_number  = [ RainData( is_ym ).DayNumber ];
      else
          daily_rain  = [ ];
          day_number  = [ ];
      end 
  end