topic Read specific raster file using python in Python Questions

Read specific raster file using python

RameshGautam — Tue, 17 Jan 2012 20:26:42 GMT

I am trying to read a specific raster image file using python, however, getting some difficulties. Can you please help me?

Here is the details about the files:

import arcpy
import sys, os, string
from arcpy.sa import *
from arcpy import env
arcpy.env.workspace = "D:\Project\Stanislaus\Timeseries\NDVI"
LstRst = arcpy.ListRasters("*", "img")
for Rst in LstRst:
    Img1 = string.find(Rst, "_stanislaus.img")
    if Img1 > -1:
        print Rst


when I print Rst, it prints all raster images with _stanislaus.img. Now, I need to read an specific image like "aug08_2010_stanislaus.img". Can you please tell me how to do that?

Thanks

Re: Read specific raster file using python

DarrenWiens2 — Sat, 11 Dec 2021 10:40:23 GMT

I'm confused. If you know the file name, why not just substitute the file name for "_stanislaus.img"?

import arcpy
 import sys, os, string
 from arcpy.sa import *
 from arcpy import env
 arcpy.env.workspace = "D:\Project\Stanislaus\Timeseries\NDVI"
 LstRst = arcpy.ListRasters("*", "img")
 for Rst in LstRst:
     Img1 = string.find(Rst, "aug08_2010_stanislaus.img")
     if Img1 > -1:
     print Rst

Or, are you looking for how to substitute a variable into the file name?

import arcpy
 import sys, os, string
 from arcpy.sa import *
 from arcpy import env
 arcpy.env.workspace = "D:\Project\Stanislaus\Timeseries\NDVI"
 LstRst = arcpy.ListRasters("*", "img")
 fileprefix = "aug08_2010"
 for Rst in LstRst:
     Img1 = string.find(Rst, fileprefix + "_stanislaus.img")
     if Img1 > -1:
     print Rst

Re: Read specific raster file using python

RameshGautam — Sat, 11 Dec 2021 10:40:26 GMT

Thanks a lot for your prompt help. When I replace "_stanislaus.img" by "aug08_2010_stanislaus.img", it can't print the file, however, as you have mentioned in second option by using fileprefix, we can read this file. Do you know why?

Thanks again for your help. I apprecaite it.

I'm confused. If you know the file name, why not just substitute the file name for "_stanislaus.img"?

import arcpy
 import sys, os, string
 from arcpy.sa import *
 from arcpy import env
 arcpy.env.workspace = "D:\Project\Stanislaus\Timeseries\NDVI"
 LstRst = arcpy.ListRasters("*", "img")
 for Rst in LstRst:
     Img1 = string.find(Rst, "aug08_2010_stanislaus.img")
     if Img1 > -1:
     print Rst

Or, are you looking for how to substitute a variable into the file name?

import arcpy
 import sys, os, string
 from arcpy.sa import *
 from arcpy import env
 arcpy.env.workspace = "D:\Project\Stanislaus\Timeseries\NDVI"
 LstRst = arcpy.ListRasters("*", "img")
 fileprefix = "aug08_2010"
 for Rst in LstRst:
     Img1 = string.find(Rst, fileprefix + "_stanislaus.img")
     if Img1 > -1:
     print Rst