gis - Better hillshading for map plots in R -


i'm working on improved hillshading topographical map plots. basic hillshade workflow documented in image() is:

require(raster) alt = getdata('alt', country='che') slope = terrain(alt, opt='slope') aspect = terrain(alt, opt='aspect') hill = hillshade(slope, aspect, 40, 270) plot(hill, col=grey(0:100/100), legend=false, main='switzerland') plot(alt, col=rainbow(25, alpha=0.35), add=true)

this image shows plot(hill..) before plot(alt..) applied:

enter image description here

the method creates solid grey under-layer of hillshades on other data layers (e.g. elevation shading) plotted semi-transparently. problem approach (a) neutral colour flat terrain (rbg (202,202,202), '#cacaca') severely shades whole model, (b) prevents multiple shade layering, such used 'swiss hillshade' approach.

i can imagine workaround converts rasters matrices , applies hillshading numerical multiplier brightness of other layers, doesn't seem elegant (although may wrong). wonder if has ideas or (preferably) experience in area? in advance.

no experience this, why not give gray undermap alpha value depends on slopes? here's try:

# before require(raster) alt = getdata('alt', country='che') slope = terrain(alt, opt='slope') aspect = terrain(alt, opt='aspect') hill = hillshade(slope, aspect, 40, 270) plot(hill, col=grey(0:100/100), legend=false, main='switzerland') plot(alt, col=rainbow(25, alpha=0.35), add=true)

enter image description here

as say, dark.

# after grayalphas <- seq(-1,1,by=.01)^2*100 grayalphas[grayalphas==100] <- 99 plot(hill, col=paste0(grey(0:100/100),sprintf("%.2d",grayalphas)), legend=false, main='switzerland')  plot(alt, col=rainbow(25, alpha=0.35), add=true)

enter image description here

i set gray alphas have parabolic shape, minimum gray value .5 , max of 99 @ gray values of 0 or 1. if choose this, you'll want tinker levels, etc, easy implement. plus you'll want put more effort did alphas, mine strictly numeric , not hex.

[edit] found nifty function adding alphas, addtrans() here in sacha epskamp's answer. keeps parabola, ranges 0 in middle 255 on extremes.

grayalphas <- seq(-1,1,length=101)^2*255 plot(hill, col=addtrans(grey(0:100/100),grayalphas), legend=false, main='switzerland') plot(alt, col=rainbow(25, alpha=0.35), add=true)

enter image description here


-

Comments

Post a Comment

Popular posts from this blog

html - Sizing a high-res image (~8MB) to display entirely in a small div (circular, diameter 100px) -

first of all, thank being patient know simple question more couple of weeks' experience figure out. i've seen folks persevere , helpful on here , hope you're willing me. here fiddle: http://jsfiddle.net/65uzf/ here 1 of places searched answers on here: ( vertically center image thumbnails fix-sized div ) tried solutions on there , didn't work out me. i've searched ~4 or 5 others, , been unsuccessful those, too. what want take large image (what linked background-image #ysypic , #tabpic divs) , display thumbnail version of image in confines of divs have set it. is, using background properties have laid out in css, portion of image displays inside of circular divs. without these properties, display defaults smaller section of image [presumably pi*(50^2) area in top left] terrible resolution. i not care if entire image doesn't display inside of circles (it won't, given difference in shape), , overflow not being displayed fine. want able see figu
Read more

java - IntelliJ - No such instance method -

using intellij's 12 ultimate, i'm running following code in debugger: java import play.api.libs.json.jsvalue; public class foo { ... public jsvalue tojson() { ... } public class fooexample { ... foo foo = new foo(); system.out.println("...); //<-- breakpoint at breakpoint, right-clicked source code , picked "evaluate expression," typing in: foo.tojson() . but following error showed up: no such instance method: play.api.libs.json.jsvalue$class.com.foo.foo.tojson () am doing wrong? foo#tojson calls scala code, if matters. edit had breakpoint after instantiation of foo. downvoted it, deserved it. edit my answer no longer valid after op modified question (i.e. moved breakpoint actual is).. a breakpoint hit before line executed. in case foo has not yet been declared or instantiated (i.e. foo constructor has not yet been called). you'll need put breakpoint @ next line (or step on current line) if want evaluate foo .
Read more

identifier - Is it possible for an html5 document to have two ids? -

i wondering if html5 element can have 2 ids, this: <button id="thing1" id="thing2" name="button1" onclick="change()">click here</button> or this: <button id="thing1, thing2" name="button1" onclick="change()">click here</button> also, if 1 of ids said 1 thing, , other id said thing, 1 used? tl;dr : no. in versions of html, each element can have 1 id. see the spec : the id attribute specifies element's unique identifier (id). [dom] the value must unique amongst ids in element's home subtree , must contain @ least 1 character. the value must not contain space characters. note: there no other restrictions on form id can take; in particular, ids can consist of digits, start digit, start underscore, consist of punctuation , etc. note: element's unique identifier can used variety of purposes, notably way link specific parts of document usin
Read more