how to convert a netcdf file into a tiff/geotiff file

i want to convert netcdf files to geotiff in matlab.any idea on how to go about it.

Answer 1

[Reyadh Albarakat]

I recommend you to use "geotiffwrite". https://www.mathworks.com/help/map/ref/geotiffwrite.html.

[Sir Felix Otieno]

ncfiles = dir('*.nc') ;
N = length(ncfiles) ;
for i = 1:N
    ncfile = ncfiles(i).name ;
      lat = ncread(ncfile,'lat') ;
      lon = ncread(ncfile,'lon') ;
      %%Change lon and lat to double
      lon = double(lon) ;
      lat = double(lat) ;
      A = ncread(ncfile,'MYD08_D3_6_AOD_550_Dark_Target_Deep_Blue_Combined_Mean') ;
      % Write the .grd data into geotiff
      R = georasterref('RasterSize',size(A),'LatitudeLimits',[min(lat),max(lat)],........
          'LongitudeLimits',[min(lon),max(lon)]);
      tiffile = strcat(ncfile,'.tif') ;
      geotiffwrite(tiffile,A,R)
  end

i have tried to use this sample but itsnot working out for me

[Reyadh Albarakat]

Dear;

I wrote this code for you. Please try it!

    % 
OutFolder = 'Your Outputs Path';
cd Your Data Path
dinfo = dir('*.nc');
nfile = length(dinfo);
filenames = {dinfo.name};

for k = 1:nfile
  file_name{k} = filenames{k};
  lat = ncread(file_name{i},'lat') ;
  lon = ncread(file_name{i},'lon') ; 
  A{k} = ncread(file_name{k},'MYD08_D3_6_AOD_550_Dark_Target_Deep_Blue_Combined_Mean');
  A{k} = (double(A{k}));
  R = georasterref('RasterSize',size(A{k}),'LatitudeLimits',[min(lat),max(lat)],'LongitudeLimits',[min(lon),max(lon)]);
  tiffile{k} = strcat(filenames{k},'.tif') ;
  geotiffwrite(fullfile(OutFolder,tiffile{k}),A{k},R);
end

[Sir Felix Otieno]

OutFolder = 'C:\Users\Student\Desktop\phelix\SOFTWARES'; cd C:\Users\Student\Desktop\phelix\SOFTWARES dinfo = dir('*.nc'); nfile = length(dinfo); filenames = {dinfo.name}; for k = 1:nfile file_name{k} = filenames{k};

lat = ncread(file_name{i},'lat') ;
lon = ncread(file_name{i},'lon') ; 

A{k} = ncread(file_name{k},'ch1');
A{k} = (double(A{k}));
R = georasterref('RasterSize',size(A{k}),'LatitudeLimits',[min(lat),max(lat)],'LongitudeLimits',[min(lon),max(lon)]);
tiffile{k} = strcat(filenames{k},'.tif') ;
geotiffwrite(fullfile(OutFolder,tiffile{k}),A{k},R);
end

Subscript indices must either be real positive integers or logicals.

Error in reference (line 9) lat = ncread(file_name{i},'lat') ; i am getting that error. you give me your e mail i send you the file i am using

Answer 2

[Sir Felix Otieno]

i have tried to use geotiffwrite but my files are in netcdf format.i dont know how i can go about that