laba4
A.S. Semin
3 апреля 2019 г
Лабораторная работа №4
Мое задание состоит в том, что надо выбрать любой показатель из раздела “Транспорт” и составить 2 по РФ. Я выбрал такой показатель, как “Объем отправленных грузов(т)”.
Загрузка пакетов
library('R.utils') # gunzip() для распаковки архивов ## Warning: package 'R.utils' was built under R version 3.5.3## Loading required package: R.oo## Loading required package: R.methodsS3## R.methodsS3 v1.7.1 (2016-02-15) successfully loaded. See ?R.methodsS3 for help.## R.oo v1.22.0 (2018-04-21) successfully loaded. See ?R.oo for help.##
## Attaching package: 'R.oo'## The following objects are masked from 'package:methods':
##
## getClasses, getMethods## The following objects are masked from 'package:base':
##
## attach, detach, gc, load, save## R.utils v2.8.0 successfully loaded. See ?R.utils for help.##
## Attaching package: 'R.utils'## The following object is masked from 'package:utils':
##
## timestamp## The following objects are masked from 'package:base':
##
## cat, commandArgs, getOption, inherits, isOpen, parse, warningslibrary('sp') # функция spplot()## Warning: package 'sp' was built under R version 3.5.3library('ggplot2') # функция ggplot()
library('RColorBrewer') # цветовые палитры
require('rgdal') # функция readOGR()## Loading required package: rgdal## Warning: package 'rgdal' was built under R version 3.5.3## rgdal: version: 1.4-3, (SVN revision 828)
## Geospatial Data Abstraction Library extensions to R successfully loaded
## Loaded GDAL runtime: GDAL 2.2.3, released 2017/11/20
## Path to GDAL shared files: C:/Users/Sanya/Documents/R/win-library/3.5/rgdal/gdal
## GDAL binary built with GEOS: TRUE
## Loaded PROJ.4 runtime: Rel. 4.9.3, 15 August 2016, [PJ_VERSION: 493]
## Path to PROJ.4 shared files: C:/Users/Sanya/Documents/R/win-library/3.5/rgdal/proj
## Linking to sp version: 1.3-1##
## Attaching package: 'rgdal'## The following object is masked from 'package:R.oo':
##
## getDescriptionlibrary('broom') # функция tidy()## Warning: package 'broom' was built under R version 3.5.3require('dplyr') # функция join()## Loading required package: dplyr##
## Attaching package: 'dplyr'## The following objects are masked from 'package:stats':
##
## filter, lag## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, unionlibrary('scales') # функция pretty_breaks()
library('mapproj')## Warning: package 'mapproj' was built under R version 3.5.3## Loading required package: maps## Warning: package 'maps' was built under R version 3.5.3library('gpclib')## General Polygon Clipper Library for R (version 1.5-5)
## Type 'class ? gpc.poly' for helplibrary('maptools')## Warning: package 'maptools' was built under R version 3.5.3## Checking rgeos availability: FALSE
## Note: when rgeos is not available, polygon geometry computations in maptools depend on gpclib,
## which has a restricted licence. It is disabled by default;
## to enable gpclib, type gpclibPermit()Объем перевезок за 2017 год
Загрузка данных
gpclibPermit()## Warning in gpclibPermit(): support for gpclib will be withdrawn from
## maptools at the next major release## [1] TRUEShapeFileURL <- "https://biogeo.ucdavis.edu/data/gadm3.6/shp/gadm36_RUS_shp.zip"
if (!file.exists('./data')) dir.create('./data')
if (!file.exists('./data/gadm36_RUS_shp.zip')) {
download.file(ShapeFileURL, destfile = './data/gadm36_RUS_shp.zip')
}
# распаковать архив
unzip('./data/gadm36_RUS_shp.zip', exdir = './data/gadm36_RUS_shp')
# посмотреть список файлов распакованного архива
dir('./data/gadm36_RUS_shp')## [1] "gadm36_RUS_0.cpg" "gadm36_RUS_0.dbf" "gadm36_RUS_0.prj"
## [4] "gadm36_RUS_0.shp" "gadm36_RUS_0.shx" "gadm36_RUS_1.cpg"
## [7] "gadm36_RUS_1.dbf" "gadm36_RUS_1.prj" "gadm36_RUS_1.shp"
## [10] "gadm36_RUS_1.shx" "gadm36_RUS_2.cpg" "gadm36_RUS_2.dbf"
## [13] "gadm36_RUS_2.prj" "gadm36_RUS_2.shp" "gadm36_RUS_2.shx"
## [16] "gadm36_RUS_3.cpg" "gadm36_RUS_3.dbf" "gadm36_RUS_3.prj"
## [19] "gadm36_RUS_3.shp" "gadm36_RUS_3.shx" "license.txt"Выбор иерархии
# прочитать данные уровней 0, 1, 2
Regions0 <- readOGR("./data/gadm36_RUS_shp/gadm36_RUS_0.shp")## OGR data source with driver: ESRI Shapefile
## Source: "C:\Users\Sanya\Desktop\R\4\data\gadm36_RUS_shp\gadm36_RUS_0.shp", layer: "gadm36_RUS_0"
## with 1 features
## It has 2 fieldsRegions1 <- readOGR("./data/gadm36_RUS_shp/gadm36_RUS_1.shp")## OGR data source with driver: ESRI Shapefile
## Source: "C:\Users\Sanya\Desktop\R\4\data\gadm36_RUS_shp\gadm36_RUS_1.shp", layer: "gadm36_RUS_1"
## with 83 features
## It has 10 fieldsRegions2 <- readOGR("./data/gadm36_RUS_shp/gadm36_RUS_2.shp")## OGR data source with driver: ESRI Shapefile
## Source: "C:\Users\Sanya\Desktop\R\4\data\gadm36_RUS_shp\gadm36_RUS_2.shp", layer: "gadm36_RUS_2"
## with 2445 features
## It has 13 fields# контурные карты для разных уровней иерархии
par(mfrow = c(1, 3))
par(oma = c(0, 0, 0, 0))
par(mar = c(0, 0, 1, 0))
plot(Regions0, main = 'adm0', asp = 1.8)
plot(Regions1, main = 'adm1', asp = 1.8)
plot(Regions2, main = 'adm2', asp = 1.8) 
par(mfrow = c(1, 1))
# убрать лишние объекты из памяти
rm(Regions0, Regions2)
# имена слотов
slotNames(Regions1)## [1] "data" "polygons" "plotOrder" "bbox" "proj4string"# слот "данные"
Regions1@data## GID_0 NAME_0 GID_1 NAME_1
## 0 RUS Russia RUS.1_1 Adygey
## 1 RUS Russia RUS.2_1 Altay
## 2 RUS Russia RUS.3_1 Amur
## 3 RUS Russia RUS.4_1 Arkhangel'sk
## 4 RUS Russia RUS.5_1 Astrakhan'
## 5 RUS Russia RUS.6_1 Bashkortostan
## 6 RUS Russia RUS.7_1 Belgorod
## 7 RUS Russia RUS.8_1 Bryansk
## 8 RUS Russia RUS.9_1 Buryat
## 9 RUS Russia RUS.10_1 Chechnya
## 10 RUS Russia RUS.11_1 Chelyabinsk
## 11 RUS Russia RUS.12_1 Chukot
## 12 RUS Russia RUS.13_1 Chuvash
## 13 RUS Russia RUS.14_1 City of St. Petersburg
## 14 RUS Russia RUS.15_1 Dagestan
## 15 RUS Russia RUS.16_1 Gorno-Altay
## 16 RUS Russia RUS.17_1 Ingush
## 17 RUS Russia RUS.18_1 Irkutsk
## 18 RUS Russia RUS.19_1 Ivanovo
## 19 RUS Russia RUS.20_1 Kabardin-Balkar
## 20 RUS Russia RUS.21_1 Kaliningrad
## 21 RUS Russia RUS.22_1 Kalmyk
## 22 RUS Russia RUS.23_1 Kaluga
## 23 RUS Russia RUS.24_1 Kamchatka
## 24 RUS Russia RUS.25_1 Karachay-Cherkess
## 25 RUS Russia RUS.26_1 Karelia
## 26 RUS Russia RUS.27_1 Kemerovo
## 27 RUS Russia RUS.28_1 Khabarovsk
## 28 RUS Russia RUS.29_1 Khakass
## 29 RUS Russia RUS.30_1 Khanty-Mansiy
## 30 RUS Russia RUS.31_1 Kirov
## 31 RUS Russia RUS.32_1 Komi
## 32 RUS Russia RUS.33_1 Kostroma
## 33 RUS Russia RUS.34_1 Krasnodar
## 34 RUS Russia RUS.35_1 Krasnoyarsk
## 35 RUS Russia RUS.36_1 Kurgan
## 36 RUS Russia RUS.37_1 Kursk
## 37 RUS Russia RUS.38_1 Leningrad
## 38 RUS Russia RUS.39_1 Lipetsk
## 39 RUS Russia RUS.40_1 Maga Buryatdan
## 40 RUS Russia RUS.41_1 Mariy-El
## 41 RUS Russia RUS.42_1 Mordovia
## 42 RUS Russia RUS.43_1 Moscow City
## 43 RUS Russia RUS.44_1 Moskva
## 44 RUS Russia RUS.45_1 Murmansk
## 45 RUS Russia RUS.46_1 Nenets
## 46 RUS Russia RUS.47_1 Nizhegorod
## 47 RUS Russia RUS.48_1 North Ossetia
## 48 RUS Russia RUS.49_1 Novgorod
## 49 RUS Russia RUS.50_1 Novosibirsk
## 50 RUS Russia RUS.51_1 Omsk
## 51 RUS Russia RUS.52_1 Orel
## 52 RUS Russia RUS.53_1 Orenburg
## 53 RUS Russia RUS.54_1 Penza
## 54 RUS Russia RUS.55_1 Perm'
## 55 RUS Russia RUS.56_1 Primor'ye
## 56 RUS Russia RUS.57_1 Pskov
## 57 RUS Russia RUS.58_1 Rostov
## 58 RUS Russia RUS.59_1 Ryazan'
## 59 RUS Russia RUS.60_1 Sakha
## 60 RUS Russia RUS.61_1 Sakhalin
## 61 RUS Russia RUS.62_1 Samara
## 62 RUS Russia RUS.63_1 Saratov
## 63 RUS Russia RUS.64_1 Smolensk
## 64 RUS Russia RUS.65_1 Stavropol'
## 65 RUS Russia RUS.66_1 Sverdlovsk
## 66 RUS Russia RUS.67_1 Tambov
## 67 RUS Russia RUS.68_1 Tatarstan
## 68 RUS Russia RUS.69_1 Tomsk
## 69 RUS Russia RUS.70_1 Tula
## 70 RUS Russia RUS.71_1 Tuva
## 71 RUS Russia RUS.72_1 Tver'
## 72 RUS Russia RUS.73_1 Tyumen'
## 73 RUS Russia RUS.74_1 Udmurt
## 74 RUS Russia RUS.75_1 Ul'yanovsk
## 75 RUS Russia RUS.76_1 Vladimir
## 76 RUS Russia RUS.77_1 Volgograd
## 77 RUS Russia RUS.78_1 Vologda
## 78 RUS Russia RUS.79_1 Voronezh
## 79 RUS Russia RUS.80_1 Yamal-Nenets
## 80 RUS Russia RUS.81_1 Yaroslavl'
## 81 RUS Russia RUS.82_1 Yevrey
## 82 RUS Russia RUS.83_1 Zabaykal'ye
## VARNAME_1
## 0 Adygea|Adygeya|Adygheya|Republic of Adygeya|Adygeyskaya A.Obl.|Respublika Adygeya
## 1 Altayskiy Kray
## 2 Amurskaya Oblast
## 3 Arcangelo|Archangel|Archangelsk|Arkhangelskaya Oblast
## 4 Astrachan|Astrakhanskaya Oblast
## 5 Bashkir|Bashkiriya|Bashkirskaya A.S.S.R.|Republic of Bashkortostan|Respublika Bashkortostan
## 6 Belgorodskaya Oblast
## 7 Bryanskaya Oblast
## 8 Buryatiya|Buryat-Mongol A.S.S.R.|Republic of Buryatia|Buryatskaya A.S.S.R.|Respublika
## 9 Cecenia|Chechenia|ChechГЄnia|Tchetchnia|Chechen-Ingush A.S.S.R.|Checheno-Ingushetia|Checheno-Ingushetia
## 10 Chelyabinskaya Oblast
## 11 Chukotka|Chukotskiy Aok
## 12 Chuvashskaya A.S.S.R.|Chuvashskaya Respublika|Chuvashiya|Chuvash Republic
## 13 Sankt-Peterburg gorsovet
## 14 Dagestanskaya A.S.S.R.|Daghestan|Republic of Dagestan|Respublika Dagestan|Dagistan
## 15 Gorno-Altayskaya A.Obl.|Respublika Altay|Oirot|Republic of Altai
## 16 Ingouchie|InguchГ©tia|Inguschetien|Ingushetia|Ingushetiya|Ingush Republic|Ingushskaya Respublika|Respublika Ingushetiya
## 17 Irkutskaya Oblast
## 18 Ivanovskaya Oblast
## 19 Kabardin A.S.S.R.|Kabardino-Balkarskaya A.S.S.R.|Kabardino-Balkariya|Kabardino-Balkarsk|Kabard|Kabardino-Balkarskaya Republic
## 20 Kaliningradskaya Oblast
## 21 Kalmykiya|Khalmg Tangch|Republic of Kalmykia|Kalmytskaya A.S.S.R.|Respublika Kalmykiya
## 22 Kaluzhskaya Oblast
## 23 Kamçatka|Kamchatskaya Kray
## 24 Karaçay-Çerkes|Karachay-Cherkessiya|Karachayevo-Cherkesskaya Respublika|Karachayevo-Cherkessiya|Karachayevo-Cherkess Republic|K
## 25 Karelian A.S.S.R.|Karelo-Finnish A.S.S.R.|Karel'skaya A.S.S.R.|Olonets|Olonetskaya G.|Kareliya|Republic of Karelia
## 26 Kemerovskaya Oblast
## 27 Khabarovskiy Kray
## 28 Khakassiya|Republic of Khakasia|Khakasskaya A.Obl.|Respublika Khakasiya|Republic of Khakasia
## 29 Khanty-Mansiysk|Khanty-Mansiyskiy A.Okr.|Khanty-Mansiyskiy A.Okr.-Yugra|Khanty-Mansiyskiy Aok
## 30 Vyatka|Vyatskaya G.|Kirovskaya Oblast
## 31 Komi A.S.S.R.|Republic of Komi|Respublika Komi
## 32 Kostromskaya Oblast
## 33 Cossacks of the Black Sea|Kuban|Kubanskaya|Yekaterinodar|Krasnodarskiy Kray
## 34 Krasnoyarskiy Kray|Yeniseisk|Yeniseyskaya G.
## 35 Kurganskaya Oblast
## 36 Kurskaya Oblast
## 37 Saint Petersburg|Sankt-Peterburgskaya G.|Leningradskaya Oblast
## 38 Lipetskaya Oblast
## 39 Magadanskaya Oblast
## 40 Mari|Mari-El|Republic of Mari El|Mariyskaya A.S.S.R.|Respublika Mariy El
## 41 Mordov|Mordvian Autonomous Republic|Mordvinia|Republic of Mordovia|Mordovian A.S.S.R.
## 42 <NA>
## 43 Mosca|Moscou|Moscow|Moskau|Moskova|Moskovskaya
## 44 Murmanskaya Oblast
## 45 Nenetskiy A.Okr.|Nenetskiy Aok
## 46 Gor'kiy|Gor'kovskaya|Gorky|Nizhegorodskaya|Nizhniy-Novgorod|Nizhegorodskaya Oblast
## 47 Kuzey Osetya|Respublika Severnaya Osetiya|Severnaya Osetiya-Alaniya|North Ossetian A.S.S.R.|Republic of North Osetia-Alania
## 48 Novgorodskaya Oblast
## 49 Novosibirskaya Oblast
## 50 Omskaya Oblast
## 51 Orlovskaya|Or'ol|Oryol|Orlovskaya Oblast
## 52 Chkalov|Orenburgskaya|Orenburgskaya Oblast
## 53 Penzenskaya Oblast
## 54 Molotov|Permskaya
## 55 KГјsten-Gebiet|Maritime Territory|Primorsk|Primorskiy Kray
## 56 Pskovskaya Oblast
## 57 Province of the Don Cossacks|Provinz des Donischen Heeres|Voyska Donskovo|Rostovskaya Oblast
## 58 Ryazanskaya Oblast|Ryazanskaya Oblast
## 59 Republic of Sakha|Yakutia-Sakha|Yakutsk|Yakut A.S.S.R.|Yakutskaya A.S.S.R.|Respublika Sakha
## 60 Sakhalinskaya Oblast
## 61 Kuybyshev|Kuybyshevskaya|Samarskaya Oblast
## 62 Saratovskaya Oblast
## 63 Smolenskaya Oblast
## 64 Stavropol'skiy Kray|Stavropolskiy Kray
## 65 Yekaterinburg|Sverdlovskaya Oblast
## 66 Tambovskaya Oblast
## 67 Kazan|Kazanskaya G.|Tatar A.S.S.R.|Tatarskaya A.S.S.R.|Republic of Tatarstan|Respublika Tatars
## 68 Tomskaya Oblast
## 69 Tul'skaya|Tulskaya Oblast
## 70 Respublika Tyva|Republic of Tuva|Tyva|Tuvinskaya A.S.S.R.|Republic of Tyva
## 71 Kalinin|Kalininskaya|Tverskaya Oblast
## 72 Tobol'sk|Tobol'skaya G.|Tyumenskaya Oblast
## 73 Udmurtiya|Udmurt Republic|Udmurtskaya A.S.S.R.|Udmurtskaya Respublika
## 74 Simbirsk|Simbirskaya G.|Ul'yanovskaya Oblast
## 75 Vladimirskaya Oblast
## 76 Stalingrad|Volgogradskaya Oblast
## 77 Vologodskaya Oblast
## 78 Voronezhskaya Oblast
## 79 Yamalo-Nenetskiy A. Okrug
## 80 Yaroslavskaya
## 81 Den jГёdiske autonome oblasten|Evrey|Jewish A.Obl.|Yahudi|Yevreyskaya A.Obl.|Evreyskaya AOb
## 82 Zabaykalsky
## NL_NAME_1
## 0 Республика Адыгея
## 1 АлтаР\271СЃРєРёР\271 РєСЂР°Р\271
## 2 Амурская область
## 3 Архангельская область
## 4 Астраханская область
## 5 Республика БаС\210кортостан
## 6 Белгородская область
## 7 Брянская область
## 8 Республика Бурятия
## 9 Республика Чечено-Р\230РЅРіСѓС\210ская
## 10 Челябинская область
## 11 ЧукотскиР\271 РђРћРє
## 12 ЧуваС\210ская Республика
## 13 Санкт-Петербург (горсовет)
## 14 Республика Дагестан
## 15 Республика АлтаР\271
## 16 Респу́блика Р\230РЅРіСѓС\210е́тия
## 17 Р\230ркутская область
## 18 Р\230вановская область
## 19 Кабардино-Балкарская Республика
## 20 Калининградская область
## 21 Республика Калмыкия
## 22 Калужская область
## 23 Камчатская РєСЂР°Р\271
## 24 Карачаево-Черкессия Республика
## 25 Республика Карелия
## 26 Кемеровская область
## 27 ХабаровскиР\271 РєСЂР°Р\271
## 28 Республика Хакасия
## 29 Ханты-МансиР\271СЃРєРёР\271 РђРћРє
## 30 Кировская область
## 31 Республика Коми
## 32 Костромская область
## 33 КраснодарскиР\271 РєСЂР°Р\271
## 34 КрасноярскиР\271 РєСЂР°Р\271
## 35 Курганская область
## 36 Курская область
## 37 Ленинградская область
## 38 Липецкая область
## 39 Магаданская область
## 40 Республика МариР\271 РР»
## 41 Республика Мордовия
## 42 <NA>
## 43 Московская область
## 44 Мурманская область
## 45 НенецкиР\271 РђРћРє
## 46 Нижегородская область
## 47 Республика Северная Осетия-Алания
## 48 Новгородская область
## 49 Новосибирская область
## 50 Омская область
## 51 Орловская область
## 52 Оренбургская область
## 53 Пензенская область
## 54 Пермская РєСЂР°Р\271
## 55 РџСЂРёРјРѕСЂСЃРєРёР\271 РєСЂР°Р\271
## 56 Псковская область
## 57 Ростовская область
## 58 Рязанская область
## 59 Республика Саха
## 60 Сахалинская область
## 61 Самарская область
## 62 Саратовская область
## 63 Смоленская область
## 64 СтавропольскиР\271 РєСЂР°Р\271
## 65 Свердловская область
## 66 Тамбовская область
## 67 Республика Татарстан
## 68 Томская область
## 69 Тульская область
## 70 Республика Тыва
## 71 Тверская область
## 72 Тюменская область
## 73 Удмуртская Республика
## 74 Ульяновская область
## 75 Владимирская область
## 76 Волгоградская область
## 77 Вологодская область
## 78 Воронежская область
## 79 Ямало-НенецкиР\271 РђРћРє
## 80 Ярославская область
## 81 EРІСЂРµР\271ская РђРћР±
## 82 ЗабаР\271кальскиР\271 РєСЂР°Р\271
## TYPE_1 ENGTYPE_1 CC_1 HASC_1
## 0 Respublika Republic <NA> RU.AD
## 1 Kray Territory <NA> RU.AL
## 2 Oblast Region <NA> RU.AM
## 3 Oblast Region <NA> RU.AR
## 4 Oblast Region <NA> RU.AS
## 5 Respublika Republic <NA> RU.BK
## 6 Oblast Region <NA> RU.BL
## 7 Oblast Region <NA> RU.BR
## 8 Respublika Republic <NA> RU.BU
## 9 Respublika Republic <NA> RU.CN
## 10 Oblast Region <NA> RU.CL
## 11 Avtonomnyy Okrug Autonomous Province 7001 RU.CK
## 12 Respublika Republic <NA> RU.CV
## 13 Gorsovet City 601 RU.SP
## 14 Respublika Republic <NA> RU.DA
## 15 Respublika Republic <NA> RU.GA
## 16 Respublika Republic 4201 RU.IN
## 17 Oblast Region <NA> RU.IK
## 18 Oblast Region <NA> RU.IV
## 19 Respublika Republic <NA> RU.KB
## 20 Oblast Region <NA> RU.KN
## 21 Respublika Republic <NA> RU.KL
## 22 Oblast Region <NA> RU.KG
## 23 Kray Territory <NA> RU.KQ
## 24 Respublika Republic <NA> RU.KC
## 25 Respublika Republic <NA> RU.KI
## 26 Oblast Region <NA> RU.KE
## 27 Kray Territory <NA> RU.KH
## 28 Respublika Republic <NA> RU.KK
## 29 Avtonomnyy Okrug Autonomous Province <NA> RU.KM
## 30 Oblast Region <NA> RU.KV
## 31 Respublika Republic <NA> RU.KO
## 32 Oblast Region <NA> RU.KT
## 33 Kray Territory <NA> RU.KD
## 34 Kray Territory 6209 RU.KX
## 35 Oblast Region <NA> RU.KU
## 36 Oblast Region <NA> RU.KS
## 37 Oblast Region <NA> RU.LN
## 38 Oblast Region <NA> RU.LP
## 39 Oblast Region 7009 RU.MG
## 40 Respublika Republic <NA> RU.ME
## 41 Respublika Republic <NA> RU.MR
## 42 Gorod City <NA> <NA>
## 43 Oblast Region <NA> RU.MS
## 44 Oblast Region <NA> RU.MM
## 45 Avtonomnyy Okrug Autonomous Province <NA> RU.NN
## 46 Oblast Region <NA> RU.NZ
## 47 Respublika Republic <NA> RU.NO
## 48 Oblast Region <NA> RU.NG
## 49 Oblast Region <NA> RU.NS
## 50 Oblast Region <NA> RU.OM
## 51 Oblast Region <NA> RU.OL
## 52 Oblast Region <NA> RU.OB
## 53 Oblast Region <NA> RU.PZ
## 54 Kray Territory <NA> RU.PE
## 55 Kray Territory <NA> RU.PR
## 56 Oblast Region <NA> RU.PS
## 57 Oblast Region <NA> RU.RO
## 58 Oblast Region <NA> RU.RZ
## 59 Respublika Republic <NA> RU.SK
## 60 Oblast Region <NA> RU.SL
## 61 Oblast Region <NA> RU.SA
## 62 Oblast Region <NA> RU.SR
## 63 Oblast Region <NA> RU.SM
## 64 Kray Territory <NA> RU.ST
## 65 Oblast Region <NA> RU.SV
## 66 Oblast Region <NA> RU.TB
## 67 Respublika Republic <NA> RU.TT
## 68 Oblast Region <NA> RU.TO
## 69 Oblast Region <NA> RU.TL
## 70 Respublika Republic <NA> RU.TU
## 71 Oblast Region <NA> RU.TV
## 72 Oblast Region <NA> RU.TY
## 73 Respublika Republic <NA> RU.UD
## 74 Oblast Region <NA> RU.UL
## 75 Oblast Region <NA> RU.VL
## 76 Oblast Region <NA> RU.VG
## 77 Oblast Region <NA> RU.VO
## 78 Oblast Region <NA> RU.VR
## 79 Avtonomnyy Okrug Autonomous Province <NA> RU.YN
## 80 Oblast Region <NA> RU.YS
## 81 Avtonomnaya Oblast Autonomous Region 6701 RU.YV
## 82 Kray Territory <NA> RU.ZB# делаем фактор из имён областей (т.е. нумеруем их)
Regions1@data$NAME_1 <- as.factor(Regions1@data$NAME_1)
Regions1@data$NAME_1## [1] Adygey Altay Amur
## [4] Arkhangel'sk Astrakhan' Bashkortostan
## [7] Belgorod Bryansk Buryat
## [10] Chechnya Chelyabinsk Chukot
## [13] Chuvash City of St. Petersburg Dagestan
## [16] Gorno-Altay Ingush Irkutsk
## [19] Ivanovo Kabardin-Balkar Kaliningrad
## [22] Kalmyk Kaluga Kamchatka
## [25] Karachay-Cherkess Karelia Kemerovo
## [28] Khabarovsk Khakass Khanty-Mansiy
## [31] Kirov Komi Kostroma
## [34] Krasnodar Krasnoyarsk Kurgan
## [37] Kursk Leningrad Lipetsk
## [40] Maga Buryatdan Mariy-El Mordovia
## [43] Moscow City Moskva Murmansk
## [46] Nenets Nizhegorod North Ossetia
## [49] Novgorod Novosibirsk Omsk
## [52] Orel Orenburg Penza
## [55] Perm' Primor'ye Pskov
## [58] Rostov Ryazan' Sakha
## [61] Sakhalin Samara Saratov
## [64] Smolensk Stavropol' Sverdlovsk
## [67] Tambov Tatarstan Tomsk
## [70] Tula Tuva Tver'
## [73] Tyumen' Udmurt Ul'yanovsk
## [76] Vladimir Volgograd Vologda
## [79] Voronezh Yamal-Nenets Yaroslavl'
## [82] Yevrey Zabaykal'ye
## 83 Levels: Adygey Altay Amur Arkhangel'sk Astrakhan' ... Zabaykal'yeСтроим сам график
# загружаем статистику с показателями по регионам
fileURL <- 'https://raw.githubusercontent.com/luibrain/luibrain1/laba4/cargos.csv'
stat.Regions <- read.csv2(fileURL, stringsAsFactors = F)
stat.Regions## Region delivered_cargos_2017
## 1 Adygey 0.8
## 2 Altay 18.4
## 3 Amur 9.8
## 4 Arkhangel'sk 11.1
## 5 Astrakhan' 17.3
## 6 Bashkortostan 28.5
## 7 Belgorod 49.9
## 8 Bryansk 10.1
## 9 Buryat 16.1
## 10 Chechnya 0.2
## 11 Chelyabinsk 67.8
## 12 Chukot 0
## 13 Chuvash 0.7
## 14 City of St. Petersburg 6.5
## 15 Dagestan 2.5
## 16 Gorno-Altay 0
## 17 Ingush 0.1
## 18 Irkutsk 50.4
## 19 Ivanovo 0.9
## 20 Kabardin-Balkar 0.5
## 21 Kaliningrad 3.9
## 22 Kalmyk 0
## 23 Kaluga 2.9
## 24 Kamchatka 0
## 25 Karachay-Cherkess 1.3
## 26 Karelia 27.6
## 27 Kemerovo 265.3
## 28 Khabarovsk 21.1
## 29 Khakass 21.2
## 30 Khanty-Mansiy 13.9
## 31 Kirov 6.2
## 32 Komi 13.7
## 33 Kostroma 2.1
## 34 Krasnodar 34.8
## 35 Krasnoyarsk 53.3
## 36 Kurgan 13.8
## 37 Kursk 20.2
## 38 Leningrad 27.1
## 39 Lipetsk 26.6
## 40 Maga Buryatdan 0
## 41 Mariy-El 1.2
## 42 Mordovia 3.1
## 43 Moscow City 4.1
## 44 Moskva 11.3
## 45 Murmansk 29.1
## 46 Nenets 0
## 47 Nizhegorod 15.6
## 48 North Ossetia 0.9
## 49 Novgorod 6.9
## 50 Novosibirsk 21
## 51 Omsk 15.7
## 52 Orel 1.4
## 53 Orenburg 37.3
## 54 Penza 2.4
## 55 Perm' 45.3
## 56 Primor'ye 13.7
## 57 Pskov 3.8
## 58 Rostov 25.3
## 59 Ryazan' 11.5
## 60 Sakha 13.6
## 61 Sakhalin 1.2
## 62 Samara 20.9
## 63 Saratov 15
## 64 Smolensk 16
## 65 Stavropol' 7.3
## 66 Sverdlovsk 0
## 67 Tambov 2.4
## 68 Tatarstan 15.8
## 69 Tomsk 2.6
## 70 Tula 11.7
## 71 Tuva 0
## 72 Tver' 1.6
## 73 Tyumen' 38.4
## 74 Udmurt 2.1
## 75 Ul'yanovsk 3.2
## 76 Vladimir 1.5
## 77 Volgograd 18.5
## 78 Vologda 20.4
## 79 Voronezh 14.6
## 80 Yamal-Nenets 11.9
## 81 Yaroslavl' 12.9
## 82 Yevrey 2.9
## 83 Zabaykal'ye 11.5stat.Regions$delivered_cargos_2017 <- as.numeric(stat.Regions$delivered_cargos_2017)
# вносим данные в файл карты
Regions1@data <- merge(Regions1@data, stat.Regions,
by.x = 'NAME_1', by.y = 'Region')
# задаём палитру
mypalette <- colorRampPalette(c('whitesmoke', 'coral3'))
spplot(Regions1, 'delivered_cargos_2017', main = 'Отправлено грузов (т)',
col.regions = mypalette(10), # цветовая шкала
# (10 градаций)
col = 'coral4', # цвет контурных линий
par.settings = list(axis.line = list(col = NA)) # без
# осей
)
rm(Regions1, stat.Regions)Муниципальные образования за 2014 год
Загрузка данных
gpclibPermit()## Warning in gpclibPermit(): support for gpclib will be withdrawn from
## maptools at the next major release## [1] TRUEfileURL <- 'https://raw.githubusercontent.com/luibrain/luibrain1/laba4/min_obr.csv'
stat.Regions <- read.csv2(fileURL, stringsAsFactors = F)Строим график
Regions <- readOGR(dsn = './data/gadm36_RUS_shp', # папка
layer = 'gadm36_RUS_1') # уровень ## OGR data source with driver: ESRI Shapefile
## Source: "C:\Users\Sanya\Desktop\R\4\data\gadm36_RUS_shp", layer: "gadm36_RUS_1"
## with 83 features
## It has 10 fieldsRegions@data$id <- Regions@data$NAME_1
Regions.points <- fortify(Regions, region = 'id')
Regions.df <- merge(Regions.points, Regions@data, by = 'id')
stat.Regions$id <- stat.Regions$Region
Regions.df <- merge(Regions.df,
stat.Regions[, c('id',
'mun_obr_2014')],
by = 'id')
names(Regions.df)## [1] "id" "long" "lat" "order"
## [5] "hole" "piece" "group" "GID_0"
## [9] "NAME_0" "GID_1" "NAME_1" "VARNAME_1"
## [13] "NL_NAME_1" "TYPE_1" "ENGTYPE_1" "CC_1"
## [17] "HASC_1" "mun_obr_2014"centroids.df <- as.data.frame(coordinates(Regions))
centroids.df$id <- Regions@data$id
colnames(centroids.df) <- c('long', 'lat', 'id')
gp <- ggplot() +
geom_polygon(data = Regions.df,
aes(long, lat, group = group,
fill = mun_obr_2014)) +
geom_path(data = Regions.df,
aes(long, lat, group = group),
color = 'coral4') +
coord_map(projection = 'gilbert', orientation = c(90, 0, 100)) +
scale_fill_distiller(palette = 'OrRd',
direction = 1,
breaks = pretty_breaks(n = 5)) +
labs(x = 'Долгота', y = 'Широта',
title = "Количество муниципальных образований")
gp