[GRASS-stats] rgrass7 read/write SpatialPolygonsDataFrames errors
Eduardo Diez
eduardodiez at gmx.com
Sat Oct 10 09:20:21 PDT 2015
This is what i got:
> vinfo0 <- execGRASS("v.info", flags="c", map=zm.gnrl)INTEGER|cat
INTEGER|ID
INTEGER|GRIDCODE
Displaying column types/names for database connection of layer <1>:
And trying to read back the original polygon (without v.clean):
> zm.fnl <- readVECT(zm.pol)Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, :
line 4 did not have 2 elements
2015-10-10 13:12 GMT-03:00 Roger Bivand <Roger.Bivand at nhh.no>:
> On Sat, 10 Oct 2015, Eduardo Diez wrote:
>
> Apparently "v.info" through execGRASS is expecting the "layer" argument to
>> be a string rather than an integer.
>>
>
> Yes, please re-run with "1"; this resulted from GRASS6 to GRASS7
> over-enthusiastic changing everything just for fun - this isn't the
> problem. Run (assuming that zm.gnrl contains a character string):
>
> vinfo0 <- execGRASS("v.info", flags="c", map=zm.gnrl)
>
> My guess is that v.info is saying that this map is broken, or suffering
> from something that vColumns does not know about.
>
> v.info -c map=
>
> should simply respond with "<storage type>|<name>" lines for the fields in
> the layer. Here, its fourth line (unlike the first three) does not have two
> elements (the type/name pairs).
>
> Does readVECT() work if you do not use v.clean? Is v.clean having a side
> effect of causing v.info -c to "say" something else? How should the input
> columns of the object to be cleaned be turned into output columns (simply
> copy across the retained areas?)?
>
> Roger
>
>
>>
>> zm.fnl <- readVECT(zm.gnrl)Error in scan(file, what, nmax, sep, dec,
>>> quote, skip, nlines, na.strings, :
>>>
>> line 4 did not have 2 elements> vinfo0 <- execGRASS("v.info",
>> flags="c", map=zm.gnrl, layer=1L,+ intern=TRUE)Error in doGRASS(cmd,
>> flags = flags, ..., parameters = parameters, echoCmd = echoCmd, :
>> Parameter <layer> does not have string value>
>> parseGRASS("v.info")Command: v.info
>>
>> Extension: .exe
>> Description: Outputs basic information about a vector map.
>> Keywords: vector, metadata, topology, extent, history, attribute columns
>> Parameters:
>> name: map, type: string, required: yes, multiple: no
>> keydesc: name, keydesc_count: 1
>> [Or data source for direct OGR access]
>> name: layer, type: string, required: no, multiple: no
>> default: 1
>> [Vector features can have category values in different layers. This
>> number determines which layer to use. When used with direct OGR access
>> this is the layer name.]
>> Flags:
>> name: h [Print history instead of info and exit] {FALSE}
>> name: c [Print types/names of table columns for specified layer
>> instead of info and exit] {FALSE}
>> name: g [Print region info in shell script style] {FALSE}
>> name: e [Print extended metadata info in shell script style] {FALSE}
>> name: t [Print topology info in shell script style] {FALSE}
>> name: help [Print usage summary] {FALSE}
>> name: verbose [Verbose module output] {FALSE}
>> name: quiet [Quiet module output] {FALSE}
>>
>>
>>
>>>
>> 2015-10-10 10:56 GMT-03:00 Roger Bivand <Roger.Bivand at nhh.no>:
>>
>> OK, thanks. So the failure is in:
>>>
>>> vinfo0 <- execGRASS("v.info", flags="c", map="<name>", layer=1L,
>>> intern=TRUE)
>>>
>>> used to return the column/field names. Please repeat with that.
>>>
>>> Roger
>>>
>>>
>>> On Fri, 9 Oct 2015, Eduardo Diez wrote:
>>>
>>> Here are the traceback and the vInfo(zm.gnrl). The latter seems to work
>>>
>>>> fine.
>>>>
>>>> zm.fnl <- readVECT(zm.gnrl)
>>>>
>>>>>
>>>>> Error in scan(file, what, nmax, sep, dec, quote, skip, nlines,
>>>> na.strings,
>>>> :
>>>> line 4 did not have 2 elements
>>>>
>>>> traceback()
>>>>>
>>>>> 8: scan(file = file, what = what, sep = sep, quote = quote, dec = dec,
>>>> nmax = nrows, skip = 0, na.strings = na.strings, quiet = TRUE,
>>>> fill = fill, strip.white = strip.white, blank.lines.skip =
>>>> blank.lines.skip,
>>>> multi.line = FALSE, comment.char = comment.char, allowEscapes =
>>>> allowEscapes,
>>>> flush = flush, encoding = encoding, skipNul = skipNul)
>>>> 7: read.table(con, header = FALSE, sep = "|")
>>>> 6: vColumns(vname)
>>>> 5: nchar(as.character(vColumns(vname)$name))
>>>> 4: .read_vect_non_plugin(vname = vname, layer = layer, type = type,
>>>> remove.duplicates = remove.duplicates, ignore.stderr =
>>>> ignore.stderr,
>>>> pointDropZ = pointDropZ, driver = driver, with_prj = with_prj,
>>>> with_c = with_c)
>>>> 3: tryCatchList(expr, classes, parentenv, handlers)
>>>> 2: tryCatch({
>>>> stopifnot(is.logical(plugin) || is.null(plugin))
>>>> stopifnot(is.logical(ignore.stderr))
>>>> if (missing(layer))
>>>> layer <- 1L
>>>> layer <- as.character(layer)
>>>> stopifnot(is.logical(with_c))
>>>> if (!is.null(driver) && driver == "GRASS")
>>>> plugin <- TRUE
>>>> if (!requireNamespace("rgdal", quietly = TRUE)) {
>>>> stop("rgdal not available")
>>>> }
>>>> if (is.null(plugin)) {
>>>> ogrD <- rgdal::ogrDrivers()$name
>>>> plugin <- "GRASS" %in% ogrD
>>>> }
>>>> if (plugin) {
>>>> res <- .read_vect_plugin(vname = vname, layer = layer,
>>>> type = type, ignore.stderr = ignore.stderr, pointDropZ =
>>>> pointDropZ,
>>>> mapset = mapset)
>>>> }
>>>> else {
>>>> res <- .read_vect_non_plugin(vname = vname, layer = layer,
>>>> type = type, remove.duplicates = remove.duplicates,
>>>> ignore.stderr = ignore.stderr, pointDropZ = pointDropZ,
>>>> driver = driver, with_prj = with_prj, with_c = with_c)
>>>> }
>>>> }, finally = {
>>>> if (get.suppressEchoCmdInFuncOption()) {
>>>> tull <- set.echoCmdOption(inEchoCmd)
>>>> }
>>>> })
>>>> 1: readVECT(zm.gnrl)
>>>>
>>>> vInfo(zm.gnrl)
>>>>>
>>>>> nodes points lines boundaries centroids areas
>>>> islands primitives map3d
>>>> 274 0 0 302 31 31
>>>> 3 333 0
>>>>
>>>> Thanks
>>>>
>>>>
>>>> 2015-10-09 11:36 GMT-03:00 Roger Bivand <Roger.Bivand at nhh.no>:
>>>>
>>>> If it is the v.info call, please run:
>>>>
>>>>>
>>>>> res <- execGRASS("v.info", flags="t", map="zm.gnrl", layer=1L,
>>>>> intern=TRUE)
>>>>>
>>>>> and report the output of:
>>>>>
>>>>> print(res)
>>>>>
>>>>> and maybe send me offlist the output of
>>>>>
>>>>> save(res, file="vinfo_odd_output.RData")
>>>>>
>>>>> Roger
>>>>>
>>>>> --
>>>>> Roger Bivand
>>>>> Norwegian School of Economics
>>>>> Helleveien 30, 5045 Bergen, Norway
>>>>> Roger.Bivand at nhh.no
>>>>>
>>>>>
>>>>> ________________________________________
>>>>> Fra: grass-stats [grass-stats-bounces at lists.osgeo.org] på vegne
>>>>> av
>>>>> Roger Bivand [Roger.Bivand at nhh.no]
>>>>> Sendt: 9. oktober 2015 4:27
>>>>> Til: Eduardo Diez
>>>>> Kopi: grass-stats at lists.osgeo.org
>>>>> Emne: Re: [GRASS-stats] rgrass7 read/write SpatialPolygonsDataFrames
>>>>> errors
>>>>>
>>>>> OK, please provide the output of traceback() run immediately after the
>>>>> failing readVECT(). That will show which command invokes scan(). Also
>>>>> try
>>>>> to run vInfo(zm.gnrl) - I think the failure is in:
>>>>>
>>>>> execGRASS("v.info", flags="t", map="zm.gnrl", layer=1L)
>>>>>
>>>>> not returning a structure that vInfo understands, and which ois then
>>>>> used
>>>>> in readVECT().
>>>>>
>>>>> Roger
>>>>>
>>>>> --
>>>>> Roger Bivand
>>>>> Norwegian School of Economics
>>>>> Helleveien 30, 5045 Bergen, Norway
>>>>> Roger.Bivand at nhh.no
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Fra: diez.edu at gmail.com [diez.edu at gmail.com] på vegne av Eduardo Diez
>>>>> [
>>>>> eduardodiez at gmx.com]
>>>>>
>>>>> Sendt: 9. oktober 2015 4:04
>>>>>
>>>>> Til: Roger Bivand
>>>>>
>>>>> Kopi: grass-stats at lists.osgeo.org
>>>>>
>>>>> Emne: Re: [GRASS-stats] rgrass7 read/write SpatialPolygonsDataFrames
>>>>> errors
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Roger, with GRASS version 7.0.1 and rgrass7_0.1-3, readVECT now doesn't
>>>>> ask for the driver to use but gives the same error as before:
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> writeVECT(SDF = mz.pol, vname = zm.pol, v.in.ogr_flags = "o")
>>>>> Over-riding projection check
>>>>> Check if OGR layer <q7468fb> contains polygons...
>>>>>
>>>>>
>>>>>
>>>>> Importing 205 features (OGR layer <q7468fb>)...
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>> Registering primitives...
>>>>> 3924 primitives registered
>>>>> 12866 vertices registered
>>>>> Number of nodes: 2884
>>>>> Number of primitives: 3924
>>>>> Number of points: 0
>>>>> Number of lines: 0
>>>>> Number of boundaries: 3924
>>>>> Number of centroids: 0
>>>>> Number of areas: -
>>>>> Number of isles: -
>>>>> -----------------------------------------------------
>>>>> Cleaning polygons
>>>>> -----------------------------------------------------
>>>>> Breaking polygons...
>>>>> Breaking polygons (pass 1: select break points)...
>>>>>
>>>>>
>>>>>
>>>>> Breaking polygons (pass 2: break at selected points)...
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>> Removing duplicates...
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>> Breaking boundaries...
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>> Removing duplicates...
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>> Cleaning boundaries at nodes...
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>> Merging boundaries...
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>> Removing dangles...
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>> Building areas...
>>>>>
>>>>>
>>>>>
>>>>> 205 areas built
>>>>> 13 isles built
>>>>> Number of nodes: 2929
>>>>> Number of primitives: 7268
>>>>> Number of points: 0
>>>>> Number of lines: 0
>>>>> Number of boundaries: 7268
>>>>> Number of centroids: 0
>>>>> Number of areas: 205
>>>>> Number of isles: 13
>>>>> -----------------------------------------------------
>>>>> Removing bridges...
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>> Registering primitives...
>>>>> 498 primitives registered
>>>>> 5828 vertices registered
>>>>> Building areas...
>>>>>
>>>>>
>>>>>
>>>>> 205 areas built
>>>>> 13 isles built
>>>>> Attaching islands...
>>>>>
>>>>>
>>>>>
>>>>> Number of nodes: 306
>>>>> Number of primitives: 498
>>>>> Number of points: 0
>>>>> Number of lines: 0
>>>>> Number of boundaries: 498
>>>>> Number of centroids: 0
>>>>> Number of areas: 205
>>>>> Number of isles: 13
>>>>> -----------------------------------------------------
>>>>> Finding centroids for OGR layer <q7468fb>...
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>> Writing centroids...
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>> 205 input polygons
>>>>> Total area: 1.25187E+006 (205 areas)
>>>>> -----------------------------------------------------
>>>>> Copying features...
>>>>>
>>>>>
>>>>>
>>>>> WARNING: Unable to delete file
>>>>>
>>>>>
>>>>>
>>>>> 'C:/Users/ediez1/AppData/Local/Temp/Rtmpqkeay4/file480c1e5e16a7/file480c11cc233b/.tmp/unknown/vector/tmp_8308/coor'
>>>>> Building topology for vector map <q7468fb at file480c11cc233b>...
>>>>> Registering primitives...
>>>>> 703 primitives registered
>>>>> 6033 vertices registered
>>>>> Building areas...
>>>>>
>>>>>
>>>>>
>>>>> 205 areas built
>>>>> 13 isles built
>>>>> Attaching islands...
>>>>>
>>>>>
>>>>>
>>>>> Attaching centroids...
>>>>>
>>>>>
>>>>>
>>>>> Number of nodes: 306
>>>>> Number of primitives: 703
>>>>> Number of points: 0
>>>>> Number of lines: 0
>>>>> Number of boundaries: 498
>>>>> Number of centroids: 205
>>>>> Number of areas: 205
>>>>> Number of isles: 13
>>>>>
>>>>> execGRASS("v.clean", flags = c("overwrite"), input = zm.pol,
>>>>>>
>>>>>> + output = zm.gnrl, tool = "rmarea", threshold = 3000)
>>>>> --------------------------------------------------
>>>>> Tool: Threshold
>>>>> Remove small areas: 3000
>>>>> --------------------------------------------------
>>>>> Copying features...
>>>>>
>>>>>
>>>>>
>>>>> Rebuilding parts of topology...
>>>>> Building topology for vector map <u10f44b at file480c11cc233b>...
>>>>> Registering primitives...
>>>>> 703 primitives registered
>>>>> 6033 vertices registered
>>>>> Building areas...
>>>>>
>>>>>
>>>>>
>>>>> 205 areas built
>>>>> 13 isles built
>>>>> Attaching islands...
>>>>>
>>>>>
>>>>>
>>>>> Attaching centroids...
>>>>>
>>>>>
>>>>>
>>>>> Number of nodes: 306
>>>>> Number of primitives: 703
>>>>> Number of points: 0
>>>>> Number of lines: 0
>>>>> Number of boundaries: 498
>>>>> Number of centroids: 205
>>>>> Number of areas: 205
>>>>> Number of isles: 13
>>>>> --------------------------------------------------
>>>>> Tool: Remove small areas
>>>>>
>>>>>
>>>>>
>>>>> 174 areas of total size 50601.7 removed
>>>>> --------------------------------------------------
>>>>> Rebuilding topology for output vector map...
>>>>> Building topology for vector map <u10f44b at file480c11cc233b>...
>>>>> Registering primitives...
>>>>> 333 primitives registered
>>>>> 4705 vertices registered
>>>>> Building areas...
>>>>>
>>>>>
>>>>>
>>>>> 31 areas built
>>>>> 3 isles built
>>>>> Attaching islands...
>>>>>
>>>>>
>>>>>
>>>>> Attaching centroids...
>>>>>
>>>>>
>>>>>
>>>>> Number of nodes: 274
>>>>> Number of primitives: 333
>>>>> Number of points: 0
>>>>> Number of lines: 0
>>>>> Number of boundaries: 302
>>>>> Number of centroids: 31
>>>>> Number of areas: 31
>>>>> Number of isles: 3
>>>>>
>>>>> zm.fnl <- readVECT(zm.gnrl)
>>>>>>
>>>>>>
>>>>>
>>>>> Error in scan(file, what, nmax, sep, dec, quote, skip, nlines,
>>>>> na.strings, :
>>>>> line 4 did not have 2 elements
>>>>>
>>>>> sessionInfo()
>>>>>>
>>>>>> R version 3.2.2 (2015-08-14)
>>>>> Platform: x86_64-w64-mingw32/x64 (64-bit)
>>>>> Running under: Windows 7 x64 (build 7601) Service Pack 1
>>>>>
>>>>>
>>>>>
>>>>> locale:
>>>>> [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United
>>>>> States.1252 LC_MONETARY=English_United States.1252
>>>>> [4] LC_NUMERIC=C LC_TIME=English_United
>>>>> States.1252
>>>>>
>>>>>
>>>>>
>>>>> attached base packages:
>>>>> [1] stats graphics grDevices utils datasets methods base
>>>>>
>>>>>
>>>>>
>>>>> other attached packages:
>>>>> [1] RPyGeo_0.9-3 RSAGA_0.94-3 plyr_1.8.3 gstat_1.0-26
>>>>> shapefiles_0.7 foreign_0.8-65 maptools_0.8-37 rgdal_1.0-7
>>>>> classInt_0.1-23
>>>>> [10] raster_2.4-20 rgrass7_0.1-3 XML_3.98-1.3 sp_1.2-0
>>>>>
>>>>>
>>>>>
>>>>> loaded via a namespace (and not attached):
>>>>> [1] Rcpp_0.12.1 intervals_0.15.1 lattice_0.20-33 FNN_1.1
>>>>> spacetime_1.1-4 class_7.3-13 zoo_1.7-12 grid_3.2.2
>>>>> [9] e1071_1.6-7 xts_0.9-7 tools_3.2.2
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Thanks
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> 2015-10-08 18:57 GMT-03:00 Roger Bivand
>>>>> <Roger.Bivand at nhh.no>:
>>>>>
>>>>>
>>>>> Sorry this has taken a long time; please try installing the development
>>>>> version available from:
>>>>>
>>>>>
>>>>>
>>>>> http://win-builder.r-project.org/8Y3Ihj451Clx
>>>>>
>>>>>
>>>>>
>>>>> (rgrass7_0.1-3.zip). The problem was the use of "_" instead of " " in
>>>>> the
>>>>> name of the OGR shapefile driver - in OGR and R it is "ESRI Shapefile",
>>>>> in
>>>>> GRASS it is "ESRI_Shapefile". Because the new long-name work-around
>>>>> worked
>>>>> on systems not needing to use this
>>>>> driver, it wasn't properly tested - and on Windows the GRASS
>>>>> standalone
>>>>> has OGR with the SQLite driver, but the Windows rgdal binary does not
>>>>> have
>>>>> SQLite. You may wish to install the R RSQLite package to move data if
>>>>> you
>>>>> use long field/column names.
>>>>>
>>>>>
>>>>>
>>>>> Before installing this development version, and after library(rgrass7),
>>>>> run example(readVECT) in the basic nc location from:
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> https://grass.osgeo.org/sampledata/north_carolina/nc_basic_spm_grass7.zip
>>>>>
>>>>>
>>>>>
>>>>> it should fail. Then restart R, install the development version of
>>>>> rgrass7
>>>>> from the first link, and repeat - it works for me and should resolve
>>>>> the
>>>>> problem.
>>>>>
>>>>>
>>>>>
>>>>> Hope this helps,
>>>>>
>>>>>
>>>>>
>>>>> Roger
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> --
>>>>>
>>>>> Roger Bivand
>>>>>
>>>>> Norwegian School of Economics
>>>>>
>>>>> Helleveien 30, 5045 Bergen, Norway
>>>>>
>>>>> Roger.Bivand at nhh.no
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Fra:
>>>>> grass-stats-bounces at lists.osgeo.org [
>>>>> grass-stats-bounces at lists.osgeo.org
>>>>> ]
>>>>> på vegne av Eduardo Diez [eduardodiez at gmx.com]
>>>>>
>>>>> Sendt: 24. september 2015 9:11
>>>>>
>>>>> Til:
>>>>> grass-stats at lists.osgeo.org
>>>>>
>>>>> Emne: [GRASS-stats] rgrass7 read/write SpatialPolygonsDataFrames errors
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Hi, I'm a user of R mostly for spatial information processing. Besides
>>>>> sp,
>>>>> maptools, etc. i also use functions from GRASS, at first through
>>>>> spgrass6
>>>>> and then rgrass7. A few days ago, I started having trouble reading and
>>>>> writing SpatialPolygonsDataFrames using
>>>>> readVECT and writeVECT. I tried both reinstalling both grass (6.4.4,
>>>>> 7.0.0 and 7.0.1) and the R packages (spgrass6 0.8-6 and rgrass7 0.1-2),
>>>>> and
>>>>> also with different R objects but i kept getting the same errors.
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> The "mult.r" object is a RasterLayer with 4 zones and the objective is
>>>>> to
>>>>> convert it to polygons and "clean" small areas with v.clean rmarea (as
>>>>> i
>>>>> don't currently know a way of doing this in plain R).
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> The first thing that came out strange was that when executing writeVECT
>>>>> it
>>>>> would ask me for the driver to use (in the way of rgdal's writeOGR, it
>>>>> didn't ask for it before). I figured out that object to be written
>>>>> should
>>>>> be projected (not in lat-lon, don't know
>>>>> why). After doing this i could write the vector and execute v.clean.
>>>>> The
>>>>> problem came again when trying to read the cleaned polygons back to R.
>>>>>
>>>>> Below i paste an extract of my R console with the errors and some
>>>>> information.
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Thanks in advance,
>>>>>
>>>>>
>>>>>
>>>>> Eduardo
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> mult.rp <- projectRaster(mult.r, crs = prj.crs, method = "ngb")
>>>>>
>>>>>>
>>>>>>
>>>>> mult.pol <- rasterToPolygons(mult.rp, dissolve = T)
>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>> initGRASS(gisBase = "c:/Program Files (x86)/GRASS GIS 7.0.1", home =
>>>>>
>>>>>>
>>>>>> tempdir(), override = T)
>>>>>
>>>>> gisdbase C:/Users/ediez1/AppData/Local/Temp/RtmpegbyiQ
>>>>>
>>>>> location file2f7c6dc070e6
>>>>>
>>>>> mapset file2f7c30eb3f52
>>>>>
>>>>> rows 1
>>>>>
>>>>> columns 1
>>>>>
>>>>> north 1
>>>>>
>>>>> south 0
>>>>>
>>>>> west 0
>>>>>
>>>>> east 1
>>>>>
>>>>> nsres 1
>>>>>
>>>>> ewres 1
>>>>>
>>>>> projection NA
>>>>>
>>>>> writeVECT(mult.pol, "pol1", v.in.ogr_flags = "o")
>>>>>
>>>>>>
>>>>>>
>>>>> Error: !is.null(driver) is not TRUE
>>>>>
>>>>> writeVECT(mult.pol, "pol1", v.in.ogr_flags = "o", driver = "ESRI
>>>>>
>>>>>>
>>>>>> Shapefile")
>>>>>
>>>>>
>>>>>
>>>>> Over-riding projection check
>>>>>
>>>>> Check if OGR layer <r848659> contains polygons...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Importing 4 features (OGR layer <r848659>)...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Registering primitives...
>>>>>
>>>>> 3005 primitives registered
>>>>>
>>>>> 15587 vertices registered
>>>>>
>>>>> Number of nodes: 2691
>>>>>
>>>>> Number of primitives: 3005
>>>>>
>>>>> Number of points: 0
>>>>>
>>>>> Number of lines: 0
>>>>>
>>>>> Number of boundaries: 3005
>>>>>
>>>>> Number of centroids: 0
>>>>>
>>>>> Number of areas: -
>>>>>
>>>>> Number of isles: -
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Cleaning polygons
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Breaking polygons...
>>>>>
>>>>> Breaking polygons (pass 1: select break points)...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Breaking polygons (pass 2: break at selected points)...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Removing duplicates...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Breaking boundaries...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Removing duplicates...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Cleaning boundaries at nodes...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Merging boundaries...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Removing dangles...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Building areas...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> 95 areas built
>>>>>
>>>>> 8 isles built
>>>>>
>>>>> Number of nodes: 2771
>>>>>
>>>>> Number of primitives: 7307
>>>>>
>>>>> Number of points: 0
>>>>>
>>>>> Number of lines: 0
>>>>>
>>>>> Number of boundaries: 7307
>>>>>
>>>>> Number of centroids: 0
>>>>>
>>>>> Number of areas: 95
>>>>>
>>>>> Number of isles: 8
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Removing bridges...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Registering primitives...
>>>>>
>>>>> 231 primitives registered
>>>>>
>>>>> 7594 vertices registered
>>>>>
>>>>> Building areas...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> 95 areas built
>>>>>
>>>>> 8 isles built
>>>>>
>>>>> Attaching islands...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Number of nodes: 144
>>>>>
>>>>> Number of primitives: 231
>>>>>
>>>>> Number of points: 0
>>>>>
>>>>> Number of lines: 0
>>>>>
>>>>> Number of boundaries: 231
>>>>>
>>>>> Number of centroids: 0
>>>>>
>>>>> Number of areas: 95
>>>>>
>>>>> Number of isles: 8
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Finding centroids for OGR layer <r848659>...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Writing centroids...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> 95 input polygons
>>>>>
>>>>> Total area: 1.28981E+006 (95 areas)
>>>>>
>>>>> -----------------------------------------------------
>>>>>
>>>>> Copying features...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> WARNING: Unable to delete file
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> 'C:/Users/ediez1/AppData/Local/Temp/Rtmpg10lvi/filec2473b0126/filec24129b63ac/.tmp/unknown/vector/tmp_10484/coor'
>>>>>
>>>>> Building topology for vector map <r848659 at filec24129b63ac>...
>>>>>
>>>>> Registering primitives...
>>>>>
>>>>> 326 primitives registered
>>>>>
>>>>> 7689 vertices registered
>>>>>
>>>>> Building areas...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> 95 areas built
>>>>>
>>>>> 8 isles built
>>>>>
>>>>> Attaching islands...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Attaching centroids...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Number of nodes: 144
>>>>>
>>>>> Number of primitives: 326
>>>>>
>>>>> Number of points: 0
>>>>>
>>>>> Number of lines: 0
>>>>>
>>>>> Number of boundaries: 231
>>>>>
>>>>> Number of centroids: 95
>>>>>
>>>>> Number of areas: 95
>>>>>
>>>>> Number of isles: 8
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> execGRASS("v.clean", flags = c("overwrite"), input = "pol1",
>>>>>
>>>>>>
>>>>>>
>>>>> + output = "clnpol", tool = "rmarea", threshold = 3000)
>>>>>
>>>>> --------------------------------------------------
>>>>>
>>>>> Tool: Threshold
>>>>>
>>>>> Remove small areas: 3000
>>>>>
>>>>> --------------------------------------------------
>>>>>
>>>>> Copying features...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Rebuilding parts of topology...
>>>>>
>>>>> Building topology for vector map <k1f9f75 at filec24129b63ac>...
>>>>>
>>>>> Registering primitives...
>>>>>
>>>>> 326 primitives registered
>>>>>
>>>>> 7689 vertices registered
>>>>>
>>>>> Building areas...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> 95 areas built
>>>>>
>>>>> 8 isles built
>>>>>
>>>>> Attaching islands...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Attaching centroids...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Number of nodes: 144
>>>>>
>>>>> Number of primitives: 326
>>>>>
>>>>> Number of points: 0
>>>>>
>>>>> Number of lines: 0
>>>>>
>>>>> Number of boundaries: 231
>>>>>
>>>>> Number of centroids: 95
>>>>>
>>>>> Number of areas: 95
>>>>>
>>>>> Number of isles: 8
>>>>>
>>>>> --------------------------------------------------
>>>>>
>>>>> Tool: Remove small areas
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> 55 areas of total size 29487.1 removed
>>>>>
>>>>> --------------------------------------------------
>>>>>
>>>>> Rebuilding topology for output vector map...
>>>>>
>>>>> Building topology for vector map <k1f9f75 at filec24129b63ac>...
>>>>>
>>>>> Registering primitives...
>>>>>
>>>>> 212 primitives registered
>>>>>
>>>>> 6840 vertices registered
>>>>>
>>>>> Building areas...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> 40 areas built
>>>>>
>>>>> 4 isles built
>>>>>
>>>>> Attaching islands...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Attaching centroids...
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Number of nodes: 136
>>>>>
>>>>> Number of primitives: 212
>>>>>
>>>>> Number of points: 0
>>>>>
>>>>> Number of lines: 0
>>>>>
>>>>> Number of boundaries: 172
>>>>>
>>>>> Number of centroids: 40
>>>>>
>>>>> Number of areas: 40
>>>>>
>>>>> Number of isles: 4
>>>>>
>>>>> cln.pol <- readVECT("clnpol")
>>>>>
>>>>>>
>>>>>>
>>>>> Error: !is.null(driver) is not TRUE
>>>>>
>>>>> cln.pol <- readVECT("clnpol", driver = "ESRI Shapefile")
>>>>>
>>>>>>
>>>>>>
>>>>> Error in scan(file, what, nmax, sep, dec, quote, skip, nlines,
>>>>> na.strings, :
>>>>>
>>>>> line 4 did not have 2 elements
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> sessionInfo()
>>>>>
>>>>>>
>>>>>>
>>>>> R version 3.2.2 (2015-08-14)
>>>>>
>>>>> Platform: x86_64-w64-mingw32/x64 (64-bit)
>>>>>
>>>>> Running under: Windows 7 x64 (build 7601) Service Pack 1
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> locale:
>>>>>
>>>>> [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United
>>>>> States.1252 LC_MONETARY=English_United States.1252
>>>>>
>>>>> [4] LC_NUMERIC=C LC_TIME=English_United
>>>>> States.1252
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> attached base packages:
>>>>>
>>>>> [1] stats graphics grDevices utils datasets methods base
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> other attached packages:
>>>>>
>>>>> [1] rgrass7_0.1-2 classInt_0.1-22 maptools_0.8-36 ggplot2_1.0.1
>>>>> XML_3.98-1.3 raster_2.4-20 sp_1.2-0
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> loaded via a namespace (and not attached):
>>>>>
>>>>> [1] Rcpp_0.12.1 magrittr_1.5 MASS_7.3-43 munsell_0.4.2
>>>>> colorspace_1.2-6 lattice_0.20-33 stringr_1.0.0 plyr_1.8.3
>>>>>
>>>>> [9] tools_3.2.2 parallel_3.2.2 rgdal_1.0-7 grid_3.2.2
>>>>> gtable_0.1.2 e1071_1.6-7 class_7.3-13 rgeos_0.3-12
>>>>>
>>>>> [17] digest_0.6.8 reshape2_1.4.1 labeling_0.3 stringi_0.5-5
>>>>> scales_0.3.0 foreign_0.8-65 proto_0.3-10
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> grass-stats mailing list
>>>>> grass-stats at lists.osgeo.org
>>>>> http://lists.osgeo.org/mailman/listinfo/grass-stats
>>>>>
>>>>>
>>>>>
>>>> --
>>> Roger Bivand
>>> Department of Economics, Norwegian School of Economics,
>>> Helleveien 30, N-5045 Bergen, Norway.
>>> voice: +47 55 95 93 55; fax +47 55 95 91 00
>>> e-mail: Roger.Bivand at nhh.no
>>>
>>>
>>
> --
> Roger Bivand
> Department of Economics, Norwegian School of Economics,
> Helleveien 30, N-5045 Bergen, Norway.
> voice: +47 55 95 93 55; fax +47 55 95 91 00
> e-mail: Roger.Bivand at nhh.no
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-stats/attachments/20151010/1320da9e/attachment-0001.html>
More information about the grass-stats
mailing list