Skip to content

Cnaim updates2 #27

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
NicolaiHandreck wants to merge 10 commits into
base: master
Choose a base branch
from
Open

Cnaim updates2 #27

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
10 commits
Select commit Hold shift + click to select a range
881b939
new
NicolaiHandreck Jul 21, 2022
3e34a1e
added tests
NicolaiHandreck Jul 21, 2022
7ffd3a2
new tests added
NicolaiHandreck Jul 25, 2022
7608288
more added
NicolaiHandreck Jul 25, 2022
00675fe
last tests added
NicolaiHandreck Jul 26, 2022
08b69f2
updated
NicolaiHandreck Aug 3, 2022
4f8a039
Update test-pof_future_transformer_30_60kv.R
NicolaiHandreck Aug 3, 2022
d5b2600
Update pof_cables_60_30kv.R
NicolaiHandreck Aug 3, 2022
25b0bc0
Made Meters, RTU and relays simpler
NicolaiHandreck Sep 15, 2022
e27e2ef
Found mistake in location factor
NicolaiHandreck Feb 21, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
319 changes: 159 additions & 160 deletions R/location_factor.R
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -86,188 +86,187 @@ location_factor <- function(placement = "Default",
sub_division = NULL) {


if (asset_type == "LV UGB" ||
asset_type == "33kV UG Cable (Non Pressurised)" ||
asset_type == "33kV UG Cable (Oil)" ||
asset_type == "33kV UG Cable (Gas)" ||
asset_type == "66kV UG Cable (Non Pressurised)" ||
asset_type == "66kV UG Cable (Oil)" ||
asset_type == "66kV UG Cable (Gas)" ||
asset_type == "132kV UG Cable (Non Pressurised)" ||
asset_type == "132kV UG Cable (Oil)" ||
asset_type == "132kV UG Cable (Gas)" ){
if (asset_type == "LV UGB" ||
asset_type == "33kV UG Cable (Non Pressurised)" ||
asset_type == "33kV UG Cable (Oil)" ||
asset_type == "33kV UG Cable (Gas)" ||
asset_type == "66kV UG Cable (Non Pressurised)" ||
asset_type == "66kV UG Cable (Oil)" ||
asset_type == "66kV UG Cable (Gas)" ||
asset_type == "132kV UG Cable (Non Pressurised)" ||
asset_type == "132kV UG Cable (Oil)" ||
asset_type == "132kV UG Cable (Gas)" ){

location_factor_asset <- 1
location_factor_asset <- 1

} else {
# Find generic term -------------------------------------------------------
asset_category <- gb_ref$categorisation_of_assets$
`Health Index Asset Category`[which(gb_ref$
categorisation_of_assets$
`Asset Register Category` ==
asset_type)]
} else {
# Find generic term -------------------------------------------------------
asset_category <- gb_ref$categorisation_of_assets$
`Health Index Asset Category`[which(gb_ref$
categorisation_of_assets$
`Asset Register Category` ==
asset_type)]

generic_term_1 <- gb_ref$generic_terms_for_assets$
`Generic Term...1`[which(gb_ref$
generic_terms_for_assets$
`Health Index Asset Category` ==
asset_category)]
generic_term_1 <- gb_ref$generic_terms_for_assets$
`Generic Term...1`[which(gb_ref$
generic_terms_for_assets$
`Health Index Asset Category` ==
asset_category)]



if (asset_category == "EHV OHL Conductor (Tower Lines)" ||
asset_category == "132kV OHL Conductor (Tower Lines)") {
generic_term_1 <- "Towers (Conductor)"
} else if (asset_category == "EHV OHL Fittings" ||
asset_category == "132kV OHL Fittings") {
generic_term_1 <- "Towers (Fittings)"
}else if (asset_category == "HV OHL Support - Poles" ||
asset_category == "EHV OHL Support - Poles" ||
asset_category == "LV OHL Support" ) {
# All the poles
if(sub_division %>% is.null())
stop("No sub division specified for the pole")
if(sub_division == "Steel")
generic_term_1 <- "Poles (Steel)"
if(sub_division == "Concrete")
generic_term_1 <- "Poles (Concrete)"
if(sub_division == "Wood")
generic_term_1 <- "Poles (Wood)"
}else if(asset_category == "EHV OHL Support - Towers" ||
asset_category == "132kV OHL Support - Tower"){
generic_term_1 <- "Towers (Structure)"
}
if (asset_category == "Overhead Line") {
stop(paste0("Asset type not implemented: ", asset_type))
}
if (asset_category == "EHV OHL Conductor (Tower Lines)" ||
asset_category == "132kV OHL Conductor (Tower Lines)") {
generic_term_1 <- "Towers (Conductor)"
} else if (asset_category == "EHV OHL Fittings" ||
asset_category == "132kV OHL Fittings") {
generic_term_1 <- "Towers (Fittings)"
}else if (asset_category == "HV OHL Support - Poles" ||
asset_category == "EHV OHL Support - Poles" ||
asset_category == "LV OHL Support" ) {
# All the poles
if(sub_division %>% is.null())
stop("No sub division specified for the pole")
if(sub_division == "Steel")
generic_term_1 <- "Poles (Steel)"
if(sub_division == "Concrete")
generic_term_1 <- "Poles (Concrete)"
if(sub_division == "Wood")
generic_term_1 <- "Poles (Wood)"
}else if(asset_category == "EHV OHL Support - Towers" ||
asset_category == "132kV OHL Support - Tower"){
generic_term_1 <- "Towers (Structure)"
}
if (asset_category == "Overhead Line") {
stop(paste0("Asset type not implemented: ", asset_type))
}

# Altitude ----------------------------------------------------------------
altitude_factor_asset_df <- dplyr::select(gb_ref$altitude_factor_lut,
c("Lower", "Upper",
generic_term_1))
if (altitude_m == "Default") {
row_no <- which(altitude_factor_asset_df$Lower == "Default")
} else if (altitude_m <= 100) {
row_no <- which(altitude_factor_asset_df$Lower == "0")
} else if (100 < altitude_m && altitude_m <= 200) {
row_no <- which(altitude_factor_asset_df$Lower == "100")
} else if (200 < altitude_m && altitude_m <= 300) {
row_no <- which(altitude_factor_asset_df$Lower == "200")
} else {
row_no <- which(altitude_factor_asset_df$Lower == "300")
}
# Altitude ----------------------------------------------------------------
altitude_factor_asset_df <- dplyr::select(gb_ref$altitude_factor_lut,
c("Lower", "Upper",
generic_term_1))
if (altitude_m == "Default") {
row_no <- which(altitude_factor_asset_df$Lower == "Default")
} else if (altitude_m <= 100) {
row_no <- which(altitude_factor_asset_df$Lower == "-Inf")
} else if (100 < altitude_m && altitude_m <= 200) {
row_no <- which(altitude_factor_asset_df$Lower == "100")
} else if (200 < altitude_m && altitude_m <= 300) {
row_no <- which(altitude_factor_asset_df$Lower == "200")
} else {
row_no <- which(altitude_factor_asset_df$Lower == "300")
}

altitude_factor <-
as.numeric(altitude_factor_asset_df[row_no, generic_term_1])
altitude_factor <-
as.numeric(altitude_factor_asset_df[row_no, ..generic_term_1])

# Corrosion ----------------------------------------------------------------
corrosion_category_factor_a <-
dplyr::select(gb_ref$corrosion_category_factor_lut,
c("Corrosion Category Index", generic_term_1))
# Corrosion ----------------------------------------------------------------
corrosion_category_factor_a <-
dplyr::select(gb_ref$corrosion_category_factor_lut,
c("Corrosion Category Index", generic_term_1))

if (corrosion_category_index == "Default") {
row_no <-
which(corrosion_category_factor_a$`Corrosion Category Index` ==
"Default")
} else if (corrosion_category_index == 1) {
row_no <- which(corrosion_category_factor_a$
`Corrosion Category Index` == "1")
} else if (corrosion_category_index == 2) {
row_no <- which(corrosion_category_factor_a$
`Corrosion Category Index` == "2")
} else if (corrosion_category_index == 3) {
row_no <- which(corrosion_category_factor_a$
`Corrosion Category Index` == "3")
} else if (corrosion_category_index == 4) {
row_no <- which(corrosion_category_factor_a$
`Corrosion Category Index` == "4")
}else if (corrosion_category_index == 5) {
row_no <- which(corrosion_category_factor_a$
`Corrosion Category Index` == "5")
}
corrosion_factor <-
as.numeric(corrosion_category_factor_a[row_no, generic_term_1])
if (corrosion_category_index == "Default") {
row_no <-
which(corrosion_category_factor_a$`Corrosion Category Index` ==
"Default")
} else if (corrosion_category_index == 1) {
row_no <- which(corrosion_category_factor_a$
`Corrosion Category Index` == "1")
} else if (corrosion_category_index == 2) {
row_no <- which(corrosion_category_factor_a$
`Corrosion Category Index` == "2")
} else if (corrosion_category_index == 3) {
row_no <- which(corrosion_category_factor_a$
`Corrosion Category Index` == "3")
} else if (corrosion_category_index == 4) {
row_no <- which(corrosion_category_factor_a$
`Corrosion Category Index` == "4")
}else if (corrosion_category_index == 5) {
row_no <- which(corrosion_category_factor_a$
`Corrosion Category Index` == "5")
}
corrosion_factor <-
as.numeric(corrosion_category_factor_a[row_no, ..generic_term_1])

# Distance from coast -----------------------------------------------------
distance_from_coast_factor_lut <- gb_ref$distance_from_coast_factor_lut
# Distance from coast -----------------------------------------------------

distance_from_coast_factor_a <-
dplyr::select(distance_from_coast_factor_lut,
c("Lower", "Upper", generic_term_1))
distance_from_coast_factor_a <-
dplyr::select(gb_ref$distance_from_coast_factor_lut,
c("Lower", "Upper", generic_term_1))

if (distance_from_coast_km == "Default") {
row_no <- which(distance_from_coast_factor_a$Lower == "Default")
} else if (distance_from_coast_km <= 1) {
row_no <- which(distance_from_coast_factor_a$Lower == "0")
} else if (1 < distance_from_coast_km && distance_from_coast_km <= 5) {
row_no <- which(distance_from_coast_factor_a$Lower == "1")
} else if (5 < distance_from_coast_km && distance_from_coast_km <= 10) {
row_no <- which(distance_from_coast_factor_a$Lower == "5")
} else if (10 < distance_from_coast_km && distance_from_coast_km <= 20) {
row_no <- which(distance_from_coast_factor_a$Lower == "10")
}else {
row_no <- which(distance_from_coast_factor_a$Lower == "20")
}
coast_factor <-
as.numeric(distance_from_coast_factor_a[row_no, generic_term_1])

# Increment constant ------------------------------------------------------
increment_constants <- gb_ref$increment_constants
inc_constant <- increment_constants[, generic_term_1]
if (distance_from_coast_km == "Default") {
row_no <- which(distance_from_coast_factor_a$Lower == "Default")
} else if (distance_from_coast_km <= 1) {
row_no <- which(distance_from_coast_factor_a$Lower == "-Inf")
} else if (1 < distance_from_coast_km && distance_from_coast_km <= 5) {
row_no <- which(distance_from_coast_factor_a$Lower == "1")
} else if (5 < distance_from_coast_km && distance_from_coast_km <= 10) {
row_no <- which(distance_from_coast_factor_a$Lower == "5")
} else if (10 < distance_from_coast_km && distance_from_coast_km <= 20) {
row_no <- which(distance_from_coast_factor_a$Lower == "10")
}else {
row_no <- which(distance_from_coast_factor_a$Lower == "20")
}
coast_factor <-
as.numeric(distance_from_coast_factor_a[row_no, ..generic_term_1])

# All factors -------------------------------------------------------------
factors <- c(coast_factor, corrosion_factor, altitude_factor)
# Increment constant ------------------------------------------------------
increment_constants <- gb_ref$increment_constants
inc_constant <- increment_constants[, ..generic_term_1]

# Location factor outdoor -------------------------------------------------
environment_indoor_outdoor <- gb_ref$environment_indoor_outdoor
# All factors -------------------------------------------------------------
factors <- c(coast_factor, corrosion_factor, altitude_factor)

if (placement == "Default") {
placement <-
environment_indoor_outdoor$
`Default 'environment' to be assumed when deriving Location Factor`[
which(
environment_indoor_outdoor$`Asset Register Category` == asset_type)]
}
# Location factor outdoor -------------------------------------------------
environment_indoor_outdoor <- gb_ref$environment_indoor_outdoor

if (placement == "Outdoor") {
if (max(factors) > 1) {
count_factor <- length(which(factors > 1))
location_factor_asset <- max(factors) +
((count_factor - 1) * inc_constant)
} else {
location_factor_asset <- min(factors)
}
} else if (placement == "Indoor") {
if (max(factors) > 1) {
count_factor <- length(which(factors > 1))
initial_location_factor <- max(factors) +
((count_factor - 1) * inc_constant)
} else {
initial_location_factor <- min(factors)
if (placement == "Default") {
placement <-
environment_indoor_outdoor$
`Default 'environment' to be assumed when deriving Location Factor`[
which(
environment_indoor_outdoor$`Asset Register Category` == asset_type)]
}

min_coast_factor <-
min(as.numeric(distance_from_coast_factor_a[, generic_term_1]))
min_corrosion_factor <-
min(as.numeric(corrosion_category_factor_a[, generic_term_1]))
min_altitude_factor <-
min(as.numeric(altitude_factor_asset_df[, generic_term_1]))
min_initial_location_factors <-
c(min_coast_factor, min_corrosion_factor, min_altitude_factor)
if (placement == "Outdoor") {
if (max(factors) > 1) {
count_factor <- length(which(factors > 1))
location_factor_asset <- max(factors) +
((count_factor - 1) * inc_constant)
} else {
location_factor_asset <- min(factors)
}
} else if (placement == "Indoor") {
if (max(factors) > 1) {
count_factor <- length(which(factors > 1))
initial_location_factor <- max(factors) +
((count_factor - 1) * inc_constant)
} else {
initial_location_factor <- min(factors)
}

if (max(min_initial_location_factors) > 1) {
count_min_factor <- length(which(min_initial_location_factors > 1))
min_initial_location_factor <- max(min_initial_location_factors) +
((count_min_factor - 1) * inc_constant)
} else {
min_initial_location_factor <- min(min_initial_location_factors)
min_coast_factor <-
min(distance_from_coast_factor_a[, ..generic_term_1])
min_corrosion_factor <-
min(corrosion_category_factor_a[, ..generic_term_1])
min_altitude_factor <-
min(altitude_factor_asset_df[, ..generic_term_1])
min_initial_location_factors <-
c(min_coast_factor, min_corrosion_factor, min_altitude_factor)

if (max(min_initial_location_factors) > 1) {
count_min_factor <- length(which(min_initial_location_factors > 1))
min_initial_location_factor <- max(min_initial_location_factors) +
((count_min_factor - 1) * inc_constant)
} else {
min_initial_location_factor <- min(min_initial_location_factors)
}
location_factor_asset <- 0.25 *
(initial_location_factor - min_initial_location_factor) +
min_initial_location_factor
}
location_factor_asset <- 0.25 *
(initial_location_factor - min_initial_location_factor) +
min_initial_location_factor
}

}
}



Expand Down
Loading