.csv to .shp in r

- 
# Load Library
library(rgdal)
library(sp)

# Load data.
# Download dataset Facebook https://data.humdata.org/organization/facebook.

data  <- read.csv(file.choose())
head(data,3)

# SpatialPointsDataFrame.

class(data) # data.frame
coordinates(data)<-~longitude+latitude # whatever the equivalent is in your 

class(data) # [1] "SpatialPointsDataFrame"...

# EPSG: 4326.

proj4string(data) = CRS("+init=epsg:4326")

# To .shp

writeOGR(data, "/home/willpecc/Escritorio", "facebook_4326", driver = "ESRI Shapefile")

Popular posts from this blog

How to fill gaps in SLC-Off Landsat 7 ETM+ Images with QGIS?

-
Image
from qgis.core import QgsProcessing from qgis.core import QgsProcessingAlgorithm from qgis.core import QgsProcessingMultiStepFeedback from qgis.core import QgsProcessingParameterRasterLayer from qgis.core import QgsProcessingParameterNumber from qgis.core import QgsProcessingParameterBoolean from qgis.core import QgsProcessingParameterRasterDestination import processing class Gapfilllandsat7(QgsProcessingAlgorithm): def initAlgorithm(self, config=None): self.addParameter(QgsProcessingParameterRasterLayer('inputband', 'Input Band', defaultValue=None)) self.addParameter(QgsProcessingParameterRasterLayer('inputmask', 'Input mask (define qué áreas se deben rellenar)', defaultValue=None)) self.addParameter(QgsProcessingParameterNumber('searchdistanceelnmerodepxelesparabuscarentodaslasdireccionesparainterpolar', 'Search distance (El número de píxeles para buscar en todas las direcciones para interpolar)', type...
Read more