I am attempting to warp a GTiff from a lat/lon projection to a higher (~2x) resolution Polar Stereographic grid. 

In the output GTiff any cells north of 64 deg lat become NaN instead of their actual values.

I am using gdalwarp thusly:

```py
  gdalwarp -t_srs 'wkt string here' -tr 5000 5000 -te -4252339.5,-5302590.0,4597660.5,2747410.25 in.tif out.tif

```

What could be causing the output raster to have the blank area?

Input extents are [-84.84, -50.19, 55.07, 79.96]. [minLon,maxLon,minLat,maxLat]

Output Coordinate system:

proj4 (+proj=stere +lat\_0=90 +lat\_ts=60 +lon\_0=-100 +k=1 +x\_0=0 +y\_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no\_defs) . 

It's standard Polar Stereographic but with Central Meridian at -100 deg and latitude of origin at 60 deg.

---

The issue was caused by the large difference in coordinate systems. (The domain was self-intersecting in one of the coord systems). It was largely solved by explicitly entering the output extent instead of having GDAL attempt to guess it.

i.e *-te -180 0 180 90* was added to the gdalwarp line.