Skip to content

Commit

Permalink
Merge pull request #17 from debboutr/master
Browse files Browse the repository at this point in the history
update get_geo_info namespace in align_rasters
  • Loading branch information
ozak committed Jan 31, 2017
2 parents 63f6402 + b7bba03 commit 9572842
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions georasters/georasters.py
Original file line number Diff line number Diff line change
Expand Up @@ -148,8 +148,8 @@ def align_rasters(raster, alignraster, how=np.ma.mean, cxsize=None, cysize=None,
how: function used to aggregate cells (if the rasters have different sizes)
It is assumed that both rasters have the same size
'''
ndv1, xsize1, ysize1, geot1, projection1, datatype1 = GetGeoInfo(raster)
ndv2, xsize2, ysize2, geot2, projection2, datatype2 = GetGeoInfo(alignraster)
ndv1, xsize1, ysize1, geot1, projection1, datatype1 = get_geo_info(raster)
ndv2, xsize2, ysize2, geot2, projection2, datatype2 = get_geo_info(alignraster)
if projection1.ExportToMICoordSys() == projection2.ExportToMICoordSys():
blocksize = (np.round(geot2[1]/geot1[1]), np.round(geot2[-1]/geot1[-1]))
mraster = gdalnumeric.LoadFile(raster)
Expand Down

0 comments on commit 9572842

Please sign in to comment.