python - ArcPy: Initialize zero valued grid Raster and iteratively increment -


i'm jumping right arcpy out friend needs automate arcgis work. problem @ hand can reduced following pseudocode:

accumulationraster = somezeroedraster each raster pair, r, s, pertaining rain , snow:       accumulationraster += con(s == 0, con(r >= 1, 1, 0), 0)

in other words, iterate through data of rain , snow of each day. perform logical operation on each cell pair, corresponding cell in accumulationraster gets incremented if no snow (s == 0) , rain (r >= 1).

having no luck far various attempts, thought i'd ask here, problem should have easy solution. full working code missing gaps pointed out below. of trivial file path set up, , above pseudocode should suffice:

the missing pieces numbered 1, 2, , 3.

  1. how initializase 0 raster w/proper dimensions.

  2. how make copy of 0 raster. if 1 simple, instead.

  3. increment accumulation raster.

import os import arcpy arcpy import env arcpy.sa import *  in_workspace = "c:/data/projectfolder" os.chdir(in_workspace) arcpy.env.workspace = in_workspace  maaneder = ["01_januar", "02_februar", "03_mars", "04_april",             "05_mai", "06_juni", "07_juli", "08_august",             "09_september", "10_oktober", "11_november", "12_desember"]  ################ 1 ################### baseaccum = #1: ?? (zero raster w/correct type & dimensions) ######################################  outputpath = "raindaysingrowthseason/" year in range(1997,2015):     month in range(1, 13):         ################ 2 ###################         monthaccum = #2: ?? (copy of basesum)         ######################################          day in range(1,32):             # e.g 1997/09_september/             fpath = str(year) + "/" + maaneder[month-1] + "/"             # e.g. 1997_09_08.asc             filesuffix = str(year) + "_" + str(month).zfill(2) + "_" + str(day).zfill(2) + ".asc"             # e.g. snoedybde/1997/09_september/sd_1997_09_08.asc             snoefile = "snoedybde/" + fpath + "sd_" + filesuffix             # e.g. rrl/1997/09_september/rrl_1997_09_08.asc             rrlfile = "rrl/" + fpath + "rrl_" + filesuffix              hassnoe = os.path.isfile(snoefile)             hasrrl = os.path.isfile(rrlfile)             if (not hassnoe or not hasrrl):                 continue              snoeraster = raster(snoefile)             rrlraster = raster(rrlfile)              ############## 3 #####################             monthaccum = #3: ??                          # a'la:                          # += con(snoeraster == 0, con(rrlraster >= 1, 1, 0), 0)             ######################################          # e.g. raindaysingrowthseason/1997         outputdir = outputpath + str(year)         if (not os.path.isdir(outputdir)):             os.mkdir(outputdir)         monthaccum.save(outputdir + "/r_" + str(year) + "_" + str(month).zfill(2))

you can use arcpy.numpyarraytoraster make 0 base raster.

if wanted raster 50 rows , 50 columns 1 meter resolution. can this.

import numpy  myarr = numpy.zeros(50, 50)  myzeroraster = arcpy.numpyarraytoraster(myarr, arcpy.point(355355, 3199277), 1, 1)

you can convert raster numpy array , accumulations using 1 numpy array another.

or if have spatial analyst, use tool create constant raster create 0 based raster , use math -> plus tool add rasters together


Comments

Post a Comment

Popular posts from this blog

javascript - Jquery show_hide, what to add in order to make the page scroll to the bottom of the hidden field once button is clicked -

i have button on website when clicked, reveals map. button can seen when page first loads up, when clicked, section holding map slides open , goes out of screen since page loads @ top (normal). want page scroll down bottom of div select when button pushed, if map out of screen. if page scrolled down, , clicks show map button, don't want page jump bottom of map. i tried figure out myself, no luck. seasoned jquery coders able give me guidance? many thanks! map = mapdiv.gmap3("get"); infobox = new infobox({ pixeloffset: new google.maps.size(-220, -310), closeboxurl: '', enableeventpropagation: true }); mapdiv.delegate('.infobox .close','click',function () { infobox.close(); }); jquery(".slidingdiv").hide(); jquery(".show_hide").show(); jquery('.show_hide').click(function(){ lastcenter=map.getcenter(); jquery(".slidingdiv&quo
Read more

javascript - Highcharts multi-color line -

Image
i looking create highchart multi-colored line. if goes down under value, become red, , when comes on value, become green . here picture example: or maybe highcharts have other charts models differentiate high or low level color? you can achieve coloring using lineargradient coloring: color: { lineargradient: { x1: 0, x2: 0, y1: 0, y2: 1 }, stops: [ [0, 'red'], [0.25, 'yellow'], [0.50, 'green'], [0.75, 'yellow'], [1, 'red'] ] } here 0 in stops relates minimum value of series , 1 maximum value, , in between percentages between values. if find percentages according values prior creating chart use correct stops values make relate specific values in data (at least approximately). see this jsfiddle example .
Read more

javascript - Enter key does not work in search box -

tried search guys sorry. input in search box doesnt work (or search) when hit enter when user clicks search button. doing wrong? here code. thanks. <script type="text/javascript"> $(document).ready(function() { $("#submit").click(function() { var url = "http://asla.org/awardssearch.html"; url += "?s=" + $('#googlecse').val(); window.location = url; }); $('#googlecse').keydown(function(event) { if (event.keycode == '13') { $("#submit").click(); } else { //do nothing :) } }); }); </script> <form class="navbar-form navbar-right" role="search"> <div class="search"> <input id="googlecse" type="text" onblur="if(this.value=='')this.value=this.defaultvalue;" onfocus="if(this.value==this.defaultvalue)t
Read more