Thursday, 9 June 2016

raster - Finding zones with nodata in ArcGIS 9.3


I have a DEM raster with some small and large areas of NoData. I also have a raster defining zones (the zones consists of a rectangular grid).


I want to create a new raster marking the zones where there exists any nodata values in the DEM raster.


My initial approach was to use the Zonal Statistics tool, using any statistics type (for example Maximum) and uncheck the "Ignore NoData" option. According to what I can read in the documentation, this should give me an output raster with NoData in the entire zones where the input raster contained any NoData pixels. (I could then perform some trivial operations to give me the desired final result).


However, this doesn't work as expected. The output raster from the Zonal Statistics tool have some holes, but there are several places where there are definitely NoData in the input raster and data in the output.


So, are there any other possibilities of achieving what I want? Or have I misunderstood how I should use the Zonal statistics tool?



Answer



Perform a zonal summary of the NoData indicator grid, IsNull([my DEM]). This indicator has 1's where NoData exist and 0's elsewhere. Thus, a sum, maximum, range, standard deviation, or average that exceeds zero will all characterize zones containing at least one NoData value. The sum actually tells you how many NoData cells are in each zone.



No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...