WarehouseMap/BuiltWH_intersect.R at main · RadicalResearchLLC/WarehouseMap · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
'%ni%' <- Negate('%in%') ## not in operator
library(leaflet)


recently_builtWH <- final_parcels |>
  filter(year_built > 2016)

earlier_projects <- plannedWarehouses #|>
##FIXME - recv'd date is most recent date and may indicate projects with recent updates
  #filter(recvd_date <= mdy('1/1/2024')) #|>
  #st_centroid()

check_wh_list <- st_intersection(recently_builtWH, earlier_projects) #|>
  #st_set_geometry(value = NULL)

check_wh_list2 <- st_intersection(recently_builtWH, earlier_projects2)

check_wh_list3 <- check_wh_list2 |>
  st_set_geometry(value = NULL)  |>
  select(sch_number)

area_of_intersection <- check_wh_list |>
  st_set_geometry(value = NULL) |>
  group_by(sch_number) |>
  summarize(count = n(), areaBuilt = sum(shape_area2),
            totalParcel = median(parcel_area), .groups = 'drop') |>
  mutate(built_area = round(areaBuilt/totalParcel, 2))

belowHalf <- area_of_intersection |>
  filter(built_area < 0.4)

alreadyBuilt <- area_of_intersection |>
  filter(built_area >= 0.4)


#recently_builtIntersect <- recently_builtWH |>
#  st_filter(check_wh_list)

#check_wh_list2 <- st_intersection(plannedWarehouses, recently_builtIntersect)

leaflet() |>
  setView(lat = 34, lng = -117.60, zoom = 9) |>
  addProviderTiles(providers$CartoDB.Positron, group = 'Basemap')  |>
  addProviderTiles(providers$Esri.WorldImagery, group = 'Imagery') |>
  addLayersControl(baseGroups = c('Basemap', 'Imagery')) |>
  addPolygons(data = recently_builtWH,
              color = 'blue',
              weight = 2,
              fillOpacity = 0.2,
              label = ~year_built
  ) |>
  addPolygons(data = earlier_projects,
              color = 'orange',
              fillOpacity = 0.2,
              label = ~sch_number) |>
  addPolygons(data = check_wh_list2,
              color = 'red',
              weight = 1,
              fillOpacity = 0.2,
              label = ~sch_number) #|>
 # addMarkers(data = check_wh_list)

## check for non-centroid parcels
leaflet() |>
  setView(lat = 34, lng = -117.60, zoom = 9) |>
  addProviderTiles(providers$CartoDB.Positron, group = 'Basemap')  |>
  addProviderTiles(providers$Esri.WorldImagery, group = 'Imagery') |>
  addLayersControl(baseGroups = c('Basemap', 'Imagery')) |>
  addPolygons(data = earlier_projects2,
              color = 'orange',
              fillOpacity = 0.2,
              label = ~sch_number) |>
  addPolygons(data = difference,
              color = 'red',
              weight = 1,
              fillOpacity = 0.2,
              label = ~sch_number) |>
  addMarkers(data = st_centroid(difference))

OddCentroidsBuilt <- check_wh_list2 |>
  filter(sch_number == '2020079023') |>
  select(project, apn, sch_number, ceqa_url) |>
  st_set_geometry(value = NULL)

approved_built_since2022 <- check_wh_list |>
  st_set_geometry(value = NULL) |>
  filter(sch_number %ni% c('2020059028', '2022020461')) |>
  select(project, apn, sch_number, ceqa_url) |> # less than 50% buildout projects
  bind_rows(OddCentroidsBuilt)

## Anti-join the built warehouses

planned_notBuilt <- plannedWarehouses |>
  anti_join(alreadyBuilt)

## Add filter list for specific projects that are built out more than 50%
## FIXME

planned_tidy <- planned_notBuilt |>
  mutate(shape_area =  round(parcel_area, -3),
         class = stage_pending_approved,
         year_chr = 'future',
         year_built = 2025,
         type = 'warehouse',
         row = row_number()) |>
  select(-parcel_area, -stage_pending_approved) |>
  mutate(county = str_c(county, ' County'))

leaflet() |>
  setView(lat = 34, lng = -117.60, zoom = 9) |>
  addProviderTiles(providers$CartoDB.Positron, group = 'Basemap')  |>
  addProviderTiles(providers$Esri.WorldImagery, group = 'Imagery') |>
  addLayersControl(baseGroups = c('Basemap', 'Imagery')) |>
  addPolygons(data = recently_builtWH,
              color = 'blue',
              fillOpacity = 0.2,
              weight = 1) |>
  addPolygons(data = planned_tidy,
              color = 'orange',
              weight = 1,
              fillOpacity = 0.2,
              label = ~sch_number)

rm(ls = planned_notBuilt, approved_built_since2022,
   check_wh_list, check_wh_list2, earlier_projects,
   recently_builtWH, recently_builtIntersect,
   earlier_projects2, difference, check_wh_list3)

rm(ls = alreadyBuilt, area_of_intersection, belowHalf)
rm(ls = OddCentroidsBuilt)
#leaflet() |>
  #addPolygons(data = planned_tidy,
    #          weight = 2,
   #           color = 'red') |>
  #addPolygons(data = recently_builtWH,
 #             weight = 2, color = 'blue') |>
#  addTiles()