I have a point shapefile, with no associated .prj file, that I have to rasterize. Each point should be transformed into a pixel of 30x30 size aligned to the pixels of a Landsat image (the point file contains accuracy assessment points for a land cover semi-automated classification).
I don't know the original coordinate system but I'm reasonably sure of which one it is. The steps I took are: 1) tried saving the shp with different coordinate system, compared the results, and decided which might be the most likely one (because in the other cases the land cover classes assigned to the points would be mostly wrong). 2) saved the shp in that coordinate system, then converted into EPSG:32632 (same as the landsat images) 3) used gdal rasterize (from Qgis' Processing Toolbox) to rasterize the points. I set "output resolution in map units per pixel" as 30 x 30 m.
The green pixel is the result of rasterization, the yellow point is the original point.
I tried redoing the whole procedure saving the original shp in different coordinate systems, but no luck. I also tried the plugin rasmover to manually align the raster, but then it occurred to me that I don't know for sure how the raster should be aligned (i mean, I guess it should be aligned to the pixel it overlaps the most, but how can I be certain?). Moreover , if I use rasmover, many of the original points end up falling outside of the pixel that is supposed to contain them, because they are located near the edge.
So the questions are, what might be causing this problem? Is it related to the original georeferencing of the two files? And is there a way to fix this problem by using QGis, or by setting gdal_rasterize so that it will align the raster correctly? I've seen answers to similar problems related to ArcGis and GRASS but not specifically QGIS.