R/GEOmap.list.R

In GEOmap: Topographic and Geologic Mapping

GEOmap.list<-function(MAP, SEL=1)
  {
########   return a GEOmap with the points as a list of strokes
    ########  the other parts of the map are unchanged
    if(missing(SEL)) SEL = 1:length(MAP$STROKES$index)

    J = list()
    NEWMAP = list(STROKES = list(nam = NULL, num = NULL, index = NULL,
                    col = NULL, style = NULL, code = NULL, LAT1 = NULL, LAT2 = NULL,
                    LON1 = NULL, LON2 = NULL), POINTS = list(lat = NULL,
                                                 lon = NULL))

    Kmap = 1
    index1 = 0

    lon = NULL
    lat = NULL

    for(i in SEL)
      {

        j1 = MAP$STROKES$index[i] + 1
        j2 = j1 + MAP$STROKES$num[i] - 1
        JEC = j1:j2
        lon = c(MAP$POINTS$lon[JEC])
        lat = c( MAP$POINTS$lat[JEC])

        J[[Kmap]] = list(lat=lat, lon=lon)
        
        NEWMAP$STROKES$nam[Kmap] =  MAP$STROKES$nam[i]
        NEWMAP$STROKES$num[Kmap] = length(lon)
        NEWMAP$STROKES$index[Kmap] = index1
        NEWMAP$STROKES$col[Kmap] =  MAP$STROKES$col[i] 
        NEWMAP$STROKES$style[Kmap] = MAP$STROKES$style[i] 
        NEWMAP$STROKES$code[Kmap] = MAP$STROKES$code[i] 
        NEWMAP$STROKES$LAT1[Kmap] = MAP$STROKES$LAT1[i] 
        NEWMAP$STROKES$LAT2[Kmap] = MAP$STROKES$LAT2[i] 
        NEWMAP$STROKES$LON1[Kmap] = MAP$STROKES$LON1[i] 
        NEWMAP$STROKES$LON2[Kmap] = MAP$STROKES$LON2[i] 
        ## NEWMAP$POINTS$lon = c(NEWMAP$POINTS$lon,  lon)
       ##   NEWMAP$POINTS$lat = c(NEWMAP$POINTS$lat,   lat)
        index1 = index1+length(lon)
        Kmap = Kmap + 1

        
      }

    NEWMAP$LL = J

    return(NEWMAP)
  }