Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- for image in rasters:
- print(image)
- name, ext = os.path.splitext(image)
- if ext == '.hdr':
- img = sp.open_image(ssrc_directory + image)
- print(image)
- im_HS = img[:,:,0]
- cols = im_HS.shape[0] # Number of column
- rows = im_HS.shape[1] # Number of lines
- bands = im_HS.shape[2] # Number of bands
- names = []
- ##mask NDVI
- names.append('maskNDVI')
- for i in range(0,cols):
- for j in range(0,rows):
- if im_HS[i,j] >= 0.6 :
- im_HS[i,j] = np.min(im_HS)
- img = sp.open_image(ssrc_directory + image)
- print(image)
- im_HS = img[:,:,:]
- cols = im_HS.shape[0] # Number of column
- rows = im_HS.shape[1] # Number of lines
- bands = im_HS.shape[2] # Number of bands
- names = []
- ##mask NDVI
- names.append('mask')
- for h in range(0,bands):
- for i in range(0,cols):
- for j in range(0,rows):
- if im_HS[i,j,0] >= 0.6 :
- im_HS[i,j,h] = np.min(im_HS)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement