Georeferences the pixels of a scan into a SpatialPointsDataFrame object.

scan_to_spatial(scan, lat, lon, k = 4/3, re = 6378, rp = 6357)

Arguments

scan

a scan (sweep) of class scan

lat

Geodetic latitude of the radar in degrees. If missing taken from scan.

lon

Geodetic longitude of the radar in degrees. If missing taken from scan.

k

Standard refraction coefficient.

re

Earth equatorial radius in km.

rp

Earth polar radius in km.

Value

a SpatialPointsDataFrame

Details

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.

Examples

# 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