NetCDF into GeoTIFF/TIFF file conversion

I have lots of NetCDF files which I have to first read in MATLAB and then convert into TIFF files. As I am totally new to MATLAB kindly help.

[José-Luis]

If you are totally new to Matlab, why use it for this? It is a rather expensive tool if your problem is limited to that. There are better (and free) tools for manipulating geographic data anyway.

GDAL comes to mind.

Answer 1

[KSSV]

First you need to read the netCDF file.....for this check ncdisp, ncinfo , ncread...Once you read the netCDf files. YOu can write the data into geotiff using geotiffwrite.

Read about:

http://in.mathworks.com/help/matlab/ref/ncread.html

http://in.mathworks.com/help/map/ref/geotiffwrite.html

[Gahana Jain]

Yes, I tried reading it but the issue is how to provide georeferencing in geotiff function. I am not able to extract the georeferencing from the nc file.

[KSSV]

this depends on what info your nc file has....attach one nc file here....let me see...

[Gahana Jain]

Hereby attaching a nc file.

[KSSV]

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

It works for all the nc file sin the folder....it creates tiff with the same as netcdf...I am not sure of the variable name though read in A..

[Mian]

Hi KSSC,

I have used your code to convert netcdf into tiff. When I verfied converted file in ArcGIS with orginal netcdf file, both files did not match each other. Then I run same code but just took transpose of "A" i.e [ A = A.']. Now, newly developed tiff is same as of netcdf.

[Preet Lal]

Hello,

I am getting error using above code

Error using map.rasterref.GeographicCellsReference.set.LatitudeLimits

Expected latlim to be of size 1x2, but it is of size 1x26896.

Error in map.rasterref.GeographicRasterReference/set.LatitudeLimits (line 964)

validateattributes(latlim, ...

Error in setSpatialReferencingProperties (line 38)

R.(name) = value;

Error in georasterref (line 111)

R = setSpatialReferencingProperties(R, ...