[prev in list] [next in list] [prev in thread] [next in thread] 

List:       mapserver-users
Subject:    [mapserver-users] Faster tileindex layer with reverse painter's algorithm
From:       "Rahkonen Jukka (MML)" <jukka.rahkonen () maanmittauslaitos ! fi>
Date:       2018-11-22 11:25:29
Message-ID: 5dab876abdd64dbfbb353b3e1196f739 () C119S212VM042 ! msvyvi ! vaha ! local
[Download RAW message or body]

Hi,

Currently raster tileindex layers are rendered according to painter's algor=
ithm. First image that it found from the tileindex is rendered first, then =
pixels from the next image are added by method "b over a" which is the defa=
ult https://mapserver.org/uk/development/rfc/ms-rfc-113.html#rfc113 until t=
he last image is processed. The end result is fine and if user wants to get=
 a result that is composed from the newest available images overall it can =
be done by reading the tileindex sorted by time "select * from tileindex or=
der by time".

Unfortunately the service gets pretty slow when the tileindex contains lots=
 of overlapping images from different times because every image that gets s=
elected by the bounding box gets also rendered even if all the pixels from =
some image will get overwritten later.  Setting a lower limit for time does=
 not resolve our problem because the newest image can be quite old in some =
places.

I would like to see suggestions about how to implement a fast pre-selection=
 that discards all the images which will not have any effect on the final r=
aster. Ideally Mapserver would do it natively if the default blend mode is =
used but SQL or function for PostGIS or Spatialite would suit me as well. I=
 have been thinking that some sort of pre-selection with the reverse painte=
r's method might do for this purpose: test the footprints of the images fro=
m top to bottom, skip images which are covered by the growing union of foot=
prints and stop when the whole BBOX is full.

-Jukka Rahkonen-


[Attachment #3 (text/html)]

<html xmlns:v="urn:schemas-microsoft-com:vml" \
xmlns:o="urn:schemas-microsoft-com:office:office" \
xmlns:w="urn:schemas-microsoft-com:office:word" \
xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" \
xmlns="http://www.w3.org/TR/REC-html40"> <head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@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;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
	{margin:0cm;
	margin-bottom:.0001pt;
	font-size:11.0pt;
	font-family:"Calibri",sans-serif;
	mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
	{mso-style-priority:99;
	color:#0563C1;
	text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
	{mso-style-priority:99;
	color:#954F72;
	text-decoration:underline;}
span.Shkpostityyli17
	{mso-style-type:personal-compose;
	font-family:"Calibri",sans-serif;
	color:windowtext;}
.MsoChpDefault
	{mso-style-type:export-only;
	font-family:"Calibri",sans-serif;
	mso-fareast-language:EN-US;}
@page WordSection1
	{size:612.0pt 792.0pt;
	margin:70.85pt 2.0cm 70.85pt 2.0cm;}
div.WordSection1
	{page:WordSection1;}
--></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]-->
</head>
<body lang="FI" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">Hi,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p>&nbsp;</o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Currently raster tileindex layers are \
rendered according to painter&#8217;s algorithm. First image that it found from the \
tileindex is rendered first, then pixels from the next image are added by method \
&#8220;b over a&#8221; which is  the default <a \
href="https://mapserver.org/uk/development/rfc/ms-rfc-113.html#rfc113"> \
https://mapserver.org/uk/development/rfc/ms-rfc-113.html#rfc113</a> until the last \
image is processed. The end result is fine and if user wants to get a result that is \
composed from the newest available images overall it can be done by reading the \
tileindex  sorted by time &#8220;select * from tileindex order by \
time&#8221;.<o:p></o:p></span></p> <p class="MsoNormal"><span \
lang="EN-US"><o:p>&nbsp;</o:p></span></p> <p class="MsoNormal"><span \
lang="EN-US">Unfortunately the service gets pretty slow when the tileindex contains \
lots of overlapping images from different times because every image that gets \
selected by the bounding box gets also rendered even if all the pixels  from some \
image will get overwritten later.&nbsp; Setting a lower limit for time does not \
resolve our problem because the newest image can be quite old in some \
places.<o:p></o:p></span></p> <p class="MsoNormal"><span \
lang="EN-US"><o:p>&nbsp;</o:p></span></p> <p class="MsoNormal"><span lang="EN-US">I \
would like to see suggestions about how to implement a fast pre-selection that \
discards all the images which will not have any effect on the final raster. Ideally \
Mapserver would do it natively if the default blend  mode is used but SQL or function \
for PostGIS or Spatialite would suit me as well. I have been thinking that some sort \
of pre-selection with the reverse painter&#8217;s method might do for this purpose: \
test the footprints of the images from top to bottom, skip images  which are covered \
by the growing union of footprints and stop when the whole BBOX is \
full.<o:p></o:p></span></p> <p class="MsoNormal"><span \
lang="EN-US"><o:p>&nbsp;</o:p></span></p> <p class="MsoNormal"><span \
lang="EN-US">-Jukka Rahkonen-<o:p></o:p></span></p> <p class="MsoNormal"><span \
lang="EN-US"><o:p>&nbsp;</o:p></span></p> </div>
</body>
</html>


[Attachment #4 (unknown)]

_______________________________________________
mapserver-users mailing list
mapserver-users@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/mapserver-users

[prev in list] [next in list] [prev in thread] [next in thread] 

Configure | About | News | Add a list | Sponsored by KoreLogic