<div>Hi Proj4 developers and users,</div><div><br></div><div>I have a polygon in kml file (epsg 4326).<br></div><div>I wanted to calculate its area in sq. km.<br></div><div>I used the following code<br></div><div><br></div><div>poly = transform(partial(pyproj.transform,  pyproj.Proj(init='epsg:4326'),pyproj.Proj(init='epsg:3857')), poly)<br></div><div>print (<a href="http://poly.area/1000000">poly.area/1000000</a>)<br></div><div><br></div><div>However, this result is more than 100 sq.km. larger than the area shown in google earth.<br></div><div>So, what's wrong with transferring into epsg 3857 and calculating area in shapely (the first code)?<br></div><div><br></div><div>Thanks for your ideas.<br></div><div><br></div><div>Isak</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div>