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.<br>
<br><div class="gmail_quote">On Fri, Nov 19, 2010 at 2:10 AM, Boss, Ken (DNR) <span dir="ltr"><<a href="mailto:Ken.Boss@state.mn.us">Ken.Boss@state.mn.us</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
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:<br>
<div class="im"><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>
</div> <SourceFilename>D:/MyFiles/aptana_workspace/fire-webpage-mgmt/data/tmp_images/interpolated.tif</SourceFilename><br>
<div class="im"> <SourceBand>1</SourceBand><br>
<NoDataValue>255</NoDataValue><br>
</div><div class="im"> <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><br>
<br>
</div>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>
<font color="#888888"><br>
--Ken<br>
</font><div><div></div><div class="h5"><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ʃaɪθənjə/ /kʊmɑr/<br>
> +91-9494447584<br>
> 17.2416N 80.1426E<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Best regards,<br>Chaitanya kumar CH.<br>/tʃaɪθənjə/ /kʊmɑr/ <br>+91-9494447584<br>17.2416N 80.1426E<br>