topic Re: Nodata areas in the raster created by conversion tool "Polygon to Raster" in Geoprocessing Questions

Nodata areas in the raster created by conversion tool "Polygon to Raster"

GiuseppeGentili — Tue, 14 Dec 2021 10:28:17 GMT

Hoping to have post this message in the right section, I reports a trouble experienced in the creation of raster from polygon feature class using the conversion tool "PolygonToRaster".
If I I set the environment extent the raster that has been created by the tool reports a "blank" area whose pixels are "nodata". If no extent is set, the raster is created with the correct representation of input polygon coverage (any "nodata" pixels ...)
The extent of the feature class is little bit grater that the output extent I have set.
I have also performed a "Check geometry" before to execute the tool.

Here is the python code that reproduces the wrong result:

import arcpy
import os

arcpy.env.overwriteOutput = True

GDB_locale = r'C:\TEMP\rcwh0\scratch.gdb'
clu_fc = r'C:\TEMP\rcwh0\scratch.gdb\rcwh0'
mca_ras = os.path.join(GDB_locale, "mca_ras")
out_cellsize = 2

desc = arcpy.Describe(clu_fc)

xmin = desc.extent.XMin
xmax = desc.extent.XMax
ymin = desc.extent.YMin
ymax = desc.extent.YMax

print "Extent of %s \nxmin: %s \nxmax: %s \nymin: %s \nymax: %s" % (clu_fc, xmin, xmax, ymin, ymax)

arcpy.env.extent = "359900 5059900 370100 5070100"

xmin = arcpy.env.extent.XMin
xmax = arcpy.env.extent.XMax
ymin = arcpy.env.extent.YMin
ymax = arcpy.env.extent.YMax

print "Env extent \nxmin: %s \nxmax: %s \nymin: %s \nymax: %s" % (xmin, xmax, ymin, ymax)

arcpy.PolygonToRaster_conversion(clu_fc, "CODE_H", mca_ras, "MAXIMUM_COMBINED_AREA", "NONE", out_cellsize)

I get the same wrong raster regardless the cell assignement parameter value set in the tool.

I have attached the file GDB in which there is the polygon feature class (rcwh0) used for the conversion, and the screen shot of the wrong raster created (in red color the nodata)

This trouble was experienced in ArcMap 10.8.1

Re: Nodata areas in the raster created by conversion tool "Polygon to Raster"

DanPatterson — Tue, 14 Dec 2021 11:00:55 GMT

Try the code variant using a "with" statement

Using environment settings in Python—ArcGIS Pro | Documentation

and use Code formatting ... the Community Version - Esri Community

so there are line numbers in your code people can reference

Re: Nodata areas in the raster created by conversion tool "Polygon to Raster"

GiuseppeGentili — Tue, 14 Dec 2021 11:28:10 GMT

Unfortunately the "with" variant does not fix the trouble.
I have got the same bad raster.

Here is the code with the variant:

with arcpy.EnvManager(cellSize=2, extent="359900 5059900 370100 5070100"): xmin = arcpy.env.extent.XMin xmax = arcpy.env.extent.XMax ymin = arcpy.env.extent.YMin ymax = arcpy.env.extent.YMax print "Env extent \nxmin: %s \nxmax: %s \nymin: %s \nymax: %s" % (xmin, xmax, ymin, ymax) arcpy.PolygonToRaster_conversion(clu_fc, "CODE_H", mca_ras, "MAXIMUM_COMBINED_AREA", "NONE", out_cellsize)

Re: Nodata areas in the raster created by conversion tool "Polygon to Raster"

DanPatterson — Tue, 14 Dec 2021 11:36:28 GMT

out_cellsize is defined elsewhere? or use

out_cellsize=cellSize from your "with" statement

What happens when you change the Cell assignment type to the default

and try a *.tif in a folder as output

Re: Nodata areas in the raster created by conversion tool "Polygon to Raster"

GiuseppeGentili — Tue, 14 Dec 2021 12:07:59 GMT

Yes, out_cellsize is set before the "with" statement.
Unfortunately, also the tif file shows the same trouble. I have attached it as zip file.

Re: Nodata areas in the raster created by conversion tool "Polygon to Raster"

DanPatterson — Tue, 14 Dec 2021 15:12:34 GMT

It keeps failing when done manually (aka, no code)

You have thousands of tiny polygons less than a fraction of a m^2, I cleaned those out and it still fails.

I would work with the polygon file to find out what the problem is. Do it in ArcToolbox, I got a warning about the spatial index, try recreating the spatial indices after removing polygons less than your desired cell size

Re: Nodata areas in the raster created by conversion tool "Polygon to Raster"

GiuseppeGentili — Tue, 14 Dec 2021 15:44:50 GMT

I have added a spatial index to the polygon feature class, having left the setting grid sizes to 0,0,0
Anyway, the converted raster again shows the nodata areas.

Without removing any tiny polygons, I have reduced the extent by the value of one cell size.
The conversion worked well !!

There aren't any tiny polygons along the "nodata" area, so I think that the tiny polygons are not responsible of this trouble. I have successfully converted other polygon coverages having thousands of tiny polygons as well.

Re: Nodata areas in the raster created by conversion tool "Polygon to Raster"

SteveLynch — Tue, 14 Dec 2021 18:10:36 GMT

I've logged an issue for the problem.

-Steve

Re: Nodata areas in the raster created by conversion tool "Polygon to Raster"

GiuseppeGentili — Wed, 15 Dec 2021 14:26:48 GMT

Many thanks for your support. I'll look forward to your good news.