I have several reclassified rasters (0 = snow cover, 1 = no snow) and I would like to use the multiply function on all rasters that have the same spatial extent. In ArcGIS 10 this function can be performed on only two rasters. Is it possible to do it on multiple rasters in Python?


Actually, this function can be performed on large numbers of rasters at once even in ArcGIS 10 (without scripting). It is equivalent to finding the local minimum, because the product is 1 if and only if the minimum equals 1 and otherwise is 0. You can also obtain a solution using Equal to Frequency, Greater than Frequency, Less than Frequency, or Popularity by comparing the frequency of 1's to the number of rasters. This latter approach yields more specific information than you call for originally: rather than indicating whether snow was ever present, it can tell you the proportion of rasters in which it is present, cell by cell. Retaining such precision in the data is practically guaranteed to produce better models than the crude binary (yes/no, 1/0, present/absent) models often employed in GIS analysis.

  • Awesome, I like your frequency recomendations, especially if the results are to contain more info. I'm assuming that these functions can be performed on multiple rasters? I'll give them a try. – dchaboya Dec 9 '11 at 15:38

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.