| #!/usr/bin/env python3 | |
| # -*- coding: utf-8 -*- | |
| """ | |
| Created on Fri May 12 19:09:58 2017 | |
| @author: aa | |
| """ | |
| import os | |
| import shutil |
| # Run it the QGIS Python command window: | |
| qgis.utils.iface.mapNavToolToolBar().move(150,40) |
I periodically benchmark my computer systems in terms of their ability to carry out data analysis and related tasks in R. To do this I use the R benchmarkme package and I persist each set of results to a file (.rds to retain R class attribute) for posterity. I also log the run here at this very Gist.
I use the run_mark.R function, which I wrote based loosely on suggestions in the benchmarkme vignette, to do the benchmark tests and afterparty.
| height: 780 |
| # Example use of the openxlsx package | |
| require(openxlsx) # For writing output in excel for collaborators | |
| df1 <- data.frame(aa = letters[3:12], bb = 21:30, stringsAsFactors = FALSE) | |
| df2 <- data.frame(cc = rnorm(20, 5, 1), dd = runif(20, 1, 5)) | |
| ## Write the data to a new excel file | |
| wb <- openxlsx::createWorkbook() | |
| addWorksheet(wb, 'Wksheet name') |
S4 is a certain specification of R object class (thus, it is useful in object oriented programming). You define an S4 class based on the data you want it to hold, and you can then create many unique copies of that S4 class that hold different data but all adhere to the same specification and can take advantage of the same functions.
The concepts of class and method are crucial to understanding object-oriented programming:
| srid,auth_name,auth_srid,srtext,proj4text | |
| 2236,EPSG,2236,PROJCS["NAD83 / Florida East (ftUS)",GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6269"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4269"]],UNIT["US survey foot",0.3048006096012192,AUTHORITY["EPSG","9003"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",24.33333333333333],PARAMETER["central_meridian",-81],PARAMETER["scale_factor",0.999941177],PARAMETER["false_easting",656166.667],PARAMETER["false_northing",0],AUTHORITY["EPSG","2236"],AXIS["X",EAST],AXIS["Y",NORTH]],+proj=tmerc +lat_0=24.33333333333333 +lon_0=-81 +k=0.999941177 +x_0=200000.0001016002 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs | |
| 2237,EPSG,2237,PROJCS["NAD83 / Florida West (ftUS)",GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS 1980",6378137,298. |
| # Allen Aven November 2015 | |
| # Function to thin a point input layer according to an input raster | |
| require(sp) | |
| require(rgdal) | |
| require(raster) | |
| require(dismo) | |
| pt_thin <- function(infile, grid=grid, n_points=5) { |