Skip to content

Instantly share code, notes, and snippets.

View oscarperpinan's full-sized avatar

Oscar Perpiñán Lamigueiro oscarperpinan

238 followers · 7 following
View GitHub Profile
@oscarperpinan
oscarperpinan / texPart.tex
Last active August 29, 2015 13:57
parent latex file of orgChapter
\part{My First Part}
\label{part:part1}
Lorem ipsum dolor sit amet, consectetur adipiscing elit. Aliquam nec
lacus nec dui semper aliquet sit amet in leo. Integer volutpat eu diam
sed iaculis.
\chapter{My First Chapter}
\label{cha:cha1}
@oscarperpinan
oscarperpinan / orgChapter.org
Last active August 29, 2015 13:57
orgmode example for writing a chapter

Code blocks

<sec:sectionA>

Code is enclosed within #+begin_src/#+end_src blocks.

x <- 1:10
@oscarperpinan
oscarperpinan / meteo.R
Created February 2, 2014 16:04
Using meteo package
## devtools::install_github('meteo', 'oscarperpinan')
library(meteo)
Sys.setenv(TZ='UTC')
st <- data.frame(name=c('Almeria','Granada','Huelva','Malaga','Caceres'),
elev=c(42, 702, 38, 29, 448))
coordinates(st) <- cbind(c(-2.46, -3.60, -6.94, -4.42, -6.37),
@oscarperpinan
oscarperpinan / stamenPoly.R
Last active December 28, 2015 09:59
Variation of https://gist.github.com/oscarperpinan/5451682 using SpatialPolygons
library(sp)
library(ggmap)
## latticeExtra must be loaded after ggmap because both ggplot2 and
## latticeExtra define a 'layer' function. We need the definition from
## latticeExtra.
library(latticeExtra)
## We only need maptools to get an example
library(maptools)
SIDS <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1],
IDvar="FIPSNO",
@oscarperpinan
oscarperpinan / californiaOpenPV.geojson
Created September 20, 2013 08:57
OpenPV systems https://openpv.nrel.gov/search?&state=CA&minsize=100&maxsize=30600&pagenum=1&nppage=2
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@oscarperpinan
oscarperpinan / siar.R
Created September 20, 2013 08:09
Code to create siar.geojson
library(sp)
library(rgdal)
setwd(tempdir())
download.file('https://raw.github.com/oscarperpinan/solar/gh-pages/data/SIAR.csv', 'siar.csv', method='wget')
siar <- read.csv('siar.csv')
summary(siar)
siarSP <- SpatialPointsDataFrame(siar[,c(6, 7)], siar[,-c(6,7)])
writeOGR(siarSP, 'siar.geojson', 'siarSP', driver='GeoJSON')
@oscarperpinan
oscarperpinan / siar.geojson
Created September 20, 2013 08:05
SIAR stations
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@oscarperpinan
oscarperpinan / solar2Local.R
Last active December 19, 2015 05:19
library(solaR)
## Download data from NREL (see help(calcG0))
NRELurl <- 'http://goo.gl/fFEBN'
dat <- read.table(NRELurl, header=TRUE, sep=',')
names(dat) <- c('date', 'hour', 'G0', 'B', 'D0', 'Ta')
##B is direct normal. We need direct horizontal.
dat$B0 <- dat$G0-dat$D0
@oscarperpinan
oscarperpinan / solaR36.R
Created June 13, 2013 10:06
library(solaR)
## Monthly averages of daily values (Wh/m2)
G0dm <- c(2.766,3.491,4.494,5.912,6.989,7.742,7.919,7.027,5.369,3.562,2.814,2.179)*1000;
Ta <- c(10, 14.1, 15.6, 17.2, 19.3, 21.2, 28.4, 29.9, 24.3, 18.2, 17.2, 15.2)
g0 <- calcG0(lat=37.2, modeRad='aguiar', dataRad=G0dm)
xyplot(g0)
@oscarperpinan
oscarperpinan / hillShade.R
Last active July 15, 2020 21:46
library(raster)
library(rasterVis)
library(maptools) ## for readShapeLines
library(colorspace) ## for terrain_hcl
## Get data from GADM and DIVA-GIS
setwd(tempdir())
download.file('http://www.gadm.org/data/shp/ETH_adm.zip', 'ETH_adm.zip')
unzip('ETH_adm.zip', exdir='.')
download.file('http://www.diva-gis.org/data/msk_alt/ETH_msk_alt.zip', 'ETH_msk_alt.zip')