Как эффективно конвертировать спутниковые netcdf файлы Sentinel3A в изображения геотифов в Matlab? - PullRequest
Новая
       71

Как эффективно конвертировать спутниковые netcdf файлы Sentinel3A в изображения геотифов в Matlab?

0 голосов
/ 13 июня 2019

Я получил файлы Sentinel 3A NetCDF онлайн, и я хочу преобразовать их в geotiff с широтой и долготой.

При чтении файлов NetCDF в Matlab я понял, что в моих файлах NetCDF диапазонов яркости есть не переменные широты / долготы, а строки и столбцы и размер. Есть еще один бесплатный файл (tie_geo_coordinates.nc), в котором переменные широта / долгота, но долгота отрицательная, а широта положительная, что должно быть наоборот. Я получил значения lat / long в приведенном ниже коде, найдя max / min каждого в tie_geo_coordinates.nc. 

file = dir('Oa14_radiance.nc') ; 
N = length(file) ;
ncfile = file.name ;
  A = ncread(ncfile,'Oa14_radiance') ;
  % Write  into geotiff
  B = imrotate(A,90);
  R = georasterref('RasterSize',size(B),...
  'LongitudeLimits',[60.4556,73.8768],'LatitudeLimits', 
    [-67.5003,-29.9244]);
  tiffile = strcat(ncfile,'.tif') ;
  geotiffwrite(tiffile,B,R);

Когда я пытаюсь использовать приведенный ниже код, при открытии в Arc изображения искажаются, а широта / долгота оказываются совершенно неправильными и не пропорциональными. У кого-нибудь есть похожая проблема?

...