Author Topic: Cascade method Rainfall disaggragation  (Read 10953 times)

ujjwalsaha

  • Team river
  • ***
  • Karma: +15/-0
    • View Profile
Cascade method Rainfall disaggragation
« on: August 25, 2012, 05:09:36 PM »
i am now working on rainfall disaggregation by cascade method.. i am not able to understand the whole process clearly. whether anybody have worked on it ? ...can anyone send a flowchart giving step by step process..also the theory in simple way so that i can understand...i will write the code in matlab....any reply will be appreciated..


Sat Kumar Tomer

  • Administrator
  • Team sea
  • *****
  • Karma: +116/-0
    • View Profile
    • ambhas
  • Institute : IISc
  • Programming language : R, Python and Matlab
Re: Cascade method Rainfall disaggragation
« Reply #1 on: August 25, 2012, 11:43:22 PM »
I have written code in Python and Matlab for the disaggregation of rainfall. The python code is available at the ambhas pythong library (http://code.google.com/p/ambhas/. Unfortunately the documentation of the Python code is not available yet. I will try to make documentation for this soon. I wrote the code in Matlab also, but I will have to find it and clean it so that others can understand. I will try to attach the Matlab code soon.

For the theory and implementation part, you can read my article and presentation. The article is attached and the presentation can be accessed from the presentation page of the ambhas website (http://ambhas.com/pre.php).

ujjwalsaha

  • Team river
  • ***
  • Karma: +15/-0
    • View Profile
Re: Cascade method Rainfall disaggragation
« Reply #2 on: August 26, 2012, 10:44:34 AM »
Thank you so much sat kumar...i will read your paper and will get back to you if i have any question...
it would be very helpful for me if you upload the matlab programme today or tomorrow ( if possible)...it will give me  guidance...

Sat Kumar Tomer

  • Administrator
  • Team sea
  • *****
  • Karma: +116/-0
    • View Profile
    • ambhas
  • Institute : IISc
  • Programming language : R, Python and Matlab
Re: Cascade method Rainfall disaggragation
« Reply #3 on: August 26, 2012, 05:17:55 PM »
You are welcome. The matlab files are attached. rain_disagg.m is the main file that you need to edit, other files contain functions used by the main file.

ujjwalsaha

  • Team river
  • ***
  • Karma: +15/-0
    • View Profile
Re: Cascade method Rainfall disaggragation
« Reply #4 on: August 26, 2012, 08:19:58 PM »
great work sat kumar...it is real help for me...

mohsen.pourreza

  • Team tributary
  • *
  • Karma: +0/-0
    • View Profile
Re: Cascade method Rainfall disaggragation
« Reply #5 on: March 16, 2013, 04:38:29 PM »
You are welcome. The matlab files are attached. rain_disagg.m is the main file that you need to edit, other files contain functions used by the main file.

Dear Dr. Sat,
During the running of your m.file I faced with a error in second line of "rain.disagg.m" file
data=dlmread('Rainfall.txt'); It needs to relate a txt file named rainfall.txt,
May you please give a rainfall file so it would be able to pass this error? Actually I need to know about format of rainfall file ..
All the best,
Mohsen

Sat Kumar Tomer

  • Administrator
  • Team sea
  • *****
  • Karma: +116/-0
    • View Profile
    • ambhas
  • Institute : IISc
  • Programming language : R, Python and Matlab
Re: Cascade method Rainfall disaggragation
« Reply #6 on: April 19, 2013, 02:41:14 PM »
The rainfall.txt file is attached.
First column is the day
second column is the month
third column is the year
fourth column is the rainfall.

Only fourth column (rainfall) is used in the program.
This file is for the daily data, but you can use the same for the hourly (or at some other frequency) also. Just put your rainfall in the forum column, entries in the other column do not matter.

Nachiketa

  • Team tributary
  • *
  • Karma: +2/-0
    • View Profile
    • Google Scholar
Re: Cascade method Rainfall disaggragation
« Reply #7 on: April 24, 2013, 11:06:34 AM »
Dear Dr.Sat,
I have read your articular and all the MATLAB code. It is really a great work. Congratulation.
I have a question that can we disaggregate Monthly rainfall to daily rainfall using your code?
If it is can you please tell what will be the arrangement in Rainfall.txt file.
Thanks in advance
Nachiketa

Sat Kumar Tomer

  • Administrator
  • Team sea
  • *****
  • Karma: +116/-0
    • View Profile
    • ambhas
  • Institute : IISc
  • Programming language : R, Python and Matlab
Re: Cascade method Rainfall disaggragation
« Reply #8 on: April 24, 2013, 12:43:59 PM »
Hi,
You can use the code in the same way.
Just put the monthly data in the rainfall column.

But be careful about one thing, that sometimes there is breakup in the relationship.
I will try to find some reference on this issue.

Nachiketa

  • Team tributary
  • *
  • Karma: +2/-0
    • View Profile
    • Google Scholar
Re: Cascade method Rainfall disaggragation
« Reply #9 on: April 24, 2013, 04:26:12 PM »
Hi Dr.Sat,
Thanks for your quick reply. As per your suggestion, we have to modify the Rainfall.txt file. But I have some specific quires which are below:
 In this file the first column is for daily. Therefore, for each year and for each month we should have 30/31 rainfall inputs. But for our case, for each month only one rainfall value is exist. For example, 10 mm rainfall for January 2008. In that case the dimension of the data file will reduce.So can we repeat this rainfall amount for 30 times in this file.
Kindly guide in this regard.
Thanks
Nachiketa

Sat Kumar Tomer

  • Administrator
  • Team sea
  • *****
  • Karma: +116/-0
    • View Profile
    • ambhas
  • Institute : IISc
  • Programming language : R, Python and Matlab
Re: Cascade method Rainfall disaggragation
« Reply #10 on: April 24, 2013, 05:14:09 PM »
Quote
The rainfall.txt file is attached.
First column is the day
second column is the month
third column is the year
fourth column is the rainfall.

Only fourth column (rainfall) is used in the program.
This file is for the daily data, but you can use the same for the hourly (or at some other frequency) also. Just put your rainfall in the forum column, entries in the other column do not matter.

Just put the monthly rainfall in the fourth column, and do the analysis. Don't repeat the rainfall.
Make sure that you have long enough series. I think around 1000 data values should be ok.

You may like to read the article attached above and the presentation (which can be downloaded from the above given link).
« Last Edit: April 24, 2013, 05:15:40 PM by Sat Kumar Tomer »

Nachiketa

  • Team tributary
  • *
  • Karma: +2/-0
    • View Profile
    • Google Scholar
Re: Cascade method Rainfall disaggragation
« Reply #11 on: April 26, 2013, 11:23:53 AM »
Dear Dr.Sat Kumar,
Many thanks for your reply.
I am trying to modify the file.
I am also reading your paper thoroughly so that i can have better understand.
best wishes
Nachiketa

mohsen.pourreza

  • Team tributary
  • *
  • Karma: +0/-0
    • View Profile
Re: Cascade method Rainfall disaggragation
« Reply #12 on: May 27, 2013, 08:43:06 PM »
Dear Dr. Kumar,
Regarding the plots in Fig 4 from your manuscript please explain how can I produce CDF of "measured rainfall". As you know, the code just plots CDF of "simulated rainfall". Then please send if there is additional code which produces it.
Finally I would like know more about RMC method, then please put some papers which describe RMC more .
All the best,
Mohsen
« Last Edit: May 27, 2013, 08:48:48 PM by mohsen.pourreza »

Sat Kumar Tomer

  • Administrator
  • Team sea
  • *****
  • Karma: +116/-0
    • View Profile
    • ambhas
  • Institute : IISc
  • Programming language : R, Python and Matlab
Re: Cascade method Rainfall disaggragation
« Reply #13 on: May 27, 2013, 09:38:15 PM »
The code only disaggregate 24 hour rainfall into 12 hour rainfall.
The measured rainfalls are given at 1, 2, 4, 8, 16 .. hour.
You need to take both the rainfall at the same time step.

I am taking both the rainfall at the 12 hour time step.
You can modify the last section of the code as:
Code: [Select]
RR_12=sum(reshape(RR_1,12,length(RR_1)/12));
figure
[f,x] = ecdf(RR_pre12);
[f1,x1] = ecdf(RR_1);
stairs(x,f,'b');
hold on
stairs(x1,f1,'r');
xlabel('Rainfall -- 12 h')
ylabel('CDF')
legend('Simulated','Measured','Location','Best')
hold off
[\code]

mohsen.pourreza

  • Team tributary
  • *
  • Karma: +0/-0
    • View Profile
Re: Cascade method Rainfall disaggragation
« Reply #14 on: May 27, 2013, 10:45:46 PM »
Thanks for your reply..I still need your guide..
1- "RR_12=sum(reshape(RR_1,12,length(RR_1)/12))" faced with error, because RR-1 is a matrix 704*1 then 704/12 is not a integer value.
2- Is "RR_pre12" means 12-hourly simulated rainfall?if so then why its dimension is 352 while we have 704 daily rainfall i.e I think it can be changed to a matrix with dimension (704*2,1)=(1408,1)?(there is two 12-hourly step in a daily or 24-hourly rainfall.)
3- you said "The measured rainfalls are given at 1, 2, 4, 8, 16 .. hour." but in the rainfall.txt file there is just daily rainfall,can I have a sample data file containing 1,2,4..hour??
All the best,
Mohsen