mapas.R

library(sf)

## Warning: package 'sf' was built under R version 3.6.1

## Linking to GEOS 3.6.1, GDAL 2.2.3, PROJ 4.9.3

library(ggplot2)

## Warning: package 'ggplot2' was built under R version 3.6.1

library(rgdal)

## Warning: package 'rgdal' was built under R version 3.6.1

## Loading required package: sp

## Warning: package 'sp' was built under R version 3.6.1

## rgdal: version: 1.4-4, (SVN revision 833)
##  Geospatial Data Abstraction Library extensions to R successfully loaded
##  Loaded GDAL runtime: GDAL 2.2.3, released 2017/11/20
##  Path to GDAL shared files: C:/Users/Samsung/Documents/R/win-library/3.6/rgdal/gdal
##  GDAL binary built with GEOS: TRUE 
##  Loaded PROJ.4 runtime: Rel. 4.9.3, 15 August 2016, [PJ_VERSION: 493]
##  Path to PROJ.4 shared files: C:/Users/Samsung/Documents/R/win-library/3.6/rgdal/proj
##  Linking to sp version: 1.3-1

library(gstat)

## Warning: package 'gstat' was built under R version 3.6.1

## Registered S3 method overwritten by 'xts':
##   method     from
##   as.zoo.xts zoo

library(RSAGA)

## Warning: package 'RSAGA' was built under R version 3.6.1

## Loading required package: shapefiles

## Warning: package 'shapefiles' was built under R version 3.6.1

## Loading required package: foreign

## 
## Attaching package: 'shapefiles'

## The following objects are masked from 'package:foreign':
## 
##     read.dbf, write.dbf

## Loading required package: plyr

## Warning: package 'plyr' was built under R version 3.6.1

library(raster)

## Warning: package 'raster' was built under R version 3.6.1

pozos <- st_read("2017HuataBampo.shp")

## Reading layer `2017HuataBampo' from data source `C:\Users\Samsung\Documents\Huatabampo\2017HuataBampo.shp' using driver `ESRI Shapefile'
## Simple feature collection with 293 features and 10 fields
## geometry type:  POINT
## dimension:      XY
## bbox:           xmin: 620553 ymin: 2956532 xmax: 642471 ymax: 2975407
## epsg (SRID):    32612
## proj4string:    +proj=utm +zone=12 +datum=WGS84 +units=m +no_defs

ggplot(pozos) +
  geom_sf()

mapview::mapview(pozos,labels=F)