Clipping Temporal raster

148
1
04-16-2024 05:29 AM
Colome
by
New Contributor II

As part of a script I am running, I am clipping the  Sentinel-2 10m Land Use/Land Cover Time Series data from ESRI (https://www.arcgis.com/home/item.html?id=cfcb7609de5f478eb7666240902d4d3d) using Python. I was using the previous version no problem, however, with this version of the data, when I do my clip, as default it clips the oldest version of the raster, whereas I actually want the latest / most recent version of the land cover. Does anyone know if its possible to specify this in Python? 

 
        in_rst_landcover = "in_rst_landcover"
        arcpy.MakeRasterLayer_management(in_raster=esri_lulc,
                                         out_rasterlayer=in_rst_landcover,
                                         envelope=in_fl_extent
                                         )
        der_rst_step1a = "step1a.tif"
        arcpy.management.Clip(in_raster=in_rst_landcover,
                              rectangle="",
                              out_raster=der_rst_step1a,
                              in_template_dataset=in_fl_extent,
                              nodata_value='',
                              clipping_geometry="ClippingGeometry",
                              maintain_clipping_extent="NO_MAINTAIN_EXTENT")
0 Kudos
1 Reply
DuncanHornby
MVP Notable Contributor

Just an idea when you run the make raster layer tool may be you should set the band number to what would be the most recent date?

0 Kudos