Question

我有一个栅格格式（lu.tif）的土地利用数据库（纬度/经度），我已经用ldal投影中的gdal对其进行了如下投影：

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 ";

问题在于它无法正确识别球形基准。在spatial_ref行中，似乎是：“ SPHEROID [\” unnamed \“，6370000,0]”。第二半径也应为6370000而不是0。

对此有任何帮助吗？

gdal无法重新格式化lcc中的栅格文件netcdf网格

0 个答案: