forked from impact-initiatives-geospatial/surveyGEER
-
Notifications
You must be signed in to change notification settings - Fork 0
/
dev_history.R
53 lines (44 loc) · 1.65 KB
/
dev_history.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
# usethis::use_build_ignore(files = "dev_history.R")
con <- DBI::dbConnect(RPostgres::Postgres(),
dbname = "global_gdb",
user = rstudioapi::askForPassword("Database user"),
password = rstudioapi::askForPassword("Database password"),
port = 5432)
library(sf);
adm1 <- st_read(con, "som_adm1")
adm1_diss <- adm1 |>
summarise() |>
st_simplify(0.1)
som_boundary <- adm1_diss
usethis::use_data(som_boundary)
library(sf)
dat <- read_csv(here::here("../../HH_HC_MSNA.csv"))
dat2 <- dat |>
select(`_uuid`,lon = `_gps_reading_longitude`,lat= `_gps_reading_latitude`,informed_consent,survey_date, end_survey,electricity_grid,
solar_light,illness_HH_count,`cooking_fuel/collected_firewood`,
`income_source/agricultural_production_sale` ,
agricultural_land ,
`employment_source/agricultural_casual`,
`employment_source/non_agricultural_casual`,
`employment_source/fishing` ) |>
filter(informed_consent=="yes")
dat_sf <- dat2 |>
st_as_sf(crs=4326,coords= c("lon","lat")) |>
st_transform(crs=32646) |>
st_jitter(amount = 150) |>
st_transform(crs=4326 )
dat_df <- dat_sf |>
st_drop_geometry() |>
cbind(st_coordinates(dat_sf)) |>
tibble() |>
rename( `_gps_reading_latitude`=Y,
`_gps_reading_longitude`=X )
dat_df |> glimpse()
dat <- dat_df
usethis::use_data(dat,overwrite=T
)
write_csv(dat,"data/example_raw_data.csv")
usethis::use_git_ignore("data/dat.rda")
usethis::use_git_ignore("data/example_raw_data.csv")
usethis::use_build_ignore("data/example_raw_data.csv")
usethis::use_build_ignore("data/dat.rda")