geospatial - R: How to extract values from contiguous raster cells that are not touched by SpatialLines? -


i've been trying extract values single attribute raster (area, in m2) overlaps lines (that is, .shp spatiallines).

the problem that, along these lines, raster goes 1 several contiguous cells in directions. using extract function values cells touched lines extracted. thus, when add extracted values lines significant amount of area (m2) lost due cells not touched line , therefore values not extracted.

i tried work around by:

step 1 - first aggregating raster lower resolution (i.e. increasing fact argument) , step 2 - rasterizing lines using aggregated raster (created in step 1) mold make sure rasterized lines thick enough cover horizontal spread of cells in original resolution raster. step 3 - resample rasterized lines (created in step 2) original resolution started with. step 4 - finally, extracted values resampled rasterized lines (created in step 3).

however, didn't quite work total area (m2) varies according fact="" value use when first aggregating raster (in step 1).

i appreciate if has dealt similar problem , can me out here. here codes i've been running try work:

# input raster file g.025 <- raster("ras.asc") g.1 <- aggregate(g.025, fact=2, fun=sum)  # input spatiallines spline1 <- readogr("/users/xxxxx.shp") spline2 <- readogr("/users/xxxxx.shp") spline3 <- readogr("/users/xxxxx.shp")  # rasterizing using low resolution raster (aggregated) c1 <- rasterize(spline1, g.1, field=spline1$type, fun=sum) c2 <- rasterize(spline2, g.1, field=spline2$type, fun=sum) c3 <- rasterize(spline3, g.1, field=spline3$type, fun=sum)  # resampling higher resolution c1 <- resample(c1, g.025) c2 <- resample(c2, g.025) c3 <- resample(c3, g.025)  # preparing extract area (m2) values raster “g.025” c1tab <- as.data.frame(c1, xy=t) c2tab <- as.data.frame(c2, xy=t) c3tab <- as.data.frame(c3, xy=t) c1tab <- c1tab[which(is.na(c1tab$layer)!=t),] c2tab <- c2tab[which(is.na(c2tab$layer)!=t),] c3tab <- c3tab[which(is.na(c3tab$layer)!=t),]  # extracting area (m2) values raster “g.025” c1tab[,4] <- extract(g.025, c1tab[,1:2]) c2tab[,4] <- extract(g.025, c2tab[,1:2]) c3tab[,4] <- extract(g.025, c3tab[,1:2]) names(c1tab)[4] <- "area_m2" names(c2tab)[4] <- "area_m2" names(c3tab)[4] <- "area_m2"  # sum total area (m2) c1_area <- sum(c1tab$area_m2) c2_area <- sum(c2tab$area_m2) c3_area <- sum(c3tab$area_m2) tot_area <- sum(c1_area, c2_area, c3_area)

thanks!

andre


Comments

Post a Comment

Popular posts from this blog

Is there a better way to structure post methods in Class Based Views -

i'm using class based views in django , never specify url in ajax calls , omit action parameter in form because know these requests go through post method of class based view. notice view - profile page make many requests depending on user updating on page - post method become convoluted , if/else mess example: def post(self, request, *args, **kwargs): # handle user changing available date if request.post.get('availabledate') != none: self.updateavailabledate(request) return jsonresponse({'result':'success'}) if request.post.get('newprojectname') != none: creative_user = creativeuserprofile.objects.get(id = request.user.id) project = project.create_project(creative_user, request.post['newprojectname']) project.save() return jsonresponse({'projectid': project.id }) if request.post.get('projectimage') != none , request.files['file']: proj
Read more

performance - Why is XCHG reg, reg a 3 micro-op instruction on modern Intel architectures? -

i'm doing micro-optimization on performance critical part of code , came across sequence of instructions (in at&t syntax): add %rax, %rbx mov %rdx, %rax mov %rbx, %rdx i thought had use case xchg allow me shave instruction , write: add %rbx, %rax xchg %rax, %rdx however, dimay found agner fog's instruction tables , xchg 3 micro-op instruction 2 cycle latency on sandy bridge, ivy bridge, broadwell, haswell , skylake. 3 whole micro-ops , 2 cycles of latency! 3 micro-ops throws off 4-1-1-1 cadence , 2 cycle latency makes worse original in best case since last 2 instructions in original might execute in parallel. now... cpu might breaking instruction micro-ops equivalent to: mov %rax, %tmp mov %rdx, %rax mov %tmp, %rdx where tmp anonymous internal register , suppose last 2 micro-ops run in parallel latency 2 cycles. given register renaming occurs on these micro-architectures, though, doesn't make sense me done way. why wouldn't register renamer
Read more

jquery - Responsive Navbar with Sub Navbar -

good day, forgive me nav vocabulary not best. if question confusing apologize. i'm working on navigation bar needs centered on page @ 768px , above. items in nav open (or perhaps toggle) subnav lives directly below. these nav should appear centered on page. i have been able work through responsive portion of nav bar, having main items align vertically , sub nav menu of each item roll out below, , gave them indent separation. working well. however, issue not responsive side, behavior @ width greater 768px. happens toggling menu's sub navigation pushes other main items down page below sub nav. here's code: $(document).ready(function () { //stick in fixed 100% height behind navbar don't wrap $('#slide-nav.navbar-inverse').after($('<div class="inverse" id="navbar-height-col"></div>')); $('#slide-nav.navbar-default').after($('<div id="navbar-height-col"></d
Read more