Tag: sf
-
California cities from CDTFA
At least based on the dates listed on https://lab.data.ca.gov/dataset/california-city-boundaries-and-identifiers — still says “03/24/25” — I’m guessing the city boundaries haven’t been updated yet. For a current project, I want the best boundaries available so let’s work directly with the source pointed out by the CIO, GIS: https://gis.data.ca.gov/datasets/CDTFA::city-and-county-boundary-line-changes/explore?layer=0&location=33.924386%2C-118.008294%2C11.67. > library(sf)> d <- read_sf(“~/Downloads/City_and_County_Boundary_Line_Changes_217120214122415134.gpkg”) It turns out…
-
A hole in Los Altos
Looking at the Census page for Los Altos, something funny caught my eye: there’s a hole in the map for Los Altos: https://data.census.gov/profile/Los_Altos_city,_California?g=160XX00US0643280 Returning to the data from https://lab.data.ca.gov/dataset/california-city-boundaries-and-identifiers, we see the same thing: > library(ggplot2)> library(sf)> cities <- read_sf(“~/Downloads/California_Cities_and_Identifiers_Blue_Version_view_-6943741225906831761.gpkg”)> losAltos <- subset(cities, CDTFA_CITY == “Los Altos”)> ggplot(losAltos) + geom_sf() What’s in that hole? >…
-
What Wittier Whittier
I was really impressed to get an almost immediate reply from the CIO, GIS to my note about Whittier appearing to some extra stuff in the geometries provided here: https://lab.data.ca.gov/dataset/california-city-boundaries-and-identifiers. They pointed to the California Department of Tax and Fee Administration as the original source of the geometries, which has been updated in the meantime:…
-
California cities (continued)
Let’s get back roughly to where we were with the R session from California cities, reading in a copy of the data pulled from Wikipedia. > library(dplyr)> library(fillpattern)> library(sf)> library(ggplot2)> cities <- read_sf(“~/Downloads/California_Cities_and_Identifiers_Blue_Version_view_-6943741225906831761.gpkg”)> stats <- read.csv(url(“https://raw.githubusercontent.com/fadend/ca_cities_data/refs/heads/main/ca_cities.csv”), stringsAsFactors=FALSE) Oops. There’s some disagreement on names: > setdiff(stats$city_name, cities$city_name)[1] “Angels Camp” “California City” “City of Industry”> setdiff(cities$city_name, stats$city_name)[1]…