ottadini

'''
Constants for use in processing scripts

For Windows paths use the raw string formatter, e.g. 

# my_path = r"C:\"
'''
# Will be different in QGIS install vs OSGEO4W install
QGIS_PLUGINS_PATH = r"C:\OSGeo4W64\apps\qgis\python\plugins"
WORKING_RASTER = "IslandOfKahoolawe.sdat"

ottadini

# Get a vtk.PolyData object for the input
pdi = self.GetPolyDataInput()
# Get a vtk.PolyData object for the output
pdo = self.GetPolyDataOutput()

numPoints = pdi.GetNumberOfPoints()

# Points for the line:
newPoints = vtk.vtkPoints()
for i in range(0, numPoints):

ottadini

# Copying the kohonen:::plot.kohcodes routine, for equivalent kmeans stars diagram
plot.centers <- function (x, main) 
{
  palette.name <- terrain.colors

  codes <- x$codes
  nvars <- ncol(codes)
  maxlegendcols <- 3
  
  codeRendering <- "segments"

ottadini

ottadini

wrap <- function(string){
  paste(strwrap(string), collapse="\n")
}

ottadini

# Cannot remember the source for this code, it isn't mine. 


genHexGrid <- function(dx, ll = c(0, 0), ur = c(1, 1))
{
    dy <- sqrt(3) * dx / 2
    
    x <- seq(ll[1], ur[1] - dx / 2, dx)
    y <- seq(ll[2], ur[2], dy)
    
    y <- rep(y, each = length(x))

ottadini

# I cannot remember the source for this code. It isn't mine. Probably somewhere on stackoverflow.




require ("kohonen")
require ("latticeExtra")
require ("deldir")
data(wines)
som.wines <- som (scale (wines), grid = somgrid(5, 5, "hexagonal"))

df <- as.data.frame (cbind (som.wines$grid$pts, som.wines$codes))
levelplot (alcohol ~ x * y, data = df, panel = panel.voronoi, aspect = "iso")

ottadini

plot.multi.dens <- function(s){
# From: http://onertipaday.blogspot.com.au/2007/09/plotting-two-or-more-overlapping.html
# the input s MUST be a numeric list
# e.g. 
# rating <- rnorm(200)
# rating2 <- rnorm(200, mean=.8)
# plot.multi.dens( list(rating, rating2)) 
    junk.x = NULL
    junk.y = NULL

ottadini

# Base graphics

# Create a PDF for a square shaped plot (e.g. for predictions vs observations regression)

pdf(file="pred-vs-meas_MLR.pdf", width=5, height=5, family="Times")
op <- par(pty='s')  # plot type = square
plot(X, Y, 
     xlab=Xlabel, ylab=Ylabel, ylim=c(0, 3), xlim=c(0, 3)) #
abline(fit, lty=2)
par(op)

ottadini