From 1bacd2467c78c721b7b16941df245c1ca49f6d4b Mon Sep 17 00:00:00 2001 From: Alex Date: Fri, 3 Apr 2026 16:41:51 -0600 Subject: [PATCH] Updated to include water rights --- 1_Data_Proc.r | 41 +++++++++++++++++++++++++++++++++++++++-- 1 file changed, 39 insertions(+), 2 deletions(-) diff --git a/1_Data_Proc.r b/1_Data_Proc.r index 0cd5469..c693ae2 100755 --- a/1_Data_Proc.r +++ b/1_Data_Proc.r @@ -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) -CROP_LEGAL %>% select(farm_unit_2) %>% group_by(farm_unit_2) %>% summarize(n()) +#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 %>% 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