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; } } }

SnapRaster not followed in Raster to Polygon

gyetman

I'm using python to convert polygons to a series of rasters and I need the results to exactly match the extents of a set of existing rasters. I tried doing this with SnapRaster, but some of the outputs have an extra row and column. Is this a bug, or have I missed something?

tile <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>image_tile_ws<SPAN class="punctuation token">,</SPAN>tile_name <SPAN class="operator token">+</SPAN> <SPAN class="string token">'.tif'</SPAN><SPAN class="punctuation token">)</SPAN> env<SPAN class="punctuation token">.</SPAN>snapRaster <SPAN class="operator token">=</SPAN> tile env<SPAN class="punctuation token">.</SPAN>extent <SPAN class="operator token">=</SPAN> Raster<SPAN class="punctuation token">(</SPAN>tile<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>extent env<SPAN class="punctuation token">.</SPAN>cellSize <SPAN class="operator token">=</SPAN> tile arcpy<SPAN class="punctuation token">.</SPAN>PolygonToRaster_conversion<SPAN class="punctuation token">(</SPAN>level0<SPAN class="punctuation token">,</SPAN>iso_field<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'l0_'</SPAN> <SPAN class="operator token">+</SPAN> tile_name<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'MAXIMUM_AREA'</SPAN><SPAN class="punctuation token">,</SPAN>weight_field<SPAN class="punctuation token">,</SPAN>env<SPAN class="punctuation token">.</SPAN>cellSize<SPAN class="punctuation token">)</SPAN>‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN></SPAN>

Vind meer berichten getagd met

Accepted answers

DanPatterson_Retired

possible... Snap Raster (Environment setting)—Help | ArcGIS for Desktop because

Tools that honor the Snap Raster environment will adjust the extent of output rasters so that they match the cell alignment of the specified snap raster.

All comments

gyetman

tile = os.path.join(image_tile_ws,tile_name + '.tif')

env.snapRaster = tile

env.extent = Raster(tile).extent

env.cellSize = tile

print('Creating raster...')

arcpy.PolygonToRaster_conversion(level0,iso_field,'l0_' + tile_name, 'MAXIMUM_AREA',weight_field,env.cellSize)

Second go at adding code, the formatting stripped all of the line breaks out. I also should mention that I'm running this with ArcGIS Desktop 10.3.1 on Windows Server r2. I've replicated the problem on a Windows 8 VM as well.

DanPatterson_Retired

possible... Snap Raster (Environment setting)—Help | ArcGIS for Desktop because

Tools that honor the Snap Raster environment will adjust the extent of output rasters so that they match the cell alignment of the specified snap raster.

gyetman

I saw that, but according to the "How the Snap Raster environment setting works" page:

"The lower left corner of the extent is snapped to a cell corner of the snap raster and the upper right corner is adjusted using the output cell size. As a result, when the output cell size is the same as the snap raster cell size, the cells in the output raster are aligned with the cells of the snap raster."

Which is why I set the cell size based on the raster that I'm trying to match.

How the Snap Raster environment setting works—Help | ArcGIS for Desktop

DanPatterson_Retired

but the extent can be adjusted in order for this to occur, hence the possibility of gaining another row/column in order for this to happen

Doe mee aan het gesprek

Aangemelde leden kunnen berichten plaatsen, updates volgen en meer. Nieuw hier? Registreer een gratis account.

Snelle Links

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