[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