[gdal-dev] MSG-CPP satellite projection to wgs84 with gdal

Tibor Szabo szabot.hu at gmail.com
Thu Jul 16 12:00:40 PDT 2020


Hi All,

I have tried to wrap a netcdf file for changing projection from satellite
view to EPSG:4326 with gdal_warp in Ubuntu and windows also with error.
Sample file:
https://z003.s3-us-west-2.amazonaws.com/SEVIR_OPER_R___MSGCPP__L2__20200712T220000_20200712T221500_0001.nc
D:\meteo_model>gdalwarp -of NetCDF -s_srs "+proj=geos +lon_0=0.000000
+lat_0=0 +h=35807.414063 +a=6378.169 +b=6356.5838" -t_srs EPSG:4326 -te 0
50 10 60
NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":cot
TEST_MOD_COT_NL.nc
ERROR 5: OSRCalcInvFlattening(): Wrong input values
ERROR 1: PROJ: proj_create_operations: Source and target ellipsoid do not
belong to the same celestial body
ERROR 6: Cannot find coordinate operations from
`PROJCRS["unknown",BASEGEOGCRS["unknown",DATUM["unknown",ELLIPSOID["unknown",6378.169,295.48806589701,LENGTHUNIT["metre",1,ID["EPSG",9001]]]],PRIMEM["Reference
meridian",0,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]]],CONVERSION["unknown",METHOD["Geostationary
Satellite (Sweep Y)"],PARAMETER["Longitude of natural
origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Satellite
Height",35807.414063,LENGTHUNIT["metre",1,ID["EPSG",9001]]],PARAMETER["False
easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False
northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'
to `EPSG:4326'

D:\meteo_model>

My gdal version:
Win:
D:\meteo_model>ogrinfo --version
GDAL 3.1.0dev, released 2019/99/99
Linux:
szabot at manicminer:~/gdal_knmi_test$ ogrinfo --version
GDAL 3.0.4, released 2020/01/28

BR
Tibor


D:\meteo_model>gdalinfo
SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc
Driver: netCDF/Network Common Data Format
Files: SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc
Size is 512, 512
Metadata:
  NC_GLOBAL#comment=none
  NC_GLOBAL#Conventions=CF-1.4
  NC_GLOBAL#history=NetCDF files based on KNMI MSG-CPP HDF5 files.
  NC_GLOBAL#institution=Royal Netherlands Meteorological Institute (KNMI)
  NC_GLOBAL#references=http://msgcpp.knmi.nl
  NC_GLOBAL#source=Royal Netherlands Meteorological Institute (KNMI)
  NC_GLOBAL#title=MSG-CPP derived Cloud Physical Products from Meteosat
Second Generation (MSG)
  NC_GLOBAL#version=0001
Subdatasets:

SUBDATASET_1_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":azidiff
  SUBDATASET_1_DESC=[1x3712x3712]
angle_of_rotation_from_solar_azimuth_to_platform_azimuth (16-bit integer)

SUBDATASET_2_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":cth
  SUBDATASET_2_DESC=[1x3712x3712] height_at_cloud_top (16-bit integer)

SUBDATASET_3_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":cldmask
  SUBDATASET_3_DESC=[1x3712x3712] cloud_area_fraction status_flag (16-bit
integer)

SUBDATASET_4_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":cot
  SUBDATASET_4_DESC=[1x3712x3712] atmosphere_optical_thickness_due_to_cloud
(16-bit integer)

SUBDATASET_5_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":cph
  SUBDATASET_5_DESC=[1x3712x3712] thermodynamic_particle_phase_at_cloud_top
status_flag (16-bit integer)

SUBDATASET_6_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":ctt
  SUBDATASET_6_DESC=[1x3712x3712] air_temperature_at_cloud_top (16-bit
integer)

SUBDATASET_7_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":cwp
  SUBDATASET_7_DESC=[1x3712x3712] atmosphere_cloud_condensed_water_content
(16-bit integer)

SUBDATASET_8_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":dcld
  SUBDATASET_8_DESC=[1x3712x3712] thickness_of_liquid_water_cloud (16-bit
integer)

SUBDATASET_9_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":dcot
  SUBDATASET_9_DESC=[1x3712x3712] atmosphere_optical_thickness_due_to_cloud
standard_error (16-bit integer)

SUBDATASET_10_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":dcwp
  SUBDATASET_10_DESC=[1x3712x3712] atmosphere_cloud_condensed_water_content
standard_error (16-bit integer)

SUBDATASET_11_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":dndv
  SUBDATASET_11_DESC=[1x3712x3712]
number_concentration_of_cloud_liquid_water_particles_in_air_at_liquid_water_cloud_top
(16-bit integer)

SUBDATASET_12_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":dreff
  SUBDATASET_12_DESC=[1x3712x3712]
effective_radius_of_cloud_condensed_water_particles_at_cloud_top
standard_error (16-bit integer)

SUBDATASET_13_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":precip
  SUBDATASET_13_DESC=[1x3712x3712] lwe_precipitation_rate (16-bit integer)

SUBDATASET_14_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":precip_ir
  SUBDATASET_14_DESC=[1x3712x3712] lwe_precipitation_rate (16-bit integer)

SUBDATASET_15_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":qa
  SUBDATASET_15_DESC=[1x3712x3712] *** TBD *** (16-bit integer)

SUBDATASET_16_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":reff
  SUBDATASET_16_DESC=[1x3712x3712]
effective_radius_of_cloud_condensed_water_particles_at_cloud_top (16-bit
integer)

SUBDATASET_17_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":satz
  SUBDATASET_17_DESC=[1x3712x3712] platform_zenith_angle (16-bit integer)

SUBDATASET_18_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":sds
  SUBDATASET_18_DESC=[1x3712x3712] surface_downwelling_shortwave_flux
(16-bit integer)

SUBDATASET_19_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":sds_cs
  SUBDATASET_19_DESC=[1x3712x3712]
surface_downwelling_shortwave_flux_assuming_clear_sky (16-bit integer)

SUBDATASET_20_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":sds_diff
  SUBDATASET_20_DESC=[1x3712x3712]
surface_diffuse_downwelling_shortwave_flux_in_air (16-bit integer)

SUBDATASET_21_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":sds_diff_cs
  SUBDATASET_21_DESC=[1x3712x3712]
surface_diffuse_downwelling_shortwave_flux_in_air_assuming_clear_sky
(16-bit integer)

SUBDATASET_22_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":sunz
  SUBDATASET_22_DESC=[1x3712x3712] solar_zenith_angle (16-bit integer)

SUBDATASET_23_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":lat
  SUBDATASET_23_DESC=[3712x3712] latitude (32-bit floating-point)

SUBDATASET_24_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":lon
  SUBDATASET_24_DESC=[3712x3712] longitude (32-bit floating-point)

SUBDATASET_25_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":time_offset
  SUBDATASET_25_DESC=[3712x3712] time_offset (16-bit integer)

SUBDATASET_26_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":version_id
  SUBDATASET_26_DESC=[1x1] version_id (64-bit floating-point)

SUBDATASET_27_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":sat_id
  SUBDATASET_27_DESC=[1x1] sat_id (64-bit floating-point)

SUBDATASET_28_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":earth_sun_distance
  SUBDATASET_28_DESC=[1x1] earth_sun_distance (64-bit floating-point)

SUBDATASET_29_NAME=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":annual_mean_solar_extraterrestrial_flux
  SUBDATASET_29_DESC=[1x1] annual_mean_solar_extraterrestrial_flux (64-bit
floating-point)
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  512.0)
Upper Right (  512.0,    0.0)
Lower Right (  512.0,  512.0)
Center      (  256.0,  256.0)

D:\meteo_model>

------------------------------------------------------------

D:\meteo_model>gdalinfo
NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":cot
ERROR 5: OSRCalcInvFlattening(): Wrong input values
Driver: netCDF/Network Common Data Format
Files: none associated
Size is 3712, 3712
Coordinate System is:
GEOGCRS["unknown",
    DATUM["unnamed",
        ELLIPSOID["Spheroid",6356584,0,
            LENGTHUNIT["metre",1,
                ID["EPSG",9001]]]],
    PRIMEM["Greenwich",0,
        ANGLEUNIT["degree",0.0174532925199433,
            ID["EPSG",9122]]],
    CS[ellipsoidal,2],
        AXIS["latitude",north,
            ORDER[1],
            ANGLEUNIT["degree",0.0174532925199433,
                ID["EPSG",9122]]],
        AXIS["longitude",east,
            ORDER[2],
            ANGLEUNIT["degree",0.0174532925199433,
                ID["EPSG",9122]]]]
Data axis to CRS axis mapping: 2,1
Origin = (-5568.748339200000373,5568.748339200000373)
Pixel Size = (3.000403200000000,-3.000403200000000)
Metadata:
  cot#coordinates=lat lon
  cot#grid_mapping=projection
  cot#long_name=Cloud Optical Thickness
  cot#offset=0
  cot#scale_factor=0.0099999998
  cot#standard_name=atmosphere_optical_thickness_due_to_cloud
  cot#units=1
  cot#_FillValue=-1
  NC_GLOBAL#comment=none
  NC_GLOBAL#Conventions=CF-1.4
  NC_GLOBAL#history=NetCDF files based on KNMI MSG-CPP HDF5 files.
  NC_GLOBAL#institution=Royal Netherlands Meteorological Institute (KNMI)
  NC_GLOBAL#references=http://msgcpp.knmi.nl
  NC_GLOBAL#source=Royal Netherlands Meteorological Institute (KNMI)
  NC_GLOBAL#title=MSG-CPP derived Cloud Physical Products from Meteosat
Second Generation (MSG)
  NC_GLOBAL#version=0001
  NETCDF_DIM_EXTRA={time}
  NETCDF_DIM_time_DEF={1,6}
  NETCDF_DIM_time_VALUES=1594628100
  projection#EPSG_code=none
  projection#grid_mapping_name=MSGnavigation
  projection#height_from_earth_center=42163972
  projection#latitude_of_projection_origin=0
  projection#longitude_of_projection_origin=0
  projection#long_name=Projection information
  projection#proj4_params=+proj=geos +lon_0=0.000000 +lat_0=0
+h=35807.414063 +a=6378.169 +b=6356.5838
  projection#scale_x=35785.832
  projection#scale_y=-35785.832
  projection#semi_major_axis=6356584
  projection#semi_minor_axis=6378169
  time#calendar=none
  time#long_name=time
  time#standard_name=time
  time#units=seconds since 1970-01-01 00:00:00
  x#long_name=x coordinate of projection
  x#standard_name=projection_x_coordinate
  x#units=km
  y#long_name=y coordinate of projection
  y#standard_name=projection_y_coordinate
  y#units=km
Geolocation:
  LINE_OFFSET=0
  LINE_STEP=1
  PIXEL_OFFSET=0
  PIXEL_STEP=1
  SRS=GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS
84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]
  X_BAND=1

X_DATASET=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":lon
  Y_BAND=1

Y_DATASET=NETCDF:"SEVIR_OPER_R___MSGCPP__L2__20200713T081500_20200713T083000_0001.nc":lat
Corner Coordinates:
Upper Left  (   -5568.748,    5568.748) (Invalid angle,Invalid angle)
Lower Left  (   -5568.748,   -5568.748) (Invalid angle,Invalid angle)
Upper Right (    5568.748,    5568.748) (Invalid angle,Invalid angle)
Lower Right (    5568.748,   -5568.748) (Invalid angle,Invalid angle)
Center      (  -0.0000000,   0.0000000) (  0d 0' 0.00"W,  0d 0' 0.00"N)
Band 1 Block=3712x3712 Type=Int16, ColorInterp=Undefined
  NoData Value=-1
  Unit Type: 1
  Offset: 0,   Scale:0.00999999977648258
  Metadata:
    coordinates=lat lon
    grid_mapping=projection
    long_name=Cloud Optical Thickness
    NETCDF_DIM_time=1594628100
    NETCDF_VARNAME=cot
    offset=0
    scale_factor=0.0099999998
    standard_name=atmosphere_optical_thickness_due_to_cloud
    units=1
    _FillValue=-1

--------------------------------------------------



>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20200716/73681d52/attachment-0001.html>


More information about the gdal-dev mailing list