Me alegro, saludos

Hola Jesús, entiendo que no hay problema en sumar los ráster de varios días ya que cada uno de ellos tiene su número de horas de insolación. Eso sí, si son dias consecutivos, serán raster iguales o muy similares.

Hola. Si que veo que es una solución viable, pero realmente lo veo un proceso pesado... Generar, por ejemplo, 200 rásters y luego hacer un sumatorio, frente a un cálculo directo de "horas de sol entre 2 fechas concretas". Supongo que no habrá más remedio que sumar uno a uno. Gracias por la respuesta!

@@jesushernandezsaseta7563 Bueno, tienes la opción de hacer un raster cada cinco días, sumar y multiplicar por 5, no creo que el error sea significativo. Por otra parte, me suena que hay extensiones que permiten poner rangos largos, pero habría que buscarlas.

Con el detalle de 0,25 horas la verdad es que no se apreciaría. Muchas gracias!!

Saludos igualmente. Vas a sketchfab.com/ y en el botón "Sign up"

Gracias a ti, saludos.

ArcGIS online es gratuito. Busca por internet "cargar mapas base en QGIS" y encontrarás rápidamente las instrucciones.

Me alegro de que te guste, saludos.

Hola, por convención, las orientaciones geográfica son siempre las mismas: 0º al norte, 90º al este, etc. No debería ser necesario cambiar nada. En cualquier caso, sugiero que hagas un ensayo y verifiques que los valores de insolación en cualquier ladera orientada al norte son mayores que en otra orientada al sur.

Hola, no hay definida una "unidad" de rugosidad. Lo que sí tenemos es que los valores tienen límites acotados entre 0 y 1; cero cuando la rugosidad es "infinita" (la suma de vectores unitarios es nula; no va a pasar nunca en un MDE) y 1 cuando es una superficie plana (vectores paralelos; puede pasar en un MDE).

Para tener datos reales, lo mejor es acudir a la estación meteorológica oficial más cercana porque miden ambas cosas y algunas más. Este tipo de simulaciones genera datos potenciales, es decir, en ausencia de nubosidad.

Dice descargar directamente esa información de la web de la estación meteorológica e importarla a Qgis?

Pues no debe ser, claro. Verifica que estás usando un MDE proyectado en el que las unidades son metros tanto en los ejes X e Y como en la altura ya que si está en grados (EPSG 4326) no funciona.

Está en UTM, 25830

@@karabaia5417 Pues lo siento pero no tengo idea de qué puede pasar.

Estos métodos son para escalas no muy detalladas ya que normalmente no hay datos microclimáticos. Sería posible estimar la radiación solar en un MDE de pocos metros y tal vez esa sea la única vía para trabajar con ese detalle. El método es simple aunque conviene estudiar la teoría (como siempre). Explico esa estimación en kzbin.info/www/bejne/a3LIgXSld5V-gqssi=I-IVGZ7rf1dnkJhi

Gracias.

Gracias.

La única forma que conozco es exportar el modelo a formato gltf y subirlo a Sketchfab.como en skfb.ly/6TDZX

No, yo tengo cada proyecto en una carpeta sin más control que evitar caracteres extraños en el "path": tildes, espacios... Usa solamente letras sin tilde, números y el carácter _

En el segundo 16 tienes las referencias que yo he usado como principios de la elaboración del video. Las explicaciones concretas y su implementación en QGIS son cosa mía, lógicamente.

Buenos días igualmente. Sólo hace falta un modelo digital de elevaciones, no imágenes. El MDE puede estar en formato geotiff o en cualquier otro, no es limitante. Lógicamente, vale cualquier MDE de cualquier zona del mundo.

Debería estar, busca por EPSG:3857 Si no lo encuentras, puedes usar cualquier otra proyección siempre que las unidades sean metros.

@ lo logre pero el resultado no me convense: los niveles del terrenos me quedan como muy achatados no parece lo real. Yo estuve en el lugar donde quiero sacar el 3d y no para mi los niveles son mucho mas bruscos

Supongo que tienes las unidades en metros incluidas las de elevación ¿no? Ahora no puedo decir más porque estoy de vacaciones y este ordenador no es el de trabajo @@MrNiqozer

Sí, claro, si es una ortoimagen bien hecha no hay problema.

Hola, en mi correo [email protected]

mappinggis.com/2018/03/como-anadir-mapas-base-en-qgis-3-0-openstreetmap-google-carto-stamen/

Sí, aunque con alguna complicación. La ortoimagen, si está georreferenciada correctamente, se carga en QGIS y el complemento la usará sin problema sobr el MDE. La exageración del relieve es una opción del complemento; de hecho, en estos modelos del terreno, la Z se suele exagerar a 1.5x o 2x. Otra cuestión diferente es la de integrar modelos 3D de Agisoft. QGIS no reconoce esos formatos (obj, stl, dae...). En ese caso tendrás que exportar el modelo a un formato intermedio (obj, por ejemplo) y transformarlo a un fichero de texto como puntos xyz (puede hacer con Meshlab); luego habrá que leerlo en QGIS como texto delimitado y trasformarlo a raster mediante interpolación.

@ que detalle de tu parte el poder responder tan extensamente, me aclaraste completamente el panorama, saludos

No lo veo posible, suponiendo que te refieras a la superficie de la lámina de agua. La rugosidad de los océanos sí se mide desde satélites por la dispersión de la señal enviada por radares pero la resolución temporal es baja (10 km aprox.) Si te refieres al fondo de cauce habría que tener la batimetría pero sí sería posible por este método si se consigue la resolución suficiente.

Gracias Alexander, muy amable. Me alegro de que aquel librito valiera para algo, ya llovió mucho desde entonces :-) Un abrazo.

Mi correo es [email protected] pero aviso de que no trabajo específicamente con datos LiDAR sino con MDE ya elaborados.

@ yo cuento con dron para obtener datos LIDAR, significa que no se podría empezar desde ahí, si no desde un archivo LAS?

@@LuisRubin Tampoco. ya que los archivos de puntos no son propiamente MDE,. Yo no trabajo con esos formatos ya que me dedico a modelizar a partir de los MDE raster ya elaborados. No creo que pueda ayudar con eso, lo siento.

No, necesariamente. El MDE es la estructura de datos pero esa estructura (raster habitualmente) puede venir de fuentes diversas: LiDAR, radar, fotogrametría...

Hola, gracias por la opinión. En principio, no sería engorroso, es cuestión de probar a ver en qué casos es interesante aplicarlo.

Hola, habría que demostrar que esa covariable tiene una correlación estrecha con la que quieres interpolar. La interpolación en sí de una variable no necesita covariables, pero en el caso climático es importante porque los datos están dispersos y no reflejan bien la variablidad espacial (apenas hay datos en zonas de montaña). Dado que la elevación correlaciona muy significativamente con la temperatura, es razonable aprovecharlo para mejorar los mapas. No tengo claro que el valor de la vivienda sea interpolable, habría que realizar análisis previos, pero una covariable podría ser la proximidad a las vías de comunicación, por ejemplo, ya que es posible que el valor sea mayor cuando estás mejor comunicado. En cualquier caso, insisto, siempre habría que analizarlo previamente.

@ entiendo. Lo otro que tr quería preguntar es sobre los parámetros del variograma. En el video tu dejas los que vienen por default pero la teoría habla de que primero hay que sacar el variograma para determinar el nugget, la meseta y demás.... En ArcGIS se puede sacar antes de proceder a hacer el kriging. Hay alguna manera de hacer esto en qgis?

@@EdgardsGranados En QGIS creo que no, aunque me suena alguna rutina de R que podía integrarse, creo. En este caso concreto de interpolación de temperaturas, la incertidumbre es tan alta y los datos tan imprecisos que no creo que merezca la pena hacer ningun análisis pero, si se quiere rigor, sí sería necesario.

@ muchas gracias por tus respuestas Angel.

Lamentablemente, no tengo idea de qué puede pasar, a mí me funciona perfectamente y lo uso bastante.

@ me di cuenta que cuando tengo una imágen satelital de fondo se queda trabado y se cierra. Le quito la imagen y funciona más rápido. Pero estoy modelando un sistema de monocultivo de mango y en los puntos de siembra le coloco la imagen. Y allí se cae por completo.

@@Hams2522 Es raro porque yo le pongo encima imágenes de varios tipos y formatos y no hay problema.

@ saludos. Si me parece super raro. En mi laptop no funciona. Y en otra computadora de mesa se cierra solo.

Hola, los filtros sirven para definir el entorno para el cual quieres calcular la rugosidad. La verdad es que nunca he probado no usarlos porque es lo conceptualmente correcto. No se si el problema que tienes está ahí. Lo que se me ocurre es que vayas mirando con cuidado los valores mínimos y máximos que te van saliendo en las capas calculadas para estar segura de que corresponden con lo que se calcula, hasta encontrar el error. Por ejemplo, es muy típico calcular los senos y cosenos con la capa en grados en vez de en radianes y eso hace que no salga nada. Igualmente, en la operación de raiz cuadrada de la suma de los cuadrados de las capas, asegúrate de poner todos los paréntesis necesarios para que se haga bien. Suerte y saludos.

@ Podria facilitarnos como realizar el filtro para rugosidad?, tengo el mismo el error...

@@10105864 No tengo una rutina específica. La tenía pero en AML, el antiguo lenguaje de Arcinfo. Por si alguien se anima a traducirla, es este: /* Cálculo del modelo de rugosidad /* Método: /* 1. Cálculo de los vectores unitarios normales a la superficie en un /* entorno de nxn celdas /* 2. Cálculo del módulo del vector suma normalizado en el intervalo 0-1 /* /* Autor: Angel M. Felicísimo /* Esta rutina puede usarse libremente; si se utiliza en algún trabajo /* sería elegante citar su procedencia &args theMDE .ventana &if %:program% <> GRID &then &return &warning Este programa debe ejecutarse en GRID &if [null %theMDE%] &then &return Uso: &r rugosidad MDE VENTANA &if not [exists %theMDE% -grid] &then &return El MDE no existe: &r rugosidad MDE VENTANA &if [null %.ventana%] &then &do &return Usando una ventana de 5x5 ... &s .ventana = 5 &end &type Calculando pendiente ... theMDP = slope (%theMDE%) / 57.2956 &type Calculando orientación ... theMDO = aspect (%theMDE%) / 57.2956 &type Calculando componentes X,Y,Z ... theXT = sin(theMDP) * cos(theMDO) theYT = sin(theMDP) * sin(theMDO) theZT = cos(theMDP) kill theMDO all &type Revisando horizontales ... XT = con(theMDP == 0, 0, theXT) YT = con(theMDP == 0, 0, theYT) ZT = con(theMDP == 0, 1, theZT) kill theXT all; kill theYT all; kill theZT all kill theMDP all &Type Calculando los vectores medios ... sumXT = focalmean(XT,rectangle,%.ventana%,%.ventana%) sumYT = focalmean(YT,rectangle,%.ventana%,%.ventana%) sumZT = focalmean(ZT,rectangle,%.ventana%,%.ventana%) kill XT all; kill YT all; kill ZT all &type Generando el modelo de rugosidad ... mdr = (sqrt(sqr(sumXT) + sqr(sumYT) + sqr(sumZT))) * 100 kill sumXT all; kill sumYT all; kill sumZT all &type Cálculo finalizado con ventana de %.ventana% x %.ventana% &type El modelo de rugosidad es MDR &return

Gracias Ramón, un saludo

Creo que el problema es que estás trabajando en grados EPSG 4326 probablemente. Usa psudomercator o cualquier otra proyección en el proyecto y te saldrá bien (creo).

me pasa lo mismo, pudiste solucionarlo? ^^

Gracias a tí por comentar, saludos.

Hola Antony, la respuesta corta es que no. Hace años, cuando trabajada con Arcinfo, lo tenía programado en AML porque, lógicamente, es muy útil tener una rutina que haga todas las operaciones. Actualmente, Arcinfo ya no funciona y la rutina tampoco :-( A ver si alguien con más conocimiento de Python se anima, yo sólo programo en lenguajes añejos :-)

Gracias. aun así, nos ayuda mucho esta detallada explicación, espero que no deje de subir información QGIS, que nos ayude a los que queremos aprender. ¡¡Saludos!!

No te preocupes, soy el primero que conoce las limitaciones de los modelos predictivos, especialmente cuando no se tienen en cuenta más que variables climáticas obviando las biológicas, edáficas... La idea es mostrar cómo se pueden hacer esos modelos.

​@ Si, era un comentario conociendo solo el título. Tengo mi propia base de rasters ya hecha con las siguientes capas: PH, PMA, TMA, Altitud, Pisos bioclimáticos, hillshide e índice topográfico de humedad. Pensaba que era un video para planificar un caso local de restauración. Un abrazo.

¿A qué viene eso Manuel?