<div dir="ltr"><div><div><div><div><div><div><div><div><div>Estimado Ulises:<br></div>A mí se me ocurre lo siguiente:<br><br></div>1) Para cada punto con elevación verificada en campo "muestrea" el valor actual de MDT (añade una columna a la capa de puntos con la cota actual del MDT). En QGis, a través de SAGA, tienes "Add grid values to points" por ejemplo.<br><br></div>2) Con una simple operación matemática entre columnas calcula la corrección que debería realizarse en cada punto para que el MDT "malo" tuviera los valores "correctos" (valor verificado - valor actual = corrección). Puedes utilizar, por ejemplo, en QGis la calculadora de campos.<br><br></div>3) A partir de la capa de puntos, y tomando el valor de la corrección, crea una nueva capa raster (capa que tendrá muy pocos valores "válidos" y muchos "nodata"). En QGis tienes "Rasterize (vector to raster)".<br><br>4) Rellena los "nodata" con algún algoritmo de interpolación (esta nueva capa raster servirá de corrección para todos los puntos del MDT "malo"). En QGis tienes un algoritmo que se llama "Fill nodata".<br><br></div>5) Utiliza la calculadora raster para sumar el MDT "malo" y la capa rellena de corrección. El resultado debería ser el MDT "bueno".<br><br></div>6) Repite el punto 1 para verificar que el valor muestreado en todos los puntos con el MDT bueno coincide exactamente con el valor verificado en campo (lo de exactamente es un decir, dependiendo del método de interpolación que se use en el punto 4, el valor en la celda podría no ser igual al valor en el punto -esto no lo he comprobado, pero mi intuición me dice que podría ser así teniendo en cuenta cómo funcionan esos algoritmos-).<br><br></div>Seguro que a alguien se le ocurre algún otro procedimiento, por ejemplo con kriging, co-kriging y cosas similares, pero teniendo en cuenta la cantidad de parámetros que necesitarías estimar en esos métodos y lo difícil que es ajustarlos bien... me quedo con la cuenta de la vieja (o en este caso del viejo, que soy yo).<br><br></div>Espero haberte ayudado,<br></div>Juan Carlos<br></div><div id="DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2"><br> <table style="border-top:1px solid #d3d4de">
<tr>
<td style="width:55px;padding-top:18px"><a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" target="_blank"><img src="https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-animated-no-repeat-v1.gif" alt="" width="46" height="29" style="width: 46px; height: 29px;"></a></td>
<td style="width:470px;padding-top:17px;color:#41424e;font-size:13px;font-family:Arial,Helvetica,sans-serif;line-height:18px">Libre de virus. <a href="https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail" target="_blank" style="color:#4453ea">www.avast.com</a> </td>
</tr>
</table>
<a href="#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2" width="1" height="1"></a></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div style="text-align:left"><span><p style="margin:0px;text-align:center"><span style="font-size:8.0pt;line-height:107%"><br></span></p></span></div><div style="text-align:left">_______________________________________________________<br></div>Juan Carlos Giménez Fernández<br><b>Grado en Ingeniería Forestal y del Medio Natural</b><br><i>Centro Universitario de Plasencia</i><br>Universidad de Extremadura<br>Avda. Virgen del Puerto nº 2<br>10600 Plasencia Cáceres (Spain)<br>Phone: <span style="color:rgb(0,0,255)">+34-927-25-70-00 (ext. 52148)</span><br>FAX: <span style="color:rgb(0,0,255)">+34-927-42-52-09</span><br>Móvil: <span style="color:rgb(0,0,255)">+34-<span>655-06-87-86</span></span><span><p style="margin:0px"><span style="font-size:8.0pt;line-height:107%"><img src="cid:2fed91a4d97429243d645f5e5dee671824e92281@zimbra"></span></p></span>_______________________________________________________<br><p style="margin:0px"><span style="font-size:8.0pt;line-height:107%"></span></p></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">El 23 de marzo de 2018, 1:54, Ulises Ibarra <span dir="ltr"><<a href="mailto:ulisesmartinibarra@gmail.com" target="_blank">ulisesmartinibarra@gmail.com</a>></span> escribió:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">
<div class="m_-5831616708324047702m_-1260165718438023363gmail-_5pbx m_-5831616708324047702m_-1260165718438023363gmail-userContent m_-5831616708324047702m_-1260165718438023363gmail-_22jv m_-5831616708324047702m_-1260165718438023363gmail-_3576" id="m_-5831616708324047702m_-1260165718438023363gmail-js_l"><p>Estimada Lista:</p><p>Una duda que me quita el sueño: Tengo un Modelo Digital de Elevaciones que se ubtuvo por medio de filtrar una nube de puntos resultado de un vuelo con drone. Aún teniendo puntos de control terrestre hay partes que quedaron altas.</p><p>Requiero ajustar en algunas franjas la elevación. Ttengo un muestreo de puntos con elevaciones verificadas al costado de algunas terracerías, que vendrían siendo Terreno Natural.</p><p>Pregunta: <br>
¿Qué se les ocurre y con qué metodo para modificar el MDT por
interpolación para que por zonas baje o suba según los puntos de
muestreo?</p><p> Tengo dos capas:<br> 1.-MDT con elevación no verificada.<br> 2.-Muestreo de puntos con elevación verificada en campo.</p><p>Cordialmente</p><span class="HOEnZb"><font color="#888888"><p>Ulises<br></p></font></span></div>
<br></div>
<br>______________________________<wbr>_________________<br>
Spanish mailing list<br>
<a href="https://lists.osgeo.org/mailman/listinfo/spanish" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/spanish</a><br>
<a href="http://es.osgeo.org" rel="noreferrer" target="_blank">http://es.osgeo.org</a><br>
<a href="http://twitter.com/osgeoes" rel="noreferrer" target="_blank">http://twitter.com/osgeoes</a><br></blockquote></div><br></div>