Davinci example import and process ASTER Level 1B imagery

gistfile1.txt

Davinci http://davinci.asu.edu/ is a NASA-funded flight code software package used as a data processing tool, developed at Arizona State University by the Mars Space Flight Facility development team. Davinci providing powerful numeric data manipulation and spectral analytical tools and has a robust suite of tools designed to mosaic, normalize, register, and blend large datasets . The core of Davinci is a free, open-source C-like interactive and scripting language that relies on a math engine developed  in C.  Davinci has a variety of  vector-oriented features that make working with large, complicated datasets much easier.  Davinci is available as a prebuilt binary for all major operating systems (Windows, Mac, and Linux).

This is an example of bringing in ASTER L1B data and processing through Davinci. The data can then be displayed in JMARS - a Java-based GIS optimised for imagery.

Follow the directions to get Davinci installed and running first: http://davinci.asu.edu/index.php?title=Download_Davinci

Access the ASTER L1B file here (Must have ASU login)
https://drive.google.com/a/asu.edu/file/d/0B0KiGz16sLxrUnpINWt2LThWelU/view?usp=sharing

In Davinci:

  a=load_aster("AST_L1B_00308072004185055_20140212133605_5109.hdf")
  IR=a.data.ir
  IR=deplaid_aster(IR)
  IR2 = rrmnoise_pca(IR,b10=0)
  #make a decorrelation stretch of everything
  dcs=do_dcs(IR2,4,3,2,sample=IR)
  #display the dcs  display(dcs)
  display(dcs)
  
  #output to JMARS to look at with other data in a GIS
  jmars(dcs,a,name=name+"_L1B",id=name,tmpdir="/Users/lela/Desktop/davinci/png/");