[GRASS-dev] Re: [GRASS GIS] #352: r.in.gdal region troubles in LL

[GRASS GIS]

#352: r.in.gdal region troubles in LL
--------------------------+-------------------------------------------------
  Reporter:  neteler      |       Owner:  grass-dev at lists.osgeo.org
      Type:  defect       |      Status:  new                      
  Priority:  major        |   Milestone:  6.4.0                    
 Component:  default      |     Version:  unspecified              
Resolution:               |    Keywords:                           
  Platform:  Unspecified  |         Cpu:  Unspecified              
--------------------------+-------------------------------------------------
Changes (by neteler):

 * cc: Agustin.Lobo at ija.csic.es (added)

Comment:

 Replying to [comment:1 glynn]:
 > Replying to [ticket:352 neteler]:
 > > gdalinfo npp.tif
 >
 > > Pixel Size = (0.250000000000200,-0.250000000000200)
 >
 > It's that extra 0.0000000000002 that's causing the problems:

 I see.

 ...
 >
 > [And which will forever cause problems. We'll just have to work around
 them as we find them. Almost everything that you "know" about coordinate
 geometry is wrong for spherical geometry. E.g. just because x1 == x2, that
 doesn't mean that x1-x0 == x2-x0.]

 The user took this map

 http://user.uni-
 frankfurt.de/~grieser/downloads/NetPrimaryProduction/npp_CruP_All_fine_geo.zip

 and passed it to R, and then redefined:

 {{{
 >b <-  npp_geotiff
 >b at data$band1[b at data$band1  <0 ]<- 0
 }}}

 and then he wrote b as npp.tif using the QGIS plugin manageR.

 (Source: http://lists.osgeo.org/pipermail/grass-user/2008-
 October/047211.html)

 Maybe there is a precision issue before, somwhere in the R or QGIS part?

 Markus

-- 
Ticket URL: <http://trac.osgeo.org/grass/ticket/352#comment:2>
GRASS GIS <http://grass.osgeo.org>