gdal неправильно переформатирует растровый файл, сетку netcdf в lcc - PullRequest
Новая
       31

gdal неправильно переформатирует растровый файл, сетку netcdf в lcc

0 голосов
/ 20 декабря 2018

У меня есть база данных землепользования (широта / долгота) в растровом формате (lu.tif), которую я уже перепроецировал с помощью gdal в проекции lcc следующим образом: 

gdalwarp -s_srs "EPSG:4326" -t_srs "+proj=lcc +lat_1=30 +lat_2=60 +lat_0=26 +lon_0=13 +x_0=0 +y_0=0 +a=6370000 +b=6370000 +units=m +no_defs" -tr 1000 1000 -te xmin amin xmax ymaxlu.tif lu_reproj.tif

Моя проекция lcc имеетсферические данные.

Я переформатирую растр в файл netcdf следующим образом: 

gdal_translate -unscale -of netCDF lu_reproj.tif lu_reproj.nc

Сгенерированный файл netcdf имеет следующие параметры: 

char lambert_conformal_conic;
  :grid_mapping_name = "lambert_conformal_conic";
  :longitude_of_central_meridian = 13.0; // double
  :false_easting = 0.0; // double
  :false_northing = 0.0; // double
  :latitude_of_projection_origin = 26.0; // double
  :standard_parallel = 30.0, 60.0; // double
  :longitude_of_prime_meridian = 0.0; // double
  :semi_major_axis = 6370000.0; // double
  :inverse_flattening = 0.0; // double
  :spatial_ref = "PROJCS[\"unnamed\",GEOGCS[\"unnamed ellipse\",DATUM[\"unknown\",SPHEROID[\"unnamed\",6370000,0]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]],PROJECTION[\"Lambert_Conformal_Conic_2SP\"],PARAMETER[\"standard_parallel_1\",30],PARAMETER[\"standard_parallel_2\",60],PARAMETER[\"latitude_of_origin\",26],PARAMETER[\"central_meridian\",13],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]]]";
  :GeoTransform = "3574465 1000 0 1034651 0 -1000 ";

Проблема заключается вчто он не распознает правильно сферические данные.В строкеatial_ref это выглядит так: "SPHEROID [\" unonym \ ", 6370000,0]".Второй радиус также должен быть 6370000, а не 0.

Любая помощь с этим?

...