How to divide polygon in R? -


i have shapefile of amazonas state (brazil) , other 6 biggest rivers in state (negro, solimões, amazonas, madeira, purus , juruá). want divide state using rivers, interfluvial areas (madeira-purus, purus-juruá, etc). want final 6 regions delimited rivers, each region different polygon. how do that?

enter image description here

i'm finding "clipping" algorithms, , give me area of rivers inside state, , that's not want.

following @jbaums comment, used gdifference, package rgeos:

intflv <- gdifference(state,rivers)

but since "state" has 1 polygon, intflv became spatialpolygons object 1 polygon, made of thousands of sub-polygons. work better in gis software, wanted object thousands of polygons, each single sub-polygon. did:

l <- list() total <- length(intflv@polygons[[1]]@polygons) (i in 1:total) {     print(paste(i,total,sep="/"))     flush.console()     p <- intflv@polygons[[1]]@polygons[[i]]     l[[length(l)+1]] <- polygons(list(p),i) } sp <- spatialpolygons(l) d <- data.frame(ids=1:total) intflv1 <- spatialpolygonsdataframe(sp,data=d) writeogr(intflv1,"shp","intflv","esri shapefile")

the result (after keeping biggest areas):

enter image description here


Comments

Post a Comment

Popular posts from this blog

android - Get AccessToken using signpost OAuth without opening a browser (Two legged Oauth) -

i using signpost oauth access data magento server. have read various tutorials on same , reach point open browser user can enter credentials. however, per requirement have automate part. hence, user should not browser page. have set-up done on server side ( magento ), hit url , returned call page. same through program in android. i have tried below, consumer = new commonshttpoauthconsumer(key, secret); provider = new commonshttpoauthprovider(oauth_init_url,access_token_url, authorize_url); try { provider.retrieverequesttoken(consumer, oauth.out_of_band); log.d("tokens" , consumer.gettoken() + " -- " + consumer.gettokensecret()); } and request tokens. dont know how bypass next step. tried directly accessing accesstoken (stupid of me) provider.retrieveaccesstoken(consumer, "mycallback://callback"); no luck, ends in oauth.signpost.exception.oauthnotauthorizedexception: authorization failed (server replied 401). ...
Read more

org.mockito.exceptions.misusing.InvalidUseOfMatchersException: mockito -

getting error org.mockito.exceptions.misusing.invaliduseofmatchersexception: invalid use of argument matchers! 1 matchers expected, 2 recorded. exception may occur if matchers combined raw values: //incorrect: somemethod(anyobject(), "raw string"); when using matchers, arguments have provided matchers. example: //correct: somemethod(anyobject(), eq("string matcher")); at when( getprogramservice .callservice(any(getprogramtype.class))) .thenreturn(jaxbresponse); please explain error , possible resolution the code posted looks fine; don't see problem it. code before or above causing exception. key here: invalid use of argument matchers! 1 matchers expected, 2 recorded. any , used above, doesn't return "special kind of null" or "special kind of getprogramtype"; don't exist. instead, returns null , side effect puts matcher on stack. whe...
Read more

google shop client API returns 400 bad request error while adding an item -

i'm working on solution add items in google shopping market java application. following tutorial real mess because of missing classes of laste release, after setting 1.5-beta maven managed make base example typecheck correctly classes. still can make simple task of adding 1 item. instead of posting whole code i'm posting part creates product feed. private product createproduct(string id) { product product = new product(); product.title = "red wool sweater"; product.content = new content( "text", "comfortable , soft, sweater keep warm on " + "cold winter nights. red , blue stripes."); product.appcontrol = new appcontrol(); product.appcontrol.addrequireddestination("productads"); product.externalid = id; product.lang = "it"; product.country = "it"; product.condition = "nuovo"; product.price = new price("eu...
Read more