Export spatial results in GIS-ready formats
Examples
# \donttest{
# Use tempdir to avoid polluting check directory
results <- run_builtin_analysis(scale = "county", year = 2016)
#>
#> ======================================================================
#> BATCH MANURESHED ANALYSIS
#> ======================================================================
#> Year: 2016
#> Scale: county
#> Nutrients: nitrogen, phosphorus
#> ----------------------------------------------------------------------
#> Checking data availability...
#> Data availability confirmed
#> Available scales:county, huc8, huc2
#> Available years forcounty:1987-2016
#> Built-in WWTP data:Available (2007-2016)
#>
#> Loading built-in NuGIS data...
#> Using cached version of nugis_county_data
#> Loaded NuGIS county data for year 2016
#> Number of spatial units: 3058
#> Loading built-in spatial boundaries...
#> Using cached version of county_boundaries
#> Loaded county boundaries
#> Number of spatial units: 3112
#> Calculating cropland threshold...
#> Cropland threshold:1235.53acres
#>
#> Processing agricultural classifications...
#> Starting complete agricultural classification for county scale...
#> Processed NuGIS data for county scale:
#> Spatial units: 3058
#> Converted P2O5 to P using factor: 0.436
#> Nitrogen classification summary:
#> Excluded: 158 units
#> Sink_Deficit: 2450 units
#> Sink_Fertilizer: 214 units
#> Source: 83 units
#> Within_County: 153 units
#> Phosphorus classification summary:
#> Excluded: 158 units
#> Sink_Deficit: 1891 units
#> Sink_Fertilizer: 372 units
#> Source: 317 units
#> Within_County: 320 units
#> Agricultural classification complete!
#> Applied threshold: 1235.53 acres
#> Agricultural classification complete
#> Spatial units processed:3112
#> Nitrogen classes:Excluded ( 212 ), Sink_Deficit ( 2450 ), Sink_Fertilizer ( 214 ), Source ( 83 ), Within_County ( 153 )
#> Phosphorus classes:Excluded ( 212 ), Sink_Deficit ( 1891 ), Sink_Fertilizer ( 372 ), Source ( 317 ), Within_County ( 320 )
#>
#> Processing WWTP data...
#> Nutrients:nitrogen, phosphorus
#> WWTP year:2016
#> Load units:kg
#> Data source:Built-in (2016)
#> Loading built-in nitrogen WWTP data for2016...
#> Using cached version of wwtp_nitrogen_combined
#> Loaded WWTP nitrogen data for year 2016
#> Number of facilities: 20846
#> Loading built-in phosphorus WWTP data for2016...
#> Using cached version of wwtp_phosphorus_combined
#> Loaded WWTP phosphorus data for year 2016
#> Number of facilities: 10148
#> Processing nitrogen WWTP facilities...
#> Filtered for positive nitrogen loads:
#> Original: 20846 facilities
#> With positive loads: 20846 facilities
#> WWTP nitrogen source classification:
#> Minor Source: 17801 facilities
#> Small Source: 1799 facilities
#> Medium Source: 729 facilities
#> Large Source: 417 facilities
#> Very Large Source: 100 facilities
#> Created spatial WWTP data with 20846 facilities
#> Looking for boundary ID column: 'FIPS'
#> Available columns in boundaries: FIPS, County, State_Name, geometry
#> Aggregating WWTP nitrogen loads by spatial boundaries...
#> Aggregation complete:
#> WWTP facilities: 20846
#> Spatial units with facilities: 2581
#> Total nitrogen load: 582738.6 tons/year
#> Processing phosphorus WWTP facilities...
#> Filtered for positive phosphorus loads:
#> Original: 10148 facilities
#> With positive loads: 10148 facilities
#> WWTP phosphorus source classification:
#> Minor Source: 6846 facilities
#> Small Source: 1839 facilities
#> Medium Source: 795 facilities
#> Large Source: 561 facilities
#> Very Large Source: 107 facilities
#> Created spatial WWTP data with 10148 facilities
#> Looking for boundary ID column: 'FIPS'
#> Available columns in boundaries: FIPS, County, State_Name, geometry
#> Aggregating WWTP phosphorus loads by spatial boundaries...
#> Aggregation complete:
#> WWTP facilities: 10148
#> Spatial units with facilities: 1966
#> Total phosphorus load: 177064.66 tons/year
#> WWTP data processing complete
#> nitrogen:20846facilities in2581spatial units
#> phosphorus:10148facilities in1966spatial units
#>
#> Integrating WWTP and agricultural data...
#> Integrating WWTP nitrogen data with agricultural classifications...
#> Using agricultural ID column: FIPS
#> Combined nitrogen classification summary:
#> Excluded: 212 units
#> Sink_Deficit: 2358 units
#> Sink_Fertilizer: 185 units
#> Source: 157 units
#> Within_County: 200 units
#> Integrating WWTP phosphorus data with agricultural classifications...
#> Using agricultural ID column: FIPS
#> Combined phosphorus classification summary:
#> Excluded: 212 units
#> Sink_Deficit: 1819 units
#> Sink_Fertilizer: 344 units
#> Source: 389 units
#> Within_County: 348 units
#> Integration complete
#>
#> Saving results...
#> Saved spatial data to: /tmp/RtmpJm2r1A/county_agricultural_2016.rds
#> File size: 1.04 MB
#> Rows: 3112, Columns: 17
#> Geometry type: POLYGON
#> CRS: EPSG:5070
#> Saved spatial data to: /tmp/RtmpJm2r1A/county_nitrogen_integrated_2016.rds
#> File size: 1.1 MB
#> Rows: 3112, Columns: 23
#> Geometry type: POLYGON
#> CRS: EPSG:5070
#> Saved centroid data to: /tmp/RtmpJm2r1A/county_nitrogen_centroids_2016.csv
#> File size: 760.52 KB
#> Rows: 3112, Columns: 24
#> Longitude range: [-124.158, -67.637]
#> Latitude range: [25.49, 48.826]
#> Saved spatial data to: /tmp/RtmpJm2r1A/county_phosphorus_integrated_2016.rds
#> File size: 1.09 MB
#> Rows: 3112, Columns: 23
#> Geometry type: POLYGON
#> CRS: EPSG:5070
#> Saved centroid data to: /tmp/RtmpJm2r1A/county_phosphorus_centroids_2016.csv
#> File size: 736.84 KB
#> Rows: 3112, Columns: 24
#> Longitude range: [-124.158, -67.637]
#> Latitude range: [25.49, 48.826]
#> Saved analysis summary to: /tmp/RtmpJm2r1A/analysis_summary_2016.rds
#> Format: RDS
#> File size: 0.76 KB
#> Results saved to:/tmp/RtmpJm2r1A
#> Files created:6
#>
#> ======================================================================
#> ANALYSIS COMPLETE
#> ======================================================================
#> Processing time:0.07minutes
#> Scale:county
#> Year:2016
#> Nutrients analyzed:nitrogen, phosphorus
#> Spatial units:3112
#> WWTP facilities:30994
#> Output directory:/tmp/RtmpJm2r1A
#> ======================================================================
output_dir <- file.path(tempdir(), "gis_outputs")
gis_files <- export_for_gis(results, output_dir)
#> Warning: Field names abbreviated for ESRI Shapefile driver
#> Warning: GDAL Error 1: /tmp/RtmpJm2r1A/gis_outputs/agricultural_results.shp does not appear to be a file or directory.
#> Created shapefile: /tmp/RtmpJm2r1A/gis_outputs/agricultural_results.shp
#> Created GeoJSON: /tmp/RtmpJm2r1A/gis_outputs/agricultural_results.geojson
#> Warning: Field names abbreviated for ESRI Shapefile driver
#> Warning: GDAL Error 1: /tmp/RtmpJm2r1A/gis_outputs/nitrogen_integrated.shp does not appear to be a file or directory.
#> Warning: Field names abbreviated for ESRI Shapefile driver
#> Warning: GDAL Error 1: /tmp/RtmpJm2r1A/gis_outputs/phosphorus_integrated.shp does not appear to be a file or directory.
# }