Collect several satellite images from the INPE/CPTEC collection
1 次查看(过去 30 天)
显示 更早的评论
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 个评论
采纳的回答
Mathieu NOE
2022-11-14
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 个评论
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Reference Applications 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!