<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=utf-8"><meta name=Generator content="Microsoft Word 12 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"MS Mincho";
panose-1:2 2 6 9 4 2 5 8 3 4;}
@font-face
{font-family:"MS Mincho";
panose-1:2 2 6 9 4 2 5 8 3 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:"\@MS Mincho";
panose-1:2 2 6 9 4 2 5 8 3 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
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=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Thanks very much Chaitanya, that was indeed my problem.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>I am using gdal_grid to create the image that I am filtering, and am setting "nodata=255" as a parameter on the interpolation algorithm. I had expected that this would cause nodata areas to be written as 255, and that 255 would be marked as the nodata value in the output image (geotiff) metadata, but apparently only the former is true, and I need to run the output through gdal_translate as you indicate in order to set the metadata.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Thanks again,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>--Ken<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><div style='border:none;border-left:solid blue 1.5pt;padding:0in 0in 0in 4.0pt'><div><div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in'><p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>From:</span></b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'> gdal-dev-bounces@lists.osgeo.org [mailto:gdal-dev-bounces@lists.osgeo.org] <b>On Behalf Of </b>Chaitanya kumar CH<br><b>Sent:</b> Thursday, November 18, 2010 10:59 PM<br><b>To:</b> Boss, Ken (DNR)<br><b>Cc:</b> gdal-dev@lists.osgeo.org<br><b>Subject:</b> Re: [gdal-dev] KernelFilteredSource and nodata<o:p></o:p></span></p></div></div><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal style='margin-bottom:12.0pt'>Ken,<br><br>KernelFilteredSource is derived from SimpleSource. It does not support setting a nodata value to the source band.<br>Make sure your source image has a nodata value. You can set it using gdal_translate with the -a_nodata option.<o:p></o:p></p><div><p class=MsoNormal>On Fri, Nov 19, 2010 at 2:10 AM, Boss, Ken (DNR) <<a href="mailto:Ken.Boss@state.mn.us">Ken.Boss@state.mn.us</a>> wrote:<o:p></o:p></p><p class=MsoNormal>Chaitanya--<br><br>Thanks much for your response.<br><br>I don't think that I described my problem effectively. Let me restate it and see if we are talking about the same thing.<br><br>Here is a .vrt file that I think should accomplish what I am after:<o:p></o:p></p><div><p class=MsoNormal><br><VRTDataset rasterYSize='690' rasterXSize='575'><br> <SRS>EPSG:26915</SRS><br> <GeoTransform>190000, 1000, 0, 4795000, 0, 1000</GeoTransform><br> <VRTRasterBand band='1' dataType='Byte'><br> <KernelFilteredSource><o:p></o:p></p></div><p class=MsoNormal> <SourceFilename>D:/MyFiles/aptana_workspace/fire-webpage-mgmt/data/tmp_images/interpolated.tif</SourceFilename><o:p></o:p></p><div><p class=MsoNormal> <SourceBand>1</SourceBand><br> <NoDataValue>255</NoDataValue><o:p></o:p></p></div><div><p class=MsoNormal style='margin-bottom:12.0pt'> <Kernel normalized="1"><br> <Size>5</Size><br> <Coefs>0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04</Coefs><br> </Kernel><br> </KernelFilteredSource><br> </VRTRasterBand><br></VRTDataset><o:p></o:p></p></div><p class=MsoNormal>I believe this file to be saying that interpolated.tif uses a value of 255 for nodata, and that the kernel is to be normalized - that is, it should:<br><br> o ignore any nodata pixels within the kernel window when operating on a valid pixel<br> o adjust the coefficients to sum to one for the number of valid pixels in the window when operating on a valid pixel<br> o do nothing at all when centered on a nodata pixel<br><br>I would expect then that any pixels that go in as nodata should come out as nodata, and that pixels with valid data values that are near nodata pixels should not be influenced by them.<br><br>However, the outputs are not what I expect. The averaging is applied to both nodata and valid pixels alike. Have I set this up correctly, and is there a bug somewhere, or (more likely) is my vrt improperly constructed?<br><br>Thanks for your help,<br><span style='color:#888888'><br>--Ken</span><o:p></o:p></p><div><div><p class=MsoNormal><br>> -----Original Message-----<br>> From: Chaitanya kumar CH [mailto:<a href="mailto:chaitanya.ch@gmail.com">chaitanya.ch@gmail.com</a>]<br>> Sent: Wednesday, November 17, 2010 11:11 PM<br>> To: Boss, Ken (DNR)<br>> Cc: <a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>> Subject: Re: [gdal-dev] KernelFilteredSource and nodata<br>><br>> Ken,<br>><br>> The current kernel filter does not ignore the boundaries of nodata pixels. But it<br>> does normalize the kernel after ignoring the nodata pixels in the kernel, if the<br>> 'normalized' attribute is set to 1.<br>><br>> You can raise a ticket to request for this feature at<br>> <a href="http://trac.osgeo.org/gdal/newticket" target="_blank">http://trac.osgeo.org/gdal/newticket</a><br>><br>>> On Thu, Nov 18, 2010 at 6:08 AM, Boss, Ken (DNR) <<a href="mailto:Ken.Boss@state.mn.us">Ken.Boss@state.mn.us</a>><br>>> wrote:<br>>> Hello list--<br>>><br>>> I am attempting to filter a raster using gdal_translate (v 1.7) and a vrt with a<br>>> KernelFilteredSource. The input raster contains large areas of "nodata" values. I<br>>> would like the filter to ignore those areas. I have tried various combinations of<br>>> <NoDataValue>, <HideNoDataValue>, <NODATA> and <Kernel normalized='1'>,<br>>> but have not been able to prevent the filter from applying itself at data/nodata<br>>> boundaries.<br>>><br>>> My current VRT and gdal_translate command lines are below. Can anyone tell me<br>>> what I am doing wrong?<br>>><br>>> Thanks,<br>>><br>>> Ken Boss<br>>> Minnesota DNR<br>>><br>>> ===== kernel_filter.vrt ===================<br>>><br>>> <VRTDataset rasterYSize='690' rasterXSize='575'><br>>> <SRS>EPSG:26915</SRS><br>>> <GeoTransform>190000, 1000, 0, 4795000, 0, 1000</GeoTransform><br>>> <VRTRasterBand band='1' dataType='Byte'><br>>> <KernelFilteredSource><br>>> <SourceFilename>interpolated.tif</SourceFilename><br>>> <SourceBand>1</SourceBand><br>>> <NoDataValue>255</NoDataValue><br>>> <HideNoDataValue>1</HideNoDataValue><br>>> <Kernel normalized='1'><br>>> <Size>5</Size><br>>> <Coefs>0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04<br>>> 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04 0.04</Coefs><br>>> </Kernel><br>>> </KernelFilteredSource><br>>> </VRTRasterBand><br>>> </VRTDataset><br>>><br>>> =====================================<br>>><br>>> gdal_translate -of GTiff -ot Byte -a_srs EPSG:26915 kernel_filter.vrt<br>>> kernel_filtered.tif<br>>><br>> _______________________________________________<br>> gdal-dev mailing list<br>> <a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>> <a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>><br>><br>><br>> --<br>> Best regards,<br>> Chaitanya kumar CH.<br>> /t<span style='font-family:"MS Mincho"'>ʃ</span>a<span style='font-family:"MS Mincho"'>ɪ</span>θənjə/ /k<span style='font-family:"MS Mincho"'>ʊ</span>m<span style='font-family:"MS Mincho"'>ɑ</span>r/<br>> +91-9494447584<br>> 17.2416N 80.1426E<o:p></o:p></p></div></div></div><p class=MsoNormal><br><br clear=all><br>-- <br>Best regards,<br>Chaitanya kumar CH.<br>/t<span style='font-family:"MS Mincho"'>ʃ</span>a<span style='font-family:"MS Mincho"'>ɪ</span>θənjə/ /k<span style='font-family:"MS Mincho"'>ʊ</span>m<span style='font-family:"MS Mincho"'>ɑ</span>r/ <br>+91-9494447584<br>17.2416N 80.1426E<o:p></o:p></p></div></div></body></html>