r - Add XY points conditionally to raster map generated by levelplot -

-

i have tricky ifelse task here. below code showing data 2 periods(future , current). data has mean, 5th , 95th confidence bounds alongside xy cordinates. want compare mean, 5th , 95th confidence bounds (ci) 2 dfs (future , current).

conditions:

1) if cis future not overlap of current , cis of futrue above current, thenpch=2.

2) if cis future not overlap of current , cis of futrue below current, thenpch=3.

3) if cis of future overlap of current, pch=4

library(raster) library(rastervis)      s <- stack(replicate(2, raster(matrix(runif(100), 3))))     current <- data.frame(coordinates(samplerandom(s, 3, sp=true)),                      c5th=c(17.643981,16.83572,9.979904),                      cmean=c(26.66364,19.74286,15.10000),c95th=c(35.68329,22.64999,20.22010))       future <- data.frame(coordinates(samplerandom(s, 3, sp=true)),                           c5th=c(17.643981,16.83572,9.979904)*2,                           cmean=c(26.66364,19.74286,15.10000)*2,c95th=c(35.68329,22.64999,20.22010)*2)

the result of above 3 conditions added map. (just attempt):

levelplot(s, margin=false, at=seq(0, 1, 0.05)) +    layer(sp.points(xy, pch=ifelse(condition, 2, 3,4), cex=2, col=1), columns=1) +   layer(sp.points(xy, pch=ifelse(condition, 2, 3,4), cex=2, col=1), columns=2)

for example, in figure below, if minimum of nfc (future) lies above maximum of afc (current), condtion 1. if maximum of nfc lies below minimum of afc, condtion 1.the plot shown below satisfies condition 3.

enter image description here

please help. at.

it easier if define complete spatialpointsdataframe object additional categorical variable defined according conditions need.

library(raster) library(rastervis)  s <- stack(replicate(2, raster(matrix(runif(1000), 3)))) ## coordinates cc <- samplerandom(s, 3, sp = true) ## data current <- data.frame(c5th=c(17.643981,16.83572,9.979904),                       cmean=c(26.66364,19.74286,15.10000),                       c95th=c(35.68329,22.64999,20.22010))  future <- data.frame(c5th=c(17.643981,16.83572,9.979904)*2,                      cmean=c(26.66364,19.74286,15.10000)*2,                      c95th=c(35.68329,22.64999,20.22010)*2)  cf <- data.frame(current, future) ## define categorical variable checking conditions need cf$class <- with(cf,                  ifelse(c5th > c95th.1, 'a',                         ifelse(c95th < c5th.1, 'b',                                ifelse(c5th < c95th.1 && c5th > c5th.1, 'c', 'd')                                )                         )                  ) cf$class <- factor(cf$class)  ## build spdf object coordinates , data pp <- spatialpointsdataframe(cc, cf)

this object can displayed spplot. can choose symbols, sizes, etc.

levelplot(s) + spplot(pp["class"],                       pch = 21:23,                       col.regions = 'gray')

Comments

Post a Comment

Popular posts from this blog

java - Could not locate OpenAL library -

-
i writing game using slick2d library, uses lwjgl. lwjgl version 2.9.3. using netbeans 8. i tried adding sounds game, example music track on main menu. problem is, exception when run game, despite game still compiling , running except sound of course: sat mar 21 01:23:57 edt 2015 error:could not locate openal library. org.lwjgl.lwjglexception: not locate openal library. @ org.lwjgl.openal.al.create(al.java:156) @ org.lwjgl.openal.al.create(al.java:156) @ org.lwjgl.openal.al.create(al.java:102) @ org.lwjgl.openal.al.create(al.java:206) @ org.newdawn.slick.openal.soundstore$1.run(soundstore.java:295) @ java.security.accesscontroller.doprivileged(native method) @ org.newdawn.slick.openal.soundstore.init(soundstore.java:292) @ org.newdawn.slick.music.<init>(music.java:156) @ org.newdawn.slick.music.<init>(music.java:75) @ view.mainmenu.init(mainmenu.java:58) i didn't include rest of error states error happens in game code. the...
Read more

c++ - Delete matches in OpenCV (Keypoints and descriptors) -

-
i want check scene image against 2 train images. that, detect features , compute descriptors of both training images. before detecting, computing , matching scene image, delete matches of train1 , train2. because these matches won't facilitate matching of scene image train1 , train2. so, match train1 train2 , vector of matches trainidx , queryidx. how can delete these matches in keypoints-vector , descriptor matrix of train1 , train2? best regards, dwi i have done below: std::vector<cv::keypoint> keypoints[2]; cv::mat descriptor[2]; std::vector< cv::dmatch > matches; /* write code generate keypoints, descriptors , matches here... keypoint[0] -> train image 1 keypoints keypoint[1] -> train image 2 keypoints descriptor[0] -> train image 1 descriptors descriptor[1] -> train image 2 descriptors matches -> matched between train image 1 , 2 */ // logic keep unmatched keypoints , corresponding descriptor...
Read more

sorting - opencl Bitonic sort with 64 bits keys -

-
i used nvidia sdk bitonic sort , works great me. 32 bits (uint) need ulong keys. have typically 2^14 keys @ time , power of 2. searched on not find kernel designed ulong. i tried modified nvidia sdk bitonic sort using ulong keys not work. kernel not crash after de call clenqueuendrangekernel got error : cl_invalid_command_queue can tell me how modify bitonic sort or radixsort, or can sort ulong keys? i running nvidia quadro 4000 opencl 1.1 cuda 6.5.20 full_profile used original nvidia sdk bitonicsort.cl i used"ulong" keys , value input , ouput instead of "uint" thanks helping i made work-around suit need. used 30 bits of key , 18 high bits of value sorting. changed comparators (2 lines): '#define mask 0xffffc000 // suit need from this: if( (*keya > *keyb) == dir ) this: if(((*keya > *keyb)||((*keya == *keyb)&&((*vala & mask)> (*valb & mask)))) == dir ) it works fine. but, still, it's work-a...
Read more