clear; clc; close all;

if 1 == 2
    filename='Bild2';
    endung='.png';
    foldername='C:\Users\y1swun\Desktop\Spektrum_aus_Bild\';
    x_min=0;
    x_max=3000;
    y_min=0.0;
    y_max=120;
    picture=imread(strcat(foldername,filename,endung));
    picture(1,:,:) = 255;
    picture(end,:,:) = 255;
    picture(:,1,:) = 255;
    picture(:,end,:) = 255;
    for ii = 1:size(picture,1)
        ids = find((picture(ii,:,1)==0)&(picture(ii,:,2)==255)&(picture(ii,:,3)==255));
        picture(ii,ids,:) = 0;
    end
    [spectrum] = f_spectrum(picture,x_min,x_max,y_min,y_max);
    spectrum(:,2) = spectrum(:,2).*1e-4.*spectrum(:,1).^2;
    spectrum(:,1) = 1e4./spectrum(:,1);
    figure(); plot(spectrum(:,1),spectrum(:,2)); xlim([3,200])
    spectrum(1:6,:) = [];
    spectrum(end-5:end,:) = [];
    
    save(strcat(foldername,'Turner2012','.txt'), 'spectrum', '-ASCII', '-tabs');
end

spectrum = dlmread('Turner2012.txt','',0,0);


h=6.626e-34;
c=3e8;
k=1.3806e-23;
c1 = 1e-10*2*pi*h*c^2;
c2 = h*c/k;
T = 280;
L = spectrum(:,1)*1e-6;

I = c1./( L.^5 .* ( exp(c2./L/T) -1 ) );


f = figure(); hold on;
plot(spectrum(:,1),3e6*I,'r-','LineWidth',2)
plot(spectrum(:,1),spectrum(:,2),'k','LineWidth',1.5)
xlim([0,50])
xlabel('Wellenlänge [µm]')
ylabel('Strahlungsfluss [mW/m²/sr/µm]')
legend(sprintf('schwarzer Strahler %.0f K (%.0f °C)',T,T-273.15),'Erde');
grid on;
set(f,'PaperPositionMode','auto')
print(f,'Turner2012.png','-dpng','-r150');