Collect several satellite images from the INPE/CPTEC collection
    5 Ansichten (letzte 30 Tage)
  
       Ältere Kommentare anzeigen
    
    Augusto Gabriel da Costa Pereira
 am 11 Nov. 2022
  
    
    
    
    
    Kommentiert: Mathieu NOE
      
 am 15 Nov. 2022
            I have this script below to download ONLY ONE satellite image from channel 11 (infrared) for 11/10/2022 at 19:20 GMT.
QUESTION 1
But what I really want now is to download ALL of his imagesfrom 00:00 GMT to 23:50 GMT.
QUESTION 2
Another code I want is to specify the times, being for the same day but from 00:00 UTC to 02:00 UTC.
% Website below with the collection of images from INPE/CPTEC
% http://satelite.cptec.inpe.br/acervo/goes16.formulario.logic 
% The script below downloads only an image from the INPE/CPTEC collection
% Specify the filename (URL).
filename = 'http://satelite.cptec.inpe.br/repositoriogoes/goes16/goes16_web/ams_ret_ch11_baixa/2022/11/S11635384_202211011920.jpg';
% Read it in to a variable in my m-file.
rgbImage = imread(filename);
% Display the image.
imshow(rgbImage);
% Save the image to a file on my local hard drive.
imwrite(rgbImage, 'satellite image.jpg'); % imwrite(rgbImage, 'satellite image.png')
Best regards,
AP
0 Kommentare
Akzeptierte Antwort
  Mathieu NOE
      
 am 14 Nov. 2022
        hello Augusto
try this code 
as far as I could find out there are pictures available with a time delta of 10 minutes 
choose the option to get the  result for question Q1 or Q2 : 
%% options 
option = 1 % :  to download ALL of his images from 00:00 GMT to 23:50 GMT
% option = 2 % :  to download  images  from 00:00 UTC to 02:00 UTC 
all the best 
full code : 
% GMT stands for Greenwich Mean Time. UTC is known as Universal Time. UTC is 0 hours ahead of GMT.
% So, when it is 11:00 am GMT  it will be  11:00 am UTC
% in short : Coordinated Universal Time is the same as in Greenwich Mean Time
% date and time intervals
day_date = '20221101';
% images are available by 10 minutes interval
%% options 
option = 1 % :  to download ALL of his images from 00:00 GMT to 23:50 GMT
% option = 2 % :  to download  images  from 00:00 UTC to 02:00 UTC 
httpsUrl = "http://satelite.cptec.inpe.br/repositoriogoes/goes16/goes16_web/ams_ret_ch11_baixa/2022/11/S11635384_";
%% loop every 10 minutes name
mm = 00; % start minutes at 00
hh = 00; % start hours at 00
iter = 0;
while 1    % indeed there are 145 pictures distant of 10 minustes in one day (from 00:00 GMT to 23:50 GMT)
    iter = iter + 1;
    if mm == 60
        mm = 00;
        hh = hh+1; % update hours 
    end
    %% option 1 
    if option == 1 && hh == 24
        break
    end
    %% option 2 :  to download  images  from 00:00 UTC to 02:00 UTC 
    if option == 2 && hh == 02 && mm == 10
        break
    end
    if hh<10
        str_hh = ['0' num2str(hh)];
    else
        str_hh = num2str(hh);
    end
    if mm<10
        str_mm = ['0' num2str(mm)];
    else
        str_mm = num2str(mm);
    end
    hh_mm_str = [str_hh str_mm];
    mm = mm+10; % update minutes 
    %%%%%%%%%%%%%%%%%%%% main code %%%%%%%%%%%%%%%%%%
    filename = ['S11635384_' day_date hh_mm_str '.jpg'];
    disp(strcat(filename, "  is processed"));
    filename_url = ['http://satelite.cptec.inpe.br/repositoriogoes/goes16/goes16_web/ams_ret_ch11_baixa/2022/11/' filename];
    % Read it in to a variable in my m-file.
    rgbImage = imread(filename_url);
%     % Display the image.
%     imshow(rgbImage);
    % Save the image to a file on my local hard drive.
    imwrite(rgbImage,filename); 
end
2 Kommentare
Weitere Antworten (0)
Siehe auch
Kategorien
				Mehr zu Reference Applications finden Sie in Help Center und File Exchange
			
	Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!

