[GRASS-dev] [GRASS GIS] #3356: v.to.db: incorrect area calculations in lat-long location

GRASS GIS trac at osgeo.org
Fri Jun 2 09:03:39 PDT 2017


#3356: v.to.db: incorrect area calculations in lat-long location
--------------------------+-----------------------------------
  Reporter:  mlennert     |      Owner:  grass-dev@…
      Type:  defect       |     Status:  new
  Priority:  normal       |  Milestone:  7.4.0
 Component:  Vector       |    Version:  svn-trunk
Resolution:               |   Keywords:  v.to.db area lat-long
       CPU:  Unspecified  |   Platform:  Unspecified
--------------------------+-----------------------------------

Comment (by hellik):

 Replying to [comment:1 hellik]:
 > Replying to [ticket:3356 mlennert]:
 > > {{{
 > > $area in field calculator in QGIS:
 > > 425.1112801
 > > }}}
 >
 > just tested it with
 >
 > {{{
 > QGIS version 2.18.9 (OTF enabled)
 > }}}
 >
 > {{{
 > 1     34.815138348493
 > }}}
 >
 > and when I do manual measure in QGIS (OTF on)
 >
 > {{{
 > 32,4 m²
 > }}}

 now tested

 {{{
 g.region -p vector=training_single at test
 projection: 3 (Latitude-Longitude)
 zone:       0
 datum:      wgs84
 ellipsoid:  wgs84
 north:      0:13:32.187577N
 south:      0:13:32.175318N
 west:       73:12:48.504873E
 east:       73:12:48.518218E
 nsres:      0:00:00.012259
 ewres:      0:00:00.013345
 rows:       1
 cols:       1
 cells:      1
 }}}

 {{{
 v.in.region output=reg_training_single
 }}}

 {{{
 v.report map=reg_training_single at test option=area units=meters
 cat|area
 1|0.155378405506674
 }}}

 then export reg_training_single to a shapefile, reprojet this shapefile to
 EPSG:32743, creating a EPSG:32743-location by ogr2ogr, import the
 reprojected shapefile into the new location

 {{{
 v.report map=myreg2_32743 at data option=area units=meters
 cat||area
 1|0.155406186822802
 }}}

 now area of the bbox of training_single in wgs84 and EPSG:32743 seems to
 be similar.

 check this bbox vector area also in qgis (OTF enabled).

 {{{
 cat     areawgs84       areareproj
 1       0.1574797419    0.157479741
 }}}

 QGIS and GRASS area seems similar and reasonable, when you're looking at
 the kml-file in real life conditions.

 these numbers are similar to the number in the original report

 {{{
 v.report training_reproj_grass op=area u=me
 cat|id|area
 1|1|0.126369981910102
 }}}

 but not similar to

 {{{
 v.report training op=area u=me
 cat|id|area
 1|1|38.2256243922775
 }}}

 or

 {{{
 $area in field calculator in QGIS:
 425.1112801
 }}}

 any idea?

--
Ticket URL: <https://trac.osgeo.org/grass/ticket/3356#comment:3>
GRASS GIS <https://grass.osgeo.org>



More information about the grass-dev mailing list