Even Rouault, Thank you! Now the program is working.
On Mon, Nov 24, 2014 at 2:15 AM, Even Rouault
wrote:
> Le dimanche 23 novembre 2014 17:04:05, Simen Langseth a écrit :
> > import gdal, numpy as np
> > from gdalconst import *
> >
> > infile = r"D:\given.tif"
> > inraster = gdal.Open(infile, G
Le dimanche 23 novembre 2014 17:04:05, Simen Langseth a écrit :
> import gdal, numpy as np
> from gdalconst import *
>
> infile = r"D:\given.tif"
> inraster = gdal.Open(infile, GA_ReadOnly)
> x_total,y_total=refraster.RasterXSize,refraster.RasterYSize
> n_chunks=10
> x_offsets=np.linspace(0,x_tota
Wong,
On Fri, Nov 21, 2014 at 12:33 PM, Wong Hua wrote:
> Hello, I am a osgeo4w user and very green in the field of geo data.
>
> I am using gdal_grid to grid and ungridded set of points (a CSV file).
>
> I ultimately want to use it with gdaltransform to change the projection type
> and gdaltrans
import gdal, numpy as np
from gdalconst import *
infile = r"D:\given.tif"
inraster = gdal.Open(infile, GA_ReadOnly)
x_total,y_total=refraster.RasterXSize,refraster.RasterYSize
n_chunks=10
x_offsets=np.linspace(0,x_total,n_chunks).astype(int)
x_offsets=zip(x_offsets[:-1],x_offsets[1:])
y_offsets=np
