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@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 -------------------------------------------------- >
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
