nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

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

GiuseppeGentili

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

output_raster.jpg

rcwh0.zip

cle444_3-1718361493524.png

cle444_2-1718361457240.png

Encontre mais publicações marcadas com

Accepted answers

All comments

DanPatterson

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

GiuseppeGentili

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)

DanPatterson

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

GiuseppeGentili

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.

mca.zip

DanPatterson

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

GiuseppeGentili

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.

SteveLynch

I've logged an issue for the problem.

-Steve

GiuseppeGentili

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

Participe da Conversa

Membros conectados podem postar, seguir atualizações e mais. Novo aqui? Registre uma conta gratuita.

Links Rápidos

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.