Просмотр плотности размещения ячеистой башни

Открыть сценарий в реальном времени

В этом примере показано, как использовать географический график плотности для просмотра плотности размещения сотовой башни в Калифорнии.

Загрузка данных о размещении сотовой башни

Загрузите таблицу данных о размещении вышек сотовой связи в рабочее пространство и просмотрите первые несколько строк. Таблица включает поля, которые идентифицируют местоположение сотовой башни по широте и долготе и идентифицируют тип башни.

load cellularTowers
head(cellularTowers)

ans=8×10 table
     ID     Latitude    Longitude          City             County       State    NEPA    SUPSTRUC    ALLSTRUC    STRUCTYPE
    ____    ________    _________    _________________    ___________    _____    ____    ________    ________    _________

    2166     37.119      -121.83     "MORGAN HILL"        SANTA CLARA     CA       N        46.9        50.6        MAST   
    2167     37.346      -121.63     "SAN JOSE"           SANTA CLARA     CA       N         2.4         2.4        PIPE   
    2168     37.158      -121.98     "REDWOOD ESTATES"    SANTA CLARA     CA       N        24.7        25.3        TOWER  
    2169     37.366      -122.14     "LOS ALTOS HILLS"    SANTA CLARA     CA       N        18.3        19.8        POLE   
    2170     37.402      -122.18     "STANFORD"           SANTA CLARA     CA       N         6.4         6.4        POLE   
    2171     37.258      -122.03     "SARATOGA"           SANTA CLARA     CA       N        10.1        11.9        B      
    2172     37.434      -121.89     "MILPITAS"           SANTA CLARA     CA       N        17.1        17.7        POLE   
    2173     37.446      -121.89     "MILPITAS"           SANTA CLARA     CA       N        19.5        19.5        B

Просмотр данных как графика географического рассеяния

Постройте график данных сотовой башни с помощью geoscatter функция. В сюжете есть четкие области вокруг Сан-Франциско, где количество башен слишком плотное, чтобы быть представленным с помощью участка рассеяния.

geoscatter(cellularTowers.Latitude, cellularTowers.Longitude, '.')
text(gca,37.75,-122.75,'San Francisco','HorizontalAlignment','right')

Просмотр данных в виде графика географической плотности

Плотная площадь башен в районе Сан-Франциско может быть показана с помощью geodensityplot.

geodensityplot(cellularTowers.Latitude, cellularTowers.Longitude)
text(gca,37.75,-122.75,'San Francisco','HorizontalAlignment','right')

Создание графика плотности с указанием радиуса

При создании графика географической плотности по умолчанию график плотности автоматически выбирает значение радиуса, используя данные широты и долготы. Используйте Radius для выбора радиуса в метрах вручную.

radiusInMeters = 50e3; % 50 km
geodensityplot(cellularTowers.Latitude, cellularTowers.Longitude,'Radius',radiusInMeters)

Использовать свойства осей для настройки прозрачности

Если установлено значение 'interp', график плотности FaceAlpha и FaceColor свойства используют Alphamap и Colormap свойства нижележащих географических осей соответственно. Изменение Alphamap изменяет отображение значений плотности на интенсивности цвета.

geodensityplot(cellularTowers.Latitude, cellularTowers.Longitude)
alphamap(normalize((1:64).^0.5,'range'))

AlphaScale свойство на географических осях может также использоваться для изменения прозрачности. Это свойство особенно полезно при попытке показать, где обнаружена какая-либо плотность, а не выделять наиболее плотные области.

figure
dp = geodensityplot(cellularTowers.Latitude, cellularTowers.Longitude)

dp = 
  DensityPlot with properties:

        FaceColor: [0 0.4470 0.7410]
        FaceAlpha: 'interp'
     LatitudeData: [1×1193 double]
    LongitudeData: [1×1193 double]
       WeightData: [1×0 double]
           Radius: 1.8291e+04

  Show all properties

gx = gca

gx = 
  GeographicAxes with properties:

     Basemap: 'streets-light'
    Position: [0.1300 0.1100 0.7750 0.8150]
       Units: 'normalized'

  Show all properties