topic Re: How to clip map sheets to a Map Gird in Data Management Questions

How to clip map sheets to a Map Gird

JosephPlessis — Fri, 03 Feb 2017 19:17:00 GMT

I have multiple geo referenced scanned maps from the 1950's, Like 60 of them.

I need to clip some white space off the edges so I can apply a fishnet.

So the main question is I have a 50k nts grid where the maps sit perfectly within each grid tile, When i try to clip the map to the grid it does nothing because it is reading the outside border of the nts grid. how would i get the clip to read each individual grid tile for a clip?

thanks

Re: How to clip map sheets to a Map Gird

DarrenWiens2 — Fri, 03 Feb 2017 19:19:31 GMT

Have you tried select (the grid polygon), then clip?

Re: How to clip map sheets to a Map Gird

JosephPlessis — Fri, 03 Feb 2017 19:36:17 GMT

I did selct a grid portion and created a layer form selcted data then proceeded to do the clip.

Your way sounds easier.

Is there a way to script this for I have a lot of maps to do.

Re: How to clip map sheets to a Map Gird

DarrenWiens2 — Fri, 03 Feb 2017 21:24:37 GMT

Yes, you could certainly do this with a Python script, looping through rasters, calculating raster centroid from extent, selecting the grid cell covering the centroid, doing the clip.

Re: How to clip map sheets to a Map Gird

JosephPlessis — Mon, 06 Feb 2017 00:06:00 GMT

is there any good resources to script this

Re: How to clip map sheets to a Map Gird

DarrenWiens2 — Mon, 06 Feb 2017 01:48:18 GMT

Start at the following link, get started on your script, and then post it back here when you need more help.

What is ArcPy?—ArcPy Get Started | ArcGIS Desktop

Basically, figure out the geoprocessing tools you need, find and use the Python syntax on each tool's help page, and then put it in a for loop to repeat. But, first you need to figure out some Python.

Re: How to clip map sheets to a Map Gird

NeilAyres — Mon, 06 Feb 2017 10:55:20 GMT

But if these are scanned map sheets, why not load them into a mosaic, then use the map polygon grid as the footprints.

This will auto magically leave behind all the stuff around the edge.

Re: How to clip map sheets to a Map Gird

JosephPlessis — Mon, 06 Feb 2017 18:52:32 GMT

I like this idea!

So after i could split aprt the mosiac using split raster?

Re: How to clip map sheets to a Map Gird

AbdullahAnter — Mon, 06 Feb 2017 19:20:12 GMT

Prefect Answer,Neil Ayres . It is it

Re: How to clip map sheets to a Map Gird

JosephPlessis — Mon, 06 Feb 2017 20:53:03 GMT

so far i have this but no return on my clip

import arcpy

# Set the current workspace
arcpy.env.workspace = "Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Mapsheets\\soil"

#variables
R1 = "102I16"
R2 = "102I09"
R3 = "102I08"
R4 = "092L13"
R5 = "092L12"
R6 = "092L11"
R7 = "092L10"
R8 = "092L08"
R9 = "092L07"
R10 = "092L06"
R11 = "092L05"
R12 = "092L04"
R13 = "092L03"
R14 = "092L02"
R15 = "092L01"
R16 = "092E16"
R17 = "092E15"
R18 = "092E14"
R19 = "092E10"
R20 = "092E09"
R21 = "092E08"
R22 = "092E07"

# Local variables:
Output_Raster ="Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Geodatabases\\Grid_points.gdb\\gridlist.tif"

gridlist= R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,R11,R12,R13,R14,R15,R16,R17,R18,R19,R20,R21,R22

grid="Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Geodatabases\\Grid_points.gdb\\GRID"

arcpy.MakeFeatureLayer_management(grid, "gridlyr")
arcpy.SelectLayerByAttribute_management("gridlyr", "NEW_SELECTION", " NTS_SNRC = 'gridlist' " )
print gridlist

# Get and print a list of TIFs from the workspace
rasters = arcpy.ListRasters("*", "TIF")
for raster in rasters:

print(raster)

while raster==gridlist:
arcpy.Clip_management("raster", Output_Raster, "gridlist", "", "", "NONE", "NO_MAINTAIN_EXTENT")

Re: How to clip map sheets to a Map Gird

NeilAyres — Tue, 07 Feb 2017 10:05:31 GMT

Joe,

can you try and format your code to make it readable...

script-formatting

The point about making a mosaic is that you don't have to do all this clipping.

No extra time and no extra disk space.

Re: How to clip map sheets to a Map Gird

JosephPlessis — Fri, 10 Dec 2021 20:05:18 GMT

THIS WORKS!

import arcpy
from arcpy import env
from arcpy.sa import *
arcpy.CheckOutExtension("3D")
arcpy.env.overwriteOutput = True

# Set the current workspace
arcpy.env.workspace = "Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Mapsheets\\soil"
Dir = "Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Mapsheets\\ClippedMaps"


GRID="Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Geodatabases\\Grid_points.gdb\\GRID"

GRIDlyr=arcpy.MakeFeatureLayer_management(GRID, "GRIDlyr")


# Get and print a list of TIFs from the workspace

rasters = arcpy.ListRasters("*")
for raster in rasters:
    print raster
    name = Dir + "\\" +"clp_"+ str(raster)
    print name
    arcpy.RasterDomain_3d(raster, "in_memory/temp", "POLYGON")

    # Process: Central Feature
    arcpy.MeanCenter_stats("in_memory/temp", "in_memory/temp1")

    # Process: Select Layer By Location
    arcpy.SelectLayerByLocation_management(GRIDlyr, "intersect", "in_memory/temp1", "", "NEW_SELECTION")

    # Clip Raster
    arcpy.Clip_management(raster,"#", name,GRIDlyr, "#", "ClippingGeometry")

    # Delete in_memory
    arcpy.Delete_management("in_memory")

    print "processing " + raster + " complete..."


    print "All processing is now finished"


‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍