xuanlongma

## TRMM 3B43 Monthly precipitation does NOT need flip (which 3B42 need)
r.trmm <- raster(nrows = 400, ncols = 1440, xmn = 0, xmx = 360, ymn =
                 -50, ymx = 50, crs="+proj=longlat +datum=WGS84")

r.trmm[] <- readBin(files.input, 'double', n = 576000, size  = 4,
           endian = 'big')

## TRMM 3B42
## r.trmm <- flip(rotate(r.trmm), 'y')

xuanlongma

pandoc --bibliography=bibtex.bib --mathjax --parse-raw -s inpud.md -o output.html

xuanlongma

If multiple input files are given, pandoc will concatenate them all (with blank lines between them) before parsing. -- from Pandoc website

Pandoc command: pandoc -s input1.md input2.md input3.md -o output.html

xuanlongma

## a tiff file 
input <- "modis.evi.tif"
output <- "modis.evi.reprojected.tif"

## GDAL command line for reprojection
proj.cmd.warp <- 'gdalwarp -t_srs \'+proj=latlong +datum=WGS84\' -r near -overwrite'

## Invoke the system GDAL command
system(command = paste(proj.cmd.warp, input, output, sep = ' '))

xuanlongma

# Script for question posted on Stack Overflow
#	Load relevant libraries
library(ggplot2)
library(raster)
library(gridExtra)


vplayout <- function(x, y) {
	viewport(layout.pos.row = x, layout.pos.col = y)
	}

xuanlongma

require(ncdf)
require(raster)

## Input: a netCDF file
file.nc <- 'rain.nc'

## Output: a GeoTIFF file
file.tiff <- 'rain.tiff'

## Import netCDF

xuanlongma

strptime('2008201', format = '%Y%j')
## or
as.Date('2008201', format = '%Y%j')