<div dir="ltr">Hi Jukka,<div><br></div><div>I have no choice but to use the standard InspireCRS84Quad tiling grid.</div><div><br></div><div>Thanks for the "-a_ullr fix"</div><div><br></div><div>Philippe</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 18, 2024 at 2:42 PM Rahkonen Jukka <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi">jukka.rahkonen@maanmittauslaitos.fi</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"><div class="msg-6343131382266751906">
<div lang="FI" style="overflow-wrap: break-word;">
<div class="m_4603536328409217765WordSection1">
<p class="MsoNormal"><span>Hi,<u></u><u></u></span></p>
<p class="MsoNormal"><span><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Can you imagine using a custom tiling schema instead of the world-wide
</span>InspireCRS84Quad<span lang="EN-US">? These commands work<br>
<br>
gdal_create -of GTiff -co COMPRESS=DEFLATE -outsize 2000 2000 -burn 255 -burn 0 -burn 0 -bands 3 -ot Byte -a_srs epsg:32631 -a_ullr 1000000 1000002 1000002 1000000 1mm.tif<br>
<br>
gdal_translate -of GPKG 1mm.tif 1mm.gpkg<br>
<br>
Notice the fixed ullr for making the image north-up.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">-Jukka Rahkonen-<br>
<br>
<br>
<br>
gdal_translate -of GPKG 1mm.tif 1mm.gpkg<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">Input file size is 2000, 2000<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">0ERROR 6: Only north-up non rotated geotransform supported<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">..ERROR 1: Raster table 1mm not correctly initialized due to missing call to SetGeoTransform()<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(225,225,225);padding:3pt 0cm 0cm">
<p class="MsoNormal"><b>Lähettäjä:</b> gdal-dev <<a href="mailto:gdal-dev-bounces@lists.osgeo.org" target="_blank">gdal-dev-bounces@lists.osgeo.org</a>>
<b>Puolesta </b>Philippe Ghesquiere via gdal-dev<br>
<b>Lähetetty:</b> torstai 18. tammikuuta 2024 15.22<br>
<b>Vastaanottaja:</b> gdal-dev <<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a>><br>
<b>Aihe:</b> [gdal-dev] GPKG raster creation : error with too small resolution<u></u><u></u></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">Dear all,<br>
<br>
I just had an error while converting an image to GPKG. This image has a very small resolution and gdal_translate pops up an error message.<br>
<br>
<u>Test image :</u><br>
gdal_create -of GTiff -co COMPRESS=DEFLATE -outsize 2000 2000 -burn 255 -burn 0 -burn 0 -bands 3 -ot Byte -a_srs epsg:32631 -a_ullr 1000000 1000000 1000002 1000002 1mm.tif<br>
<br>
<u>Conversion to GPKG :</u><br>
gdal_translate -of GPKG -co TILING_SCHEME=InspireCRS84Quad -co ZOOM_LEVEL_STRATEGY=UPPER 1mm.tif 1mm.gpkg<br>
Input file size is 2000, 2000<br>
ERROR 1: Could not find an appropriate zoom level<br>
<br>
<u>Analysis</u><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">The error is due to <a href="https://github.com/OSGeo/gdal/blob/master/ogr/ogrsf_frmts/gpkg/ogrgeopackagedatasource.cpp#L3169" title="Follow link" target="_blank">
ogrgeopackagedatasource.cpp</a> (line 3169), and it does make sense : the zoom level is too high (Z>25).<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">However, wouldn't it be possible to modify gdal_translate behavior ?<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">For example, we could :<u></u><u></u></p>
</div>
<div>
<ul type="disc">
<li class="MsoNormal">
let GDAL compute the appropriate zoom level : -co ZOOM_LEVEL_STRATEGY=AUTO/LOWER/UPPER,<u></u><u></u></li><li class="MsoNormal">
optionally specify a maximum zoom level : -co MAX_ZOOM_LEVEL=18<u></u><u></u></li></ul>
<div>
<p class="MsoNormal">What do you think ?<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Philippe<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<p class="MsoNormal"><span style="font-size:7pt">The information in this e-mail is confidential. The contents may not be disclosed or used by anyone other than the addressee. Access to this e-mail by anyone else is unauthorised.<br>
If you are not the intended recipient, please notify Airbus immediately and delete this e-mail.<br>
Airbus cannot accept any responsibility for the accuracy or completeness of this e-mail as it has been sent over public networks. If you have any concerns over the content of this message or its Accuracy or Integrity, please contact Airbus immediately.<br>
All outgoing e-mails from Airbus are checked using regularly updated virus scanning software but you should take whatever measures you deem to be appropriate to ensure that this message and any attachments are virus free.</span><u></u><u></u></p>
</div>
</div>
</div></blockquote></div>
<font style="font-size: 9px;">The information in this e-mail is confidential. The contents may not be disclosed or used by anyone other than the addressee. Access to this e-mail by anyone else is unauthorised.<br>If you are not the intended recipient, please notify Airbus immediately and delete this e-mail.<br>Airbus cannot accept any responsibility for the accuracy or completeness of this e-mail as it has been sent over public networks. If you have any concerns over the content of this message or its Accuracy or Integrity, please contact Airbus immediately.<br>All outgoing e-mails from Airbus are checked using regularly updated virus scanning software but you should take whatever measures you deem to be appropriate to ensure that this message and any attachments are virus free.</font>