tag:blogger.com,1999:blog-15095239960504663252024-03-13T23:25:20.494-03:00Mapas GeológicosBlog que explica cómo hacer figuras (mapas, gráficos, perfiles, bloques 3D) con GMT (Generic Mapping Tools) con Windows.Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.comBlogger82125tag:blogger.com,1999:blog-1509523996050466325.post-73756482112591725502022-09-10T18:46:00.007-03:002022-09-16T14:42:12.387-03:00Webinar Gratuito «Elaboración de Mapas con GENERIC MAPPING TOOLS»<p></p><div class="separator" style="clear: both; text-align: center;"><a href="https://blogger.googleusercontent.com/img/a/AVvXsEinPJmu2XfG2oWSL8PsBm9c85POWa2dUiQxqQVXLcpngZ9B1MOCpEYCMuQZ4lEhOfkjR3VKTetAi5rnyyiT3d-rWAEVDpAIziYa4mASHVrNp7mFOPF_PFy31ybZBZrIIF7W8gj1CmYc9Ka7QZSypk12cSuG_ILtbtA_OYLh10RqbJxniVodFWTohf2YNQ" style="margin-left: 1em; margin-right: 1em;"><img alt="" data-original-height="945" data-original-width="960" height="630" src="https://blogger.googleusercontent.com/img/a/AVvXsEinPJmu2XfG2oWSL8PsBm9c85POWa2dUiQxqQVXLcpngZ9B1MOCpEYCMuQZ4lEhOfkjR3VKTetAi5rnyyiT3d-rWAEVDpAIziYa4mASHVrNp7mFOPF_PFy31ybZBZrIIF7W8gj1CmYc9Ka7QZSypk12cSuG_ILtbtA_OYLh10RqbJxniVodFWTohf2YNQ=w640-h630" width="640" /></a></div><div class="separator" style="clear: both; text-align: center;"><a href="https://blogger.googleusercontent.com/img/a/AVvXsEhUrXPBzIlRQj51FzN943uaqL9wI3QoPVWcsHXWPejSb-QFTupwVViJXMWE7N76704Ixz9J8rVxfJ2IuHjbWpJIb-l6kv8TnpnictHQElB3p8ZPKl6Y1o1uE8mzMMOix_bY3yAATs7QrCBMaWe0Xselzaf1pBPa4m9y7R688IHDeR20AMEI6oTLXy3qJw" style="margin-left: 1em; margin-right: 1em;"><div class="separator" style="clear: both; text-align: center;"><br /></div></a></div><p></p><h3 style="line-height: 1.38; margin-bottom: 6pt; margin-top: 20pt; text-align: left;"><span style="font-family: Arial; font-size: medium; font-weight: 400; white-space: pre-wrap;">Repositorio con material y link al video:</span></h3><div><a href="https://github.com/Esteban82/CAC-2022/blob/main/README.md " style="font-family: Arial; white-space: pre-wrap;"><span style="font-size: medium;">https://github.com/Esteban82/CAC-2022/blob/main/README.md </span></a></div><div style="line-height: 1.38; margin-bottom: 6pt; margin-top: 20pt; text-align: left;"><span><a name='more'></a></span></div><h1 dir="ltr" style="line-height: 1.38; margin-bottom: 6pt; margin-top: 20pt;"><span style="font-family: Arial; font-size: 21pt; font-weight: 400; white-space: pre-wrap;"><strike>Inscripción:</strike></span></h1><div><span style="font-family: Arial; font-size: 21pt; font-weight: 400; white-space: pre-wrap;"><strike><br /></strike></span></div><p style="background-color: white; border-color: rgb(232, 232, 232); box-sizing: border-box; color: #666666; font-family: Roboto; font-size: 15px; margin: 0px 0px 20px;"><strike><img alt="🔸" height="16" loading="lazy" src="https://static.xx.fbcdn.net/images/emoji.php/v9/t72/1/16/1f538.png" style="border: 0px rgb(232, 232, 232); box-sizing: border-box; height: auto; max-width: 100%; vertical-align: middle;" width="16" /> Link de inscripción<a href="https://www.eventbrite.com.ar/e/elaborando-mapas-con-generic-mapping-tools-tickets-415324917007" style="background: transparent; border-color: rgb(232, 232, 232); box-sizing: border-box; color: #537383;">: https://www.eventbrite.com.ar/e/elaborando-mapas-con-generic-mapping-tools-tickets-415324917007</a></strike></p><div class="separator" style="clear: both; text-align: center;"><p style="background-color: white; border-color: rgb(232, 232, 232); box-sizing: border-box; color: #666666; font-family: Roboto; font-size: 15px; margin: 0px 0px 20px; text-align: start;"><strike><img alt="🔸" height="16" loading="lazy" src="https://static.xx.fbcdn.net/images/emoji.php/v9/t72/1/16/1f538.png" style="background-color: transparent; border: 0px rgb(232, 232, 232); box-sizing: border-box; height: auto; max-width: 100%; vertical-align: middle;" width="16" /> Cierre de inscripción: miiércoles 14 de septiembre a las 13 horas.</strike></p><div><br /></div></div><h1 dir="ltr" style="line-height: 1.38; margin-bottom: 6pt; margin-top: 20pt;"></h1><h1 dir="ltr" style="line-height: 1.38; margin-bottom: 6pt; margin-top: 20pt;"><span style="background-color: transparent; color: black; font-family: Arial; font-size: 18pt; font-style: normal; font-variant: normal; font-weight: 400; text-decoration: none; vertical-align: baseline; white-space: pre;">Material para ver antes del webinar: </span></h1><div class="separator" style="clear: both; text-align: center;"><span id="docs-internal-guid-f4c8257c-7fff-7677-7dd6-ed21fbeb9bed"><p dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt;"><span style="font-family: Arial; font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline; white-space: pre-wrap;">Se invita a ver a los siguientes links (previamente idealmente para ver mapas y figuras más complejas): </span></p><ul style="margin-bottom: 0px; margin-top: 0px; padding-inline-start: 48px;"><li aria-level="1" dir="ltr" style="font-family: Arial; font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; list-style-type: disc; vertical-align: baseline; white-space: pre;"><p dir="ltr" role="presentation" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt; text-align: left;"><a href="https://www.instagram.com/mapasgeologicos/" style="text-decoration-line: none;"><span style="color: #1155cc; font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; text-decoration-line: underline; text-decoration-skip-ink: none; vertical-align: baseline; white-space: pre-wrap;">https://www.instagram.com/mapasgeologicos/</span></a><span style="font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline; white-space: pre-wrap;"> : Instragram oficial con ejemplos de mapas y animaciones realizados con GMT.</span></p></li><li aria-level="1" dir="ltr" style="font-family: Arial; font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; list-style-type: disc; vertical-align: baseline; white-space: pre;"><p dir="ltr" role="presentation" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt; text-align: left;"><a href="https://www.generic-mapping-tools.org/" style="font-size: 11pt; text-decoration-line: none;"><span style="color: #1155cc; font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; text-decoration-line: underline; text-decoration-skip-ink: none; vertical-align: baseline; white-space: pre-wrap;">https://www.generic-mapping-tools.org/</span></a><span style="font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline; white-space: pre-wrap;">: Sitio oficial de GMT.</span></p></li></ul><h1 dir="ltr" style="line-height: 1.38; margin-bottom: 6pt; margin-top: 20pt; text-align: left;"><span style="font-family: Arial; font-size: 19pt; font-variant-east-asian: normal; font-variant-numeric: normal; font-weight: 400; vertical-align: baseline; white-space: pre-wrap;">Programas a utilizar:</span></h1><ul style="margin-bottom: 0px; margin-top: 0px; padding-inline-start: 48px;"><li aria-level="1" dir="ltr" style="font-family: Arial; font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; list-style-type: disc; vertical-align: baseline; white-space: pre;"><p dir="ltr" role="presentation" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt; text-align: left;"><span style="font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline; white-space: pre-wrap;">Git for Windows: </span><a href="https://gitforwindows.org/" style="text-decoration-line: none;"><span style="color: #1155cc; font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; text-decoration-line: underline; text-decoration-skip-ink: none; vertical-align: baseline; white-space: pre-wrap;">https://gitforwindows.org/</span></a></p></li><li aria-level="1" dir="ltr" style="font-family: Arial; font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; list-style-type: disc; vertical-align: baseline; white-space: pre;"><p dir="ltr" role="presentation" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt; text-align: left;"><span style="font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline; white-space: pre-wrap;">GMT 6.4: </span><a href="https://www.generic-mapping-tools.org/download/" style="text-decoration-line: none;"><span style="color: #1155cc; font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; text-decoration-line: underline; text-decoration-skip-ink: none; vertical-align: baseline; white-space: pre-wrap;">https://www.generic-mapping-tools.org/download/</span></a></p></li></ul><div style="text-align: left;"><br /></div><p dir="ltr" style="line-height: 1.38; margin-bottom: 0pt; margin-top: 0pt; text-align: left;"><span style="font-family: Arial; font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline; white-space: pre-wrap;">Ver link para su instalación: </span><a href="https://youtu.be/1bPMIN7noTI" style="text-decoration-line: none;"><span style="color: #1155cc; font-family: Arial; font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; text-decoration-line: underline; text-decoration-skip-ink: none; vertical-align: baseline; white-space: pre-wrap;">https://youtu.be/1bPMIN7noTI</span></a><span style="font-family: Arial; font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline; white-space: pre-wrap;">. </span></p><span style="font-family: Arial; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline; white-space: pre-wrap;"><div style="font-size: 11pt; text-align: left;"><span style="font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline;">Nota: La instalación del </span><span style="font-size: 11pt; font-style: italic; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline;">software </span><span style="font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline;">Visual Studio Code no es necesaria.</span></div><div style="font-size: 11pt; text-align: left;"><span style="font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline;"><br /></span></div><div style="font-size: 11pt; text-align: left;"><span style="font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline;"><br /></span></div><div style="text-align: left;"><span style="font-size: medium; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline;">Referencias:</span></div><div style="font-size: 11pt; text-align: left;"><span style="font-size: 11pt; font-variant-east-asian: normal; font-variant-numeric: normal; vertical-align: baseline;"><a href="https://centroargentinodecartografia.org/2022/09/webinar-gratuito-elaboracion-de-mapas-con-generic-mapping-tools/" style="font-family: "Times New Roman"; font-size: medium; text-align: center; white-space: normal;">https://centroargentinodecartografia.org/2022/09/webinar-gratuito-elaboracion-de-mapas-con-generic-mapping-tools/</a></span></div></span></span></div>Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0tag:blogger.com,1999:blog-1509523996050466325.post-15637909364109603652022-09-01T09:04:00.006-03:002022-09-01T09:04:58.188-03:00Curso Postgrado - Buenos Aires. Octubre-Noviembre 2022 - Presencial.<div style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; text-align: center;"><span style="color: red;"><span style="font-size: x-large;"><b>Geoinformática aplicada a la cartografía multitemática</b></span></span></div><p align="center" class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px; text-align: center;"><span style="color: navy; font-family: Verdana, "sans-serif"; font-size: 24pt;"></span></p><div style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small;"><p align="center" class="MsoNormal" style="margin: 0px; text-align: center;"><i><b><span style="color: navy; font-family: Verdana, "sans-serif"; font-size: 18pt;">Curso de Postgrado y Extensión</span></b></i><i><b><span style="color: navy; font-family: Verdana, "sans-serif"; font-size: 13.5pt;"></span></b></i><span style="color: navy; font-family: Verdana, "sans-serif"; font-size: 10pt;"></span></p><p align="center" class="MsoNormal" style="margin: 0px; text-align: center;"><span style="color: navy; font-family: Verdana, "sans-serif"; font-size: 10pt;"></span></p><p align="center" class="MsoNormal" style="margin: 0px; text-align: center;"><span style="color: navy;">Por:</span><b><span style="color: navy; font-size: 13.5pt;"> Dr. Federico D. ESTEBAN</span></b></p><p align="center" class="MsoNormal" style="margin: 0px; text-align: center;"><span style="font-size: medium;">Departamento de Ciencias Geológicas - FCEN - Universidad de Buenos Aires<span style="color: navy; font-family: Arial, "sans-serif";"></span></span></p><p align="center" class="MsoNormal" style="margin: 0px; text-align: center;"><span style="color: navy; font-family: Arial, "sans-serif"; font-size: 10pt;"></span></p></div><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><b><br /></b></p><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><span style="color: black;"><b><br /></b></span></p><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><b style="color: black; font-family: Arial, "sans-serif";">Fecha de Inicio: </b><span style="color: black; font-family: Arial, "sans-serif";">Martes 11 de octubre, 14 horas. </span><br /></p><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><b style="color: black; font-family: Arial, "sans-serif";">Fechas de Cursada:</b><span style="color: black; font-family: Arial, "sans-serif";"> 11 de octubre al 17 de noviembre (6 semanas).</span><span style="color: black;"><span style="font-family: Arial, "sans-serif"; font-size: 10pt;"><br /></span></span></p><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><b style="color: black;"><span style="font-family: Arial, "sans-serif";">Días y horario de clases:</span></b><span style="color: black; font-family: Arial, "sans-serif";"> Martes y Jueves de 14 a 18 hs.</span><br /></p><span style="background-color: white; font-family: Arial, Helvetica, sans-serif; font-size: small;"></span><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><span style="color: black;"><b><span style="font-family: Arial, "sans-serif";">Duración del curso:</span></b><span style="font-family: Arial, "sans-serif";"> 12 clases (Teórico - Práctic</span><span style="font-family: Arial, "sans-serif";">o</span><span style="font-family: Arial, "sans-serif";">), </span><span style="font-family: Arial, "sans-serif"; font-size: 10pt;"><br /></span></span></p><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><span style="color: black;"><b>Carga horaria total</b>: 48 hs.</span></p><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><span style="color: black;"><span style="font-family: Arial, "sans-serif"; font-size: 10pt;"></span><span style="font-family: Arial, "sans-serif"; font-size: 14pt;"></span></span></p><span style="background-color: white; font-family: Arial, Helvetica, sans-serif; font-size: small;"></span><div style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small;"><p class="MsoNormal" style="margin: 0px;"><b style="font-family: arial, sans-serif;">Dictado en:</b><span style="font-family: arial, sans-serif;"> Departamento de Ciencias Geológicas --- (Pabellón 2, Ciudad Universitaria, Buenos Aires). Google maps </span><a data-saferedirecturl="https://www.google.com/url?q=https://www.google.com/maps/place/Departamento%2Bde%2BGeolog%25C3%25ADa/@-34.5418281,-58.4426073,15z/data%3D!4m5!3m4!1s0x0:0x204de95207c75963!8m2!3d-34.5418281!4d-58.4426073&source=gmail&ust=1662120093167000&usg=AOvVaw3wLDnjS4iIiEGiddSLkYlY" href="https://www.google.com/maps/place/Departamento+de+Geolog%C3%ADa/@-34.5418281,-58.4426073,15z/data=!4m5!3m4!1s0x0:0x204de95207c75963!8m2!3d-34.5418281!4d-58.4426073" style="color: #1155cc; font-family: arial, sans-serif;" target="_blank">link</a><span style="font-family: arial, sans-serif;">.</span><br /></p></div><span style="background-color: white; font-family: Arial, Helvetica, sans-serif; font-size: small;"><b>Modalidad: Presencial<br />Forma de Evaluación</b>: Desarrollo de un proyecto final.<br /></span><span style="background-color: white; font-family: Arial, Helvetica, sans-serif; font-size: small;"><b>Cupo máximo:</b> 20 alumnas/os.<br /></span><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><span style="color: black;"><b>Destinado a</b>: Estudiantes</span> de postgrado, profesores/as y profesionales de geología, geofísica y disciplinas afines dentro de las Ciencias de la Tierra.</p><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><b>Sitio web</b>: <a data-saferedirecturl="https://www.google.com/url?q=https://github.com/Esteban82/FCEN-2022&source=gmail&ust=1662120093167000&usg=AOvVaw0CKZ2zNQIOVyyHXZEP9w9u" href="https://github.com/Esteban82/FCEN-2022" style="color: #1155cc;" target="_blank">https://github.com/<wbr></wbr>Esteban82/FCEN-2022</a></p><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><b><span style="font-family: Arial, "sans-serif";"></span></b></p><b style="background-color: white; font-family: Arial, Helvetica, sans-serif; font-size: small;"><b>Programa: </b><span style="color: #222222; font-weight: normal;"><a data-saferedirecturl="https://www.google.com/url?q=https://github.com/Esteban82/FCEN-2022/blob/main/Programa.md&source=gmail&ust=1662120093167000&usg=AOvVaw07c9QXv-fuQIrLYiZI_JkO" href="https://github.com/Esteban82/FCEN-2022/blob/main/Programa.md" style="color: #1155cc;" target="_blank">https://github.com/<wbr></wbr>Esteban82/FCEN-2022/blob/main/<wbr></wbr>Programa.md</a></span></b><br /><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><span style="color: black;"><b>Destinado a</b>: Estudiantes</span> de postgrado, profesores/as y profesionales de geología, geofísica y disciplinas afines dentro de las Ciencias de la Tierra.</p><div><p class="MsoNormal" style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small; margin: 0px;"><b style="background-color: transparent; font-family: inherit;">Importante</b><span style="background-color: transparent; font-family: inherit;">: En el curso se trabajará con computadoras. Es muy recomendable que cada uno pueda asistir con una notebook para trabajar.</span></p></div><span style="font-size: x-small;"><a href="https://github.com/Esteban82/FCEN-2022#dirigido-a"></a></span><div style="background-color: white; font-family: Arial, Helvetica, sans-serif; font-size: small;"><b><span style="font-size: medium;"><br /></span></b></div><div style="background-color: white; font-family: Arial, Helvetica, sans-serif; font-size: small;"><b><span style="font-size: medium;">Inscripción:</span></b></div><b><div style="background-color: white; font-family: Arial, Helvetica, sans-serif; font-size: small;"><span style="color: #222222; font-weight: normal;"><br /></span></div><div style="background-color: white; font-family: Arial, Helvetica, sans-serif; font-size: small;"><span style="color: #222222; font-weight: normal;">Ver información, requisitos y aranceles para inscripción formal en </span><a data-saferedirecturl="https://www.google.com/url?q=https://exactas.uba.ar/ensenanza/cursos-de-posgrado/&source=gmail&ust=1662120093167000&usg=AOvVaw0inwvPQjV2ieT5BvAiBsPP" href="https://exactas.uba.ar/ensenanza/cursos-de-posgrado/" rel="noreferrer" style="color: #1155cc; font-weight: normal;" target="_blank">https://exactas.uba.ar/<wbr></wbr>ensenanza/cursos-de-posgrado/</a><br /></div></b><div style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small;"><b>Informes por aranceles</b>: Secretaría de Posgrado , Pabellón 2, Planta Baja, Horario: 9-13 y 14 a 19.00 hs. Tel: 5285-8060/61/56.</div><p><b style="color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small;">Informes y Pre-Inscripción, consultas por email al profesor: </b><a href="mailto:federico.esteban@gmail.com" style="color: #1155cc; font-family: Arial, Helvetica, sans-serif; font-size: small;" target="_blank">federico.esteban@<wbr></wbr>gmail.com</a> </p>Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0tag:blogger.com,1999:blog-1509523996050466325.post-53730940822925781132021-11-04T11:43:00.003-03:002021-11-04T11:46:46.342-03:00Animación. Vuelo de Ushuaia a La Quiaca<p> </p><div class="separator" style="clear: both; text-align: center;"><iframe allowfullscreen="" class="BLOG_video_class" height="412" src="https://www.youtube.com/embed/ovQE_LnWn6U" width="768" youtube-src-id="ovQE_LnWn6U"></iframe></div><p></p><span><a name='more'></a></span><p>Simulación de vuelo Ushuaia (-68.304444/-57.807222) a La Quiaca (-65.6/-22.105556) a una altura de 160 km con frames cada 5 km. <br /></p><p><b>Nota: Usar GMT 6.3 o posterior. </b><br /></p><p></p><p></p><p>Adaptado de <a href="https://docs.generic-mapping-tools.org/dev/animations/anim04.html">animación 4</a>.<br /></p><p>Datos: <a href="https://docs.generic-mapping-tools.org/dev/datasets/remote-data.html#global-earth-day-night-images">Blue Marble</a> </p><p>Imagen de Fondo: <a href="https://cdni.rt.com/actualidad/public_images/2019.12/article/r/1104/90/jpeg/5dfd168159bf5b27737f86fb.jpg">1.jpg</a><br /></p><p><br /></p><p></p><p>Script completo:<br /></p><p><br /></p><p></p><blockquote>#!/usr/bin/env bash<br /><br /># 1. Crear archivos para la animacion<br />cat << 'EOF' > pre.sh<br />gmt begin<br /> gmt project -C-68.304444/-57.807222 -E-65.6/-22.105556 -G5 -Q > flight_path.txt<br /> gmt image 1.jpg -D0/0+w24c -X0 -Y0<br />gmt end<br />EOF<br />cat << 'EOF' > main.sh<br />gmt begin<br /> gmt grdimage -JG${MOVIE_COL0}/${MOVIE_COL1}/${MOVIE_WIDTH}+du+z160+a3+t55+v36 -Rg @earth_day -Xc -Y-5c<br /> gmt coast -N1/thinner -N2/faint<br />gmt end<br />EOF<br /># 2. Crear la animacion<br />gmt movie main.sh -Chd -Nanim04 -Tflight_path.txt -Sbpre.sh -H2 -Vi -M0,png -Fmp4</blockquote><br /><p></p>Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0Buenos Aires, CABA, Argentina-34.6036844 -58.3815591-62.913918236178844 -93.5378091 -6.2934505638211533 -23.225309099999997tag:blogger.com,1999:blog-1509523996050466325.post-57264945160721102262021-10-18T09:36:00.001-03:002021-10-18T09:36:07.022-03:00Animación. Orbitando sobre las Islas Malvinas<p></p><div class="separator" style="clear: both; text-align: center;"><iframe allowfullscreen="" class="BLOG_video_class" height="366" src="https://www.youtube.com/embed/Tp6vyOZo7Uk" width="684" youtube-src-id="Tp6vyOZo7Uk"></iframe></div><span><a name='more'></a></span><br /><a href="https://www.dropbox.com/s/cji2kcanyjx4v12/Malvinas.sh?dl=0">Bash script</a><br /><p></p>Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0Buenos Aires, CABA, Argentina-34.6036844 -58.3815591-62.913918236178844 -93.5378091 -6.2934505638211533 -23.225309099999997tag:blogger.com,1999:blog-1509523996050466325.post-3025906069994857532021-09-26T20:16:00.005-03:002022-01-27T11:01:48.284-03:00Webinar "Programación y cartografía: Enséñale a tu PC a hacer los mapas que queres"<p></p><div class="separator" style="clear: both; text-align: center;"><iframe allowfullscreen="" class="BLOG_video_class" height="480" src="https://www.youtube.com/embed/xfEEWGu_jBo" width="670" youtube-src-id="xfEEWGu_jBo"></iframe></div><br /><div class="separator" style="clear: both; text-align: center;"></div><p></p><p></p><a name='more'></a><p></p><p></p><p>Ver repositorio en github para mas información: <a href="https://github.com/Esteban82/AcademiaGISArgentina" target="_blank">Link</a>.<br /></p><p>Webinar <span title="Editado">gratuito con <a class="notranslate" href="https://www.instagram.com/academiagisargentina/" tabindex="0">@academiagisargentina</a>.
Se realizó una demostración del uso del software Generic Mapping Tools
(GMT) mediante el uso de línea de comandos y scripts. En la misma se
elaboraron en vivo mapas políticos (con líneas de costas, ríos, lagos,
límite de países y provincias), físicos y satelitales con distintas
proyecciones geográficas (Mercator, Mollweide, UTM, estereográfica,
etc).<br /><br />La
disertación estuve dirigida a un público general y a la comunidad
educativa interesada en la elaboración de mapas de gran calidad gráfica.
No se requieren conocimientos previos de informática ni programación.<br /><br />En
el siguiente <a href="https://www.youtube.com/watch?v=1bPMIN7noTI">link</a> pueden ver como instalar los programas usados (Git for
windows y GTM, el Visual Studio Code no es necesario). <br /></span></p><br />Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0Buenos Aires, CABA, Argentina-34.6036844 -58.3815591-62.913918236178844 -93.5378091 -6.2934505638211533 -23.225309099999997tag:blogger.com,1999:blog-1509523996050466325.post-8836337912774313052021-08-23T12:35:00.001-03:002021-08-23T12:35:53.756-03:00Próximo curso. Asociación Geológica Argentina. Septiembre - Diciembre 2021<p></p><div class="separator" style="clear: both; text-align: center;"><a href="https://lh3.googleusercontent.com/-hMnvtpJnLvo/YSO9o6-xHhI/AAAAAAAAaUY/3xBd18wJaoItOZFYYgppKxQIpo3vithfwCLcBGAsYHQ/flyer_curso_geoinform%25C3%25A1tica_2021.png" style="margin-left: 1em; margin-right: 1em;"><img alt="" data-original-height="1152" data-original-width="2048" height="360" src="https://lh3.googleusercontent.com/-hMnvtpJnLvo/YSO9o6-xHhI/AAAAAAAAaUY/3xBd18wJaoItOZFYYgppKxQIpo3vithfwCLcBGAsYHQ/w640-h360/flyer_curso_geoinform%25C3%25A1tica_2021.png" width="640" /></a></div><p></p><p><span style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small;">Programa detallado del curso: </span><a href="https://www.dropbox.com/s/he7ahi2dff99upt/Programa_Curso_AGA_2021.pdf?dl=0" style="font-family: Arial, Helvetica, sans-serif; font-size: small;">link</a><span style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small;">.</span></p><div style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small;">Cualquier duda pueden escribir a <a href="mailto:secretariageologicaargentina@gmail.com" style="color: #1155cc;" target="_blank">secretariageologicaargentina@<wbr></wbr>gmail.com</a>. </div><div style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small;"><br /></div><div style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small;">Formulario de pre-inscripción: <br /></div><div style="background-color: white; color: #222222; font-family: Arial, Helvetica, sans-serif; font-size: small;"><b><span style="color: #990000;"><a data-saferedirecturl="https://www.google.com/url?q=https://docs.google.com/forms/d/e/1FAIpQLSerYPdTj2469rA0MVaZkEy4jq4RitDPKWd5vexXYhnsmhrOGQ/viewform&source=gmail&ust=1629818166257000&usg=AFQjCNHkbBmpzK3THtJD8xdYpLjSJVJmLQ" href="https://docs.google.com/forms/d/e/1FAIpQLSerYPdTj2469rA0MVaZkEy4jq4RitDPKWd5vexXYhnsmhrOGQ/viewform" style="color: #1155cc;" target="_blank">https://docs.google.com/forms/<wbr></wbr>d/e/<wbr></wbr>1FAIpQLSerYPdTj2469rA0MVaZkEy4<wbr></wbr>jq4RitDPKWd5vexXYhnsmhrOGQ/<wbr></wbr>viewform</a></span></b></div>Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0tag:blogger.com,1999:blog-1509523996050466325.post-84720343768317647722021-08-13T10:58:00.008-03:002021-08-13T11:19:51.270-03:00Animación goles de Messi en el planeta Tierra (2004-2021).<p><br /></p><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-qfaGRsENf4g/YQBAYvA4qDI/AAAAAAAAaJ4/bzc5IJJsQp0dC7lzJZ6sxp0ZlogSCKMlQCLcBGAsYHQ/s1920/GolesMessi_World_SinMarca.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1080" data-original-width="1920" height="360" src="https://1.bp.blogspot.com/-qfaGRsENf4g/YQBAYvA4qDI/AAAAAAAAaJ4/bzc5IJJsQp0dC7lzJZ6sxp0ZlogSCKMlQCLcBGAsYHQ/w640-h360/GolesMessi_World_SinMarca.png" width="640" /></a></div><p></p><div class="separator" style="clear: both; text-align: center;"><iframe allowfullscreen="" class="BLOG_video_class" height="369" src="https://www.youtube.com/embed/BFxq9xkrliU" width="657" youtube-src-id="BFxq9xkrliU"></iframe></div><p><span><span></span></span></p><a name='more'></a>Animación de los goles de Messi simulando sismos. <br /><p></p><p><a href="https://www.dropbox.com/s/153o51lx762rl9x/GolesMessi_World.sh?dl=0">Script Bash</a></p><p> </p><p>Descripción: Próximamente <br /></p><p></p><p></p><p><span></span></p><!--more--><p></p>Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com2Buenos Aires, CABA, Argentina-34.6036844 -58.3815591-62.913918236178844 -93.5378091 -6.2934505638211533 -23.225309099999997tag:blogger.com,1999:blog-1509523996050466325.post-43435834061924258862021-05-04T00:27:00.005-03:002021-11-24T14:38:18.416-03:00Animacion 14. Mecanismos focales 1981-2021<p> </p><div class="separator" style="clear: both; text-align: center;"><div class="separator" style="clear: both; text-align: center;"><a href="https://lh3.googleusercontent.com/-G3gEtqR9NAQ/YJC8E-5_68I/AAAAAAAAZSc/VL4v0cgjDHwSRiVrknYv5UFYRV_Hd49YQCLcBGAsYHQ/anim14_1981-2021.png" style="margin-left: 1em; margin-right: 1em;"><img alt="" data-original-height="1080" data-original-width="1920" height="360" src="https://lh3.googleusercontent.com/-G3gEtqR9NAQ/YJC8E-5_68I/AAAAAAAAZSc/VL4v0cgjDHwSRiVrknYv5UFYRV_Hd49YQCLcBGAsYHQ/w640-h360/anim14_1981-2021.png" width="640" /></a></div><br /></div><div class="separator" style="clear: both; text-align: center;"></div><span></span><p></p><span><a name='more'></a></span><p><span>Se muestra como crear una animación con mecanismos focales entre 1981 y 2021 con datos del Global CMT. La animación muestra en mapa y en sección los mecanismos en la zona de subducción plana del centro de Chile y el oeste de Argentina. </span></p><p><span>Script en: <a href="https://docs.generic-mapping-tools.org/dev/animations/anim14.html.">https://docs.generic-mapping-tools.org/dev/animations/anim14.html</a></span></p><p><span>Video en Youtube:</span></p><p><span></span></p><div class="separator" style="clear: both; text-align: center;"><iframe allowfullscreen="" class="BLOG_video_class" height="305" src="https://www.youtube.com/embed/02kEykwrTIc" width="320" youtube-src-id="02kEykwrTIc"></iframe></div><p><br /></p><p>GMT Catalog Citation Information<br />Dziewonski, A. M., T.-A. Chou and J. H. Woodhouse, Determination of earthquake source parameters from waveform data for studies of global and regional seismicity, J. Geophys. Res., 86, 2825-2852, 1981. doi:10.1029/JB086iB04p02825<br />Ekström, G., M. Nettles, and A. M. Dziewonski, The global CMT project 2004-2010: Centroid-moment tensors for 13,017 earthquakes, Phys. Earth Planet. Inter., 200-201, 1-9, 2012. doi:10.1016/j.pepi.2012.04<br /><span><span>https://www.globalcmt.org/</span></span></p><p></p><p><span> </span></p>Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0Buenos Aires, CABA, Argentina-34.6036844 -58.3815591-62.913918236178844 -93.5378091 -6.2934505638211533 -23.225309099999997tag:blogger.com,1999:blog-1509523996050466325.post-59341817503067444742021-03-04T00:24:00.006-03:002021-03-04T00:24:47.510-03:00EJ3.9. Gradiente de colores para figuras<p></p><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-7IXkdVitJww/YEBNEkgGygI/AAAAAAAAZJA/nCR4V9JBA6IfH99qrKKV-niAvC11q5l9wCLcBGAsYHQ/s2048/Color_Gradient_Fill.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="2042" data-original-width="2048" height="638" src="https://1.bp.blogspot.com/-7IXkdVitJww/YEBNEkgGygI/AAAAAAAAZJA/nCR4V9JBA6IfH99qrKKV-niAvC11q5l9wCLcBGAsYHQ/w640-h638/Color_Gradient_Fill.png" width="640" /></a></div><p></p><p></p><span><a name='more'></a></span><p></p><p>Se muestra como hacer para crear un gradiente de color para el fondo de figuras usando makecpt, grdmath y grdedit. </p><p><u>Scripts con sintaxis moderna: </u><br /></p><p><span> </span><a href="https://www.dropbox.com/s/vujo4qyaa5j2gcv/EJ3.9_Color_Gradient_Fill.sh?dl=0">shell</a></p><p><span> </span><a href="https://www.dropbox.com/s/is2z5vtkasepq49/EJ3.9_Color_Gradient_Fill.bat?dl=0">bat</a><br /> </p><p>1. <b>makecpt</b>: Se utiliza para crear una paleta de colores (CPT) continuo. En este caso con los colores blanco a gris y para valores de 0 a 1.<br /></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt makecpt -Cwhite,50 -T0/1 -Z </span></span></p><p></p><p> </p><p>2. <b>gradmath</b>: Permite crear superficies 2d para colorear con el CPT. Crear grilla con grdmath. Se utiliza una resolución de 0.01 unidades (-I) lo suficientemente chico para evitar un efecto de pixelado. <br /></p><p><u>Nota</u>: Todas las grillas tienen valores entre 0 y 1 (-R0/1/0/1). <br /></p><p> </p><p></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdmath -I0.01 Y = %GRD%</span></span></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdmath -I0.01 X = %GRD%</span></span></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdmath -I0.01 X Y MUL = %GRD%</span></span></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdmath -I0.01 X Y ADD NORM = %GRD%<br /></span></span></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdmath -I0.01 X Y ADD NORM = %GRD%<br />gmt grdedit %GRD% -Er<br /></span></span></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdmath -I0.01 0.5 0.5 CDIST NORM = %GRD%</span></span></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdmath -I0.01 0.5 0.5 CDIST NORM NEG 1 ADD = %GRD%<br /></span></span></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdmath -I0.01 1 1 CDIST NORM = %GRD%<br /></span></span></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdmath -I0.01 1 0 CDIST NORM = %GRD%</span></span><br /></p><p> <br /></p><p><b>3. grdimage</b>: Se utliza para crear el fondo degradado. </p><p><span style="background-color: #cccccc;">gmt grdimage %GRD% -C</span><br /></p><p><b>4. Graficar encima</b>: Una vez graficado el fondo de la figura, se puedan graficar datos. Hay que tener el cuidado de modificar la región de la figura y de mantener sus dimensiones. <br /></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt histogram -Bxaf -Byaf @v3206_06.txt -R-6000/0/0/30 -Gorange -W1p -Z1 -T250</span></span><br /></p><p><br /></p>Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0Buenos Aires, CABA, Argentina-34.6036844 -58.3815591-62.913918236178844 -93.5378091 -6.2934505638211533 -23.225309099999997tag:blogger.com,1999:blog-1509523996050466325.post-74925590419724417042020-11-26T00:00:00.004-03:002020-11-27T08:46:33.410-03:00Repositorio Curso "Geoinformática aplicada a la cartografía multitemática" <p>Les comparto el link del repositorio del Curso "Geoinformática aplicada a la cartografía
multitemática" brindado por la Asociación Geológica Argentina (AGA) en
agosto y septiembre de 2020.</p>
<p>En el mismo hay ejercicios (muchos de ellos explicados en este blog), guías, datos y links a las clases en youtube (10 clases de 2 hs cada una).</p>
<p>Link: <a href="https://github.com/Esteban82/AGA-Geoinformatica-2020">https://github.com/Esteban82/AGA-Geoinformatica-2020</a></p><p><br /></p><p> <br /></p>Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com1tag:blogger.com,1999:blog-1509523996050466325.post-83468256153589307522020-10-21T14:47:00.009-03:002022-11-24T16:05:20.641-03:00Estilo bibliográfico de la RAGA para Mendeley y Zotero.<div dir="ltr" style="text-align: left;" trbidi="on">
Como agregar el estilo bibliográfico (csl) de la Revista de la Asociación Geológica Argentina (RAGA) en Mendeley.</div><div dir="ltr" style="text-align: left;" trbidi="on"><span><a name='more'></a></span></div><div dir="ltr" style="text-align: left;" trbidi="on">
</div><div dir="ltr" style="text-align: left;" trbidi="on"></div><div dir="ltr" style="text-align: left;" trbidi="on">Los archivos Citation Style Languague (<a href="https://en.wikipedia.org/wiki/Citation_Style_Language">csl</a>) sirven para dar formato a la citas y referencias bibliográficas en software como Mendeley o Zotero (<a href="https://es.wikipedia.org/wiki/Gestor_de_referencias_bibliogr%C3%A1ficas">gestores de referencias bibliográficas</a>) según los requisitos de cada <i>journal</i>. <br /><span face="sans-serif" style="background-color: white; color: #202122; font-size: 14px;"><br />El </span><b style="background-color: white; color: #202122; font-family: sans-serif; font-size: 14px;">Citation Style Language</b><span face="sans-serif" style="background-color: white; color: #202122; font-size: 14px;"> (</span><b style="background-color: white; color: #202122; font-family: sans-serif; font-size: 14px;">CSL</b><span face="sans-serif" style="background-color: white; color: #202122; font-size: 14px;">) es un lenguaje abierto basado en </span><a href="https://es.wikipedia.org/wiki/Extensible_Markup_Language" style="background: none rgb(255, 255, 255); color: #0645ad; font-family: sans-serif; font-size: 14px; text-decoration-line: none;" title="Extensible Markup Language">XML</a><span face="sans-serif" style="color: #202122;"><span style="background-color: white; font-size: 14px;">.</span></span></div><div dir="ltr" style="text-align: left;" trbidi="on">
<br />
En el siguiente link esta disponible el formato de la Revista de la Asociación Geológica Argentina (<a href="http://ppct.caicyt.gov.ar/raga">RAGA</a>). <b>Por ahora es una versión Beta </b>(cualquier error por favor avisar).</div><div dir="ltr" style="text-align: left;" trbidi="on"><u><br /></u>
</div><div dir="ltr" style="text-align: left;" trbidi="on"><u>Link archivo csl</u>:
<a href="https://csl.mendeley.com/styles/544939841/revista-de-la-asociacion-geologica-argentina">https://csl.mendeley.com/styles/544939841/revista-de-la-asociacion-geologica-argentina</a></div><div dir="ltr" style="text-align: left;" trbidi="on"><br />
<br />
<b>A. Como agregar el estilo bibliográfico en Mendeley</b><br />
<br />
1. Seleccionar <b>View</b>, <b>Citation Style</b>, <b>More Styles...</b><br />
<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-sch4ypxjzgI/XR1pXjig2mI/AAAAAAAAUCc/nohNwlArZZIYN1VZ3BjDPq-BtcCky53_gCLcBGAs/s1600/1.jpeg" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="328" data-original-width="579" height="226" src="https://1.bp.blogspot.com/-sch4ypxjzgI/XR1pXjig2mI/AAAAAAAAUCc/nohNwlArZZIYN1VZ3BjDPq-BtcCky53_gCLcBGAs/s400/1.jpeg" width="400" /></a></div>
<br />
2. En la ventana que se abre, ir a la pestaña <b>Get More Styles</b> y copiar el link del archivo csl en el campo <b>Download Style</b>, y hacer clic en <b>Download</b>.</div><div dir="ltr" style="text-align: left;" trbidi="on"><br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-EFn1wMXJ61s/XR1pXjTcUCI/AAAAAAAAUCg/j7LAEOUxpYsgitbW1BaII6ENMQ99TwmSgCEwYBhgL/s1600/2.jpeg" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="501" data-original-width="670" height="298" src="https://1.bp.blogspot.com/-EFn1wMXJ61s/XR1pXjTcUCI/AAAAAAAAUCg/j7LAEOUxpYsgitbW1BaII6ENMQ99TwmSgCEwYBhgL/s400/2.jpeg" width="400" /></a></div>
<br />
3. En la pestaña debe aparecer el estilo de la "<i>Revista de la Asociación Geológica Argentina</i>". En <b>Citation and Bibliography Language</b> seleccionar Español.<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/--HIGuOgbBb0/XR1pXSejMrI/AAAAAAAAUCk/VqPuvjr9g30-E8qFP8WRZ_W-g6cCK2SiwCEwYBhgL/s1600/3.jpeg" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="498" data-original-width="683" height="291" src="https://1.bp.blogspot.com/--HIGuOgbBb0/XR1pXSejMrI/AAAAAAAAUCk/VqPuvjr9g30-E8qFP8WRZ_W-g6cCK2SiwCEwYBhgL/s400/3.jpeg" width="400" /></a></div>
<br /></div><div dir="ltr" style="text-align: left;" trbidi="on"><b><br /></b></div><div dir="ltr" style="text-align: left;" trbidi="on">PD: 24-Nov-2022</div><div dir="ltr" style="text-align: left;" trbidi="on"><br /></div><div dir="ltr" style="text-align: left;" trbidi="on"><b>B. Como agregar el estilo bibliográfico en Zotero</b></div><div dir="ltr" style="text-align: left;" trbidi="on"><b><br /></b></div><div dir="ltr" style="text-align: left;" trbidi="on">1. Seleccionar <b>Editar</b>, <b>Preferencias</b>.</div><div dir="ltr" style="text-align: left;" trbidi="on"><br /></div><div dir="ltr" style="text-align: left;" trbidi="on"><br /></div><div class="separator" style="clear: both; text-align: center;"><a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgh2bZgCxdI1M2gsMcWcxCWaW8F85Y37tjr6fgiKZ41YAwHoX1nr0pu0yyPtjSXTy7ZlB6vKq4fo2gaXzAq57RJS-aPqJBUGcVh_dNt9fag-QhQC3qQYX7tB4Uv_wQV8PVBi3jKpVEsWVGqWezLswAc24JscQ8KM3sjMob7m-Umtb4rQaySPQJuK2EjnA/s1247/1_CLS.jpeg" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="930" data-original-width="1247" height="478" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgh2bZgCxdI1M2gsMcWcxCWaW8F85Y37tjr6fgiKZ41YAwHoX1nr0pu0yyPtjSXTy7ZlB6vKq4fo2gaXzAq57RJS-aPqJBUGcVh_dNt9fag-QhQC3qQYX7tB4Uv_wQV8PVBi3jKpVEsWVGqWezLswAc24JscQ8KM3sjMob7m-Umtb4rQaySPQJuK2EjnA/w640-h478/1_CLS.jpeg" width="640" /></a></div><br /><div dir="ltr" style="text-align: left;" trbidi="on">2. En la ventana que se abre, seleccionar la pestaña <b>Citar</b>, y luego hacer click en el "<b>+</b>".</div><div dir="ltr" style="text-align: left;" trbidi="on"><br /></div><div class="separator" style="clear: both; text-align: center;"><a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhOg7sr3h6miVqTgInV89d3I6a87UewHOPYQhlnqkt7q4v3X5Ke0M4n2dR1TTY9As5Nu06jzuyX-96u3aSUInCQ4TRdqMsgUOccdS9ZOTZ-dizSYFzSygwTXUBrGXie6IbibecZJv4rbGJWq_lSR5mgx4LwkcOC7MoVsXLPezgEsWdcOoEH5KM02gwVbg/s1245/2_AbrirCls.jpeg" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="928" data-original-width="1245" height="478" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhOg7sr3h6miVqTgInV89d3I6a87UewHOPYQhlnqkt7q4v3X5Ke0M4n2dR1TTY9As5Nu06jzuyX-96u3aSUInCQ4TRdqMsgUOccdS9ZOTZ-dizSYFzSygwTXUBrGXie6IbibecZJv4rbGJWq_lSR5mgx4LwkcOC7MoVsXLPezgEsWdcOoEH5KM02gwVbg/w640-h478/2_AbrirCls.jpeg" width="640" /></a></div><div><br /></div>3. En la ventana que se abre, elegir el archivo <b>csl</b> descargado.<div><br /><div class="separator" style="clear: both; text-align: center;"><a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiGs9Z7YUc803D2MHoUZ1CkxNh5zTiiH7O_4fakDLm20km96EM8OQofsT_7OkaATE28AXWBAEr2sOQGFfB514njaAduVj_-RijdZWkgCgrvC0qrWZ9NBi7iI1k5A2EGxgObX2M7XebM_-d3u55Kulv0IcdNyEjVb0YQgdk0z6rgw3dNZs4HGAZANnYZzg/s1242/3_AbrirCLS.jpeg" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="929" data-original-width="1242" height="478" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEiGs9Z7YUc803D2MHoUZ1CkxNh5zTiiH7O_4fakDLm20km96EM8OQofsT_7OkaATE28AXWBAEr2sOQGFfB514njaAduVj_-RijdZWkgCgrvC0qrWZ9NBi7iI1k5A2EGxgObX2M7XebM_-d3u55Kulv0IcdNyEjVb0YQgdk0z6rgw3dNZs4HGAZANnYZzg/w640-h478/3_AbrirCLS.jpeg" width="640" /></a></div><br /><br /><div dir="ltr" style="text-align: left;" trbidi="on"><br /></div><div dir="ltr" style="text-align: left;" trbidi="on">Mas informacion:</div><div dir="ltr" style="text-align: left;" trbidi="on">https://en.wikipedia.org/wiki/Citation_Style_Language</div><div dir="ltr" style="text-align: left;" trbidi="on">https://github.com/zotero/bundled-styles</div>
</div>Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com4Argentina-38.416097 -63.616671999999987-66.6670613196018 -98.772922 -10.165132680398195 -28.460421999999987tag:blogger.com,1999:blog-1509523996050466325.post-60297830445462107752020-10-12T00:00:00.028-03:002021-03-01T20:42:40.134-03:0023.1. Armónicos esféricos. Topografía lunar. <div dir="ltr" style="text-align: left;" trbidi="on">
<div dir="ltr"><div class="separator" style="clear: both; text-align: center;"><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-mA9Tf97SHik/X4Ibp5fhw-I/AAAAAAAAY0c/8LrTQW2eMeoaqKBGZj35YrU_dWAeCEbmQCLcBGAsYHQ/s1966/MoonTopo%252BEspectro_1440_0.125.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1942" data-original-width="1966" height="632" src="https://1.bp.blogspot.com/-mA9Tf97SHik/X4Ibp5fhw-I/AAAAAAAAY0c/8LrTQW2eMeoaqKBGZj35YrU_dWAeCEbmQCLcBGAsYHQ/w640-h632/MoonTopo%252BEspectro_1440_0.125.png" width="640" /></a></div></div></div><div dir="ltr"><span><a name='more'></a></span></div><div dir="ltr">En ciencias de la tierra frecuentemente se utilizan los armónicos esféricos para modelar distintos parámetros (topografía, geoide, campo gravitatorio/magnetométrico) a escala planetaria. En este ejemplo se muestra como crear un mapa topográfico lunar a partir de la tabla de coeficientes armónicos de Wieczorek (2015).</div><div dir="ltr"> </div><div dir="ltr"><a href="https://www.dropbox.com/s/srkba7rtq5t2wju/MoonTopo%2BEspectro.bat?dl=0" target="_blank">Script</a><br /></div><div dir="ltr"> </div><div dir="ltr"><br /></div><div dir="ltr"><b>0. Datos:</b> </div><div dir="ltr"><ul style="text-align: left;"><li>MoonTopo2600p.shape: Armónicos esféricos de la luna normalizados 4π hasta el orden 2600. <a href="https://zenodo.org/record/3870924">Link zenodo</a>.<b> <br /></b></li></ul><p><b>1. Armónicos esféricos:</b> Tabla de datos conformado por 4 columnas con los datos de grado (L), orden (M) y coeficientes de coseno (C) y de seno (S) hasta el grado 2600 que totalizan casi 3,4 millones de registros (<span class="qv3Wpe" id="cwos">3.383.901).<br /></span></p><p><span class="qv3Wpe" id="cwos"><b>Nota</b>: El grado 0 (orden 0) describe una esfera con el radio promedio lunar (1737 km).</span></p><span style="background-color: #cccccc;"><span style="font-family: courier;"></span></span><p><span class="qv3Wpe" id="cwos"></span></p><table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto;"><tbody><tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-GjBlAepLKGw/X4EYmkcfPkI/AAAAAAAAYzQ/n-Q1mTNKUogr7x5fdP9wkjg7qld2IKZUgCLcBGAsYHQ/s421/MoonTopo_Armonicos_Esfericos.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="364" data-original-width="421" height="276" src="https://1.bp.blogspot.com/-GjBlAepLKGw/X4EYmkcfPkI/AAAAAAAAYzQ/n-Q1mTNKUogr7x5fdP9wkjg7qld2IKZUgCLcBGAsYHQ/w320-h276/MoonTopo_Armonicos_Esfericos.png" width="320" /></a></td></tr><tr><td class="tr-caption" style="text-align: center;"><span style="font-size: x-small;">Tabla de coeficientes de armónicos esféricos hasta el grado 5 y orden 5.</span><br /></td></tr></tbody></table><span class="qv3Wpe" id="cwos"></span><p></p><p><b>2. Elipsoide lunar: </b>Con el siguiente comando se establece el elipsoide de luna:<br /></p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmtset PROJ_ELLIPSOID Moon</span></span><p><b>3. Crear grilla</b>: con <a href="https://docs.generic-mapping-tools.org/6.1/sph2grd.html" target="_blank">sph2grd </a>se convierte la tabla de armónicos esféricos (%SPH%) en una grilla con los siguientes argumentos:</p><ul style="text-align: left;"><li>-G%CUT%: Nombre de la grilla creada.</li><li>-R: Región geográfica (o "selenográfica") del mapa.</li><li>-Ng: Indica normalización geodésica de los coeficientes (suma de los productos sobre la superficie = 4π).</li><li>-I: Resolución de la grilla creada en grados.</li><li>-F: Filtrado de los armónicos esféricos según su grado (l) utilizados para crear la grilla. <br /></li></ul><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt sph2grd -R "%SPH%" -G%CUT% -Ng -I0.25 -F1/1/1410/1440</span></span></p><p> </p><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-W3D6pV491bg/X4J2I7esmdI/AAAAAAAAY0w/mAXlmc566lgKLEAV_se5UvNSsRUSQs2gACLcBGAsYHQ/s1776/MoonTopo_1440_0.125.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1130" data-original-width="1776" height="255" src="https://1.bp.blogspot.com/-W3D6pV491bg/X4J2I7esmdI/AAAAAAAAY0w/mAXlmc566lgKLEAV_se5UvNSsRUSQs2gACLcBGAsYHQ/w400-h255/MoonTopo_1440_0.125.png" width="400" /></a></div><p></p><p>Es importante seleccionar cuidadosamente un filtro para los armónicos esféricos y la resolución de la grilla. </p><p><b>A. Filtrado de armónicos esféricos</b>: El filtro (-F) es opcional y puede ser del tipo <i>Cosine band-pass</i> (definiendo 4 parámetros) o <i>Gaussian band-pass</i> (2 parámetros). En los ejemplos siguientes se aplicó el <i>Cosine band-pass</i>.<br /></p><p><b>i. Filtrar orden 0:</b> Como el modelo corresponde a la forma de la luna, entonces si se incluye el orden 0 la grilla resultante corresponde a la distancia al centro de la luna (o radio lunar). Por lo tanto, si se filtra el grado 0, entonces se obtiene las diferencias de alturas con respecto al radio lunar promedio (el cual es utilizado como <a href="https://en.wikipedia.org/wiki/Topography_of_the_Moon#Reference_elevation" target="_blank">"nivel del mar" para los mapas topográficos de la luna</a>).<br /></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt sph2grd -R "%SPH%" -G%CUT% -Ng -I0.25 -F0/0/710/720</span></span><br /><br /></p><p></p><div class="separator" style="clear: both; text-align: center;"><b><a href="https://1.bp.blogspot.com/-O9-wO_oxuKs/X4Ev5VsNQpI/AAAAAAAAYzo/ZvnhLUyitLcRjrP6UcL2jhW5WTdUNN_UwCLcBGAsYHQ/s1776/MoonTopo_0-720_0.25.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1130" data-original-width="1776" height="255" src="https://1.bp.blogspot.com/-O9-wO_oxuKs/X4Ev5VsNQpI/AAAAAAAAYzo/ZvnhLUyitLcRjrP6UcL2jhW5WTdUNN_UwCLcBGAsYHQ/w400-h255/MoonTopo_0-720_0.25.png" width="400" /> </a></b></div><div class="separator" style="clear: both; text-align: center;"><b> </b></div><div class="separator" style="clear: both; text-align: left;"><b>ii. Filtrar ordenes mayores</b>: El mayor grado (L) de los armónicos esféricos determina la resolución espacial (λ; detalles de los rasgos) que corresponde a grados y km respectivamente según las fórmulas:</div><div class="separator" style="clear: both; text-align: left;"><b></b></div><div class="separator" style="clear: both; text-align: left;"><b></b></div><div class="separator" style="clear: both; text-align: left;"><b><br /></b></div><div class="separator" style="clear: both; text-align: left;">λ grados = 180/L </div><div class="separator" style="clear: both; text-align: left;">λ km = Radio*π / L (que equivale a 20.000 km/L para la Tierra y <span class="qv3Wpe" id="cwos">5.400 km para la Luna).</span><br /></div><div class="separator" style="clear: both; text-align: left;"> </div><div class="separator" style="clear: both; text-align: left;">Filtrar las órdenes mayores es muy recomendable. Utilizar todos los armónicos puede consumir
gran tiempo de procesamiento que luego pueden no observarse en la
figura final. Las siguientes figuras todas tiene la misma resolución de grilla (-I0.125) pero distinta resolución espacial (a partir de distintos filtros aplicados).<b><br /></b></div><div class="separator" style="clear: both; text-align: center;"><b> </b></div><b></b></div><div dir="ltr"><b></b></div><div dir="ltr"><b></b></div><div dir="ltr"><b></b></div><div dir="ltr"><b></b><p></p><p><b></b></p><table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto;"><tbody><tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-jMM06KPZjR4/X4H4X_DtXWI/AAAAAAAAY0E/CVUAgqg8tVA-iuEvYqoPAjipNPS3WJ-YwCLcBGAsYHQ/s1966/MoonTopo%252BEspectro_90_0.125.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1942" data-original-width="1966" height="395" src="https://1.bp.blogspot.com/-jMM06KPZjR4/X4H4X_DtXWI/AAAAAAAAY0E/CVUAgqg8tVA-iuEvYqoPAjipNPS3WJ-YwCLcBGAsYHQ/w400-h395/MoonTopo%252BEspectro_90_0.125.png" width="400" /></a></td></tr><tr><td class="tr-caption" style="text-align: center;"><span style="font-size: x-small;">Filtrado hasta el grado 90 (-F1/1/85/90).</span><br /></td></tr></tbody></table><div class="separator" style="clear: both; text-align: center;"><b></b></div><b><br /><table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto;"><tbody><tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-UCRxO6XsXV4/X4H4Va0uiDI/AAAAAAAAYz0/eZqECfIF9jwXToYCm52XeTpWLJazWChIwCLcBGAsYHQ/s1966/MoonTopo%252BEspectro_180_0.125.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1942" data-original-width="1966" height="395" src="https://1.bp.blogspot.com/-UCRxO6XsXV4/X4H4Va0uiDI/AAAAAAAAYz0/eZqECfIF9jwXToYCm52XeTpWLJazWChIwCLcBGAsYHQ/w400-h395/MoonTopo%252BEspectro_180_0.125.png" width="400" /></a></td></tr><tr><td class="tr-caption" style="text-align: center;"><span style="font-size: x-small;">Filtrado hasta el grado 180 (-F1/1/180/180).</span></td></tr></tbody></table><br /><table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto;"><tbody><tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-ZJ0i7aISxkQ/X4H4Wc2CarI/AAAAAAAAYz4/hP2l4Clo8RInRoYD2jPlJq3BVnf3oNtrACLcBGAsYHQ/s1966/MoonTopo%252BEspectro_360_0.125.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1942" data-original-width="1966" height="395" src="https://1.bp.blogspot.com/-ZJ0i7aISxkQ/X4H4Wc2CarI/AAAAAAAAYz4/hP2l4Clo8RInRoYD2jPlJq3BVnf3oNtrACLcBGAsYHQ/w400-h395/MoonTopo%252BEspectro_360_0.125.png" width="400" /></a></td></tr><tr><td class="tr-caption" style="text-align: center;"><span style="font-size: x-small;">Filtrado hasta el grado 360 (-F1/1/35<span></span>0/360).</span></td></tr></tbody></table><br /><table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto;"><tbody><tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-ZEJEFxRx9oI/X4H4XgWCoBI/AAAAAAAAY0A/XqucRn3CiO4Bt-Ekb2lnwgXR4Sf5dU5iACLcBGAsYHQ/s1966/MoonTopo%252BEspectro_720_0.125.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1942" data-original-width="1966" height="395" src="https://1.bp.blogspot.com/-ZEJEFxRx9oI/X4H4XgWCoBI/AAAAAAAAY0A/XqucRn3CiO4Bt-Ekb2lnwgXR4Sf5dU5iACLcBGAsYHQ/w400-h395/MoonTopo%252BEspectro_720_0.125.png" width="400" /></a></td></tr><tr><td class="tr-caption" style="text-align: center;"><span style="font-size: x-small;">Filtrado hasta el grado 90 (-F1/1/710/720).</span></td></tr></tbody></table><br /><table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto;"><tbody><tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-7guqsn_09ow/X4H4WoOhJNI/AAAAAAAAYz8/cPmLDCST9r08OZZfPfq68KglyRJBXOBTwCLcBGAsYHQ/s1966/MoonTopo%252BEspectro_1440_0.125.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1942" data-original-width="1966" height="395" src="https://1.bp.blogspot.com/-7guqsn_09ow/X4H4WoOhJNI/AAAAAAAAYz8/cPmLDCST9r08OZZfPfq68KglyRJBXOBTwCLcBGAsYHQ/w400-h395/MoonTopo%252BEspectro_1440_0.125.png" width="400" /></a></td></tr><tr><td class="tr-caption" style="text-align: center;"><span style="font-size: x-small;">Filtrado hasta el grado 1440 (-F1/1/1420/1440).</span></td></tr></tbody></table> <br /></b><p></p><p><b>B. Resolución de la grilla</b>: Se debe elegir cuidadosamente. Un valor grande puede resultar que la resolución espacial no se observe correctamente, hasta incluso se vea un efecto de pixelado. Un valor muy chico puede resultar en excesivos tiempos de procesamiento y con un detalle excesivo que no se observa en la figura final.</p><p><br /></p><table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto;"><tbody><tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-NszTO516k3Y/X3_FD2svCKI/AAAAAAAAYyQ/lm1Q_FOmEkALkIxTQiMvdw3-LtcH8KUmwCLcBGAsYHQ/s1776/MoonTopo_G720_1_1.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1130" data-original-width="1776" height="254" src="https://1.bp.blogspot.com/-NszTO516k3Y/X3_FD2svCKI/AAAAAAAAYyQ/lm1Q_FOmEkALkIxTQiMvdw3-LtcH8KUmwCLcBGAsYHQ/w400-h254/MoonTopo_G720_1_1.png" width="400" /></a></td></tr><tr><td class="tr-caption" style="text-align: center;"><span style="font-size: x-small;">Resolución de 1 grado de arco (-I1).</span></td></tr></tbody></table><p></p><table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody><tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-vKDmFV4t0OY/X3_FErIIsQI/AAAAAAAAYyU/MwdHg_BVm-E9GrIRUvnmFiK4ZgwFIezMQCLcBGAsYHQ/s1776/MoonTopo_G720_2_0.5.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1130" data-original-width="1776" height="254" src="https://1.bp.blogspot.com/-vKDmFV4t0OY/X3_FErIIsQI/AAAAAAAAYyU/MwdHg_BVm-E9GrIRUvnmFiK4ZgwFIezMQCLcBGAsYHQ/w400-h254/MoonTopo_G720_2_0.5.png" width="400" /></a></td></tr><tr><td class="tr-caption" style="text-align: center;"><span style="font-size: x-small;">Resolución de 30 minutos de arco (-I0.5).</span></td><td class="tr-caption" style="text-align: center;"><br /></td><td class="tr-caption" style="text-align: center;"><br /></td><td class="tr-caption" style="text-align: center;"><br /></td><td class="tr-caption" style="text-align: center;"><br /></td><td class="tr-caption" style="text-align: center;"><br /></td><td class="tr-caption" style="text-align: center;"><br /></td><td class="tr-caption" style="text-align: center;"><br /></td></tr></tbody></table></div><div dir="ltr"><table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody><tr><td class="tr-caption" style="text-align: center;"><br /></td></tr></tbody></table></div></div><div dir="ltr" style="text-align: left;" trbidi="on"></div><div dir="ltr" style="text-align: left;" trbidi="on"></div><table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto;"><tbody><tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-zxBhk5HRqNQ/X3_Ph_ojJeI/AAAAAAAAYy8/vcMA2iFgqEwdvJ-4BS0UJwCk8mOwDKiXwCLcBGAsYHQ/s1776/MoonTopo_G720_3_0.25.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1130" data-original-width="1776" height="255" src="https://1.bp.blogspot.com/-zxBhk5HRqNQ/X3_Ph_ojJeI/AAAAAAAAYy8/vcMA2iFgqEwdvJ-4BS0UJwCk8mOwDKiXwCLcBGAsYHQ/w400-h255/MoonTopo_G720_3_0.25.png" width="400" /></a></td></tr><tr><td class="tr-caption" style="text-align: center;"><span style="font-size: x-small;">Resolución de 15 minutos de arco (-I0.25).</span></td></tr></tbody></table><br /><table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto;"><tbody><tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-nN-IxFvlNKs/X3_PhKaYM4I/AAAAAAAAYy4/77qsoUTViCQHH0x1-ZAClGB5f61oDv9mgCLcBGAsYHQ/s1776/MoonTopo_G720_4_0.125.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1130" data-original-width="1776" height="255" src="https://1.bp.blogspot.com/-nN-IxFvlNKs/X3_PhKaYM4I/AAAAAAAAYy4/77qsoUTViCQHH0x1-ZAClGB5f61oDv9mgCLcBGAsYHQ/w400-h255/MoonTopo_G720_4_0.125.png" width="400" /></a></td></tr><tr><td class="tr-caption" style="text-align: center;"><span style="font-size: x-small;">Resolución de 7,5 minutos de arco (-I0.125).</span></td></tr></tbody></table><div dir="ltr" style="text-align: left;" trbidi="on"> <div dir="ltr"><p><b>Referencias</b>:</p><ul style="text-align: left;"><li><span class="ng-binding">Wieczorek,
Mark A. (2015). Spherical harmonic model of the shape of Earth's Moon:
MoonTopo2600p [Data set]. Zenodo. http://doi.org/10.5281/zenodo.3870924</span><span class="ng-binding"> </span></li><li><span class="ng-binding"> </span>Wieczorek, Mark A. Gravity and Topography of the
Terrestrial Planets, Treatise on Geophysics, 2nd edition, Oxford,
153-193, doi:10.1016/B978-0-444-53802-4.00169-X, 2015.</li></ul></div>
</div>
Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0tag:blogger.com,1999:blog-1509523996050466325.post-9224952255836165792020-09-28T00:10:00.001-03:002020-09-28T00:10:49.157-03:007.6. Anomalias gravimétricas con sombreado topografico<p> Graficar anomalías gravimétricas con efecto de sombreado topográfico.</p><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-EKj3P24AXCY/X26xLaLO4vI/AAAAAAAAX3U/NLdurYNG2FYbvc_Lw_Y0dSXb8xjXJJt8gCLcBGAsYHQ/s2048/EJ7.7_Mapa_DEM%252BGrav_AireLibre_2.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1256" data-original-width="2048" height="392" src="https://1.bp.blogspot.com/-EKj3P24AXCY/X26xLaLO4vI/AAAAAAAAX3U/NLdurYNG2FYbvc_Lw_Y0dSXb8xjXJJt8gCLcBGAsYHQ/w640-h392/EJ7.7_Mapa_DEM%252BGrav_AireLibre_2.png" width="640" /></a></div><p></p><a name="more"></a><p></p><p><span></span></p><a name='more'></a><p></p><p>Cuando se grafican grillas de anomalías gravimétricas (o cualquier otra grilla en general) puede ser difícil correlacionar estas anomalías con el relieve a pesar de incluir distintos elementos en el mapa (escala, coordenadas, límites, etc). Una forma de lograrlo (y sin sobrecargar el mapa) es agregar el efecto de sombreado (hill-shade) de una grilla topográfica de alta resolución.</p><p>Datos:</p><ul style="text-align: left;"><li>grav_29.1.nc (%FAA%): Grilla de anomalías de Aire Libre de Sandwell y Smith (versión 29.1) con resolución de 1 minuto de arco (<a href="https://topex.ucsd.edu/pub/global_grav_1min/grav_29.1.nc">link</a>).<br /></li></ul><ul style="text-align: left;"><li>GMRTv3_1.grd (%DEM%): Grilla topográfica de alta resolución (<i>ca</i>. 8 x 7 segundos de arco). </li></ul><div><div><p><a href="https://www.dropbox.com/s/93l1edxq7d2c21w/EJ7.6_Mapa_Dem%2BGrav_AireLibre.bat?dl=0" target="_blank">Script</a><br /></p><p><b>0. Mapa de anomalías gravimétricas</b>: Con efecto de sombreado calculado a partir de la misma grilla de gravimétrica adaptada a la región de Cuyo (ver <a href="http://7.1">7.1</a>).</p><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-W2exQzaJGaM/X3FSPYLezSI/AAAAAAAAX44/hlJEyqknNfIb_a99rDUtuw-8ueMLEivzwCLcBGAsYHQ/s2048/EJ7.6_Mapa_DEM%252BGrav_AireLibre_1.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1256" data-original-width="2048" height="392" src="https://1.bp.blogspot.com/-W2exQzaJGaM/X3FSPYLezSI/AAAAAAAAX44/hlJEyqknNfIb_a99rDUtuw-8ueMLEivzwCLcBGAsYHQ/w640-h392/EJ7.6_Mapa_DEM%252BGrav_AireLibre_1.png" width="640" /></a></div><br /><p></p><p><b>1. Grilla de Sombreado</b>: Con grdcut se recorta la grilla topográfica a la región del mapa y luego con grdgradient se cree la grilla (%SHADOW%) para agregar el efecto de sombreado (ver <a href="https://mapasgeologicos.blogspot.com/2016/03/52-mapa-topografico-ii.html" target="_blank">5.2</a> para más detalles).<br /></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdcut "%DEM%" -G%CUT2% -R</span></span></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdgradient %CUT2% -A270 -G%SHADOW% -Ne0.8</span></span></p><br /><b>2. Remuestrear:</b> Para poder graficar una grilla con efecto de sombreado es necesario que ambas grillas cubran la misma región y tengan la misma resolución (esto garantiza que para cada nodo de una grilla haya uno y sólo un un nodo en la otra). Esto se puede lograr con <a href="https://docs.generic-mapping-tools.org/latest/grdsample.html" target="_blank">grdsample</a> con los siguientes argumentos:<ul style="text-align: left;"><li>%FAA%: Grilla de anomalía de aire libre.<br /></li><li> -R%CUT%: Al indicar el nombre de la grilla previamente recortada, GMT obtiene la información de la (sub-)región y resolución de la grilla-</li><li>-G%CUT%: Nombre de la grilla creada.</li></ul><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdsample %FAA% -G%CUT% -R%CUT2%</span></span><span style="font-family: courier;"><span style="background-color: #cccccc;"> </span><br /></span></div><div><p><b>3. Imagen:</b> con grdimage se combina la grilla de anomalías de aire (%CUT) con la grilla para el efecto de sombreado (%SHADOW%).<span style="font-family: courier;"><br /></span></p><p><span style="background-color: #cccccc;"><span style="font-family: courier;">gmt grdimage -R -J -O -K %CUT% -C%color% >> %OUT% -I%SHADOW%</span></span></p><p><span style="font-family: courier;"></span></p><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-EKj3P24AXCY/X26xLaLO4vI/AAAAAAAAX3U/NLdurYNG2FYbvc_Lw_Y0dSXb8xjXJJt8gCLcBGAsYHQ/s2048/EJ7.7_Mapa_DEM%252BGrav_AireLibre_2.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1256" data-original-width="2048" height="187" src="https://1.bp.blogspot.com/-EKj3P24AXCY/X26xLaLO4vI/AAAAAAAAX3U/NLdurYNG2FYbvc_Lw_Y0dSXb8xjXJJt8gCLcBGAsYHQ/w307-h187/EJ7.7_Mapa_DEM%252BGrav_AireLibre_2.png" width="307" /></a><a href="https://1.bp.blogspot.com/-W2exQzaJGaM/X3FSPYLezSI/AAAAAAAAX44/hlJEyqknNfIb_a99rDUtuw-8ueMLEivzwCLcBGAsYHQ/s2048/EJ7.6_Mapa_DEM%252BGrav_AireLibre_1.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1256" data-original-width="2048" height="188" src="https://1.bp.blogspot.com/-W2exQzaJGaM/X3FSPYLezSI/AAAAAAAAX44/hlJEyqknNfIb_a99rDUtuw-8ueMLEivzwCLcBGAsYHQ/w310-h188/EJ7.6_Mapa_DEM%252BGrav_AireLibre_1.png" width="310" /></a></div><br /><br /><p></p></div></div>Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0tag:blogger.com,1999:blog-1509523996050466325.post-17270181435520453452020-06-15T00:00:00.036-03:002020-07-24T15:07:05.648-03:0022.1 Animación Tierra Girando<div>Animaciones tierra girando.<br /></div><div class="separator" style="clear: both; text-align: center;"><br /></div><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-TQceJXreLfs/XuO6iQxJAGI/AAAAAAAAXLY/WMdUqMsrqNQQEwUgHXv8GsEnok0Rbt5AwCK4BGAsYHg/s1500/TierraGirando.gif" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1500" data-original-width="1500" height="310" src="https://1.bp.blogspot.com/-TQceJXreLfs/XuO6iQxJAGI/AAAAAAAAXLY/WMdUqMsrqNQQEwUgHXv8GsEnok0Rbt5AwCK4BGAsYHg/w320-h320/TierraGirando.gif" width="310" /></a><a href="https://1.bp.blogspot.com/-E0CMGM3XGjw/XuZaTysaOEI/AAAAAAAAXMU/ifCZxHXMRKQl7k14m1lno7TnRyB0RmQBQCK4BGAsYHg/s600/TierraGirando2.gif" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="600" data-original-width="600" height="310" src="https://1.bp.blogspot.com/-E0CMGM3XGjw/XuZaTysaOEI/AAAAAAAAXMU/ifCZxHXMRKQl7k14m1lno7TnRyB0RmQBQCK4BGAsYHg/w320-h320/TierraGirando2.gif" width="310" /></a><a href="https://1.bp.blogspot.com/-TQceJXreLfs/XuO6iQxJAGI/AAAAAAAAXLY/WMdUqMsrqNQQEwUgHXv8GsEnok0Rbt5AwCK4BGAsYHg/s1500/TierraGirando.gif" style="margin-left: 1em; margin-right: 1em;"></a><a href="https://1.bp.blogspot.com/-TQceJXreLfs/XuO6iQxJAGI/AAAAAAAAXLY/WMdUqMsrqNQQEwUgHXv8GsEnok0Rbt5AwCK4BGAsYHg/s1500/TierraGirando.gif" style="margin-left: 1em; margin-right: 1em;"></a></div><div></div><span><a name='more'></a></span><div><div>Para crear una animación en GMT se requieren 2 scripts. El primero o principal (<i>mainscript</i>) se usa para crear la figura
base. El segundo script (animacion) creará distintas figuras (o fotogramas/<i>frames</i>) a
partir del principal con ligeras modificaciones que finalmente combinará creando la animación.</div><div><br /></div><div><u>Nota</u>: La creación de la animación puede llevar mucho tiempo. Se sugiere empezar con animaciones de pocas figuras y baja calidad para ver que se obtiene el resultado deseado.<br /></div><div><br /></div></div><div>Requerimientos:</div><ol style="text-align: left;"><li>GMT 6: El comando <a href="https://docs.generic-mapping-tools.org/latest/movie.html#a">movie</a> no esta en versiones anteriores.</li></ol><div>Opcional:<br /></div><ol style="text-align: left;"><li>GraphicsMagick: Necesario <a href="http://www.graphicsmagick.org/">instalarlo</a> para crear GIF.<br /></li></ol><div><br /></div><div>Scripts:</div><div><ol style="text-align: left;"><li><a href="https://www.dropbox.com/s/h0qa7crtss2wc1h/EJ1.2_Modern_Hemisferico.bat?dl=0">EJ1.2_Modern_Hemisferico</a></li><li><a href="https://www.dropbox.com/s/bj8ejrjia3zf9uu/EJ22.1_Animacion.bat?dl=0">EJ22.1_Animacion</a></li></ol></div><div><br /></div><div><u>Nota</u>: Los <i>scripts</i> deben estar escritos con <a href="https://docs.generic-mapping-tools.org/6.0/cookbook/introduction.html#modern-and-classic-mode" target="_blank">Sintaxis Moderna</a>.<br /></div><div><br /></div><div><b>1. Script principal</b>: Se usa el <i>script</i> del mapa hemisférico (ver 1.2 para más detalles) como base para crear todas las figuras que se utilizarán para la animación. Para adaptarlo, sólo hay que reemplazar el valor que se desea que varía por la variable %MOVIE_COL0%. En este caso, para que lograr el efecto de la tierra girando hay que cambiar el parámetro de la proyección que indica en qué meridiano se centra la proyección (que en el original era -65).<span style="background-color: #eeeeee;"></span></div><div><br /><span style="background-color: #eeeeee;"></span><span style="background-color: #eeeeee;"></span></div><div style="text-align: left;"><pre><span style="background-color: #eeeeee;"><span style="font-family: "courier";"><strike>SET PROJ=G-65/-30/90/15c</strike><br />SET PROJ=G<b>%MOVIE_COL0%</b>/-30/90/14.9c</span><br /><br /></span></pre></div><div><br /></div><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-Ko5RyEI_RYQ/Xs7bWlUYl2I/AAAAAAAAXIs/OQAPKaU28p8f_UajSYvwBhNDZt4j_DRcACK4BGAsYHg/Globe1.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1755" data-original-width="1755" height="320" src="https://1.bp.blogspot.com/-Ko5RyEI_RYQ/Xs7bWlUYl2I/AAAAAAAAXIs/OQAPKaU28p8f_UajSYvwBhNDZt4j_DRcACK4BGAsYHg/s320/Globe1.png" /></a></div><div><br /></div><div><div class="separator" style="clear: both; text-align: center;"><br /></div></div><div><br /></div><div><b>2. Animación:</b> tiene dos partes fundamentales: Lista de Valores y Movie.<br /></div><div><br /></div><div><b>2.1. Lista de Valores</b>: En esta parte se utiliza el comando <a href="https://docs.generic-mapping-tools.org/latest/gmtmath.html">gmtmath</a> para generar una lista de valores que se le asignará a la variable %MOVIE_COL0% del script principal (y que se guarda en "temp_col0.txt"). En este caso, el archivo tendrá valores de 0 a 359 con incrementos de a 1 (360 valores en total).</div><div><br /></div><div><u>Nota</u>: También es posible crear este archivo a mano con todos los valores que uno quisiera.<br /></div><div style="text-align: left;"><pre><span style="background-color: #9e9e9e;"><span style="font-family: "courier";"><span style="background-color: #eeeeee;">gmt math -T0/359/1 -o1 T = "temp_col0.txt"</span></span><br /><br /><br /></span><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-EgLCF2zq-rI/XuKC-cpMu4I/AAAAAAAAXK4/jSzCq6gNjn8kK_ZBR-eIr6qqbl5Vv8IewCK4BGAsYHg/s367/TierraGirando_ScreenCapture_1.PNG" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="367" data-original-width="203" height="320" src="https://1.bp.blogspot.com/-EgLCF2zq-rI/XuKC-cpMu4I/AAAAAAAAXK4/jSzCq6gNjn8kK_ZBR-eIr6qqbl5Vv8IewCK4BGAsYHg/s320/TierraGirando_ScreenCapture_1.PNG" /></a></div></pre></div><div><br /></div><div><b>2.2: Movie</b>: En esta parte se utiliza <a href="https://docs.generic-mapping-tools.org/latest/movie.html#a">movie</a> para crear todos las figuras a partir del script principal y, además, se indica las características que tendrá la animación. En este ejemplo se utilizan los siguientes opciones:</div><div><ul style="text-align: left;"><li>"EJ1.2_Modern_Hemisferico.bat": <span>Nombre del script principal.</span></li><li><span>-N: Nombre del archivo de salida (y de la subcarpeta donde se crearan todos los archivos necesarios para el video).</span></li><li><span>-A: Permite opcionalmente crear un gif.<br /></span></li><li><span>-F: Permite opcionalmente crear un video en formato mp4.</span></li><li><span>-C: Define el tamaño del lienzo que tendrá el video. En este caso, 15 x 15 cm y 100 dpu (dots-per-unit).</span></li><li><span>-T</span><span><span>"temp_col0.txt"</span>: Nombre del archivo con la lista de valores. <br /></span></li></ul></div><div></div><div style="text-align: left;"><pre style="text-align: left;"><span style="background-color: #9e9e9e;"><span style="font-family: "courier";"><span style="background-color: #eeeeee;">gmt movie "EJ1.2_Modern_Hemisferico.bat" -N"TierraGirando" -T"temp_col0.txt" -C15cx15cx100 -Fmp4 -D12 -Agif -V3</span></span><br /><br /><b><br /></b><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-IB3gT6X7oO8/XuZfxeORtjI/AAAAAAAAXMw/CjDfTTFs9O8isfD2sD3QwaCwR-K7Q-RtwCK4BGAsYHg/s1500/TierraGirando.gif" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1500" data-original-width="1500" height="320" src="https://1.bp.blogspot.com/-IB3gT6X7oO8/XuZfxeORtjI/AAAAAAAAXMw/CjDfTTFs9O8isfD2sD3QwaCwR-K7Q-RtwCK4BGAsYHg/s320/TierraGirando.gif" /></a><div class="separator" style="clear: both; text-align: center;"><iframe allowfullscreen="" class="BLOG_video_class" height="266" src="https://www.youtube.com/embed/b96nDjucJzQ" width="320" youtube-src-id="b96nDjucJzQ"></iframe></div><div class="separator" style="clear: both; text-align: center;"><br /></div><div class="separator" style="clear: both; text-align: center;"><br /></div><div class="separator" style="clear: both; text-align: center;"><br /></div></div></span></pre><p><b>2.3. Reusar figuras: </b>El argumento -V de movie muestra el comando de <a href="https://ffmpeg.org/">ffmpeg</a> que utilizó GMT para crear la animación. Esto es útil para hacer una nueva animación a partir de las mismas figuras pero con distintos parámetros (como framerate).</p><p></p><pre style="text-align: left;"><pre><span style="font-family: "courier";"><span style="background-color: #9e9e9e;"><span style="background-color: #eeeeee;">ffmpeg -loglevel warning -f image2 <b>-framerate 12</b> -y -i "TierraGirando/TierraGirando_%%03d.png"<br /> -vcodec libx264 -pix_fmt yuv420p TierraGirando.mp4</span></span></span></pre><pre><span style="background-color: #9e9e9e;"></span></pre><span style="background-color: #9e9e9e;"><br /></span><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-U-tt9YBbVIg/XuUM3Vz5OuI/AAAAAAAAXL0/4lV4-qki9d84LuvsOsHQcl_mnong7j1xgCK4BGAsYHg/s653/TierraGirando_ScreenCapture_2.PNG" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="230" data-original-width="653" height="224" src="https://1.bp.blogspot.com/-U-tt9YBbVIg/XuUM3Vz5OuI/AAAAAAAAXL0/4lV4-qki9d84LuvsOsHQcl_mnong7j1xgCK4BGAsYHg/w640-h224/TierraGirando_ScreenCapture_2.PNG" width="640" /></a></div></pre></div><div><br /></div><div><b>2.4. Apagar/Hibernar PC</b>: El script Animacion puede tardar mucho tiempo en realizar todo el proceso (el de este ejemplo estuvo 3:27 hs) por lo que es conveniente dejar la pc trabajando sóla. Los siguientes comandos ubicados al final del script permiten apagar (o hibernar) la PC (en Windows) inmediatamente una vez finalizada el proceso.<br /></div><div></div><div></div><div style="text-align: left;"><pre><span style="background-color: #eeeeee;"><span style="font-family: "courier";">shutdown /s<br />shutdonw /h<br /></span></span><br /></pre><p style="text-align: left;"><b>3. Animación 2</b>: En cambio para crear una animación con la tierra girando de polo sur a polo norte hay que utilizar los siguientes comandos:<br /></p><b>1. Script Principal:<span style="background-color: #eeeeee;"><span style="font-family: "courier";"></span></span></b><span style="background-color: #eeeeee;"><span style="font-family: "courier";"></span><span style="font-family: "courier";"></span></span><span style="background-color: #eeeeee;"><span style="font-family: "courier";"></span></span><br /></div><div style="text-align: left;"><span style="background-color: #eeeeee;"><span style="font-family: "courier";"><br /></span></span></div><div style="text-align: left;"><span style="background-color: #eeeeee;"><span style="font-family: "courier";">SET PROJ=G-65/%MOVIE_COL0%/90/14.9c</span></span><b><br /></b></div><div style="text-align: left;"><b><br /></b></div><div style="text-align: left;"><b>2. Lista de Valores:</b><pre style="text-align: left;"><span style="background-color: #eeeeee;"><span style="font-family: "courier";">gmt math -T-90/269/1 -o1 T = "temp_col0.txt"<br /><br /><br /></span><div class="separator" style="clear: both; text-align: center;"><a href="https://1.bp.blogspot.com/-CsO3wTX1dVA/XuZo5F9IQBI/AAAAAAAAXNY/eHxV5UX2VEM-sWuDvI4o8Y6VfgAr3hZagCK4BGAsYHg/s1500/TierraGirando2.gif" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1500" data-original-width="1500" height="320" src="https://1.bp.blogspot.com/-CsO3wTX1dVA/XuZo5F9IQBI/AAAAAAAAXNY/eHxV5UX2VEM-sWuDvI4o8Y6VfgAr3hZagCK4BGAsYHg/s320/TierraGirando2.gif" /></a><div class="separator" style="clear: both; text-align: center;"></div><div class="separator" style="clear: both; text-align: center;"><pre style="text-align: left;"><span style="background-color: #9e9e9e;"><div class="separator" style="clear: both; text-align: center;"><iframe allowfullscreen="" class="BLOG_video_class" height="367" src="https://www.youtube.com/embed/rqtN-QTkYrw" width="414" youtube-src-id="rqtN-QTkYrw"></iframe></div></span></pre></div></div><br /></span></pre></div><div style="text-align: left;"><pre></pre></div>Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0Buenos Aires, CABA, Argentina-34.6036844 -58.3815591-62.913918236178844 -93.5378091 -6.2934505638211533 -23.225309099999997tag:blogger.com,1999:blog-1509523996050466325.post-31599225895016619092020-04-13T00:00:00.000-03:002020-05-05T14:26:18.742-03:009.1. Superficie de tendencia. Anomalia regional y local<div dir="ltr" style="text-align: left;" trbidi="on">
<div class="separator" style="clear: both; text-align: center;">
</div>
Filtrar una grilla con el objetivo de separar el aporte regional de la local. Para eso la grilla se ajusta a una superficie de tendencia definida por un polinomio de bajo grado.<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-kPaILGJOghA/XpDMqoIelxI/AAAAAAAAW6o/PCf4QkUf-WMwYHO6bY6Z85ThIUgz3E2QwCLcBGAsYHQ/s1600/EJ8.6_Mapa_Grav_AireLibre_T10.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="810" data-original-width="1600" height="320" src="https://1.bp.blogspot.com/-kPaILGJOghA/XpDMqoIelxI/AAAAAAAAW6o/PCf4QkUf-WMwYHO6bY6Z85ThIUgz3E2QwCLcBGAsYHQ/s640/EJ8.6_Mapa_Grav_AireLibre_T10.png" width="640" /></a></div>
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-2ghEA0B3nEc/XpDMvF_oXTI/AAAAAAAAW6s/GYABFe7XgJ8KsjFHj1b5LHM4ofjAKLKRwCLcBGAsYHQ/s1600/EJ8.6_Mapa_Grav_AireLibre_D10.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="810" data-original-width="1600" height="320" src="https://1.bp.blogspot.com/-2ghEA0B3nEc/XpDMvF_oXTI/AAAAAAAAW6s/GYABFe7XgJ8KsjFHj1b5LHM4ofjAKLKRwCLcBGAsYHQ/s640/EJ8.6_Mapa_Grav_AireLibre_D10.png" width="640" /></a></div>
<div style="text-align: center;">
</div>
<a name='more'></a><b>Introducción:</b> Las anomalías gravitatorias usualmente provienen de distintas fuentes. Según sus dimensiones se pueden definir anomalías regionales (de gran longitud de onda) y anomalías locales (o residuales, de corta longitud de onda). Una forma de separarlas es aproximar las anomalías a una superficie suave descripta por un polinomio de grado bajo. En este ejemplo se muestra como separar (y graficar) las anomalías regionales y locales. <br />
<br />
<a href="https://www.dropbox.com/s/72tpm42xzgt2qim/EJ9.1_Filtrado_Superficie-Tendencia.bat?dl=0">Script</a><br />
<br />
<div dir="ltr">
<b>0. Datos</b>:<br />
<ul style="text-align: left;">
<li>Anomalías de Aire Libre: <a href="ftp://topex.ucsd.edu/pub/global_grav_1min/grav.img.23.1">grav.img.29.1 </a>(en formato img). Ver <a href="http://mapasgeologicos.blogspot.com.ar/p/blog-page_70.html">Base de Datos</a>.</li>
</ul>
</div>
<b>1A. Mapa Base</b>: Se crea con los siguientes comandos (ver <a href="http://mapasgeologicos.blogspot.com/2016/02/5-mapa-topografico.html">7.1</a> para más detalles):<br />
<br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt img2grd "E:\Facultad\Datos_Geofisicos\Gravimetria\Sandwell-Smith\Datos\grav_29.1.img" -R -G%CUT% -T1 -I1 -E -S0.1 </span></span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt grdgradient %CUT% -A0/270 -G%SHADOW% -Ne0.5</span></span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt grdimage -R -J -O -K %CUT% -C%color% -I%SHADOW% >> %OUT%</span></span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt pscoast -R -J -O -K -Df -W1/thinnest >> %OUT%</span></span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt psbasemap -R -J -O -K -Bxaf -Byaf >> %OUT%</span></span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt psscale -O -K -Dx15.3/0.3+w8.1/0.618c+e -C%color% -Ba+l"Anomal\355as Aire Libre (mGal)" -I >> %OUT%</span></span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt psbasemap -R -J -O -K -Ln0.08/0.075+c-52:00+w100k+f+l >> %OUT%</span></span></span><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-OqGKHXwf7bc/Xo3dJvItK7I/AAAAAAAAW5Y/NdZOseOi_H4PLOLeqdewGNRznjDkuznlwCLcBGAsYHQ/s1600/EJ8.6_Mapa_Grav_AireLibre_0.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="810" data-original-width="1600" height="200" src="https://1.bp.blogspot.com/-OqGKHXwf7bc/Xo3dJvItK7I/AAAAAAAAW5Y/NdZOseOi_H4PLOLeqdewGNRznjDkuznlwCLcBGAsYHQ/s400/EJ8.6_Mapa_Grav_AireLibre_0.png" width="400" /></a></div>
<br />
<b>1A. Mapa Base. Curvas de Nivel</b>: Agrego curvas de nivel (ver <a href="https://mapasgeologicos.blogspot.com/2016/03/53-mapa-topografico-iii-curvas-de-nivel.html">5.3</a> para más detalles).<br />
<br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt grdcontour -R -J -O -K %CUT% >> %OUT% -C20 -Q100 -W0.25,black,-</span></span><br />
<br />
<div style="text-align: center;">
<a href="https://1.bp.blogspot.com/-bm4e-CgoDg0/XozZ8ynDUyI/AAAAAAAAW5M/xRvnabWG7ww30wQj9v7JEQjhvXY8SCrEACLcBGAsYHQ/s1600/EJ8.6_Mapa_Grav_AireLibre_T0.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="200" src="https://1.bp.blogspot.com/-bm4e-CgoDg0/XozZ8ynDUyI/AAAAAAAAW5M/xRvnabWG7ww30wQj9v7JEQjhvXY8SCrEACLcBGAsYHQ/s400/EJ8.6_Mapa_Grav_AireLibre_T0.png" width="400" /></a> </div>
<br />
<b>2. Grilla de tendencia:</b> con <a href="https://docs.generic-mapping-tools.org/latest/grdtrend.html">grdtrend</a> se lee la grilla %CUT% y se ajusta al polinomio de bajo grado:<br />
<br />
<div style="text-align: center;">
<b>m<sub>1</sub> + m<sub>2</sub>x + m<sub>3</sub>y + m<sub>4</sub>xy + m<sub>5</sub>x<sup>2</sup> + m<sub>6</sub>y<sup>2</sup> + m<sub>7</sub>x<sup>3</sup> + m<sub>8</sub>x<sup>2</sup>y + m<sub>9</sub>xy<sup>2</sup> + m<sub>10</sub>y<sup>3</sup></b></div>
<ul style="text-align: left;"></ul>
<br />
<b>2A. Superficie de Tendencia: </b>La superficie se según las opciones:<br />
<ul style="text-align: left;">
<li>-N: grado del polinomio (m<sub>x</sub>) entre 1 y 10. </li>
<li>-T: Calcula la superficie de tendencia.</li>
</ul>
<u>Nota</u>: Es necesario conocer la geología de la zona para determinar que superficie (y por lo tanto grado del polinomio) se utiliza para describir la anomalía regional.<br />
<ul style="text-align: left;"></ul>
<br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt grdtrend %CUT% -T%CUT% -N10</span></span> <br />
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-ED428ptMUm4/XpDOAX7r5XI/AAAAAAAAW68/AS-GSWfPiV4hF80OcKrG280DCm2RZfzKACLcBGAsYHQ/s1600/EJ8.6_Mapa_Grav_AireLibre_T01.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="810" data-original-width="1600" height="200" src="https://1.bp.blogspot.com/-ED428ptMUm4/XpDOAX7r5XI/AAAAAAAAW68/AS-GSWfPiV4hF80OcKrG280DCm2RZfzKACLcBGAsYHQ/s400/EJ8.6_Mapa_Grav_AireLibre_T01.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">N1: Plano horizontal.</td></tr>
</tbody></table>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-cwNHqPmTJ_Q/XpDOGR4TFnI/AAAAAAAAW7A/AQFVsTbwGaIY7epND6ykp_TXOgPgZBi_QCLcBGAsYHQ/s1600/EJ8.6_Mapa_Grav_AireLibre_T03.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="810" data-original-width="1600" height="200" src="https://1.bp.blogspot.com/-cwNHqPmTJ_Q/XpDOGR4TFnI/AAAAAAAAW7A/AQFVsTbwGaIY7epND6ykp_TXOgPgZBi_QCLcBGAsYHQ/s400/EJ8.6_Mapa_Grav_AireLibre_T03.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">N3: Plano inclinado (bilinear).</td></tr>
</tbody></table>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"></td></tr>
<tr><td class="tr-caption" style="text-align: center;"></td></tr>
</tbody></table>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-qUVlZ-c23Mk/XpDOORRwULI/AAAAAAAAW7E/rnnTcrLZNq8qh7A_z5k2cbECQekTRXmgwCLcBGAsYHQ/s1600/EJ8.6_Mapa_Grav_AireLibre_T06.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="810" data-original-width="1600" height="200" src="https://1.bp.blogspot.com/-qUVlZ-c23Mk/XpDOORRwULI/AAAAAAAAW7E/rnnTcrLZNq8qh7A_z5k2cbECQekTRXmgwCLcBGAsYHQ/s400/EJ8.6_Mapa_Grav_AireLibre_T06.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">N6: Superficie cuadrática.</td></tr>
</tbody></table>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-HxfCGsyhl2M/XrGLCn1fGRI/AAAAAAAAXEU/KaJdd2BwcA0o8ZNvM6qdvqUKQEfsOmfVwCLcBGAsYHQ/s1600/EJ9.1_Filtrado_Superficie-Tendencia_5_T10.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="810" data-original-width="1600" height="201" src="https://1.bp.blogspot.com/-HxfCGsyhl2M/XrGLCn1fGRI/AAAAAAAAXEU/KaJdd2BwcA0o8ZNvM6qdvqUKQEfsOmfVwCLcBGAsYHQ/s400/EJ9.1_Filtrado_Superficie-Tendencia_5_T10.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">N10: Superficie bicúbica.</td></tr>
</tbody></table>
<br />
<b>2A. Superficie Residual: </b>En cambio con -D se calcula la diferencia de la grilla original con respecto a la grilla de tendencia.<br />
<ul style="text-align: left;">
<li>-D: Calcula la diferencia entre la superfice original y la de tendencia. </li>
</ul>
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt grdtrend %CUT% -D%CUT% -N10</span></span><br />
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-75mFsnaGWGE/XpDXYPEZ5RI/AAAAAAAAW7Y/PGy-y42NcXQ--J_JvKUa0xPPijXhI_AeACLcBGAsYHQ/s1600/EJ8.6_Mapa_Grav_AireLibre_D6.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="810" data-original-width="1600" height="200" src="https://1.bp.blogspot.com/-75mFsnaGWGE/XpDXYPEZ5RI/AAAAAAAAW7Y/PGy-y42NcXQ--J_JvKUa0xPPijXhI_AeACLcBGAsYHQ/s400/EJ8.6_Mapa_Grav_AireLibre_D6.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">N6: Superficie cuadrática.</td></tr>
</tbody></table>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-BcHMMadQ4fM/XpDXZWQbVhI/AAAAAAAAW7c/IyML7Vcd8PoYAckhGVwsmtjTkjFT5Sh6gCLcBGAsYHQ/s1600/EJ8.6_Mapa_Grav_AireLibre_D10.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="810" data-original-width="1600" height="201" src="https://1.bp.blogspot.com/-BcHMMadQ4fM/XpDXZWQbVhI/AAAAAAAAW7c/IyML7Vcd8PoYAckhGVwsmtjTkjFT5Sh6gCLcBGAsYHQ/s400/EJ8.6_Mapa_Grav_AireLibre_D10.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">N10: Superficie bicúbica. </td></tr>
</tbody></table>
<blockquote>
<div>
</div>
</blockquote>
</div>
Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0tag:blogger.com,1999:blog-1509523996050466325.post-37987488102189670572019-12-11T00:00:00.000-03:002020-02-03T23:28:43.984-03:0021.2. Zona Wadati-Benioff II.<div dir="ltr" style="text-align: left;" trbidi="on">
Graficar sismos y mecanismos focales en un perfil de distancia versus profundidad en cualquier dirección.<br />
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div style="text-align: center;">
<a href="https://1.bp.blogspot.com/-Ucx2atWIS7M/Xehq4juiGKI/AAAAAAAAVxs/uhxM8gugxbA6ahJxcsHrtbt284xS956MACEwYBhgL/s1600/EJ21.2_Mapa_Sismicidad_II.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="840" data-original-width="1600" height="336" src="https://1.bp.blogspot.com/-o2dKnL2vxYY/Xd9AEKu0gzI/AAAAAAAAVt8/ofM3ciFL4CEGeEGLSGvYzJyjo_Xgt0xXwCLcBGAsYHQ/s640/EJ21.2_Perfil_Wadati-Benioff_II.png" width="640" /> <img border="0" data-original-height="1002" data-original-width="1600" height="250" src="https://1.bp.blogspot.com/-Ucx2atWIS7M/Xehq4juiGKI/AAAAAAAAVxs/uhxM8gugxbA6ahJxcsHrtbt284xS956MACEwYBhgL/s400/EJ21.2_Mapa_Sismicidad_II.png" width="400" /></a> </div>
<a name='more'></a><b>Introducción</b>:<br />
Esta ejemplo es una adaptación de la entrada <a href="https://mapasgeologicos.blogspot.com/2019/11/211-zona-wadati-benioff-i-perfil.html">21.1</a>. Se explican las adaptaciones para hacer perfiles siguiendo círculos máximos en cualquier orientación (no sólo E-W o N-S). Aquí sólo se explican las diferencias. Se filtran los eventos (sismos y mecanismos focales) ubicados hasta 100 km de cada lado del perfil.<br />
<br />
<b>Datos</b>:<br />
<ul style="text-align: left;">
<li>DEM: Grilla GMRT 3.6 (ver <a href="https://mapasgeologicos.blogspot.com/2017/01/anexo-4-descargar-grilla-gmrt.html">anexo 4</a>).</li>
<li>Sismos USGS: Archivos query_*.csv (<a href="https://www.dropbox.com/s/0gmbmhms6umfdt7/USGS.zip?dl=0">USGS.zip</a>; ver ANSS en <a href="https://mapasgeologicos.blogspot.com/p/blog-page_70.html">base de datos</a>).</li>
<li>Mecanismos Focales: Archivos CMT_*.txt (comprimidos en <a href="https://www.dropbox.com/s/z427zknrdsd1g7p/CMT_1976-2019.zip?dl=0">CMT_1976-2019.zip</a>; ver Global CMT en <a href="https://mapasgeologicos.blogspot.com/p/blog-page_70.html">base de datos</a>).</li>
</ul>
<br />
<a href="https://www.dropbox.com/s/1k0ye11bg6wwi9a/EJ21.2_Mapa_Sismicidad_II.bat?dl=0">Script Mapa</a><br />
<a href="https://www.dropbox.com/s/k9unli530aepqyb/EJ21.2_Perfil_Wadati-Benioff_II.bat?dl=0">Script Perfil</a><br />
<br />
<b>0. Mapa Base</b>: con los sismos y mecanismos focales ploteados y perfil dibujado (ver <a href="https://mapasgeologicos.blogspot.com/2016/02/4-mapa-tectonico.html">4.1</a> y <a href="https://mapasgeologicos.blogspot.com/2017/10/171-mapa-de-calor-de-sismicidad.html">17.1</a> para detalles de los comandos).<br />
<br />
<span style="background-color: #999999;"><span style="font-family: "courier new" , "courier" , monospace;">echo -74 -29 A > "temp_perfil"<br />echo -64 -33 A' >> "temp_perfil"</span></span><br />
<b><br /></b><a href="https://www.blogger.com/null" name="more"></a> <br />
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-bbtWjO1z-kA/Xehq4BkogZI/AAAAAAAAVxo/8dj0l-x66ZQoBsg4zhkb4qRbDea8uEaIQCLcBGAsYHQ/s1600/EJ21.2_Mapa_Sismicidad_II_0.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1002" data-original-width="1600" height="250" src="https://1.bp.blogspot.com/-bbtWjO1z-kA/Xehq4BkogZI/AAAAAAAAVxo/8dj0l-x66ZQoBsg4zhkb4qRbDea8uEaIQCLcBGAsYHQ/s400/EJ21.2_Mapa_Sismicidad_II_0.png" width="400" /></a></div>
<b><br /></b><a href="https://www.blogger.com/null" name="more"></a> <br />
<b>1. Mapa Base 1</b>: Con <a href="http://gmt.soest.hawaii.edu/doc/latest/gmtselect.html">select</a> se filtran los eventos (sismos y mecanismos focales) que se localizan hasta 100 km a cada lado del perfil A-A' con los siguientes argumentos:<br />
<ul style="text-align: left;">
<li>-L"temp_perfil": archivo con los puntos que definen el perfil.</li>
<li>+d100k: Distancia de 100 km perpendicular al perfil.</li>
<li>+p: Solo selecciona los puntos hacia los lados del perfil, excluyendo los que están más allá de los extremos.</li>
<li>-fg: Indica que los datos son coordenadas geográficas (y no cartesianos).</li>
</ul>
<br />
<span style="background-color: #999999;"><span style="font-family: "courier new" , "courier" , monospace;">gmt select "query_*.csv" -h1 -i2,1,3 > "temp_sismos" -L"temp_perfil"+d</span><span style="font-family: "courier new" , "courier" , monospace;"><span style="font-family: "courier new" , "courier" , monospace;">100k</span>+p -fg<br />gmt select "CMT_*.txt" > "temp_CMT" -L"temp_perfil"+d100k+p -fg </span></span><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-Ucx2atWIS7M/Xehq4juiGKI/AAAAAAAAVxs/uhxM8gugxbA6ahJxcsHrtbt284xS956MACEwYBhgL/s1600/EJ21.2_Mapa_Sismicidad_II.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1002" data-original-width="1600" height="250" src="https://1.bp.blogspot.com/-Ucx2atWIS7M/Xehq4juiGKI/AAAAAAAAVxs/uhxM8gugxbA6ahJxcsHrtbt284xS956MACEwYBhgL/s400/EJ21.2_Mapa_Sismicidad_II.png" width="400" /></a></div>
<br />
<div style="text-align: center;">
</div>
<br />
<b>2. Perfil topográfico y en profundidad</b>: Se utilizan los comandos explicados en <a href="https://mapasgeologicos.blogspot.com/2016/05/111-perfil-topografico.html">11.1</a><a href="https://www.blogger.com/null">. Perfil Topográfico</a>
para crear el perfil, extraer los datos del dem, calcular cuanto mide
(en km), graficar el perfil con la exageración vertical y su
orientación.<br />
<br />
<br />
<b>3. Perfil en profundidad (inferior)</b>: Se adjunta un perfil en profundidad hasta los 300 km de profundidad. <br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;"><br /></span></span>
<span style="background-color: #999999;"><span style="font-family: "courier new" , "courier" , monospace;">SET DepthMin=0</span></span><br />
<span style="background-color: #999999;"><span style="font-family: "courier new" , "courier" , monospace;">SET DepthMax=300 </span></span><br />
<span style="background-color: #999999;"><span style="font-family: "courier new" , "courier" , monospace;"><br /></span></span>
<span style="background-color: #999999;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psbasemap -R -J -O -K -Bxaf+l"Longitud (\260)" >> %OUT% --MAP_FRAME_AXES=Sn</span></span><br />
<span style="background-color: #999999;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psbasemap -R -J -O -K -Byaf+l"Profundidad (km)" >> %OUT% --FONT_ANNOT_PRIMARY=8,Helvetica,blue --MAP_FRAME_AXES=wE </span></span><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-fpVKHm5lGGY/XehtUSBtrQI/AAAAAAAAVx8/oHxOW9CfkcYpuP33I-js1ZlIlTeW1J67ACLcBGAsYHQ/s1600/EJ21.2_Perfil_Wadati-Benioff_II_0.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="844" data-original-width="1600" height="210" src="https://1.bp.blogspot.com/-fpVKHm5lGGY/XehtUSBtrQI/AAAAAAAAVx8/oHxOW9CfkcYpuP33I-js1ZlIlTeW1J67ACLcBGAsYHQ/s400/EJ21.2_Perfil_Wadati-Benioff_II_0.png" width="400" /></a></div>
<br />
<br />
<b>4. Filtrar y Proyectar sismos al círculo máximo: </b>con <a href="http://gmt.soest.hawaii.edu/doc/latest/project.html">project</a> se filtran los sismos y proyectan al círculo máximo (perfil) con los siguientes argumentos:<br />
<br />
<ul style="text-align: left;">
<li>-C: Inicio del perfil (definido por las variables Long1 y Lat1).</li>
<li>-E: Fin del perfil (Long2 y Lat2).</li>
<li>-W: Filtra eventos según la distancia perpendicular al perfil (q hasta ± 100 km en este caso).</li>
<li>-Lw: Se excluyen los eventos más alla de los extremos del perfil.</li>
<li>-S: Ordena los eventos del mas cercano al inicio del perfil al más lejano.</li>
<li>-i2,1,3: Columnas con los datos de longitud (2), latitud (1) y profundidad (3) en los archivos "query_*.csv"</li>
<li>-Q: Indica unidades de grado para columnas de Longitud y Latitud.</li>
</ul>
<span style="background-color: #999999;"><span style="font-family: "courier new" , "courier" , monospace;">gmt project "query_*.csv" -C%Long1%/%Lat1% -E%Long2%/%Lat2% -W-%Dist_Perfil%/%Dist_Perfil% -Lw -S -Q -h -i2,1,3 > "temp_sismos_project.txt"</span></span><br />
<br />
Los datos se guardan en el archivo "temp_sismos_project.txt" que contiene 7 columnas. Las 3 primeras corresponden a Longitud, Latitud y Profundidad de cada sismo según el archivo de entrada. Para cada evento se agregan columnas con 4 datos: <br />
<ul style="text-align: left;">
<li>p: Distancia al inicio del perfil (definido por -C) medido sobre el perfil.</li>
<li>q: Distancia perpendicular al perfil.</li>
<li>Long Proj: Longitud del punto más cercano sobre el círculo máximo.</li>
<li>Lat Proj: Latitud del punto más cercano sobre el círculo máximo.</li>
</ul>
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-1UDzWo4G6KM/Xe6Yt8i5SPI/AAAAAAAAVyc/VUVvHBSuvOg5IXmrrrIQage6IIcEnfXVwCLcBGAsYHQ/s1600/Captura_Seis_Project.PNG" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="405" data-original-width="527" height="306" src="https://1.bp.blogspot.com/-1UDzWo4G6KM/Xe6Yt8i5SPI/AAAAAAAAVyc/VUVvHBSuvOg5IXmrrrIQage6IIcEnfXVwCLcBGAsYHQ/s400/Captura_Seis_Project.PNG" width="400" /></a></div>
<br />
<br />
<b>6. Graficar sismos en distancia del perfil versus profundidad</b>: con psxy se grafican los datos del archivo "temp_sismos_project". Para graficarlos se seleccionan las columnas 3 (p) y profundidad (2).<br />
<span style="background-color: #999999;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psxy -R -J -O -K "temp_sismos_project" >> %OUT% -C"temp_seis" -Sc0.05c -i3,2,2 </span></span><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-A7yxX9YaO-s/Xe-7QxZzgII/AAAAAAAAVzE/UJlradZ5_5kdcMsVnvABkVQ3U1vANajggCLcBGAsYHQ/s1600/EJ21.2_Perfil_Wadati-Benioff_II_1.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="844" data-original-width="1600" height="336" src="https://1.bp.blogspot.com/-A7yxX9YaO-s/Xe-7QxZzgII/AAAAAAAAVzE/UJlradZ5_5kdcMsVnvABkVQ3U1vANajggCLcBGAsYHQ/s640/EJ21.2_Perfil_Wadati-Benioff_II_1.png" width="640" /></a></div>
<br />
<b>7. Filtrar Mecanismos Focales:</b> Con los mismos comandos del punto 1, se filtran
los datos de los eventos a 100 km de cada lado del perfil.<br />
<br />
<span style="background-color: #999999;"><span style="font-family: "courier new" , "courier" , monospace;">gmt select "Mecanismos_Focales\CMT_*.txt" > "temp_CMT" -L"temp_perfil"+d%Dist_Perfil%k+p -fg </span></span><br />
<br />
<b>8. Proyectar mecanismos focales</b>:<a href="https://docs.generic-mapping-tools.org/latest/supplements/seis/pscoupe.html">pscoupe</a> permite graficar los mecanimos focales previamente filtrados (temp_CMT) en un perfil con los siguientes argumentos:<br />
<ul style="text-align: left;">
<li>-Aa: Perfil (longitud y latitud inicial y final), la inclinación del
perfil (90 = vertical), extensión hacia los lados del perfil (10 mil
km, para incluir a todos los mecanismos previamente filtrados) y profundidad mínima y máxima. </li>
<li>-Sd0.15/0: <i>beachballs</i> de un tamaño de 0,15 (para los eventos de magnitud 5
en pulgadas por default). El segundo valor corresponde al tamaño de la
fuente usado para el nombre del evento. Con 0 no se escribe. </li>
<li>-Q: No crea los archivos de información de los mecanismos focales.</li>
<li>-G: Color para los cuadrantes compresivos.</li>
<li>-M: Tamaño fijo.</li>
</ul>
<span style="background-color: #999999;">gmt pscoupe -R -J -O -K "temp_CMT" >> %OUT% -Sd0.15/0 -Gred -M -Aa%Long1%/%Lat1%/%Long2%/%Lat2%/90/10000/%DepthMin%/%DepthMax%f -W.1 -Q </span><br />
<span style="background-color: #999999;"><br /></span>
<span style="background-color: #999999;"></span><br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-2maGaqZwfi0/XfBVrIALu9I/AAAAAAAAVzo/G72Gx809sXAuJlSDGpADLTD9KMk38RvbgCLcBGAsYHQ/s1600/EJ21.2_Perfil_Wadati-Benioff_II_2.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="844" data-original-width="1600" height="336" src="https://1.bp.blogspot.com/-2maGaqZwfi0/XfBVrIALu9I/AAAAAAAAVzo/G72Gx809sXAuJlSDGpADLTD9KMk38RvbgCLcBGAsYHQ/s640/EJ21.2_Perfil_Wadati-Benioff_II_2.png" width="640" /></a></div>
<br />
<ul style="text-align: left;">
</ul>
</div>
Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0Buenos Aires, Argentina-34.640856073580395 -58.346582648706942-36.305057073580393 -60.928369648706941 -32.976655073580396 -55.764795648706944tag:blogger.com,1999:blog-1509523996050466325.post-43417138557214429842019-11-12T00:00:00.000-03:002020-02-03T23:31:05.664-03:0021.1. Zona Wadati-Benioff I. Perfil Topografico + Sismos/Mecanismos Focales.<div dir="ltr" style="text-align: left;" trbidi="on">
<div dir="ltr">
<div class="separator" style="clear: both; text-align: center;">
</div>
Graficar sismos y mecanismos focales en un perfil E-W. <br />
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
<img border="0" data-original-height="835" data-original-width="1600" height="332" src="https://1.bp.blogspot.com/-YIGOhgpSfR4/XcGJXrLWrRI/AAAAAAAAVmU/u2r523wV6841kz_3t4vGPNQmuWEOmTuwwCEwYBhgL/s640/EJ21.1_Perfil_Wadati-Benioff_Longitud-Profundidad_3.png" width="640" /> </div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-u8DX96X_I3o/Xbzq2vmx5rI/AAAAAAAAVlQ/RxqnS25VBQ8rfacovtleZUdFHXZAioltgCLcBGAsYHQ/s1600/EJ21.1_Mapa_Sismicidad_1.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1002" data-original-width="1600" height="250" src="https://1.bp.blogspot.com/-u8DX96X_I3o/Xbzq2vmx5rI/AAAAAAAAVlQ/RxqnS25VBQ8rfacovtleZUdFHXZAioltgCLcBGAsYHQ/s400/EJ21.1_Mapa_Sismicidad_1.png" width="400" /></a></div>
<br />
<a name='more'></a><b>Introducción</b>:En las zonas de convergencia activa, se produce la subducción de una placa tectónica por debajo de otra. En la placa que desciende se producen gran cantidad de sismos (y mecanismos focales) que son detectados en la superficie. Al graficarlos en perfiles perpendiculares a la zona de subducción, estos se alinean y conforman una zona de sismicidad denominada zona de <a href="https://es.wikipedia.org/wiki/Zona_de_Wadati-Benioff">Wadati-Benioff</a>. A continuación se muestra un ejemplo de Sudamérica a los 32 ºS.<br />
<br />
0. Datos:</div>
<ul style="text-align: left;">
<li>DEM: Grilla GMRT 3.6 (ver <a href="https://mapasgeologicos.blogspot.com/2017/01/anexo-4-descargar-grilla-gmrt.html">anexo 4</a>).</li>
<li>Sismos USGS: Archivos query_*.csv (<a href="https://www.dropbox.com/s/0gmbmhms6umfdt7/USGS.zip?dl=0">USGS.zip</a>; ver ANSS en <a href="https://mapasgeologicos.blogspot.com/p/blog-page_70.html">base de datos</a>).</li>
<li>Mecanismos Focales: Archivos CMT_*.txt (comprimidos en <a href="https://www.dropbox.com/s/z427zknrdsd1g7p/CMT_1976-2019.zip?dl=0">CMT_1976-2019.zip</a>; ver Global CMT en <a href="https://mapasgeologicos.blogspot.com/p/blog-page_70.html">base de datos</a>).</li>
</ul>
<div dir="ltr">
</div>
<div dir="ltr">
<br />
<a href="https://www.dropbox.com/s/pooqoqsw9sy858g/EJ21.1_Mapa_Sismicidad.bat?dl=0">Script Mapa</a><br />
<a href="https://www.dropbox.com/s/spnf5macf2p53i9/EJ21.1_Perfil_Wadati-Benioff.bat?dl=0">Script Perfil</a><br />
<b><br /></b>
<b><a href="https://www.blogger.com/null" name="more"></a>0. Mapa Base</b>: con los sismos y mecanismos focales ploteados y perfil dibujado (ver <a href="https://mapasgeologicos.blogspot.com/2016/02/4-mapa-tectonico.html">4.1</a> y <a href="https://mapasgeologicos.blogspot.com/2017/10/171-mapa-de-calor-de-sismicidad.html">17.1</a> para detalles de los comandos).</div>
<div dir="ltr">
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-uIU-1YGuD5U/Xbzqv3xCP5I/AAAAAAAAVlM/b15blWxRfiIhNx8PkINbS2jCnpePEtOmgCLcBGAsYHQ/s1600/EJ21.1_Mapa_Sismicidad_0.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1002" data-original-width="1600" height="250" src="https://1.bp.blogspot.com/-uIU-1YGuD5U/Xbzqv3xCP5I/AAAAAAAAVlM/b15blWxRfiIhNx8PkINbS2jCnpePEtOmgCLcBGAsYHQ/s400/EJ21.1_Mapa_Sismicidad_0.png" width="400" /></a></div>
<br />
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<b>1. Mapa base 1</b>: Se filtran que se localizan en la región que se extiende 1 grado hacia el norte y el sur del perfil (entre -33 y -31 ) que van a ser graficados en el perfil perpendicular A-A'. Se utiliza <a href="http://gmt.soest.hawaii.edu/doc/latest/gmtselect.html">select</a> con el argumento -R donde se indica la región de interés.<br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">SET REGION2=-74/-64/-33/-31</span></span><br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt select -R%REGION2% "query_*.csv" -h1 -i2,1,3 > "temp_sismos"<br />gmt select -R%REGION2% "CMT_*.txt" > "temp_CMT"<br /><br />gmt psxy -R%REGION% -J -T -O -K >> %OUT%<br />gmt psxy -R -J -O -K "temp_sismos" >> %OUT% -Sc0.05 -C"temp_seis"<br />gmt psmeca -R -J -O -K "temp_CMT" >> %OUT% -Sd0.12/0 -Gorange </span></span><br />
<br /></div>
<div dir="ltr">
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-u8DX96X_I3o/Xbzq2vmx5rI/AAAAAAAAVlQ/RxqnS25VBQ8rfacovtleZUdFHXZAioltgCLcBGAsYHQ/s1600/EJ21.1_Mapa_Sismicidad_1.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1002" data-original-width="1600" height="250" src="https://1.bp.blogspot.com/-u8DX96X_I3o/Xbzq2vmx5rI/AAAAAAAAVlQ/RxqnS25VBQ8rfacovtleZUdFHXZAioltgCLcBGAsYHQ/s400/EJ21.1_Mapa_Sismicidad_1.png" width="400" /></a></div>
<br />
<b><br /></b>
<b>2. Perfil topográfico (superior)</b>: Se utilizan los comandos explicados en <a href="https://mapasgeologicos.blogspot.com/2016/05/111-perfil-topografico.html">11.1</a><a href="https://www.blogger.com/null">. Perfil Topográfico</a> para crear el perfil, extraer los datos del dem, calcular cuanto mide (en km), graficar el perfil con la exageración vertical y su orientación.<br />
<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-krnFYxdC-jw/XcGJW4gFf0I/AAAAAAAAVmA/SUQCRcC6NkIO5H7nIVn_PeD0N38LjRE1gCLcBGAsYHQ/s1600/EJ21.1_Perfil_Wadati-Benioff_Longitud-Profundidad_0.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="260" data-original-width="1600" height="61" src="https://1.bp.blogspot.com/-krnFYxdC-jw/XcGJW4gFf0I/AAAAAAAAVmA/SUQCRcC6NkIO5H7nIVn_PeD0N38LjRE1gCLcBGAsYHQ/s400/EJ21.1_Perfil_Wadati-Benioff_Longitud-Profundidad_0.png" width="400" /></a></div>
<br />
<div class="separator" style="clear: both; text-align: center;">
<b><br /></b></div>
<b>3. Perfil en profundidad (inferior)</b>: Se adjunta un perfil en profundidad hasta los 300 km de profundidad. <br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;"><br /></span></span>
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">SET DepthMin=0</span></span><br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">SET DepthMax=300 </span></span><br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;"><br /></span></span>
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psbasemap -R -J -O -K -Bxaf+l"Longitud (\260)" >> %OUT% --MAP_FRAME_AXES=Sn</span></span><br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psbasemap -R -J -O -K -Byaf+l"Profundidad (km)" >> %OUT% --FONT_ANNOT_PRIMARY=8,Helvetica,blue --MAP_FRAME_AXES=wE </span></span><br />
<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-V3U_MedNo98/XcGJW6cvn_I/AAAAAAAAVmE/A9jDURDq0p4CncV55xcELbGXRs3y_i3UQCEwYBhgL/s1600/EJ21.1_Perfil_Wadati-Benioff_Longitud-Profundidad_1.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="835" data-original-width="1600" height="207" src="https://1.bp.blogspot.com/-V3U_MedNo98/XcGJW6cvn_I/AAAAAAAAVmE/A9jDURDq0p4CncV55xcELbGXRs3y_i3UQCEwYBhgL/s400/EJ21.1_Perfil_Wadati-Benioff_Longitud-Profundidad_1.png" width="400" /></a></div>
<br />
<br />
<b>4. Filtrar datos</b>: Con los mismos comandos del punto 2, se filtran los datos de los sismos y de los mecanismos focales según la región geográfica entre los 31 y 33 latitud sur.<br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">SET REGION2=%LongMin%/%LongMax%/-33/-31<br />gmt select -R%REGION2% "query_*.csv" -h1 -i2,1,3 > "temp_sismos"<br />gmt select -R%REGION2% "CMT_*.txt" > "temp_CMT"</span></span><br />
<br />
<b>5. Plotear sismos en profundidad versus longitud</b>: con psxy se grafican los sismos previamente filtrados (temp_sismos). Es
importante indicar las columnas con los datos de longitud y profundidad
(-i0,2). Se utiliza una CPT con valores entre 0 y 300 para pintar los sismos según su profundidad. <br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt makecpt > "temp_seis.cpt" -Crainbow -T0/300 -Z -I </span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt psxy -R -J -O -K "temp_sismos" >> %OUT% -W.1 -C"temp_seis" -Sc0.05c -i0,2,2</span></span><br />
<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-X0SE19NevB8/XcGJXBJe0DI/AAAAAAAAVmI/zesXZ41IiYkh0PCGSVDXmmtNjLCMci7KwCEwYBhgL/s1600/EJ21.1_Perfil_Wadati-Benioff_Longitud-Profundidad_2.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="835" data-original-width="1600" height="207" src="https://1.bp.blogspot.com/-X0SE19NevB8/XcGJXBJe0DI/AAAAAAAAVmI/zesXZ41IiYkh0PCGSVDXmmtNjLCMci7KwCEwYBhgL/s400/EJ21.1_Perfil_Wadati-Benioff_Longitud-Profundidad_2.png" width="400" /></a></div>
</div>
<div dir="ltr">
<br /></div>
<div dir="ltr">
<b>5. Mecanismos focales</b>: <a href="https://docs.generic-mapping-tools.org/latest/supplements/seis/pscoupe.html">pscoupe</a> permite graficar los mecanimos focales previamente filtrados (temp_CMT) en un perfil con los siguientes argumentos:<br />
<ul style="text-align: left;">
<li>-Aa: Perfil (longitud y latitud inicial y final), la inclinación del perfil (90 = vertical), extensión hacia los lados del perfil (10 mil km, para incluir a todos los mecanismos) y profundidad mínima y máxima. </li>
<li>-Sd0.15/0: <i>beachballs</i> de un tamaño de 0,15 (para los eventos de magnitud 5
en pulgadas por default). El segundo valor corresponde al tamaño de la
fuente usado para el nombre del evento. Con 0 no se escribe. </li>
<li>-Q: No crea los archivos de información de los mecanismos focales.</li>
<li>-G: Color para los cuadrantes compresivos.</li>
<li>-M: Tamaño fijo.</li>
</ul>
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt pscoupe -R -J -O -K "temp_CMT" >> %OUT% -Sd0.15/0 -Gred -M -Aa%LongMin%/%LatMin%/%LongMax%/%LatMax%/90/10000/%DepthMin%/%DepthMax%f -Q</span></span><br />
<br />
<br />
<div style="text-align: center;">
<a href="https://1.bp.blogspot.com/-YIGOhgpSfR4/XcGJXrLWrRI/AAAAAAAAVmM/jSapp3HMaNEmGcJwk5EGliUMYGjJPNPMgCEwYBhgL/s1600/EJ21.1_Perfil_Wadati-Benioff_Longitud-Profundidad_3.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="835" data-original-width="1600" height="207" src="https://1.bp.blogspot.com/-YIGOhgpSfR4/XcGJXrLWrRI/AAAAAAAAVmM/jSapp3HMaNEmGcJwk5EGliUMYGjJPNPMgCEwYBhgL/s400/EJ21.1_Perfil_Wadati-Benioff_Longitud-Profundidad_3.png" width="400" /></a></div>
<b><br /></b>
<b>6. Otros Perfiles Longitudinales</b>:<br />
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-kBLKf6k_a84/XcTNXM6Jj9I/AAAAAAAAVnM/Wr2AzAAhcZIwqdBFJWSemk22l5venh0JwCLcBGAsYHQ/s1600/EJ21.1_Perfil_Wadati-Benioff_-28.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="835" data-original-width="1600" height="207" src="https://1.bp.blogspot.com/-kBLKf6k_a84/XcTNXM6Jj9I/AAAAAAAAVnM/Wr2AzAAhcZIwqdBFJWSemk22l5venh0JwCLcBGAsYHQ/s400/EJ21.1_Perfil_Wadati-Benioff_-28.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Perfil a los 28 S entre 74 y 62 L y hasta 700 km de profundidad.</td></tr>
</tbody></table>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-fORtlcd7asY/XcTNXAsgyfI/AAAAAAAAVnQ/eAdzc0-AiGcm_Fp-Sf_xIisEMMAWl8KxQCLcBGAsYHQ/s1600/EJ21.1_Perfil_Wadati-Benioff_-30.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="835" data-original-width="1600" height="207" src="https://1.bp.blogspot.com/-fORtlcd7asY/XcTNXAsgyfI/AAAAAAAAVnQ/eAdzc0-AiGcm_Fp-Sf_xIisEMMAWl8KxQCLcBGAsYHQ/s400/EJ21.1_Perfil_Wadati-Benioff_-30.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Perfil a los 28 S entre 74 y 42 L y hasta 300 km de profundidad.</td></tr>
</tbody></table>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-NT41SmDRZgU/XcTNX68TgjI/AAAAAAAAVnY/09nJ4O6pzGs_wGXcvQXt_4UGUDfmZMSFwCLcBGAsYHQ/s1600/EJ21.1_Perfil_Wadati-Benioff_-34.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="835" data-original-width="1600" height="207" src="https://1.bp.blogspot.com/-NT41SmDRZgU/XcTNX68TgjI/AAAAAAAAVnY/09nJ4O6pzGs_wGXcvQXt_4UGUDfmZMSFwCLcBGAsYHQ/s400/EJ21.1_Perfil_Wadati-Benioff_-34.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Perfil a los 28 S entre 74 y 42 L y hasta 300 km de profundidad.</td></tr>
</tbody></table>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-lfKjA8x8GtE/XcTNXQuv9rI/AAAAAAAAVnU/opSXxwgrwp8pOpcCmRsYS6o-MPvj8r9ggCLcBGAsYHQ/s1600/EJ21.1_Perfil_Wadati-Benioff_-28A.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="835" data-original-width="1600" height="207" src="https://1.bp.blogspot.com/-lfKjA8x8GtE/XcTNXQuv9rI/AAAAAAAAVnU/opSXxwgrwp8pOpcCmRsYS6o-MPvj8r9ggCLcBGAsYHQ/s400/EJ21.1_Perfil_Wadati-Benioff_-28A.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Perfil a los 28 S entre 74 y 42 L y hasta 300 km de profundidad.</td></tr>
</tbody></table>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
</div>
</div>
Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com2Argentina-40.008154981029357 -62.2658601859938-66.714981481029355 -103.5744541859938 -13.301328481029358 -20.957266185993802tag:blogger.com,1999:blog-1509523996050466325.post-57303199733947823112019-10-23T00:00:00.000-03:002020-04-11T17:38:12.398-03:0011.5. Multiples perfiles geofísicos en mapa<div dir="ltr" style="text-align: left;" trbidi="on">
Graficar múltiples perfiles en un mapa.<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<b><a href="https://1.bp.blogspot.com/-3VSlyQOn0YM/XafIMd-lfdI/AAAAAAAAVZo/l28j2EN4m8YJIXeLIlcBxIEnuJEqUEJwACLcBGAsYHQ/s1600/EJ11.5_Perfiles%252BMapa_5000-1-200.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1600" data-original-width="1586" height="640" src="https://1.bp.blogspot.com/-3VSlyQOn0YM/XafIMd-lfdI/AAAAAAAAVZo/l28j2EN4m8YJIXeLIlcBxIEnuJEqUEJwACLcBGAsYHQ/s640/EJ11.5_Perfiles%252BMapa_5000-1-200.png" width="633" /></a></b></div>
<b></b><br />
<a name='more'></a><a href="https://www.blogger.com/null"></a><b>Introducción:</b>En el ejemplo anterior (<a href="https://mapasgeologicos.blogspot.com/2019/10/114-perfil-geofisico-en-mapa.html">11.4</a>) se muestra como graficar un perfil de anomalías sobre un mapa a partir de una linea. En este, se muestra como graficar múltiples perfiles perpendiculares a un perfil guía.<br />
<br />
<a href="https://www.dropbox.com/s/7ti88dc0ogsb2rv/EJ11.5_Perfiles%2BMapa.bat?dl=0">Script</a><br />
<br />
<div dir="ltr">
<b>0. Datos</b>:<br />
<ul style="text-align: left;">
<li>DEM: <a href="https://www.dropbox.com/s/wkch1cvb37aplf7/GMRTv3_5.grd?dl=0">GMRTv3_5.grd</a></li>
<li>Anomalías de Aire Libre: <a href="ftp://topex.ucsd.edu/pub/global_grav_1min/grav.img.23.1">grav.img.27.1 </a>(en formato img). Ver <a href="http://mapasgeologicos.blogspot.com.ar/p/blog-page_70.html">Base de Datos</a>.</li>
<li>Trinchera: Datos de Plates Project (Trench_PlateProject.txt). Ver <a href="http://mapasgeologicos.blogspot.com.ar/p/blog-page_70.html">Base de Datos</a>.</li>
</ul>
</div>
<b>1. Mapa Base</b>: Se crea con los siguientes comandos (ver <a href="http://mapasgeologicos.blogspot.com/2016/02/5-mapa-topografico.html">5.1</a> y <a href="http://mapasgeologicos.blogspot.com/2016/03/52-mapa-topografico-ii.html">5.2</a> para más detalles). <br />
<br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt grdcut "GMRTv3_5.grd" -G%CUT% -R </span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt makecpt -Cibcso -T-8400/0/10 -Z -N > %color%<br />gmt makecpt -Celevation -T0/6050/50 -Z >> %color%</span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt grdgradient %CUT% -A270 -G%SHADOW% -Ne0.5</span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt grdimage -R -J -O -K %CUT% -C%color% -I%SHADOW% >> %OUT%</span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt psscale -R -J -O -K -DJRM+o0.3c/0+w15/0.618c -C%color% -Ba1+l"Elevaciones (km)" -I >> %OUT% -W0.001</span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt psxy -R -J -O -K >> %OUT% "Trench_PlateProject.txt" -W1,red</span></span><br />
<br />
<div style="text-align: center;">
<a href="https://1.bp.blogspot.com/-uutf-qfztFc/XaKW7aqWbcI/AAAAAAAAVTY/RGm1dXnDpEokqZqI1m-2RPk1NECQKmnxgCLcBGAsYHQ/s1600/EJ11.5_Perfiles%252BMapa_0.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1600" data-original-width="1586" height="400" src="https://1.bp.blogspot.com/-uutf-qfztFc/XaKW7aqWbcI/AAAAAAAAVTY/RGm1dXnDpEokqZqI1m-2RPk1NECQKmnxgCLcBGAsYHQ/s400/EJ11.5_Perfiles%252BMapa_0.png" width="396" /></a> </div>
<br />
<b>2. Datos Gravimétricos</b>: Con <a href="http://gmt.soest.hawaii.edu/doc/latest/supplements/img/img2grd.html">img2grd</a> se convierten los datos gravimétricos (anomalías de Aire Libre) para su posterior extracción. Ver el punto 4 de <a href="https://mapasgeologicos.blogspot.com/2016/04/8-mapa-gravimetrico-anomalias-de-aire.html">7.1</a> para más detalles.<br />
<br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt img2grd %GRA% -R -G%TGRA% -T1 -I1 -E -S0.1</span></span><br />
<br />
<b>3. Definir Perfil Guía</b>: Definir una línea N-S según el meridiano 55 ºW entre 10 ºS y 55ºS. Los datos se guardan en el archivo "temp_line" con los siguientes comandos. La unión de los puntos se realiza siguiendo un <a href="https://es.wikipedia.org/wiki/Gran_c%C3%ADrculo">círculo máximo</a>.<br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">REM Perfil: Crear archivo para dibujar perfil (Long Lat)</span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">echo -55 -10 > "temp_line"</span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">echo -55 -55 >> "temp_line" </span></span><br />
<br />
<b>4. Extraer datos de los perfiles</b>: con <a href="https://docs.generic-mapping-tools.org/latest/grdtrack.html?highlight=grdtrack">grdtrack</a> se crean los perfiles perpendiculares al perfil guía ("temp_line") con las siguientes opciones. Los datos de los perfiles se guardan en el archivo "temp_data"<br />
<ul style="text-align: left;">
<li>-G%TGRA%: Grilla con datos de anomalías de Aire Libre.</li>
<li>-C5000k/50k/200k: Perfiles perpendiculares equiespaciados cada 200 km según el perfil guía, cada uno de 5000 km de longitud y con datos cada 50 km. </li>
</ul>
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt grdtrack "temp_line" -G%TGRA% -C5000k/50k/200k > "temp_data"</span></span><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div style="text-align: left;">
<b>5. Graficar perfiles en el mapa</b>: <a href="https://docs.generic-mapping-tools.org/latest/pswiggle.html">pswiggle</a> permite graficar los perfiles con las siguientes opciones:<br />
<ul style="text-align: left;">
<li>-i0,1,4: Selecciona las columnas con los datos del archivo "temp_data"</li>
<li>-F: Agrega recuadro blanco debajo de la escala (+gwhite), con un borde (+p) y efecto sombreado (+s).</li>
<li>-Gred@50+p: Pinta de rojo las anomalías <b>positivas</b> con 50% de transparencia (@).</li>
<li>-Gblue@50+n: Pinta de azul las anomalías <b>negativas</b> con 50% de transparencia (@).</li>
<li>-Z500: Escala de las anomalías.</li>
<li>-D: Dibuja una escala de 100 mGal abajo a la derecha (en este ejemplo).</li>
<li>-T: Dibuja el perfil.</li>
<li>-W: Dibuja el borde de las anomalías.</li>
</ul>
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt pswiggle -R -J -O -K "temp_data" -i0,1,4 >> %OUT% -Gred@50+p -Gblue@50+n -Z500 -T -W -DjRB+o0.5/0.5+w100+lmGal -F+gwhite+p+s</span></span><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
</div>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-Z4Fgx7z9_g0/XafT4v2htwI/AAAAAAAAVaA/dH4-S2CQzHsgdQYV8Z7RtPVPAjAIhSqLgCLcBGAsYHQ/s1600/EJ11.5_Perfiles%252BMapa_2000-100-1000.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1586" height="400" src="https://1.bp.blogspot.com/-Z4Fgx7z9_g0/XafT4v2htwI/AAAAAAAAVaA/dH4-S2CQzHsgdQYV8Z7RtPVPAjAIhSqLgCLcBGAsYHQ/s400/EJ11.5_Perfiles%252BMapa_2000-100-1000.png" width="395" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">-C2000/100/100</td></tr>
</tbody></table>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-V19TDAHWkIk/XafT2wys0pI/AAAAAAAAVZ8/RLrFXOIx_FwrPU4HYmMP8hYwW3KtKAGgwCLcBGAsYHQ/s1600/EJ11.5_Perfiles%252BMapa_2000-10-500.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1586" height="400" src="https://1.bp.blogspot.com/-V19TDAHWkIk/XafT2wys0pI/AAAAAAAAVZ8/RLrFXOIx_FwrPU4HYmMP8hYwW3KtKAGgwCLcBGAsYHQ/s400/EJ11.5_Perfiles%252BMapa_2000-10-500.png" width="396" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">-C2000k/10k/500k</td></tr>
</tbody></table>
<br />
<div class="" style="clear: both; text-align: center;">
</div>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-vWTWVq1NJSc/XafT2UuTGtI/AAAAAAAAVZ4/MPmqOSTkgwk1OuTyJg1s-LVRyDIrz1CMwCLcBGAsYHQ/s1600/EJ11.5_Perfiles%252BMapa_2000-10-200.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1586" height="400" src="https://1.bp.blogspot.com/-vWTWVq1NJSc/XafT2UuTGtI/AAAAAAAAVZ4/MPmqOSTkgwk1OuTyJg1s-LVRyDIrz1CMwCLcBGAsYHQ/s400/EJ11.5_Perfiles%252BMapa_2000-10-200.png" width="396" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">-C2000/10/200</td></tr>
</tbody></table>
<div class="separator" style="clear: both; text-align: center;">
</div>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-1KpGwDi6sh4/XafT8hRjlzI/AAAAAAAAVaE/453ul6iNSe8y09uU6sl_0utiCgVqTu0IACLcBGAsYHQ/s1600/EJ11.5_Perfiles%252BMapa_5000-100-1000.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1586" height="400" src="https://1.bp.blogspot.com/-1KpGwDi6sh4/XafT8hRjlzI/AAAAAAAAVaE/453ul6iNSe8y09uU6sl_0utiCgVqTu0IACLcBGAsYHQ/s400/EJ11.5_Perfiles%252BMapa_5000-100-1000.png" width="396" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">-C5000k/100k/1000k</td></tr>
</tbody></table>
<div class="separator" style="clear: both; text-align: center;">
</div>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-SqCwLMZL6oU/XafT8_JLQxI/AAAAAAAAVaI/kzLSgHJ440oQtC5-D3FK75bGxt8IVg6qQCLcBGAsYHQ/s1600/EJ11.5_Perfiles%252BMapa_5000-10-200.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1586" height="400" src="https://1.bp.blogspot.com/-SqCwLMZL6oU/XafT8_JLQxI/AAAAAAAAVaI/kzLSgHJ440oQtC5-D3FK75bGxt8IVg6qQCLcBGAsYHQ/s400/EJ11.5_Perfiles%252BMapa_5000-10-200.png" width="396" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">-C5000k/50k/200k</td></tr>
</tbody></table>
<br />
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<b>6. Argentina-Chile</b>: Ejemplos acotando la región a Argentina-Chile.<br />
<br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">SET REGION=-88/-42/-56/-21</span></span><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-ut-sxOsWgp8/Xa8p-xALz_I/AAAAAAAAVeM/Tv3KW8nVEe4KSMwRu6he-JVRQrHrWR2lgCLcBGAsYHQ/s1600/EJ11.5_Perfiles%252BMapa_Argentina_1.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1358" data-original-width="1600" height="338" src="https://1.bp.blogspot.com/-ut-sxOsWgp8/Xa8p-xALz_I/AAAAAAAAVeM/Tv3KW8nVEe4KSMwRu6he-JVRQrHrWR2lgCLcBGAsYHQ/s400/EJ11.5_Perfiles%252BMapa_Argentina_1.png" width="400" /></a></div>
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/--7-fa5MFJnw/Xa8p_IEq6bI/AAAAAAAAVeU/gikbHbtyXUoDIqWCDtxyBX3n7DFwsMccgCLcBGAsYHQ/s1600/EJ11.5_Perfiles%252BMapa_Argentina_2.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1358" data-original-width="1600" height="338" src="https://1.bp.blogspot.com/--7-fa5MFJnw/Xa8p_IEq6bI/AAAAAAAAVeU/gikbHbtyXUoDIqWCDtxyBX3n7DFwsMccgCLcBGAsYHQ/s400/EJ11.5_Perfiles%252BMapa_Argentina_2.png" width="400" /></a></div>
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-uYDp-eMLu88/Xa8p-3dcZAI/AAAAAAAAVeQ/2beXT40Knd8og49abD6Ng52Q8jEbQ9JSQCLcBGAsYHQ/s1600/EJ11.5_Perfiles%252BMapa_Argentina_3.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1358" data-original-width="1600" height="338" src="https://1.bp.blogspot.com/-uYDp-eMLu88/Xa8p-3dcZAI/AAAAAAAAVeQ/2beXT40Knd8og49abD6Ng52Q8jEbQ9JSQCLcBGAsYHQ/s400/EJ11.5_Perfiles%252BMapa_Argentina_3.png" width="400" /></a></div>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-4HFSYeYDvNk/Xa8p_7dkjDI/AAAAAAAAVeY/nVzIaNyhrnUMYz2HaQZiMuZoUHiacpTPwCLcBGAsYHQ/s1600/EJ11.5_Perfiles%252BMapa_Argentina_4.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1358" data-original-width="1600" height="338" src="https://1.bp.blogspot.com/-4HFSYeYDvNk/Xa8p_7dkjDI/AAAAAAAAVeY/nVzIaNyhrnUMYz2HaQZiMuZoUHiacpTPwCLcBGAsYHQ/s400/EJ11.5_Perfiles%252BMapa_Argentina_4.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Sólo se pinta la anomalía (sin -T ni -W).</td></tr>
</tbody></table>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
</div>
</div>
Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0Buenos Aires, CABA, Argentina-34.6036844 -58.381559100000004-34.8127694 -58.704282600000006 -34.3945994 -58.0588356tag:blogger.com,1999:blog-1509523996050466325.post-16755789476190777172019-10-01T00:25:00.002-03:002020-04-11T17:39:44.639-03:0011.4. Perfil Geofísico en mapa<div dir="ltr" style="text-align: left;" trbidi="on">
<div dir="ltr">
Dibujar perfil gravimétrico de anomalías de aire libre en un mapa.<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-NAbQhaZGs1I/XZF5b1dZMEI/AAAAAAAAVK8/RjUZzzXEMiw5vapIE46awUcjUtzxigokgCEwYBhgL/s1600/EJ11.4_PerfilGeofisico%252BMapa_3_Z300_%252BT%252BW.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="588" data-original-width="1600" height="234" src="https://1.bp.blogspot.com/-NAbQhaZGs1I/XZF5b1dZMEI/AAAAAAAAVK8/RjUZzzXEMiw5vapIE46awUcjUtzxigokgCEwYBhgL/s640/EJ11.4_PerfilGeofisico%252BMapa_3_Z300_%252BT%252BW.png" width="640" /></a></div>
<a name='more'></a><b>Introducción</b>: En ocasiones puede ser útil graficar datos de un perfil en un mapa. Se explica como extraer los datos de anomalías de Aire Libre a lo largo de un perfil y graficarlos sobre un mapa batimétrico-topográfico.<br />
<br />
<u>Nota</u>: Es el primer ejemplo con la versión GMT6. Puede haber ligeros cambios en la sintaxis en algunas herramientas. <br />
<br />
<u>Versión GMT</u>: 6.0.0rc4<br />
<br />
<a href="https://www.dropbox.com/s/390tb01gt1poz7o/EJ11.4_PerfilGeofisico%2BMapa.bat?dl=0">Script</a><br />
<br /></div>
<div dir="ltr">
<b>0. Datos</b>:<br />
<ul style="text-align: left;">
<li>DEM: GEBCO_2014</li>
<li>Anomalías de Aire Libre: <a href="ftp://topex.ucsd.edu/pub/global_grav_1min/grav.img.23.1">grav.img.23.1 </a>(en formato img). Ver <a href="http://mapasgeologicos.blogspot.com.ar/p/blog-page_70.html">Base de Datos</a>.</li>
</ul>
</div>
<div dir="ltr">
<b>1. Mapa Base</b>: Se crea con los siguientes comandos (ver <a href="http://mapasgeologicos.blogspot.com/2016/02/5-mapa-topografico.html">5.1</a> y <a href="http://mapasgeologicos.blogspot.com/2016/03/52-mapa-topografico-ii.html">5.2</a> para más detalles).<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: left;">
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">grdcut %DEM% -G%CUT% -R%REGION%</span><span style="font-family: "courier new" , "courier" , monospace;"> </span></span></span></div>
<div class="separator" style="clear: both; text-align: left;">
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;"><span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;"></span></span></span></span></span></div>
<div class="separator" style="clear: both; text-align: left;">
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">makecpt -Crelief -G-5850/2500 -T-5850/2500/50 -Z -D >%color%<br />gmt grdgradient %CUT% -A90 -G%SHADOW% -Ne0.5h<br />gmt grdimage -R -J -O -K %CUT% -C%color% -I%SHADOW% >> %OUT%</span></span></div>
<div style="text-align: center;">
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;"></span></span><br />
<a href="https://1.bp.blogspot.com/-KWlPDVIf1tI/XZLBdQjNxoI/AAAAAAAAVLE/zYLMYrJGgLgt3yC20o-r0gALr61IkOx2ACLcBGAsYHQ/s1600/EJ11.4_PerfilGeofisico%252BMapa.png"></a><a href="https://1.bp.blogspot.com/-KWlPDVIf1tI/XZLBdQjNxoI/AAAAAAAAVLE/zYLMYrJGgLgt3yC20o-r0gALr61IkOx2ACLcBGAsYHQ/s1600/EJ11.4_PerfilGeofisico%252BMapa.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="588" data-original-width="1600" height="233" src="https://1.bp.blogspot.com/-KWlPDVIf1tI/XZLBdQjNxoI/AAAAAAAAVLE/zYLMYrJGgLgt3yC20o-r0gALr61IkOx2ACLcBGAsYHQ/s640/EJ11.4_PerfilGeofisico%252BMapa.png" width="640" /> </a></div>
<br />
<b>2. Definir Perfil</b>: Crear un archivo (temp_line) con las coordenadas del perfil. La unión de los puntos se realiza siguiendo un <a href="https://es.wikipedia.org/wiki/Gran_c%C3%ADrculo">círculo máximo</a>.<br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">REM Perfil: Crear archivo para dibujar perfil (Long Lat)</span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">echo -76 -32 > "temp_line"</span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">echo -46 -32 >> "temp_line" </span></span><br />
<br />
<b>3. Datos Gravimétricos</b>: Se crea el archivo "temp_data" con los datos de las anomalias de aire libre. Ver <a href="https://mapasgeologicos.blogspot.com/2016/05/111-perfil-topografico.html">11.1</a> para la explicación de los comandos y <a href="https://mapasgeologicos.blogspot.com/2016/04/8-mapa-gravimetrico-anomalias-de-aire.html">7.1</a> para el comando img2grd.<br />
<br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt img2grd %GRA% -R -G%TGRA% -T1 -I1 -E -S0.1<br />gmt sample1d "temp_line" -I0.2k > "temp_sample1d"<br />gmt mapproject "temp_sample1d" -G+k > "temp_track" <br />gmt grdtrack "temp_track" -G%TGRA% > "temp_data"<br />gmtinfo "temp_data"</span></span><u></u><br />
<div style="text-align: left;">
<u><br /></u>
<b>3. Graficar perfil en el mapa</b>: El comando <a href="https://docs.generic-mapping-tools.org/latest/pswiggle.html">pswiggle</a> permite graficar el perfil con las siguientes opciones:</div>
<ul style="text-align: left;">
<li>-Gred+p: Pinta de rojo las anomalías <b>positivas</b>.</li>
<li>-Gblue+n: Pinta de azul las anomalías <b>negativas</b>.</li>
<li>-Z300: Escala de las anomalías.</li>
<li>-D: Dibuja una escala de 100 mGal abajo a la derecha (en este ejemplo).</li>
<li>-T: Dibuja el perfil.</li>
<li>-W: Dibuja el borde de las anomalías</li>
<li>-t50: transparencia (50%).</li>
<li>-i0,1,3: Selecciona las columnas donde están los datos en el archivo "temp_data" de longitud, latitud y anomalía de Aire Libre. </li>
</ul>
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt pswiggle -R -J -O -K "temp_data" -i0,1,3 >> %OUT% -Gred+p -Gblue+n -Z300 -DjRB+o0.1/0.1+w100+lmGal -t50 -T -W</span></span><u><br /></u><br />
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-qqoJpN3FLQo/XZF0ONuYS2I/AAAAAAAAVKg/UAas3Pjqllkqss7zqlRxIF3w-R0M8gWjACLcBGAsYHQ/s1600/EJ11.4_PerfilGeofisico%252BMapa_Z300.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="588" data-original-width="1600" height="234" src="https://1.bp.blogspot.com/-qqoJpN3FLQo/XZF0ONuYS2I/AAAAAAAAVKg/UAas3Pjqllkqss7zqlRxIF3w-R0M8gWjACLcBGAsYHQ/s640/EJ11.4_PerfilGeofisico%252BMapa_Z300.png" width="640" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">pswiggle sin -T ni -W.</td></tr>
</tbody></table>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-h9fzi5nx3HM/XZF0NndE8uI/AAAAAAAAVKs/pRifkKJ3sYMu4I4eBFGSk1IZwXqfKnjRQCEwYBhgL/s1600/EJ11.4_PerfilGeofisico%252BMapa_Z300%252BT.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="588" data-original-width="1600" height="234" src="https://1.bp.blogspot.com/-h9fzi5nx3HM/XZF0NndE8uI/AAAAAAAAVKs/pRifkKJ3sYMu4I4eBFGSk1IZwXqfKnjRQCEwYBhgL/s640/EJ11.4_PerfilGeofisico%252BMapa_Z300%252BT.png" width="640" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">pswiggle con -T y sin -W.</td></tr>
</tbody></table>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-NAbQhaZGs1I/XZF5b1dZMEI/AAAAAAAAVK4/Che0dmbG7uolenRsaMwtUVzjXsMvXtFzwCLcBGAsYHQ/s1600/EJ11.4_PerfilGeofisico%252BMapa_3_Z300_%252BT%252BW.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="588" data-original-width="1600" height="234" src="https://1.bp.blogspot.com/-NAbQhaZGs1I/XZF5b1dZMEI/AAAAAAAAVK4/Che0dmbG7uolenRsaMwtUVzjXsMvXtFzwCLcBGAsYHQ/s640/EJ11.4_PerfilGeofisico%252BMapa_3_Z300_%252BT%252BW.png" width="640" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">pswiggle con -T y -W.</td></tr>
</tbody></table>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-z4oHK2vJrE0/XZF0M-5DUyI/AAAAAAAAVKk/gzze5QLll_ofQdl3Fdd-2S6szPoMB3XOACEwYBhgL/s1600/EJ11.4_PerfilGeofisico%252BMapa.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="633" data-original-width="1600" height="252" src="https://1.bp.blogspot.com/-z4oHK2vJrE0/XZF0M-5DUyI/AAAAAAAAVKk/gzze5QLll_ofQdl3Fdd-2S6szPoMB3XOACEwYBhgL/s640/EJ11.4_PerfilGeofisico%252BMapa.png" width="640" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Con proyección UTM 20 sur (U-20/15c).</td></tr>
</tbody></table>
</div>
</div>
Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0tag:blogger.com,1999:blog-1509523996050466325.post-27647068324970825082019-09-10T00:00:00.000-03:002020-04-07T12:20:18.288-03:005.6. Mapa batimétrico-topográfico. Combinar CPT.<div dir="ltr" style="text-align: left;" trbidi="on">
Como combinar paletas de colores (CPT) para utilizarlos en mapas topobatimetricos.<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-Wr2PQTD6FTQ/XWSXgsfjeEI/AAAAAAAAU30/QAUkQZywA7kQyzd_9CHdFigJJXkt4i2xQCLcBGAs/s1600/EJ5.6_CPT_Combinado_1.png"><img border="0" data-original-height="1600" data-original-width="1411" height="640" src="https://1.bp.blogspot.com/-Wr2PQTD6FTQ/XWSXgsfjeEI/AAAAAAAAU30/QAUkQZywA7kQyzd_9CHdFigJJXkt4i2xQCLcBGAs/s640/EJ5.6_CPT_Combinado_1.png" width="563" /></a></div>
<br />
<div style="text-align: center;">
</div>
<br />
<a name='more'></a>Se explica como crear una nueva Paleta de Color (CPT) a partir de la combinación de otras dos, por ejemplo para hacer un mapas batimetrico-topográfico.<br />
<div dir="ltr" style="text-align: left;" trbidi="on">
<br />
<br />
Este ejemplo es una adaptación del explicado en <a href="http://mapasgeologicos.blogspot.com/2016/02/5-mapa-topografico.html">5.1. Mapa Topográfico</a>. En <a href="https://mapasgeologicos.blogspot.com/2016/12/54-mapas-topografico-paleta-de-color-cpt.html">5.4</a> se explican el comando makecpt. En este ejemplo sólo se explican como combinar los CPT.</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<br /></div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<a href="https://www.dropbox.com/s/a5o1u479u2b5dhl/EJ5.6_CPT_Combinado.bat?dl=0">Script</a></div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<br />
<b>Datos</b>: <a href="https://www.dropbox.com/s/wkch1cvb37aplf7/GMRTv3_5.grd?dl=0">GMRT3_5.grd</a> <br />
<br />
1. Para combinar los CPT utilizamos <a href="http://gmt.soest.hawaii.edu/doc/5.4.5/makecpt.html">makecpt </a>dos veces. Para su correcta combinación es importante cumplir los siguientes puntos:<br />
<ul style="text-align: left;">
<li>Primer CPT: </li>
<ul>
<li><b>Valores menores</b>: En este caso de -8400 a 0 dado que corresponde a la batimetría.</li>
<li>Agregar <b>-N</b>: Omite que se escriban los para colores valores para <i>background</i>, <i>foreground</i> y NaN.</li>
<li>Redireccionar con "<b>></b>". Crea el archivo de salida (%color%). </li>
</ul>
</ul>
<ul style="text-align: left;">
<li>Segundo CPT: </li>
<ul>
<li><b>Valores mayores</b>: En este caso de 0 (donde terminó el primer CPT) hasta 6500. </li>
<li>Redireccionar con "<b>>></b>": Agrega los nuevos valores al mismo archivo (%color%). </li>
</ul>
</ul>
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<br />
1. Combinación 1:<br />
<span style="background-color: #cccccc;"><br /></span>
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt makecpt -Cabyss -T-8400/0/10 -Z -N > %color%</span></span><br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt makecpt -Cdem2 -T0/6050/50 -Z >> %color%</span></span><br />
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-Wr2PQTD6FTQ/XWSXgsfjeEI/AAAAAAAAU30/QAUkQZywA7kQyzd_9CHdFigJJXkt4i2xQCLcBGAs/s1600/EJ5.6_CPT_Combinado_1.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1411" height="400" src="https://1.bp.blogspot.com/-Wr2PQTD6FTQ/XWSXgsfjeEI/AAAAAAAAU30/QAUkQZywA7kQyzd_9CHdFigJJXkt4i2xQCLcBGAs/s400/EJ5.6_CPT_Combinado_1.png" width="352" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Combinación abyss y dem2.</td></tr>
</tbody></table>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-deN9mG5_Ub8/XWSXiCH7qpI/AAAAAAAAU38/p8Yv5oRL3AAW1Q6x64LqxHN3KlrFHEyrQCLcBGAs/s1600/EJ5.6_CPT_Combinado.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1411" height="400" src="https://1.bp.blogspot.com/-deN9mG5_Ub8/XWSXiCH7qpI/AAAAAAAAU38/p8Yv5oRL3AAW1Q6x64LqxHN3KlrFHEyrQCLcBGAs/s400/EJ5.6_CPT_Combinado.png" width="350" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Combinación arendal y wiki.</td></tr>
</tbody></table>
<br />
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-JOmYeuyeKCE/XWSXhXHPzwI/AAAAAAAAU34/9VY37VeuI4c2MJMB1RXUu58l7INMT6qXwCLcBGAs/s1600/EJ5.6_CPT_Combinado_2.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1411" height="400" src="https://1.bp.blogspot.com/-JOmYeuyeKCE/XWSXhXHPzwI/AAAAAAAAU34/9VY37VeuI4c2MJMB1RXUu58l7INMT6qXwCLcBGAs/s400/EJ5.6_CPT_Combinado_2.png" width="352" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Combinación ibcso y elevation.</td></tr>
</tbody></table>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-y4iBvNL7rvU/XWSXkasBOLI/AAAAAAAAU4A/LboBYgSkbPoAiroFJSQOlRIqDYqHKV3MACLcBGAs/s1600/EJ5.6_CPT_Combinado_3.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1411" height="400" src="https://1.bp.blogspot.com/-y4iBvNL7rvU/XWSXkasBOLI/AAAAAAAAU4A/LboBYgSkbPoAiroFJSQOlRIqDYqHKV3MACLcBGAs/s400/EJ5.6_CPT_Combinado_3.png" width="352" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Combinación bathy y gray.</td></tr>
</tbody></table>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-OS_leKPj38M/XWSXo4innDI/AAAAAAAAU4E/PsyIwOdPipE8GH7LD7U0BZKJZDviU4EpACLcBGAs/s1600/EJ5.6_CPT_Combinado_3b.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1411" height="400" src="https://1.bp.blogspot.com/-OS_leKPj38M/XWSXo4innDI/AAAAAAAAU4E/PsyIwOdPipE8GH7LD7U0BZKJZDviU4EpACLcBGAs/s400/EJ5.6_CPT_Combinado_3b.png" width="352" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Combinación bathy y gray (invertida).</td></tr>
</tbody></table>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-Ec-dVhz7w1k/XWSXq0WrYdI/AAAAAAAAU4I/F1XE59SMruYGfDuYAmvlmvaaGohtaN03QCLcBGAs/s1600/EJ5.6_CPT_Combinado_4.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1411" height="400" src="https://1.bp.blogspot.com/-Ec-dVhz7w1k/XWSXq0WrYdI/AAAAAAAAU4I/F1XE59SMruYGfDuYAmvlmvaaGohtaN03QCLcBGAs/s400/EJ5.6_CPT_Combinado_4.png" width="352" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Combinación ocean y copper.</td></tr>
</tbody></table>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-sVdW2e6eKos/XWSXr1pDk-I/AAAAAAAAU4M/EfSHk4kocyUzVzQl2ry58kMFH4tCSit9gCLcBGAs/s1600/EJ5.6_CPT_Combinado_5.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1411" height="400" src="https://1.bp.blogspot.com/-sVdW2e6eKos/XWSXr1pDk-I/AAAAAAAAU4M/EfSHk4kocyUzVzQl2ry58kMFH4tCSit9gCLcBGAs/s400/EJ5.6_CPT_Combinado_5.png" width="352" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Combinación oslo y bilbao.</td></tr>
</tbody></table>
<a href="https://www.blogger.com/null" name="more"></a></div>
</div>
Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com2tag:blogger.com,1999:blog-1509523996050466325.post-24774443010989089082019-05-28T00:00:00.000-03:002019-05-28T00:00:07.214-03:0013.6. Bloque 3D. Animación.<div dir="ltr" style="text-align: left;" trbidi="on">
<div dir="ltr" style="text-align: left;" trbidi="on">
<div dir="ltr">
Animación mostrando un bloque 3D de la profundidad del moho en Sudamérica.<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://3.bp.blogspot.com/-OYt7qHVkUd4/XNyjbtDGqEI/AAAAAAAAThU/odaZjxpO9l4kixP1bY2K3Li_FIGAuWMrwCLcBGAs/s1600/3d_Moho_RedBlue_10.gif" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="993" data-original-width="1600" height="247" src="https://3.bp.blogspot.com/-OYt7qHVkUd4/XNyjbtDGqEI/AAAAAAAAThU/odaZjxpO9l4kixP1bY2K3Li_FIGAuWMrwCLcBGAs/s400/3d_Moho_RedBlue_10.gif" width="400" /> </a></div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<a href="https://www.blogger.com/null" name="more"></a><br />
<a href="https://www.blogger.com/null" name="more"></a><br />
<a name='more'></a></div>
<div dir="ltr">
<span style="font-family: inherit;">Esta animación permite mostrar las zonas profundas del moho ocultas por la superficie. Se arma combinando distintos bloques 3D en un gif.</span><br />
<br />
<span style="font-family: inherit;"><a href="https://www.dropbox.com/s/5osbdzbpiejnneb/3d_Moho_RB.bat?dl=0">Script</a> </span></div>
<div dir="ltr">
<br /></div>
<span style="font-family: inherit;"><b>0. Datos</b>: Se utilizó la grilla de profundidad del moho (Más información en <a href="http://mapasgeologicos.blogspot.com.ar/p/blog-page_70.html">Base de Datos</a>):</span><br />
<ul style="text-align: left;">
<li><span style="font-family: inherit;"><b>Profundidad Moho</b>: Tabla de datos depthtomoho.xyz (<a href="https://igppweb.ucsd.edu/~gabi/crust1/depthtomoho.xyz.zip">aquí comprimido en archivo zip</a>) del Modelo <a href="https://igppweb.ucsd.edu/~gabi/crust1.html">CRUST 1.0</a>. Profundidad de km (positivo hacia arriba). </span></li>
</ul>
<div dir="ltr">
<span style="font-family: inherit;"><b>1. Crear grilla: </b>Se utiliza <a href="http://gmt.soest.hawaii.edu/doc/5.4.5/surface.html">surface</a> para crear la grilla Moho.nc (ver. <a href="https://mapasgeologicos.blogspot.com/2018/07/85-mapas-crust-10.html">8.5</a> para más detalles).</span><br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt surface "depthtomoho.xyz" -I6m -R%REGION% -G"Moho.nc" -fg -r </span></span><br />
<br />
<span style="font-family: inherit;"><b>2. Crear CPT personalizada:</b> con makecpt se crea la paleta de color (CPT) que varía gradualmente entre <i>lightred</i> y <i>lightblue</i> entre los valores de 75 y 0 km de profundidad, respectivamente (ver <a href="https://mapasgeologicos.blogspot.com/2017/05/55-mapa-topografico-paleta-de-color.html">5.5</a> para más detalles).</span></div>
<div dir="ltr">
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt makecpt -Clightred,lightblue -T-75/0 -Z > %color%</span></span></div>
<div dir="ltr">
<br />
<span style="font-family: inherit;"><b>3. Efecto Sombreado:</b> Con <a href="http://gmt.soest.hawaii.edu/doc/5.4.5/grdgradient.html">grdgradient</a> se crea una grilla para hacer el efecto de sombreado (ver <a href="https://mapasgeologicos.blogspot.com/2016/03/52-mapa-topografico-ii.html">5.2</a> para más detalles). La </span><br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt grdgradient %CUT% -A%A% -G%SHADOW% -Ne0.5</span></span><br />
<br />
<span style="font-family: inherit;"><b>4. Bloques 3D</b>: <span style="font-family: sans-serif;">Para este ejemplo se crearon 35 figuras. La primera corrresponde a un bloque 3D entero, mientras que para el resto fue nivelado a distintas profundidades entre 10 y 74 km (cada 2 km).</span></span><br />
<span style="font-family: inherit;"><br />
</span></div>
<div dir="ltr">
<span style="font-family: inherit;"><b>A. Bloque 3D entero</b>: Con grdview se crea el bloque 3D (mas detalles en <a href="https://mapasgeologicos.blogspot.com/2016/06/12-bloque-3d.html">13.1</a>).</span><br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt grdview %CUT% -R -J -JZ -O -K -p -Qi300 -C%color% -I%SHADOW% >> %OUT% -N%BASE%+glightgray -Wf0.5 -BnSwEZ+b -Baf -Bzaf10+l"Altura (km)" </span></span><br />
<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-rNy__7gvGso/XOU9mIIQGJI/AAAAAAAATl8/dJfNkXOE4IcW8-lL9ZzKoGCEkSURJJnDwCLcBGAs/s1600/3d_Moho_RB_08.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="993" data-original-width="1600" height="198" src="https://1.bp.blogspot.com/-rNy__7gvGso/XOU9mIIQGJI/AAAAAAAATl8/dJfNkXOE4IcW8-lL9ZzKoGCEkSURJJnDwCLcBGAs/s320/3d_Moho_RB_08.png" width="320" /></a></div>
</div>
<div dir="ltr">
<b><br />
</b> <span style="font-family: inherit;"><b>B. Bloque 3D Nivelado</b>: </span><br />
<span style="font-family: inherit;">Se define la variable Above que es el valor de profundidad (en positivo) a la cual se quiere nivelar el bloque 3D (también se utiliza para definir el nombre de la imagen de salida)<span style="font-family: sans-serif;">. </span>Luego, con <a href="http://gmt.soest.hawaii.edu/doc/5.4.5/grdclip.html">grdclip </a>se nivela la grilla con el siguiente argumento:</span><br />
<ul style="text-align: left;">
<li><span style="font-family: inherit;">-Sa-%above%</span><span style="font-family: inherit;">/-%Above%: Setea todos los datos por encima del valor de la variable %Above% a ese mismo valor. (Se agrega el símbolo "-" porque los datos están en negativo).</span></li>
</ul>
</div>
<div dir="ltr">
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">SET Above=66 </span></span></div>
<div dir="ltr">
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt grdclip "Moho.nc" -G%CUT% -Sa-%Above%/-%Above%</span></span></div>
<div dir="ltr">
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: right;"><tbody>
<tr><td style="text-align: center;"><a href="https://3.bp.blogspot.com/-XVw4Rc0yDP4/XOU-XG9VMjI/AAAAAAAATmI/6hQTPe42fD0KUylKg_xOhjLXpQvm1vQqACLcBGAs/s1600/3d_Moho_RB_32.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="993" data-original-width="1600" height="195" src="https://3.bp.blogspot.com/-XVw4Rc0yDP4/XOU-XG9VMjI/AAAAAAAATmI/6hQTPe42fD0KUylKg_xOhjLXpQvm1vQqACLcBGAs/s320/3d_Moho_RB_32.png" width="320" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Nivelado a 32 km de profundidad.</td></tr>
</tbody></table>
<br />
<table cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://2.bp.blogspot.com/-dFiZl-zjBEI/XOU-XKFI8dI/AAAAAAAATmM/bht5tyG0jCYvFbocdUBxbeiiehhJWqQEQCLcBGAs/s1600/3d_Moho_RB_20.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="993" data-original-width="1600" height="194" src="https://2.bp.blogspot.com/-dFiZl-zjBEI/XOU-XKFI8dI/AAAAAAAATmM/bht5tyG0jCYvFbocdUBxbeiiehhJWqQEQCLcBGAs/s320/3d_Moho_RB_20.png" width="320" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Nivelado a 20 km de profundidad.</td></tr>
</tbody></table>
</div>
<div dir="ltr">
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: left;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-LznuhUUMpnQ/XOU-XvhEtxI/AAAAAAAATmQ/GXGN2i_h9NsUkhQtuv50d85XlcXpOI8sQCLcBGAs/s1600/3d_Moho_RB_56.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="993" data-original-width="1600" height="195" src="https://1.bp.blogspot.com/-LznuhUUMpnQ/XOU-XvhEtxI/AAAAAAAATmQ/GXGN2i_h9NsUkhQtuv50d85XlcXpOI8sQCLcBGAs/s320/3d_Moho_RB_56.png" width="320" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Nivelado a 44 km de profundidad.</td></tr>
</tbody></table>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: right;"><tbody>
<tr><td style="text-align: center;"><a href="https://2.bp.blogspot.com/-sEY0V9QfX8Y/XOU-XJucq_I/AAAAAAAATmE/nScrn4p0dUQ4EHmSVQ1zyYsq0MWXL3s9QCLcBGAs/s1600/3d_Moho_RB_44.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="993" data-original-width="1600" height="196" src="https://2.bp.blogspot.com/-sEY0V9QfX8Y/XOU-XJucq_I/AAAAAAAATmE/nScrn4p0dUQ4EHmSVQ1zyYsq0MWXL3s9QCLcBGAs/s320/3d_Moho_RB_44.png" width="320" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Nivelado a 56 km de profundidad.</td></tr>
</tbody></table>
<br />
<span style="font-family: inherit;"><b>5. Animación</b>: La herramienta <a href="http://www.graphicsmagick.org/convert.html">convert</a> del programa GraphicsMagik permite combinar los 35 bloques 3D en un gif con las siguientes opciones:</span><br />
<br />
<span style="font-family: inherit;"> -delay: intervalo entre imágenes en centésimas de segundos.</span><br />
<span style="font-family: inherit;"> -loop: cuantas veces se repite (0 = infinitamente).</span><br /><br /> </div>
<div dir="ltr">
<div class="separator" style="clear: both; text-align: left;">
<span style="font-family: "courier new" , "courier" , monospace;"><span style="font-family: inherit;"><span style="background-color: #cccccc;">gm convert -delay 10 -loop 0 3d_*.png 3d_Moho_RedBlue_10.gif</span></span></span></div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
Alternativamente se puede utilizar el sitio <a href="https://gifmaker.me/">gifmaker</a>. <span style="font-family: "courier new" , "courier" , monospace;"><span style="font-family: inherit;"><span style="background-color: #cccccc;"><br /></span></span></span></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<br />
<div style="text-align: center;">
<a href="https://1.bp.blogspot.com/-OYt7qHVkUd4/XNyjbtDGqEI/AAAAAAAAThY/np0EHd-iHbUjO4v1YZVgSjjrhtDkj9xdwCPcBGAYYCw/s1600/3d_Moho_RedBlue_10.gif" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="993" data-original-width="1600" height="247" src="https://1.bp.blogspot.com/-OYt7qHVkUd4/XNyjbtDGqEI/AAAAAAAAThY/np0EHd-iHbUjO4v1YZVgSjjrhtDkj9xdwCPcBGAYYCw/s400/3d_Moho_RedBlue_10.gif" width="400" /></a></div>
<span style="font-family: inherit;"><b> </b></span><br />
<span style="font-family: inherit;"><b> 6. Otras Animaciones</b>: en gris con distinto intervalo entre imágenes.</span><br />
<br />
<span style="font-family: "Courier New", Courier, monospace;"><span style="background-color: #cccccc;">gmt makecpt -Cdarkgray -T-75/0 > %color%</span></span><br />
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://2.bp.blogspot.com/-PdqkSikM5vk/XNyWgY9baBI/AAAAAAAATg8/xQV-k87pkUo7HOBU_ZOiLHeTKjJXi0fuQCEwYBhgL/s1600/3d_Moho_Gris_10.gif" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="993" data-original-width="1600" height="247" src="https://2.bp.blogspot.com/-PdqkSikM5vk/XNyWgY9baBI/AAAAAAAATg8/xQV-k87pkUo7HOBU_ZOiLHeTKjJXi0fuQCEwYBhgL/s400/3d_Moho_Gris_10.gif" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Animación con intervalos de 0.1 segundos.</td></tr>
</tbody></table>
<div class="separator" style="clear: both; text-align: center;">
</div>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://4.bp.blogspot.com/-aIn322DMx8Y/XNyWmMY2owI/AAAAAAAAThA/eynsUTnVYL8dzhBBNoxiITGVj75wz40zgCEwYBhgL/s1600/3d_Moho_Gris_15.gif" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="993" data-original-width="1600" height="247" src="https://4.bp.blogspot.com/-aIn322DMx8Y/XNyWmMY2owI/AAAAAAAAThA/eynsUTnVYL8dzhBBNoxiITGVj75wz40zgCEwYBhgL/s400/3d_Moho_Gris_15.gif" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Animación con intervalos de 0.15 segundos.</td><td class="tr-caption" style="text-align: center;"><br /></td><td class="tr-caption" style="text-align: center;"><br /></td><td class="tr-caption" style="text-align: center;"><br /></td></tr>
</tbody></table>
</div>
</div>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://4.bp.blogspot.com/-oYTovk1BKeE/XNyY-nRvVeI/AAAAAAAAThI/JMyoYpYqWRc1CPtNRGGtzOQFVPw9edRQACLcBGAs/s1600/3d_Moho_Gris_20.gif" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="993" data-original-width="1600" height="248" src="https://4.bp.blogspot.com/-oYTovk1BKeE/XNyY-nRvVeI/AAAAAAAAThI/JMyoYpYqWRc1CPtNRGGtzOQFVPw9edRQACLcBGAs/s400/3d_Moho_Gris_20.gif" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Animación con intervalos de 0.2 segundos.</td></tr>
</tbody></table>
</div>
Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0tag:blogger.com,1999:blog-1509523996050466325.post-22044600444287500882019-03-28T00:00:00.000-03:002019-08-08T23:41:16.015-03:002.3. Mapa Satelital. Blue Marble<div dir="ltr" style="text-align: left;" trbidi="on">
Mapas continentales a partir de imagen satelital (Blue Marble).<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
<img border="0" data-original-height="1409" data-original-width="1600" height="280" src="https://1.bp.blogspot.com/-hGiuZMoFyao/XJrl3qrM_II/AAAAAAAAS-c/ZuZCoVTdEbYgl-T0gpA4nhlfPvjIYpITACLcBGAs/s320/EJ2.3_Satelital_BlueMarble_W-65_4.png" width="320" /><a href="https://1.bp.blogspot.com/-8OQ5czYHGM8/XJcAZsjTIfI/AAAAAAAAS9Y/VOBWlVIAk_cMKvLVN1h7WSPuvSFtrCBrACPcBGAYYCw/s1600/EJ2.3_Satelital_BlueMarble_G-65_1.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1600" data-original-width="1600" height="320" src="https://1.bp.blogspot.com/-8OQ5czYHGM8/XJcAZsjTIfI/AAAAAAAAS9Y/VOBWlVIAk_cMKvLVN1h7WSPuvSFtrCBrACPcBGAYYCw/s320/EJ2.3_Satelital_BlueMarble_G-65_1.png" width="320" /></a></div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<a href="https://www.blogger.com/null" name="more"></a><br />
<a name='more'></a><br />
Se explica como hacer mapamundis o hemisféricos a partir de la imagen satelital Blue Marble: Next Generation de la NASA.<br />
<div dir="ltr" style="text-align: left;" trbidi="on">
<br />
<a href="https://www.dropbox.com/s/8e3g3scdbvyh9ws/EJ2.3_ImagenSatelital_BlueMarble.bat?dl=0">Script </a></div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<br /></div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<u><b>Datos</b></u>: <br />
<b>Imagen Satelital</b>: Utilizo la siguiente imagen de la NASA (<a href="https://neo.sci.gsfc.nasa.gov/view.php?datasetId=BlueMarbleNG-TB">link</a>) en formato geotiff de 0,1º de resolución. Esta imagen contiene datos en los polos y asigna colores azules a los océanos según la profundidad.<br />
<ul style="text-align: left;">
<li><a href="https://neo.sci.gsfc.nasa.gov/servlet/RenderData?si=526294&cs=rgb&format=TIFF&width=3600&height=1800">BlueMarbleNG-TB_2004-06-01_rgb_3600x1800.tiff</a></li>
</ul>
<br />
<b>1. Imagen a RGB</b>: Con mirrone se descompone la imagen en grillas RGB (ver 2.2 para más detalles).<br />
<ul style="text-align: left;">
<li><a href="https://www.dropbox.com/s/48h5k789q9u6ukb/BlueMarble_r.grd?dl=0">BlueMarble_r.grd</a></li>
<li><a href="https://www.dropbox.com/s/hbl4mt2qerppe19/BlueMarble_g.grd?dl=0">BlueMarble_g.grd</a></li>
<li><a href="https://www.dropbox.com/s/o49r1jxhnj0c5pn/BlueMarble_b.grd?dl=0">BlueMarble_b.grd</a></li>
</ul>
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<b>2. Mapas</b>: grdimage permite recombinar la grillas RGB para distintas proyecciones (-J) y regiones geográficas (-R). </div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">SET PROJ=W-65/15c</span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">SET REGION=-110/30/-90/20</span></span><br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt grdimage -R -J -O -K >> %OUT% "BlueMarble_r.grd" "BlueMarble_g.grd" "BlueMarble_b.grd" </span></span><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://4.bp.blogspot.com/-5WAnKaNnbQQ/XJwiaUkF5BI/AAAAAAAAS_A/Yf9qhKcQSbQD1en0YUbh182_Yk5LNxX0QCLcBGAs/s1600/EJ2.3_Satelital_BlueMarble_W-65_0.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1408" data-original-width="1600" height="351" src="https://4.bp.blogspot.com/-5WAnKaNnbQQ/XJwiaUkF5BI/AAAAAAAAS_A/Yf9qhKcQSbQD1en0YUbh182_Yk5LNxX0QCLcBGAs/s400/EJ2.3_Satelital_BlueMarble_W-65_0.png" width="400" /></a></div>
<br />
<b><br />
</b><b>3. Límites y marco:</b> se agregan los límites de los países, línea de costa y marco del mapa<br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;"></span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;"></span></span><span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;">gmt pscoast -R -J -O -K -N1/0.2,- >> %OUT%<br />
gmt pscoast -R -J -O -K -Df -W1/ >> %OUT%<br />
gmt psbasemap -R -J -O -K -B0 >> %OUT%</span></span><br />
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-hGiuZMoFyao/XJrl3qrM_II/AAAAAAAAS-c/ZuZCoVTdEbYgl-T0gpA4nhlfPvjIYpITACLcBGAs/s1600/EJ2.3_Satelital_BlueMarble_W-65_4.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1409" data-original-width="1600" height="350" src="https://1.bp.blogspot.com/-hGiuZMoFyao/XJrl3qrM_II/AAAAAAAAS-c/ZuZCoVTdEbYgl-T0gpA4nhlfPvjIYpITACLcBGAs/s400/EJ2.3_Satelital_BlueMarble_W-65_4.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Mollweide mostrando Sudamérica, Atlántico y parte de África). </td></tr>
</tbody></table>
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<b>4</b><b>A. Mapamundis</b>: Otros ejemplos de mapas con distintas proyecciones (-J) y regiones geográficas (-R).</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><br />
<a href="https://1.bp.blogspot.com/-j2GcA7kqk4s/XJVUr5leQjI/AAAAAAAAS8I/LCU9a9jaTwkgDx3tkTulOdlnzBwp3hGEwCEwYBhgL/s1600/EJ2.3_Satelital_BlueMarble_W.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="804" data-original-width="1600" height="200" src="https://1.bp.blogspot.com/-j2GcA7kqk4s/XJVUr5leQjI/AAAAAAAAS8I/LCU9a9jaTwkgDx3tkTulOdlnzBwp3hGEwCEwYBhgL/s400/EJ2.3_Satelital_BlueMarble_W.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Mollweide centrado en Greenwich (W).</td></tr>
</tbody></table>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-LjVJy3QUtvk/XJVUpRHC41I/AAAAAAAAS8A/mw1mw9I9WikXpY6nPTa38_I3zfbTfOmdgCEwYBhgL/s1600/EJ2.3_Satelital_BlueMarble_W-65.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="804" data-original-width="1600" height="200" src="https://1.bp.blogspot.com/-LjVJy3QUtvk/XJVUpRHC41I/AAAAAAAAS8A/mw1mw9I9WikXpY6nPTa38_I3zfbTfOmdgCEwYBhgL/s400/EJ2.3_Satelital_BlueMarble_W-65.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Mollweide centrado en Sudamérica (W-65).</td></tr>
</tbody></table>
<div class="separator" style="clear: both; text-align: center;">
</div>
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://4.bp.blogspot.com/-f2VlAJPEVcs/XJVVvUCygBI/AAAAAAAAS8Q/uExnU-6pTUsqibqvs0BZzNieCz6q8UkSACLcBGAs/s1600/EJ2.3_Satelital_BlueMarble_N-65.png" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="814" data-original-width="1600" height="201" src="https://4.bp.blogspot.com/-f2VlAJPEVcs/XJVVvUCygBI/AAAAAAAAS8Q/uExnU-6pTUsqibqvs0BZzNieCz6q8UkSACLcBGAs/s400/EJ2.3_Satelital_BlueMarble_N-65.png" width="400" /></a> </td><td style="text-align: center;"><br /></td><td style="text-align: center;"><br /></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Robinson centrado en Sudamérica (N-65).</td></tr>
</tbody></table>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><br />
<a href="https://3.bp.blogspot.com/-CAnzwd0vY7w/XJVYDdptSKI/AAAAAAAAS8c/d0ZXHI0aPUk4GA5jcfs2bdNCCCv8pHmHgCLcBGAs/s1600/EJ2.3_Satelital_BlueMarble_Ks-65.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="804" data-original-width="1600" height="200" src="https://3.bp.blogspot.com/-CAnzwd0vY7w/XJVYDdptSKI/AAAAAAAAS8c/d0ZXHI0aPUk4GA5jcfs2bdNCCCv8pHmHgCLcBGAs/s400/EJ2.3_Satelital_BlueMarble_Ks-65.png" width="400" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Eckert VI centrado en Sudamérica (Ks-65).</td></tr>
</tbody></table>
<br />
<b>4B. Mapas hemisféricos:</b><br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://3.bp.blogspot.com/-8OQ5czYHGM8/XJcAZsjTIfI/AAAAAAAAS9Q/gXEj72tqFXg_KDqdCLECWWn0H7zMl0CkACLcBGAs/s1600/EJ2.3_Satelital_BlueMarble_G-65_1.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1600" height="320" src="https://3.bp.blogspot.com/-8OQ5czYHGM8/XJcAZsjTIfI/AAAAAAAAS9Q/gXEj72tqFXg_KDqdCLECWWn0H7zMl0CkACLcBGAs/s320/EJ2.3_Satelital_BlueMarble_G-65_1.png" width="320" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Proyección estereográfica (-S-65/0/90)</td></tr>
</tbody></table>
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><br />
<a href="https://1.bp.blogspot.com/-D86Veuajy3E/XJcAZdpQYYI/AAAAAAAAS9M/UAdu-ShOXl0A678XGptmu2In-wWE768bACLcBGAs/s1600/EJ2.3_Satelital_BlueMarble_S0.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><span id="goog_140537764"></span><span id="goog_140537765"></span><img border="0" data-original-height="1600" data-original-width="1600" height="320" src="https://1.bp.blogspot.com/-D86Veuajy3E/XJcAZdpQYYI/AAAAAAAAS9M/UAdu-ShOXl0A678XGptmu2In-wWE768bACLcBGAs/s320/EJ2.3_Satelital_BlueMarble_S0.png" width="320" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Proyección ortográfica (-G-65/0/90)</td><td class="tr-caption" style="text-align: center;"><br /></td></tr>
</tbody></table>
</div>
</div>
<br />
<b>PD 08-08-2019</b>: Con la versión de GMT 6.0 es posible hacer mapas con proyección hemisférica NO centrada en el ecuador.<br />
<br />
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://1.bp.blogspot.com/-Ld2-iKsvvNI/XUzcz4u-III/AAAAAAAAUm4/ff662uMSeYUssLni29zAa12-1W1ECBQqgCLcBGAs/s1600/EJ2.3_Satelital_BlueMarble_G-65_2.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" data-original-height="1600" data-original-width="1600" height="320" src="https://1.bp.blogspot.com/-Ld2-iKsvvNI/XUzcz4u-III/AAAAAAAAUm4/ff662uMSeYUssLni29zAa12-1W1ECBQqgCLcBGAs/s320/EJ2.3_Satelital_BlueMarble_G-65_2.png" width="320" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">Proyección ortográfica (-G-65/-30/90)</td></tr>
</tbody></table>
<br />
<br />
<br />
<br />
<br /></div>
Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com3tag:blogger.com,1999:blog-1509523996050466325.post-52774404358616227532019-02-23T00:39:00.001-03:002019-02-23T00:39:25.117-03:004.2. Vectores y movimientos de placas<div dir="ltr" style="text-align: left;" trbidi="on">
<div class="separator" style="clear: both; text-align: center;">
</div>
<span id="goog_2039143337"></span><span id="goog_2039143338"></span>Se explica como dibujar vectores que indican el desplazamiento <u>relativo</u> entre las placas de Sudamérica y de Scotia. Se utilizan dos métodos: 1) de forma personalizada y 2) mediante el sitio <a href="https://www.unavco.org/">UNAVCO</a>.<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://4.bp.blogspot.com/-7xMI1HyYqdU/XHC5z_8dJoI/AAAAAAAASuM/JcUtJY02YSQtKOEY1vTdCNee7tFmEcQcgCLcBGAs/s1600/4.2_MapaVectores.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="763" data-original-width="1600" height="304" src="https://4.bp.blogspot.com/-7xMI1HyYqdU/XHC5z_8dJoI/AAAAAAAASuM/JcUtJY02YSQtKOEY1vTdCNee7tFmEcQcgCLcBGAs/s640/4.2_MapaVectores.png" width="640" /></a></div>
<br /><a name='more'></a><br />
<b><a href="https://www.dropbox.com/s/bb2efognh6fysdl/4.2_MapaVectores.bat?dl=0">Script</a> </b><br />
<br />
<b>0. Mapa previo</b>: Mapa previo a partir de una grilla con el límite transforme dibujado de Plates Project (ver 4.1).<br />
<br />
<span style="background-color: #cccccc;">gmt grdimage -R -J -O -K %CUT% -C%color% -I%SHADOW% >> %OUT%<br />
gmt psxy -R -J -O -K "transform.gmt" -A >> %OUT% -W0.4,30 -Sf1.5/0.25+l+s -Gblack</span><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-O1vNVbkAAL0/XHC5mdictaI/AAAAAAAASuU/sihLJMb-4tgWt9BKRjq5hwZgWgU5GhqIwCEwYBhgL/s1600/4.2_MapaVectores_0.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="763" data-original-width="1600" height="190" src="https://1.bp.blogspot.com/-O1vNVbkAAL0/XHC5mdictaI/AAAAAAAASuU/sihLJMb-4tgWt9BKRjq5hwZgWgU5GhqIwCEwYBhgL/s400/4.2_MapaVectores_0.png" width="400" /></a></div>
<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
</div>
<b>1. Vectores personalizados 1:</b> Los datos se ingresan manualmente.<b><br />
</b><br />
<div dir="ltr" style="text-align: left;" trbidi="on">
<b>1A. Datos</b>: longitud, latitud (del origen del vector), acimut (en grados desde el norte), largo del vector (en cm). Los datos se guardan en el archivo temp_vector.<br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">echo -68,-55,45,2c > temp_vector<br />
echo -68,-55,50,3c >> temp_vector</span></span><br />
<span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #999999;"><br />
</span></span> <br />
<b>1B. Graficar vectores</b>: con psxy se grafican con los siguientes argumentos:</div>
<ul style="text-align: left;">
<li>-SV+e: dibujar cabeza de vector de 0,3 cm de tamaño. </li>
<li>-W: propiedades de la línea.</li>
<li>-G: color del relleno de la cabeza del vector.</li>
</ul>
<div dir="ltr" style="text-align: left;" trbidi="on">
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;"><span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psxy -J -R -O -K >> %OUT% "temp_vector" -SV0.3+e -W0.4p,black -Gorange</span></span></span></span> </div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://2.bp.blogspot.com/-hJRMJKCFm3s/XHC5m9RphSI/AAAAAAAASuY/9O3Ua7hfXLkttn84I-OJQK2XAxTuc4FjgCEwYBhgL/s1600/4.2_MapaVectores_1.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="763" data-original-width="1600" height="190" src="https://2.bp.blogspot.com/-hJRMJKCFm3s/XHC5m9RphSI/AAAAAAAASuY/9O3Ua7hfXLkttn84I-OJQK2XAxTuc4FjgCEwYBhgL/s400/4.2_MapaVectores_1.png" width="400" /></a></div>
<br />
<br />
<b>2. Vectores personalizados 2: </b>echo y psxy se pueden combinar en único comando mediante una pleca. Con este método se pueden indicar los parámetros de cada vector por separado. Se utilizan los siguientes <a href="http://gmt.soest.hawaii.edu/doc/5.4.5/psxy.html#vector-attributes">modificadores de los vectores</a>: <br />
<ul style="text-align: left;">
<li>+bl: Mitad izquierda (<b>l</b>eft) del vector al inicio (<b>b</b>eginning) de la línea.</li>
<li>+er: Mitad derecha (<b>r</b>igth) del vector al final (<b>e</b>nd) de la línea. </li>
</ul>
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">echo -66,-53.5,85,4c | gmt psxy -J -R -O -K >> %OUT% -SV1+bl -W1.5p,red -Gred<br />
echo -66,-55.0,85,4c | gmt psxy -J -R -O -K >> %OUT% -SV1+er -W1.5p,red -Gred</span></span><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://3.bp.blogspot.com/-UXrC2z8g0vc/XHC6UzLUM6I/AAAAAAAASug/2Os-Qd4iuQIqD6NrdqwSd2oBKY8hPNHqACLcBGAs/s1600/4.2_MapaVectores_2.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="763" data-original-width="1600" height="190" src="https://3.bp.blogspot.com/-UXrC2z8g0vc/XHC6UzLUM6I/AAAAAAAASug/2Os-Qd4iuQIqD6NrdqwSd2oBKY8hPNHqACLcBGAs/s400/4.2_MapaVectores_2.png" width="400" /></a></div>
<br />
<br />
<b>3. Vectores UNAVCO: </b>Los datos se descargan desde el sitio web (ver Anexo 5).</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<br />
<b>3A. Datos:</b> Utilizo dos conjuntos de datos. <br />
<ul style="text-align: left;">
<li>Datos a partir de Modelos de Inversión (<a href="https://www.dropbox.com/s/kdxeir35b0lmkcy/UNAVCO_MODEL.txt?dl=0">UNAVCO_Model.txt</a>).</li>
<li>Datos a partir de Polos de Euler: (<a href="https://www.dropbox.com/s/n064k998vtya6hl/UNAVCO_EULER.txt?dl=0">UNAVCO_EULER.txt</a>).</li>
</ul>
<div class="separator" style="clear: both; text-align: center;">
<a href="http://1.bp.blogspot.com/-qg1Rb3lzVho/XGwQbXdpv3I/AAAAAAAASsQ/EdueRUzWZCcejgKXBgrHzpi0jeU_JpW7gCK4BGAYYCw/s1600/Captura.PNG" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="120" src="https://1.bp.blogspot.com/-qg1Rb3lzVho/XGwQbXdpv3I/AAAAAAAASsQ/EdueRUzWZCcejgKXBgrHzpi0jeU_JpW7gCK4BGAYYCw/s400/Captura.PNG" width="400" /></a></div>
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<b>3B. Graficar: </b>se utiliza <a href="http://gmt.soest.hawaii.edu/doc/5.4.5/supplements/meca/psvelo.html#psvelo">psvelo</a> con los siguientes argumentos:</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<ul style="text-align: left;">
<li>-W: atributos de la flecha. Por defecto línea negra, continua.</li>
<li>-Se: Selecciona el significado de las columnas de los datos de entrada. Para los de UNAVCO se utiliza e. Los valores indican la escala de los vectores (0,25), el límite de confidencia del 95% (0.95) y el tamaño de fuente para el nombre del sitio (0 para que NO anote el nombre del sitio). </li>
<li>-G: Color de la cabeza del vector.</li>
</ul>
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<span style="background-color: #cccccc;">gmt psvelo -J -R -O -K >> %OUT% "UNAVCO_MODEL.txt" -W -Se0.25/0.95/0 -Gwhite</span><br />
<span style="background-color: #cccccc;">gmt psvelo -J -R -O -K >> %OUT% "UNAVCO_EULER.txt" -W -Se0.25/0.95/0 -Ggreen</span><br />
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://4.bp.blogspot.com/-7xMI1HyYqdU/XHC5z_8dJoI/AAAAAAAASuc/yuSHv_UoXn46_ORBs_v5pIvgozgP7qT_wCEwYBhgL/s1600/4.2_MapaVectores.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="763" data-original-width="1600" height="190" src="https://4.bp.blogspot.com/-7xMI1HyYqdU/XHC5z_8dJoI/AAAAAAAASuc/yuSHv_UoXn46_ORBs_v5pIvgozgP7qT_wCEwYBhgL/s400/4.2_MapaVectores.png" width="400" /></a></div>
<br /></div>
<div dir="ltr" style="text-align: left;" trbidi="on">
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<br /></div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<b>Referencias</b>:<br />
UNAVCO Plate Motion Calculator (<a href="https://www.unavco.org/software/geodetic-utilities/plate-motion-calculator/plate-motion-calculator.html">link</a>).<br />
UNAVCO Plate Tectonics and GPS 'Reference Frames' (<a href="https://www.unavco.org/software/visualization/GPS-Velocity-Viewer/GPS-Velocity-Viewer-frames.html">link</a>).<br />
<div class="separator" style="clear: both; text-align: center;">
</div>
</div>
</div>
Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com2tag:blogger.com,1999:blog-1509523996050466325.post-43234366498247448462019-02-10T12:42:00.002-03:002019-02-20T00:23:09.137-03:00Anexo 5. Descargar datos de movimiento de placas (UNAVCO)<div dir="ltr" style="text-align: left;" trbidi="on">
El sitio UNAVCO tiene una herramienta <a href="https://www.unavco.org/software/geodetic-utilities/plate-motion-calculator/plate-motion-calculator.html#overview">online</a> (Plate Motion Calculator) que permite calcular el movimiento de placas tectónicas en cualquier punto de la Tierra. Los cálculos se pueden hacer a partir de 1) modelos inversión global y de 2) polos de Euler.<br />
<br />
<br />
<a name='more'></a><br />
<br />
<div dir="ltr" style="text-align: left;" trbidi="on">
Se explica como obtener los vectores de desplazamiento <u>relativo</u> de la Placa Scotia (placa móvil) con respecto a la Sudamericana (placa fija) en el punto 54 ºS y 63 ºO. Para ambos métodos, la forma de ingresar el sitio geográfico y el formato de salida es el mismo.<br />
<br />
<a href="https://www.unavco.org/software/geodetic-utilities/plate-motion-calculator/plate-motion-calculator.html#overview">Link calculadora </a><br />
<br />
<b>1. Sitio Geográfico: </b>Indicar longitud y latitud, altura del punto a analizar.<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://3.bp.blogspot.com/-cLLx8jGzVsk/WvWyTdKS4aI/AAAAAAAAQXE/8PDuPKygIUA5J6TOtMXG-SCWPWXnIpquwCLcBGAs/s1600/UNAVCO_Coordenadas.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="336" data-original-width="597" height="225" src="https://3.bp.blogspot.com/-cLLx8jGzVsk/WvWyTdKS4aI/AAAAAAAAQXE/8PDuPKygIUA5J6TOtMXG-SCWPWXnIpquwCLcBGAs/s400/UNAVCO_Coordenadas.png" width="400" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<br />
Alternativamente, se pueden ingresar múltiples sitios.<br />
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://2.bp.blogspot.com/-_YaN45l7IB8/XFhLED1mo2I/AAAAAAAASng/_xtivpiM86QchIDAb5L-ekr0XOMxZtYTwCEwYBhgL/s1600/Captura%2Bde%2Bpantalla%2Bde%2B2019-02-04%2B00-58-53.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="183" data-original-width="783" height="92" src="https://2.bp.blogspot.com/-_YaN45l7IB8/XFhLED1mo2I/AAAAAAAASng/_xtivpiM86QchIDAb5L-ekr0XOMxZtYTwCEwYBhgL/s400/Captura%2Bde%2Bpantalla%2Bde%2B2019-02-04%2B00-58-53.png" width="400" /></a></div>
<br />
<b>2. Formato del archivo de salida</b>: Elegir la opción <i>GMT <a href="http://gmt.soest.hawaii.edu/doc/5.4.5/supplements/meca/psvelo.html">psvelo</a></i>. El resultado se muestra en
una nueva ventana y se guarda en un archivo de texto plano.</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-l11HVNHAG2g/WvWyUNhjZ_I/AAAAAAAAQXk/qEp06KjpScUJEFIs5bSmD1A9hSTIr1DXwCEwYBhgL/s1600/UNAVCO_Submit.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="138" data-original-width="970" height="56" src="https://1.bp.blogspot.com/-l11HVNHAG2g/WvWyUNhjZ_I/AAAAAAAAQXk/qEp06KjpScUJEFIs5bSmD1A9hSTIr1DXwCEwYBhgL/s400/UNAVCO_Submit.png" width="400" /></a></div>
<br /></div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<br /></div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<b>3A. Modelos de inversión global</b>: El sitio cuenta con 18 modelos (a la fecha).<br />
<br />
<u>Nota</u>: Es importante verificar qué placas tectónicas se definen
para cada modelo, especialmente si se analizan placas
menores. Para este ejemplo sólo en los modelos NUVEL 1A, HS3-NUVEL1A, GSRM
v1.2 y v2.1, MORVEL están definidas las placas de Sudamérica (SA) y de
Scotia (ST). Por eso en el resto de los modelos los resultados son erróneos (ver 3C). <br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<span id="goog_1902270423"></span><span id="goog_1902270424"></span></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-Pir8RfY5-Qs/XFZoHwMcNcI/AAAAAAAASl8/Nv7MOrMKbgAs05YLGDzBq-po8WIp5PLrgCLcBGAs/s1600/Captura%2Bde%2Bpantalla%2Bde%2B2019-02-03%2B01-02-06.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="703" data-original-width="720" height="624" src="https://1.bp.blogspot.com/-Pir8RfY5-Qs/XFZoHwMcNcI/AAAAAAAASl8/Nv7MOrMKbgAs05YLGDzBq-po8WIp5PLrgCLcBGAs/s640/Captura%2Bde%2Bpantalla%2Bde%2B2019-02-03%2B01-02-06.png" width="640" /></a></div>
<br />
<b>3B. Elegir Modelo/s</b>: (1 en particular o todos los modelos), placa tectónica móvil
(Scotia en este ejemplo), y referencia (placa fija; Sudamérica).<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://2.bp.blogspot.com/-Df2nbKO4ZlY/WvWyT05QZMI/AAAAAAAAQXg/sDO-4RlybsQbgffh7dhD0bpx27vDCkb5gCEwYBhgL/s1600/UNAVCO_Placas.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="332" data-original-width="973" height="136" src="https://2.bp.blogspot.com/-Df2nbKO4ZlY/WvWyT05QZMI/AAAAAAAAQXg/sDO-4RlybsQbgffh7dhD0bpx27vDCkb5gCEwYBhgL/s400/UNAVCO_Placas.png" width="400" /></a></div>
<br />
<u>Nota</u>: Para obtener movimientos absolutos de las placas se debe elegir la opción NNR <i>no-net-rotation</i> en <i>Reference</i>.<br />
<br />
<div class="separator" style="clear: both; text-align: left;">
<b>3C. Datos de salida</b>: Conformado por longitud y latitud del sitio (columnas 1 y 2), velocidad hacia
el este y hacia el norte y sus respectivas incertezas
(columnas 3 a 6) de la velocidad, correlación entre los componentes (col 7). </div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://4.bp.blogspot.com/-SL6jVr-it0E/WvWyTRNfOCI/AAAAAAAAQXc/-F9zzR0KjnYbZK0SsZkiz0FY1smeH2F-wCEwYBhgL/s1600/UNAVCO_Output.png"></a><a href="https://4.bp.blogspot.com/-SL6jVr-it0E/WvWyTRNfOCI/AAAAAAAAQXc/-F9zzR0KjnYbZK0SsZkiz0FY1smeH2F-wCEwYBhgL/s1600/UNAVCO_Output.png" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="639" data-original-width="951" height="267" src="https://4.bp.blogspot.com/-SL6jVr-it0E/WvWyTRNfOCI/AAAAAAAAQXc/-F9zzR0KjnYbZK0SsZkiz0FY1smeH2F-wCEwYBhgL/s400/UNAVCO_Output.png" width="400" /> </a></div>
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<br />
<b>4A.
Polos de Euler:</b> La calculadora de UNAVCO también permite calcular el
movimiento relativo para el caso de que uno tenga el polos de Euler.<b> </b>Utilizo los de <a href="https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2007GL031699">Smalley 2007</a> quienes presentan el siguiente polo de Euler para el movimiento de Scotia con respecto a Sudamérica.<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://2.bp.blogspot.com/-9YNcZBJCZyY/XFedLDx6W5I/AAAAAAAASmw/ZmPm-X-ILRMeFGeiZGwIm0ODHagYq94mQCLcBGAs/s1600/Captura%2Bde%2Bpantalla%2Bde%2B2019-02-03%2B23-01-02.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="83" data-original-width="567" height="55" src="https://2.bp.blogspot.com/-9YNcZBJCZyY/XFedLDx6W5I/AAAAAAAASmw/ZmPm-X-ILRMeFGeiZGwIm0ODHagYq94mQCLcBGAs/s400/Captura%2Bde%2Bpantalla%2Bde%2B2019-02-03%2B23-01-02.png" width="400" /></a></div>
<br />
<br />
<b>4B. Ingreso del polo de Euler: </b>En las 2 primeras opciones se elige "user specified". En la tercera se ingresan los parámetros del polo (longitud, latitud y tasa de rotación en grados por <a href="https://es.wikipedia.org/wiki/Ma_(unidad_de_tiempo)">Ma</a>).<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://4.bp.blogspot.com/-QA50z536TD0/XFerF81RqoI/AAAAAAAASm8/5fq71OZJRLcRAKyOdBhJQZRZC07NVuT_gCLcBGAs/s1600/Captura%2Bde%2Bpantalla%2Bde%2B2019-02-03%2B01-53-24.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="588" data-original-width="780" height="301" src="https://4.bp.blogspot.com/-QA50z536TD0/XFerF81RqoI/AAAAAAAASm8/5fq71OZJRLcRAKyOdBhJQZRZC07NVuT_gCLcBGAs/s400/Captura%2Bde%2Bpantalla%2Bde%2B2019-02-03%2B01-53-24.png" width="400" /> </a></div>
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
<br />
<b>4C. Datos de salida:</b><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-QQhUoRyrK58/XFerHXVCMZI/AAAAAAAASnQ/XdB4guDaAV0CXzwntuABzJZtdfDpK4qOwCEwYBhgL/s1600/Captura%2Bde%2Bpantalla%2Bde%2B2019-02-03%2B01-54-21.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="111" data-original-width="831" height="50" src="https://1.bp.blogspot.com/-QQhUoRyrK58/XFerHXVCMZI/AAAAAAAASnQ/XdB4guDaAV0CXzwntuABzJZtdfDpK4qOwCEwYBhgL/s400/Captura%2Bde%2Bpantalla%2Bde%2B2019-02-03%2B01-54-21.png" width="400" /></a></div>
<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
</div>
</div>
<div dir="ltr" style="text-align: left;" trbidi="on">
</div>
Referencias:<br />
<ul style="text-align: left;">
<li><span style="font-size: x-small;">UNAVCO Plate Motion Calculator (<a href="https://www.unavco.org/software/geodetic-utilities/plate-motion-calculator/plate-motion-calculator.html">link</a>).</span></li>
<span style="font-size: x-small;">
</span>
<li><span style="font-size: x-small;">UNAVCO Plate Tectonics and GPS 'Reference Frames' (<a href="https://www.unavco.org/software/visualization/GPS-Velocity-Viewer/GPS-Velocity-Viewer-frames.html">link</a>).</span></li>
<span style="font-size: x-small;">
</span>
<li><span style="font-size: x-small;">
Smalley,
R.Jr.,
I. W. D. Dalziel,
M. G. Bevis,
E. Kendrick,
D. S. Stamps,
E. C. King,
F. W. Taylor,
E. Lauría,
A. Zakrajsek, and
H. Parra (<span class="pubYear">2007</span>), <span class="articleTitle">Scotia arc kinematics from GPS geodesy</span>, <i>Geophys. Res. Lett.</i>, <span class="vol">34</span>, L21308, doi:<a class="accessionId" href="https://doi.org/10.1029/2007GL031699" target="_blank" title="Link to external resource: 10.1029/2007GL031699">10.1029/2007GL031699</a>.
</span></li>
</ul>
</div>
Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0tag:blogger.com,1999:blog-1509523996050466325.post-28860736759244970712018-12-10T01:00:00.000-03:002018-12-13T13:28:15.500-03:0020.2. Bloque 3-D Columnas Apiladas<div dir="ltr" style="text-align: left;" trbidi="on">
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
Dibujar columnas 3-D en un mapa con perspectiva que indiquen la cantidad de 1º, 2º, 3º y 4º puesto de Copas Américas ganadas por país. <br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://4.bp.blogspot.com/-4HeE5WSrd78/XAWfAhLIbUI/AAAAAAAASI4/QeNuLT0Ci2cwaB41BwBHWeIYlLYRGRqxQCLcBGAs/s1600/Copa_America.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1386" data-original-width="1600" height="553" src="https://4.bp.blogspot.com/-4HeE5WSrd78/XAWfAhLIbUI/AAAAAAAASI4/QeNuLT0Ci2cwaB41BwBHWeIYlLYRGRqxQCLcBGAs/s640/Copa_America.png" width="640" /></a></div>
<br />
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<br />
<a name='more'></a><a href="https://www.blogger.com/null"></a>En el ejemplo anterior (<a href="https://mapasgeologicos.blogspot.com/2018/01/201-mapas-con-barras-3-d-copas-americas.html">20.1</a>) se muestra como dibujar columnas 3-D. Se explica como dibujar columnas apiladas en un mapa con perspectiva. Se dibujar 4 columnas correspondientes a la cantidad 1º, 2º, 3º y 4º puesto de las <a href="https://es.wikipedia.org/wiki/Copa_Am%C3%A9rica">Copa América</a>. <br />
<b>0. Datos</b>: Archivo <a href="https://es.wikipedia.org/wiki/Valores_separados_por_comas">csv</a>
conformado por 9 columnas: longitud, latitud, 1º puesto, 2º puesto, 3º
puesto, 4º puesto, 1º + 2º, 1º + 2 + 3, 1 +2 Copas Américas ganadas,
participaciones y país. La longitud y latitud corresponden al <a href="https://es.wikipedia.org/wiki/Centroide">centroide</a> de cada país (<a href="https://www.dropbox.com/s/gjcxoqn3edyb9c5/Centroids_Countrys.txt?dl=0">Centroids_Countrys.txt</a>, <a href="https://developers.google.com/public-data/docs/canonical/countries_csv">link original</a>) y definen la posición donde se dibujarán las columnas. La suma de posiciones (columnas 6, 7 y 8) se utilizan como base para apilar las columnas.<br />
<br />
<br />
Link: <a href="https://www.dropbox.com/s/v7zarglvb6v91up/CopaAmerica.csv?dl=0">CopaAmerica.csv</a><br />
<br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">REM Nombre archivo de entrada</span></span><br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">SET INPUT="CopaAmerica.csv"</span></span><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="http://2.bp.blogspot.com/-aaoBmHNAiYU/XAmjlnW8-sI/AAAAAAAASK0/5m68xfdfipUVgWRYxmfahkI0DzbNoPUOACK4BGAYYCw/s1600/Captura_Datos_CSV.PNG" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="165" src="https://2.bp.blogspot.com/-aaoBmHNAiYU/XAmjlnW8-sI/AAAAAAAASK0/5m68xfdfipUVgWRYxmfahkI0DzbNoPUOACK4BGAYYCw/s400/Captura_Datos_CSV.PNG" width="400" /></a></div>
<br />
<br />
<b>1. Mapa Base</b>: mapa en perspectiva con áreas húmedas y secas pintadas, límites de los países, título de la figura y eje vertical.<br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt pscoast -R -J -O -K -JZ -p >> %OUT% -Df -G200 </span></span><br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt pscoast -R -J -O -K -JZ -p >> %OUT% -Sdodgerblue2 -C200 -Df</span></span><br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt pscoast -R -J -JZ -O -K -p -Df >> %OUT% -N1</span></span><br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psbasemap -R -J -JZ -O -K -p >> %OUT% -BWSneZ -Ba</span></span><br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psbasemap -R -J -JZ -O -K -p >> %OUT% -B+t"Copas Am\351ricas Ganadas"</span></span><br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psbasemap -R -J -JZ -O -K -p >> %OUT% -BWSneZ+b -Bzafg+l"Cantidad"</span></span><br />
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://2.bp.blogspot.com/-qE71ZSmnGGU/W_4S4FPbWgI/AAAAAAAASGI/hUhHkmDuZ2IxbZcUn4z2IPL3Tp7rqWoZwCEwYBhgL/s1600/Copa_America_v2_0.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1386" data-original-width="1600" height="277" src="https://2.bp.blogspot.com/-qE71ZSmnGGU/W_4S4FPbWgI/AAAAAAAASGI/hUhHkmDuZ2IxbZcUn4z2IPL3Tp7rqWoZwCEwYBhgL/s320/Copa_America_v2_0.png" width="320" /> </a> </div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
<b>2. Columnas apiladas</b>: con <a href="http://gmt.soest.hawaii.edu/doc/5.4.4/psxyz.html">psxyz</a> se dibujan las columnas centradas en las coordenadas (columnas 0 y 1) que se extiende desde una base hasta una altura (valores 3 y 4) con los parámetros:</div>
<ul style="text-align: left;">
<li><a href="http://gmt.soest.hawaii.edu/doc/5.4.2/psxyz.html#s">-S</a>o0.5c: dibuja columnas en 3-D con una base cuadrada de 0,5 cm de lado.</li>
<li>-G: color de la columna.</li>
<li>-W: propiedades de las líneas.</li>
<li>-i0,1,6,2: columnas con los datos de entrada (longitud, latitud, valor base, valor altura).</li>
</ul>
<div class="separator" style="clear: both; text-align: left;">
<b>a. 1º puesto</b>: Columnas 0, 1, 2. En este caso, al omitirse el valor de la base se asume que corresponde a 0. </div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psxyz -R -J -JZ -O -K -p >> %OUT% %INPUT% -So0.5c -Gblue -Wthinner -i0,1,2</span></span></span></div>
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://2.bp.blogspot.com/-wHh38Uo4IL8/W_4S4GBaUTI/AAAAAAAASGM/AlPVbgpIiYIJBL3_uVfBG3svRvrj4FtUQCEwYBhgL/s1600/Copa_America_v2_1.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1386" data-original-width="1600" height="277" src="https://2.bp.blogspot.com/-wHh38Uo4IL8/W_4S4GBaUTI/AAAAAAAASGM/AlPVbgpIiYIJBL3_uVfBG3svRvrj4FtUQCEwYBhgL/s320/Copa_America_v2_1.png" width="320" /> </a></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
<b>b. 2º Puesto</b>: Columnas 0, 1, 6 (1º + 2º puestos) y 2 (1º puestos).</div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psxyz -R -J -JZ -O -K -p >> %OUT% %INPUT% -So0.5cb -Ggreen -Wthinner -i0,1,6,2</span></span> </div>
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://2.bp.blogspot.com/-IFhuM8poeVs/W_4S5MMs47I/AAAAAAAASGQ/xaMswOWnjzo2fXFM4rlf8yzOi-ltfLxSQCEwYBhgL/s1600/Copa_America_v2_2.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1386" data-original-width="1600" height="277" src="https://2.bp.blogspot.com/-IFhuM8poeVs/W_4S5MMs47I/AAAAAAAASGQ/xaMswOWnjzo2fXFM4rlf8yzOi-ltfLxSQCEwYBhgL/s320/Copa_America_v2_2.png" width="320" /></a></div>
<span style="font-family: "courier new" , "courier" , monospace;"> </span><br />
<br />
<b>c. 3º puesto</b>: Columnas 0, 1, 7 (1º + 2º + 3º puestos) y 6 (1º + 2º puestos).<br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psxyz -R -J -JZ -O -K -p >> %OUT% %INPUT% -So0.5cb -Gyellow -Wthinner -i0,1,6,7</span></span><br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://4.bp.blogspot.com/-vyjhXf_pNew/W_4S5e89weI/AAAAAAAASGU/QbdmaFDgaqoouI2i9PdcstosgpGR4oOfwCEwYBhgL/s1600/Copa_America_v2_3.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1386" data-original-width="1600" height="277" src="https://4.bp.blogspot.com/-vyjhXf_pNew/W_4S5e89weI/AAAAAAAASGU/QbdmaFDgaqoouI2i9PdcstosgpGR4oOfwCEwYBhgL/s320/Copa_America_v2_3.png" width="320" /></a></div>
<br />
<b>d. 4º puesto</b>: Columnas 0, 1, 8 (1º + 2º + 3º + 4º puestos) y 7 (1º + 2º + 3º puestos) .<br />
<br />
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt psxyz -R -J -JZ -O -K -p >> %OUT% %INPUT% -So0.5cb -Gred -Wthinner -i0,1,8,7</span></span><br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://1.bp.blogspot.com/-BWJ72eCR48A/W_4S5mQKpCI/AAAAAAAASGY/DvETgZ_lL-w8RBAEjNus6rSyWQQ6dZbQQCEwYBhgL/s1600/Copa_America_v2_4.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1386" data-original-width="1600" height="277" src="https://1.bp.blogspot.com/-BWJ72eCR48A/W_4S5mQKpCI/AAAAAAAASGY/DvETgZ_lL-w8RBAEjNus6rSyWQQ6dZbQQCEwYBhgL/s320/Copa_America_v2_4.png" width="320" /></a></div>
<b><br /></b>
<b>3. Agregar valores</b>:<span style="font-family: "courier new" , "courier" , monospace;"> </span>Para escribir en la base la suma de los 4 primeros puestos de campeonatos en la base de las columnas, primero se utiliza <a href="http://gmt.soest.hawaii.edu/doc/5.4.2/gmtconvert.html">gmtconvert</a> para extraer los datos de las columnas 0, 1 y 8. Estos se pasan a <a href="http://gmt.soest.hawaii.edu/doc/5.4.2/pstext.html">pstext</a> (mediante el símbolo pleca) con los siguientes argumentos:<br />
<ul style="text-align: left;">
<li>-Gwhite@30: fondo blanco con 30% de transparencia.</li>
<li>-D0/-0.8c: Desplaza el texto 8 mm hacia abajo (con respecto a la posición geográfica definidas por las columnas 0 y 1).</li>
<li>-F: Fuente del texto (tamaño, tipo y borde de letra, justificación).<br /> </li>
</ul>
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">gmt convert %INPUT% -o0,1,8 | gmt pstext -R%REGION% -J -O -K -p -Gwhite@30 -D0/-0.8c -F+f20p,Helvetica-Bold,firebrick=thinner+jCM >> %OUT% </span></span><br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://3.bp.blogspot.com/-mN2Rvl8o6eQ/XAWfAuowZDI/AAAAAAAASI0/fhz_89CfGTY25Q-1MoCJkrP6oVGu9s67QCEwYBhgL/s1600/Copa_America_5.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1376" data-original-width="1600" height="275" src="https://3.bp.blogspot.com/-mN2Rvl8o6eQ/XAWfAuowZDI/AAAAAAAASI0/fhz_89CfGTY25Q-1MoCJkrP6oVGu9s67QCEwYBhgL/s320/Copa_America_5.png" width="320" /></a></div>
<br />
<br />
<div class="separator" style="clear: both; text-align: left;">
<b>4. Leyenda</b>: Con echo primero se crea el archivo auxiliar que luego se dibuja con <a href="http://gmt.soest.hawaii.edu/doc/5.4.4/pslegend.html">pslegend</a>.</div>
<div class="separator" style="clear: both; text-align: left;">
</div>
<div class="separator" style="clear: both; text-align: left;">
</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">echo H 10 Times-Roman Copas Am\351ricas > "temp_leyenda"<br />echo N 1 >> "temp_leyenda"<br />echo S 0.25c r 0.5c blue 0.25p 0.75c 1 Puesto >> "temp_leyenda"<br />echo S 0.25c r 0.5c green 0.25p 0.75c 2 Puesto >> "temp_leyenda"</span></span><br />
<div class="separator" style="clear: both; text-align: left;">
<span style="background-color: #cccccc;"><span style="font-family: "courier new" , "courier" , monospace;">echo S 0.25c r 0.5c yellow 0.25p 0.75c 3 Puesto >> "temp_leyenda"<br />echo S 0.25c r 0.5c red 0.25p 0.75c 4 Puesto >> "temp_leyenda"<br /><br />gmt pslegend -R -J -JZ -O -K -p >> %OUT% "temp_leyenda" -DjLB+o0.2i+w1.35i/0+jBL -F+glightgrey+pthinner+s-4p/-6p/grey20@40</span></span></div>
<div class="separator" style="clear: both; text-align: center;">
</div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://2.bp.blogspot.com/-4HeE5WSrd78/XAWfAhLIbUI/AAAAAAAASI4/qbS-rnE7jTQl36m4t3rAMlUHmKddh3XVACEwYBhgL/s1600/Copa_America.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" data-original-height="1386" data-original-width="1600" height="277" src="https://2.bp.blogspot.com/-4HeE5WSrd78/XAWfAhLIbUI/AAAAAAAASI4/qbS-rnE7jTQl36m4t3rAMlUHmKddh3XVACEwYBhgL/s320/Copa_America.png" width="320" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<br /></div>
Federico Estebanhttp://www.blogger.com/profile/09230004561888756222noreply@blogger.com0