<div dir="ltr">Thanks for that information Regina. Your right, ST_Colormap does retain the projection information. That's interesting as I have been able to retain the spatial information in a jpeg using gdal_translate without creating the jpegw file.<div><br></div><div>Thanks</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Apr 26, 2016 at 12:48 AM, Regina Obe <span dir="ltr"><<a href="mailto:lr@pcorp.us" target="_blank">lr@pcorp.us</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" link="blue" vlink="purple"><div><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">David,<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">I could be mistaken, but I always thought that for jpegs the spatial info is stored in a .jpegw file and since ST_ASJPEG only outputs the jpg part it wouldn't keep the spatial ref info.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Are you saying you loose it at the ST_ColorMap part or the ST_AsJPEG.  I would expect it maintained in ST_ColorMap and lost at ST_AsJPEG call.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Thanks,<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Regina<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p><p class="MsoNormal" style="margin-left:.5in"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> postgis-devel [mailto:<a href="mailto:postgis-devel-bounces@lists.osgeo.org" target="_blank">postgis-devel-bounces@lists.osgeo.org</a>] <b>On Behalf Of </b>David Haynes II<br><b>Sent:</b> Friday, April 22, 2016 12:51 PM<br><b>To:</b> PostGIS Development Discussion <<a href="mailto:postgis-devel@lists.osgeo.org" target="_blank">postgis-devel@lists.osgeo.org</a>><br><b>Subject:</b> [postgis-devel] ST_Colormap<u></u><u></u></span></p><p class="MsoNormal" style="margin-left:.5in"><u></u> <u></u></p><div><p class="MsoNormal" style="margin-left:.5in">Hello<u></u><u></u></p><div><p class="MsoNormal" style="margin-left:.5in"><u></u> <u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">Working on an application for clipping rasters (<a href="http://data.terrapop.org/terraclip" target="_blank">data.terrapop.org/terraclip</a>) and running into an issue with the ST_Colormap function. It seems that whenever we run the ST_Colormap function we lose the spatial reference information for the image.This is the sample of the query that we are running. The query is successful, but the resulting images can be built into a vrt with gdal. I'm using PostgreSQL 9.3.7  and POSTGIS 2.1.7.<u></u><u></u></p></div><div><div><p class="MsoNormal" style="margin-left:.5in"><u></u> <u></u></p></div><div><div><p class="MsoNormal" style="margin-left:.5in">WITH projection as<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">(<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">SELECT ST_SRID(rast) as srid<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">FROM gis_rasters.apricothar<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">Limit 1<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">),polygon as<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">(<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">SELECT <a href="http://sgl.id" target="_blank">sgl.id</a> as sample_geog_level_id, <a href="http://gi.id" target="_blank">gi.id</a> as geog_instance_id, gi.label as geog_instance_label,<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">gi.code as geog_instance_code, ST_Transform(bound.geog::geometry, prj.srid) as geom<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">FROM projection prj, sample_geog_levels sgl<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">inner join geog_instances gi on <a href="http://sgl.id" target="_blank">sgl.id</a> = gi.sample_geog_level_id<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">inner join boundaries bound on bound.geog_instance_id = <a href="http://gi.id" target="_blank">gi.id</a><u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">WHERE <a href="http://sgl.id" target="_blank">sgl.id</a> = 226<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">),raster_clip as<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">(<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">SELECT ST_Clip(r.rast, 1, p.geom, TRUE) AS rast <u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">FROM polygon p inner join gis_rasters.apricothar  r on ST_Intersects(r.rast,p.geom)<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">)<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">select ST_AsJPEG(ST_Colormap(rast, 1, 'greyscale'))  as img<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">from raster_clip, projection prj<u></u><u></u></p></div></div><div><p class="MsoNormal" style="margin-left:.5in"><u></u> <u></u></p></div><p class="MsoNormal" style="margin-left:.5in">-- <u></u><u></u></p><div><div><div><div><p class="MsoNormal" style="margin-left:.5in">David Haynes, Ph.D.<u></u><u></u></p><div><p class="MsoNormal" style="margin-left:.5in">Research Associate Terra Populus<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in">Minnesota Population Center<u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in"><a href="tel:612-626-3933" value="+16126263933" target="_blank">612-626-3933</a><u></u><u></u></p></div><div><p class="MsoNormal" style="margin-left:.5in"><a href="http://www.terrapop.org" target="_blank">www.terrapop.org</a><u></u><u></u></p></div></div></div></div></div></div></div></div></div><br>_______________________________________________<br>
postgis-devel mailing list<br>
<a href="mailto:postgis-devel@lists.osgeo.org">postgis-devel@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/postgis-devel" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/postgis-devel</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr">David Haynes, Ph.D.<div>Research Associate Terra Populus</div><div>Minnesota Population Center</div><div>612-626-3933</div><div><a href="http://www.terrapop.org" target="_blank">www.terrapop.org</a></div></div></div></div></div>
</div>