gis - Plot specific point data on each layer in stack with rasterVis in R -

-

i have 2 rasters stacked together:

library(rastervis) r1 <- raster(system.file("external/test.grd", package="raster"))  r2 <- r1 / 2 r.stack <- stack(r1, r2)

since highlight areas later each specific layer in stack, create 2 point datasets based on raster values:

pts1 <- rastertopoints(r1, spatial=t) idx <- which(as.data.frame(pts1)[, 1] >= 400) pts1 <- pts1[idx, 1]  pts2 <- rastertopoints(r2, spatial=t) idx <- which(as.data.frame(pts2)[, 1] >= 400) pts2 <- pts2[idx, 1]

now, plot raster stack levelplot rastervis in r. overlay r1 pts1 , r2 pts2.

however, add 1 point dataset, used both layers:

levelplot(r.stack) + layer(sp.points(pts1, pch=20, cex=0.1, col="black"))

how can use specific point datasets specific layers while still using raster stacks?

i avoid making own subplot plotting each layer its' specific point dataset alone followed using print.trellis. tried it, result inferior compared levelplot raster stacks.

any idea on how achieve that?

with panel.number function can subset data according panel located:

pts <- list(pts1, pts2)  levelplot(r.stack) +     layer(sp.points(pts[[panel.number()]],                     pch=20, cex=0.1, col="black"))

Comments

Post a Comment

Popular posts from this blog

yii2 - Yii 2 Running a Cron in the basic template -

-
i have seen instructions on running cron in advanced template, cant figure out how on basic template. have following controller in basic controllers folder. <?php namespace app\controllers; use yii\console\controller; /** * cron controller */ class croncontroller extends controller { public function actionindex() { echo "cron service runnning"; } public function actionmail($to) { echo "sending mail " . $to; } } i have navigated root of application , tried these comands yii cron php yii cron im getting unknown comand "cron" i know old question failed find actual answer this. if @ code yii exec file you'll notice requires /common/config/aliases.php file console/config/main.php. can naturally change these or copy on these advanced template. hope solves problem wanting same thing.
Read more

asp.net - 'System.Web.HttpContext' does not contain a definition for 'GetOwinContext' Mystery -

-
i realize question might seem trivial some, it's these types of things find myself fighting quite bit , want make sense of despite seeming losing battle in .net (for me anyway). so, if following: using system.web; ... applicationuser user = system.web.httpcontext.getowincontext().getusermanager<applicationusermanager>().findbyid(system.web.httpcontext.user.identity.getuserid()); that produces error in title , red getowincontext() , error cannot resolve symbol 'getowincontext()' however, if following (remove system.web in front of httpcontext ), works expected (or @ least no errors): using system.web; ... applicationuser user = httpcontext.getowincontext().getusermanager<applicationusermanager>().findbyid(system.web.httpcontext.user.identity.getuserid()); however, if (same line that's working using system.web commented out): //using system.web; ... applicationuser user = httpcontext.getowincontext().getusermana...
Read more

mercurial graft feature, can it copy? -

-
there new feature in mercurial 3, called 'graft' (graft local). job of moving change-set different branch. is there way "copy" change-set branch? i interested in in moving stuff qa branch production branch, still need code reside in qa. want copy. found , have tried few things, maybe it'll become clear more tinkering , reading... does have better way of 'cherry picking' change-sets move cross branch (or trunk)??? graft copies changesets, doesn't move them. hg graft (emphasis mine): this command uses mercurial's merge logic copy individual changes other branches without merging branches in history graph. known 'backporting' or 'cherry-picking' . default, graft copy user, date, , description source changesets.
Read more