<div dir="ltr">Thanks Janne. <div><br></div><div>We did try drawing a limiting box inside the earth extents to get all valid data points within, and this works but cuts out data in the edges of the earth. See the following image: </div>
<div><br></div><div>
<img src="cid:ii_140781440ec30d01" alt="Inline image 1"><br></div><div>We still need the data in between the drawn dotted line blue box and the edge of the earth. The data values in the red shaded off earth areas are very large (2.14 E9) and these are the values that throw the error. If you do discard the values in the light red, that is fine but it is throwing an exception instead, preventing us from proceeding forward. </div>
<div><br></div><div>Thanks for the support!</div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Aug 13, 2013 at 3:12 AM, <span dir="ltr"><<a href="mailto:support.mn@elisanet.fi" target="_blank">support.mn@elisanet.fi</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
<br>
we just discarded all points outside any reasonable area.<br>
Another possibility might be to draw a limiting line (or box)<br>
and move all outside points to that..<br>
<br>
Hope that helps?<br>
<br>
Regards: Janne.<br>
<br>
----------------------------------------------------------------------------------------<br>
<br>
Kathryn Jablonski - NOAA Affiliate [<a href="mailto:kathryn.jablonski@noaa.gov">kathryn.jablonski@noaa.gov</a>] kirjoitti:<br>
<div><div class="h5">> I am trying to convert geographic coordinates (in degrees east/north lat<br>
> lon) to projected coordinates in meters for geostationary full disc data<br>
> (GVAR Goes East data from CLASS). To do this, I tried using PyProj/ Proj4<br>
> and am running into errors due to the off earth pixels in the corners.<br>
><br>
> This is the proj4 string given:<br>
> projection_coords = Proj("+proj=geos +h=35774290 +a= 6378137 +b= 6378137<br>
> +lon_0=-75 +units=meters +no_defs")<br>
> ll_x, ll_y = projection_coords_geos( LL_y_deg, LL_x_deg, inverse = False,<br>
> errcheck=True)<br>
> 'x=%9.3f y=%11.3f' % (ll_x,ll_y)<br>
> ur_x,ur_y = projection_coords_geos( UR_y_deg, UR_x_deg, inverse =<br>
> False, errcheck=True)<br>
> 'x=%9.3f y=%11.3f' % (ur_x,ur_y)<br>
><br>
> For full disc geos data I am receiving incorrect (1e30) values for the<br>
> extents and was wondering if anyone has ever run into the same error or<br>
> have a suggestion to correct this.<br>
><br>
> Also, for remapping GOES data, should the 'sweep axis' be including in the<br>
> Proj string and set to 'y'?<br>
><br>
> Thanks in advance for any suggestions!<br>
><br>
</div></div>> _______________________________________________<br>
> Proj mailing list<br>
> <a href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a><br>
> <a href="http://lists.maptools.org/mailman/listinfo/proj" target="_blank">http://lists.maptools.org/mailman/listinfo/proj</a><br>
<br>
<br>
_______________________________________________<br>
Proj mailing list<br>
<a href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a><br>
<a href="http://lists.maptools.org/mailman/listinfo/proj" target="_blank">http://lists.maptools.org/mailman/listinfo/proj</a><br>
</blockquote></div><br></div>