<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">Hello<br>
<br>
There is a very small change betwwen 2154 and 102110 <span
class="Apple-style-span" style="border-collapse: separate;
color: rgb(0, 0, 0); font-family: 'Lucida Console'; font-size:
13px; font-style: normal; font-variant: normal; font-weight:
normal; letter-spacing: normal; line-height: 15px; orphans: 2;
text-align: -webkit-left; text-indent: 0px; text-transform:
none; white-space: pre-wrap; widows: 2; word-spacing: 0px;
-webkit-border-horizontal-spacing: 0px;
-webkit-border-vertical-spacing: 0px;
-webkit-text-decorations-in-effect: none;
-webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;
background-color: rgb(225, 226, 229); ">
<pre tabindex="0" class="GEWYW5YBFEB" id="rstudio_console_output" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-all; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; ">+towgs84=0,0,0,0,0,0,0
</pre>
</span><br>
If I use the small example:<br>
<br>
library(sp)<br>
library(rgdal)<br>
<br>
xlla=700000<br>
ylla=6500000<br>
ncolsa=1000<br>
nrowsa=1000<br>
cellsizea=1<br>
<br>
tableau <- data.frame(IDENT = 1, ncols = 1000, row.names = 1)<br>
<br>
crds <- cbind(x=c(xlla,xlla+ncolsa*cellsizea,
xlla+ncolsa*cellsizea, xlla, xlla), y=c(ylla, ylla,
ylla+nrowsa*cellsizea, ylla+nrowsa*cellsizea, ylla))<br>
<br>
Pa <- Polygon (crds)<br>
Psa <- Polygons (list ( Pa ), 1)<br>
pola <- SpatialPolygons (list ( Psa ),
proj4string=CRS("+init=epsg:2154"))<br>
SPDFa <- SpatialPolygonsDataFrame(pola, tableau)<br>
<br>
print(SPDFa)<br>
<br>
writeOGR(SPDFa, dsn="C:\\ajeter",layer="toto",driver="ESRI
Shapefile",overwrite_layer=TRUE)<br>
<br>
xxx <- readOGR(dsn="C:\\ajeter",layer="toto")<br>
print(xxx)<br>
<br>
For SPDFa: there is this projection:<br>
<span class="Apple-style-span" style="border-collapse: separate;
color: rgb(0, 0, 0); font-family: 'Lucida Console'; font-size:
13px; font-style: normal; font-variant: normal; font-weight:
normal; letter-spacing: normal; line-height: 15px; orphans: 2;
text-align: -webkit-left; text-indent: 0px; text-transform:
none; white-space: pre-wrap; widows: 2; word-spacing: 0px;
-webkit-border-horizontal-spacing: 0px;
-webkit-border-vertical-spacing: 0px;
-webkit-text-decorations-in-effect: none;
-webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;
background-color: rgb(225, 226, 229); ">
<pre tabindex="0" class="GEWYW5YBFEB" id="rstudio_console_output" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-all; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; ">Slot "proj4string":
CRS arguments:
+init=epsg:2154 +proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
And for xxx there is this projection:
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Console'; font-size: 13px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: 15px; orphans: 2; text-align: -webkit-left; text-indent: 0px; text-transform: none; white-space: pre-wrap; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; background-color: rgb(225, 226, 229); "><pre tabindex="0" class="GEWYW5YBFEB" id="rstudio_console_output" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-a
ll; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; ">Slot "proj4string":
CRS arguments:
+proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +units=m +no_defs
I don't see the same argument in the CRS, and in 2154, there is this small part:
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Console'; font-size: 13px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: 15px; orphans: 2; text-align: -webkit-left; text-indent: 0px; text-transform: none; white-space: pre-wrap; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; background-color: rgb(225, 226, 229); "><pre tabindex="0" class="GEWYW5YBFEB" id="rstudio_console_output" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-a
ll; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; ">+towgs84=0,0,0,0,0,0,0
I think there is no big problem of shift but when I create a file and I put it in other GIS (Qgis), it is not recognize as 2154.
Last comment, with the shapefile, the prj don't give the same things betwwen 2154 and 102110
2154 PROJCS["Lambert_Conformal_Conic",GEOGCS["GCS_GRS 1980(IUGG, 1980)",DATUM["D_unknown",SPHEROID["GRS80",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["standard_parallel_1",49],PARAMETER["standard_parallel_2",44],PARAMETER["latitude_of_origin",46.5],PARAMETER["central_meridian",3],PARAMETER["false_easting",700000],PARAMETER["false_northing",6600000],UNIT["Meter",1]]
102110 PROJCS["RGF93_Lambert_93",GEOGCS["GCS_RGF93",DATUM["D_RGF_1993",SPHEROID["GRS_1980",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["standard_parallel_1",49],PARAMETER["standard_parallel_2",44],PARAMETER["latitude_of_origin",46.5],PARAMETER["central_meridian",3],PARAMETER["false_easting",700000],PARAMETER["false_northing",6600000],UNIT["Meter",1]]
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Console'; font-size: 13px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: 15px; orphans: 2; text-align: -webkit-left; text-indent: 0px; text-transform: none; white-space: pre-wrap; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; background-color: rgb(225, 226, 229); "><pre tabindex="0" class="GEWYW5YBFEB" id="rstudio_console_output" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-a
ll; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Console'; font-size: 13px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: 15px; orphans: 2; text-align: -webkit-left; text-indent: 0px; text-transform: none; white-space: pre-wrap; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; background-color: rgb(225, 226, 229); "><pre tabindex="0" class="GEWYW5YBFEB" id="rstudio_console_output" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; outline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none;
border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-all; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Console'; font-size: 13px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: 15px; orphans: 2; text-align: -webkit-left; text-indent: 0px; text-transform: none; white-space: pre-wrap; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; background-color: rgb(225, 226, 229); "><pre tabindex="0" class="GEWYW5YBFEB" id="rstudio_console_output" style="font-family: 'Lucida Console'; font-size: 10pt !important; outline-style: none; o
utline-width: initial; outline-color: initial; border-top-style: none; border-right-style: none; border-bottom-style: none; border-left-style: none; border-width: initial; border-color: initial; white-space: pre-wrap !important; word-break: break-all; margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; -webkit-user-select: text; line-height: 1.2; ">GCS_RGF93 don't appear in R!
</pre></span></pre></span></pre></span>
Thanks for your help
Frédéric
</pre></span></pre></span></pre>
</span><br>
<div class="moz-signature">
<meta http-equiv="content-type" content="text/html;
charset=windows-1252">
<title></title>
<font face="Arial" size="2,5"><b>
Frédéric Pons
</b></font><br>
<div class="moz-signature">
<div class="moz-signature"> <font color="#666666"
face="Arial" size="1"><b>
Expert hydraulique sur les inondations et aléas côtiers
<br>
</b></font><font color="#666666" face="Arial" size="1"><b>
DREC/Service Risques Inondations Littoraux et
Hydraulique
</b></font><font color="#666666" face="Arial" size="1"><b>
- Tél.: (33)4 42 24 76 68
</b></font><br>
<font color="#666666" face="Arial" size="1"><b>
Direction Territoriale Méditerranée<br>
</b></font>
<img alt="" src="cid:part1.07020408.01000805@cerema.fr"
height="45" width="183"><br>
<div class="moz-signature"><font color="#ee7f00"
face="Arial" size="1">Centre d’études et d’expertise sur
les risques, l’environnement, la mobilité et
l’aménagement
</font><br>
<a href="http://www.cerema.fr" class="moz-signature"><font
color="#ee7f00" face="Arial" size="4">www.cerema.fr</font></a><br>
<div class="moz-signature"> <font color="#666666"
face="Arial" size="1"> </font>
</div>
</div>
</div>
</div>
</div>
Le 22/09/2015 16:44, > Hodgess, Erin (par Internet) a écrit :<br>
</div>
<blockquote
cite="mid:FF9DB805FC41CC4E95825A50F68063028C5250FE@challenger.uhd.campus"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
<div style="direction: ltr;font-family: Tahoma;color:
#000000;font-size: 10pt;">Hello!<br>
<br>
I ran through your example, and brought the OGR file back into R
via the readOGR function.<br>
<br>
> xxx <- readOGR(dsn="C:\\ajeter",layer="toto")<br>
OGR data source with driver: ESRI Shapefile <br>
Source: "C:\ajeter", layer: "toto"<br>
with 1 features<br>
It has 2 fields<br>
> str(xxx)<br>
Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5
slots<br>
..@ data :'data.frame': 1 obs. of 2 variables:<br>
.. ..$ IDENT: num 1<br>
.. ..$ ncols: num 1000<br>
..@ polygons :List of 1<br>
.. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots<br>
.. .. .. ..@ Polygons :List of 1<br>
.. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5
slots<br>
.. .. .. .. .. .. ..@ labpt : num [1:2] 700500 6500500<br>
.. .. .. .. .. .. ..@ area : num 1e+06<br>
.. .. .. .. .. .. ..@ hole : logi FALSE<br>
.. .. .. .. .. .. ..@ ringDir: int 1<br>
.. .. .. .. .. .. ..@ coords : num [1:5, 1:2] 700000 700000
701000 701000 700000 ...<br>
.. .. .. ..@ plotOrder: int 1<br>
.. .. .. ..@ labpt : num [1:2] 700500 6500500<br>
.. .. .. ..@ ID : chr "0"<br>
.. .. .. ..@ area : num 1e+06<br>
..@ plotOrder : int 1<br>
..@ bbox : num [1:2, 1:2] 700000 6500000 701000 6501000<br>
.. ..- attr(*, "dimnames")=List of 2<br>
.. .. ..$ : chr [1:2] "x" "y"<br>
.. .. ..$ : chr [1:2] "min" "max"<br>
..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot<br>
.. .. ..@ projargs: chr "+proj=lcc +lat_1=49 +lat_2=44
+lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80
+units=m +no_defs"<br>
> <br>
If you put that big @projagrs into google.com, this is what you
get:<br>
<pre style="color: rgb(0, 0, 0); font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; widows: 1; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; white-space: pre-wrap;">Proj4js.defs["EPSG:2154"] = "+proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs";
Is that what you want, please?
Thanks,
Erin
</pre>
<br>
<div style="font-family: Times New Roman; color: #000000;
font-size: 16px">
<hr tabindex="-1">
<div style="direction: ltr;" id="divRpF137580"><font
color="#000000" face="Tahoma" size="2"><b>From:</b>
R-sig-Geo [<a class="moz-txt-link-abbreviated" href="mailto:r-sig-geo-bounces@r-project.org">r-sig-geo-bounces@r-project.org</a>] on behalf of
PONS Frederic - CEREMA/DTerMed/DREC/SRILH
[<a class="moz-txt-link-abbreviated" href="mailto:Frederic.Pons@cerema.fr">Frederic.Pons@cerema.fr</a>]<br>
<b>Sent:</b> Tuesday, September 22, 2015 9:25 AM<br>
<b>To:</b> <a class="moz-txt-link-abbreviated" href="mailto:r-sig-geo@r-project.org">r-sig-geo@r-project.org</a><br>
<b>Subject:</b> [R-sig-Geo] Problem in EPSG 2154 =>
102110<br>
</font><br>
</div>
<div>
<div class="moz-cite-prefix">Dear R-users,<br>
<br>
I have a problem when I create a spatialpolygondataframe<br>
<br>
I want to have an EPSG 2154 but after the export the epsg
is 102110.<br>
<br>
It seems to be the "same" CRS but it is not what I want.<br>
<br>
A small example to help me<br>
Best regards<br>
<br>
xlla=700000<br>
ylla=6500000<br>
ncolsa=1000<br>
nrowsa=1000<br>
cellsizea=1<br>
<br>
tableau <- data.frame(IDENT = 1, ncols = 1000,
row.names = 1)<br>
<br>
crds <- cbind(x=c(xlla,xlla+ncolsa*cellsizea,
xlla+ncolsa*cellsizea, xlla, xlla), y=c(ylla, ylla,
ylla+nrowsa*cellsizea, ylla+nrowsa*cellsizea, ylla))<br>
<br>
Pa <- Polygon (crds)<br>
Psa <- Polygons (list ( Pa ), 1)<br>
pola <- SpatialPolygons (list ( Psa ),
proj4string=CRS("+init=epsg:2154"))<br>
SPDFa <- SpatialPolygonsDataFrame(pola, tableau)<br>
<br>
writeOGR(SPDFa, dsn="C:\\ajeter",layer="toto",driver="ESRI
Shapefile",overwrite_layer=TRUE)<br>
<br>
<div class="moz-signature"><font face="Arial" size="2,5"><b>Frédéric
Pons </b></font><br>
<div class="moz-signature">
<div class="moz-signature"><font color="#666666"
face="Arial" size="1"><b>Expert hydraulique sur
les inondations et aléas côtiers
<br>
</b></font><font color="#666666" face="Arial"
size="1"><b>DREC/Service Risques Inondations
Littoraux et Hydraulique
</b></font><font color="#666666" face="Arial"
size="1"><b>- Tél.: (33)4 42 24 76 68
</b></font><br>
<font color="#666666" face="Arial" size="1"><b>Direction
Territoriale Méditerranée<br>
</b></font><img alt=""
src="cid:part3.09050607.01080106@cerema.fr"
height="45" width="183"><br>
<div class="moz-signature"><font color="#ee7f00"
face="Arial" size="1">Centre d’études et
d’expertise sur les risques, l’environnement, la
mobilité et l’aménagement
</font><br>
<a moz-do-not-send="true"
href="http://www.cerema.fr"
class="moz-signature" target="_blank"><font
color="#ee7f00" face="Arial" size="4">www.cerema.fr</font></a><br>
</div>
</div>
</div>
</div>
</div>
<br>
</div>
</div>
</div>
</blockquote>
<br>
</body>
</html>