<div dir="ltr">I am trying to re-project some landsat 8 bands in from one location into another location with the same Ordnance Survey Projection (epsg:27700) using r.proj. They should overlap but currently do not. Am i using the correct module? The error message is:<br>
<br>--------------------------------------------------------------------------<br>(Tue Apr 29 11:30:48 2014) <br>r.proj input=LC82040242013363LGN00_refl1 location=Landsat8 mapset=Landsat dbase=\\msastore\decsgroup\GIS Admin support share output=LC82040242013363LGN00_refl1<br>
Input Projection Parameters: +proj=tmerc +lat_0=23<br>+lon_0=-96 +k=1 +x_0=0 +y_0=0 +no_defs +a=6377563.396<br>+rf=299.3249646 +towgs84=368.000,-120.000,425.000<br>Input Unit Factor: 1<br>Output Projection Parameters: +proj=tmerc +lat_0=49<br>
+lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +no_defs<br>+a=6377563.396 +rf=299.3249646<br>+towgs84=368.000,-120.000,425.000<br>Output Unit Factor: 1<br>WARNING: pj_transform() failed: latitude or longitude exceeded limits<br>
WARNING: pj_transform() failed: latitude or longitude exceeded limits<br>WARNING: pj_transform() failed: latitude or longitude exceeded limits<br>----------------------------------------------------------------------<br><br>
<br>r.proj parameters<br><br><br>From Projection<br><br><br><br>projection: 99 (Transverse Mercator)<br>zone: 0<br>datum: osgb36<br>ellipsoid: airy<br>north: 5848830<br>south: 5605470<br>west: 296070<br>
east: 536430<br>nsres: 30<br>ewres: 30<br>rows: 8112<br>cols: 8012<br>cells: 64993344<br><br><br><br>to Projection<br><br>projection: 99 (Transverse Mercator)<br>zone: 0<br>datum: osgb36<br>
ellipsoid: airy<br>north: 1230000<br>south: 0<br>west: 0<br>east: 660000<br>nsres: 50<br>ewres: 50<br>rows: 24600<br>cols: 13200<br>cells: 324720000<br clear="all"><div><br>
-- <br><div dir="ltr"><b>Paul J. Shapley</b><span style="font-family:times new roman,serif"><font><span style="font-family:arial,helvetica,sans-serif"></span></font><br></span><b><br></b></div>
</div></div>