alex/CREP_Policy_Interaction_Analysis
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Updated to include water rights

Browse Source
This commit is contained in:
Alex 2026-04-03 16:41:51 -06:00
parent de4e1359dc
commit 1bacd2467c
1 changed files with 39 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

41
1_Data_Proc.r
View File

@ -3,10 +3,47 @@ library(fixest)
library(janitor)
###########
CROP_LEGAL <- read_csv("Data/Created/Crop_Legal_Int.csv") %>% clean_names() %>% select(-fid,-fid_2,-div,-district,-master_id,-master_id_2,-polygon_id,-shape_area,-int_area,-comments,-district_2,-path,-comments_2,-shape_leng_2,-shape_area_2,-master_id_2,-gwzone_2,-comments_2,-shape_leng_2,shape_area_2,-path)
#From hydrobase Water Rights -Net Amounts search, in Div 3 and for IRRIGATION
#https://dwr.state.co.us/Tools/WaterRights/NetAmounts?submitButton=Submit&SelectedGeoValue=waterDivisionDiv&SelectedWaterDivisionId=3&SelectedStructureId=2&SelectedUsageTypeId=1
WATER_RIGHTS <- read_csv("Data/Downloaded/CDSS_WaterRights_NetAmounts_Search_202604031353.csv") %>% clean_names()
#WATER_RIGHTS %>% filter(wdid=='2013695') %>% select(net_apex_absolute,net_apex_conditional,net_absolute,net_conditional,seasonal_limits)
WATER_RIGHTS <- WATER_RIGHTS %>% mutate(rights_absolute=ifelse(net_absolute>net_conditional,net_absolute,net_conditional),rights_apex=ifelse(net_apex_absolute>net_apex_conditional,net_apex_absolute,net_apex_conditional) ) %>% filter(decreed_units=='C')
WATER_RIGHTS <-WATER_RIGHTS %>% mutate(adjudication_date=year(as.POSIXct(adjudication_date, format="%m/%d/%Y %I:%M:%S %p")),previous_adj_date=year(as.POSIXct(previous_adj_date, format="%m/%d/%Y %I:%M:%S %p")))
WATER_RIGHTS <- WATER_RIGHTS %>% select(wdid,rights_absolute,rights_apex,adjudication_date,previous_adj_date)
WATER_RIGHTS %>% mutate(cal_year=ifelse(adjudication_date<=2002,2002,NA)) %>% filter(!is.na(cal_year)) %>% group_by(wdid,cal_year) %>% summarize(rights_absolute=sum(rights_absolute),rights_apex=sum(rights_apex)) %>% ungroup
CROP_LEGAL
PARCEL_YEAR_WDID <- CROP_LEGAL %>% select(parcel_id,cal_year,gw_id1,gw_id2,gw_id3,gw_id4,gw_id5,gw_id6,gw_id7,gw_id8,gw_id9,gw_id10,gw_id11,gw_id12,gw_id13,gw_id14,gw_id15,gw_id16,gw_id17,gw_id18,gw_id19,gw_id20) %>% pivot_longer(-c("parcel_id","cal_year"),values_to="wdid") %>% select(-name) %>% filter(!is.na(wdid))
WELL_LIST <- PARCEL_YEAR_WDID %>% pull(wdid) %>% unique
WATER_RIGHTS <- WATER_RIGHTS %>% filter(wdid %in% WELL_LIST)
WATER_RIGHTS %>% group_by(wdid) %>% summarize(rights_absolute=sum(rights_absolute),rights_apex=sum(rights_apex)) %>% filter(rights_absolute==0,rights_apex==0)
WATER_RIGHTS %>% filter(rights_absolute==0,rights_apex==0)
read_csv("https://dwr.state.co.us/Rest/GET/api/v2/waterrights/transaction/?format=csv&wdid=2005165",skip=2)
WATER_RIGHTS %>% filter(wdid=='2013695')
WATER_RIGHTS %>% mutate(
WATER_RIGHTS %>% select(wdid,rights_absolute,rights_apex) %>% filter(rights_absolute==0,rights_apex==0)
WATER_RIGHTS %>% filter(wdid=='2013695')%>% select(wdid,rights_absolute,rights_apex)
WATER_RIGHTS %>% select(rights_absolute,rights_apex)
WATER_RIGHTS
WATER_RIGHTS <- WATER_RIGHTS %>% select(wdid,rights,previous_adj_date,adjudication_date) %>% mutate(adjudication_date=year(as.POSIXct(adjudication_date, format="%m/%d/%Y %I:%M:%S %p")),previous_adj_date=year(as.POSIXct(adjudication_date, format="%m/%d/%Y %I:%M:%S %p")))
WATER_RIGHTS %>% filter(adjudication_date>2002) %>% print(n=100)
WATER_RIGHTS %>% filter(wdid==2013837) %>% select(decreed_units,adjudication_date)
WATER_RIGHTS %>% group_by(decreed_units) %>% summarize(n())
CROP_LEGAL %>% select(farm_unit_2) %>% group_by(farm_unit_2) %>% summarize(n())
CROP_LEGAL
CROP_LEGAL %>% select(farm_unit_2) %>% unique
colnames(CROP_LEGAL )
CROP_LEGAL
nrow(CROP_LEGAL)-nrow(TEST
TEST <- CROP_LEGAL %>% group_by(parcel_id) %>% filter(calc_area==max(calc_area)) %>% ungroup
TEST %>% group_by(cal_year) %>% summarize(n())
CROP_LEGAL %>% select(layer) %>% gsub(" Parcels","")
CROP_LEGAL %>% filter(name=='3-200004-0')
/home/alex/Git/Publications/CREP_Policy_Map/Data/Processed