halfeatenpie / hydrosanity Goto Github PK

better to use lattice than custom-build grid plots.

Some analyses should transform rainfall by sqrt.
This can be changed by the user by editing the plot call.

make regular
lag (a sub-period)
difference between two / diff
smooth
add series
fill gaps from another series
arbitrary transform like sqrt
change units (noted in another issue)

What steps will reproduce the problem?
1. install.packages("hydrosanity",dep=T)
2. library(hydrosanity)
3. hydrosanity()

What is the expected output? What do you see instead?

> library(hydrosanity)
Loading required package: playwith
Loading required package: lattice 
Loading required package: cairoDevice
Loading required package: gWidgetsRGtk2
Loading required package: gWidgets     
Loading required package: grid         
Hydrosanity: an interface for exploring hydrological time series. 
Version 0.8.76. (c) 2007 Felix Andrews <[email protected]>          
  GUI based on Rattle (c) 2006 [email protected]       
Type "hydrosanity()" to start the graphical user interface.       
Warning message:                                                  
In namespaceImportFrom(self, asNamespace(ns)) :                   
  replacing previous import: update                               
> hydrosanity()                                                   
Error in hydrosanity() : could not find function "gladeXMLNew"    

What version of the product are you using? On what operating system?

Version 0.8.76. / Ubuntu 9.04 / R2.10a

write tooltips for the GUI.
rattle has them.

Sorry, First post.

It has been several years since I got this program to load. Here is my 
understanding of what has happened since it was created:

Libglade has been completely replaced by GtkBuilder.  

The hydrosanity.glade file for the interface can be converted to a Gtk XML file 
with  libglade-convert or gtk-builder-convert, this throws several warnings 
about duplicate IDs, but no serious issues.

My understanding is that the software needs to be rewritten or redeveloped 
slightly for GtkBuilder. 

I have more rainfall data than has ever been collected in digital form on a 
small West African Republic. I also have the 10GB TRMM dataset which I feel 
would be best understood with good software if available.  

Given the spatial functions, the powerful time series capibilities, and 
graphical feedback to boot, I want to use hydrosanity as a base for my own 
software development. 

Step One: 

Using Ubuntu 10.4.1 LTS 

Step Two:

Used cran and ubuntu mirrors to install the usual suspects, 
i.e. 
sudo apt-get install r-cran-* 

Step Three:

sudo apt-get install subversion
svn checkout http://hydrosanity.googlecode.com/svn/trunk/ hydrosanity-read-only

Set Four:

sudo apt-get install glade
Not sure, but I think this gave me all the GtkBuilder goodies

Find the hydrosanity.glade file in "hydrosanity-read-only"

convert it with gtk-builder-convert

Step Five:

Try to figure out

http://developer.gnome.org/gtk/stable/gtk-migrating-GtkBuilder.html

Step Six:

within R I (~/freeman) ran:
>install.packages(pkgs="/home/freeman/Downloads/rattle-read-only/" repos=NULL)

After installing required libraries, calling 
>library(hydrosanity)
>hydrosanity() 

I get:

Error in checkPtrType(object, "GtkBuilder") : 
  object of class character isn't a GtkBuilder 

!!!!!

Which is better than file.choose V. choose.file or just can't find function 
gladeXMLNew

What steps will reproduce the problem?
1. Open MODSIM07.hydrosanity project
2. Go to Areal rainfall tab, select two or more gauges from list at left
3. Under "Areal rainfall from, gauges, weighted by area" choose either
"Thiessen Polygon" or "Calculate areal rainfall series"

What is the expected output? What do you see instead?

When choosing "Thiessen polygon" this is the error:

A command has failed. The error was:

cannot coerce type 'S4' to vector of type 'any'

The error occurred in: 

as.vector(x, mode)

The original command was: 

{
    tmp.poly <- hsp$catchment@polygons[[1]]@Polygons[[1]]
    tmp.subPolys <- arealSubPolygons(tmp.locs, boundary =
coordinates(tmp.poly))
}

When choosing "Calculate areal rainfall series" this is the error:

A command has failed. The error was:

cannot coerce type 'S4' to vector of type 'any'

The error occurred in: 

as.vector(x, mode)

The original command was: 

{
    tmp.poly <- hsp$catchment@polygons[[1]]@Polygons[[1]]
    tmp.subPolys <- arealSubPolygons(tmp.locs, boundary =
coordinates(tmp.poly))
    tmp.subAreas <- sapply(tmp.subPolys@polygons, getPolygonAreaSlot)
    tmp.areaFrac <- tmp.subAreas/tmp.poly@area
    tmp.polyNames <- getSpPPolygonsIDSlots(tmp.subPolys)
    tmp.data <- lapply(hsp$data[tmp.polyNames], window, hsp$timePeriod[1], 
        hsp$timePeriod[2])
    tmp.data <- lapply(tmp.data, quick.disaccumulate.timeblob)
    tmp.sync <- sync.timeblobs(tmp.data)
    tmp.synctime <- tmp.sync$Time
    tmp.sync <- tmp.sync[-1]
}

I think the error is caused by arealSubPolygons function 

What version of the product are you using? On what operating system?

Output from SessionInfo:

R version 2.8.1 (2008-12-22) 
i386-pc-mingw32 

locale:
LC_COLLATE=English_Australia.1252;LC_CTYPE=English_Australia.1252;LC_MONETARY=En
glish_Australia.1252;LC_NUMERIC=C;LC_TIME=English_Australia.1252

attached base packages:
[1] stats     graphics  grDevices datasets  grid      utils     methods  
base     

other attached packages:
 [1] gpclib_1.4-3         tripack_1.2-11       sp_0.9-32           
hydrosanity_0.8.76  
 [5] playwith_0.9-43      gWidgetsRGtk2_0.0-51 gWidgets_0.0-35     
cairoDevice_2.9-2   
 [9] lattice_0.17-20      rcom_2.0-4           rscproxy_1.0-12     

loaded via a namespace (and not attached):
[1] gridBase_0.4-3 reshape_0.8.2  RGtk2_2.12.9   tools_2.8.0   


Please provide any additional information below.

stl does a nice decomposition of seasonality, trend and remainder.

all metadata should be editable in one place.
copy to/from a spreadsheet (i.e. as tab-separated).
maybe edit as csv text.

unify side-bar icon view and main table?

simplest case is a single grid with long-term mean areal rainfall.
then use that to construct weightings for the gauges, and estimate areal
rainfall (daily) series as a weighted sum. (complications with missing values)

plot ccf (lag) over time, multiple series.
could be a shaded grid or a bubble plot for each series.

would be nice to be able to click on one year to show the pair of time
series for that year.

separate plots into separate log

Have permission from BOM to distribute rain gauge data.
Find out about stream flow data.
Catchment shapefile should be fine.

Thiessen polygon method: compute scale factors for rainfall sites, and sum.

Also need to convert streamflow to mm for water balance stuff.

write a vignette.
etc. documentation

Need to formally attach units to time series.
Allow converting units from dataset tab.

is.numeric on POSIXt or Date is FALSE in R >= 2.6

convert negative data values to NA

need a more standard time series format.

zoo is pretty standard, and can handle irregular and regular series. can be
coerced to ts.

need to handle intersect / union of multiple time series.

need to handle quality codes - store as numeric and convert to factor when
needed?

how to keep quality codes with multiple time series? can't do it unless
store series separately.

need to adapt for representing time periods - end() / window()