move some of the detection and landmarking code to plot_utils

Ben Klein
1 parent b205c6b6
Showing 3 changed files with 123 additions and 143 deletions
docs/docs/api_docs/c_api/functions.md
openbr/core/plot.cpp
share/openbr/plotting/plot_utils.R
@@ -587,7 +587,7 @@ Returns names and parameters for the requested objects. Each object is newline s
  
 ---
  
-## I
+## br_plot
  
 Renders recognition performance figures for a set of **.csv** files created by [br_eval](#br_eval).
  
@@ -232,8 +232,8 @@ bool Plot(const QStringList &amp;files, const File &amp;destination, bool show)
     p.file.write(qPrintable(QString(plot).arg("CMC", toRList(optMap["cmcOptions"]), "FALSE")));
     p.file.write("plot <- plotSD(sdData=SD)\nplot\n");
     p.file.write("plot <- plotBC(bcData=BC)\nplot\n");
-    p.file.write("plot <- plotERR(errData=ERR)\nplot\n\n");
-    p.file.write("plotEERSamples(imData=IM, gmData=GM");
+    p.file.write("plot <- plotERR(errData=ERR)\nplot\n");
+    p.file.write("plotEERSamples(imData=IM, gmData=GM)\n\n");
  
     return p.finalize(show);
 }
@@ -295,23 +295,19 @@ bool PlotDetection(const QStringList &amp;files, const File &amp;destination, bool show)
     QString plot = "plot <- plotLine(lineData=%1, options=list(%2), flipY=%3, geometry=%4)\nplot\n";
     foreach (const QString &type, QStringList() << "Discrete" << "Continuous") {
         optMap["rocOptions"].set("title", type);
-        p.file.write(qPrintable(QString(plot).arg(type + "ROC", toRList(optMap["rocOptions"]), "FALSE", plotType)));
+        p.file.write(qPrintable(QString(plot).arg(type + "ROC", toRList(optMap["rocOptions"]), "FALSE", "\"" + plotType + "\"")));
     }
  
     foreach (const QString &type, QStringList() << "Discrete" << "Continuous") {
         optMap["prOptions"].set("title", type);
-        p.file.write(qPrintable(QString(plot).arg(type + "PR", toRList(optMap["prOptions"]), "FALSE", plotType)));
+        p.file.write(qPrintable(QString(plot).arg(type + "PR", toRList(optMap["prOptions"]), "FALSE", "\"" + plotType + "\"")));
     }
-
-    p.file.write(qPrintable(QString("qplot(X, data=Overlap, geom=\"histogram\", position=\"identity\", xlab=\"Overlap\", ylab=\"Frequency\")") +
-                            QString(" + theme_minimal() + scale_x_continuous(minor_breaks=NULL) + scale_y_continuous(minor_breaks=NULL) + theme(axis.text.y=element_blank(), axis.ticks=element_blank(), axis.text.x=element_text(angle=-90, hjust=0))") +
-                            (p.major.size > 1 ? (p.minor.size > 1 ? QString(" + facet_grid(%2 ~ %1, scales=\"free\")").arg(p.minor.header, p.major.header) : QString(" + facet_wrap(~ %1, scales = \"free\")").arg(p.major.header)) : QString()) +
-                            QString(" + theme(aspect.ratio=1, legend.position=\"bottom\")\n\n")));
+    p.file.write("plot <- plotOverlap(overlapData=Overlap)\nplot\n");
  
     p.file.write(qPrintable(QString("ggplot(AverageOverlap, aes(x=%1, y=%2, label=round(X,3)), main=\"Average Overlap\") + geom_text() + theme_minimal()").arg(p.minor.size > 1 ? p.minor.header : "'X'", p.major.size > 1 ? p.major.header : "'Y'") +
                             QString("%1%2\n\n").arg(p.minor.size > 1 ? "" : " + xlab(NULL)", p.major.size > 1 ? "" : " + ylab(NULL)")));
  
-    p.file.write(qPrintable(QString("ggplot(AverageOverlap, aes(x=%1, y=%2, fill=X)) + geom_tile() + scale_fill_continuous(\"Average Overlap\") + theme_minimal()").arg(p.minor.size > 1 ? p.minor.header : "'X'", p.major.size > 1 ? p.major.header : "'Y'") +
+    p.file.write(qPrintable(QString("ggplot(AverageOverlap, aes(x=%1, y=%2, fill=X)) + geom_tile() + scale_fill_continuous(\"Average Overlap\", guide=FALSE) + theme_minimal()").arg(p.minor.size > 1 ? p.minor.header : "'X'", p.major.size > 1 ? p.major.header : "'Y'") +
                             QString("%1%2\n\n").arg(p.minor.size > 1 ? "" : " + xlab(NULL)", p.major.size > 1 ? "" : " + ylab(NULL)")));
  
     return p.finalize(show);
@@ -322,66 +318,10 @@ bool PlotLandmarking(const QStringList &amp;files, const File &amp;destination, bool sho
     qDebug("Plotting %d landmarking file(s) to %s", files.size(), qPrintable(destination));
     RPlot p(files, destination);
     p.file.write("\nformatData(type=\"landmarking\")\n\n");
-
-    p.file.write(qPrintable(QString("\nreadData <- function(data) {\n\texamples <- list()\n"
-                                    "\tfor (i in 1:nrow(data)) {\n"
-                                    "\t\tpath <- data[i,1]\n"
-                                    "\t\tvalue <- data[i,2]\n"
-                                    "\t\tfile <- unlist(strsplit(path, \"[.]\"))[1]\n"
-                                    "\t\text <- unlist(strsplit(path, \"[.]\"))[2]\n"
-                                    "\t\tif (ext == \"jpg\" || ext == \"JPEG\" || ext == \"jpeg\" || ext == \"JPG\") {\n"
-                                    "\t\t\timg <- readJPEG(path)\n"
-                                    "\t\t} else if (ext == \"PNG\" || ext == \"png\") {\n"
-                                    "\t\t\timg <- readPNG(path)\n"
-                                    "\t\t} else if (ext == \"TIFF\" || ext == \"tiff\" || ext == \"TIF\" || ext == \"tif\") { \n"
-                                    "\t\t\timg <- readTIFF(path)\n"
-                                    "}else {\n"
-                                    "\t\t\tnext\n"
-                                    "\t\t}\n"
-                                    "\t\texample <- list(file = file, value = value, image = img)\n"
-                                    "\t\texamples[[i]] <- example\n"
-                                    "\t}\n"
-                                    "\treturn(examples)\n"
-                                    "}\n")));
-
-    p.file.write(qPrintable(QString("\nlibrary(jpeg)\n"
-                                    "library(png)\n"
-                                    "library(grid)\n")));
-
-    p.file.write(qPrintable(QString("\nplotImage <- function(image, title=NULL, label=NULL) { \n"
-                                    "\tp <- qplot(1:10, 1:10, geom=\"blank\") + annotation_custom(rasterGrob(image$image), xmin=-Inf, xmax=Inf, ymin=-Inf, ymax=Inf) + theme(axis.line=element_blank(), axis.title.y=element_blank(), axis.text.x=element_blank(), axis.text.y=element_blank(), line=element_blank(), axis.ticks=element_blank(), panel.background=element_blank()) + labs(title=title) + xlab(label)\n"
-                                    "\treturn(p)"
-                                    "}\n")));
-
-    p.file.write(qPrintable(QString("\nsample <- readData(Sample) \n"
-                                    "rows <- sample[[1]]$value\n"
-                                    "algs <- unique(Box$%1)\n"
-                                    "algs <- algs[!duplicated(algs)]\n"
-                                    "print(plotImage(sample[[1]],\"Sample Landmarks\",sprintf(\"Total Landmarks: %s\",sample[[1]]$value))) \n"
-                                    "if (nrow(EXT) != 0 && nrow(EXP)) {\n"
-                                    "\tfor (j in 1:length(algs)) {\n"
-                                    "\ttruthSample <- readData(EXT[EXT$. == algs[[j]],])\n"
-                                    "\tpredictedSample <- readData(EXP[EXP$. == algs[[j]],])\n"
-                                    "\t\tfor (i in 1:length(predictedSample)) {\n"
-                                    "\t\t\tmultiplot(plotImage(predictedSample[[i]],sprintf(\"%s\\nPredicted Landmarks\",algs[[j]]),sprintf(\"Average Landmark Error: %.3f\",predictedSample[[i]]$value)),plotImage(truthSample[[i]],\"Ground Truth\\nLandmarks\",\"\"),cols=2)\n"
-                                    "\t\t}\n"
-                                    "\t}\n"
-                                    "}\n").arg(p.major.size > 1 ? p.major.header : (p.minor.header.isEmpty() ? p.major.header : p.minor.header))));
-
-    p.file.write(qPrintable(QString("\n"
-                 "# Code to format error table\n"
-                 "StatBox <- summarySE(Box, measurevar=\"Y\", groupvars=c(\"%1\",\"X\"))\n"
-                 "OverallStatBox <- summarySE(Box, measurevar=\"Y\", groupvars=c(\"%1\"))\n"
-                 "mat <- matrix(paste(as.character(round(StatBox$Y, 3)), round(StatBox$ci, 3), sep=\" \\u00b1 \"),nrow=rows,ncol=length(algs),byrow=FALSE)\n"
-                 "mat <- rbind(mat, paste(as.character(round(OverallStatBox$Y, 3)), round(OverallStatBox$ci, 3), sep=\" \\u00b1 \"))\n"
-                 "mat <- rbind(mat, as.character(round(NormLength$Y, 3)))\n"
-                 "colnames(mat) <- algs\n"
-                 "rownames(mat) <- c(seq(0,rows-1),\"Aggregate\",\"Average IPD\")\n"
-                 "ETable <- as.table(mat)\n").arg(p.major.size > 1 ? p.major.header : (p.minor.header.isEmpty() ? p.major.header : p.minor.header))));
-
-    p.file.write(qPrintable(QString("\n"
-                       "print(textplot(ETable))\n"
-                       "print(title(\"Landmarking Error Rates\"))\n")));
+    p.file.write(qPrintable(QString("algs <- uniqueBox$%1)\n").arg(p.major.size > 1 ? p.major.header : (p.minor.header.isEmpty() ? p.major.header : p.minor.header))));
+    p.file.write("algs <- algs[!duplicated(algs)]\n");
+    p.file.write("plotLandmarkSamples(samples=sample, expData=EXP, extData=EXT)\n");
+    p.file.write("plotLandmarkTables(tableData=Box)\n");
  
     p.file.write(qPrintable(QString("ggplot(Box, aes(Y,%1%2))").arg(p.major.size > 1 ? QString(", colour=%1").arg(p.major.header) : QString(),
                                                                     p.minor.size > 1 ? QString(", linetype=%1").arg(p.minor.header) : QString()) +
@@ -3,6 +3,9 @@ library(&quot;ggplot2&quot;)
 library("gplots")
 library("reshape")
 library("scales")
+library("jpeg")
+library("png")
+library("grid")
  
 # Code to format FAR values
 far_names <- list('0.001'="FAR = 0.1%", '0.01'="FAR = 1%")
@@ -53,6 +56,24 @@ plotTable &lt;- function(tableData=NULL, name=NULL, labels=NULL) {
     }
 }
  
+plotLandmarkTables <- function(tableData=NULL) {
+    if(majorSize > 1) {
+        var <- majorHeader
+    } else {
+        if(minorHeader == "") var <- majorHeader else var <- minorHeader
+    }
+    StatBox <- summarySE(tableData, measurevar="Y", groupvars=c(var,"X"))
+    OverallStatBox <- summarySE(tableData, measurevar="Y", groupvars=c(var))
+    mat <- matrix(paste(as.character(round(StatBox$Y, 3)), round(StatBox$ci, 3), sep=" \u00b1 "), nrow=rows, ncol=length(algs), byrow=FALSE)
+    mat <- rbind(mat, paste(as.character(round(OverallStatBox$Y, 3)), round(OverallStatBox$ci, 3), sep=" \u00b1 "))
+    mat <- rbind(mat, as.character(round(NormLength$Y, 3)))
+    colnames(mat) <- algs
+    rownames(mat) <- c(seq(0, rows-1), "Aggregate","Average IPD")
+    ETable <- as.table(mat)
+    print(textplot(ETable))
+    print(title("Landmarking Error Rates"))
+}
+
 plotLine <- function(lineData=NULL, options=NULL, flipY=FALSE, geometry="line") {
     textSize <- if("textSize" %in% names(options)) as.numeric(options$textSize) else 12
     p <- qplot(X, if(flipY) 1-Y else Y, data=lineData, main=options$title, geom=geometry, size=if("size" %in% names(options)) I(as.numeric(options$size)) else I(.5), colour=if(majorSize > 1) factor(eval(parse(text=majorHeader))) else NULL, linetype=if(minorSize > 1) factor(eval(parse(text=minorHeader))) else NULL, xlab=options$xTitle, ylab=options$yTitle) + theme_minimal()
@@ -138,6 +159,20 @@ plotERR &lt;- function(errData=NULL) {
     return(p)
 }
  
+plotOverlap <- function(overlapData=NULL) {
+    p <- qplot(X, data=overlapData, geom="histogram", position="identity", xlab="Overlap", ylab="Frequency")
+    p <- p + theme_minimal() + scale_x_continuous(minor_breaks=NULL) + scale_y_continuous(minor_breaks=NULL) + theme(axis.text.y=element_blank(), axis.ticks=element_blank(), axis.text.x=element_text(angle=-90, hjust=0))
+    if(majorSize > 1) {
+        if(minorSize) {
+            p <- p + facet_grid(facets=as.formula(paste(minorHeader, "~", majorHeader)), scales="free")
+        } else {
+            p <- p + facet_wrap(facets=as.formula(paste("~", majorHeader)), scales="free")
+        }
+    }
+    p <- p + theme(aspect.ratio=1, legend.position="bottom")
+    return(p)
+}
+
 formatData <- function(type="eval") {
     if (type == "eval") {
 	    # Split data into individual plots
@@ -182,7 +217,6 @@ formatData &lt;- function(type=&quot;eval&quot;) {
         ContinuousPR <<- data[grep("ContinuousPR",data$Plot),-c(1)]
         Overlap <<- data[grep("Overlap",data$Plot),-c(1)]
         AverageOverlap <<- data[grep("AverageOverlap",data$Plot),-c(1)]
-
     } else if (type == "landmarking") {
         # Split data into individual plots
         Box <<- data[grep("Box",data$Plot),-c(1)]
@@ -194,10 +228,13 @@ formatData &lt;- function(type=&quot;eval&quot;) {
         EXP <<- data[grep("EXP",data$Plot),-c(1)]
         EXP$X <<- as.character(EXP$X)
         NormLength <<- data[grep("NormLength",data$Plot),-c(1)]
+        sample <<- readImageData(Sample)
+        rows <<- sample[[1]]$value
     }
 }
  
 summarySE <- function(data=NULL, measurevar, groupvars=NULL, na.rm=FALSE, conf.interval=0.95, .drop=TRUE) {
+    # derived from http://www.cookbook-r.com/Manipulating_data/Summarizing_data/
 	require(plyr)
  
 	length2 <- function (x, na.rm=FALSE) {
@@ -246,85 +283,88 @@ multiplot &lt;- function(..., plotlist=NULL, cols) {
 plotEERSamples <- function(imData=NULL, gmData=NULL) {
     if(nrow(imData) == 0) return()
  
-    library(jpeg)
-    library(png)
-    library(grid)
-
-    for (i in 1:nrow(imData)) {
-        score <- imData[i,1]
-        files <- imData[i,2]
-        alg <- imData[i,3]
-        files <- unlist(strsplit(files, "[:]"))
-
-        ext1 <- unlist(strsplit(files[2], "[.]"))[2]
-        ext2 <- unlist(strsplit(files[4], "[.]"))[2]
-        if (ext1 == "jpg" || ext1 == "JPEG" || ext1 == "jpeg" || ext1 == "JPG") {
-            img1 <- readJPEG(files[2])
-        } else if (ext1 == "PNG" || ext1 == "png") {
-            img1 <- readPNG(files[2])
-        } else if (ext1 == "TIFF" || ext1 == "tiff" || ext1 == "TIF" || ext1 == "tif") {
-            img1 <- readTIFF(files[2])
-        } else {
-            next
-        }
-        if (ext2 == "jpg" || ext2 == "JPEG" || ext2 == "jpeg" || ext2 == "JPG") {
-            img2 <- readJPEG(files[4])
-        } else if (ext2 == "PNG" || ext2 == "png") {
-            img2 <- readPNG(files[4])
-        } else if (ext2 == "TIFF" || ext2 == "tiff" || ext2 == "TIF" || ext2 == "tif") {
-            img2 <- readTIFF(files[4])
-        } else {
-            next
-        }
-        name1 <- files[1]
-        name2 <- files[3]
+    printImages <- function(images, label) {
+        for (i in 1:nrow(images)) {
+            score <- images[i,1]
+            files <- images[i,2]
+            alg <- images[i,3]
+            files <- unlist(strsplit(files, "[:]"))
+
+            ext1 <- unlist(strsplit(files[2], "[.]"))[2]
+            ext2 <- unlist(strsplit(files[4], "[.]"))[2]
+            if (ext1 == "jpg" || ext1 == "JPEG" || ext1 == "jpeg" || ext1 == "JPG") {
+                img1 <- readJPEG(files[2])
+            } else if (ext1 == "PNG" || ext1 == "png") {
+                img1 <- readPNG(files[2])
+            } else if (ext1 == "TIFF" || ext1 == "tiff" || ext1 == "TIF" || ext1 == "tif") {
+                img1 <- readTIFF(files[2])
+            } else {
+                next
+            }
+            if (ext2 == "jpg" || ext2 == "JPEG" || ext2 == "jpeg" || ext2 == "JPG") {
+                img2 <- readJPEG(files[4])
+            } else if (ext2 == "PNG" || ext2 == "png") {
+                img2 <- readPNG(files[4])
+            } else if (ext2 == "TIFF" || ext2 == "tiff" || ext2 == "TIF" || ext2 == "tif") {
+                img2 <- readTIFF(files[4])
+            } else {
+                next
+            }
+            name1 <- files[1]
+            name2 <- files[3]
  
-        g1 <- rasterGrob(img1, interpolate=TRUE)
-        g2 <- rasterGrob(img2, interpolate=TRUE)
+            g1 <- rasterGrob(img1, interpolate=TRUE)
+            g2 <- rasterGrob(img2, interpolate=TRUE)
  
-        plot1 <- qplot(1:10, 1:10, geom="blank") + annotation_custom(g1, xmin=-Inf, xmax=Inf, ymin=-Inf, ymax=Inf) + theme(axis.line=element_blank(), axis.text.x=element_blank(), axis.text.y=element_blank(), axis.ticks=element_blank(), panel.background=element_blank()) + labs(title=alg) + ylab(unlist(strsplit(files[2], "[/]"))[length(unlist(strsplit(files[2], "[/]")))]) + xlab(name1)
-        plot2 <- qplot(1:10, 1:10, geom="blank") + annotation_custom(g2, xmin=-Inf, xmax=Inf, ymin=-Inf, ymax=Inf) + theme(axis.line=element_blank(), axis.text.x=element_blank(), axis.text.y=element_blank(), axis.ticks=element_blank(), panel.background=element_blank()) + labs(title=paste("Impostor score =", score)) + ylab(unlist(strsplit(files[4], "[/]"))[length(unlist(strsplit(files[4], "[/]")))]) + xlab(name2)
+            plot1 <- qplot(1:10, 1:10, geom="blank") + annotation_custom(g1, xmin=-Inf, xmax=Inf, ymin=-Inf, ymax=Inf) + theme(axis.line=element_blank(), axis.text.x=element_blank(), axis.text.y=element_blank(), axis.ticks=element_blank(), panel.background=element_blank()) + labs(title=alg) + ylab(unlist(strsplit(files[2], "[/]"))[length(unlist(strsplit(files[2], "[/]")))]) + xlab(name1)
+            plot2 <- qplot(1:10, 1:10, geom="blank") + annotation_custom(g2, xmin=-Inf, xmax=Inf, ymin=-Inf, ymax=Inf) + theme(axis.line=element_blank(), axis.text.x=element_blank(), axis.text.y=element_blank(), axis.ticks=element_blank(), panel.background=element_blank()) + labs(title=paste(label, " score =", score)) + ylab(unlist(strsplit(files[4], "[/]"))[length(unlist(strsplit(files[4], "[/]")))]) + xlab(name2)
  
-        multiplot(plot1, plot2, cols=2)
+            multiplot(plot1, plot2, cols=2)
+        }
     }
+    printImages(imData, "Impostor")
+    printImages(gmData, "Genuine")
+}
  
-    # Print genuine matches below the EER
-    for (i in 1:nrow(gmData)) {
-        score <- gmData[i,1]
-        files <- gmData[i,2]
-        alg <- gmData[i,3]
-        files <- unlist(strsplit(files, "[:]"))
-
-        ext1 <- unlist(strsplit(files[2], "[.]"))[2]
-        ext2 <- unlist(strsplit(files[4], "[.]"))[2]
-        if (ext1 == "jpg" || ext1 == "JPEG" || ext1 == "jpeg" || ext1 == "JPG") {
-            img1 <- readJPEG(files[2])
-        } else if (ext1 == "PNG" || ext1 == "png") {
-            img1 <- readPNG(files[2])
-        } else if (ext1 == "TIFF" || ext1 == "tiff" || ext1 == "TIF" || ext1 == "tif") {
-            img1 <- readTIFF(files[2])
-        } else {
-            next
+plotLandmarkSamples <- function(samples=NULL, expData=NULL, extData=NULL) {
+    print(plotImage(samples[[1]], "Sample Landmarks", sprintf("Total Landmarks: %s", samples[[1]]$value)))
+    if (nrow(EXT) != 0 && nrow(EXP)) {
+        for (j in 1:length(algs)) {
+            truthSample <- readData(EXT[EXT$. == algs[[j]],])
+            predictedSample <- readData(EXP[EXP$. == algs[[j]],])
+            for (i in 1:length(predictedSample)) {
+                multiplot(plotImage(predictedSample[[i]], sprintf("%s\nPredicted Landmarks", algs[[j]]), sprintf("Average Landmark Error: %.3f", predictedSample[[i]]$value)), plotImage(truthSample[[i]], "Ground Truth\nLandmarks", ""), cols=2)
+            }
         }
-        if (ext2 == "jpg" || ext2 == "JPEG" || ext2 == "jpeg" || ext2 == "JPG") {
-            img2 <- readJPEG(files[4])
-        } else if (ext2 == "PNG" || ext2 == "png") {
-            img2 <- readPNG(files[4])
-        } else if (ext2 == "TIFF" || ext2 == "tiff" || ext2 == "TIF" || ext2 == "tif") {
-            img2 <- readTIFF(files[4])
-        } else {
+    }
+}
+
+readImageData <- function(data) {
+    examples <- list()
+    for (i in 1:nrow(data)) {
+        path <- data[i,1]
+        value <- data[i,2]
+        file <- unlist(strsplit(path, "[.]"))[1]
+        ext <- unlist(strsplit(path, "[.]"))[2]
+        if (ext == "jpg" || ext == "JPEG" || ext == "jpeg" || ext == "JPG") {
+            img <- readJPEG(path)
+        } else if (ext == "PNG" || ext == "png") {
+            img <- readPNG(path)
+        } else if (ext == "TIFF" || ext == "tiff" || ext == "TIF" || ext == "tif") {
+            img <- readTIFF(path)
+        }else {
             next
         }
-        name1 <- files[1]
-        name2 <- files[3]
-
-        g1 <- rasterGrob(img1, interpolate=TRUE)
-        g2 <- rasterGrob(img2, interpolate=TRUE)
-
-        plot1 <- qplot(1:10, 1:10, geom="blank") + annotation_custom(g1, xmin=-Inf, xmax=Inf, ymin=-Inf, ymax=Inf) + theme(axis.line=element_blank(), axis.text.x=element_blank(), axis.text.y=element_blank(), axis.ticks=element_blank(), panel.background=element_blank()) + labs(title=alg) + ylab(unlist(strsplit(files[2], "[/]"))[length(unlist(strsplit(files[2], "[/]")))]) + xlab(name1)
-        plot2 <- qplot(1:10, 1:10, geom="blank") + annotation_custom(g2, xmin=-Inf, xmax=Inf, ymin=-Inf, ymax=Inf) + theme(axis.line=element_blank(), axis.text.x=element_blank(), axis.text.y=element_blank(), axis.ticks=element_blank(), panel.background=element_blank()) + labs(title=paste("Genuine score =", score)) + ylab(unlist(strsplit(files[4], "[/]"))[length(unlist(strsplit(files[4], "[/]")))]) + xlab(name2)
-
-        multiplot(plot1, plot2, cols=2)
+        example <- list(file = file, value = value, image = img)
+        examples[[i]] <- example
     }
+    return(examples)
+}
  
+plotImage <- function(image, title=NULL, label=NULL) {
+    p <- qplot(1:10, 1:10, geom="blank") + annotation_custom(rasterGrob(image$image), xmin=-Inf, xmax=Inf, ymin=-Inf, ymax=Inf)
+    p <- p + theme(axis.line=element_blank(), axis.title.y=element_blank(), axis.text.x=element_blank(), axis.text.y=element_blank(), line=element_blank(), axis.ticks=element_blank(), panel.background=element_blank())
+    p <- p + labs(title=title) + xlab(label)
+    return(p)
 }
+