<?xml version="1.0" encoding="utf-8"?><rss version="2.0" xmlns:content="http://purl.org/rss/1.0/modules/content/" xmlns:dc="http://purl.org/dc/elements/1.1/">
<channel>
<title>Das GEO-Forum - R 7-Parameter Helmert Trasfromation</title>
<link>http://forum.diegeodaeten.de/</link>
<description>DieGeodaeten.de ist ein geodätisches Portal, welches von Vermessungsingenieuren der HS Neubrandenburg ins Leben gerufen wurde. Neben Neuigkeiten aus den Bereichen Geodäsie und Geoinformatik werden Buchempfehlungen oder Downloads angeboten.</description>
<language>de</language>
<item>
<title>R 7-Parameter Helmert Trasfromation</title>
<content:encoded><![CDATA[<p>Hallo,</p>
<p>Ich würde gerne in R eine Transformation von Bessel und Lambert nach longlat und WGS84 machen. Dabei habe ich folgenden Cod:</p>
<p>(Mit dem SP library)</p>
<pre><code>
sp_pts1$long&lt;-as.numeric(lt[[n]][[1]]$x)
sp_pts1$lat&lt;-as.numeric(lt[[n]][[1]]$y)
       
coordinates(sp_pts1) &lt;- c(&quot;long&quot;, &quot;lat&quot;)

wkt&lt;-'PROJCS[&quot;Lambert_Austria&quot;,
              GEOGCS[&quot;GCS_Bessel_1841&quot;,
                     DATUM[&quot;D_Bessel_1841&quot;,
                           SPHEROID[&quot;Bessel_1841&quot;,6377397.155,299.1528128]],
                     PRIMEM[&quot;Greenwich&quot;,0],
                     UNIT[&quot;Degree&quot;,0.0174532925199433]],
              PROJECTION[&quot;Lambert_Conformal_Conic&quot;],
              PARAMETER[&quot;False_Easting&quot;,400000],
              PARAMETER[&quot;False_Northing&quot;,400000],
              PARAMETER[&quot;Central_Meridian&quot;,13.3333333333333],
              PARAMETER[&quot;Standard_Parallel_1&quot;,46],
              PARAMETER[&quot;Standard_Parallel_2&quot;,49],
              PARAMETER[&quot;Central_Parallel&quot;,48],
              UNIT[&quot;Meter&quot;,1]]'
           
proj4string(sp_pts1)&lt;-CRS(SRS_string=wkt) 
      
sp_pts1_cart&lt;-spTransform(sp_pts1, CRS(&quot;+proj=cart +ellps=intl&quot;))
sp_pts1_helmert&lt;-spTransform(sp_pts1_cart,crs(&quot;+proj=helmert +convention=coordinate_frame
                                                       +x=-586 +y=-89  +z=-468 
                                                       +rx=5.1 +ry=1.4 +rz=5.4 
                                                       +s=-1.1&quot;))
sp_pts1_cart2&lt;-spTransform(sp_pts1_cart,crs(&quot;+proj=cart +ellps=WGS84&quot;))
sp_pts1_longlat&lt;-spTransform(sp_pts1_cart2,crs(&quot;+proj=longlat +datum=WGS84&quot;)) 
</code></pre><p>Nur er geht nicht und zeigt mir immer folgenden Fehler an:</p>
<pre><code>
ERROR 1: PROJ: proj_crs_get_coordinate_system: Object is not a SingleCRS
</code></pre><p>Ich bin schon am verzweifeln <img src="http://forum.diegeodaeten.de/images/smilies/confused.gif" alt=":confused:" title="hÃ¤?!" /> Habt ihr eine Idee?</p>
]]></content:encoded>
<link>http://forum.diegeodaeten.de/index.php?id=6659</link>
<guid>http://forum.diegeodaeten.de/index.php?id=6659</guid>
<pubDate>Wed, 30 Sep 2020 23:48:17 +0000</pubDate>
<category>Geodäsie/Vermessung</category><dc:creator>Benjamin</dc:creator>
</item>
</channel>
</rss>
