<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
I am simply trying to open a raster image, manipulate in NumPy and then spit<BR>it back out in GDAL.<BR><BR>I am having trouble setting null data...<BR>I tried: ds.GetRasterBand(1).SetNoDataValue( -9999 )<BR>but nothing seems to work.<BR><BR>Any help is appreciated!<BR><BR><BR>Below is my code:<BR><BR><BR>#! /usr/bin/env python <BR><BR>import os, sys<BR>use_numeric = True<BR><BR>try:<BR>from osgeo import ogr, gdal<BR>from osgeo.gdalconst import *<BR>import numpy<BR>os.chdir(r'L:\users\gv\numPyGDAL\ospy_data4\ospy_data4')<BR>use_numeric = False<BR><BR>except ImportError:<BR>import ogr, gdal<BR>from gdalconst import *<BR>import Numeric<BR>os.chdir(r'L:\users\gv\numPyGDAL\ospy_data4\ospy_data4')<BR><BR><BR># register all of the drivers<BR>gdal.AllRegister()<BR><BR># open the image<BR>ds = gdal.Open('sequest', GA_ReadOnly)<BR>if ds is None:<BR>print 'Could not open image'<BR>sys.exit(1)<BR><BR>band = ds.GetRasterBand(1))# 1-based index<BR># read data and add the value to the string<BR>data = band.ReadAsArray()<BR>print data<BR>dims = data.shape <BR>nx=dims[1] <BR>ny=dims[0]<BR>nxarray=numpy.zeros(nx)<BR><BR>dst_format = 'HFA' <BR>dst_datatype = gdal.GDT_Int16<BR>dst_options = ['COMPRESS=LZW'] <BR>dst_file='sequest_final.img' <BR>dst_xsize = nx <BR>dst_ysize = ny <BR>dst_nbands = 1 <BR><BR>geoTransform = ds.GetGeoTransform()<BR>proj = ds.GetProjection()<BR><BR>driver = gdal.GetDriverByName( dst_format ) <BR>dst_ds = driver.Create(dst_file, dst_xsize, dst_ysize, dst_nbands, <BR>dst_datatype, dst_options)<BR><BR>for j in range(ny):<BR>nxarray = data[j,:]<BR>nxarray.shape = (1,nx)<BR>print nxarray<BR>dst_ds.GetRasterBand(1).WriteArray(nxarray,0,j)<BR>dst_ds.SetGeoTransform(geoTransform)<BR>dst_ds.SetProjection(proj)<BR>dst_ds = None <br /><hr />Hotmail: Free, trusted and rich email service. <a href='http://clk.atdmt.com/GBL/go/201469228/direct/01/' target='_new'>Get it now.</a></body>
</html>