<div dir="ltr"><div dir="ltr"><div>Hi Karsten, </div><div><br></div><div>An option is to run gdallocationinfo, get the col (P) and row (L) value and then feed that into </div><div><br></div><div>gdal_translate -srcwin P L col -outsize 1 1  in.tif out.tif ##I think that's right</div><div><br></div><div><div>You otherwise have to do the "grid logic", which (while easy ...) I'd rather like to see exposed explicitly in GDAL as functionality and I'm exploring this in various ways atm.  (I think you can't use -projwin as it seems to need at least a half pixel width and height to count as non zero so it's still a bit tricky to get the right values). <br></div><div><br></div></div><div>Cheers, Mike</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jan 27, 2023 at 5:01 PM karsten <<a href="mailto:karsten@terragis.net">karsten@terragis.net</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><u></u>



<div>
<div dir="ltr" align="left"><font size="2" face="Arial"><span><span><font color="#000080">Hi 
All,</font></span></span></font></div>
<div>
<div><span><span><font color="#000080"></font></span></span> </div>
<div><span><span><font color="#000080">in an application I have a user click on a map and I need to 
programmatically clip a raster layer </font></span></span></div>
<div><span><span><font color="#000080">to the one raster cell for those coordinates at the click 
location. </font></span></span></div>
<div><font color="#000080"><span><span>I was looking e.g. at <a href="https://gdal.org/programs/gdal_rasterize.html" target="_blank">https://gdal.org/programs/gdal_rasterize.html</a> to 
run as a command <br>line tool, but so far could not figure out if I can some 
can "feed" that a vector point to the tool as an 
input</span></span></font></div>
<div><font color="#000080"><span><span>without first writing a shape file with the point 
onto disk. </span></span></font></div>
<div><font color="#000080"><span><span>
<div><span><span><font color="#000080"><br>Does somebody have any pointers to a command line tool (or a 
python example) I could use to achieve this 
?</font></span></span><span><span></span></span><font color="#000080"><span><span> </span></span></font></div></span></span></font></div>
<div><font color="#000080"><span><span></span></span></font> </div>
<div><font color="#000080"><span>Thanks<span> <br></span></span><span style="font-size:11pt;font-family:Verdana,"sans-serif"">Karsten</span></font></div>
<div><span style="font-size:11pt;font-family:Verdana,"sans-serif""></span> </div></div>
<div align="left">
<div><span style="font-size:11pt;font-family:Verdana,"sans-serif"">Karsten 
Vennemann<br></span><span style="font-size:9pt;font-family:Verdana,"sans-serif"">Principal<br><br><span style="color:rgb(0,50,100)">Terra GIS LTD<br>2119 Boyer Ave E 
<br>Seattle, WA  98112<br>USA <br><a title="http://www.terragis.net/
blocked::http://www.terragis.net/" href="http://www.terragis.net/" target="_blank"><font title="http://www.terragis.net/" color="#0000ff">www.terragis.net</font></a></span></span></div></div></div>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature">Michael Sumner<br>Software and Database Engineer<br>Australian Antarctic Division<br>Hobart, Australia<br>e-mail: <a href="mailto:mdsumner@gmail.com" target="_blank">mdsumner@gmail.com</a></div></div>