Georeferences the pixels of a scan into a SpatialPointsDataFrame object.
scan_to_spatial(scan, lat, lon, k = 4/3, re = 6378, rp = 6357)
scan | a scan (sweep) of class scan |
---|---|
lat | Geodetic latitude of the radar in degrees. If missing taken from |
lon | Geodetic longitude of the radar in degrees. If missing taken from |
k | Standard refraction coefficient. |
re | Earth equatorial radius in km. |
rp | Earth polar radius in km. |
a SpatialPointsDataFrame
Beam altitude accounts for the curvature of the earth, using beam_height. Distance from the radar over the earth's surface is calculated using beam_distance.
# load example scan: data(example_scan) # convert to a SpatialPointsDataFrame: scan_to_spatial(example_scan)#> class : SpatialPointsDataFrame #> features : 172800 #> extent : -239848.1, 239848.1, -239848.1, 239848.1 (xmin, xmax, ymin, ymax) #> crs : +proj=aeqd +lat_0=56.3675003051758 +lon_0=12.8516998291016 +units=m +ellps=WGS84 #> variables : 9 #> names : azim, range, distance, height, VRADH, DBZH, ZDR, RHOHV, PHIDP #> min values : 0, 500, 499.980668562476, 213.377999126911, -23.8809346854687, -30, -15.9058823511004, 0, -178.588235259056 #> max values : 359, 240000, 239848.089357779, 5695.96400438249, 23.6928958296776, 55, 7.90588282048702, 0.99959997474798, 178.588244199753