Skip to content

Extending sp's SpatialPoints with altitude information

6 messages · Eamon Caddigan, Michael Sumner, Chris English

Eamon Caddigan
# 
Hi all,

I am looking to extend the SpatialPoints and SpatialPointsDataFrame classes
to include altitude information. Somebody recommended I check in with this
group for pointers before getting too far along on this.

I've looked around and haven't seen any other efforts to tackle this
specific problem, although I have seen the sp vignette about extending
classes. Am I missing something? Is this a bad idea?

Thanks,
Eamon Caddigan
Michael Sumner
# 
On Fri, 15 Apr 2016 at 10:36 Eamon Caddigan <eamon.caddigan at gmail.com>
wrote:
It's possible for Points in sp, but not for polygons or lines.

library(sp)
x <- matrix(rnorm(27), ncol = 3)
SpatialPoints(x)

But, the general support outside this in sp is very limited once you break
out of the plane.

I'm working on tools to make this much easier, and to translate from sp to
other forms and back. If you'd like to expand on what you're doing I'm
happy to help. It's possible to do this easily for lines and polygons, and
translate to rgl and ggplot2 and other forms, but it's best to use
relational tables IMO and that gets you outside the playing field pretty
quickly.

See here for some early thoughts, it's still very much in-progress and
design:

http://mdsumner.github.io/2015/12/28/gis3d.html

http://mdsumner.github.io/2016/03/03/polygons-R.html

Cheers, Mike.
Eamon Caddigan
#
On Thu, Apr 14, 2016 at 9:30 PM, Michael Sumner <mdsumner at gmail.com> wrote:
Thanks Mike!

I had no idea that coordinates() handled three-dimensional matrices
already. In that case, it looks like I don't have to do anything at all
(for now).
I look forward to seeing what comes out of this!

-Eamon
Chris English
# 
Mike,

Reading this discussion lead me to finally visiting your github and reading
gris read.md. And of course installing it. I understand mostly where you
are going in the Needs section and were I to recommend a new, fast friend,
it would be Sandro Santilli, strk at keybit.net (http://strk.keybit.net), who
agonized over (I presume) and implemented topology in PostGIS. While this
is not R, the decisions and rigor are informative. Perhaps you are fast
friends already.

I really appreciate allowing mixed objects in gris that seems informed by
spatstat rather more than sp and friends where one feels shackled to a
given representation of line or point or polygon & etc..  In an analysis of
eyetracking data which was initially informed by a mapping perspective (sp,
trajectories, trips), I have essentially had to rewrite sp to to dispense
with sp and many constraints imposed. Those constraints were theoretically
valid but constraining none the less when trying to implement some system
of proposed algorithms that took no heed of sp's object constraints.

When the eyetracking thing is finally put to bed and I get an rgl
compatible laptop that is not my wife's laptop (as things currently stand),
and when I'm a tad better at r, I'd be happy to help with the evolution of
this 3d.

Chris

On Fri, Apr 15, 2016 at 10:06 PM, Eamon Caddigan <eamon.caddigan at gmail.com>
wrote:

  
  


  


      

    

    

      
      

        


  

        

      Michael Sumner
    

    

      
      #
    

  


  

      
On Sat, 16 Apr 2016 at 07:04 chris english <englishchristophera at gmail.com>
wrote:

            

      
      
      
    

        
Thanks for the link, I know that PostGIS has some of this but I haven't
explored it. I've mostly gone on what is in "PostGIS in Action" - I need to
get familiar with that db - I just don't see it as "the solution" for R -
we can have a light-weight toolkit in R itself, and leverage heavy external
libraries only when desired.

            

      
      
      
    

        
Mostly I'm informed by the use of www.Manifold.net and the now-defunct
Eonfusion.  MapInfo .tab and .mif formats also allow mixed topology, but I
never used MapInfo. Eonfusion was so radical it's very hard to describe
briefly, but essentially all storage is via relational tables - vertices,
primitives, and objects. All tables can have any attributes (or none) and
all operations have default choices of which attributes are used (i.e. x,
y, z, time or lon, lat, z, time - or whatever you want to call them). It's
very natural once you are familiar with it. Attribute flexibility on
vertices means  you can manually set whatever you want. I.e. calculate then
triangulate in theta/phi, then process which triangles are "on top", and
voila you have a viewshed analysis. Manifold is stuck in 2D for now but the
editing tools, use of selections, natural use of layers not tied to files,
constrained triangulations, and general slickness and price put it way
ahead of other options IMO even before the new Radian engine is released.

            

      
      
      
    

        
That sounds great, I'm keen to find out what you are doing. That is another
very good point in the "spatial is not special" spectrum. "Spatial" is
really everything, maps are just one kind of "projection", long-lat, x-y-z,
and time are inherent to many things but not everything. Topology exists in
all kinds of  applications. Date time metadata on numeric is a "projection"
for example.

I'd like to see a really general framework where geometry and topology are
the basis and things like sp (and now sfr) can be built on it. I see it as
inevitable that dplyr or its successor is where that's going to happen -
seamless back-ending with a data base is just one reason -  and with ggvis
to replace ggplot2 it will be the go-to tool for visualization and
user-interaction with spatial data. Hopefully the ongoing modernization of
PROJ.4 will also assist in that being just a choice in a
"projection-engine" family.

At the moment I'm concentrating on  https://github.com/mdsumner/spbabel
 which shows some of the ways to nest spatial data in data frames  (via
tidyr) , in a two-level way analogous to sp objects and an "inside-out" way
that's a better match for gris (nested data frames and the ggplot2::fortify
approach can't do this but it's a requirement for topological structures
including triangulation).  It's all in-dev and unstable for now.

Finally, all of my work in this area relies on the tools in dplyr - for
table manipulations - and RTriangle - for constrained triangulation.

Cheers, Mike.



Chris

            

      
      
      
    

        

      
      
    

            

      
      
      
    

        
Dr. Michael Sumner
Software and Database Engineer
Australian Antarctic Division
203 Channel Highway
Kingston Tasmania 7050 Australia

  
  


  


      

    

    

      
      

        


  

        

      Chris English
    

    

      
      #
    

  


  

      
Eonfusion sounds interestingly reminiscent of Ted Nelson's Project Xanadu,
though finding a nearly working version is a very rare bird in the wild.

Chris

        
On Sat, Apr 16, 2016 at 1:15 AM, Michael Sumner <mdsumner at gmail.com> wrote: