Zev Ross President, ZevRoss Spatial Analysis
|
|
- Chester Roberts
- 5 years ago
- Views:
Transcription
1 SPATIAL ANALYSIS IN R WITH SF AND RASTER A quick refresher on the coordinate Zev Ross President, ZevRoss Spatial Analysis reference system
2
3 Coordinate reference system A place on the earth is specified by a latitude and longitude or X/Y coordinates The coordinates are based on a mathematical model of the shape of the earth Some mathematical formulas involve a transformation from the 3-D globle to a 2-dimensional map
4 Projected vs unprojected CRS An unprojected CRS uses latitude and longitude coordinates and references the earth as a three-dimensional object A projected CRS uses X and Y coordinates as a two-dimensional representation of the earth
5 Your geographic files have a CRS - but it's not always defined Vector and raster spatial data was created based on a specific CRS Usually the spatial file has metadata that tells you the CRS Sometimes there is no metadata defining the CRS
6 Both sf and raster will read the CRS if it exists in the metadata st_crs() prints out a vector object's CRS crs() prints out a raster object's CRS
7 An example: This shapefile has a defined CRS > shape1 <- st_read("shape1.shp") > st_crs(shape1) $epsg [1] 4326 $proj4string [1] "+proj=longlat +ellps=wgs84 +no_defs" attr(,"class") [1] "crs" You can tell that this is an unprojected CRS because the definition (called proj4string) starts with "+proj=longlat" referring to longitude and latitude.
8 Define CRS with EPSG or proj4string You can use either an EPSG or proj4string to define the CRS. A CRS might have both but might only need one The EPSG code is a numeric representation of a CRS (e.g, 4326) The proj4string is a full set of parameters spelled out in a string (e.g., "+proj=longlat +ellps=wgs84 +no_defs")
9 An example: This shapefile does not have a defined CRS > shape2 <- st_read("shape2.shp") > st_crs(shape2) $epsg [1] NA $proj4string [1] NA attr(,"class") [1] "crs"
10 If the CRS is not defined Do background research to find out the CRS Then tell R what the CRS with st_crs()
11 Defining the crs with st_crs() Define with the proj4string > st_crs(shape2) <- "+proj=longlat +ellps=wgs84 +no_defs" Define with the EPSG code > st_crs(shape2) <- 4236
12 For a raster define the CRS with crs() An example where the CRS is defined > crs(singleband) # CRS arguments: # +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 # +lon_0=-96 +x_0=0 +y_0=0 +ellps=grs80 # +towgs84=0,0,0,0,0,0,0 +units=m +no_defs An example where the CRS is not defined > crs(multiband) # CRS arguments: NA
13 For a raster define the CRS with crs() Here we define it > crs(multiband) <- "+proj=utm +zone=18+datum=wgs84 +units=m +no_defs+ellps=wgs84 +towgs84=0,0,0"
14 Change the CRS with st_transform() or projectraster() Use st_transform() to change the CRS for vectors Use projectraster() to change the CRS for rasters
15 Transform the CRS for vector shapes > shape1_prj <- st_transform(shape1, crs = 32618) > shape1_prj <- st_transform(shape1, crs ="+proj=utm +zone=18 +ellps=wgs84 +datum=wgs84 +units=m +no_defs")
16 Transform CRS of one layer to match another layer > shape1_prj <- st_transform(shape1, crs = crs(singleband, astext = TRUE)) Note that you need astext = TRUE to force the crs() function from raster to output the CRS as a string
17 Transform the CRS for a raster > singleband_prj <- projectraster(singleband, + crs = "+proj=utm +zone=18 +ellps=wgs84 +datum=wgs84 +units=m +no_defs") > singleband_prj <- projectraster(singleband, crs = "+init=epsg:32618") Note that to use an EPSG with projectraster the syntax is "+init=epsg:32618"
18 SPATIAL ANALYSIS IN R WITH SF AND RASTER Let's practice!
19 SPATIAL ANALYSIS IN R WITH SF AND RASTER Slicing, dicing and simplifying your Zev Ross President, ZevRoss Spatial Analysis vectors
20 Using dplyr with spatial data > head(trees) # Simple feature collection with 6 features and 2 fields # geometry type: POINT # dimension: XY # bbox: xmin: ymin: xmax: # epsg (SRID): 4326 # proj4string: +proj=longlat +ellps=wgs84 +no_defs # Source: local data frame [6 x 3] # Groups: hood [2] # A tibble: 6 x 3 # tree_id hood geometry # <dbl> <fctr> <simple_feature> # BK09 <POINT( > # BK09 <POINT( > # BK09 <POINT( > # BK17 <POINT( > # BK17 <POINT( > # BK17 <POINT( >
21 Count trees by neighborhood Use the count() function from dplyr and arrange in descending order. > cnt <- count(trees, hood) %>% + arrange(desc(n)) > cnt # Simple feature collection with 188 features and 2 fields # geometry type: MULTIPOINT # dimension: XY # bbox: xmin: ymin: xmax: # epsg (SRID): 4326 # proj4string: +proj=longlat +ellps=wgs84 +no_defs # A tibble: 188 x 3 # hood n geometry # <fctr> <int> <simple_feature> # 1 SI <MULTIPOINT(-...> # 2 BK <MULTIPOINT(-...> # 3 SI <MULTIPOINT(-...> # 4 QN <MULTIPOINT(-...> #... with 184 more rows
22 Drop geometry by setting it to NULL # Without a pipe > tree_cnt <- count(trees, hood) > tree_cnt <- st_set_geometry(tree_cnt, NULL) # With a pipe > tree_cnt <- count(trees, hood) %>% + st_set_geometry(null) > head(tree_cnt) # A tibble: 3 x 2 # hood n # <fctr> <int> # 1 BK09 3 # 2 BK17 3 # 3 BK19 3
23 Join spatial and non-spatial data - setup > tree_sm <- tree[1:3,] > head(tree_sm) # tree_id hood geometry # QN76 POINT( # MN32 POINT( # QN70 POINT( Fake data for our example. Goal is to get the val variable into the tree data. > dat <- data.frame(hood = c("qn76", "MN32", "QN70"), + val = c(1, 2, 3))
24 Join spatial and non-spatial data with inner_join() # tree_sm is spatial, dat is non-spatial > res <- inner_join(tree_sm, dat, by = "hood") # The dataset now how the "val" variable > head(res) # tree_id hood val geometry # QN76 1 POINT( # MN32 2 POINT( # QN70 3 POINT(
25 Simplify your vectors
26 Vector data can be more detailed than needed Administrative boundaries and rivers, for example, can have far more detail than required Simplifying will speed up computations
27 The size of our example data before simplification Size in memory > library(pryr) > object_size(county) 489 kb Number of vertices # "cast" our polygons to bundles of points -- "MULTIPOINT" then count > pts <- st_cast(county$geometry, "MULTIPOINT") > sum(sapply(pts, length)) [1] 57886
28 Simplify with st_simplify() Tolerance controls simplification. Bigger numbers mean more simplification. Units are the units of the CRS. > boro_simp <- st_simplify(boro, dtolerance = 500)
29 Visually there is barely a difference Non-Simplified Simplified
30 16x smaller and 33x fewer vertices Original object size Simplified object size > library(pryr) > object_size(boro) 489 kb # Cast code left off here > sum(sapply(pts, length)) [1] > library(pryr) > object_size(boro_simp) 29.4 kb # Cast code left off here > sum(sapply(pts, length)) [1] 1770
31 SPATIAL ANALYSIS IN R WITH SF AND RASTER Let's practice!
32 SPATIAL ANALYSIS IN R WITH SF AND RASTER Converting sf objects between sp and raw Zev Ross President, ZevRoss Spatial Analysis coordinates
33 sp has had a long and useful live sp was created more than a decade ago Many package make use of sp objects
34 Points often exist in non-spatial data frames of coordinates To use spatial functionality, you need to convert a data frame of coordinates to sf objects
35 Convert sf objects to sp with as() The class for sp objects is Spatial. # Our sf object > polys <- st_read("polygons.shp") > class(polys) [1] "sf" "data.frame" # Convert to Spatial object > polys_sp <- as(polys, Class = "Spatial") > class(polys_sp) [1] "SpatialPolygonsDataFrame" attr(,"package") [1] "sp"
36 Convert from sp to sf with st_as_sf() # Convert from sp to sf object > polys_sf <- st_as_sf(polys_sp) > class(polys_sf) [1] "sf" "data.frame"
37 Coordinates to an sf points object with st_as_sf() The coords argument specifies the coordinate columns and must be in longitude, latitude or X, Y order. # Simple dataframe with coordinates > pts <- data.frame(id = 1:2, lon = c(-73, -72),lat = c(41, 42)) # Convert to an sf object > pts <- st_as_sf(pts, coords = c("lon", "lat")) > pts # Simple feature collection with 2 features and 1 field # geometry type: POINT # dimension: XY # bbox: xmin: -73 ymin: 41 xmax: -72 ymax: 42 # epsg (SRID): NA # proj4string: NA # ID geometry # 1 1 POINT (-73 41) # 2 2 POINT (-72 42)
38 You can specify the CRS with the crs argument Specify the CRS with either a proj4string or a EPSG code > # WGS 84 with a proj4string > st_as_sf(pts, coords = c("lon", "lat"), + crs = "+proj=longlat +ellps=wgs84 +datum=wgs84 +no_defs") With long/lat you can probably use EPSG = 4326 (WGS 84) # WGS 84 with an EPSG code > st_as_sf(pts, coords = c("lon", "lat"), crs = 4326)
39 Write your points to a CSV with coordinates Writing to other spatial formats like shapefiles is easy with st_write() If you use st_write() with points and specify a CSV, coordinates won't be included Use a "hidden" argument, layer_options to write coordinates > st_write(pts, "pts.csv", layer_options = "GEOMETRY=AS_XY")
40 SPATIAL ANALYSIS IN R WITH SF AND RASTER Let's practice!
41 SPATIAL ANALYSIS IN R WITH SF AND RASTER Manipulating raster layers Zev Ross President, ZevRoss Spatial Analysis
42 Two key functions Reduce raster resolution with aggregate() Reclassify values with reclassify()
43 Example: Elevation data for area around Ithaca, NY > elevation <- raster("elevation.tif") > ncell(elevation) [1] > res(elevation) [1] > file.size("elevation.tif") [1] # ~ 25 megabytes > plot(elevation)
44 Reduce raster resolution with aggregate() Use the fact argument to specify the factor of aggregation > dem_low <- aggregate(dem, fact = 20) Use fun to specify the function to do aggregation > dem_low <- aggregate(dem, fact = 20, fun = mean)
45 Significant reduction in file size > dem_low <- aggregate(dem, fact = 20, fun = mean) > ncell(dem_low) [1] > res(dem_low) [1] > writeraster(dem_low, "elevation-small.tif") > file.size("elevation-small.tif") [1] # 0.11 megabytes
46 Aggregated raster
47 Reclassify raster values with reclassify() Here we use a 2-column matrix to change values of 5 to 100 > new_vals <- cbind(5, 100) > new_rast <- reclassify(old_raster, rcl = new_vals) Here we use a 3-column matrix to change values between 1 and 3 to NA > new_vals <- cbind(1, 3, NA) > new_rast <- reclassify(old_raster, rcl = new_vals)
48 Use a multi-row matrix to change many values at once Reclassify values into three groups (example from the help) # Values between 0 and 0.25 become 1 and so on > m <- c(0.00, 0.25, 1, 0.25, 0.50, 2, 0.50, 1.00, 3) > rclmat <- matrix(m, ncol = 3, byrow = TRUE) > rc <- reclassify(r, rclmat)
49 Example with the elevation data data Create a raster with just four categories of elevation > m <- c( 0, 300, 1, 300, 400, 2, 400, 500, 3, 500, 650, 4) > rclmat <- matrix(m, nrow = 4, byrow = TRUE) > rc <- reclassify(dem_low, rcl = rclmat)
50
51 SPATIAL ANALYSIS IN R WITH SF AND RASTER Let's practice!
Zev Ross President, ZevRoss Spatial Analysis
SPATIAL ANALYSIS IN R WITH SF AND RASTER Welcome! Zev Ross President, ZevRoss Spatial Analysis Packages we will use in this course Two key packages sf for vectors raster for grids Additional packages discussed
More informationZev Ross President, ZevRoss Spatial Analysis
SPATIAL ANALYSIS IN R WITH SF AND RASTER Buffers and centroids Zev Ross President, ZevRoss Spatial Analysis Use a projected coordinate reference system for spatial analysis As a general rule your layers
More informationSpatial data. Spatial data in R
Spatial data Spatial data in R What do we want to do? What do we want to do? What do we want to do? Today s plan Geographic data types in R Projections Important large-scale data sources Geographic Data
More informationWorking with Geospatial Data in R. Introducing sp objects
Introducing sp objects Data frames aren t a great way to store spatial data > head(ward_sales) ward lon lat group order num_sales avg_price 1 1-123.3128 44.56531 0.1 1 159 311626.9 2 1-123.3122 44.56531
More informationUsing R as a GIS. Alan R. Pearse. 31 July 2017
Using R as a GIS Alan R. Pearse 31 July 2017 Contact Email: ar.pearse@qut.edu.au Contents Data sources A note on syntax What is R? Why do GIS in R? Must-have R packages Overview of spatial data types in
More informationsp sf library(rgdal) library(sp) iris_paris_sp <- readogr("data/parisiris2017.shp", stringsasfactors = FALSE)
sp sf library(rgdal) library(sp) iris_paris_sp
More informationSpatial data analysis with
Spatial data analysis with Spatial is special Complex: geometry and attributes The earth is not flat Size: lots and lots of it, multivariate, time series Special plots: maps First Law of Geography: nearby
More informationSpatial Data in R. Release 1.0. Robert Hijmans
Spatial Data in R Release 1.0 Robert Hijmans Sep 16, 2018 CONTENTS 1 1. Introduction 1 2 2. Spatial data 3 2.1 2.1 Introduction............................................. 3 2.2 2.2 Vector data..............................................
More informationCreating Geo model from Live HANA Calculation View
Creating Geo model from Live HANA Calculation View This document provides instructions for how to prepare a calculation view on a live HANA system with location dimensions for use in the SAP Analytics
More informationSpatial data with R. some applications to weather and agriculture. Joe Wheatley. Dublin R 21 Feb Biospherica Risk
Spatial data with R some applications to weather and agriculture Joe Wheatley Biospherica Risk Dublin R 21 Feb 2013 mm 30 20 10 0 Packages sp vector data raster grid data rgdal input/output, coordinate
More informationPaul Ramsey PostGIS Gotchas. The Things You Didn't Know That Lead to Confusion and Dismay. Wednesday, November 18, 15
Paul Ramsey PostGIS Gotchas The Things You Didn't Know That Lead to Confusion and Dismay Happy GIS day! What do you call the day after GIS day? PostGIS Day is tomorrow Paul Ramsey
More informationSpatial Data Analysis with R
Spatial Data Analysis with R Release 0.1 Robert Hijmans April 12, 2018 CONTENTS 1 1. Introduction 1 2 2. Spatial data 3 2.1 2.1 Introduction............................................. 3 2.2 2.2 Vector
More informationCalculate Burn Severity Patch Sizes with SDMtools in R Anna Talucci 19 May 2017
Calculate Burn Severity Patch Sizes with SDMtools in R Anna Talucci 19 May 2017 Overview This analysis is focused on generating some descriptive statisitcs about patch size for burn severity classes across
More informationHow to convert coordinate system in R
How to convert coordinate system in R Dasapta Erwin Irawan 17 June 2014 Affiliation: Applied Geology Research Division, Faculty of Earth Sciences and Technology, Institut Teknologi Bandung Faculty of Agriculture
More informationSpatial Reference Systems Transformations with Boost.Geometry
Spatial Reference Systems Transformations with Boost.Geometry Adam Wulkiewicz Software Engineer at MySQL Spatial reference system A spatial reference system (SRS) or coordinate reference system (CRS) is
More information4/7/2009. Model: Abstraction of reality following formal rules e.g. Euclidean space for physical space
Model: Abstraction of reality following formal rules e.g. Euclidean space for physical space At different levels: mathematical model (Euclidean space) conceptual design model (ER model) data model (design)
More informationGeoapplications development Control work 1 (2017, Fall)
Page 1 Geoapplications development Control work 1 (2017, Fall) Author: Antonio Rodriges, Oct. 2017 http://rgeo.wikience.org/ Surname, name, patronymic: Group: Date: Signature: Select all correct statements.
More informationPackage angstroms. May 1, 2017
Package angstroms May 1, 2017 Title Tools for 'ROMS' the Regional Ocean Modeling System Version 0.0.1 Helper functions for working with Regional Ocean Modeling System 'ROMS' output. See
More informationSome methods for the quantification of prediction uncertainties for digital soil mapping: Universal kriging prediction variance.
Some methods for the quantification of prediction uncertainties for digital soil mapping: Universal kriging prediction variance. Soil Security Laboratory 2018 1 Universal kriging prediction variance In
More informationIntroduction to the raster package (version 2.6-7)
Introduction to the raster package (version 2.6-7) Robert J. Hijmans November 12, 2017 1 Introduction This vignette describes the R package raster. A raster is a spatial (geographic) data structure that
More informationPackage meteoforecast
Type Package Title Numerical Weather Predictions Version 0.51 Date 2017-04-01 Package meteoforecast April 2, 2017 Access to several Numerical Weather Prediction services both in raster format and as a
More informationPackage velox. R topics documented: December 1, 2017
Type Package Title Fast Raster Manipulation and Extraction Version 0.2.0 Date 2017-11-30 Author Philipp Hunziker Package velox December 1, 2017 Maintainer Philipp Hunziker BugReports
More informationThree-dimensional Analyses
Chapter 6 Three-dimensional Analyses Contents Figures 6.1 Three-dimensional home range......................... 110 6.2 Three-dimensional exploration of digital elevation models (DEMs)... 115 6.1 Example
More informationUnderstanding and Using Geometry, Projections, and Spatial Reference Systems in ArcGIS Annette Locke, Rob Juergens
Understanding and Using Geometry, Projections, and Spatial Reference Systems in ArcGIS Annette Locke, Rob Juergens Introduction We present fundamental concepts necessary for the correct and efficient use
More informationGeographic Information Systems. using QGIS
Geographic Information Systems using QGIS 1 - INTRODUCTION Generalities A GIS (Geographic Information System) consists of: -Computer hardware -Computer software - Digital Data Generalities GIS softwares
More informationPackage nngeo. September 29, 2018
Type Package Title k-nearest Neighbor Join for Spatial Data Version 0.2.2 Package nngeo September 29, 2018 K-nearest neighbor search for projected and non-projected 'sf' spatial layers. Nearest neighbor
More informationFrom data source to data view: A practical guide to uploading spatial data sets into MapX
From data source to data view: A practical guide to uploading spatial data sets into MapX Thomas Piller UNEP/GRID Geneva I Table of contents 1. Adding a new data source to MapX... 1 1.1 Method 1: upload
More informationfile:///c:/users/c_harmak/appdata/local/temp/arc8f36/tmp308d.tmp.htm
Page 1 of 6 FireDistricts_CoB Shapefile Tags Bradenton, boundary, fire districts Summary The best current representation of the City of Bradenton Florida's fire districts and sub-districts. Description
More informationPackage redlistr. May 11, 2018
Package redlistr May 11, 2018 Title Tools for the IUCN Red List of Ecosystems and Species Version 1.0.1 A toolbox created by members of the International Union for Conservation of Nature (IUCN) Red List
More informationExercise 1-1: Using GPS track data to create a field boundary
Exercise 1-1: Using GPS track data to create a field boundary Learning objectives: Add QGIS plugins Create a point vector file from a text file Convert GPS tracking points to a field boundary Data folder:
More informationSPECIFICATION FOR GMT VECTOR DATA FORMAT FOR OGR COMPATABILITY. DRAFT v0.2. Brent Wood
SPECIFICATION FOR GMT VECTOR DATA FORMAT FOR OGR COMPATABILITY DRAFT v0.2 Brent Wood pcreso@pcreso.com Background. The National Institute for Water and Atmospheric Research (NIWA) in New Zealand is funding
More informationRaster Data. James Frew ESM 263 Winter
Raster Data 1 Vector Data Review discrete objects geometry = points by themselves connected lines closed polygons attributes linked to feature ID explicit location every point has coordinates 2 Fields
More informationPackage ExceedanceTools
Type Package Package ExceedanceTools February 19, 2015 Title Confidence regions for exceedance sets and contour lines Version 1.2.2 Date 2014-07-30 Author Maintainer Tools
More informationSpatial data and QGIS
Spatial data and QGIS Xue Jingbo IT Center 2017.08.07 A GIS consists of: Spatial Data. Computer Hardware. Computer Software. Longitude Latitude Disease Date 26.870436-31.909519 Mumps 13/12/2008 26.868682-31.909259
More informationUsing Free and Open Source GIS to Automatically Create Standards- Based Spatial Metadata
Using Free and Open Source GIS to Automatically Create Standards- Based Spatial Metadata Claire Ellul University College London Overview The Problem with Metadata Automation Results Further Work The Problem
More informationSpatial Ecology Lab 6: Landscape Pattern Analysis
Spatial Ecology Lab 6: Landscape Pattern Analysis Damian Maddalena Spring 2015 1 Introduction This week in lab we will begin to explore basic landscape metrics. We will simply calculate percent of total
More informationMiscellaneous Code. Chapter Remove or search for duplicated GPS locations in a data frame. Contents
Chapter 10 Miscellaneous Code Contents 10.1 Remove or search for duplicated GPS locations in a data frame..... 172 10.2 Need to convert back to a matrix to be able to export the data or manipulate the
More informationSpatial Data Models. Raster uses individual cells in a matrix, or grid, format to represent real world entities
Spatial Data Models Raster uses individual cells in a matrix, or grid, format to represent real world entities Vector uses coordinates to store the shape of spatial data objects David Tenenbaum GEOG 7
More informationElec_ISO_RTO_Regions
Page 1 of 8 Elec_ISO_RTO_Regions Shapefile Tags independent system operator, iso, regional transmission organization, rto, utilitiescommunication, economy Summary The S&P Global Platts Independent System
More informationCreating Your Own Ag Database
Appendix B Creating Your Own Ag Database B.1 Creating Your Own Database (Empty Map Set) B.2 Importing Data via Add New Layers B.3 Importing Data via the Map Analysis Tool B.4 Importing Data via the File
More informationElec_ISO_LMP_PricingPoints
Page 1 of 7 Elec_ISO_LMP_PricingPoints Shapefile Tags locational marginal pricing, lmp, independent system operator, iso, nodal, trade, market Summary The S&P Global Platts ISO Nodal Pricing Points geospatial
More informationUsing R as a Geographical Information System
Using R as a Geographical Information System F. Vanwindekens Gembloux June 2016 Contents 1 Shapefiles 2 1.1 Definitions............................... 2 1.2 Loading data............................. 2
More information5 Extract the information of location from the geometry column of PostgreSQL table
5 Extract the information of location from the geometry column of PostgreSQL table Open QGIS and load PostGIS layer buildings and the raster layer Tai_wide_G (optional just to show the basemap). 1 Click
More informationQGIS Application - Bug report #418 QGIS fails to read undefined projection from user datum in shape.prj file
QGIS Application - Bug report #418 QGIS fails to read undefined projection from user datum in shape.prj file 2006-12-03 02:38 AM - neteler-itc-it - Status: Closed Priority: Low Assignee: Magnus Homann
More informationAnalysing Spatial Data in R: Representing Spatial Data
Analysing Spatial Data in R: Representing Spatial Data Roger Bivand Department of Economics Norwegian School of Economics and Business Administration Bergen, Norway 31 August 2007 Object framework To begin
More informationGeological mapping using open
Geological mapping using open source QGIS MOHSEN ALSHAGHDARI -2017- Abstract Geological mapping is very important to display your field work in a map for geologist and others, many geologists face problems
More informationGlobal_Price_Assessments
Page 1 of 6 Global_Price_Assessments Shapefile Thumbnail Not Available Tags assessments, platts, crude oil, global, crude, petroleum Summary The S&P Global Platts Global Price Assessment geospatial data
More informationPractical guidance on mapping and visualisation of crime and social data in QGIS
Practical guidance on mapping and visualisation of crime and social data in QGIS Lesson 4: Mapping of aggregated crime data in QGIS This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike
More informationMonitoring Vegetation in a Changing Climate [D5P3a] Gregory Duveiller August 3, 2018
Monitoring Vegetation in a Changing Climate [D5P3a] Gregory Duveiller August 3, 218 Introduction This document outlines the steps that will be taken in the practicals of Monitoring Vegetation in a Changing
More informationplot(seq(0,10,1), seq(0,10,1), main = "the Title", xlim=c(1,20), ylim=c(1,20), col="darkblue");
R for Biologists Day 3 Graphing and Making Maps with Your Data Graphing is a pretty convenient use for R, especially in Rstudio. plot() is the most generalized graphing function. If you give it all numeric
More informationPackage spatial.tools
Package spatial.tools February 20, 2015 Maintainer Jonathan Asher Greenberg License GPL (>= 2) Title R functions for working with spatial data. Type Package LazyLoad yes Author
More informationSpotfire Extension for OpenSpirit User s Guide
Spotfire Extension for OpenSpirit User s Guide Software Release 1.0.0 August 2012 Important Information SOME TIBCO SOFTWARE EMBEDS OR BUNDLES OTHER TIBCO SOFTWARE. USE OF SUCH EMBEDDED OR BUNDLED TIBCO
More informationGraphic Display of Vector Object
What is GIS? GIS stands for Geographic Information Systems, although the term Geographic Information Science is gaining popularity. A GIS is a software platform for storing, organizing, viewing, querying,
More informationLab 2. Vector and raster data.
Lab 2. Vector and raster data. The goal: To learn about the structure of the vector and raster data types. Objective: Create vector and raster datasets and visualize them. Software for the lab: ARCINFO,
More informationPackage geojsonsf. R topics documented: January 11, Type Package Title GeoJSON to Simple Feature Converter Version 1.3.
Type Package Title GeoJSON to Simple Feature Converter Version 1.3.0 Date 2019-01-11 Package geojsonsf January 11, 2019 Converts Between GeoJSON and simple feature objects. License GPL-3 Encoding UTF-8
More informationMaking Maps: Salamander Species in US. Read in the Data
Anything Arc can do, R can do better OR How to lose Arc in 10 days (part 1/n) UA Summer R Workshop: Week 3 Nicholas M. Caruso Christina L. Staudhammer 14 June 2016 Making Maps: Salamander Species in US
More informationFractionation_Facilities
Page 1 of 8 Fractionation_Facilities Shapefile Thumbnail Not Available Tags ngl, natural gas liquids, gas processing, fractionators Summary The S&P Global Platts Fractionation facilities geospatial data
More informationPackage gfcanalysis. August 29, 2016
Package gfcanalysis August 29, 2016 Version 1.4 Date 2015-11-20 Title Tools for Working with Hansen et al. Global Forest Change Dataset Maintainer Alex Zvoleff Depends R (>=
More informationUnderstanding and Using Geometry, Projections, and Spatial Reference Systems in ArcGIS. Rob Juergens, Melita Kennedy, Annette Locke
Understanding and Using Geometry, Projections, and Spatial Reference Systems in ArcGIS Rob Juergens, Melita Kennedy, Annette Locke Introduction We want to give you a basic understanding of geometry and
More informationLECTURE 2 SPATIAL DATA MODELS
LECTURE 2 SPATIAL DATA MODELS Computers and GIS cannot directly be applied to the real world: a data gathering step comes first. Digital computers operate in numbers and characters held internally as binary
More informationArcCatalog or the ArcCatalog tab in ArcMap ArcCatalog or the ArcCatalog tab in ArcMap ArcCatalog or the ArcCatalog tab in ArcMap
ArcGIS Procedures NUMBER OPERATION APPLICATION: TOOLBAR 1 Import interchange file to coverage 2 Create a new 3 Create a new feature dataset 4 Import Rasters into a 5 Import tables into a PROCEDURE Coverage
More informationWriting functions with the raster package
Writing functions with the raster package Robert J. Hijmans November 12, 2017 1 Introduction The raster package has a number of low-level functions (e.g. to read and write files) that allow you to write
More informationPackage seg. February 15, 2013
Package seg February 15, 2013 Version 0.2-4 Date 2013-01-21 Title A set of tools for residential segregation research Author Seong-Yun Hong, David O Sullivan Maintainer Seong-Yun Hong
More informationVector-Based GIS Data Processing. Chapter 6
Vector-Based GIS Data Processing Chapter 6 Vector Data Model Feature Classes points lines polygons Layers limited to one class of data Figure p. 186 Vector Data Model Shapefiles ArcView non-topological
More informationReview of Cartographic Data Types and Data Models
Review of Cartographic Data Types and Data Models GIS Data Models Raster Versus Vector in GIS Analysis Fundamental element used to represent spatial features: Raster: pixel or grid cell. Vector: x,y coordinate
More informationData Assembly, Part II. GIS Cyberinfrastructure Module Day 4
Data Assembly, Part II GIS Cyberinfrastructure Module Day 4 Objectives Continuation of effective troubleshooting Create shapefiles for analysis with buffers, union, and dissolve functions Calculate polygon
More informationClass #2. Data Models: maps as models of reality, geographical and attribute measurement & vector and raster (and other) data structures
Class #2 Data Models: maps as models of reality, geographical and attribute measurement & vector and raster (and other) data structures Role of a Data Model Levels of Data Model Abstraction GIS as Digital
More informationPackage qualmap. R topics documented: September 12, Type Package
Type Package Package qualmap September 12, 2018 Title Opinionated Approach for Digitizing Semi-Structured Qualitative GIS Data Version 0.1.1 Provides a set of functions for taking qualitative GIS data,
More informationPackage osrm. November 13, 2017
Package osrm November 13, 2017 Type Package Title Interface Between R and the OpenStreetMap-Based Routing Service OSRM Version 3.1.0 Date 2017-11-13 An interface between R and the OSRM API. OSRM is a routing
More informationSmart GIS Course. Developed By. Mohamed Elsayed Elshayal. Elshayal Smart GIS Map Editor and Surface Analysis. First Arabian GIS Software
Smart GIS Course Developed By Mohamed Elsayed Elshayal Elshayal Smart GIS Map Editor and Surface Analysis First Arabian GIS Software http://www.freesmartgis.blogspot.com/ http://tech.groups.yahoo.com/group/elshayalsmartgis/
More informationIntroduction to GIS. Geographic Information Systems SOCR-377 9/24/2015. R. Khosla Fall Semester The real world. What in the world is GIS?
Introduction to GIS Geographic Information Systems SOCR-377 What in the world is GIS? GIS is simultaneously the telescope, the microscope, the computer and the Xerox machine of regional analysis and synthesis
More informationR-ArcGIS Scripting. Data Sources City of Toronto: RD
Tutorial Overview R is an open-source statistical computing language that offers a large suite of data analysis and statistical tools, and is currently the de facto standard for statistical data analysis
More informationFinal Project. Victor October 21, 2017
Final Project Victor October 21, 2017 Memo on final graphic Executive summary The terrorist attack of 9/11 sounds like the epic act of terrorism and often makes on feels like America is the country with
More informationTechnical Specifications
1 Contents INTRODUCTION...3 ABOUT THIS LAB...3 IMPORTANCE OF THIS MODULE...3 EXPORTING AND IMPORTING DATA...4 VIEWING PROJECTION INFORMATION...5...6 Assigning Projection...6 Reprojecting Data...7 CLIPPING/SUBSETTING...7
More informationExercise 03 Creating and Editing Shapefiles Assigned Feb. 2, 2018 Due Feb. 9, 2018
Exercise 03 Creating and Editing Shapefiles Assigned Feb. 2, 2018 Due Feb. 9, 2018 On the class website I've posted an exercise_03_data.zip file which contains a USGS 7.5' quad map of Laramie (as laramie_quad_usgs_1963.tiff)
More informationPackage seg. September 8, 2013
Package seg September 8, 2013 Version 0.3-2 Date 2013-09-08 Title A set of tools for measuring spatial segregation Author Seong-Yun Hong, David O Sullivan Maintainer Seong-Yun Hong
More informationSpatial Data Analysis: Introduction to Raster Processing: Part-3
Spatial Data Analysis: Introduction to Raster Processing: Part-3 Background Geospatial data is becoming increasingly used to solve numerous real-life problems (check out some examples here.) In turn, R
More informationSection 13: data.table
Section 13: data.table Ed Rubin Contents 1 Admin 1 1.1 Announcements.......................................... 1 1.2 Last section............................................ 1 1.3 This week.............................................
More informationGeoreferencing & Spatial Adjustment
Georeferencing & Spatial Adjustment Aligning Raster and Vector Data to the Real World Rotation Differential Scaling Distortion Skew Translation 1 The Problem How are geographically unregistered data, either
More informationGRASS GIS - Introduction
GRASS GIS - Introduction What is a GIS A system for managing geographic data. Information about the shapes of objects. Information about attributes of those objects. Spatial variation of measurements across
More informationAlgorithms for GIS csci3225
Algorithms for GIS csci3225 Laura Toma Bowdoin College Spatial data types and models Spatial data in GIS satellite imagery planar maps surfaces networks point cloud (LiDAR) Spatial data in GIS satellite
More informationQGIS Tutorials Documentation
QGIS Tutorials Documentation Release 0.1 Nathaniel Roth November 30, 2016 Contents 1 Installation 3 1.1 Basic Installation............................................. 3 1.2 Advanced Installation..........................................
More informationPackage geomerge. July 31, 2018
Type Package Title Geospatial Data Integration Version 0.3.1 Date 2018-07-31 Author Karsten Donnay and Andrew M. Linke Maintainer Karsten Donnay Package geomerge July 31, 2018 Geospatial
More informationCSE 512 Course Project Operation Requirements
CSE 512 Course Project Operation Requirements 1. Operation Checklist 1) Geometry union 2) Geometry convex hull 3) Geometry farthest pair 4) Geometry closest pair 5) Spatial range query 6) Spatial join
More informationgeotools: Exporting cartography data from Stata to GIS systems
2018 Canadian Stata Conference Morris J. Wosk Centre for Dialogue, Vancouver, BC geotools: Exporting cartography data from Stata to GIS systems Sergiy Radyakin sradyakin@worldbank.org Development Economics
More informationQuadratics Functions: Review
Quadratics Functions: Review Name Per Review outline Quadratic function general form: Quadratic function tables and graphs (parabolas) Important places on the parabola graph [see chart below] vertex (minimum
More informationObjectives Learn how to work with projections in GMS, and how to combine data from different coordinate systems into the same GMS project.
v. 10.4 GMS 10.4 Tutorial Working with map projections in GMS Objectives Learn how to work with projections in GMS, and how to combine data from different coordinate systems into the same GMS project.
More informationOn Grid: Tools and Techniques to Place Reality Data in a Geographic Coordinate System
RC21940 On Grid: Tools and Techniques to Place Reality Data in a Geographic Coordinate System Seth Koterba Principal Engineer ReCap Autodesk Ramesh Sridharan Principal Research Engineer Infraworks Autodesk
More informationNew Media in Landscape Architecture: Advanced GIS
New Media in Landscape Architecture: Advanced GIS - Projections and Transformations - Version 10.2, English ANHALT UNIVERSITY OF APPLIED SCIENCES Hochschule Anhalt Author: Dr. Matthias Pietsch Tutorial-Version:
More informationWell Unknown ID AKA EPSG: 3857
Well Unknown ID AKA EPSG: 3857 Pamela Kanu November 2016 WGS 1984 WEB MERCATOR ALIASES: AUXILIARY SPHERE, WKID: 3857, WKID: 102100, WKID: 102113, SHERICAL MERCATOR, WGS 84/PSEUDO-MERCATOR, OPEN LAYERS:
More informationGIS Workshop Spring 2016
1/ 14 GIS Geographic Information System. An integrated collection of computer software and data used to view and manage information about geographic places, analyze spatial relationships, and model spatial
More informationTerminals. Shapefile. Thumbnail Not Available. Tags barge, rail, truck, tanker, ports, terminals, crude, refined, products
Page 1 of 11 Terminals Shapefile Thumbnail Not Available Tags barge, rail, truck, tanker, ports, terminals, crude, refined, products Summary The S&P Global Platts Oil Terminals geospatial data layer was
More informationHow to perform a quality check of a new dataset. QGIS Tutorials and Tips
How to perform a quality check of a new dataset QGIS Tutorials and Tips ZanSea zansea@suza.ac.tz 1 Objective GIS datasets can come from many different sources: From a Website. From a USB key given by a
More informationGeoreferencing & Spatial Adjustment 2/13/2018
Georeferencing & Spatial Adjustment The Problem Aligning Raster and Vector Data to the Real World How are geographically unregistered data, either raster or vector, made to align with data that exist in
More informationLecturer 2: Spatial Concepts and Data Models
Lecturer 2: Spatial Concepts and Data Models 2.1 Introduction 2.2 Models of Spatial Information 2.3 Three-Step Database Design 2.4 Extending ER with Spatial Concepts 2.5 Summary Learning Objectives Learning
More informationCoverage data model. Vector-Based Spatial Analysis: Tools Processes. Topological Data Model. Polygons Files. Geographic Information Systems.
GEOG4340 Geographic Information Systems Lecture Four 2013winter Vector-Based Spatial Analysis: Tools Processes Reading materials: Chapter 6 of Intro GIS by J. R. Jensen and R.R. Jensen Cheng. Q., Earth
More informationThe Problem. Georeferencing & Spatial Adjustment. Nature of the problem: For Example: Georeferencing & Spatial Adjustment 2/4/2014
Georeferencing & Spatial Adjustment Aligning Raster and Vector Data to a GIS The Problem How are geographically unregistered data, either raster or vector, made to align with data that exist in geographical
More informationGIS in agriculture scale farm level - used in agricultural applications - managing crop yields, monitoring crop rotation techniques, and estimate
Types of Input GIS in agriculture scale farm level - used in agricultural applications - managing crop yields, monitoring crop rotation techniques, and estimate soil loss from individual farms or agricultural
More informationThe Problem. Georeferencing & Spatial Adjustment. Nature Of The Problem: For Example: Georeferencing & Spatial Adjustment 9/20/2016
Georeferencing & Spatial Adjustment Aligning Raster and Vector Data to the Real World The Problem How are geographically unregistered data, either raster or vector, made to align with data that exist in
More informationPackage tiler. June 9, 2018
Version 0.2.0 Package tiler June 9, 2018 Title Create Geographic and Non-Geographic Map Tiles Creates geographic map tiles from geospatial map files or nongeographic map tiles from simple image files.
More information