<div dir="ltr"><div><br></div><div>I am trying to convert geographic coordinates (in degrees east/north lat lon) to projected coordinates in meters for geostationary full disc data (GVAR Goes East data from CLASS).  To do this, I tried using PyProj/ Proj4 and am running into errors due to the off earth pixels in the corners. <br>
</div><div><br></div><div>This is the proj4 string given:</div><div> projection_coords = Proj(&quot;+proj=geos +h=35774290 +a= 6378137 +b= 6378137 +lon_0=-75 +units=meters +no_defs&quot;)</div><div>ll_x, ll_y = projection_coords_geos( LL_y_deg, LL_x_deg, inverse = False, errcheck=True)</div>
<div>    &#39;x=%9.3f y=%11.3f&#39; % (ll_x,ll_y)</div><div>    ur_x,ur_y = projection_coords_geos( UR_y_deg, UR_x_deg, inverse = False, errcheck=True)</div><div> &#39;x=%9.3f y=%11.3f&#39; % (ur_x,ur_y)</div><div><br></div>
<div>For full disc geos data I am receiving incorrect (1e30) values for the extents and was wondering if anyone has ever run into the same error or have a suggestion to correct this. </div><div><br></div><div>Also, for remapping GOES data, should the &#39;sweep axis&#39; be including in the Proj string and set to &#39;y&#39;? </div>
<div><br></div><div>Thanks in advance for any suggestions!</div></div>