Regional Water Planning Areas

Texas Regional Water PLanning Areas - polygons

Download data

Data for regional water planning areas polygons from https://www.twdb.texas.gov/mapping/gisdata.asp

Manually download the shapefile

Read the file in and display

fgdb <- paste0(path, "TWDB_RWPAs_2014.dbf")

# Read the feature class
temp <- sf::st_read(fgdb)

Reading layer `TWDB_RWPAs_2014' from data source 
  `/home/ajackson/Dropbox/Rprojects/Curated_Data_Files/WaterPlanningAreas/TWDB_RWPAs_2014.dbf' 
  using driver `ESRI Shapefile'
Simple feature collection with 16 features and 6 fields
Geometry type: MULTIPOLYGON
Dimension:     XY
Bounding box:  xmin: 373011.9 ymin: 412841.6 xmax: 1618159 ymax: 1594957
Projected CRS: Custom

temp <- sf::st_transform(temp, googlecrs)
# sf::st_crs(temp) <- googlecrs

summary(temp)

   REG_NAME            LETTER               AREA             PERIMETER      
 Length:16          Length:16          Min.   :5.587e+09   Min.   : 408299  
 Class :character   Class :character   1st Qu.:3.079e+10   1st Qu.:1050271  
 Mode  :character   Mode  :character   Median :3.557e+10   Median :1306886  
                                       Mean   :4.375e+10   Mean   :1307408  
                                       3rd Qu.:5.389e+10   3rd Qu.:1601777  
                                       Max.   :1.037e+11   Max.   :1976957  
    Label             Label_2                   geometry 
 Length:16          Length:16          MULTIPOLYGON :16  
 Class :character   Class :character   epsg:4326    : 0  
 Mode  :character   Mode  :character   +proj=long...: 0  
                                                         
                                                         
                                                         

leaflet(temp) %>% 
  leaflet::setView(lng = -95.3103, lat = 29.7752, zoom = 8 ) %>%   
  addTiles() %>% 
  addPolygons(weight=1) %>% 
  addStaticLabels(label=temp$Label,
                  style = list("color" = "red", "font-weight" = "bold"))

# saveRDS(temp, paste0(path, "Water_Planning_Areas.Rds"))