Outdoor Radio Coverage with Deygout MethodΒΆ

This is an example of path loss prediction with Deygout method with srtm data.

Here a comparison of the coverage for 3 different heights 50,100 ,150 meters.

f = 474.166 MHz Saint Pern

  • ../_images/sphx_glr_plot_cover_001.png
  • ../_images/sphx_glr_plot_cover_002.png
  • ../_images/sphx_glr_plot_cover_003.png

Out:

SRTMDownloader - server= dds.cr.usgs.gov, directory=srtm/version2_1/SRTM3.
N48W002.hgt.zip
 Computation time : 16.65 seconds
 Computation time : 16.86 seconds
 Computation time : 15.45 seconds

from pylayers.gis.ezone import *
import pandas as pd
import matplotlib.pyplot as plt
import numpy as np
import time
#
# Saint Pern Transmitter
#

latgw = 48.28693
longw = -1.95677
tilename = enctile(longw,latgw)
dem=DEM(enctile(longw,latgw))
dem.dwlsrtm()
xgw,ygw = dem.m(longw,latgw)
ez = Ezone(enctile(longw,latgw))
ez.hgts=dem.hgts
ez.rebase('srtm')
pa = np.array([longw,latgw])
for Ht in [50,100,150]:
    tic = time.time()
    L1 = ez.cover(pc=pa, Ht=Ht, Hr=1.5, Nr=200, Rmax=20000, method='deygout',
            source='srtm',fGHz=0.474166)
    toc = time.time()
    print(" Computation time : {:.2f} seconds".format(toc-tic))
    tc = plt.tripcolor(L1[0],-L1[1].flatten(),
                              shading='gouraud',
                              cmap='jet',
                              vmax=-60,
                              vmin=-120,
                              alpha = 1,
                              edgecolors='k',
                              linewidth=0.0)

    plt.title("Lat {:.2f}, Lon {:.2f} , Height{:.2f}".format(latgw,longw,Ht),fontsize=18)
    cb = plt.colorbar()
    cb.set_label("Loss (dB)")
    plt.axis('equal')
plt.show()

Total running time of the script: ( 0 minutes 58.258 seconds)

Gallery generated by Sphinx-Gallery