when viewing rotated rasters in QGIS - the overall raster is rotated (as designed), though each grid cell is not. (the nodata value doesn't appear to pick up as the default in QGIS either - -2.9999e+30 instead of -3e+30, but that might be a separate question).
.
and a closeup
is this common for QGIS? viewing the same in 'the other GIS GUI' appears to rotate each cell.
i can add all the code - but the basic part is
def writeRast(gridbound, outTiff, grid, cellsize=1000, wkt=''): ''' lots of other stuff''' geotransform=(ulx, xres, xtheta, uly, ytheta, -yres) drv = gdal.GetDriverByName('GTiff') ds = drv.Create(outTiff, ncol, nrow, 1 ,gdal.GDT_Float32) # Open the file band = ds.GetRasterBand(1) band.SetNoDataValue(-3e30) ds.SetGeoTransform(geotransform) # Specify its coordinates if wkt != '': targproj = osr.SpatialReference() targproj.ImportFromWkt(wkt) band.WriteArray(grid) is there some setting in QGIS that would display each cell rotated? or something i'm missing in writing the raster with GDAL?
tried with gdal 1.11.0 and 2.0.on mac OSX, QGIS v 2.8.3.
أكثر...
.and a closeup
is this common for QGIS? viewing the same in 'the other GIS GUI' appears to rotate each cell.
i can add all the code - but the basic part is
def writeRast(gridbound, outTiff, grid, cellsize=1000, wkt=''): ''' lots of other stuff''' geotransform=(ulx, xres, xtheta, uly, ytheta, -yres) drv = gdal.GetDriverByName('GTiff') ds = drv.Create(outTiff, ncol, nrow, 1 ,gdal.GDT_Float32) # Open the file band = ds.GetRasterBand(1) band.SetNoDataValue(-3e30) ds.SetGeoTransform(geotransform) # Specify its coordinates if wkt != '': targproj = osr.SpatialReference() targproj.ImportFromWkt(wkt) band.WriteArray(grid) is there some setting in QGIS that would display each cell rotated? or something i'm missing in writing the raster with GDAL?
tried with gdal 1.11.0 and 2.0.on mac OSX, QGIS v 2.8.3.
أكثر...