Re: [R-sig-Geo] removing internal boundaries from a spatstat window that was created from three adjacent counties
Adrian, Rolf, et al--
Thanks so much; worked like a charm. Any value of r below 1 (meter, I
assume?) seemed to serve well.
Missing FAQ aside, your book is phenomenal. I'm working through it in
concert with my current project and learning a ton.
--Chris
Christopher W. Ryan, MD, MS
cryanatbinghamtondotedu
https://www.linkedin.com/in/ryancw
Early success is a terrible teacher. Youre essentially being rewarded
for a lack of preparation, so when you find yourself in a situation
where you must prepare, you cant do it. You dont know how.
--Chris Hadfield, An Astronaut's Guide to Life on Earth
Adrian Baddeley wrote:
> Dear Christopher
>
> Try
> newW <- closing(sremsWindow, r)
>
> where 'r' is a very small distance (slightly greater than the largest gap
> that you want to rub out).
>
> See closing.owin
>
> (Dang, that should have been a FAQ in our book)
>
> A
>
> Prof Adrian Baddeley DSc FAA
> Department of Mathematics and Statistics
> Curtin University, Perth, Western Australia
>
>
> From: Rolf Turner <r.tur...@auckland.ac.nz>
> Sent: Tuesday, 16 February 2016 6:03 AM
> To: Christopher W. Ryan
> Cc: r-sig-geo@r-project.org; Adrian Baddeley; Ege Rubak
> Subject: Re: [R-sig-Geo] removing internal boundaries from a spatstat window
> that was created from three adjacent counties
>
> On 13/02/16 11:59, Christopher W. Ryan wrote:
>
>> In spatstat, I want to create a single window from three adjacent
>> counties in New York State, USA. My original data is a shapefile,
>> "cty036" showing all the New York State county boundaries. Here's what
>> I've done so far:
>>
>> == begin code =
>> library(shapefiles)
>> library(spatstat)
>> library(maptools)
>> library(sp)
>> library(RColorBrewer)
>> library(rgdal)
>>
>> nysArea <- readOGR("C:/DATA/GeographicData", "cty036")
>>
>> nysAreaUTM <- spTransform(nysArea, CRS("+proj=utm +zone=18 +ellps=GRS80
>> +units=m +no_defs"))
>>
>> sremsAreaUTM <- subset(nysAreaUTM, NAME=="Broome" | NAME=="Tioga" |
>> NAME=="Chenango")
>>
>> sremsWindow <- as(sremsAreaUTM, "owin")
>> end code =
>>
>>
>> This works pretty well, producing a 3-county owin object. But the
>> dividing lines between the counties are shown in the window, whenever I
>> plot it or plot subsequent point patterns that use the window. In
>> essence, in plots it looks like 3 polygons instead of one big one. I'd
>> prefer not to have the inter-county boundaries be visible--I'd rather
>> have just one big polygon for the whole area. How can I remove them? Or
>> should I create the window differently from the get-go?
>>
>> Thanks.
>
> I think it probable that your owin object has come out as a multiple
> polygon. Look at length(sremsWindow$bdry) --- this is probably 3 (or more).
>
> I would guess that the internal boundaries actually consist of *pairs*
> of closely juxtaposed lines --- which look like single lines when plotted.
>
> Have you read the spatstat vignette "shapefiles"? That might give you
> some guidance as to how to proceed.
>
> I would have thought the automatic repair process that spatstat now uses
> would fix up this problem. What version of spatstat are you using?
>
> You *might* be able to fix things up by doing (something like):
>
> W1 <- owin(poly=sremsWindow$bdry[[1]])
> W2 <- owin(poly=sremsWindow$bdry[[2]])
> W3 <- owin(poly=sremsWindow$bdry[[3]])
>
> W <- union.owin(W1,W2,W3)
>
> But if my guess about the internal boundaries actually being pairs of
> line is correct, this may not work.
>
> It's hard to say without having your sremsWindow object to experiment on.
>
> Or perhaps you need to extract the three counties separately as owin
> objects and the apply union.owin() to the three results. (Again, if my
> guess is correct, this may not work.)
>
> Adrian or Ege may be able to propose a more efficient/effective solution.
>
> cheers,
>
> Rolf Turner
>
> --
> Technical Editor ANZJS
> Department of Statistics
> University of Auckland
> Phone: +64-9-373-7599 ext. 88276
>
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Re: [R-sig-Geo] removing internal boundaries from a spatstat window that was created from three adjacent counties
On 13/02/16 11:59, Christopher W. Ryan wrote:
In spatstat, I want to create a single window from three adjacent
counties in New York State, USA. My original data is a shapefile,
"cty036" showing all the New York State county boundaries. Here's what
I've done so far:
== begin code =
library(shapefiles)
library(spatstat)
library(maptools)
library(sp)
library(RColorBrewer)
library(rgdal)
nysArea <- readOGR("C:/DATA/GeographicData", "cty036")
nysAreaUTM <- spTransform(nysArea, CRS("+proj=utm +zone=18 +ellps=GRS80
+units=m +no_defs"))
sremsAreaUTM <- subset(nysAreaUTM, NAME=="Broome" | NAME=="Tioga" |
NAME=="Chenango")
sremsWindow <- as(sremsAreaUTM, "owin")
end code =
This works pretty well, producing a 3-county owin object. But the
dividing lines between the counties are shown in the window, whenever I
plot it or plot subsequent point patterns that use the window. In
essence, in plots it looks like 3 polygons instead of one big one. I'd
prefer not to have the inter-county boundaries be visible--I'd rather
have just one big polygon for the whole area. How can I remove them? Or
should I create the window differently from the get-go?
Thanks.
I think it probable that your owin object has come out as a multiple
polygon. Look at length(sremsWindow$bdry) --- this is probably 3 (or more).
I would guess that the internal boundaries actually consist of *pairs*
of closely juxtaposed lines --- which look like single lines when plotted.
Have you read the spatstat vignette "shapefiles"? That might give you
some guidance as to how to proceed.
I would have thought the automatic repair process that spatstat now uses
would fix up this problem. What version of spatstat are you using?
You *might* be able to fix things up by doing (something like):
W1 <- owin(poly=sremsWindow$bdry[[1]])
W2 <- owin(poly=sremsWindow$bdry[[2]])
W3 <- owin(poly=sremsWindow$bdry[[3]])
W <- union.owin(W1,W2,W3)
But if my guess about the internal boundaries actually being pairs of
line is correct, this may not work.
It's hard to say without having your sremsWindow object to experiment on.
Or perhaps you need to extract the three counties separately as owin
objects and the apply union.owin() to the three results. (Again, if my
guess is correct, this may not work.)
Adrian or Ege may be able to propose a more efficient/effective solution.
cheers,
Rolf Turner
--
Technical Editor ANZJS
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 88276
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
[R-sig-Geo] removing internal boundaries from a spatstat window that was created from three adjacent counties
In spatstat, I want to create a single window from three adjacent
counties in New York State, USA. My original data is a shapefile,
"cty036" showing all the New York State county boundaries. Here's what
I've done so far:
== begin code =
library(shapefiles)
library(spatstat)
library(maptools)
library(sp)
library(RColorBrewer)
library(rgdal)
nysArea <- readOGR("C:/DATA/GeographicData", "cty036")
nysAreaUTM <- spTransform(nysArea, CRS("+proj=utm +zone=18 +ellps=GRS80
+units=m +no_defs"))
sremsAreaUTM <- subset(nysAreaUTM, NAME=="Broome" | NAME=="Tioga" |
NAME=="Chenango")
sremsWindow <- as(sremsAreaUTM, "owin")
end code =
This works pretty well, producing a 3-county owin object. But the
dividing lines between the counties are shown in the window, whenever I
plot it or plot subsequent point patterns that use the window. In
essence, in plots it looks like 3 polygons instead of one big one. I'd
prefer not to have the inter-county boundaries be visible--I'd rather
have just one big polygon for the whole area. How can I remove them? Or
should I create the window differently from the get-go?
Thanks.
--Chris
--
Christopher W. Ryan, MD, MS
cryanatbinghamtondotedu
https://www.linkedin.com/in/ryancw
Early success is a terrible teacher. You’re essentially being rewarded
for a lack of preparation, so when you find yourself in a situation
where you must prepare, you can’t do it. You don’t know how.
--Chris Hadfield, An Astronaut's Guide to Life on Earth
___
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
