[R-sig-Geo] Problem with writeOGR()

[Michael Sumner]

My guess is that it's the logical column.

Try

gpsori1$UTM.Ch <- as.character(gpsori1$UTM.Ch)
writeOGR(gpsori1, dsn="MATA20090729",
layer="MATA20090729gps1",driver="ESRI Shapefile")


On Thu, Nov 5, 2009 at 9:35 AM, Agustin Lobo <aloboaleu at gmail.com> wrote:
> Thanks,
> No POSIXct fields:
>> summary(gpsori1)
> Object of class SpatialPointsDataFrame
> Coordinates:
>              min     max
> coords.x1  453020  457145
> coords.x2 4623500 4625866
> Is projected: NA
> proj4string : [NA]
> Number of points: 2096
> Data attributes:
>       No            UTM.Zone   UTM.Ch           UTM.East
>  Min.   :   1.0   Min.   :31   Mode:logical   Min.   :453020
>  1st Qu.: 524.8   1st Qu.:31   TRUE:2096      1st Qu.:453601
>  Median :1048.5   Median :31   NA's:0         Median :454140
>  Mean   :1048.5   Mean   :31                  Mean   :454418
>  3rd Qu.:1572.2   3rd Qu.:31                  3rd Qu.:455120
>  Max.   :2096.0   Max.   :31                  Max.   :457145
>
>   UTM.North          Altitude            Date            Time
>  Min.   :4623500   Min.   :1648   2009/07/29:2095   08:49:04:   3
>  1st Qu.:4624886   1st Qu.:1795   2009/07/30:   1   08:49:05:   2
>  Median :4625170   Median :1821                     08:49:06:   2
>  Mean   :4625068   Mean   :1835                     08:49:07:   2
>  3rd Qu.:4625338   3rd Qu.:1885                     08:49:08:   2
>  Max.   :4625866   Max.   :1952                     08:49:09:   2
>                                                    (Other) :2083
>> str(gpsori1)
> Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
>  ..@ data       :'data.frame': 2096 obs. of  8 variables:
>  .. ..$ No       : int [1:2096] 1 2 3 4 5 6 7 8 9 10 ...
>  .. ..$ UTM.Zone : int [1:2096] 31 31 31 31 31 31 31 31 31 31 ...
>  .. ..$ UTM.Ch   : logi [1:2096] TRUE TRUE TRUE TRUE TRUE TRUE ...
>  .. ..$ UTM.East : int [1:2096] 453245 453245 453245 453245 453245 453245
> 453245 453245 453245 453245 ...
>  .. ..$ UTM.North: int [1:2096] 4625338 4625338 4625338 4625338 4625338
> 4625338 4625338 4625338 4625338 4625338 ...
>  .. ..$ Altitude : num [1:2096] 1648 1648 1648 1648 1648 ...
>  .. ..$ Date     : Factor w/ 2 levels "2009/07/29","2009/07/30": 1 1 1 1 1 1
> 1 1 1 1 ...
>  .. ..$ Time     : Factor w/ 1452 levels "08:23:05","08:23:07",..: 1 2 3 4 5
> 6 7 8 9 10 ...
>  ..@ coords.nrs : num(0)
>  ..@ coords     : num [1:2096, 1:2] 453245 453245 453245 453245 453245 ...
>  .. ..- attr(*, "dimnames")=List of 2
>  .. .. ..$ : NULL
>  .. .. ..$ : chr [1:2] "coords.x1" "coords.x2"
>  ..@ bbox       : num [1:2, 1:2] 453020 4623500 457145 4625866
>  .. ..- attr(*, "dimnames")=List of 2
>  .. .. ..$ : chr [1:2] "coords.x1" "coords.x2"
>  .. .. ..$ : chr [1:2] "min" "max"
>  ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
>  .. .. ..@ projargs: chr NA
>> class(gpsori1$Date)
> [1] "factor"
>> class(gpsori1$Time)
> [1] "factor"
>
> After your message I thought it could be because of the factors, but the
> following does not work either:
>
>> gpsori1  <-
>> read.table("MATA20090729/MATA20090729_1.csv",sep=",",header=T,stringsAsFactors
>> =F)> coordinates(gpsori1) <- cbind(gpsori1$UTM.East, gpsori1$UTM.North)
>> writeOGR(gpsori1, dsn="MATA20090729",
>> layer="MATA20090729gps1",driver="ESRI Shapefile")
> Error in writeOGR(gpsori1, dsn = "MATA20090729", layer = "MATA20090729gps1",
>  :
>  unknown data type
>
>> summary(gpsori1)
> Object of class SpatialPointsDataFrame
> Coordinates:
>              min     max
> coords.x1  453020  457145
> coords.x2 4623500 4625866
> Is projected: NA
> proj4string : [NA]
> Number of points: 2096
> Data attributes:
>       No            UTM.Zone   UTM.Ch           UTM.East
>  Min.   :   1.0   Min.   :31   Mode:logical   Min.   :453020
>  1st Qu.: 524.8   1st Qu.:31   TRUE:2096      1st Qu.:453601
>  Median :1048.5   Median :31   NA's:0         Median :454140
>  Mean   :1048.5   Mean   :31                  Mean   :454418
>  3rd Qu.:1572.2   3rd Qu.:31                  3rd Qu.:455120
>  Max.   :2096.0   Max.   :31                  Max.   :457145
>   UTM.North          Altitude        Date               Time
>  Min.   :4623500   Min.   :1648   Length:2096        Length:2096
>  1st Qu.:4624886   1st Qu.:1795   Class :character   Class :character
>  Median :4625170   Median :1821   Mode  :character   Mode  :character
>  Mean   :4625068   Mean   :1835
>  3rd Qu.:4625338   3rd Qu.:1885
>  Max.   :4625866   Max.   :1952
>
> Perhaps the NA for proj4string ?
>
> Agus
>
>
> Michael Sumner wrote:
>>
>> Hi Agustin, you will need to determine which columns in the data frame
>> are causing problems.
>>
>> Can you run this and let us know what you see?
>>
>> summary(gpsori1)
>>
>> For example, you cannot write POSIXct columns to SHP and would need to
>> convert to text, or numeri:
>>
>> library(rgdal)
>> d <- data.frame(x = 1:10, y = 1:10, z = 1:10)
>> coordinates(d) <- ~x+y
>>
>> ## OK
>> writeOGR(d, ".", "test", "ESRI Shapefile")
>> unlink("test.shp")
>>
>> d$z <- Sys.time() + 1:10
>> ## not OK
>> writeOGR(d, ".", "test", "ESRI Shapefile")
>>
>>
>> Regards, Mike.
>>
>>
>>
>> On Thu, Nov 5, 2009 at 7:11 AM, Agustin Lobo <alobolistas at gmail.com>
>> wrote:
>>>
>>> I'm puzzled with this:
>>>
>>>> class(gpsori1)
>>>
>>> [1] "SpatialPointsDataFrame"
>>> attr(,"package")
>>> [1] "sp"
>>>>
>>>> writeOGR(gpsori1, dsn="CAST20090907",
>>>> layer="CAST20090907gps1",driver="ESRI Shapefile")
>>>
>>> Error in writeOGR(gpsori1, dsn = "CAST20090907", layer =
>>> "CAST20090907gps1",
>>>  :
>>>  unknown data type
>>>
>>> Am I doing something wrong or is there a recent change that I'm not aware
>>> of?
>>>
>>> Thanks
>>>
>>> Agus
>>>
>>>
>>> _______________________________________________
>>> R-sig-Geo mailing list
>>> R-sig-Geo at stat.math.ethz.ch
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>>
>>>
>>
>
> --
> Dr. Agustin Lobo
> Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
> LLuis Sole Sabaris s/n
> 08028 Barcelona
> Spain
> Tel. 34 934095410
> Fax. 34 934110012
> email: Agustin.Lobo at ija.csic.es
> http://www.ija.csic.es/gt/obster
>