[gdal-dev] ERROR 1: Too many points (440 out of 441) failed to transform, unable to compute output bounds.

[Hailey Eckstrand]

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.Info at 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20110210/b588e861/attachment.html