Созданный SpatialPointsDataFrame не показывает столбцы x, y - PullRequest
Новая
       27

Созданный SpatialPointsDataFrame не показывает столбцы x, y

0 голосов
/ 20 апреля 2020

Я использовал несколько методов для создания SpatialPointsDataFrame из CSV. Все они работают, но после его завершения координаты x и y больше не отображаются в соответствующих столбцах. Даже когда я следовал тому же коду в наборе данных из пакета (SP, meuse), он показывает их в столбцах, а мой - нет. Мне нужны x и y в их столбцах с заголовками для последующего анализа, где мне нужно указать на эти столбцы. Что я делаю не так? 

> bird
          x       y      site use  Crops   Scrub Improved  Water Unimproved   Urban Wetlands Woodland
1  428677.1 2969392    Site 1   0   6.40  771.33   743.15   3.63      53.93   14.35   338.13   529.95
2  424642.0 2965926    Site 2   0 156.91  820.81   323.45  40.38     187.83  155.89   314.61   515.41
3  422386.4 2965106    Site 3   0  40.93  303.08   393.80  10.57     474.57  692.58   247.64   297.64
4  426822.2 2963802    Site 4   0 141.61  967.26   447.54 584.06       0.00    7.02   218.33    95.03
5  429614.9 2959814    Site 6   0   0.00  921.14   671.71  12.38     258.03   16.77   180.78   400.08
6  429856.1 2957563    Site 8   0  17.45  296.79  1007.86 110.46     319.02  343.52   160.26   205.52
7  432358.2 2971414  Site CC1   0 172.92 1459.80   267.44   2.37       0.80    0.00   116.63   440.96
8  431888.7 2974011 Site Off1   0 452.20 1053.80   646.38   2.06     234.44    0.00    37.81    34.22
9  422079.4 2957859 Site Off2   0  54.23  318.36   326.06   6.93     107.65 1174.05   118.01   355.52
10 426358.3 2970949    Site 1   1 760.33  902.48   364.73   1.09      14.69    0.00   166.90   250.67
11 424933.7 2969391    Site 2   1 282.57 1229.56   167.22   1.00       6.03   22.64   118.53   633.32
12 423355.2 2962024    Site 3   1  43.74  262.90   747.57 119.38      83.27  416.18   632.81   240.45
13 428249.8 2966626    Site 4   1   0.00 1181.28   817.34   5.02       7.95   10.68   117.63   320.98
14 429865.0 2963104    Site 6   1 167.77  925.92   864.37   1.63       0.00    6.27   299.17   195.76
15 427126.8 2958142    Site 8   1  17.12  407.58  1050.60  23.85     214.19  124.91   141.23   481.37
16 426679.6 2973693  Site CC1   1 356.26  835.65   769.40  13.74      84.45   15.56    55.53   330.26
17 432716.1 2957642 Site Off1   1 302.42  760.49   563.46   2.90      72.93  269.00    63.37   426.35
18 423757.6 2960735 Site Off2   1 120.46  363.95   997.97  91.36     233.23  205.37   498.81    49.45
> coords<-bird[,c("x","y")]
> databird<-bird[,3:12]
> proj
CRS arguments:
 +proj=utm +zone=17 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +towgs84=0,0,0 
> crs<-proj
> birdframe<-SpatialPointsDataFrame(coords=coords, data=databird, proj4string=proj)
> head(birdframe)
    site use  Crops  Scrub Improved  Water Unimproved  Urban Wetlands Woodland
1 Site 1   0   6.40 771.33   743.15   3.63      53.93  14.35   338.13   529.95
2 Site 2   0 156.91 820.81   323.45  40.38     187.83 155.89   314.61   515.41
3 Site 3   0  40.93 303.08   393.80  10.57     474.57 692.58   247.64   297.64
4 Site 4   0 141.61 967.26   447.54 584.06       0.00   7.02   218.33    95.03
5 Site 6   0   0.00 921.14   671.71  12.38     258.03  16.77   180.78   400.08
6 Site 8   0  17.45 296.79  1007.86 110.46     319.02 343.52   160.26   205.52
> birdframe
class       : SpatialPointsDataFrame 
features    : 18 
extent      : 422079.4, 432716.1, 2957563, 2974011  (xmin, xmax, ymin, ymax)
crs         : +proj=utm +zone=17 +ellps=WGS84 +datum=WGS84 +units=m +no_defs +towgs84=0,0,0 
variables   : 10
names       :      site, use,  Crops,  Scrub, Improved,  Water, Unimproved,   Urban, Wetlands, Woodland 
min values  :    Site 1,   0,      0,  262.9,   167.22,      1,          0,       0,    37.81,    34.22 
max values  : Site Off2,   1, 760.33, 1459.8,   1050.6, 584.06,     474.57, 1174.05,   632.81,   633.32 
> class(birdframe)
[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"
> head(meuse)
       x      y cadmium copper lead zinc  elev       dist   om ffreq soil lime landuse dist.m   log.zn
1 181072 333611    11.7     85  299 1022 7.909 0.00135803 13.6     1    1    1      Ah     50 6.929517
2 181025 333558     8.6     81  277 1141 6.983 0.01222430 14.0     1    1    1      Ah     30 7.039660
3 181165 333537     6.5     68  199  640 7.800 0.10302900 13.0     1    1    1      Ah    150 6.461468
4 181298 333484     2.6     81  116  257 7.655 0.19009400  8.0     1    2    0      Ga    270 5.549076
5 181307 333330     2.8     48  117  269 7.480 0.27709000  8.7     1    2    0      Ah    380 5.594711
6 181390 333260     3.0     61  137  281 7.791 0.36406700  7.8     1    2    0      Ga    470 5.638355
> meuse
class       : SpatialPointsDataFrame 
features    : 155 
extent      : 178605, 181390, 329714, 333611  (xmin, xmax, ymin, ymax)
crs         : NA 
variables   : 15
names       :      x,      y, cadmium, copper, lead, zinc,  elev,     dist, om, ffreq, soil, lime, landuse, dist.m,           log.zn 
min values  : 178605, 329714,     0.2,     14,   37,  113,  5.18,        0,  1,     1,    1,    0,      Aa,     10, 4.72738781871234 
max values  : 181390, 333611,    18.1,    128,  654, 1839, 10.52, 0.880389, 17,     3,    3,    1,       W,   1000, 7.51697722460432 
> class(meuse)
[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"

Посмотрите, как у meuse есть столбцы x и y в столбцах 1 и 2, а у моего нет? Что я делаю не так?

...