<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
James Cheshire wrote:
<blockquote
cite="mid:76801fef0811250250r3a6b6879j3a911c02d000c40a@mail.gmail.com"
type="cite">Dear adegenet users,<br>
<br>
When I print the Monmonier object class it only lists the coordinates
of the the first three points in each run as follows:<br>
<br>
$call:monmonier(xy = Coordinates, dist = Dist, cn = nbDelaunay, nrun =
20, <br>
scanthres = TRUE)<br>
<br>
# Object content #<br>
Class: monmonier<br>
$nrun (number of successive runs): 20<br>
$run1 ... $run20: runs of the algorithm<br>
$threshold (minimum difference between neighbours): 0.4003<br>
$xy: spatial coordinates<br>
$cn: connection network<br>
<br>
# Runs content #<br>
# Run 1<br>
# First direction<br>
Class: list<br>
$path:<br>
x y<br>
Point_1 361546.4 280420.0<br>
Point_2 376145.2 284840.4<br>
Point_3 387056.9 281490.3<br>
...<br>
<br>
<br>
Is it possible to print or extract all the coordinates used to plot the
Monmonier boundaries? My reason for doing this is so that I can export
them into a Geographical Information System software package. <br>
<br>
Many thanks for your help in advance, it is greatly appreciated.<br>
<br>
James Cheshire<br>
<br>
PhD Student<br>
University College London<br>
<br>
</blockquote>
Hello James, <br>
<br>
yes, it is possible (see ?monmonier, "value" section for a description
of monmonier objects). Here is an example:<br>
<br>
## <br>
<tt>library(adegenet)<br>
load(system.file("files/mondata2.rda",package="adegenet"))<br>
cn2 <- chooseCN(mondata2$xy,type=1,ask=FALSE)<br>
mon2 <- monmonier(mondata2$xy,dist(mondata2$x2),cn2,threshold=2,
nrun=3)<br>
plot(mon2,mondata2$x2,method="greylevel",add.arr=FALSE,bwd=6,col=c("red","blue","green"),csize=.5)<br>
<br>
mon2$run1 # gives the content of first boundary<br>
mon2$run1$dir1 # gives the content of first boundary, 1st direction<br>
mon2$run1$dir1$path # gives the content of first boundary, 1st
direction, coordinates</tt><br>
##<br>
<br>
Another simpler way maybe using Peter Solymos's coords.monmonier (see
?coords.monmonier):<br>
##<br>
<tt>> coords.monmonier(mon2)<br>
<br>
$run1<br>
$run1$dir1<br>
x.hw y.hw first second<br>
Point_1 5.585655 3.353934 16 50<br>
Point_2 6.175788 3.745792 16 37<br>
Point_3 6.050270 4.321076 37 50<br>
Point_4 6.456283 4.801207 36 37<br>
Point_5 7.097451 5.369001 21 36<br>
Point_6 6.246300 5.912191 26 36<br>
Point_7 5.534326 6.166065 26 34<br>
<br>
$run1$dir2<br>
x.hw y.hw first second<br>
Point_1 5.585655 3.353934 16 50<br>
Point_2 5.304805 3.512757 3 50<br>
Point_3 4.769548 3.836670 3 15<br>
[...]</tt><br>
##<br>
<br>
Best, <br>
<br>
Thibaut.<br>
<br>
<div class="moz-signature">-- <br>
<title>Titre du document</title>
<meta http-equiv="Content-Type" content="text/html; ">
<meta name="GENERATOR" content="hevea 1.07">
<!--HEVEA command line is: /usr/bin/hevea signature.tex -->
<!--HTMLHEAD--><!--ENDHTML--><!--PREFIX <ARG ></ARG>--><!--CUT DEF section 1 -->######################################
<br>
Dr Thibaut JOMBART <br>
CNRS UMR 5558 - Laboratoire de Biométrie et Biologie Evolutive <br>
Universite Lyon 1 <br>
43 bd du 11 novembre 1918 <br>
69622 Villeurbanne Cedex <br>
Tél. : 04.72.43.29.35 <br>
Fax : 04.72.43.13.88 <br>
<a class="moz-txt-link-abbreviated" href="mailto:jombart@biomserv.univ-lyon1.fr">jombart@biomserv.univ-lyon1.fr</a> <br>
<a href="http://biomserv.univ-lyon1.fr/%7Ejombart/">http://biomserv.univ-lyon1.fr/%7Ejombart/</a><br>
<a href="http://adegenet.r-forge.r-project.org/">http://adegenet.r-forge.r-project.org/</a><br>
</div>
</body>
</html>