function ndsi, ref_065, ref_160, missing

   ndsi = (ref_065 - ref_160) / (ref_065 + ref_160)

   idx = where(ref_065 eq missing or ref_160 eq missing,cc)
   if (cc gt 0) then begin
      ndsi[idx] = missing
   endif

   return, ndsi
end