Hace un par de meses comencé con el análisis de datos geoespaciales.

En agosto de 2017 fui a un congreso llamado “Agente de cambio” donde se habló, entre otros temas, de la industria aeroespacial, tuve la oportunidad de tener una conferencia mediante Skype con una astronauta mexicana de la NASA que se encontraba en Polonia en una misión análoga a Marte, también tuve el gusto de conocer a Guillermo Miguel del Castillo Hoffman, estudiante de ingeniería en computación de la UAM Azcapotzalco y que en 2016 fue seleccionado por la NASA para realizar una estancia en el Ames Research Center, sin hacer preámbulo, lo que Hoffman planteó al final de su plática me motivó a realizar análisis de información satelital, Hoffman dijo; “En la NASA hace falta que se analicen grandes volúmenes de datos satelitales”, así que al día siguiente del congreso empecé a investigar cómo podría realizar consultas a los satélites de la NASA mediante el lenguaje estadístico R, encontré un sitio en internet llamado SkyWatch en la cual tuve que pedir al administrador del sitio una llave única para que pudiera descargar imágenes satelitales y datasets climáticos y atmosféricos, me otorgaron la llave en menos de 24 horas, entonces mediante Microsoft R services realicé mi primer consulta a un satélite de la NASA (fui muy feliz) llamado Landsat-8, el código de la consulta es el siguiente:

setwd("C:/Users/DataWarehouse/Desktop/Data_Science/skywatch")

library(SkyWatchr)
library(httr)
library(htmlTable)
library(sp)

api_key <- readLines('API_KEY_SKYWATCH.txt')
options(SkyWatchr.apikey = api_key)
geodata <- querySW(api_key = api_key, time_period = "2015-09", longitude_latitude = "-114.3651,43.6827", instrument_satellite = "Landsat-8", wavelength_band = "coastal aerosol")

 

 

Anuncios