<HTML><HEAD></HEAD>
<BODY dir=ltr>
<DIV dir=ltr>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: 'Calibri'; COLOR: #000000">
<DIV>Your graphic shows the earth profile as an ellipse, but your proj4 string
sets a and b to the semi-major axis of WGS84 (i.e. a sphere). Why?
-Noel</DIV>
<DIV> </DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: 'Calibri'; COLOR: #000000">Noel Zinn,
Principal, Hydrometronics LLC<BR>+1-832-539-1472 (office), +1-281-221-0051
(cell)<BR>noel.zinn@hydrometronics.com (email)<BR>http://www.hydrometronics.com
(website)<BR></DIV>
<DIV
style="FONT-SIZE: small; FONT-FAMILY: 'Calibri'; FONT-WEIGHT: normal; COLOR: #000000; FONT-STYLE: normal; TEXT-DECORATION: none; DISPLAY: inline">
<DIV style="FONT: 10pt tahoma">
<DIV> </DIV>
<DIV style="BACKGROUND: #f5f5f5">
<DIV style="font-color: black"><B>From:</B> <A title=kathryn.jablonski@noaa.gov
href="mailto:kathryn.jablonski@noaa.gov">Kathryn Jablonski - NOAA Affiliate</A>
</DIV>
<DIV><B>Sent:</B> Tuesday, August 13, 2013 9:38 AM</DIV>
<DIV><B>To:</B> <A title=proj@lists.maptools.org
href="mailto:proj@lists.maptools.org">PROJ.4 and general Projections
Discussions</A> </DIV>
<DIV><B>Subject:</B> Re: [Proj] Proj4 Projection Coordinates for Geostationary
satellitedata</DIV></DIV></DIV>
<DIV> </DIV></DIV>
<DIV
style="FONT-SIZE: small; FONT-FAMILY: 'Calibri'; FONT-WEIGHT: normal; COLOR: #000000; FONT-STYLE: normal; TEXT-DECORATION: none; DISPLAY: inline">
<DIV dir=ltr>Thanks Janne.
<DIV> </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> </DIV>
<DIV><IMG alt="Inline image 1"
src="cid:0D47AF6E28124565B1B92D72971B6C82@NoelHmPC"><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> </DIV>
<DIV>Thanks for the support!</DIV>
<DIV> </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="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">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>
<DIV> </DIV></DIV>
<P>
<HR>
_______________________________________________<BR>Proj mailing
list<BR>Proj@lists.maptools.org<BR>http://lists.maptools.org/mailman/listinfo/proj</DIV></DIV></DIV></BODY></HTML>