Using ArcGIS 10.1 with Spatial and 3D Analysts, and Python 2.7
I have a series of reclassified rasters in which every cell is either one or zero. I need to programmatically create a single raster which shows when each cell first had a value of one. The rasters contain the date in the filename, as well as in a date field in the raster. (More detailed description of the data below)
I am completely stymied. I've tried using the Plus function in Spatial Analyst, but I lose the date. I have very little experience working with rasters, so I'm hoping that the problem is my unfamiliarity with the terminology, and that someone can point me in the right direction.
If I can figure out how to do it ‘by hand’ in ArcMap, I’m confident I can create the scripts to automate the process. Of course, a ready to go script would be very nice :^)
The data is 35 years of sea-ice concentration shapefiles from the Canadian Ice Service. Clipped tifs were generated and reclassified so that each cell in the raster has a value of either one, if the ice concentration is below 50% (break-up), or zero if it is not (i.e. concentration greater than 50% or land). I then added the date, parsed from the filename, using Create and CalculateField_management. The objective is to create annual maps showing the progression of the ice break-up to use in biological analyses.