<div>Dear Prof .Roger,</div>
<div>&nbsp; There is still a little problem with your method, the following are the programs:</div>
<p>a&lt;-readShapePoly(&quot;d:/deleting/a/guichi2.shp&quot;)&nbsp; <br>W &lt;-as(as(a, &quot;SpatialPolygons&quot;), &quot;owin&quot;)&nbsp; # polygons window<br>plot(W)</p>
<p>cases&lt;-coordinates(readShapePoints(&#39;d:/deleting/a/cases.shp&#39;))&nbsp; #points<br>cases&lt;-data.frame(cases)<br>names(cases)&lt;-c(&quot;x&quot;,&quot;y&quot;)<br>cases[1:2,]</p>
<p>plot(W)<br>points(cases$x,cases$y)&nbsp; #no errors,83 points were inside the&nbsp;window</p>
<p>W&lt;-as.owin(W)<br>points&lt;-ppp(cases$x,cases$y,win=W)&nbsp; <font color="#ff0000"><strong>#????????????<br></strong></font><strong><font color="#ff0000">#error</font></strong>:83 points were rejected as lying outside the specified window in: ppp(cases$x, cases$y, win = W) 
<br><strong><font color="#000099">In fact,</font></strong> 83 points were inside the specified window,what&#39;s wrong?</p>
<div>(This&nbsp;packages:<br>spatstat &nbsp; &nbsp; mgcv maptools &nbsp;foreign &nbsp; &nbsp; &nbsp; sp<br>&quot;1.11-6&quot; &quot;1.3-24&quot; &quot;0.6-12&quot; &quot;0.8-20&quot; &quot;0.9-13&quot;)<br>&nbsp;BTW, i put the data in the attachment so that you can check it.
</div>
<div>Thanks.<br>&nbsp;</div>
<div><span class="gmail_quote">On 6/1/07, <b class="gmail_sendername">Roger Bivand</b> &lt;<a href="mailto:Roger.Bivand@nhh.no">Roger.Bivand@nhh.no</a>&gt; wrote:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">On Fri, 1 Jun 2007, zhijie zhang wrote:<br><br>&gt; Dear friends,<br>&gt;&nbsp;&nbsp; I have a polygon, which is the shapefile format. I want to use it as the
<br>&gt; window of points in ppp(x,y,window) of spatstat package, how should i do?<br>&gt; The following programs can&#39;t work,:<br>&gt; a&lt;-readShapePoly(&quot;d:/deleting/a/a.shp&quot;)&nbsp;&nbsp;# read the boundary data<br>
&gt; a2&lt;-Polygon(coordinates(a), hole=as.logical(NA))<br>&gt; ppp(x,y,a2)<br><br>Using the shapefile shipped with maptools:<br><br>library(maptools)<br>library(spatstat)<br>xx &lt;- readShapePoly(system.file(&quot;shapes/sids.shp&quot;, package=&quot;maptools&quot;)[1])
<br>xx1 &lt;- xx[1,] # and taking just the first county<br>plot(xx1)<br>W &lt;- as(as(xx1, &quot;SpatialPolygons&quot;), &quot;owin&quot;)<br>plot(W)<br>points(rpoispp(1000, win=W))<br><br>(This for:<br><br>other attached packages:
<br>spatstat&nbsp;&nbsp;&nbsp;&nbsp; mgcv maptools&nbsp;&nbsp;foreign&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sp<br>&quot;1.11-5&quot; &quot;1.3-24&quot; &quot;0.6-12&quot; &quot;0.8-20&quot; &quot;0.9-13&quot;)<br><br>Hope this helps,<br><br>Roger<br><br><br>&gt;&nbsp;&nbsp;THANKS.<br>&gt;<br>
<br>--<br>Roger Bivand<br>Economic Geography Section, Department of Economics, Norwegian School of<br>Economics and Business Administration, Helleveien 30, N-5045 Bergen,<br>Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
<br>e-mail: <a href="mailto:Roger.Bivand@nhh.no">Roger.Bivand@nhh.no</a><br><br></blockquote></div><br><br clear="all"><br>-- <br>With Kind Regards,<br><br>oooO:::::::::<br>(..):::::::::&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<br>:\.(:::Oooo::&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>::\_)::(..)::&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>:::::::)./:::&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>::::::(_/::::&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>:::::::::::::<br>[***********************************************************************]
<br>Zhi Jie,Zhang ,PHD&nbsp;&nbsp;<br>Tel:86-21-54237149&nbsp;&nbsp;&nbsp;&nbsp;<br>Dept. of Epidemiology,School of Public Health,Fudan University&nbsp;&nbsp; <br>Address:No. 138 Yi Xue Yuan Road,Shanghai,China&nbsp;&nbsp;<br>Postcode:200032&nbsp;&nbsp; <br><a href="mailto:Email:epistat@gmail.com">
Email:epistat@gmail.com</a> <br>Website: <a href="http://www.statABC.com">www.statABC.com</a><br>[***********************************************************************]<br>oooO:::::::::<br>(..):::::::::&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<br>:\.(:::Oooo::&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>::\_)::(..)::&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>:::::::)./:::&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>::::::(_/::::&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>:::::::::::::