5

I am trying to change the pixel value of a tif file (containing projection information) but after the pixel value is changed it is losing projection information.

I am using the Python imaging library for this and the code is like this 

import Image
im = Image.open("some tif file")
im.putpixel((2,3),0)

So please help me in solving this ..

whuber
  • 69,783
  • 15
  • 186
  • 281
Vicky
  • 1,852
  • 3
  • 22
  • 28

1 Answers1

6

You need a Python library that supports the GeoTIFF format. Libraries that support only the TIFF file discard the projection and geotransform metadata when they are saved. GDAL is the most used open source Python library for reading/writing spatially aware rasters. So, for your example:

from osgeo import gdal

# Load file, and access the band and get a NumPy array
src = gdal.Open('myfile.tif', gdal.GA_Update)
band = src.GetRasterBand(1)
ar = band.ReadAsArray()

# Assign a new raster pixel value
ar[2,3] = 0

# Save/close the file
band.WriteArray(ar)
del ar, band, src
Mike T
  • 42,095
  • 10
  • 126
  • 187