<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p>cf improvement <a class="moz-txt-link-freetext" href="https://github.com/OSGeo/gdal/pull/9101">https://github.com/OSGeo/gdal/pull/9101</a></p>
<p><a class="moz-txt-link-freetext" href="https://docs.ogc.org/is/17-083r2/17-083r2.html#76">https://docs.ogc.org/is/17-083r2/17-083r2.html#76</a> mentions "TileMatrices
for more detailed scales is possible if they follow the same
pattern." . The current limitation to 25 was to make some
computations easier, but supporting up to 30 (or 29 when
MatrixWidth=2 at zoom level = 0) is now possible.<br>
</p>
<div class="moz-cite-prefix">Le 18/01/2024 à 15:35, Rahkonen Jukka
via gdal-dev a écrit :<br>
</div>
<blockquote type="cite"
cite="mid:DB9PR09MB6887AC9FB67AB80EA90B72D0FD712@DB9PR09MB6887.eurprd09.prod.outlook.com">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<meta name="Generator"
content="Microsoft Word 15 (filtered medium)">
<style>@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}span.Shkpostityyli18
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}div.WordSection1
{page:WordSection1;}ol
{margin-bottom:0cm;}ul
{margin-bottom:0cm;}</style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="WordSection1">
<p class="MsoNormal"><span style="mso-fareast-language:EN-US">Hi,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-US">I am not sure how standard the level 25 is for
the
</span><span lang="EN-US">InspireCRS84Quad, the INSPIRE
technical guidance stops at 17. But I believe that you know
your own use case and have considered also that EPSG:3857 at
1 mm resolution is good for you</span><span
style="mso-fareast-language:EN-US" lang="EN-US">.<br>
<br>
The number of tiles at level 25 is 2^50=</span><span
style="color:black" lang="EN-US"> 1.1259E+15. TileCol and
TileRow would get rather large numbers but I suppose that it
is not an issue.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-US">-Jukka Rahkonen-<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-US"
lang="EN-US"><o:p> </o:p></span></p>
<div
style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b>Lähettäjä:</b> Philippe Ghesquiere
<a class="moz-txt-link-rfc2396E" href="mailto:philippe.ghesquiere@airbus.com"><philippe.ghesquiere@airbus.com></a>
<br>
<b>Lähetetty:</b> torstai 18. tammikuuta 2024 16.06<br>
<b>Vastaanottaja:</b> Rahkonen Jukka
<a class="moz-txt-link-rfc2396E" href="mailto:jukka.rahkonen@maanmittauslaitos.fi"><jukka.rahkonen@maanmittauslaitos.fi></a><br>
<b>Kopio:</b> gdal-dev <a class="moz-txt-link-rfc2396E" href="mailto:gdal-dev@lists.osgeo.org"><gdal-dev@lists.osgeo.org></a><br>
<b>Aihe:</b> Re: [gdal-dev] GPKG raster creation : error
with too small resolution<o:p></o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Hi Jukka,<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I have no choice but to use the
standard InspireCRS84Quad tiling grid.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Thanks for the "-a_ullr fix"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Philippe<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Thu, Jan 18, 2024 at 2:42 PM
Rahkonen Jukka <<a
href="mailto:jukka.rahkonen@maanmittauslaitos.fi"
moz-do-not-send="true" class="moz-txt-link-freetext">jukka.rahkonen@maanmittauslaitos.fi</a>>
wrote:<o:p></o:p></p>
</div>
<blockquote
style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<div>
<div>
<div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Hi,<o:p></o:p></p>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><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.</span><o:p></o:p></p>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span
lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span
lang="EN-US">-Jukka Rahkonen-<br>
<br>
<br>
<br>
gdal_translate -of GPKG 1mm.tif 1mm.gpkg</span><o:p></o:p></p>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span
lang="EN-US">Input file size is 2000, 2000</span><o:p></o:p></p>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span
lang="EN-US">0ERROR 6: Only north-up non rotated
geotransform supported</span><o:p></o:p></p>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span
lang="EN-US">..ERROR 1: Raster table 1mm not
correctly initialized due to missing call to
SetGeoTransform()</span><o:p></o:p></p>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span
lang="EN-US"> </span><o:p></o:p></p>
<div
style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b>Lähettäjä:</b>
gdal-dev <<a
href="mailto:gdal-dev-bounces@lists.osgeo.org"
target="_blank" moz-do-not-send="true"
class="moz-txt-link-freetext">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" moz-do-not-send="true"
class="moz-txt-link-freetext">gdal-dev@lists.osgeo.org</a>><br>
<b>Aihe:</b> [gdal-dev] GPKG raster creation :
error with too small resolution<o:p></o:p></p>
</div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">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><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">The error is
due to
<a
href="https://github.com/OSGeo/gdal/blob/master/ogr/ogrsf_frmts/gpkg/ogrgeopackagedatasource.cpp#L3169"
target="_blank" title="Follow link"
moz-do-not-send="true">
ogrgeopackagedatasource.cpp</a> (line 3169),
and it does make sense : the zoom level is too
high (Z>25).<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">However,
wouldn't it be possible to modify gdal_translate
behavior ?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">For example,
we could :<o:p></o:p></p>
</div>
<div>
<ul type="disc">
<li class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;mso-list:l0 level1 lfo3">
let GDAL compute the appropriate zoom level :
-co ZOOM_LEVEL_STRATEGY=AUTO/LOWER/UPPER,<o:p></o:p></li>
<li class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;mso-list:l0 level1 lfo3">
optionally specify a maximum zoom level : -co
MAX_ZOOM_LEVEL=18<o:p></o:p></li>
</ul>
<div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">What do you
think ?<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Philippe<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"
style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span
style="font-size:7.0pt">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><o:p></o:p></p>
</div>
</div>
</div>
</blockquote>
</div>
<p class="MsoNormal"><span style="font-size:7.0pt">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><o:p></o:p></p>
</div>
<br>
<fieldset class="moz-mime-attachment-header"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
gdal-dev mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/gdal-dev">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a>
</pre>
</blockquote>
<pre class="moz-signature" cols="72">--
<a class="moz-txt-link-freetext" href="http://www.spatialys.com">http://www.spatialys.com</a>
My software is free, but my time generally not.</pre>
</body>
</html>