This script produces the figures and tables of the manuscript. Each Figur or Table correspond to a section.
The LUCAS 2018 Copernicus Data are downloaded and loaded https://jeodpp.jrc.ec.europa.eu/ftp/jrc-opendata/LUCAS/LUCAS_2018_Copernicus.zip along with the original table containing the 2018 survey from Eurostat website.
library(raster)
## Loading required package: sp
url_lucas_copernicus<-'https://jeodpp.jrc.ec.europa.eu/ftp/jrc-opendata/LUCAS/LUCAS_2018_Copernicus.zip'
url_lucas_survey<-'https://ec.europa.eu/eurostat/cache/lucas/EU_2018_20200213.CSV'
# Download LUCAS Copernicus 2018 survey
download.file(url_lucas_copernicus,destfile=file.path(work.dir,'data',basename(url_lucas_copernicus)))
unzip(file.path(work.dir,'data',basename(url_lucas_copernicus)),exdir=file.path(work.dir,'data'))
# Load the data
copernicus_polygon<-shapefile(file.path(work.dir,'data','LUCAS_2018_Copernicus_polygons.shp')) # read the polygon geometry
lucas_table<-read.csv(file.path(work.dir,'data','LUCAS_2018_Copernicus_attributes.csv')) # read the table
spdf<-merge(x = copernicus_polygon, y = lucas_table, by.x = "POINT_ID",by.y = "POINT_ID") # merge the
# load the survey data
download.file(url_lucas_survey,destfile=file.path(work.dir,'data',basename(url_lucas_survey)))
survey<-read.csv(file.path(work.dir,'data',basename(url_lucas_survey)))
spdf
## class : SpatialPolygonsDataFrame
## features : 63287
## extent : -10.18257, 34.06011, 34.59811, 69.95665 (xmin, xmax, ymin, ymax)
## crs : +proj=longlat +datum=WGS84 +no_defs
## variables : 120
## names : POINT_ID, ID, YEAR, NUTS0, NUTS1, NUTS2, NUTS3, TH_LAT, TH_LONG, OFFICE_PI, EX_ANTE, SURVEY_DATE, CAR_LATITUDE, CAR_EW, CAR_LONGITUDE, ...
## min values : 26461768, 22534, 2018, AT, AT1, AT11, AT111, 34.598316652, -10.18181545, 0, 0, 01/05/18, -4.167545, East, -9.411686, ...
## max values : 65021668, 438280, 2018, UK, UKN, UKN1, UKN16, 69.956514696, 34.059950365, 0, 0, 31/10/18, 69.956338, West, 174.1853, ...
colnames(survey)
## [1] "POINT_ID" "NUTS0" "NUTS1"
## [4] "NUTS2" "NUTS3" "TH_LAT"
## [7] "TH_LONG" "OFFICE_PI" "EX_ANTE"
## [10] "SURVEY_DATE" "CAR_LATITUDE" "CAR_EW"
## [13] "CAR_LONGITUDE" "GPS_PROJ" "GPS_PREC"
## [16] "GPS_ALTITUDE" "GPS_LAT" "GPS_EW"
## [19] "GPS_LONG" "OBS_DIST" "OBS_DIRECT"
## [22] "OBS_TYPE" "LC1" "LC1_SPEC"
## [25] "LC1_PERC" "LC2" "LC2_SPEC"
## [28] "LC2_PERC" "LU1" "LU1_TYPE"
## [31] "LU1_PERC" "LU2" "LU2_TYPE"
## [34] "LU2_PERC" "PARCEL_AREA_HA" "TREE_HEIGHT_SURVEY"
## [37] "TREE_HEIGHT_MATURITY" "FEATURE_WIDTH" "LNDMNG_PLOUGH"
## [40] "LM_PLOUGH_SLOPE" "LM_PLOUGH_DIRECT" "LM_STONE_WALLS"
## [43] "CROP_RESIDUES" "LM_GRASS_MARGINS" "GRAZING"
## [46] "SPECIAL_STATUS" "LC_LU_SPECIAL_REMARK" "CPRN_CANDO"
## [49] "CPRN_LC" "CPRN_LC1N" "CPRNC_LC1E"
## [52] "CPRNC_LC1S" "CPRNC_LC1W" "CPRN_LC1N_BRDTH"
## [55] "CPRN_LC1E_BRDTH" "CPRN_LC1S_BRDTH" "CPRN_LC1W_BRDTH"
## [58] "CPRN_LC1N_NEXT" "CPRN_LC1E_NEXT" "CPRN_LC1S_NEXT"
## [61] "CPRN_LC1W_NEXT" "CPRN_URBAN" "CPRN_IMPERVIOUS_PERC"
## [64] "INSPIRE_PLCC1" "INSPIRE_PLCC2" "INSPIRE_PLCC3"
## [67] "INSPIRE_PLCC4" "INSPIRE_PLCC5" "INSPIRE_PLCC6"
## [70] "INSPIRE_PLCC7" "INSPIRE_PLCC8" "EUNIS_COMPLEX"
## [73] "GRASSLAND_SAMPLE" "GRASS_CANDO" "WM"
## [76] "WM_SOURCE" "WM_TYPE" "WM_DELIVERY"
## [79] "EROSION_CANDO" "SOIL_STONES_PERC" "BIO_SAMPLE"
## [82] "SOIL_BIO_TAKEN" "BULK0_10_SAMPLE" "SOIL_BLK_0_10_TAKEN"
## [85] "BULK10_20_SAMPLE" "SOIL_BLK_10_20_TAKEN" "BULK20_30_SAMPLE"
## [88] "SOIL_BLK_20_30_TAKEN" "STANDARD_SAMPLE" "SOIL_STD_TAKEN"
## [91] "ORGANIC_SAMPLE" "SOIL_ORG_DEPTH_CANDO" "PHOTO_POINT"
## [94] "PHOTO_NORTH" "PHOTO_EAST" "PHOTO_SOUTH"
## [97] "PHOTO_WEST"
LUCAS 2018 points totaling 337 854 points.
##
## 1 2 8
## 1 63364 18924 121673
## 2 0 2488 8671
## 3 0 5844 17050
## 4 0 2 23
## 5 0 0 10
## 6 0 0 2
## 7 0 0 99803
##
## Yes
## In situ < 100 m 63287
Building the Copernicus polygon geometry The Figure was done in QGIS.
Example of information collected by the Copernicus protocol
## POINT_ID ID YEAR NUTS0 NUTS1 NUTS2 NUTS3 TH_LAT TH_LONG OFFICE_PI
## 35003 45223358 258396 2018 DE DE8 DE80 DE80J 53.2925 13.01502 0
## EX_ANTE SURVEY_DATE CAR_LATITUDE CAR_EW CAR_LONGITUDE GPS_PROJ GPS_PREC
## 35003 0 19/07/18 53.29306 East 13.01448 WGS84 4
## GPS_ALTITUDE GPS_LAT GPS_EW GPS_LONG OBS_DIST OBS_DIRECT
## 35003 62 53.29281 East 13.01484 36 On the point
## OBS_TYPE OBS_RADIUS LETTER_GROUP LC1 LC1_LABEL LC1_SPEC
## 35003 In situ < 100 m NA B B32 Rape and turnip rape 8
## LC1_SPEC_LABEL LC1_PERC LC2 LC2_LABEL LC2_SPEC LC2_SPEC_LABEL LC2_PERC
## 35003 Not relevant > 75 % 8 Not relevant 8 Not relevant 0
## LU1 LU1_LABEL LU1_TYPE
## 35003 U111 Agriculture (excluding fallow land and kitchen gardens) 8
## LU1_TYPE_LABEL LU1_PERC LU2 LU2_LABEL LU2_TYPE LU2_TYPE_LABEL LU2_PERC
## 35003 Not relevant > 90 % 8 Not relevant 8 Not relevant 0
## PARCEL_AREA_HA TREE_HEIGHT_MATURITY TREE_HEIGHT_SURVEY FEATURE_WIDTH
## 35003 > 10 ha Not relevant Not relevant Not relevant
## LM_STONE_WALLS CROP_RESIDUES LM_GRASS_MARGINS GRAZING
## 35003 No No > 1 m width No signs of grazing
## SPECIAL_STATUS LC_LU_SPECIAL_REMARK CPRN_CANDO CPRN_LC
## 35003 Hunting No remark Yes B3
## CPRN_LC_LABEL CPRN_LC1N CPRNC_LC1E CPRNC_LC1S CPRNC_LC1W
## 35003 Non-permanent industrial crops 30 51 51 45
## CPRN_LC1N_BRDTH CPRN_LC1E_BRDTH CPRN_LC1S_BRDTH CPRN_LC1W_BRDTH
## 35003 15 888 888 50
## CPRN_LC1N_NEXT CPRN_LC1S_NEXT CPRN_LC1E_NEXT CPRN_LC1W_NEXT CPRN_URBAN
## 35003 E2 8 8 E2 No
## CPRN_IMPERVIOUS_PERC INSPIRE_PLCC1 INSPIRE_PLCC2 INSPIRE_PLCC3
## 35003 5 888 888 888
## INSPIRE_PLCC4 INSPIRE_PLCC5 INSPIRE_PLCC6 INSPIRE_PLCC7 INSPIRE_PLCC8
## 35003 888 888 888 888 888
## EUNIS_COMPLEX GRASSLAND_SAMPLE GRASS_CANDO WM
## 35003 Not relevant 0 Not relevant No visible water management
## WM_SOURCE WM_TYPE WM_DELIVERY EROSION_CANDO SOIL_STONES_PERC
## 35003 Not relevant Not relevant Not relevant Not relevant Not relevant
## BIO_SAMPLE SOIL_BIO_TAKEN BULK0_10_SAMPLE SOIL_BLK_0_10_TAKEN
## 35003 0 8 0 Not relevant
## BULK10_20_SAMPLE SOIL_BLK_10_20_TAKEN BULK20_30_SAMPLE
## 35003 0 Not relevant 0
## SOIL_BLK_20_30_TAKEN STANDARD_SAMPLE SOIL_STD_TAKEN ORGANIC_SAMPLE
## 35003 Not relevant 0 Not relevant 0
## SOIL_ORG_DEPTH_CANDO SOIL_TAKEN SOIL_CROP PHOTO_POINT PHOTO_NORTH
## 35003 Not relevant NA NA Photo taken Photo taken
## PHOTO_SOUTH PHOTO_EAST PHOTO_WEST TRANSECT REVISIT TH_GPS_DIST
## 35003 Photo taken Photo taken Photo taken NA 4 36.75829
## FILE_PATH_GISCO_NORTH
## 35003 https://gisco-services.ec.europa.eu/lucas/photos/2018/DE/452/233/45223358N.jpg
## FILE_PATH_GISCO_SOUTH
## 35003 https://gisco-services.ec.europa.eu/lucas/photos/2018/DE/452/233/45223358S.jpg
## FILE_PATH_GISCO_EAST
## 35003 https://gisco-services.ec.europa.eu/lucas/photos/2018/DE/452/233/45223358E.jpg
## FILE_PATH_GISCO_WEST
## 35003 https://gisco-services.ec.europa.eu/lucas/photos/2018/DE/452/233/45223358W.jpg
## FILE_PATH_GISCO_POINT
## 35003 https://gisco-services.ec.europa.eu/lucas/photos/2018/DE/452/233/45223358P.jpg
## CPRN_LC_SAME_LC1 LUCAS_CORE_INTERSECT COPERNICUS_CLEANED
## 35003 TRUE TRUE TRUE
HIST POLYGONS AREAS (in Hectare)
## Loading required package: rgeos
## rgeos version: 0.5-5, (SVN revision 640)
## GEOS runtime version: 3.8.0-CAPI-1.13.1
## Linking to sp version: 1.4-4
## Polygon checking: TRUE
## Loading required package: ggplot2
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum Unknown based on GRS80 ellipsoid in CRS definition
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
## Warning: Removed 2 rows containing missing values (geom_bar).
Distribution of level-2 land cover classes of the resulting
## Loading required package: stringr
## Class Label Number
## 1 Roofed build-up areas A1 12
## 2 Artificial non-built up areas A2 374
## 3 Other artificial areas A3 21
## 4 Cereals B1 12774
## 5 Root crops B2 877
## 6 Non-permanent industrial crops B3 2435
## 7 Dry pulses, vegetables and flowers B4 767
## 8 Fodder crops B5 2757
## 9 Permanent crops: fruit trees B7 817
## 10 Other permanent crops B8 1148
## 11 Broadleaved woodland C1 8481
## 12 Coniferous woodland C2 5996
## 13 Mixed woodland C3 4484
## 14 Shrubland with sparse tree cover D1 1308
## 15 Shrubland without tree cover D2 1546
## 16 Grassland with sparse tree/shrub cover E1 2078
## 17 Grassland without tree/shrub cover E2 13053
## 18 Spontaneously re-vegetated surfaces E3 2608
## 19 Rocks and stones F1 35
## 20 Sand F2 30
## 21 Lichens and moss F3 3
## 22 Other bare soil F4 1503
## 23 Inland water bodies G1 1
## 24 Inland running water G2 9
## 25 Inland wetlands H1 164
## 26 Coastal wetlands H2 6
## 27 TOTAL 63287
Map of LUCAS Copernicus polygons (N=63287) surveyed in 2018 per level-1 land cover class over the EU28.
## Regions defined for each Polygons
## Coordinate system already present. Adding new coordinate system, which will replace the existing one.
Histogram distribution of LC1 level-3 classes for COPERNICUS_CLEANED
## Warning: Use of `df$LC1` is discouraged. Use `LC1` instead.
LUCAS 208 dataset
## Loading required package: xtable
## [1] 120
## [1] "POINT_ID" "ID" "YEAR"
## [4] "NUTS0" "NUTS1" "NUTS2"
## [7] "NUTS3" "TH_LAT" "TH_LONG"
## [10] "OFFICE_PI" "EX_ANTE" "SURVEY_DATE"
## [13] "CAR_LATITUDE" "CAR_EW" "CAR_LONGITUDE"
## [16] "GPS_PROJ" "GPS_PREC" "GPS_ALTITUDE"
## [19] "GPS_LAT" "GPS_EW" "GPS_LONG"
## [22] "OBS_DIST" "OBS_DIRECT" "OBS_TYPE"
## [25] "OBS_RADIUS" "LETTER_GROUP" "LC1"
## [28] "LC1_LABEL" "LC1_SPEC" "LC1_SPEC_LABEL"
## [31] "LC1_PERC" "LC2" "LC2_LABEL"
## [34] "LC2_SPEC" "LC2_SPEC_LABEL" "LC2_PERC"
## [37] "LU1" "LU1_LABEL" "LU1_TYPE"
## [40] "LU1_TYPE_LABEL" "LU1_PERC" "LU2"
## [43] "LU2_LABEL" "LU2_TYPE" "LU2_TYPE_LABEL"
## [46] "LU2_PERC" "PARCEL_AREA_HA" "TREE_HEIGHT_MATURITY"
## [49] "TREE_HEIGHT_SURVEY" "FEATURE_WIDTH" "LM_STONE_WALLS"
## [52] "CROP_RESIDUES" "LM_GRASS_MARGINS" "GRAZING"
## [55] "SPECIAL_STATUS" "LC_LU_SPECIAL_REMARK" "CPRN_CANDO"
## [58] "CPRN_LC" "CPRN_LC_LABEL" "CPRN_LC1N"
## [61] "CPRNC_LC1E" "CPRNC_LC1S" "CPRNC_LC1W"
## [64] "CPRN_LC1N_BRDTH" "CPRN_LC1E_BRDTH" "CPRN_LC1S_BRDTH"
## [67] "CPRN_LC1W_BRDTH" "CPRN_LC1N_NEXT" "CPRN_LC1S_NEXT"
## [70] "CPRN_LC1E_NEXT" "CPRN_LC1W_NEXT" "CPRN_URBAN"
## [73] "CPRN_IMPERVIOUS_PERC" "INSPIRE_PLCC1" "INSPIRE_PLCC2"
## [76] "INSPIRE_PLCC3" "INSPIRE_PLCC4" "INSPIRE_PLCC5"
## [79] "INSPIRE_PLCC6" "INSPIRE_PLCC7" "INSPIRE_PLCC8"
## [82] "EUNIS_COMPLEX" "GRASSLAND_SAMPLE" "GRASS_CANDO"
## [85] "WM" "WM_SOURCE" "WM_TYPE"
## [88] "WM_DELIVERY" "EROSION_CANDO" "SOIL_STONES_PERC"
## [91] "BIO_SAMPLE" "SOIL_BIO_TAKEN" "BULK0_10_SAMPLE"
## [94] "SOIL_BLK_0_10_TAKEN" "BULK10_20_SAMPLE" "SOIL_BLK_10_20_TAKEN"
## [97] "BULK20_30_SAMPLE" "SOIL_BLK_20_30_TAKEN" "STANDARD_SAMPLE"
## [100] "SOIL_STD_TAKEN" "ORGANIC_SAMPLE" "SOIL_ORG_DEPTH_CANDO"
## [103] "SOIL_TAKEN" "SOIL_CROP" "PHOTO_POINT"
## [106] "PHOTO_NORTH" "PHOTO_SOUTH" "PHOTO_EAST"
## [109] "PHOTO_WEST" "TRANSECT" "REVISIT"
## [112] "TH_GPS_DIST" "FILE_PATH_GISCO_NORTH" "FILE_PATH_GISCO_SOUTH"
## [115] "FILE_PATH_GISCO_EAST" "FILE_PATH_GISCO_WEST" "FILE_PATH_GISCO_POINT"
## [118] "CPRN_LC_SAME_LC1" "LUCAS_CORE_INTERSECT" "COPERNICUS_CLEANED"
Copernicus survey with relation to observation type (OBS_TYPE), observation direction (OBS_DIRECT) and parcel area (PAR-CEL_AREA_HA)
## 1 2 8 TOTAL
## 1 63364 18924 121673 203961
## 2 0 2488 8671 11159
## 3 0 5844 17050 22894
## 4 0 2 23 25
## 5 0 0 10 10
## 6 0 0 2 2
## 7 0 0 99803 99803
## TOTAL 63364 27258 247232 337854
## 1 2 8 TOTAL
## 1 63057 23641 237089 323787
## 2 251 2834 8191 11276
## 3 56 781 1917 2754
## 8 0 2 35 37
## TOTAL 63364 27258 247232 337854
## Warning in rbind(t(tab), TOTAL = margin.table(tab.sum, margin = 2)): number of
## columns of result is not a multiple of vector length (arg 2)
## [,1] [,2] [,3] [,4] [,5] [,6]
## PARCEL_AREA_HA.Var1 "1" "2" "3" "4" "5" "8"
## PARCEL_AREA_HA.Freq " 18346" " 18598" " 19531" "104538" "176814" " 27"
## TOTAL "63364" "27258" "247232" "337854" "63364" "27258"
##
## 1 2 3 4 5 8
## 18346 18598 19531 104538 176814 27
## PARCEL_AREA_HA Number
## 1 1 18346
## 2 2 18598
## 3 3 19531
## 4 4 104538
## 5 5 176814
## 6 8 27
## 7 TOTAL 337854
## 1 2 3 4 5 8 TOTAL
## 1 12488 13864 15294 78354 83961 0 203961
## 2 80 176 299 3824 6780 0 11159
## 3 1640 1271 1044 4844 14095 0 22894
## 4 1 0 0 2 3 19 25
## 5 0 0 0 0 3 7 10
## 6 0 0 0 0 1 1 2
## 7 4137 3287 2894 17514 71971 0 99803
## TOTAL 18346 18598 19531 104538 176814 27 337854
## Class Number
## 1 U111 - Agriculture (excluding fallow land and kitchen gardens) 33603
## 2 U112 - Fallow land 1779
## 3 U113 - Kitchen garden 57
## 4 U120 - Forestry 15990
## 5 U130 - Aquaculture and fishing 1
## 6 U140 - Mining and quarrying 44
## 7 U150 - Other primary production 19
## 8 U210 - Energy production 13
## 9 U221 - Manufacturing of food, beverages and tobacco products 1
## 10 U223 - Coal, oil and metal processing 2
## 11 U224 - Production of non-metal mineral goods 3
## 12 U225 - Chemical and allied industries and manufacturing 1
## 13 U226 - Machinery and equipment 1
## 14 U227 - Wood based products 4
## 15 U311 - Railway transport 12
## 16 U312 - Road transport 163
## 17 U313 - Water transport 3
## 18 U314 - Air transport 12
## 19 U315 - Transport via pipelines 1
## 20 U316 - Telecommunication 2
## 21 U317 - Logistics and storage 17
## 22 U318 - Protection infrastructures 16
## 23 U319 - Electricity, gas and thermal power distribution 20
## 24 U321 - Water supply and treatment 8
## 25 U322 - Waste treatment 10
## 26 U330 - Construction 31
## 27 U341 - Commerce 38
## 28 U342 - Financial, professional and information services 2
## 29 U350 - Community services 99
## 30 U361 - Amenities, museums, leisure 326
## 31 U362 - Sport 167
## 32 U370 - Residential 377
## 33 U411 - Abandoned industrial areas 23
## 34 U412 - Abandoned commercial areas 2
## 35 U413 - Abandoned transport areas 4
## 36 U414 - Abandoned residential areas 30
## 37 U415 - Other abandoned areas 459
## 38 U420 - Semi-natural and natural areas not in use 5086
## 39 TOTAL 58426