Parallel Raster Extraction

This is my job submission script for the raster extraction. #!/bin/bash -l #PBS -l walltime=48:00:00,nodes=1:ppn=24,mem=61gb #PBS -m abe #PBS -M myaddress@email.com module load gdal/2.0.1 module load proj module load R module load ompi module load udunits module load geos export OMP_NUM_THREADS=1 mpirun -v -hostfile $PBS_NODEFILE -np 1 R –slave CMD BATCH ./extract.R /dev/stdout This is […]

Continue reading..Parallel Raster Extraction

Weekend GIS Fun: Part 2

See part 1 for the relevant shape file and description of what’s going on! In the last post, we identified the map codes of interest to subset the shape file. Here’s some code for working with it in R. library(rgdal) bedrock <- readOGR(dsn=”~/Arrowhead_South_Bedrock_Geology”,”Arrowhead_South_Bedrock_Geology”) basalt = subset(bedrock,map_label==”Mnd” | map_label==”Mna”) writeOGR(basalt,dsn=”.”,layer=”basalt”,driver=”ESRI Shapefile”) So now we have a […]

Continue reading..Weekend GIS Fun: Part 2