Usluge zasnovane na lokaciji - LV6

rm(list=ls())
#!!! vAŽNO!!!Prije izvedbe p=ma, OBAVEZNO instalirati
# sljedede dodatne Knjižnice:
# leaflet, sp, sf, raster, mapview, rgdal
#Položaji točaka
#Koordinate mjesta dobiveneppm1211 alata:
# https://www.mapdevelopers.com/geocode_tool.php
name <- LETTERS[1:10]
# imena tocaka
longitude <- c(14.4258566, 14.3029538, 14.688725, 14.7607367, 14.4087218, 14.4636243, 14.8020393, 14.631247, 14.5395314, 14.7508013)
latitude <- c(45.337594, 45.332782, 45.1740698, 44.7562894, 44.9604914, 44.5380381, 45.3991761, 45.0195884, 45.1588133, 44.9709329)
#Geografske sirine i geografske duzine kao koordinate tocaka
#u podatkovnom okviru 'stations'
stations <- as.data.frame(cbind(longitude, latitude))
# simulacije Kolicina padavina po tockama
set.seed(0)
precip <- round((runif(length(latitude))*10)^3)
psize <- 1 + precip/500
plot(stations, cex=psize, pch=20, col='red', main='Padavine')
# Dodavanje imena stanicama text(stations, name, pos=4)
# Dodavanje legende
breaks <- c(100, 250, 500, 1000)
legend.psize <- 1+breaks/500
legend("topright", legend=breaks, pch=20, pt.cex=legend.psize, col='red', bg='gray')

# crtanje linija
plot(stations, main='Precipitation')
polygon(stations, col='blue', border='light blue')
lines(stations, lwd=3, col='red')
points(stations, cex=2, pch=20)
points(stations, cex=psize, pch=20, col='red', main='Padavine')
# Prikaz podataka na podlozi ogenstreetMap (OSM)
library(leaflet)
map.point <- leaflet() %>%
  addTiles() %>%
  # use the default base map which is openStreetMap tiles
  addMarkers(lng=stations$longitude, lat=stations$latitude)
map.point
# MEĐUTIM, navedenipc.:~1 ip„§ uldiejs nisuREgs.,Wni!
# Kao takvi nisu niti definirani u R-u!


#Kreiranje i prikaz prostornih podataka
library(sp)
# Dodatne R knjiznice 'sp', 'raster', and 'rgdal' su osnova za rad s
# prostornim podatcima!
pts <- SpatialPoints(stations) # stations -> spatialPoints

# spatialPoiwts provjera
class (pts)
showDefault(pts)

#Defirricija referentnog kooridnatnog sustava
crdref <- CRS('+proj=longlat +datum.WGS84')
#spatial_Points s referentnim koordi_natnim sustavom
pts <- SpatialPoints(stations, proj4string=crdref)

# slučajna opažanja kiše (simulacija)
precipvalue <- runif(nrow(stations), min=0, max=100)
df <- data.frame(ID=1:nrow(stations), precip=precipvalue)
ptsdf <- SpatialPointsDataFrame(pts, data=df) # kise dodana u

print(ptsdf) # Provjera podatkovnog okvira
plot(ptsdf) #Graficki prikaz tocaka


# SpatialPoints in R
longitude <- c(-116.7, -120.4, -116.7, -113.5, -115.5, -120.8, -119.5, -113.7, -113.7, 110.7)
latitude <- c(45.3, 42.6, 38.9, 42.1, 35.7, 38.9, 36.2, 39, 41.6, 36.9)
lonlat <- cbind(longitude, latitude)

library(sp)
# 'sp', 'raster', and 'rgdal' are the essential additional R libraries # that manage spatial data!
pts <- SpatialPoints(lonlat) # lonlat -> SpatialPoints

# spatialPoints examination
class (pts)
showDefault(pts)

crdref <- CRS('+proj=longlat +datum=WGS84') #defines the reference co-or
pts <- SpatialPoints(lonlat, proj4string=crdref) #SpatialPoints with co-

# Random precipitation observations
precipvalue <- runif(nrow(lonlat), min=0, max=100)
df <- data.frame(ID=1:nrow(lonlat), precip=precipvalue)
ptsdf <- SpatialPointsDataFrame(pts, data=df) # Precip. assigned to spat

print(ptsdf) # Inspection of spatial point data frame
library(leaflet)
map.point <- leaflet() %>%
  addTiles() %>% # Koristena osnovna OSM podloga
  addPolylines(lng=stations$longitude, lat=stations$latitude)
map.point
library(raster)
#spatialLines and spatialPolygons
lon <- c(-116.8, -114.2, -112.9, -111.9, -114.2, -115.4, -117.7)
lat <- c(41.3, 42.9, 42.4, 39.8, 37.6, 38.3, 37.6)
lonlat <- cbind(lon, lat)
lns <- spLines(lonlat, crs=crdref)
print(lns)

pols <- spPolygons(lonlat, crs=crdref)
print(pols)

plot(pols, axes=TRUE, las=1)
plot(pols, border='blue', col='yellow', lwd=3, add=TRUE)
points(pts, col='red', pch=20, cex=3)

library(leaflet)
map.poly <- leaflet() %>%
  addTiles() %>% # use the default base map which is OpenStreetMap
  addPolygons(lng=lon, lat=lat)
map.poly

rm(list=ls())
library(raster)
r <- raster(ncol=10, nrow=10, xmx=-80, xmn=-150, ymn=20, ymx=60)

values(r) <- runif(ncell(r))
print(r)
values(r) <- 1:ncell(r)
plot(r)
lon <- c(-116.8, -114.2, -112.9, -111.9, -114.2, -115.4, -117.7)
lat <- c(41.3, 42.9, 42.4, 39.8, 37.6, 38.3, 37.6)
lonlat <- cbind(lon, lat)
pols <- spPolygons(lonlat, crs='+proj=longlat +datum=WGS84')
points(lonlat, col='red', pch=20, cex=3)
plot(pols, border='blue', lwd=2, add=TRUE)
library(raster)
filename <- system.file("external/lux.shp", package="raster")
print(filename)

s <- shapefile(filename)
print(s)
plot(s)

f <- system.file("external/rlogo.grd", package="raster")
print(f)

r1 <- raster(f)
print(r1)
plot(r1)
r2 <- raster(f, band=2)
print(r2)
plot(r2)

b <- brick(f)
plot(b)
s <- stack(f)
plot(s)

#Postavljanje referentnog koordinatnog sustava i proiekcije
library(rgdal)
epsg <- make_EPSG()
i <- grep("France", epsg$note, ignore.case=TRUE) # first three
epsg[i[1:3], ]
library(raster)
library(rgdal)
f <- system.file("external/lux.shp", package="raster")
p <- shapefile(f)
print(p)
plot(p)
#proyiera referentnog koordinatnog sustava, ako je nepoznat

#Kreiranje vlastite SHP baze
#Postaviti setwd() na mapu ukojoj se nalazi datoteka 'polozaji.csv
setwd('C:/Users/Krepana Krava/Documents/Faks/lokacije/kod')

#Položaji točaka
stations <- read.csv('polozaji_HR.csv', header = FALSE, sep =',')
colnames(stations) <- c('latitude', 'longitude')
#Geografske širine i geografske dužine kao
#koordinate točaka
#u podatkovnom okviru stations

# Prikaz podataka na podlozi openstreetMap (osm)
library(leaflet)
map.point <- leaflet() %>%
  addTiles() %>% # use the default base
  #map which is openstreetmap tiles
  addMarkers(lng=stations$longitude, lat=stations$latitude)
map.point
crdref <- CRS('+proj=longlat +datum=WGS84')
#spatialPoints s referntnim koordinatnim sustavom
pts <- SpatialPoints(stations, proj4string=crdref)
# slučajna opazanja kiše (simulacija)
precipvalue <- runif(nrow(stations), min=0, max=100)
df <- data.frame(ID=1:nrow(stations), precip=precipvalue)
ptsdf <- SpatialPointsDataFrame(pts, data=df) # Kol. kišel
#dodana u podatkovni okvir
print(ptsdf) # Provjera podatkovnog okvira
plot(ptsdf) #Grafički prikaz točaka

library(leaflet)
map.point <- leaflet() %>%
  addCircleMarkers(lng=stations$longitude, lat=stations$latitude)%>%
  addTiles() %>% # Korištena osnovna
  #05Mpp41~
  addPolylines(lng=stations$longitude, lat=stations$latitude)
map.point

library(leaflet)
map.point <- leaflet() %>%
  addTiles() %>% # Korištena osnovna
  #osm podloga
  addPolylines(lng=stations$longitude, lat=stations$latitude)
map.point