Merge `gtfs2gps_dt` scripts into a single function about gtfs2gps

[rafapereirabr]

Hi Joao. I've created the function test_gtfs_freq.R to test whether a GTFS feed is frequency based or whether it's a simple feed with a detailed stop_times.txt file. This should help solving this issue

[Joaobazzo]

After merging the functions:

• see if it's possible to keep the same units throughout the code, we keep converting "km" <-> "m" sometimes

[Joaobazzo]

There are basically two scripts,gtfs2gps_single_dt.R and mod_updates.R, which will be pulled today @pedro-andrade-inpe . Here's the approach

gtfs2gps_single_dt.R

gtfs2gps_dt_single <- function(gtfszip, filepath, spatial_resolution = 15, week_days = TRUE){
- 1) read GTFS
- 2) convert into sf object
# Core function

corefun <- function(shapeid){
- 3) spatial interpolation
# test data type

if(test_gtfs_freq(gtfs_data)=="frequency"){
- 4.1) update newstoptimes based on update_freq(all_tripids)
else
- 4.2) update newstoptimes based on update_dt(all_tripids)
}

5) write object
}

6) future.apply for all shapeids

return(output)
}

mod_updates.R

update_freq <- function(tripid,new_stoptimes){
1) build 'id' vector all speed between stop sequences
2) speed estimation
3) update departure times
4) get freq info for that trip

# start frequencies functions
update_newstoptimes_freq(starttimes,freq_temp){

update_departure_stoptimes(i,dt_list){
5.1) update departure stoptimes
return(dt_list)

5.2) function to update stoptimes of each departure (update_departure_stoptimes)
5.3) apply function and return the stop times of all departures from that period (update_departure_stoptimes)
return(departure_stoptimes)

6) update newstoptimes frequency
return(new_stoptimes)

the simple data frame behaves similarly

update_dt <- function(tripid,new_stoptimes){
1) build 'id' vector all speed between stop sequences
2) speed estimation
3) update departure times
4) recalculate time stamps
return(new_stoptimes)

[pedro-andrade-inpe]

This issue was already solved by #36, #48, and 4f37503.