When I run this example, I get 4 output bands which corresponds to each of the input bands:
http://www.harrisgeospati...sholdRasterTask.html

If you replace your file with the sample dataset, do you get an error with this code? Is there additional code that is causing this message?

Sorry, I did not try the code on the sample dataset, because I already overwrote the program file with new code. To solve the problem, I took a single band and used the task, BinaryGTThresholdRaster. However, the new code is not portable, and will create problems when I change data sets.

Now I am trying to figure out how to use the pixel state for band statistics. In the help file:

> maskedData = maskedRaster.GetData(BAND=0, PIXEL_STATE=pixelState)
>pos = WHERE(pixelState EQ 0, count)
>IF (count NE N_Elements(maskedData)) THEN $
> maskedData = maskedData[pos]

I get that I only want pixels with the state = 0, but the example makes no sense to me. I assume the 2nd statement counts the number of pixels with state = 0, but how would the count not being equal to N_elements help to identify which pixels to exclude?

Thanks

Please ignore post on 20 March 2019

hi

Infrastructure Engineer

https://www.fieldengineer.com/skills/what-is-an-infrastructure-engineer