Ryan Hill ryan-hill
ryan-hill
/ float2int.py
Created
October 14, 2015 00:35
— forked from bryanluman/float2int.py
Converting a large DEM with rasterio
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| import numpy as np | |
| import rasterio | |
| """ | |
| 2014-02-13 | |
| Bryan Luman | |
| Use it however you like at your own risk | |
| Problem: | |
| I have a huge DEM converted from LiDAR LAS points. I'd like to make it slightly |
ryan-hill
/ Lookup.R
Created
October 1, 2015 19:10
— forked from mhweber/Lookup.R
This code takes a lookup table and applies it to a data frame, updating only values for records that occur in the lookup table using indexing and match
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| # Create data frame 1 | |
| x = c("ID1","ID2","ID3","ID4","ID5") | |
| y = c("C1","C2","C3","C4","C5") | |
| d1 = data.frame("SiteID" = x, "Value" = y) | |
| d1 | |
| # Create lookup table | |
| x = c("ID2","ID5") | |
| y = c("C5","C2") | |
| lookup = data.frame("SiteID" = x, "Value" = y) | |
| lookup |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| import pysal as ps | |
| import pandas as pd | |
| ''' | |
| Arguments | |
| --------- | |
| dbfile : DBF file - Input to be imported | |
| upper : Condition - If true, make column heads upper case | |
| ''' | |
| def dbf2DF(dbfile, upper=True): #Reads in DBF files and returns Pandas DF | |
| db = ps.open(dbfile) #Pysal to open DBF |
ryan-hill
/ allPairwisePredictorCombos.R
Last active
September 18, 2015 20:33
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #predData is dataframe of predictor variable where 1st column is a unique site ID | |
| vect = names(predData)[2:length(predData)] | |
| combos = cbind(combn(vect, 2)[1,], combn(vect, 2)[2,]) | |
| #Useful for comparing pairwise combo models | |
| combos = data.frame(combos); combos$AUC <- 0.0 | |
| names(combos) <- c('Pred1','Pred2','AUC') |
ryan-hill
/ setRasterNAto0.py
Created
September 11, 2015 19:02
Code snippet to convert NA values within conterminous USA (NLCD data as mask) to 0s
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| path = 'xyz/' | |
| mask_grid = path + 'nlcd2006.tif' | |
| arcpy.env.mask = mask_grid | |
| PCH = path + 'PCH.tif' | |
| PCH = Con(IsNull(PCH),0,PCH) | |
| PCH.save(path + 'PCH2.tif') |
NewerOlder