<div dir="ltr">On Tue, Dec 3, 2013 at 5:23 AM, James David Smith <span dir="ltr"><<a href="mailto:james.david.smith@gmail.com" target="_blank">james.david.smith@gmail.com</a>></span> wrote:<br><div class="gmail_extra">
<div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">Hi all,<div><br></div>
<div>Apologies for cross-posting this to the RPostrgeSQL and PostGIS mailing lists, but I'm not sure where the problem lies.</div><div><br></div><div>I am using the RPostgreSQL package of R to connect to my installation of PostgreSQL/PostGIS. I want to retrieve some coordinates from my database and run the below command. Is this just that R doesn't know what to do with Geometry's, or is it something I should be more concerned about? I ask, as I'm having some wider issues with my workflow and this may be related.</div>


<div><br></div><div>Thanks</div><div><br></div><div>James</div><div><br></div><div><pre style="line-height:15px;font-size:14px;background-color:rgb(225,226,229);margin-bottom:0px;font-family:Consolas,'Lucida Console',monospace;margin-top:0px;border:none;outline:none;white-space:pre-wrap!important">
<span style="white-space:pre-wrap;color:blue">> </span><span style="color:blue">dbGetQuery(con, "SELECT st_transform(st_setsrid(st_makepoint(soose::numeric, soosn::numeric),27700),4326) FROM stage WHERE ssid = '707187161010202'")
</span>                        </pre><pre style="line-height:15px;font-size:14px;background-color:rgb(225,226,229);margin-bottom:0px;font-family:Consolas,'Lucida Console',monospace;margin-top:0px;border:none;outline:none;white-space:pre-wrap!important">
st_transform
1 0101000020E6100000F71184EF3961B3BF830420AF33C64940
<span style="color:rgb(197,6,11)">Warning message:
</span><span style="color:rgb(197,6,11)">In postgresqlExecStatement(conn, statement, ...) :
  RS-DBI driver warning: (unrecognized PostgreSQL field type geometry (id:480912) in column 0)</span></pre><pre style="line-height:15px;font-size:14px;background-color:rgb(225,226,229);margin-bottom:0px;font-family:Consolas,'Lucida Console',monospace;margin-top:0px;border:none;outline:none;white-space:pre-wrap!important">
<span style="color:rgb(197,6,11)"><br></span></pre></div></div></blockquote><div class="gmail_quote"><br></div>James,<div><br></div><div>R knows how to handle geometries using the sp package, but can't read the data in using RPostgreSQL. You would need to use readOGR() in the rgdal package. But what data do you actually need in R, and what do you intend to do with it? If you're not going to be mapping it or doing spatial statistics on it, there may be a way to get the data you want using RPostgreSQL.</div>
<div><br></div><div>Best,</div><div>--Lee </div></div><div><br></div>-- <br>Lee Hachadoorian<br>Asst Professor of Geography, Dartmouth College<br><a href="http://freecity.commons.gc.cuny.edu/" target="_blank">http://freecity.commons.gc.cuny.edu/</a>
</div></div>