I have nc files. from the metadata ,the projection is cylindrical and the resolution is 25 km:
indrical NC_GLOBAL#ease_resolution=25 NC_GLOBAL#history=none NC_GLOBAL#institution=SMOS CATDS Processing Chain NC_GLOBAL#netcdf_version_id=3.6.2 NC_GLOBAL#product_version=1.0 NC_GLOBAL#references=DDI CATDS - Ref. CAT-DDI-CT-00020-CG NC_GLOBAL#source=CATDS L3 SM Filtrering Processor NC_GLOBAL#title=L3-SM P11p : Preliminary Daily Retrieval Map Corner Coordinates: Upper Left ( 0.0, 0.0) Lower Left ( 0.0, 586.0) Upper Right ( 1383.0, 0.0) Lower Right ( 1383.0, 586.0) Center ( 691.5, 293.0) I want to reproject it to WGS84 with a resolution of 0.25 * 0.25 degree. I used this code:
import os for doy in range(1,367): inputFile='D:\\358\\2010SMOSl3_sm'+str(doy)+'.nc' inputDataset='NETCDF:\"'+inputFile+'\":Soil_Moisture' outputFile='D:\\SMOS\\RE01\\MIR_CLF31D\\2010_sm'+str(doy)+'.img' gdalcmd='gdalwarp -of "ENVI" -ot Float32 -srcnodata -32768 -dstnodata -9999 -te -180 -90 180 90 -overwrite '+inputDataset+' '+outputFile os.system( gdalcmd) but I got an error:
ERROR 1: Unable to compute a transformation between pixel/line and georeferenced coordinates for NETCDF
:\SMOS\RE01\MIR_CLF31D\SM-D-RE-renamed- daily-2010-from day 14 to day 358\2010SMOSl3_sm307.nc:Soil_Moisture. There is no affine transformation and no GCPs.
أكثر...
indrical NC_GLOBAL#ease_resolution=25 NC_GLOBAL#history=none NC_GLOBAL#institution=SMOS CATDS Processing Chain NC_GLOBAL#netcdf_version_id=3.6.2 NC_GLOBAL#product_version=1.0 NC_GLOBAL#references=DDI CATDS - Ref. CAT-DDI-CT-00020-CG NC_GLOBAL#source=CATDS L3 SM Filtrering Processor NC_GLOBAL#title=L3-SM P11p : Preliminary Daily Retrieval Map Corner Coordinates: Upper Left ( 0.0, 0.0) Lower Left ( 0.0, 586.0) Upper Right ( 1383.0, 0.0) Lower Right ( 1383.0, 586.0) Center ( 691.5, 293.0) I want to reproject it to WGS84 with a resolution of 0.25 * 0.25 degree. I used this code:
import os for doy in range(1,367): inputFile='D:\\358\\2010SMOSl3_sm'+str(doy)+'.nc' inputDataset='NETCDF:\"'+inputFile+'\":Soil_Moisture' outputFile='D:\\SMOS\\RE01\\MIR_CLF31D\\2010_sm'+str(doy)+'.img' gdalcmd='gdalwarp -of "ENVI" -ot Float32 -srcnodata -32768 -dstnodata -9999 -te -180 -90 180 90 -overwrite '+inputDataset+' '+outputFile os.system( gdalcmd) but I got an error:
ERROR 1: Unable to compute a transformation between pixel/line and georeferenced coordinates for NETCDF
:\SMOS\RE01\MIR_CLF31D\SM-D-RE-renamed- daily-2010-from day 14 to day 358\2010SMOSl3_sm307.nc:Soil_Moisture. There is no affine transformation and no GCPs. أكثر...