Dylan Beaudette wrote:
> Hi,
>
> Here is another approach, using something like a 'concave hull' -- i.e. alpha
> shapes:
>
> http://casoilresource.lawr.ucdavis.edu/drupal/node/919
>
Thanks Dylan and thanks Roman for your advice. (convexhull.xy was what I
wanted).
signature.asc
Description: O
to above and uploaded to this URI
http://sociologi.cjb.net/temp/to-polygons-help.RData, so
print(load(url("http://sociologi.cjb.net/temp/to-polygons-help.RData";)))
plot(sdn)
sapply(sorted.list.of.clusters, function(x)
{plot(clust.geo.test[which(my.cluster == x),], add = T)})
On Wed, Apr 21, 2010 at 03:14:46PM +0200, Roger Bivand wrote:
> On Wed, 21 Apr 2010, Hans Ekbrand wrote:
[...]
>> Well, hclust was useful, once I understood how cutree works. What
>> would be the benefit of dnearneigh(), is it faster?
>>
>
> For larger data sets
On Wed, Apr 21, 2010 at 01:42:05PM +0200, Roger Bivand wrote:
> So you do not want hclust at all, really. Look at dnearneigh() in spdep,
> setting a 100m bound. Then use n.comp.nb() to see which points belong to
> which graph component, using perhaps plot.nb with colours to distinguish
> the
On Wed, Apr 21, 2010 at 09:59:51AM +0200, Hans Ekbrand wrote:
[...]
> head(clust.geo.test.tree$height, 70)
> [1] 0.00 0.00 0.00 0.00 0.00 0.00
> 0.00 0.00 0.00 0.00
> [11] 0.00 0.00 0.00 0.0
On Tue, Apr 20, 2010 at 11:13:22PM +0200, Hans Ekbrand wrote:
> Roger Bivand wrote:
> > On Tue, 20 Apr 2010, Hans Ekbrand wrote:
> >
> >> I have just read about clustering on wikipedia, and learnt that what I
> >> want is:
> >>
> >> Agglomera
Roger Bivand wrote:
> On Tue, 20 Apr 2010, Hans Ekbrand wrote:
>
>> I have just read about clustering on wikipedia, and learnt that what I
>> want is:
>>
>> Agglomerative hierarchical clustering, with complete linkage
>
> library(cluster)
> ?hclust
>
&g
" +proj=utm +zone=33 +ellps=GRS80 +units=m +no_defs"
head(fi...@coords, 2)
East North
[1,] 315359.9 6393110
[2,] 325862.4 6405239
dist(head(fi...@coords, 2))
1
2 16044.10
Is this in meters?
TIA
--
Hans Ekbrand (http://sociologi.cjb.net)
signature.asc
this should be fine:
>
> plot(yta)
> plot(vagar, add = TRUE)
> sdn.reproj <- spTransform(sdn, CRS(proj4string(vagar)))
Michael, you just made my day!
Thanks!
--
Hans Ekbrand (http://sociologi.cjb.net)
A. Because it breaks the logical sequence of discussion
Q. Why is top p
Hi list,
I'm a newcomer to GIS and r-sig-geo. I have exported some shapefiles
from ArcGIS and succeeded in plotting them with R. But one shapefile has
a different coordinate system, I think, and it does not turn up in the
graph when added to the others. In ArcGIS, this layer had an "unknown"
coord
10 matches
Mail list logo
