summary sma

From DavinciWiki
Jump to: navigation, search

Description

Create a summary table of results for one pixel from sma output.

Arguments and Return Values

Arguments: takes a standard sma_output structure and a position in the mixed cube

Return Value: a table displaying the grouped concentrations

Usage

Syntax: summary_sma(sma_output, x, y [, xlow=INT][, xhigh=INT][ ,ylow=INT][, yhigh=INT][, group=BOOL][, output=BOOL])

sma_output = output structure from sma

x = x-coordinate

y = y-coordinate

xlow=low value for x-range

xhigh=high value for x-range

ylow=low value for y-range

yhigh=high value for y-range

group=use the .grouped component of the sma_output structure

output=return a structure with the sorted abundances, error, and endmemeber names

Examples

dv> final=sma(mixspec,library,wave=250,wave2=1305,exclude=2//9//44, group=1)
Division by zero, 130049 times
struct, 16 elements                                                
   forcedlib: struct, 1 elements                                     
       label: "none"     
   excluded: 3x1x1 array of int, bsq format [12 bytes]               
   bb: 180x70x1 array of float, bsq format [50,400 bytes]            
   bberror: 180x70x1 array of float, bsq format [50,400 bytes]       
   normconc: 180x70x47 array of float, bsq format [2,368,800 bytes]  
   measured: 180x70x73 array of float, bsq format [3,679,200 bytes]  
   modeled: 180x70x73 array of float, bsq format [3,679,200 bytes]   
   conc: 180x70x50 array of float, bsq format [2,520,000 bytes]      
   rms: 180x70x1 array of float, bsq format [50,400 bytes]           
   nsamples: 71                                                      
   endlib: struct, 4 elements                                       
       data: 44x1x73 array of float, bsq format [12,848 bytes]
       group: Text Buffer with 44 lines of text
           1: Quartz
           2: K-Feldspar
           3: Plagioclase
           4: Plagioclase
           5: Plagioclase
           6: Plagioclase
           7: Plagioclase
           8: Plagioclase
           9: Orthopyroxene
           10: Orthopyroxene
       label: Text Buffer with 44 lines of text
           1: Quartz BUR-4120
           2: Microcline BUR-3460
           3: Albite WAR-0235
           4: Oligoclase BUR-060D
           5: Andesine WAR-0024
           6: Labradorite BUR-3080A
           7: Bytownite WAR-1384
           8: Anorthite BUR-340
           9: Bronzite NMNH-93527
           10: Enstatite HS-9.4B
       xaxis: 1x1x73 array of double, bsq format [584 bytes]
   error: 180x70x50 array of float, bsq format [2,520,000 bytes]     
   grouped: struct, 3 elements                                          
       grouped_conc: 180x70x20 array of float, bsq format [1,008,000 bytes]  
       grouped_error: 180x70x20 array of float, bsq format [1,008,000 bytes]
       grouped_labels: Text Buffer with 20 lines of text
           1: Quartz
           2: K-Feldspar
           3: Plagioclase
           4: Orthopyroxene
           5: Pigeonite
           6: Clinopyroxene
           7: Olivine
           8: Mica
           9: Serpentine
           10: Sheet-silicate
   spectral_range: "Channels 3-73"                                   
   wave1: 250                                                        
   wave2: 1305                                                       

Example: show the modeled concentrations for pixel [100,32] of mixspec
dv> summary_sma(final,100,32)
*********************Summary for Pixel Coordinate (1,1) ******************* 
              Endmember             Abundance(%)       Normalized for BB(%) 
           Quartz BUR-4120        0.33 +/-  0.34          0.62 +/-  0.63
           Albite WAR-0235        1.66 +/-  0.88          3.10 +/-  1.65
     Labradorite BUR-3080A        1.26 +/-  4.52          2.36 +/-  8.46
        Bytownite WAR-1384       13.76 +/-  3.87         25.75 +/-  7.23
         Enstatite HS-9.4B        0.57 +/-  1.30          1.06 +/-  2.43
   Hypersthene NMNH-B18247        0.82 +/-  3.33          1.54 +/-  6.22
        Augite NMHN-122302       15.94 +/-  2.57         29.82 +/-  4.80
              KI 3362 Fo60        0.84 +/-  1.29          1.57 +/-  2.41
              KI 3373 Fo35        0.62 +/-  1.89          1.16 +/-  3.53
         saponite <0.2 mic        1.01 +/-  0.97          1.88 +/-  1.81
           Swy-1 < 0.2 mic        1.58 +/-  3.03          2.95 +/-  5.66
             02-011 Opal A        0.38 +/-  0.98          0.72 +/-  1.83
 Magnesiohornblende WAR-03        1.94 +/-  0.79          3.63 +/-  1.47
              Anhydrite S9        2.30 +/-  0.29          4.31 +/-  0.54
    Gypsum (Satin spar) S6        1.74 +/-  0.73          3.26 +/-  1.36
               Calcite C40        2.18 +/-  0.22          4.08 +/-  0.40
              Dolomite C20        2.21 +/-  0.22          4.13 +/-  0.42
 Crystalline stilbite (zeo        4.31 +/-  1.47          8.06 +/-  2.75
___________________________________________________________________________
                  Sum (includes bb) =  100.22
                Blackbody Abundance =   46.76 +/-  2.66
                       RMS Error(%) =   0.127584 

DavinciWiki Mini-Nav Bar

Contents


Contact Developers

  • davinci-dev [AT] mars.asu.edu

All other topics

  • See navigation on the left

Major Sub-Functions

Related Functions

Recent Library Changes

Created On: 01-30-2008
Modified On: 02-03-2011

Personal tools