<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto">Hi,<div><br><div>You can also correct this 5x5 problem by using gdal grid with Nearest neighbour interpolation. You could use the point data to do this.</div><div><br></div><div>I also post a small presentation that can help go from vector to raster with point data. Sometimes, keeping it simple is better.</div><div><br></div><div><a href="https://www.slideshare.net/shencoop/qgis-raster-to-point">https://www.slideshare.net/shencoop/qgis-raster-to-point</a></div><div><br></div><div>If you do decide to make a shapefile with the csv, make sure to store the coordinates in a text field. Placing the coordinates in a $x and $y object fields will lead to rounding. Any rounding would corrupt the raster therefore you need to original coordinates for this to work. This may explain your problems also.</div><div><br></div><div>I used this method to change the geoids of a thousands of rasters because the program I was using was designed to change the z value for texte files only and not rasters (NRCAN GPS-H).<br><br><div dir="ltr">Nicolas Cadieux<div>Ça va bien aller!</div></div><div dir="ltr"><br><blockquote type="cite">Le 13 mai 2020 à 11:25, chris hermansen <clhermansen@gmail.com> a écrit :<br><br></blockquote></div><blockquote type="cite"><div dir="ltr"><meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr"><div dir="ltr">Maria and list,<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, May 12, 2020 at 7:00 PM Priv.-Doz. Dr. Maria Shinoto <<a href="mailto:maria.shinoto@zaw.uni-heidelberg.de">maria.shinoto@zaw.uni-heidelberg.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi, <br>
<br>
Thanks for helping. -- Thanks to you and Chris Hermansen I got a result, but it could be better. <br>
<br>
For the records, a short explanation:<br>
<br>
*****<br>
Well, I checked the properties, jgd2011 is in Meters, the raster is said to be 5m. In the official Japanese viewer, which creates a beautiful raster image without white pixels, the pixels are exactly 5m*5m. <br>
<br>
Today I tried the export to .xyz since the shapefile looked ugly, and after realising that the Japanese xyz is indeed yxz, everything looked fine, and I could store it in a Geopackage. But the grid is now 5,276m * 6,146m. But it fits well on top of the basemap. The basemap is of the same special Japanese GML format, but QGIS could read it all without problem. I do not understand why QGIS does not read the point data from the GML fille, but that is an aside, I am amazed by what QGIS actually can do. <br>
<br>
From the Geopackage I could rasterize. It is as Chris Hermansen said, thanks. Unfortunately, I did not get it done from the shapefiles, they always looked weird or like nothing, even with identical settings. But the geopackages from xyz tiles are fine. <br>
<br>
For resolution, I chose georeferenced units as Chris suggested, and since the measurement tool got some different length, I put it to 5,276m by 6,146m. A 5m by 5m resolution created a weird layer with horizontally expanding white pixels.<br>
<br>
It seems that tweaking with the resolution might lead to an even better result, but for the time being, it is OK as it is. <br>
*****<br>
<br>
<br></blockquote><div>Upon reflection I think the basic problem here is that the point data should be interpolated to create a raster if you want a precise 5x5m resolution.<br></div><div><br></div><div>For this, rather than use the Raster > Rasterize tool, the approach should be:</div><div><ol><li>open the processing toolbox Processing > Toolbox</li><li>in the toolbox open Interpolation > TIN interpolation</li><li>in the TIN Interpolation screen:</li><ol><li>select the Vector layer</li><li>select the Interpolation attribute</li><li>click the + to add to the vector layer panel</li><li>choose the interpolation method - probably best to use cubic<br></li><li>click on the ... next to extent and set it to the layer extent</li><li>set the pixel size to 5.0 and 5.0</li><li>click Run</li></ol></ol><div>This way you won't have the odd sizes you mentioned. This may give you a smoother surface in the end as well.<br></div></div></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr">Chris Hermansen · clhermansen "at" gmail "dot" com<br><br>C'est ma façon de parler.</div></div></div>
</div></blockquote></div></div></body></html>