clean isolated pixel

Question

Hi,

I'm looking for a way to remove isolated pixels in order to speed up vectorization. In other words, I would like to set the value of its 8 neighbours to any pixel having a different value. I tried Majority but I don't get what I want. This will be coded aftr with Arcobjects. I would be grateful to any idea.
Thanks,

Olivier

OlivierOlivier · Answer

Mark, Curtis, thanks for your kind help,  I'll try your ideas, hope it'll work.

Olivier

curtvprice · Answer

Olivier,I'm looking for a way to remove isolated pixels

I'd try Boundary Clean or Thin (Generalization tools) before you get too fancy.

If you decide you truly need to zap single pixel zones, you can find them the Region Group tool and zap them using Nibble, a tool specifically designed to fill in "holes." 
from arcpy.sa import *
# set all single-cell zones to NoData
catmask = SetNull(Lookup(RegionGroup(catraster),"COUNT") == 1, catraster)
# fill them in with the Nibble tool
cleanras = Nibble(catraster, catmask)

Hope this helps...

MarkBoucher · Answer

Cell statistics compares the value of the cell it is processing with other layers values in the same location. By using and shift function, you can place the surrounding 8 cells "under" the center cell and do some comparisons that you couldn't do otherwise. For example, if you wanted to know if the grid (pixel) were isolated and surrounded by cells that all had the same value, you could use con() to test if the values of the surrounding cells all equaled their average and if the cell had a different value than the average.

I'm sure there are other SA functions that could do this or something like it more eloquently. If not, this might be a way to attack the problems.

Cheers!