line - R from SpatialPointsDataFrame to SpatialLines -

-

i'have spatialpointsdataframe load 

pst<-readogr("/data_spatial/coast/","points_coast")

and spatiallines in output, have find somthing 

coord<-as.data.frame(coordinates(pst)) slo1<-line(coord)  sli1<-lines(list(slo1),id="coastline") coastline <- spatiallines(list(sli1)) class(coastline)

it seems work when try plot(coastline) , have line should not there ... a bad coastline

some 1 can me ? the shapefile here !

i have looked @ shapefile. there id column, if plot data, seems id not ordered north-south or something. lines created because point order not perfect, connecting points next each other in table, far each other in terms of space. try figure out correct ordering of data calculating distances between points , ordering on distance.

a workaround remove lines longer distance, e.g. 500 m.. first, find out distance between consecutive coordinates larger distance: breaks. take subset of coordinates between 2 breaks , lastly create lines subset. end coastline consisting of several (breaks-1) segments , without erroneous ones.

# read data library(rgdal) pst<-readogr("/data_spatial/coast/","points_coast") coord<-as.data.frame(coordinates(pst)) colnames(coord) <- c('x','y')  # determine distance between consective coordinates linelength = linelength(as.matrix(coord),sum=f) # 'id' of long lines, plus first , last item of dataset breaks = c(1,which(linelength>500),nrow(coord))  # check position of breaks breaks = c(1,which(linelength>500),nrow(coord))  # plot extent of coords , check breaks plot(coord,type='n') points(coord[breaks,], pch=16,cex=1)  # create vector filled lines of each subset ll <- vector("list", length(breaks)-1) (i in 1: (length(breaks)-1)){     subcoord = coord[(breaks[i]+1):(breaks[i+1]),]  # check if subset contains more 2 coordinates     if (nrow(subcoord) >= 2){         slo1<-line(subcoord)         sli1<-lines(list(slo1),id=paste0('section',i))         ll[[i]] = sli1     }  }  # remove invalid lines nulls = which(unlist(lapply(ll,is.null))) ll = ll[-nulls] lin = spatiallines(ll) # add result plot lines(lin,col=2)  # write shapefile df = data.frame(row.names=names(lin),id=1:length(names(lin))) lin2 = spatiallinesdataframe(sl=lin, data=df) proj4string(lin2) <- proj4string(pst) writeogr(obj=lin2, layer='coastline', dsn='/data_spatial/coast', driver='esri shapefile')

coastline


Comments

Post a Comment

Popular posts from this blog

php - Magento - Deleted Base url key -

-
i working on cleaning url: http://mydomain.com/shop/shop/deals i going through settings in configuration tab in admin panel on magento , change secure , unsecure base urls remove first base url "shop" web page read: http://mydomain.com/shop/deals so changed unsecure base urls mydomain.com/shop mydomain.com , secure 12.345.67.89/shop 12.345.67.89. wrong thing do. i cannot reload page or admin console. how access these files through ftp , change correctly there?? when make changes in admin panel changes made in database not in file. have make changes in database. 1.go database , select core_config_data table under table find web/unsecure/base_url , web/secure/base_url . make changes here. don't forget clear cache in var/cache folder
Read more

javascript - Tooltipster plugin not firing jquery function when button or any click even occur -

-
i want use tooltipster plugin click element in body show tooltip , want have close button inside tooltip close it. however use jquery click() function handle won't fired. tried solution in post. worked when tooltip triggered hover event. tooltipster plugin not firing jquery function when link in in tooltip clicked the original solution using hover in jsfiddle trigger: 'hover', instead using click show tooltip trigger: 'click', http://jsfiddle.net/bcqyl/7/ it won't fire click event in jquery block. can captured normal javascript function only. checked should code in tooltipster plugin locked , captured "click()" event inside tips when in click trigger mode. i tried other event onchange event of radio button fired when using code in original solution @evan right. tooltip doesn't exist in dom , can't bind/delegate action (even if go way $(document) after domready). if you're after close button inside of too...
Read more

java - WrongTypeOfReturnValue exception thrown when unit testing using mockito -

-
my test list<person> mylist; @test public void testisvalidperson() { mylist = new arraylist<person>(); mylist.add(new person("tom")); when(persondao.get(person)).thenreturn(mylist); when((persondao.get(person)).isempty()).thenreturn(false);//------exception thrown boolean result = service.isvalid("tom"); assertfalse(result); } method tested: public boolean isvalid(string person){ persondao = new persondao(); person personobj = new person(person); return (persondao.get(person).isempty())?false : true; } exception thrown: org.mockito.exceptions.misusing.wrongtypeofreturnvalue: boolean cannot returned get() get() should return list *** if you're unsure why you're getting above error read on. due nature of syntax above problem might occur because: 1. exception *might* occur in wrongly written multi-threaded tests. please refer mockito faq on limitations of concurrency testing. 2. sp...
Read more