Hfrhyu

Is it tax fraud for an individual to declare non-taxable revenue as taxable income? (US tax laws)

What's the output of a record cartridge playing an out-of-speed record

How can bays and straits be determined in a procedurally generated map?

Why doesn't Newton's third law mean a person bounces back to where they started when they hit the ground?

Can an x86 CPU running in real mode be considered to be basically an 8086 CPU?

Pattern match does not work in bash script

Can a Warlock become Neutral Good?

Why don't electron-positron collisions release infinite energy?

Arthur Somervell: 1000 Exercises - Meaning of this notation

Why does Kotter return in Welcome Back Kotter?

Mage Armor with Defense fighting style (for Adventurers League bladeslinger)

Show that if two triangles built on parallel lines, with equal bases have the same perimeter only if they are congruent.

Why "Having chlorophyll without photosynthesis is actually very dangerous" and "like living with a bomb"?

Is it legal for company to use my work email to pretend I still work there?

How to find program name(s) of an installed package?

Modeling an IP Address

How old can references or sources in a thesis be?

The use of multiple foreign keys on same column in SQL Server

Have astronauts in space suits ever taken selfies? If so, how?

What does it mean to describe someone as a butt steak?

What typically incentivizes a professor to change jobs to a lower ranking university?

Why is 150k or 200k jobs considered good when there's 300k+ births a month?

How did the USSR manage to innovate in an environment characterized by government censorship and high bureaucracy?

Is a conference paper whose proceedings will be published in IEEE Xplore counted as a publication?

Reclassifying raster with 33 bands in ArcMap?

Display a raster stack with several bands in ArcGIS?Bands disappear after Raster Calculation ArcGISCreating Quantile classes when reclassifying raster in QGIS?Reclassifying multiple rasters in ArcGIS Desktop?Why is ArcPy script looping through bands slow?Is it possible to read in specific bands from a multi-band raster with gdal or rasterio?Reclassify Raster in ArcGIS Pro 2 by Multiple BandsRetaining integer values in raster when reclassifying values in QGIS?Identify spectral bands in a multispectral raster datasetGet colortable from three-band RGB raster (brick) with R

.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;

0

1

I have a file geodatabase raster dataset with 33 bands. Each band is comprised of pixels with values 1-31. This represents 33 years of landuse cover data (with 31 different landuse classes). I'm trying to move from this format to a single TIFF file that has 1 band comprised of pixels with the value 1-33, corresponding to the year that the pixel switched from forest to non-forest. Its a huge file (100GB), so I need to reduce the file size before I can do any analysis.

My plan is to reclassify each band/year into 0 (forest) and 1 (non-forest), then sum all the bands (using the raster calculator, I think?). However, the reclassify tool is only reclassifying a single band.

How do I do this for all the bands.

arcmap reclassify multi-band

share|improve this question

edited Apr 2 at 21:37

PolyGeo♦

53.9k1781245

asked Apr 2 at 21:10

Jacy HydeJacy Hyde

183

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  If result of sum will say 2, what that means for you?
  
  – FelixIP
  Apr 3 at 1:41
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  @FelixIP I could use it to calculate the year of deforestation. The dataset goes through 2017, so for each pixel, it would be 2017- (final raster value) = year deforested
  
  – Jacy Hyde
  Apr 3 at 13:56
  

  
  
  
  

add a comment | 

0

1

I have a file geodatabase raster dataset with 33 bands. Each band is comprised of pixels with values 1-31. This represents 33 years of landuse cover data (with 31 different landuse classes). I'm trying to move from this format to a single TIFF file that has 1 band comprised of pixels with the value 1-33, corresponding to the year that the pixel switched from forest to non-forest. Its a huge file (100GB), so I need to reduce the file size before I can do any analysis.

My plan is to reclassify each band/year into 0 (forest) and 1 (non-forest), then sum all the bands (using the raster calculator, I think?). However, the reclassify tool is only reclassifying a single band.

How do I do this for all the bands.

arcmap reclassify multi-band

share|improve this question

edited Apr 2 at 21:37

PolyGeo♦

53.9k1781245

asked Apr 2 at 21:10

Jacy HydeJacy Hyde

183

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  If result of sum will say 2, what that means for you?
  
  – FelixIP
  Apr 3 at 1:41
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  @FelixIP I could use it to calculate the year of deforestation. The dataset goes through 2017, so for each pixel, it would be 2017- (final raster value) = year deforested
  
  – Jacy Hyde
  Apr 3 at 13:56
  

  
  
  
  

add a comment | 

I have a file geodatabase raster dataset with 33 bands. Each band is comprised of pixels with values 1-31. This represents 33 years of landuse cover data (with 31 different landuse classes). I'm trying to move from this format to a single TIFF file that has 1 band comprised of pixels with the value 1-33, corresponding to the year that the pixel switched from forest to non-forest. Its a huge file (100GB), so I need to reduce the file size before I can do any analysis.

My plan is to reclassify each band/year into 0 (forest) and 1 (non-forest), then sum all the bands (using the raster calculator, I think?). However, the reclassify tool is only reclassifying a single band.

How do I do this for all the bands.

arcmap reclassify multi-band

share|improve this question

edited Apr 2 at 21:37

PolyGeo♦

53.9k1781245

asked Apr 2 at 21:10

Jacy HydeJacy Hyde

183

share|improve this question

edited Apr 2 at 21:37

PolyGeo♦

53.9k1781245

asked Apr 2 at 21:10

Jacy HydeJacy Hyde

183

share|improve this question

edited Apr 2 at 21:37

PolyGeo♦

53.9k1781245

asked Apr 2 at 21:10

Jacy HydeJacy Hyde

183

edited Apr 2 at 21:37

PolyGeo♦

53.9k1781245

edited Apr 2 at 21:37

PolyGeo♦

53.9k1781245

asked Apr 2 at 21:10

Jacy HydeJacy Hyde

183

asked Apr 2 at 21:10

Jacy HydeJacy Hyde

183

1 Answer
1

active

oldest

votes

0

Big task it is. So set environment settings - cell size, snap raster and extent, to ones of that raster.
Run script below from ArcMap or ArcCatalog:

import arcpy
from arcpy.sa import *
## change next 2 lines
ForestCode=184
multiBand="F:/AERIALS/Rectified/image0051.tif"
## output raster name
outRaster = arcpy.GetParameterAsText(0)

nBands=int(arcpy.GetRasterProperties_management(multiBand, "BANDCOUNT")[0])
outList=[]
for i in range(nBands):
 arcpy.MakeRasterLayer_management(multiBand,"AFTER","","",i+1)
 after=arcpy.Raster("AFTER")
 if i>0:
 deforest=Con((before==ForestCode) & (after != ForestCode),i)
 outList.append(deforest)
 before=after*1
 arcpy.AddMessage("Band %i processed" %i) 
arcpy.AddMessage("Computing last year of deforestation ...")
deforest=CellStatistics(outList, "MAXIMUM", "DATA")
deforest.save(outRaster)

I strongly suggest to test it on a tiny subset first, i.e. clip bit of that monster in the area where deforestation certainly happened.

Update:

Replace 184 by value that represents forest in your raster.
Script compares consecutive pair of bands - before and after - to find when last deforestation happened.

If you think that 'before' is not defined during first reference, have a better look at the code.

share|improve this answer

edited yesterday

answered Apr 3 at 19:43

FelixIPFelixIP

16.6k11642

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  thanks for answering. I'm going to be honest, I have absolutely no idea what this code is doing. I have a few specific questions, but would you mind explaining what's happening here in general? I also specifically don't understand the line ForestCode= 184. What is this line doing/ what is 184? also, can you explain what's happening with before==ForestCode?-- "before" isn't defined prior to this line.
  
  – Jacy Hyde
  yesterday
  

  
  
  
  

add a comment | 

Your Answer

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "79"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);



);

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);

Sign up using Google

Sign up using Facebook

Sign up using Email and Password

Post as a guest

Name

Email

Required, but never shown

StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f317558%2freclassifying-raster-with-33-bands-in-arcmap%23new-answer', 'question_page');

);

Post as a guest

Name

Email

Required, but never shown

0

Big task it is. So set environment settings - cell size, snap raster and extent, to ones of that raster.
Run script below from ArcMap or ArcCatalog:

import arcpy
from arcpy.sa import *
## change next 2 lines
ForestCode=184
multiBand="F:/AERIALS/Rectified/image0051.tif"
## output raster name
outRaster = arcpy.GetParameterAsText(0)

nBands=int(arcpy.GetRasterProperties_management(multiBand, "BANDCOUNT")[0])
outList=[]
for i in range(nBands):
 arcpy.MakeRasterLayer_management(multiBand,"AFTER","","",i+1)
 after=arcpy.Raster("AFTER")
 if i>0:
 deforest=Con((before==ForestCode) & (after != ForestCode),i)
 outList.append(deforest)
 before=after*1
 arcpy.AddMessage("Band %i processed" %i) 
arcpy.AddMessage("Computing last year of deforestation ...")
deforest=CellStatistics(outList, "MAXIMUM", "DATA")
deforest.save(outRaster)

I strongly suggest to test it on a tiny subset first, i.e. clip bit of that monster in the area where deforestation certainly happened.

Update:

Replace 184 by value that represents forest in your raster.
Script compares consecutive pair of bands - before and after - to find when last deforestation happened.

If you think that 'before' is not defined during first reference, have a better look at the code.

share|improve this answer

edited yesterday

answered Apr 3 at 19:43

FelixIPFelixIP

16.6k11642

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  thanks for answering. I'm going to be honest, I have absolutely no idea what this code is doing. I have a few specific questions, but would you mind explaining what's happening here in general? I also specifically don't understand the line ForestCode= 184. What is this line doing/ what is 184? also, can you explain what's happening with before==ForestCode?-- "before" isn't defined prior to this line.
  
  – Jacy Hyde
  yesterday
  

  
  
  
  

add a comment | 

0

Big task it is. So set environment settings - cell size, snap raster and extent, to ones of that raster.
Run script below from ArcMap or ArcCatalog:

import arcpy
from arcpy.sa import *
## change next 2 lines
ForestCode=184
multiBand="F:/AERIALS/Rectified/image0051.tif"
## output raster name
outRaster = arcpy.GetParameterAsText(0)

nBands=int(arcpy.GetRasterProperties_management(multiBand, "BANDCOUNT")[0])
outList=[]
for i in range(nBands):
 arcpy.MakeRasterLayer_management(multiBand,"AFTER","","",i+1)
 after=arcpy.Raster("AFTER")
 if i>0:
 deforest=Con((before==ForestCode) & (after != ForestCode),i)
 outList.append(deforest)
 before=after*1
 arcpy.AddMessage("Band %i processed" %i) 
arcpy.AddMessage("Computing last year of deforestation ...")
deforest=CellStatistics(outList, "MAXIMUM", "DATA")
deforest.save(outRaster)

I strongly suggest to test it on a tiny subset first, i.e. clip bit of that monster in the area where deforestation certainly happened.

Update:

Replace 184 by value that represents forest in your raster.
Script compares consecutive pair of bands - before and after - to find when last deforestation happened.

If you think that 'before' is not defined during first reference, have a better look at the code.

share|improve this answer

edited yesterday

answered Apr 3 at 19:43

FelixIPFelixIP

16.6k11642

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  thanks for answering. I'm going to be honest, I have absolutely no idea what this code is doing. I have a few specific questions, but would you mind explaining what's happening here in general? I also specifically don't understand the line ForestCode= 184. What is this line doing/ what is 184? also, can you explain what's happening with before==ForestCode?-- "before" isn't defined prior to this line.
  
  – Jacy Hyde
  yesterday
  

  
  
  
  

add a comment | 

Big task it is. So set environment settings - cell size, snap raster and extent, to ones of that raster.
Run script below from ArcMap or ArcCatalog:

import arcpy
from arcpy.sa import *
## change next 2 lines
ForestCode=184
multiBand="F:/AERIALS/Rectified/image0051.tif"
## output raster name
outRaster = arcpy.GetParameterAsText(0)

nBands=int(arcpy.GetRasterProperties_management(multiBand, "BANDCOUNT")[0])
outList=[]
for i in range(nBands):
 arcpy.MakeRasterLayer_management(multiBand,"AFTER","","",i+1)
 after=arcpy.Raster("AFTER")
 if i>0:
 deforest=Con((before==ForestCode) & (after != ForestCode),i)
 outList.append(deforest)
 before=after*1
 arcpy.AddMessage("Band %i processed" %i) 
arcpy.AddMessage("Computing last year of deforestation ...")
deforest=CellStatistics(outList, "MAXIMUM", "DATA")
deforest.save(outRaster)

I strongly suggest to test it on a tiny subset first, i.e. clip bit of that monster in the area where deforestation certainly happened.

Update:

Replace 184 by value that represents forest in your raster.
Script compares consecutive pair of bands - before and after - to find when last deforestation happened.

If you think that 'before' is not defined during first reference, have a better look at the code.

share|improve this answer

edited yesterday

answered Apr 3 at 19:43

FelixIPFelixIP

16.6k11642

import arcpy
from arcpy.sa import *
## change next 2 lines
ForestCode=184
multiBand="F:/AERIALS/Rectified/image0051.tif"
## output raster name
outRaster = arcpy.GetParameterAsText(0)

nBands=int(arcpy.GetRasterProperties_management(multiBand, "BANDCOUNT")[0])
outList=[]
for i in range(nBands):
 arcpy.MakeRasterLayer_management(multiBand,"AFTER","","",i+1)
 after=arcpy.Raster("AFTER")
 if i>0:
 deforest=Con((before==ForestCode) & (after != ForestCode),i)
 outList.append(deforest)
 before=after*1
 arcpy.AddMessage("Band %i processed" %i) 
arcpy.AddMessage("Computing last year of deforestation ...")
deforest=CellStatistics(outList, "MAXIMUM", "DATA")
deforest.save(outRaster)

share|improve this answer

edited yesterday

answered Apr 3 at 19:43

FelixIPFelixIP

16.6k11642

answered Apr 3 at 19:43

FelixIPFelixIP

16.6k11642

answered Apr 3 at 19:43

FelixIPFelixIP

16.6k11642