My GDAL package version: GDAL 1.9dev, released 2011/01/18 My input file is a NETCDF file which contains global average temperatures. This is the gdalinfo spew:
Driver: netCDF/Network Common Data Format Files: gfdl_cm2_1-sresa2-tas-run1-2041-2070.nc Size is 144, 88 Coordinate System is `' Origin = (0.000000000000000,88.988767426589433) Pixel Size = (2.500000000000000,-2.022471986967942) Metadata: NC_GLOBAL#title=GFDL CM2.1, SRESA2 (run1) SRES A2 experiment output for IPCC AR4 and US CCSP NC_GLOBAL#institution=NOAA GFDL (US Dept of Commerce / NOAA / Geophysical Fluid Dynamics Laboratory, Princeton, NJ, USA) NC_GLOBAL#source=GFDL_CM2.1 (2004): atmosphere: AM2.1 (am2p13fv, M45L24); ocean: OM3.1 (mom4p1p7_om3p5, tripolar360x200L50); sea ice: SIS; land: LM2; infrastructure: FMS preK release NC_GLOBAL#contact=gfdl.climate.model.i...@noaa.gov NC_GLOBAL#project_id=IPCC Fourth Assessment and US CCSP Projects NC_GLOBAL#table_id=Table A1 (20 September 2004) NC_GLOBAL#experiment_id=SRES A2 experiment NC_GLOBAL#realization=1 NC_GLOBAL#cmor_version=9.600000e-01 NC_GLOBAL#Conventions=CF-1.0 NC_GLOBAL#history=Wed Feb 9 15:43:49 2011: ncks -d lat,1,88 gfdl_cm2_1-sresa2-tas-run1-2041-2070.nc gfdl_cm2_1-sresa2-tas-run1-2041-2070.nc input/atmos.200101-210012.t_ref.nc At 10:21:57 on 01/28/2005, CMOR rewrote data to comply with CF standards and IPCC Fourth Assessment and US CCSP Projects requirements NC_GLOBAL#references=The GFDL Data Portal (http://nomads.gfdl.noaa.gov/) provides access to NOAA/GFDL's publicly available model input and output data sets. From this web site one can view and download data sets and documentation, including those related to the GFDL CM2.1 model experiments run for the IPCC's 4th Assessment Report and the US CCSP. NC_GLOBAL#comment=GFDL experiment name = CM2.1U-H2_SresA2_W1. PCMDI experiment name = SRESA2 (run1). Initial conditions for this experiment were taken from 1 January 2001 of the climate of the 20th century experiment named CM2.1U-D4_1860-2000-AllForc_H2. Several time varying forcing agents used during the first 100 years of this CM2.1U-H2_SresA2_W1 experiment were based on output of the MOZART global atmospheric chemistry climate model when driven with the SRES A2 emission scenarios. For the CM2.1U-H2_SresA2_W1 experiment volcanic aerosols were zero, land surface cover types were time invariant, and solar irradiance varied seasonally based upon late 20th century averages but with no interannual variation. NC_GLOBAL#gfdl_experiment_name=CM2.1U-H2_SresA2_W1 NC_GLOBAL#NCO=20110209 tas#standard_name=air_temperature tas#long_name=Surface Air Temperature tas#units=K tas#cell_methods=time: mean tas#coordinates=height tas#original_name=t_ref lat#standard_name=latitude lat#long_name=latitude lat#units=degrees_north lat#axis=Y lat#bounds=lat_bnds lon#standard_name=longitude lon#long_name=longitude lon#units=degrees_east lon#axis=X lon#bounds=lon_bnds Corner Coordinates: Upper Left ( 0.0000000, 88.9887674) Lower Left ( 0.0000000, -88.9887674) Upper Right ( 360.000, 88.989) Lower Right ( 360.000, -88.989) Center ( 180.0000000, -0.0000000) Band 1 Block=144x1 Type=Float32, ColorInterp=Undefined NoData Value=9.96920996838686905e+36 Metadata: NETCDF_VARNAME=tas NETCDF_DIMENSION_timeofyear=1 this is the command I am trying to perform: gdalwarp -s_srs "+proj=longlat" -t_srs "+proj=somerc +lat_0=47.5 +lon_0=-97.0 +y_0=3150000.0 +x_0=3900000.0 +a=6371229 +b=6371229" -tr 50000 -50000 -te 575000.000 575000.000 7275000.000 5775000.000 -srcnodata "1e20" -dstnodata "1e20" -r bilinear -overwrite NETCDF:gfdl_cm2_1-sresa2-tas-run1-2041-2070.nc:tas test.tif as per this thread: http://tinyurl.com/4enmf2o* *I have also tried: gdalwarp --config CENTER_LONG 180 -s_srs "+proj=longlat" -t_srs "+proj=somerc +lat_0=47.5 +lon_0=-97.0 +y_0=3150000.0 +x_0=3900000.0 +a=6371229 +b=6371229" -tr 50000 -50000 -te 575000.000 575000.000 7275000.000 5775000.000 -srcnodata "1e20" -dstnodata "1e20" -r bilinear -overwrite NETCDF:gfdl_cm2_1-sresa2-tas-run1-2041-2070.nc:tas test.tif and with a negative center_long gdalwarp --config CENTER_LONG -180 -s_srs "+proj=longlat" -t_srs "+proj=somerc +lat_0=47.5 +lon_0=-97.0 +y_0=3150000.0 +x_0=3900000.0 +a=6371229 +b=6371229" -tr 50000 -50000 -te 575000.000 575000.000 7275000.000 5775000.000 -srcnodata "1e20" -dstnodata "1e20" -r bilinear -overwrite NETCDF:gfdl_cm2_1-sresa2-tas-run1-2041-2070.nc:tas test.tif But the stays the same: ERROR 1: Too many points (440 out of 441) failed to transform, unable to compute output bounds. I am not sure if it makes a different, but I have slightly altered the input file by removing the top 2 rows of latitude using the command: ncks -d lat,1,88 gfdl_cm2_1-sresa2-tas-run1-2041-2070.nc Any help with this error would be very appreciated, Hailey
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
