anchor.diff:         call error("image_setup", "unexpected value of io_bytes (should be 4 or 8)")
       end select
   
  +    ! Set up filters
  +    if(mp_exists_keyword(handle, path, 'use_filters')) then
  +       call mp_read_keyword(handle, path, 'use_filters',img%use_filters)
  +       if(img%use_filters) then
  +          if(img%use_exact_nu) call error("image_setup", "cannot use filters in monochromatic mode")
  +          ! TODO: also don't allow if using raytracing
  +          ! TODO: set up filters here
  +       end if
  +    else
  +       img%use_filters = .false.
  +    end if
  +
  +
     end subroutine image_setup
   
     subroutine find_sed_bin(img,x_image,y_image,ir)
  
positive.code:        if(img%compute_sed) then
            call find_sed_bin(img,x_image,y_image,ir)
            if(ir >= 1 .and. ir <= img%n_ap) then
               img%sed(inu,ir,im,io,:) = img%sed(inu,ir,im,io,:) + img%tmp_stokes * p%energy
               if(img%uncertainties) then
                  img%sed2(inu,ir,im,io,:) = img%sed2(inu,ir,im,io,:) + img%tmp_stokes**2._dp * p%energy**2._dp
                  img%sedn(inu,ir,im,io,:) = img%sedn(inu,ir,im,io,:) + 1._dp
               end if
            end if
         end if
      end if
  
    end subroutine image_bin
  
    subroutine image_bin_raytraced(img,p,x_image,y_image,im,spectrum)
  
  
negative.code:      real(dp), allocatable :: log10_emissivity(:,:)
    end type dust_helper
  
    type filter
  
       real(dp) :: n_wav
       real(dp), allocatable :: nu(:), tr(:)
  
    end type filter
  
    type image
  
       ! IMAGE PARAMETERS
  