add earth curvature and refraction correction to cartesian_to_polar()

adokter · Apr 1, 2019 · 820e85f956faa0b8ff351bfe7ad6379a68f46193 · 820e85f
1 parent 791af7b
commit 820e85f956faa0b8ff351bfe7ad6379a68f46193
Unified Split

Showing with 2 additions and 2 deletions.

+2 −2 R/project_as_ppi.R
diff --git a/R/project_as_ppi.R b/R/project_as_ppi.R
@@ -195,8 +195,8 @@ proj_to_wgs <- function(x, y, proj4string) {
   return(res)
 }
 
-cartesian_to_polar <- function(coords, elev = 0) {
-  range <- sqrt(coords[, 1]^2 + coords[, 2]^2) / cos(elev)
+cartesian_to_polar <- function(coords, elev = 0, k = 4/3, lat = 35, re = 6378, rp = 6357) {
+  range <- beam_range(sqrt(coords[, 1]^2 + coords[, 2]^2),elev,k=k, lat=lat, re=re, rp=rp)
   azim <- (0.5 * pi - atan2(coords[, 2], coords[, 1])) %% (2 * pi)
   data.frame(range = range, azim = azim * 180 / pi)
 }