<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=Content-Type content="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.StileMessaggioDiPostaElettronica17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.apple-style-span
{mso-style-name:apple-style-span;}
.MsoChpDefault
{mso-style-type:export-only;}
@page Section1
{size:612.0pt 792.0pt;
margin:70.85pt 2.0cm 2.0cm 2.0cm;}
div.Section1
{page:Section1;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=IT link=blue vlink=purple>
<div class=Section1>
<p class=MsoNormal><span lang=EN-US>Thanks for All help give me.<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>I spoke with my professor and now I must do
some </span><span class=apple-style-span><span lang=EN-US style='font-size:
9.0pt;font-family:"Arial","sans-serif";color:black'>changes. After subtract the
column Z with clomun DEM_SGRD to obtain column H. I wish to create a new shape
with only this H column. I tried and tried but when I write the shape file
there is always an error.<o:p></o:p></span></span></p>
<p class=MsoNormal><span class=apple-style-span><span lang=EN-US
style='font-size:9.0pt;font-family:"Arial","sans-serif";color:black'>Sorry for
all disturb.<o:p></o:p></span></span></p>
<p class=MsoNormal><span class=apple-style-span><span lang=EN-US
style='font-size:9.0pt;font-family:"Arial","sans-serif";color:black'><o:p> </o:p></span></span></p>
<p class=MsoNormal><span class=apple-style-span><span lang=EN-US
style='font-size:9.0pt;font-family:"Arial","sans-serif";color:black'>I don’t
want to be ill</span></span><span lang=EN-US>-mannered, to use improperly your
time, but before to save the shape is It possible to save only all records with
H value > of 2. (the statistical analysis need this clean)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Thanks<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Ale<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>> prova <- readOGR(".",
"prova")<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>OGR data source with driver: ESRI Shapefile
<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Source: ".", layer:
"prova"<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>with 42 rows and 2
columns<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>> summary(prova)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Object of class SpatialPointsDataFrame<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Coordinates:<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
min max<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>coords.x1 267980.1 267998.1<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>coords.x2 4147799.4 4147818.8<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Is projected: NA <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>proj4string : [NA]<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Number of points: 42<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Data attributes:<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
Z
DEM_SGRD <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> Min. :1415
Min. :1412 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 1st Qu.:1419 1st
Qu.:1413 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> Median :1421 Median
:1414 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> Mean :1423
Mean :1414 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 3rd Qu.:1424 3rd
Qu.:1415 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> </span>Max.
:1437 Max. :1417 <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>> prova$H <- prova$Z - prova$DEM_SGRD <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><span lang=EN-US>> summary(prova)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Object of class SpatialPointsDataFrame<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Coordinates:<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
min max<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>coords.x1 267980.1 267998.1<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>coords.x2 4147799.4 4147818.8<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Is projected: NA <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>proj4string : [NA]<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Number of points: 42<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>Data attributes:<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US>
Z
DEM_SGRD
H <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> Min. :1415
Min. :1412 Min. : 3.169 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 1st Qu.:1419 1st
Qu.:1413 1st Qu.: 5.612 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> Median :1421 Median
:1414 Median : 7.320 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> Mean :1423
Mean :1414 Mean : 9.265 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> 3rd Qu.:1424 3rd
Qu.:1415 3rd Qu.:10.982 <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> </span>Max.
:1437 Max. :1417 Max.
:21.806 <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>> prova2 <- prova$H<o:p></o:p></p>
<p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p>
<p class=MsoNormal><span lang=EN-US>> summary(prova2)<o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> Min. 1st Qu.
Median Mean 3rd Qu. Max. <o:p></o:p></span></p>
<p class=MsoNormal><span lang=EN-US> </span>3.169
5.612 7.320 9.265 10.980 21.810 <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>> writeOGR(prova2, ".", "prova2",
driver="ESRI Shapefile")<o:p></o:p></p>
<p class=MsoNormal>Errore in writeOGR(prova2, ".",
"prova2", driver = "ESRI Shapefile") : <o:p></o:p></p>
<p class=MsoNormal> obj of wrong class<o:p></o:p></p>
</div>
</body>
</html>