Re: [R-sig-Geo] Merging Shape files

[Dexter Locke]

Within st_join() see the options like "left = TRUE".

You may also want to look at st_intersects()

Here are other helpful reference websites:
https://r-spatial.github.io/sf/reference/geos_binary_pred.html

https://cran.r-project.org/web/packages/sf/vignettes/sf3.html

-Dexter


On Tue, Mar 17, 2020 at 10:29 AM El Mechry El Koudouss <
eelkoudo...@fordham.edu> wrote:

> Dear readers,
> I downloaded a shapefile with data on crime incidents in Washington DC in
> 2016 (available here:
> https://opendata.dc.gov/datasets/crime-incidents-in-2016) and another
> shapefile with data on parks and recreation areas, also in Washington DC
> (available here:
> https://opendata.dc.gov/datasets/parks-and-recreation-areas).
> I was able to read both shapefiles using st_read() from the sf package. My
> question is, how can I merge the two data sets? I tried st_join() but it
> simply resulted in empty columns from the parks dataset getting added to
> the crime data. Any guidance would be much appreciated.
> library(sf)
> parks <- st_read("Parks_and_Recreation_Areas.shp")
> crime <- st_read("Crime_Incidents_in_2016.shp")
> df <- st_join(crime, parks)
>
> --
> El Mechry, El Koudouss (Meshry)
> Graduate Research Assistant
> Center for International Policy Studies
> Fordham University
> Department of Economics
> Website: www.meshry.com
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

[[alternative HTML version deleted]]

___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] Extract VIIRS data in R

[Miluji Sb]

Dear all,

Hope everyone is keeping safe.

I am trying to extract/read VIIRS nighttime lights data but the output
seems rather strange.

I have uploaded the file here
 so as
not exceed the size limit of the email.

## Code ##
library(ncdf4)
library(rgdal)

netcdf_file <- c("VNP02DNB_NRT.A2020069.1048.001.nc")
nl <- brick(netcdf_file, lvar=0, values=TRUE,
varname="observation_data/DNB_observations")

## Detail ##
class  : RasterLayer
dimensions : 3232, 4064, 13134848  (nrow, ncol, ncell)
resolution : 1, 1  (x, y)
extent : 0.5, 4064.5, 0.5, 3232.5  (xmin, xmax, ymin, ymax)
crs: NA
source : C:/Users/shour/Desktop/VNP02DNB_NRT.A2020069.1048.001.nc
names  : DNB.observations.at.pixel.locations
zvar   : observation_data/DNB_observations

The spatial resolution is supposed to be 750m but this shows 1°. What am I
doing wrong? Any help will be greatly appreciated. Thank you.

Sincerely,

Millu

[[alternative HTML version deleted]]

___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] Merging Shape files

[El Mechry El Koudouss]

Dear readers,
I downloaded a shapefile with data on crime incidents in Washington DC in
2016 (available here:
https://opendata.dc.gov/datasets/crime-incidents-in-2016) and another
shapefile with data on parks and recreation areas, also in Washington DC
(available here: https://opendata.dc.gov/datasets/parks-and-recreation-areas).
I was able to read both shapefiles using st_read() from the sf package. My
question is, how can I merge the two data sets? I tried st_join() but it
simply resulted in empty columns from the parks dataset getting added to
the crime data. Any guidance would be much appreciated.
library(sf)
parks <- st_read("Parks_and_Recreation_Areas.shp")
crime <- st_read("Crime_Incidents_in_2016.shp")
df <- st_join(crime, parks)

-- 
El Mechry, El Koudouss (Meshry)
Graduate Research Assistant
Center for International Policy Studies
Fordham University
Department of Economics
Website: www.meshry.com

[[alternative HTML version deleted]]

___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo