>
> ... but basically I came to the conclusion that when people talk about an
> area of something on Earth they really mean the area of a projection (am i
> right about that?) and so i need to specify a projection before getting a
> meaningful area. So i ended up with this code:
>
> wkt = "POLYGON ((lon1 lat1, lon2 lat2, lon3 lat3))"
> poly = ogr.CreateGeometryFromWkt(wkt)
> sr = osr.SpatialReference()
> sr.SetWellKnownGeogCS("NAD83")
> sr.SetProjCS("SRS_PT_TRANSVERSE_MERCATOR")
> poly.AssignSpatialReference(sr)
> poly.GetArea() <https://lists.osgeo.org/mailman/listinfo/gdal-dev>
>
> I might be dense, but it seems you are assigning sr as a specific
transverse mercator, i.e. UTM zone 18N, or a custom one centered on the
centroid of your polygon. i.e. you don't have sufficient arguments.
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail&utm_term=icon>
Virus-free.
www.avast.com
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail&utm_term=link>
<#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
_______________________________________________
gdal-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
