maps - Clipping rasters in R -


i building map northeastern u.s.a. map background needs either altitude map or mean annual temperature map. have 2 rasters worldclim.org give me these variables need clip them extent of states interested in. suggestions on how this. have far:

   #load libraries library (sp) library (rgdal) library (raster) library (maps) library (mapproj)   #load data state<- data (statemapenv) elevation<-raster("alt.bil") meantemp<-raster ("bio_1.asc")  #build raw map nestates<- c("maine", "vermont", "massachusetts", "new hampshire" ,"connecticut",   "rhode island","new york","pennsylvania", "new jersey",   "maryland", "delaware", "virginia", "west virginia")  map(database="state", regions = nestates, interior=t,  lwd=2) map.axes()  #add site localities sites<-read.csv("sites.csv", header=t) lat<-sites$latitude lon<-sites$longitude  map(database="state", regions = nestates, interior=t, lwd=2) points (x=lon, y=lat, pch=17, cex=1.5, col="black") map.axes() library(maps)                                                                  #add axes  main map map.scale(x=-73,y=38, relwidth=0.15, metric=t,  ratio=f)  #create inset map   # next, create new graphics space in lower-right hand corner.  numbers proportional distances within graphics window (xmin,xmax,ymin,ymax) on scale of 0 1.   # "plt" key parameter adjust     par(plt = c(0.1, 0.53, 0.57, 0.90), new = true)    # think key command http://www.stat.auckland.ac.nz/~paul/rgraphics/examples-map.r     plot.window(xlim=c(-127, -66),ylim=c(23,53))    # fill box white     polygon(c(0,360,360,0),c(0,0,90,90),col="white")    # draw map     map(database="state", interior=t, add=true, fill=false)     map(database="state", regions=nestates, interior=true, add=true, fill=true, col="grey")

the elevation , meantemp objects ones need clipped area extent of nestates object. input

the function crop in raster package allows use extent object or object extent can calculated cut(subset) object. package example:

r <- raster(nrow=45, ncol=90) r[] <- 1:ncell(r) e <- extent(-160, 10, 30, 60) rc <- crop(r, e)

if wanted cut in more detailed manner maybe use shp of states , on function in sp package.


Comments

Post a Comment

Popular posts from this blog

.htaccess - First slash is removed after domain when entering a webpage in the browser -

sometimes when hit address of webpage in browser, removes first slash after domain. examples: www.mywebsite.com/scripts/script_6.php -> trying load www.mywebsite.comscripts/script_6.php www.mywebsite.com/test.php -> trying load www.mywebsite.comtest.php after entering again, loads page. a .htaccess settings must wrong: rewritebase / errordocument 404 http://www.mywebsite.com/oops.php errordocument 401 /notallowed.php rewriteengine on # turn on rewriting engine #gzip <ifmodule mod_deflate.c> addoutputfilterbytype deflate text/text text/html text/plain text/xml text/css application/x- javascript application/javascript </ifmodule> #end gzip rewritecond %{http_host} ^mywebsite\.com rewriterule ^(.*)$ http://www.mywebsite.com$1 [r=permanent,l] directoryindex index6.php rewriterule ^termsofuse/?$ termsofuse.php [nc,l] # handle requests "adpolicy" rewriterule ^privacy-policy/?$ privacy-policy.php [nc,l] # handle requ...
Read more

Automatically create pages in phpfox -

i have array of information of places , wanna use them create multiple pages in phpfox site. first tried manually insert data in phpfox database (in tables of phpfox_pages , phpfox_pages_text ). the page shown in site when opening it's link, site says: the page looking cannot found. do know other tables should manipulate make work? or know plugin phpfox same? thanks, solved ( @purefan guide): "pages" has user if take @ phpfox_user table, you'll found out how fetch new records on that. there 2 ambiguous fields in table ( password & password_salt ). can use script create values fields: <?php function generaterandomstring() { $length = 164; $characters = '-_0123456789abcdefghijklmnopqrstuvwxyzabcdefghijklmnopqrstuvwxyz'; $randomstring = ''; ($i = 0; $i < $length; $i++) { $randomstring .= $characters[rand(0, strlen($characters) - 1)]; } return $randomstring; } $ssalt = '';...
Read more

c# - Farseer ContactListener is not working -

i'm using farseer in xna project, have trouble contactlistener. created class contactlistener these 2 error messages , don't know how fix problems. the type or namespace name 'contactlistener' not found (are missing using directive or assembly reference?) the type or namespace name 'contactimpulse' not found (are missing using directive or assembly reference?) what wrong contactlistener class? class mycontactlistener: contactlistener { void begincontact(contact contact) { /* handle begin event */ } void endcontact(contact contact) { /* handle end event */ } void presolve(contact contact, ref manifold oldmanifold) { fixture fixturea = contact.fixturea; fixture fixtureb = contact.fixtureb; if (fixtureb.collisioncategories == category.cat10) { contact.enabled = false; } } void postsolve(contact contact, ref contactimpulse impulse) { /* handle post-solve event */ } }...
Read more